Is there any reason to expect the density of primes $p$ such that $2p+1$ is also a prime where $p=3\bmod4$ holds would be different from case of $p=1\bmod4$?


What if $2p+1$ is replaced by $2p-1$ and what is the conjectured density and what is known so far for both cases $2p+1$ and $2p-1$?



Density of Sophie Germain $3\bmod 4$ primes
As with most questions about finding primes in linear forms, nothing is known. We expect that the number of Sophie Germain primes $\leq x$ equivalent to 3 modulo 4 is $$\sim \frac{2x}{\log^2 x} \prod\_{p > 2} (1-2/p) (1-1/p)^{-2},$$ and for the same estimate to hold with 3 replaced by 1. The same should hold true with $2p+1$ replaced by $2p-1$. 


The question of finding Sophie Germian primes equivalent to 3 modulo 4 can be formulated as a sieve theory problem. Start with the set $$A : = \{n \leq x : n \equiv 3 \mod 4\}.$$ Fix a parameter $z \leq \sqrt{x}$. We are interested in the $S(A,z)$, the number of $n \in A$ such that $n(2n+1)$ has no prime factors $\leq z$. For instance, when $z = \sqrt{x}$ then $S(A,z)$ in the number of Sophie Germain primes in $(z,x]$ (plus 1 because of the integer 1).


A [sieve](https://terrytao.wordpress.com/2015/01/29/254a-supplement-5-the-linear-sieve-and-chens-theorem-optional/) will tell us that if $z = x^{o(1)}$, then $$S(A,z) \sim \frac{x}{2} \prod\_{2 < p \leq z} \left(1- \frac{2}{p}\right).$$ If $z = x^{\delta}$ for some fixed $\delta >0$, then a sieve yields $$S(A,z) \gg \frac{x}{\log^2 x}.$$


Current sieves do not allow one to take $z = \sqrt{x}$ and actually do not predict the correct asymptotic formula. For instance, the number of primes $\leq x$ is $\sim x/ \log x$, while a sieve would predict $$x \prod\_{p < \sqrt{x}} (1 - \frac{1}{p}) \sim x 2 e^{- \gamma}/ \log x.$$ This discrepancy is related to the fact that if a prime $p$ divides $n\leq x$ and $p > \sqrt{x}$, then it is impossible that a different prime $q > \sqrt{x}$ also divides $n$. 


To formulate a heuristic, we make the assumption that each element of $A$ is independently prime with probability $1 / \log x$, in accordance with the prime number theorem. Then we consider $\mathbb{P}(n , 2n+ 1 \text{ is prime})$. If these events were independent, then the probability would be $\frac{1}{\log^2x}$. But there are some minor local obstructions to independence that have to be accounted for. The probability that $p>2$ divides $n(2n+1)$ is $1-2/p$, while the probability that $p | n$ times the probability that $p | 2n+1$ is $(1-1/p)^2$. Thus we see we should add the correction faction $$\prod\_{p>2} (1-2/p) (1-1/p)^{-2}.$$ Note that the above product is convergent. For $p =2$, the correction factor is $(1-1/2)^{-2}$, since we assumed all the elements of $A$ are odd. 


Thus we expect the number of Sophie Germain primes equivalent to 3 modulo 4 to be $$\frac{2x}{\log^2 x} \prod\_{p > 2} (1-2/p) (1-1/p)^{-2}.$$ This is known as the twin prime constant and is in the same spirit of the [singular series](http://mathworld.wolfram.com/k-TupleConjecture.html) arising from the the Hardy-Littlewood prime $k$-tuples conjecture. The only difference is that in this question, we seek primes in linear forms $\{a\_i n + b\_i\}\_{i=1}^k$ where $a\_i$ is not required to be 1. State of the art work on finding primes in linear forms can be found [here.](https://arxiv.org/pdf/1405.2593v2.pdf)


