Innumber theory,Brun's theorem states that the sum of thereciprocals of thetwin primes (pairs ofprime numbers which differ by 2)converges to a finite value known asBrun's constant, usually denoted byB2 (sequenceA065421 in theOEIS). Brun's theorem wasproved byViggo Brun in 1919, and it has historical importance in the introduction ofsieve methods.
The convergence of the sum of reciprocals of twin primes follows from bounds on thedensity of the sequence of twin primes.Let denote the number ofprimesp ≤x for whichp + 2 is also prime (i.e. is the number of twin primes with the smaller at mostx). Then, we have
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 asmall set. In explicit terms, the sum
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 anirrational number only if there are infinitely many twin primes.
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 thePentium 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 andPatrick Demichel used all twin primes up to 1016 to give the estimate[2] thatB2 ≈ 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) thatB2 < 2.1754 (assuming theextended Riemann hypothesis). Then in 2025, Lachlan Dunn showedB2 < 2.1609, assuming the generalised Riemann hypothesis.[3] It has been shown unconditionally thatB2 < 2.347.[4]
There is also aBrun's constant for prime quadruplets. Aprime 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 byB4, is the sum of the reciprocals of all prime quadruplets:
with value:
This constant should not be confused with theBrun's constant forcousin primes, as prime pairs of the form (p, p + 4), which is also written asB4. Wolf derived an estimate for the Brun-type sumsBn of 4/n.
Let (sequenceA005597 in theOEIS) be thetwin prime constant. Then it isconjectured that
In particular,
for every and all sufficiently largex.
Many special cases of the above have been proved. Jie Wu proved that for sufficiently largex,
The digits of Brun's constant were used in a bid of $1,902,160,540 in theNortel patent auction. The bid was posted byGoogle and was one of three Google bids based onmathematical constants.[5] Furthermore, academic research on the constant ultimately resulted in thePentium FDIV bug becoming a notablepublic relations fiasco forIntel.[6][7]