Brun's theorem

From HandWiki
Short description: Theorem that the sum of the reciprocals of the twin primes converges
The convergence to Brun's constant.

In number theory, Brun's theorem states that the sum of the reciprocals of the twin primes (pairs of prime numbers which differ by 2) converges to a finite value known as Brun's constant, usually denoted by B2 (sequence A065421 in the OEIS). Brun's theorem was proved by Viggo Brun in 1919, and it has historical importance in the introduction of sieve methods.

Asymptotic bounds on twin primes

The convergence of the sum of reciprocals of twin primes follows from bounds on the density of the sequence of twin primes. Let [math]\displaystyle{ \pi_2(x) }[/math] denote the number of primes px for which p + 2 is also prime (i.e. [math]\displaystyle{ \pi_2(x) }[/math] is the number of twin primes with the smaller at most x). Then, we have

[math]\displaystyle{ \pi_2(x) = O\!\left(\frac {x(\log\log x)^2}{(\log x)^2} \right)\!. }[/math]

That is, twin primes are less frequent than prime numbers by nearly a logarithmic factor. This bound gives the intuition that the sum of the reciprocals of the twin primes converges, or stated in other words, the twin primes form a small set. In explicit terms, the sum

[math]\displaystyle{ \sum\limits_{ p \, : \, p + 2 \in \mathbb{P} } {\left( {\frac{1}{p} + \frac{1}{{p + 2}}} \right)} = \left( {\frac{1}{3} + \frac{1}{5}} \right) + \left( {\frac{1}{5} + \frac{1}{7}} \right) + \left( {\frac{1}{{11}} + \frac{1}{{13}}} \right) + \cdots }[/math]

either has finitely many terms or has infinitely many terms but is convergent: its value is known as Brun's constant.

If it were the case that the sum diverged, then that fact would imply that there are infinitely many twin primes. Because the sum of the reciprocals of the twin primes instead converges, it is not possible to conclude from this result that there are finitely many or infinitely many twin primes. Brun's constant could be an irrational number only if there are infinitely many twin primes.

Numerical estimates

The series converges extremely slowly. Thomas Nicely remarks that after summing the first billion (109) terms, the relative error is still more than 5%.[1]

By calculating the twin primes up to 1014 (and discovering the Pentium FDIV bug along the way), Nicely heuristically estimated Brun's constant to be 1.902160578.[1] Nicely has extended his computation to 1.6×1015 as of 18 January 2010 but this is not the largest computation of its type.

In 2002, Pascal Sebah and Patrick Demichel used all twin primes up to 1016 to give the estimate[2] that B2 ≈ 1.902160583104. Hence,

Year B2 set of twin
primes below #
by
1976 1.902160540 1 × 1011 Brent
1996 1.902160578 1 × 1014 Nicely
2002 1.902160583104 1 × 1016 Sebah and Demichel

The last is based on extrapolation from the sum 1.830484424658... for the twin primes below 1016. Dominic Klyve showed conditionally (in an unpublished thesis) that B2 < 2.1754 (assuming the extended Riemann hypothesis). It has been shown unconditionally that B2 < 2.347.[3]

There is also a Brun's constant for prime quadruplets. A prime quadruplet is a pair of two twin prime pairs, separated by a distance of 4 (the smallest possible distance). The first prime quadruplets are (5, 7, 11, 13), (11, 13, 17, 19), (101, 103, 107, 109). Brun's constant for prime quadruplets, denoted by B4, is the sum of the reciprocals of all prime quadruplets:

[math]\displaystyle{ B_4 = \left(\frac{1}{5} + \frac{1}{7} + \frac{1}{11} + \frac{1}{13}\right) + \left(\frac{1}{11} + \frac{1}{13} + \frac{1}{17} + \frac{1}{19}\right) + \left(\frac{1}{101} + \frac{1}{103} + \frac{1}{107} + \frac{1}{109}\right) + \cdots }[/math]

with value:

B4 = 0.87058 83800 ± 0.00000 00005, the error range having a 99% confidence level according to Nicely.[1]

This constant should not be confused with the Brun's constant for cousin primes, as prime pairs of the form (pp + 4), which is also written as B4. Wolf derived an estimate for the Brun-type sums Bn of 4/n.

Further results

Let [math]\displaystyle{ C_2=0.6601\ldots }[/math] (sequence A005597 in the OEIS) be the twin prime constant. Then it is conjectured that

[math]\displaystyle{ \pi_2(x) \sim 2C_2\frac{x}{(\log x)^2}. }[/math]

In particular,

[math]\displaystyle{ \pi_2(x) \lt (2C_2+\varepsilon)\frac{x}{(\log x)^2} }[/math]

for every [math]\displaystyle{ \varepsilon\gt 0 }[/math] and all sufficiently large x.

Many special cases of the above have been proved. Most recently, Jie Wu proved that for sufficiently large x,

[math]\displaystyle{ \pi_2(x) \lt 4.5\,\frac{x}{(\log x)^2} }[/math]

where 4.5 corresponds to [math]\displaystyle{ \varepsilon \approx 3.18 }[/math] in the above.

In popular culture

The digits of Brun's constant were used in a bid of $1,902,160,540 in the Nortel patent auction. The bid was posted by Google and was one of three Google bids based on mathematical constants.[4] Furthermore, academic research on the constant ultimately resulted in the Pentium FDIV bug becoming a notable public relations fiasco for Intel.[5][6]

See also

Notes

  1. 1.0 1.1 1.2 Nicely, Thomas R. (18 January 2010). "Enumeration to 1.6*10^15 of the twin primes and Brun's constant". Some Results of Computational Research in Prime Numbers (Computational Number Theory). http://www.trnicely.net/twins/twins2.html. 
  2. Sebah, Pascal; Gourdon, Xavier. "Introduction to twin primes and Brun's constant computation". CiteSeerX 10.1.1.464.1118.
  3. Klyve, Dominic. "Explicit bounds on twin primes and Brun's Constant". https://collections.dartmouth.edu/archive/object/dcdis/dcdis-klyve2007?ctx=dcdis#?length=12&start=0&view=list&rdat_only_u=no&rdat_u=yes&col=dcdis&oc_0=main-title&od_0=a&sv=Brun%27s+constant. 
  4. Damouni, Nadia (1 July 2011). "Dealtalk: Google bid "pi" for Nortel patents and lost". Reuters. https://www.reuters.com/article/us-dealtalk-nortel-google-idUSTRE76104L20110702. 
  5. "Pentium FDIV flaw FAQ". http://www.trnicely.net/pentbug/pentbug.html. 
  6. Price, D. (1995). "Pentium FDIV flaw-lessons learned". IEEE Micro 15 (2): 86–88. doi:10.1109/40.372360. https://ieeexplore.ieee.org/document/372360. 

References

External links