Euclid showed by a deliciously simple argument that the number of primes is infinite. In a completely different manner, Euler confirmed the same result. Euler’s conclusion followed from his demonstration that the sum of the reciprocals of the primes diverges:

Obviously, this could not happen if there were only finitely many primes.

**Prime Pairs**

Although the prime numbers become sparser among the natural numbers as their size increases, pairs of primes differing by just two continue to appear, albeit with decreasing frequency, as we move to larger and larger numbers. Thus we have prime pairs such as , and .

The question of whether the set of prime pairs is infinite remains unresolved. If we could prove that the sum

diverged, the infinity of the prime pairs would follow. However,* this series actually converges*. In 1919, Norwegian mathematician Viggo Brun proved that

where , now known as Brun’s Constant, has a value of about 1.9. The value of is not known exactly and, indeed, it is not known whether it is rational or irrational (of course, irrationality would imply that the sum is infinite).

The prime number theorem shows that , the number of primes less than or equal to , is asymptotic to :

Let denote the number of primes for which is also prime. Then

Thus, twin primes are fewer than primes by approximately a logarithmic factor. Brun’s Theorem follows from this result. In other words,

is either a finite series, or a convergent infinite series with limit denoted by . However, the convergence of the series does not tell us whether it is a finite or infinite series, so it does not answer the main question:* Are there infinitely many prime pairs*?

Direct computations of with increasing numbers of twin primes converge very slowly (Fig. 1), and more sophisticated extrapolations are used to estimate the ultimate value. Fig. 2 shows that extrapolation indicates a limiting value marginally larger than 1.9. The intersection of the extrapolating line with the vertical axis is Brun’s constant if the twin prime conjecture is valid. According to this approach, the direct estimation of does not exceed 1.9 until the value , far beyond the computational range.

**The Pentium Bug**

In 1995, Thomas Nicely of Lynchburg College, Virginia, showed that . To be on the safe side, Nicely performed all his calculations twice, using two algorithms on two separate computers. He found discrepancies and, after eliminating other possibilities, pinned down the anomalous numerical results to the floating point unit in the brand new Intel Pentium chip, introduced in 1993. The Pentium gave incorrect values for the reciprocals of the twin primes:

with errors in the tenth decimal place. This is a nice example of how problems in number theory can provide excellent test cases for computer hardware and software.

The cause was discovered to be missing entries in a lookup table used for division in the floating point unit of the chip. Intel recalled the defective processors and estimated that the total cost associated with their replacement would be $ 475 million. Although the errors in the faulty chips were rare, the story generated widespread debate, and caused widespread criticism of how Intel had reacted to the problem.

**Sources**

Cipra, Barry A. (1995). “How number theory got the best of the Pentium chip”. *Science*. **267** (5195): 175. doi:10.1126/science.267.5195.175.

Nicely, Thomas (1996): Enumeration to of the Twin Primes and Brun’s Constant, *Virginia J. Sci*., Vol. **46**, p. 195-204

Sebah, Pascal and Xavier Gourdon, 2002: Introduction to twin primes and Brun’s constant computation.

HTML: http://numbers.computation.free.fr/Constants/constants.html

PDF: http://numbers.computation.free.fr/Constants/Primes/twin.pdf

## 1 Response to “Brun’s Constant and the Pentium Bug”