Arxiv:2102.10542V1 [Math.GM] 21 Feb 2021

February 21, 2021

Solution of the Basel problem in the framework of distribution theory

Andreas Aste

Department of Physics, University of Basel, Klingelbergstrasse 82, CH-4056 Basel, Switzerland E-Mail: [email protected]

Abstract A simple proof of Euler’s formula which states that the sum of the reciprocals of all natural numbers squared equals π2/6 is presented based on the distribution theory introduced by Laurent Schwartz. Additional identities are obtained as a byproduct of the derivation. arXiv:2102.10542v1 [math.GM] 21 Feb 2021

Mathematics Subject Classiﬁcation MSC 2010: 40A25, 46F05 Keywords: Basel problem; Zeta function; distribution theory; generalized functions; test functions; summation of series 1 Introduction

P∞ 1 2 The so-called Basel problem to determine the sum ζ(2) = n=1 n2 = π /6 was ﬁrst posed in 1644 by Pietro Mengoli, an Italian mathematician and clergyman from Bologna, and solved by the Swiss mathematician Leonhard Euler (*1707 in Basel, †1783 in Saint Petersburg) in 1735. Several ways have been found in the meantime to calculate ζ(2) (see [2] and references therein). A further simple method to derive Euler’s result using the the theory of distributions and test functions which is based on elementary arguments like translational invariance is presented in this letter.

Distribution theory [1], which represents a mathematical discipline in its own right, is of fundamental signiﬁcance for a rigourous treatment of quantum ﬁeld theories in classical spacetime [3, 4]. It is also hoped that the stunning exercise presented in this letter serves as an incentive for graduate students with some basic knowledge of distribution theory to study the subject of generalized functions and their applications in theoretical physics in greater detail.

2 Calculating ζ(2)

0 We consider the distribution ∆0 ∈ D (R) deﬁned by the formal expression

∞ X inx −3ix −2ix −x ix 2ix 3ix ∆0(x) = e = ... + e + e + e + 1 + e + e + e + ..., (1) n=−∞ which acts on (smooth) test functions (with compact support) ϕ ∈ D(R) as a linear and, in the sense of distributions, continuous functional according to

∞ ∞ ∞ Z N Z X inx X inx ∆0[ϕ] := e ϕ(x)dx = lim e ϕ(x)dx . (2) N→∞ n=−∞−∞ n=−N−∞

0 In fact, ∆0 is well-defined by equation (2) as a distribution in D (R), the dual space of D(R), and equation (2) highlights the meaning of the formal definition (1) of ∆0 as a generalized function [6]. Note that a more intuitive representation of ∆0 as an alternative infinite sum of Dirac delta distributions is motivated in the appendix.

By deﬁnition, ∆0 is a periodic distribution invariant under a translation T2π, i.e. formally

∞ ∞ X in(x+2π) X inx (T2π∆0)(x) = ∆0(x + 2π) = e = e = ∆0(x) , (3) n=−∞ n=−∞ or in distributional notation

(T2π∆0)[ϕ] = ∆0[T2πϕ] ∀ϕ ∈ D(R) , where (T2πϕ)(x) = ϕ(x − 2π) , (4) and ∆0 is symmetric ∆0(x) = ∆0(−x) . (5) ix Now since ∆0(x) is invariant with respect to a multiplication with e , i.e.

∞ ∞ ix ix X inx X inx e ∆0(x) = e e = e , (6) n=−∞ n=−∞

1 ∆0 must vanish as a distribution on R\{2πn | n ∈ Z}, since only for x = 2πn with n ∈ Z one ix has a trivial factor e = 1; therefore the distributional support of ∆0 must be contained in a corresponding discrete set supp ∆0 ⊆ {2πn | n ∈ Z} . (7) For a moment, the following considerations are restricted to the open interval I = (0, 2π). Calculating the ﬁrst antisymmetric antiderivative ∆1 of ∆0 with x ∈ I x Z ∞ inx ∞ 0 0 X e X sin(nx) ∆1(x) = lim ∆0(x )dx = −i + x = 2 + x (8) &0 n n n=−∞ n=1 n6=0 with ∆1(x) = −∆1(−x) , (9)

∆1 must be constant on I, since its derivative ∆0 vanishes there. This also implies that the Fourier sum in equation (8) represents a linear function on I. Calculating the mean value µI,1 of ∆1 on I according to 2π− 1 Z µI,1 = lim ∆1(x)dx , (10) 2π &0 inx the oscillatory terms ∼ e in equation (8) do not contribute to µI,1 and one is left with 2π 1 Z µ = xdx = π . (11) I,1 2π 0

Finally turning to the antiderivative of ∆1 on I

x ∞ Z X einx 1 X cos(nx) 1 ∆ (x) = ∆ (x0)dx0 = − + x2 = −2 + x2 , (12) 2 1 n2 2 n2 2 Z 0 n∈ \0 n=1 one arrives at an expression containing a series that converges absolutely to a continuous function on I. However, since the distributional derivative of ∆2 is ∆1 which is constant, i.e., π on I, ∆2 must be of the form

∆2(x) = πx + γ , x ∈ I (13) with an integration constant γ. This constant can be calculated by considering the average value of ∆2 on I 2π 2π 1 Z dx 2π2 1 Z µ = x2 = = (πx + γ)dx = π2 + γ , (14) I,2 2π 2 3 2π 0 0 hence γ = −π2/3, an ﬁnally Euler’s famous result ∞ ∞ X 1 π2 X 1 π2 ∆ (0) = −2 = γ = − → = (15) 2 n2 3 n2 6 n=1 n=1 follows.

As an exercise, the reader may verify that by considering additional antiderivatives of ∆2 like ∆4, ∆6 et cetera, further values of the Euler-Riemann zeta function like ∞ ∞ X 1 π4 X 1 π6 ζ(4) = = , ζ(6) = = ,... (16) n4 90 n6 945 n=1 n=1 follow directly from strategy outlined above.

2 A Explicit representation of ∆0 as an inﬁnite sum of Dirac delta distributions

0 R We consider the following sequence {δN }N∈N0 ⊂ D ( ) of distributions [6] represented by the functions  N N  P inx 2 2 X inx  e |x| < π δN (x) = Θ(π − x ) e = n=−N (17) n=−N 0 |x| ≥ π with supp δN = [−π, π], where ( 1 x > 0 Θ(x) = (18) 0 x ≤ 1

−iNx −i(N−1)x −ix ix iNx is the Heaviside function. With δN (x) = e + e + ... + e + 1 + e + ... + e and ix i(N+1)x −iNx e δN (x) = δN (x) + e − e (19) for x ∈ (−π, π) one immediately obtains the compact representation

ei(N+1)x − e−iNx ei(N+1/2)x − e−i(N+1/2)x sin((N + 1/2)x) δ (x) = = = , x ∈ (−π, π)\{0} N eix − 1 eix/2 − e−ix/2 sin(x/2) (20) and from the deﬁnition (17) one has δN (0) = 2N + 1 which removes the singularity appearing at x = 0 in the representation (20). Only the term ei0x = 1 in deﬁnition (17) contributes to the integral ∞ π Z Z δN (x)dx = δN (x)dx = 2π . (21) −∞ −π 0 R For illustrative purposes, the graph of δ50 is depicted in Fig. 1. In fact, {δN }N∈N0 ⊂ D ( ) is a δ-sequence converging to 2π times the Dirac delta distribution δ for N → ∞. Applying δN on a (smooth) test function ϕ ∈ D(R) (with compact support) leads to

∞ π Z Z sin((N + 1/2)x) δ [ϕ] = δ (x)ϕ(x)dx = ϕ(x)dx N N sin(x/2) −∞ −π

π π Z sin((N + 1/2)x) x/2 Z sin((N + 1/2)x) x/2 = 2 ϕ(x)dx = 2 β(x) ϕ(x)dx (22) x sin(x/2) x sin(x/2) −π −π where a smooth bump function β ∈ D(R) with the properties β(x) = 1 for |x| ≤ π and β(x) = 0 for |x| ≥ 3π/2 was introduced which does not change the integral above. Since one has ( 1 = x/2 |x| ∈ (0, 3π/2] σ(x) = si(x/2) sin(x/2) , σ ∈ C∞([−3π/2, 3π/2]) , (23) 1 x = 0 i.e. since σ is a smooth function on the interval [−3π/2, 3π/2], alsoϕ ˜(x) = β(x)σ(x)ϕ(x) is smooth and has compact support:ϕ ˜ ∈ D(R). Furthermore,ϕ ˜(0) = ϕ(0) holds.

Now, equation (22) becomes, with x0 = (N + 1/2)x in the limit N → ∞ in the sense of distributions

∞ π (N+1/2)π Z Z sin((N + 1/2)x) Z sin(x0) δ [ϕ] = δ (x)ϕ(x)dx = 2 ϕ˜(x)dx = 2 ϕ˜(x0/(N+1/2))dx0 N N x x0 −∞ −π −(N+1/2)π

3 100

(x) 40 50

-20

-4-3-2-101234 x

Figure 1: The graph of δ50 deﬁned by equation (17).

∞ Z sin(x0) −−−−→N→∞ 2 ϕ˜(0)dx0 = 2πϕ˜(0) = 2πϕ(0) = 2πδ[ϕ] . (24) x0 −∞ The normalization of the δ-distribution follows from equation (21), i.e., as a byproduct of the derivation presented above the integral

(N+1/2)π ∞ Z sin(x) Z sin(x) lim dx = dx = π (25) N→∞ x x −(N+1/2)π −∞ is obtained. Neglecting the cutoﬀ in deﬁnition (17) leads to the periodic distributional identity

∞ ∞ ∞ X inx X X ∆0(x) = e = 2π δ(x − 2πn) or ∆0[ϕ] = 2π ϕ(2πn) . (26) n=−∞ n=−∞ n=−∞

One readily expresses the antisymmetric antiderivative of ∆0 by the help of the ﬂoor function b·c and the ceiling function d·e

x x ∆ (x) = π + , (27) 1 2π 2π which simpliﬁes to ∆1(x) = π sign(x) (28) on the open interval (−2π, 2π), and the symmetric antiderivative of ∆1 is represented by the continuous function x x x x π2 X einx x2 ∆ (x) = πx + −2π2 − = − + . (29) 2 2π 2π 2π 2π 3 n2 2 n∈Z\{0}

4 References

[1] Schwartz, L.: Généralisation de la notion de fonction, de dérivation,de transformation de Fourier et applications mathématiques et physiques. Ann. Univ. Grenoble. Sect. Sci. Math. Phys. (N.S.) 21 (1945), pp. 57–74 (1945).

[2] Riemenschneider, O.: Uber¨ einige elementare analytische Berechnungen von ζ(2). Vari- ationen ¨uber ein Thema von Leonhard Euler. Mitt. Math. Ges. Hamburg XXXXVI, pp. 53-69 (2016).

[3] Streater, R.F., Wightman, A.S.: PCT, Spin, Statistics and All That. Benjamin- Cummings Publishing Company, 1964.

[4] Scharf, G.: Finite Quantum Electrodynamics: The Causal Approach. Dover Books on Physics, 2014.

[5] Epstein H., Glaser V.: The role of locality in perturbation theory. Annales Poincar´ePhys. Theor. A19, pp. 211-295 (1973).

[6] Constantinescu, F.: Distributions and Their Applications in Physics. Pergamon Press, 1980.