Basel problem


The Basel problem is a problem in mathematical analysis with relevance to number theory, concerning an infinite sum of inverse squares. It was first posed by Pietro Mengoli in 1650 and solved by Leonhard Euler in 1734,[1] and read on 5 December 1735 in The Saint Petersburg Academy of Sciences.[2] Since the problem had withstood the attacks of the leading mathematicians of the day, Euler's solution brought him immediate fame when he was twenty-eight. Euler generalised the problem considerably, and his ideas were taken up more than a century later by Bernhard Riemann in his seminal 1859 paper "On the Number of Primes Less Than a Given Magnitude", in which he defined his zeta function and proved its basic properties. The problem is named after Basel, hometown of Euler as well as of the Bernoulli family who unsuccessfully attacked the problem.

The Basel problem asks for the precise summation of the reciprocals of the squares of the natural numbers, i.e. the precise sum of the infinite series:

The sum of the series is approximately equal to 1.644934.[3] The Basel problem asks for the exact sum of this series (in closed form), as well as a proof that this sum is correct. Euler found the exact sum to be and announced this discovery in 1735. His arguments were based on manipulations that were not justified at the time, although he was later proven correct. He produced an accepted proof in 1741.

The solution to this problem can be used to estimate the probability that two large random numbers are coprime. Two random integers in the range from 1 to , in the limit as goes to infinity, are relatively prime with a probability that approaches , the reciprocal of the solution to the Basel problem.[4]

Euler's approach edit

Euler's original derivation of the value   essentially extended observations about finite polynomials and assumed that these same properties hold true for infinite series.

Of course, Euler's original reasoning requires justification (100 years later, Karl Weierstrass proved that Euler's representation of the sine function as an infinite product is valid, by the Weierstrass factorization theorem), but even without justification, by simply obtaining the correct value, he was able to verify it numerically against partial sums of the series. The agreement he observed gave him sufficient confidence to announce his result to the mathematical community.

To follow Euler's argument, recall the Taylor series expansion of the sine function

Dividing through by   gives

The Weierstrass factorization theorem shows that the left-hand side is the product of linear factors given by its roots, just as for finite polynomials. Euler assumed this as a heuristic for expanding an infinite degree polynomial in terms of its roots, but in fact it is not always true for general  .[5] This factorization expands the equation into:


If we formally multiply out this product and collect all the x2 terms (we are allowed to do so because of Newton's identities), we see by induction that the x2 coefficient of sin x/x is [6]


But from the original infinite series expansion of sin x/x, the coefficient of x2 is 1/3! = −1/6. These two coefficients must be equal; thus,


Multiplying both sides of this equation by −π2 gives the sum of the reciprocals of the positive square integers.


This method of calculating   is detailed in expository fashion most notably in Havil's Gamma book which details many zeta function and logarithm-related series and integrals, as well as a historical perspective, related to the Euler gamma constant.[7]

Generalizations of Euler's method using elementary symmetric polynomials edit

Using formulae obtained from elementary symmetric polynomials,[8] this same approach can be used to enumerate formulae for the even-indexed even zeta constants which have the following known formula expanded by the Bernoulli numbers:


For example, let the partial product for   expanded as above be defined by  . Then using known formulas for elementary symmetric polynomials (a.k.a., Newton's formulas expanded in terms of power sum identities), we can see (for example) that


and so on for subsequent coefficients of  . There are other forms of Newton's identities expressing the (finite) power sums   in terms of the elementary symmetric polynomials,   but we can go a more direct route to expressing non-recursive formulas for   using the method of elementary symmetric polynomials. Namely, we have a recurrence relation between the elementary symmetric polynomials and the power sum polynomials given as on this page by


which in our situation equates to the limiting recurrence relation (or generating function convolution, or product) expanded as


Then by differentiation and rearrangement of the terms in the previous equation, we obtain that


Consequences of Euler's proof edit

By the above results, we can conclude that   is always a rational multiple of  . In particular, since   and integer powers of it are transcendental, we can conclude at this point that   is irrational, and more precisely, transcendental for all  . By contrast, the properties of the odd-indexed zeta constants, including Apéry's constant  , are almost completely unknown.

The Riemann zeta function edit

The Riemann zeta function ζ(s) is one of the most significant functions in mathematics because of its relationship to the distribution of the prime numbers. The zeta function is defined for any complex number s with real part greater than 1 by the following formula:


Taking s = 2, we see that ζ(2) is equal to the sum of the reciprocals of the squares of all positive integers:


Convergence can be proven by the integral test, or by the following inequality:


This gives us the upper bound 2, and because the infinite sum contains no negative terms, it must converge to a value strictly between 0 and 2. It can be shown that ζ(s) has a simple expression in terms of the Bernoulli numbers whenever s is a positive even integer. With s = 2n:[9]


A proof using Euler's formula and L'Hôpital's rule edit

The normalized sinc function   has a Weierstrass factorization representation as an infinite product:


The infinite product is analytic, so taking the natural logarithm of both sides and differentiating yields


(by uniform convergence, the interchange of the derivative and infinite series is permissible). After dividing the equation by   and regrouping one gets


We make a change of variables ( ):


Euler's formula can be used to deduce that

or using the corresponding hyperbolic function:



Now we take the limit as   approaches zero and use L'Hôpital's rule thrice. By Tannery's theorem applied to  , we can interchange the limit and infinite series so that   and by L'Hôpital's rule


A proof using Fourier series edit

Use Parseval's identity (applied to the function f(x) = x) to obtain


for n ≠ 0, and c0 = 0. Thus,





as required.

Another proof using Parseval's identity edit

Given a complete orthonormal basis in the space   of L2 periodic functions over   (i.e., the subspace of square-integrable functions which are also periodic), denoted by  , Parseval's identity tells us that


where   is defined in terms of the inner product on this Hilbert space given by


We can consider the orthonormal basis on this space defined by   such that  . Then if we take  , we can compute both that


by elementary calculus and integration by parts, respectively. Finally, by Parseval's identity stated in the form above, we obtain that


Generalizations and recurrence relations edit

Note that by considering higher-order powers of   we can use integration by parts to extend this method to enumerating formulas for   when  . In particular, suppose we let


so that integration by parts yields the recurrence relation that


Then by applying Parseval's identity as we did for the first case above along with the linearity of the inner product yields that


Proof using differentiation under the integral sign edit

It's possible to prove the result using elementary calculus by applying the differentiation under the integral sign technique to an integral due to Freitas:[10]


While the primitive function of the integrand cannot be expressed in terms of elementary functions, by differentiating with respect to   we arrive at

which can be integrated by substituting   and decomposing into partial fractions. In the range   the definite integral reduces to


The expression can be simplified using the arctangent addition formula and integrated with respect to   by means of trigonometric substitution, resulting in


The integration constant   can be determined by noticing that two distinct values of   are related by

because when calculating   we can factor   and express it in terms of   using the logarithm of a power identity and the substitution  . This makes it possible to determine  , and it follows that


This final integral can be evaluated by expanding the natural logarithm into its Taylor series:


The last two identities imply


Cauchy's proof edit

While most proofs use results from advanced mathematics, such as Fourier analysis, complex analysis, and multivariable calculus, the following does not even require single-variable calculus (until a single limit is taken at the end).

For a proof using the residue theorem, see here.

History of this proof edit

The proof goes back to Augustin Louis Cauchy (Cours d'Analyse, 1821, Note VIII). In 1954, this proof appeared in the book of Akiva and Isaak Yaglom "Nonelementary Problems in an Elementary Exposition". Later, in 1982, it appeared in the journal Eureka,[11] attributed to John Scholes, but Scholes claims he learned the proof from Peter Swinnerton-Dyer, and in any case he maintains the proof was "common knowledge at Cambridge in the late 1960s".[12]

The proof edit

The inequality
is shown pictorially for any  . The three terms are the areas of the triangle OAC, circle section OAB, and the triangle OAB. Taking reciprocals and squaring gives

The main idea behind the proof is to bound the partial (finite) sums

between two expressions, each of which will tend to π2/6 as m approaches infinity. The two expressions are derived from identities involving the cotangent and cosecant functions. These identities are in turn derived from de Moivre's formula, and we now turn to establishing these identities.

Let x be a real number with 0 < x < π/2, and let n be a positive odd integer. Then from de Moivre's formula and the definition of the cotangent function, we have


From the binomial theorem, we have


Combining the two equations and equating imaginary parts gives the identity


We take this identity, fix a positive integer m, set n = 2m + 1, and consider xr = rπ/2m + 1 for r = 1, 2, ..., m. Then nxr is a multiple of π and therefore sin(nxr) = 0. So,


for every r = 1, 2, ..., m. The values xr = x1, x2, ..., xm are distinct numbers in the interval 0 < xr < π/2. Since the function cot2 x is one-to-one on this interval, the numbers tr = cot2 xr are distinct for r = 1, 2, ..., m. By the above equation, these m numbers are the roots of the mth degree polynomial


By Vieta's formulas we can calculate the sum of the roots directly by examining the first two coefficients of the polynomial, and this comparison shows that


Substituting the identity csc2 x = cot2 x + 1, we have


Now consider the inequality cot2 x < 1/x2 < csc2 x (illustrated geometrically above). If we add up all these inequalities for each of the numbers xr = rπ/2m + 1, and if we use the two identities above, we get


Multiplying through by (π/2m + 1)2
, this becomes


As m approaches infinity, the left and right hand expressions each approach π2/6, so by the squeeze theorem,


and this completes the proof.

Proof assuming Weil's conjecture on Tamagawa numbers edit

A proof is also possible assuming Weil's conjecture on Tamagawa numbers.[13] The conjecture asserts for the case of the algebraic group SL2(R) that the Tamagawa number of the group is one. That is, the quotient of the special linear group over the rational adeles by the special linear group of the rationals (a compact set, because   is a lattice in the adeles) has Tamagawa measure 1:


To determine a Tamagawa measure, the group   consists of matrices

with  . An invariant volume form on the group is

The measure of the quotient is the product of the measures of   corresponding to the infinite place, and the measures of   in each finite place, where   is the p-adic integers.

For the local factors,

where   is the field with   elements, and   is the congruence subgroup modulo  . Since each of the coordinates   map the latter group onto   and  , the measure of   is  , where   is the normalized Haar measure on  . Also, a standard computation shows that  . Putting these together gives  .

At the infinite place, an integral computation over the fundamental domain of   shows that  , and therefore the Weil conjecture finally gives

On the right-hand side, we recognize the Euler product for  , and so this gives the solution to the Basel problem.

This approach shows the connection between (hyperbolic) geometry and arithmetic, and can be inverted to give a proof of the Weil conjecture for the special case of  , contingent on an independent proof that  .

Other identities edit

See the special cases of the identities for the Riemann zeta function when   Other notably special identities and representations of this constant appear in the sections below.

Series representations edit

The following are series representations of the constant:[14]


There are also BBP-type series expansions for ζ(2).[14]

Integral representations edit

The following are integral representations of  [15][16][17]


Continued fractions edit

In van der Poorten's classic article chronicling Apéry's proof of the irrationality of  ,[18] the author notes as "a red herring" the similarity of a continued fraction for Apery's constant, and the following one for the Basel constant:

where  . Another continued fraction of a similar form is:[19]
where  .

See also edit

References edit

  • Weil, André (1983), Number Theory: An Approach Through History, Springer-Verlag, ISBN 0-8176-3141-0.
  • Dunham, William (1999), Euler: The Master of Us All, Mathematical Association of America, ISBN 0-88385-328-0.
  • Derbyshire, John (2003), Prime Obsession: Bernhard Riemann and the Greatest Unsolved Problem in Mathematics, Joseph Henry Press, ISBN 0-309-08549-7.
  • Edwards, Harold M. (2001), Riemann's Zeta Function, Dover, ISBN 0-486-41740-9.

Notes edit

  1. ^ Ayoub, Raymond (1974), "Euler and the zeta function", Amer. Math. Monthly, 81 (10): 1067–86, doi:10.2307/2319041, JSTOR 2319041
  2. ^ E41 – De summis serierum reciprocarum
  3. ^ Sloane, N. J. A. (ed.), "Sequence A013661", The On-Line Encyclopedia of Integer Sequences, OEIS Foundation
  4. ^ Vandervelde, Sam (2009), "Chapter 9: Sneaky segments", Circle in a Box, MSRI Mathematical Circles Library, Mathematical Sciences Research Institute and American Mathematical Society, pp. 101–106
  5. ^ A priori, since the left-hand-side is a polynomial (of infinite degree) we can write it as a product of its roots as
    Then since we know from elementary calculus that  , we conclude that the leading constant must satisfy  .
  6. ^ In particular, letting   denote a generalized second-order harmonic number, we can easily prove by induction that   as  .
  7. ^ Havil, J. (2003), Gamma: Exploring Euler's Constant, Princeton, New Jersey: Princeton University Press, pp. 37–42 (Chapter 4), ISBN 0-691-09983-9
  8. ^ Cf., the formulae for generalized Stirling numbers proved in: Schmidt, M. D. (2018), "Combinatorial Identities for Generalized Stirling Numbers Expanding f-Factorial Functions and the f-Harmonic Numbers", J. Integer Seq., 21 (Article 18.2.7)
  9. ^ Arakawa, Tsuneo; Ibukiyama, Tomoyoshi; Kaneko, Masanobu (2014), Bernoulli Numbers and Zeta Functions, Springer, p. 61, ISBN 978-4-431-54919-2
  10. ^ Freitas, F. L. (2023), "Solution of the Basel problem using the Feynman integral trick", arXiv:2312.04608 [math.CA]
  11. ^ Ransford, T J (Summer 1982), "An Elementary Proof of  " (PDF), Eureka, 42 (1): 3–4
  12. ^ Aigner, Martin; Ziegler, Günter M. (2001), Proofs from THE BOOK (2nd ed.), Springer, p. 32, ISBN 9783662043158; this anecdote is missing from later editions of this book, which replace it with earlier history of the same proof.
  13. ^ Vladimir Platonov; Andrei Rapinchuk (1994), Algebraic groups and number theory, translated by Rachel Rowen, Academic Press|
  14. ^ a b Weisstein, Eric W., "Riemann Zeta Function \zeta(2)", MathWorld
  15. ^ Connon, D. F. (2007), "Some series and integrals involving the Riemann zeta function, binomial coefficients and the harmonic numbers (Volume I)", arXiv:0710.4022 [math.HO]
  16. ^ Weisstein, Eric W., "Double Integral", MathWorld
  17. ^ Weisstein, Eric W., "Hadjicostas's Formula", MathWorld
  18. ^ van der Poorten, Alfred (1979), "A proof that Euler missed ... Apéry's proof of the irrationality of ζ(3)" (PDF), The Mathematical Intelligencer, 1 (4): 195–203, doi:10.1007/BF03028234, S2CID 121589323, archived from the original (PDF) on 2011-07-06
  19. ^ Berndt, Bruce C. (1989), Ramanujan's Notebooks: Part II, Springer-Verlag, p. 150, ISBN 978-0-387-96794-3

External links edit

  • An infinite series of surprises by C. J. Sangwin
  • From ζ(2) to Π. The Proof. step-by-step proof
  • Remarques sur un beau rapport entre les series des puissances tant directes que reciproques (PDF), English translation with notes of Euler's paper by Lucas Willis and Thomas J. Osler
  • Ed Sandifer, How Euler did it (PDF)
  • James A. Sellers (February 5, 2002), Beyond Mere Convergence (PDF), retrieved 2004-02-27
  • Robin Chapman, Evaluating ζ(2) (fourteen proofs)
  • Visualization of Euler's factorization of the sine function
  • Johan W Ästlund (December 8, 2010), Summing inverse squares by Euclidean geometry (PDF)
    • Why is pi here? And why is it squared? A geometric answer to the Basel problem on YouTube (animated proof based on the above)