Next Article in Journal
The Spectrum of Teleparallel Gravity
Next Article in Special Issue
Dimension and Dimensional Reduction in Quantum Gravity
Previous Article in Journal
Explaining Defects of the Universal Vacua with Black Holes-Hedgehogs and Strings
Previous Article in Special Issue
Geometric Operators in the Einstein–Hilbert Truncation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Role of Riemann’s Zeta Function in Mathematics and Physics †,‡

Institut für Theoretische Physik, Universität Tübingen, Auf der Morgenstelle 14, D-72076 Tübingen, Germany
This article is dedicated to Martin Reuter on the occasion of his 60th birthday, as well as to a mutual hero of ours: Bernhard Riemann.
This paper is based on a talk given at the conference “Quantum Fields—From Fundamental Concepts to Phenomenological Questions”, Mainz, Germany, 26–28 September 2018.
Submission received: 28 January 2019 / Revised: 28 February 2019 / Accepted: 11 March 2019 / Published: 14 March 2019

Abstract

:
In particular, Riemann’s impact on mathematics and physics alike is demonstrated using methods originating from the theory of numbers and from quantum electrodynamics, i.e., from the behavior of an electron in a prescribed external electromagnetic field. More specifically, we employ Riemann’s zeta function to regularize the otherwise infinite results of the so-called Heisenberg–Euler Lagrangian. As a spin-off, we also calculate some integrals that are useful in mathematics and physics.

1. Usefulness of Riemann’s Functional Equation for the Zeta Function in Physics

Riemann introduced in his talk in Berlin in 1859 [1] (p. 147)
Γ s 2 π s 2 ζ ( s ) = 1 d x ψ ( x ) x s 2 1 + x s 2 1 2 1 s ( 1 s ) .
where ψ ( x ) is related to one of Jacobi’s θ functions. Notice that there is no change of the right-hand side under s ( 1 s ) . π s 2 Γ s 2 ζ ( s ) has simple poles at s = 0 (from Γ ) and s = 1 (from ζ ). To remove these poles, we multiply by 1 2 s ( s 1 ) . This is the reason why Riemann defines
ξ ( s ) = 1 2 s ( s 1 ) π s 2 Γ s 2 ζ ( s ) ,
which is an entire function ( ζ is a meromorphic function). Obviously, we have
ξ ( s ) = ξ ( 1 s )
together with the symmetrical form of the functional equation, which was proved by Riemann for all complex s:
Γ s 2 π s 2 ζ ( s ) = Γ 1 s 2 π 1 2 ( 1 s ) ζ ( 1 s ) .
Notice that the right-hand side is obtained from the left-hand side by replacing s by 1 s .
Before we continue with Equation (3), we make use of two important formulae due to Euler and Legendre.
Legendre’s duplication formula: Γ functions of argument 2 s can be expressed in terms of Γ functions of smaller arguments.
Γ ( 2 s ) = ( 2 π ) 1 2 2 2 s 1 2 Γ ( s ) Γ s + 1 2 = 1 π 2 2 s 1 Γ ( s ) Γ s + 1 2 .
When we replace s s 2 , we obtain:
Γ ( s ) = 1 π 2 s 1 Γ s 2 Γ s + 1 2 or π 2 s 1 Γ ( s ) = Γ s 2 Γ s + 1 2 .
Euler’s reflection formula for the Γ function:
Γ ( s ) Γ ( 1 s ) = π sin π s .
Replace s s + 1 2 = s 2 + 1 2 :
Γ s + 1 2 Γ 1 s + 1 2 = π sin π s 2 + π 2 or Γ s + 1 2 Γ 1 s 2 = π cos π s 2 .
Now multiply Equation (3) by Γ s + 1 2 to obtain
π s 2 Γ s 2 Γ s + 1 2 ζ ( s ) = π 1 s 2 Γ s + 1 2 Γ 1 s 2 ζ ( 1 s ) .
Then, we get from Equations (5) and (7) the following equations:
π s 2 π 2 s 1 Γ ( s ) ζ ( s ) = π 1 s 2 π cos π s 2 ζ ( 1 s ) , ζ ( 1 s ) = 2 ( 2 π ) s cos π s 2 Γ ( s ) ζ ( s ) .
Here, we replace s 1 s thereby obtaining
ζ ( s ) = 2 ( 2 π ) 1 s cos π ( 1 s ) 2 Γ ( 1 s ) ζ ( 1 s )
or ζ ( s ) = 2 s π s 1 sin π s 2 Γ ( 1 s ) ζ ( 1 s ) .
The latter equation is of great importance in the following sections.

2. Correction of the Classical Electromagnetic Lagrangian by Vacuum Electrons

In 1936, Werner Heisenberg and Hans Euler [2] wrote down the first effective Lagrangian in quantum field theory, which incorporates a quantum correction to the classical Lagrangian of a constant electromagnetic field; this correction is due to the polarization of the quantum vacuum (Dirac’s idea), i.e., the effect of an external constant electromagnetic field on the motion of the vacuum electrons. To simplify matters, we only consider a constant magnetic field in z direction. For this special case, the modified Lagrangian takes the form—in Schwinger’s representation:
L ( B ) = L ( 0 ) + L ( 1 ) , L ( 0 ) = 1 2 B 2 , L ( 1 ) ( B ) = 1 8 π 2 0 d s s 3 e i m 2 s ( e B s ) cot ( e B s ) + 1 3 ( e B s ) 2 1 .
The integral was explicitly calculated for the first time in [3] by dimensional regularization and thereafter in [4] by the so-called zeta-function regularization. The findings of these two different methods agree exactly, whereby the result obtained by the zeta-function regularization is finite without the usual subtraction of divergent counterterms. The result turns out to be
L ( 1 ) ( B ) = 1 32 π 2 3 m 4 + 4 ( e B ) 2 1 3 4 ζ ( 1 ) + 4 m 2 ( e B ) ( ln 2 π 1 ) 2 m 4 ln 2 e B m 2 4 m 2 ( e B ) ln 2 e B m 2 4 3 ( e B ) 2 ln 2 e B m 2 16 ( e B ) 2 1 1 + m 2 2 e B d x ln Γ ( x ) .
For those values of the field strength, i.e., for strong fields e B m 2 1 , the integral over the logarithm of the gamma function only yields a constant.
With b = e B m 2 , we obtain
b 2 1 1 + 1 2 b d x ln Γ ( x ) b 2 1 1 + 1 2 b d x ln Γ ( 1 ) + d d x ln Γ ( x ) | x = 1 ( x 1 ) = b 2 Ψ ( 1 ) 1 1 + 1 2 b d x ( x 1 ) = 1 8 Ψ ( 1 ) = 1 8 C .
Here, C is Euler’s number, C = ln γ = 0.57721566490 , and the digamma function
Ψ ( x ) = d d x ln Γ ( x ) = Γ ( x ) Γ ( x ) at x = 1 is   given   by Ψ ( 1 ) = ln γ .
Therefore, by only considering the dominant forms for large magnetic field strength, we obtain for the asymptotic form of the one-loop effective Lagrangian in spinor quantum electrodynamics (QED)
lim e B m 2 L ( 1 ) ( B ) = 1 32 π 2 4 ( e B ) 2 1 3 4 ζ ( 1 ) 4 3 ( e B ) 2 ln 2 e B m 2 = α B 2 6 π ln e B m 2 + 12 ζ ( 1 ) 1 + ln 2 , α = e 2 4 π .
On the other hand, we find in Ritus’ paper [5] under Formula (60) the expression
lim e B m 2 L ( 1 ) ( B ) = α B 2 6 π ln e B γ π m 2 + 6 π 2 ζ ( 2 ) .
Since Equations (15) and (16) are just two different representations of the same strong-field Lagrangian L ( 1 ) ( B ) , we have the equality
ln 2 e B m 2 1 + 12 ζ ( 1 ) = ln 2 e B m 2 2 π γ + 6 π 2 ζ ( 2 )
or
1 + 12 ζ ( 1 ) = ln ( 2 π γ ) + 6 π 2 ζ ( 2 ) .
This important equation is in fact a direct consequence of a variant of the famous functional equation of Riemann’s zeta function [1]
ζ ( s ) = 2 s π s 1 sin π s 2 Γ ( 1 s ) ζ ( 1 s ) .
Note that the result in Equation (17), which arises from long, complicated field theoretic calculations, follows from solving a physics problem—not from analytical theory of numbers, namely by studying the behavior of vacuum electrons in presence of a constant external strong classical magnetic field. (Later, we meet ζ ( 1 ) and ζ ( 2 ) again in connection with the Glaisher–Kinkelin constant.) Thus far, we consider nonlinear spinor QED where spin 1 2 particles with mass m are coupled to an external constant magnetic field. The corresponding effective Lagrangian is given by Equation (12).
Now, we study charged spinless particles with mass m, associated with a complex scalar field, which interact with a constant magnetic field. Here, the starting point is given by the Heisenberg–Euler effective Lagrangian, which in Schwinger’s proper time representation reads:
L s c a l a r ( 1 ) ( B ) = 1 16 π 2 0 1 s 3 e i m 2 s e B s sin ( e B s ) 1 6 ( e B s ) 2 1 d s .
Without going into detail, we just repeat the former results for spinor particles obtained for strong magnetic fields, when performing the calculation in both the ζ function and proper time regularization. Here are the results for scalar QED:
lim e B m 2 L s c a l a r ( 1 ) ( B ) = α B 2 24 π ln 2 e B m 2 + 12 ζ ( 1 ) 1 + ln 2
= α B 2 24 π ln 2 e B m 2 + ln γ π + 6 π 2 ζ ( 2 )
This brings us back to Equation (17), which is equivalent to Riemann’s functional equation for the zeta function. This functional equation is evidently independent of the masses involved, be they fermionic or scalar. Were it not for the factors 1 24 π instead of 1 6 π and ln 2 e B m 2 instead of ln e B m 2 in the spinor case, we could have guessed Equations (20) and (21). However, arriving from the proper-time integrals of Equation (12) or Equation (19) at Equations (15), (16), (20), and (21) is a highly challenging undertaking.
Finally, let us mention that the results of Equations (15) and (20) can be used in the Callan–Szymanzik renormalization equation to calculate the β ζ ( α ) function for spinor and scalar QED to result in
β ζ ( α ) = 2 3 α 2 π ( spinor ) , β ζ ( α ) = 1 6 α π ( scalar ) .
This confirms the correctness of the results in Equation (22) calculated otherwise.

3. Proof of Equation (17) from the Functional Equation of Riemann’s Zeta Function

Since Equation (17) contains derivatives of the ζ ( s ) function at s = 1 and s = 2 , we need the derivative of Equation (18). Writing 2 s π s 1 = e s ln 2 e ( s 1 ) ln π , we obtain
ζ ( s ) = ( ln 2 + ln π ) 2 s π s 1 sin π s 2 Γ ( 1 s ) ζ ( 1 s ) + π 2 2 s π s 1 cos π s 2 Γ ( 1 s ) ζ ( 1 s ) 2 s π s 1 sin π s 2 ( Γ ( 1 s ) ζ ( 1 s ) + Γ ( 1 s ) ζ ( 1 s ) ) .
For s = 2 , we find
ζ ( 2 ) = π 2 2 2 π 3 cos ( π ) Γ ( 3 ) ζ ( 3 ) , Γ ( 3 ) = 2 , ζ ( 2 ) = ζ ( 3 ) 4 π 2
or
ζ ( 3 ) = 4 π 2 ζ ( 2 ) ,
a result we use below.
For s = 1 , we employ ζ ( 2 ) = π 2 6 (Euler’s Basel problem) and Γ ( 2 ) = 1 C = 1 ln γ , Γ ( 2 ) = 1 , so that
ζ ( 1 ) = ln 2 × 1 12 ln π × 1 12 + ( 1 ln γ ) × 1 12 + 1 2 × 1 π 2 ζ ( 2 )
or
ζ ( 1 ) = 1 12 1 ln ( 2 π γ ) + 6 π 2 ζ ( 2 ) ,
that can also be written as
1 + 12 ζ ( 1 ) = ln ( 2 π γ ) + 6 π 2 ζ ( 2 ) ,
which is exactly the equation that followed from the two regularization methods, which produced the physically motivated results Equations (15) and (16).

4. Asymptotic Expansions of ln Γ ( x + 1 ) and ln Γ 1 ( x + 1 )

On the way to calculating the constant L 1 (for the definition cf. Equation (33)), we start with the asymptotic expansion of ln Γ ( x + 1 ) [6]:
ln ( x ! ) = x = 1 x ln x = L 0 + x + 1 2 ln x x + B 2 1 × 2 x + B 4 3 × 4 x 3 + B 6 5 × 6 x 5 + = ln Γ ( x + 1 ) .
This is called the Moivre–Stirling formula. Actually, it was discovered by Moivre (1667–1754) with the aid of Euler’s (1707–1783) summation formula. Stirling (1692–1770) “only” showed, using Wallis’ (1616–1703) product formula, that the constant L 0 is given by
L 0 = ln 2 π = 0.918938533 .
The Bernoulli numbers in Equation (26) are
B 1 = 1 2 , B 2 n + 1 = 0 , n = 1 , 2 , 3 , , B 2 = 1 6 , B 4 = 1 30 , B 6 = 1 42 , B 8 = 1 30 .
Useful integrals of ln Γ ( x + 1 ) are due to Raabe (1801–1859):
x x + 1 d t ln Γ ( t ) = x 1 x d t ln Γ ( t + 1 ) = x ln x x + ln 2 π .
In particular,
0 1 d x ln Γ ( x + 1 ) = ln 2 π 1 , 0 1 d x ln Γ ( x ) = ln 2 π .
The generalized Γ function Γ 1 ( x ) shows up in the value of the integral
0 x d t ln Γ ( t ) = ln Γ 1 ( x ) x 2 ( x 1 ) + x ln 2 π .
One could use Equation (30) as defining equation for ln Γ 1 ( x ) . Putting x = 1 2 , we obtain from Equation (30)
0 1 2 d x ln Γ ( x ) = ln Γ 1 ( 1 2 ) + 1 8 + 1 2 ln 2 π = 3 2 L 1 + 5 24 ln 2 + 1 4 ln π ,
where we use
ln Γ 1 ( 1 2 ) = 3 2 L 1 1 8 1 24 ln 2 ,
which can be looked up in Ref. [6]. The constant L 1 remains to be determined, which is the main goal of this paper.
Next, we introduce the asymptotic expansion of ln Γ 1 ( x + 1 ) [6]:
ln 1 1 × 2 2 × 3 3 x x = k = 1 x k ln k = : ln Γ 1 ( x + 1 ) = L 1 + x 2 ( x + 1 ) + 1 12 ln x 1 4 x 2 B 4 2 × 3 × 4 x 2 + B 6 4 × 5 × 6 x 4 + B 8 6 × 7 × 8 x 6 + .
As pointed out in the title of this paper, the constant L 1 is determined by comparing different regularization methods arising from physics arguments.
The generalized Γ function of the first kind satisfies the relations
ln Γ 1 ( x + 1 ) = 0 x d t ln Γ ( t + 1 ) + x 2 ( x + 1 ) x ln 2 π ,
ln Γ 1 ( x ) = 0 x d t ln Γ ( t ) + x 2 ( x 1 ) x ln 2 π .
Γ 1 ( x ) satisfies the functional equation
Γ 1 ( x + 1 ) = x x Γ 1 ( x )
with the constraints Γ 1 ( 0 ) = Γ 1 ( 1 ) = Γ 1 ( 2 ) = 1 .
The resemblance to the usual Γ function (of the zeroth kind Γ 0 Γ ) becomes obvious in observing that Γ 1 for integer values of the argument reduces to
Γ 1 ( n + 1 ) = 1 1 × 2 2 × 3 3 n n , n > 0 .
The generalized Γ function of the kth kind satisfies
Γ k + 1 ( x + 1 ) = x x k + 1 Γ k + 1 ( x ) , ln Γ k + 1 ( x + 1 ) = x k + 1 ln x + ln Γ k + 1 ( x ) , k = 0 : ln Γ 1 ( x + 1 ) = x ln x + ln Γ 1 ( x ) .
The analog of Raabe’s formula is given by
x 1 x d t ln Γ 1 ( t + 1 ) = x 2 2 ln x x 2 4 + L 1 1 12 .
In particular,
0 1 d x ln Γ 1 ( x + 1 ) = L 1 1 3 ,
0 1 d x ln Γ 1 ( x ) = L 1 1 12 .
Using x = 1 2 in Equation (34) and Equation (32) in
ln Γ 1 3 2 = 1 2 ln 1 2 + ln Γ 1 ( 1 2 )
we obtain
ln Γ 1 3 2 = 13 24 ln 2 + 3 2 L 1 1 8 ,
so that
0 1 2 d x ln Γ ( x + 1 ) = ln Γ 1 3 2 3 8 + 1 4 ln 2 + 1 4 ln π = 7 24 ln 2 + 3 2 L 1 1 2 + 1 4 ln π .

5. Numerical Value of L 1 from Γ 1 Function Regularization

Here, we make substantial use of the results contained in the rather elaborate paper of the authors of Ref. [7], from which we can extract the Γ 1 function regularization for L ( 1 ) ( B ) limiting ourselves to strong fields, e B / m 2 1 .
Here is the result:
lim e B m 2 L ( 1 ) ( B ) = α B 2 6 π e B m 2 + ln 2 12 L 1 .
Comparing the right-hand side with the results from the ζ function regularization in Equation (15) or from Ritus’ formula in Equation (16), the equality of the three different regularization methods teaches us
L 1 = 1 12 ζ ( 1 )
or
L 1 = C 12 + 1 12 ln 2 π 1 2 π 2 ζ ( 2 ) .
Upon using the constants
ζ ( 1 ) = 0.165421 , ζ ( 2 ) = 0.93754
and Euler’s constant
ln γ = C = Ψ ( 1 ) = 0.57721566490 ,
we obtain
L 1 = 0.248754477 .
This number can be inserted into the text wherever we meet the constant L 1 .
Given the explicit expression in Equation (46), it is worthwhile looking at the so-called Glaisher–Kinkelin constant A. Here is one of the many representations:
A = e 1 12 ζ ( 1 ) from ln A = 1 12 ζ ( 1 ) = L 1 .
When we take the value of L 1 from Equation (45), we obtain:
A = ( 2 π ) 1 12 e π 2 6 ln γ ζ ( 2 ) 1 2 π 2 = 1.2824271291
Here are two more representations:
0 1 2 ln Γ ( x ) d x = 3 2 ln A + 5 24 ln 2 + 1 4 ln π
yields
A = 2 5 36 π 1 6 exp 2 3 0 1 2 ln Γ ( x ) d x .
0 1 2 ln Γ ( x + 1 ) d x = 3 2 ln A 7 24 ln 2 + 1 4 ln π 1 2
results in
A = 2 7 36 π 1 6 exp 1 3 + 2 3 0 1 2 ln Γ ( x + 1 ) d x .
A final remark should be made concerning the importance of ζ ( 3 ) , i.e., the infinite sum of all the inverse cubes,
ζ ( 3 ) = 1 + 1 2 3 + 1 3 3 + 1 4 3 + 1 5 3 + .
No complex numbers are needed because ζ ( 3 ) lies in the original Euler domain of convergence, i.e., x > 1 . No wonder that Euler became interested in the fundamentals of its value, e.g. is it rational, transcendental, etc.? Remarkably, it took until the late 20th century before a major breakthrough was achieved in 1979 by the 61-year-old French mathematician Roger Apéry, who was able to prove the irrationality of ζ ( 3 ) . Today, due to his proof, the constant ζ ( 3 ) is known as Apéry constant. Its numerical value is given by
ζ ( 3 ) 1.202056903159 .
Hence, everywhere we meet ζ ( 3 ) in the foregoing text, one can use this value. However, ζ ( 3 ) is not only of great interest to researchers working in number theory, but is also of immense value for people doing calculations in QED, e.g. the second- and third-order (in α ) radiative corrections of the electron (and muon) anomalous magnetic moment, which is one of the best-measured and -calculated numbers in all of physics. Furthermore, ζ ( 3 ) is needed in the second-order radiative corrections to the so-called triangle anomaly, which is essential for the understanding of the π 0 decay into two photons.
Finally, while knowledge of the numerical value of Apéry’s constant ζ ( 3 ) = 4 π 2 ζ ( 2 ) 1.202056903159 is absolutely necessary for computing the just-mentioned elementary particle processes, the constants ζ ( 1 ) 0.165421 and ζ ( 2 ) 0.93754 —tied together by the relation in Equation (25)—are likewise of utmost importance for evaluating effects arising from the effective Lagrangian in QED, i.e., from the Heisenberg–Euler non-linear contribution. An example of this is the impact of the QED birefringent quantum vacuum structure on the propagation of light in the universe in the neighborhood of neutron stars, which constitute the laboratory for studying physical processes in superstrong magnetic fields.

Funding

This research received no external funding.

Acknowledgments

Part of this article was used in a talk given at a Topical Workshop on “Quantum Fields—From Fundamental Concepts to Phenomenological Questions”, 26–28 September 2018, at the Mainz Institute for Theoretical Physics. My sincere thanks to the organizers and the institute staff; their help on several occasions is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Riemann, B. Über die Anzahl der Primzahlen unter einer gegebenen Grösse. Ges. Math. Werke und Wissenschaftlicher Nachlaß 1859, 2, 145–155. [Google Scholar]
  2. Heisenberg, W.; Euler, H. Folgerungen aus der Diracschen Theorie des Positrons. Z. Phys. 1936, 98, 714. [Google Scholar] [CrossRef]
  3. Dittrich, W. One-loop effective potentials in quantum electrodynamics. J. Phys. A 1976, 9, 1171. [Google Scholar] [CrossRef]
  4. Dittrich, W.; Reuter, M. Effective Lagrangian in Quantum Electrodynamics; Lecture Notes in Physics; No. 220; Springer: Berlin, Germany, 1985. [Google Scholar]
  5. Ritus, V.I. Lagrangian of an intense electromagnetic field and quantum electrodynamics at short distances. Sov. J. Exp. Theor. Phys. 1976, 42, 774. [Google Scholar]
  6. Bendersky, L. Sur la fonction gamma géneralisée. Acta Math. 1933, 61, 263–322. [Google Scholar] [CrossRef]
  7. Dittrich, W.; Tsai, W.; Zimmermann, K.H. Evaluation of Some Integrals Following from L1, the Constant of the Asymptotic Expansion of ln Γ1(x + 1), Originating from Physics (QED). Phys. Rev. D 1979, 19, 2829. [Google Scholar]

Share and Cite

MDPI and ACS Style

Dittrich, W. The Role of Riemann’s Zeta Function in Mathematics and Physics †,‡. Universe 2019, 5, 79. https://0-doi-org.brum.beds.ac.uk/10.3390/universe5030079

AMA Style

Dittrich W. The Role of Riemann’s Zeta Function in Mathematics and Physics †,‡. Universe. 2019; 5(3):79. https://0-doi-org.brum.beds.ac.uk/10.3390/universe5030079

Chicago/Turabian Style

Dittrich, Walter. 2019. "The Role of Riemann’s Zeta Function in Mathematics and Physics †,‡" Universe 5, no. 3: 79. https://0-doi-org.brum.beds.ac.uk/10.3390/universe5030079

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop