Suppose $C$ is a random subset of $\mathbb{N}\setminus\{1, 2\}$, such that $\forall n \in \mathbb{N}\setminus\{1, 2\}$, $P(n \in C) = \frac{1}{\ln(n)}$ and the events of different numbers belonging to $C$ are independent (Such random subset is called Cramer’s Prime Model, as it was first introduced by Harald Cramer as probabilistic model for distribution of prime numbers). What is the probability that $\exists N \in \mathbb{N}$ such that $\forall n > N$, $2n \in C + C$, where $+$ stands for pointwise summation?
This question was inspired by the famous Goldbach conjecture from number theory:
$\exists N \in \mathbb{N}\,|\,\forall n > N$, $2n \in P + P$, where $P$ is the set of all prime numbers.
That conjecture was first stated in 1742, remains open until now and is considered to be quite hopeless. However, maybe, it is known, what is the probability of similar statement being true in Cramers Model…
This probability is easy to find. For any $k$ from $3$ to $n-1$ the probability that $k$ or $2n-k$ are not in $C$ is $$1-\frac1{\ln k\ln(2n-k)}$$ And the probability, that $n$ is not in $C$ is $$1-\frac1{\ln n}$$ So $$p(2n\in C+C)=1-\left(1-\frac1{\ln n}\right)\prod_{k=3}^{n-1}\left(1-\frac1{\ln k\ln(2n-k)}\right)$$