Divergence of the sum of the reciprocals of the primes

The sum of the reciprocals of all prime numbers diverges, that is:

\sum_{p\text{ prime }}\frac1p = \frac12 + \frac13 + \frac15 + \frac17 + \frac1{11} + \frac1{13} + \frac1{17} +\cdots = \infty.

This was proved by Leonhard Euler in 1737, and strengthens Euclid's 3rd-century-BC result that there are infinitely many prime numbers.

There is a variety of proofs of Euler's result, including a lower bound for the partial sums stating that

\sum_{\scriptstyle p\text{ prime }\atop \scriptstyle p\le n}\frac1p \ge \ln \ln (n+1) - \ln\frac{\pi^2}6

for all natural numbers n. The iterated natural logarithm indicates that the divergence might be very slow, which is indeed the case, see Meissel–Mertens constant.

Contents

The harmonic series

First, we describe how Euler originally discovered the result. He was considering the harmonic series

\sum_{n=1}^\infty \frac{1}{n} = 
1 + \frac{1}{2} + \frac{1}{3} + \frac{1}{4} +
\cdots.

He had already used the following "product formula" to show the existence of infinitely many primes.

\sum_{n=1}^\infty \frac{1}{n} = \prod_{p} \frac{1}{1-p^{-1}}
=\prod_{p} \left( 1+\frac{1}{p}+\frac{1}{p^2}+\cdots \right).

(Here, the product is taken over all primes p; in the following, a sum or product taken over p always represents a sum or product taken over a specified set of primes, unless noted otherwise.)

Such infinite products are today called Euler products. The product above is a reflection of the fundamental theorem of arithmetic. (Multiply out the right side as you would like to do.) Of course, the above "equation" is not necessary because the harmonic series is known (by other means) to diverge. This type of formal manipulation was common at the time, when mathematicians were still experimenting with the new tools of calculus.[citation needed]

Euler noted that if there were only a finite number of primes, then the product on the right would clearly converge, contradicting the divergence of the harmonic series. (In modern language, we now say that the existence of infinitely many primes is reflected by the fact that the Riemann zeta function has a simple pole at s = 1.)

First proof

Euler took the above product formula and proceeded to make a sequence of audacious leaps of logic. First, he took the natural logarithm of each side, then he used the Taylor series expansion for ln(1 − x) as well as the sum of a geometric series:


\begin{align}
\ln \left( \sum_{n=1}^\infty \frac{1}{n}\right) & {} = \ln \left( \prod_p \frac{1}{1-p^{-1}}\right)
= \sum_p \ln \left( \frac{1}{1-p^{-1}}\right) = \sum_p - \ln(1-p^{-1}) \\
& {} = \sum_p \left( \frac{1}{p} + \frac{1}{2p^2} + \frac{1}{3p^3} + \cdots \right) \\
& {} = \left( \sum_{p}\frac{1}{p} \right) + \sum_p \frac{1}{p^2} \left( \frac{1}{2} + \frac{1}{3p} + \frac{1}{4p^2} + \cdots \right) \\
& {} < \left( \sum_p \frac{1}{p} \right) + \sum_p \frac{1}{p^2} \left( 1 + \frac{1}{p} + \frac{1}{p^2} + \cdots \right) \\
& {} = \left( \sum_p \frac{1}{p} \right) + \left( \sum_p \frac{1}{p(p-1)} \right) \\
& {} = \left( \sum_p \frac{1}{p} \right) + C
\end{align}

for a fixed constant C < 1. Since the sum of the reciprocals of the first n positive integers is asymptotic to ln(n), (i.e. their ratio approaches one as n approaches infinity), Euler then concluded

\frac{1}{2} + \frac{1}{3} + \frac{1}{5} + \frac{1}{7} + \frac{1}{11} + \cdots = \ln \ln (+ \infty).

It is almost certain that Euler meant that the sum of the reciprocals of the primes less than n is asymptotic to ln(ln(n)) as n approaches infinity. It turns out this is indeed the case; Euler had reached a correct result by questionable means.

Simpler version of the above proof


\begin{align}
& {} \quad \ln \left( \sum_{n=1}^\infty \frac{1}{n}\right) = \ln \left( \prod_p \frac{1}{1-p^{-1}}\right) = \sum_p \ln \left( \frac{p}{p-1}\right) = \sum_p \ln\left(1+\frac{1}{p-1}\right)
\end{align}

Since

 e^x = 1 + x + \frac{x^2}{2!} + \frac{x^3}{3!} + \cdots,

we get ex > 1 + x and x > ln(1 + x).

So

 \sum_p \ln\left(1+\frac{1}{p-1}\right) < \sum_p \frac{1}{p - 1}

Hence \sum_p \frac{1}{p-1} diverges. But 1/(pi − 1) < 1/pi−1 where pi is the ith prime. Hence \sum_p \frac{1}{p} diverges.

Second proof

The following proof by contradiction is due to Paul Erdős.

Let pi denote the ith prime number. Assume that the sum of the reciprocals of the primes converges, i.e.

\sum_{i=1}^\infty {1\over p_i} < \infty\,.

Then there exists a positive integer k such that

\sum_{i=k+1}^\infty {1\over p_i} < {1\over 2}.\qquad(1)

For a positive integer x let Mx denote the set of those n in {1, 2, . . ., x} which are not divisible by any prime greater than pk. We will now derive an upper and a lower estimate for the number |Mx| of elements in Mx. For large x, these bounds will turn out to be contradictory.

Upper estimate: Every n in Mx can be written as n = r m2 with positive integers m and r, where r is square-free. Since only the k primes p1 . . ., pk can show up (with exponent 1) in the prime factorization of r, there are at most 2k different possibilities for r. Furthermore, there are at most √x possible values for m. This gives us the upper estimate

|M_x| \le  2^k\sqrt{x}\,.\qquad(2)

Lower estimate: The remaining x − |Mx| numbers in the set difference {1, 2, . . ., x} \ Mx are all divisible by a prime greater than pk. Let Ni,x denote the set of those n in {1, 2, . . ., x} which are divisible by the ith prime pi. Then

\{1,2,\ldots,x\}\setminus M_x\subset\bigcup_{i=k+1}^\infty N_{i,x}\;.

Since the number of integers in Ni,x is at most x/pi (actually zero for pi > x), we get

x-|M_x| \le \sum_{i=k+1}^\infty |N_{i,x}|< \sum_{i=k+1}^\infty {x\over p_i}\,.

Using (1), this implies

{x\over 2} < |M_x|\,.\qquad(3)

Contradiction: For every integer x ≥ 22k + 2, the estimates (2) and (3) cannot hold simultaneously.

Third proof

Here is another proof that actually gives a lower estimate for the partial sums; in particular, it shows that these sums grow at least as fast as ln(ln(n)). The proof is an adaptation of the product expansion idea of Euler. In the following, a sum or product taken over p always represents a sum or product taken over a specified set of primes.

The proof rests upon the following four inequalities:

  • Every positive integer i can be uniquely expressed as the product of a square-free integer and a square. This gives the inequality
 \sum_{i=1}^n{\frac{1}{i}} \le \prod_{p \le n}{\biggl(1+\frac{1}{p}\biggr)}\sum_{k=1}^n{\frac{1}{k^2}}\,,
where for every i between 1 and n the (expanded) product contains to the square-free part of i and the sum contains to the square part of i (see fundamental theorem of arithmetic).

\ln(n+1)
= \int_1^{n+1}\frac{dx}x
=\sum_{i=1}^n\underbrace{\int_i^{i+1}\frac{dx}x}_{{}<1/i}
<\sum_{i=1}^n{\frac{1}{i}}.
  • Let n ≥ 2. The upper bound (using a telescoping sum) for the partial sums (convergence is all we really need)

 \sum_{k=1}^n{\frac{1}{k^2}}
< 1+\sum_{k=2}^n\underbrace{\biggl(\frac1{k-1/2}-\frac1{k+1/2}\biggr)}_{= 1/(k^2-1/4) > 1/k^2}=1+\frac23-\frac1{n+1/2} < \frac53.

Combining all these inequalities, we see that


\ln(n+1) 
< \sum_{i=1}^n\frac{1}{i}
\le \prod_{p \le n}{\biggl(1+\frac{1}{p}\biggr)}\sum_{k=1}^n{\frac{1}{k^2}}
< \frac53\prod_{p \le n}{\exp\biggl(\frac{1}{p}\biggr)}
= \frac53\exp\biggl(\sum_{p \le n}{\frac{1}{p}}\biggr).

Dividing through by 5/3 and taking the natural logarithm of both sides gives

 \ln \ln(n+1) - \ln\frac53 < \sum_{p \le n}{\frac{1}{p}}

as desired. 

Using


\sum_{k=1}^\infty{\frac{1}{k^2}}=\frac{\pi^2}6

(see Basel problem), the above constant ln (5/3) = 0.51082... can be improved to ln(π2/6) = 0.4977...; in fact it turns out that

\lim_{n \to \infty } \biggl( 
\sum_{p \leq n} \frac{1}{p}  - \ln(\ln(n)) \biggr)= M,

where M = 0.261497... is the Meissel–Mertens constant (somewhat analogous to the much more famous Euler–Mascheroni constant).

Fourth proof

From Dusart's inequality (see PNT), we get

 p_n <  n \ln n + n \ln \ln n \quad\mbox{for } n \ge 6.

Then


\sum_{n=1}^\infty \frac1{ p_n}
\ge \sum_{n=6}^\infty \frac1{ p_n}
\ge \sum_{n=6}^\infty \frac1{ n \ln n + n \ln \ln n}
\ge \sum_{n=6}^\infty \frac1{2n \ln n}
= \infty

by the integral test for convergence. This shows that the series on the left diverges.

See also

References

External links


Wikimedia Foundation. 2010.

Look at other dictionaries:

  • Proof that the sum of the reciprocals of the primes diverges — In the third century BC, Euclid proved the existence of infinitely many prime numbers. In the 18th century, Leonhard Euler proved a stronger statement: the sum of the reciprocals of all prime numbers diverges. Here, we present a number of proofs… …   Wikipedia

  • Contributions of Leonhard Euler to mathematics — The 18th century Swiss mathematician Leonhard Euler (1707–1783) is among the most prolific and successful mathematicians in the history of the field. His seminal work had a profound impact in numerous areas of mathematics and he is widely… …   Wikipedia

  • Harmonic series (mathematics) — In mathematics, the harmonic series is the divergent infinite series: Its name derives from the concept of overtones, or harmonics in music: the wavelengths of the overtones of a vibrating string are 1/2, 1/3, 1/4, etc., of the string s… …   Wikipedia

  • Integral test for convergence — In mathematics, the integral test for convergence is a method used to test infinite series of non negative terms for convergence. An early form of the test of convergence was developed in India by Madhava in the 14th century, and by his followers …   Wikipedia

  • Leonhard Euler — Infobox Scientist name = Leonhard Euler|box width = 300px |200px image width = 200px caption = Portrait by Johann Georg Brucker birth date = birth date|df=yes|1707|4|15 birth place = Basel, Switzerland death date = 18 September (O.S 7 September)… …   Wikipedia

  • List of mathematical proofs — A list of articles with mathematical proofs:Theorems of which articles are primarily devoted to proving them: See also: *Bertrand s postulate and a proof *Estimation of covariance matrices *Fermat s little theorem and some proofs *Gödel s… …   Wikipedia

  • Convergent series — redirects here. For the short story collection, see Convergent Series (short story collection). In mathematics, a series is the sum of the terms of a sequence of numbers. Given a sequence , the nth partial sum Sn is the sum of the first n terms… …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”

We are using cookies for the best presentation of our site. Continuing to use this site, you agree with this.