How does one systematically check if a given configuration of prime numbers $p_1, p_2, ... p_n$ is the densest possible configuration of primes in the range $[p_1, p_n]$? (The densest configurations seem to be called “prime constellations”).
I am writing a paper for my high school science class (not math). It deals with interesting aspects of prime numbers. In this context I have two questions:
1) It is fairly clear that the configuration of a prime triplet must either be (p, p+2, p+6) or (p, p+4, p+6), as the „obvious“ shorter choice (p, p+2, p+4) is an equally spaced sequence of three elements and therefore one of the elements will be divisible by three (the exception of course is a sequence with “3” as its first element). While this simple example of three primes is clear I would like to know how to continue? What is the rule/algorithm for identifying longer prime constellations and how many constellations of the same length exist?
2) If a configuration of primes is given as $p, p+k_1, p+k_2, … p+k_n$, how can I check if this is the densest possible configuration, a so-called “prime constellation“ in the light of the above?
Thank you for your help. I would appreciate if the explanations are hands-on with examples that are fairly easy to understand.
$n,n+2,n+4$ can't be all primes more than once because whatever is $n \bmod 3$ one of $n,n+2,n+4$ will be $\equiv 0 \bmod 3$.
If $p > b_k+1$ there is the solution $a_p=1$.
Thus it suffices to look at the primes $p\le b_k+1$ which means there is no obstruction exactly when for some $A$ all of $A,A+b_1,\ldots,A+b_k$ are coprime with $(b_k+1)!$
The natural generalization of the twin prime conjecture is that when there is no obstruction $n,n+b_1,\ldots,n+b_k$ are all primes infinitely many times.
The conjectured asymptotic (from the random model for the primes) for the number of such $n$ is $$\pi_b(x) \sim \sum_{n \le x} \prod_{p \le n^r} \frac{\sum_{a =1,\ p\, \nmid\, a (a+b_1)\ldots (a+b_k)}^{p-1} 1}{p} \sim C(b) \frac{x}{\ln^{k+1} x}$$ where for $p > b_k+1, \frac{\sum_{a =1,\ p\, \nmid\, a (a+b_1)\ldots (a+b_k)}^{p-1} 1}{p} = 1-\frac{b_k+1}{p}$
and $r =e^{-\gamma} \approx 0.56$ is the constant such that $\prod_{p \le n^r} \frac{p-1}{p} \sim \frac{1}{\ln n}$ (making the predicted $C(b)$ constants compatible with the PNT)
The random model for the primes is merely saying that when choosing $n$ randomly in $[1,x]$, for $p \le n^r$, then $n \bmod p$ can be considered uniformly distributed and more important the events $n \bmod p$ and $n \bmod q$ can be considered independent. Thus, under this model, the probability that $n$ is prime can be considered to be $\prod_{p \le n^r} \frac{p-1}{p}$ which by Mertens's theorem is $\sim \frac1{\ln n}$.