Misplaced Pages

Riemann zeta function

Article snapshot taken from Wikipedia with creative commons attribution-sharealike license. Give it a read and then ask your questions in the chat. We can research this topic together.
(Redirected from Euler product formula) Analytic function in mathematics
The Riemann zeta function ζ(z) plotted with domain coloring.
The pole at z = 1 {\displaystyle z=1} and two zeros on the critical line.

The Riemann zeta function or Euler–Riemann zeta function, denoted by the Greek letter ζ (zeta), is a mathematical function of a complex variable defined as ζ ( s ) = n = 1 1 n s = 1 1 s + 1 2 s + 1 3 s + {\displaystyle \zeta (s)=\sum _{n=1}^{\infty }{\frac {1}{n^{s}}}={\frac {1}{1^{s}}}+{\frac {1}{2^{s}}}+{\frac {1}{3^{s}}}+\cdots } for Re ( s ) > 1 {\displaystyle \operatorname {Re} (s)>1} , and its analytic continuation elsewhere.

The Riemann zeta function plays a pivotal role in analytic number theory and has applications in physics, probability theory, and applied statistics.

Leonhard Euler first introduced and studied the function over the reals in the first half of the eighteenth century. Bernhard Riemann's 1859 article "On the Number of Primes Less Than a Given Magnitude" extended the Euler definition to a complex variable, proved its meromorphic continuation and functional equation, and established a relation between its zeros and the distribution of prime numbers. This paper also contained the Riemann hypothesis, a conjecture about the distribution of complex zeros of the Riemann zeta function that many mathematicians consider the most important unsolved problem in pure mathematics.

The values of the Riemann zeta function at even positive integers were computed by Euler. The first of them, ζ(2), provides a solution to the Basel problem. In 1979 Roger Apéry proved the irrationality of ζ(3). The values at negative integer points, also found by Euler, are rational numbers and play an important role in the theory of modular forms. Many generalizations of the Riemann zeta function, such as Dirichlet series, Dirichlet L-functions and L-functions, are known.

Definition

Bernhard Riemann's article On the number of primes below a given magnitude

The Riemann zeta function ζ(s) is a function of a complex variable s = σ + it, where σ and t are real numbers. (The notation s, σ, and t is used traditionally in the study of the zeta function, following Riemann.) When Re(s) = σ > 1, the function can be written as a converging summation or as an integral:

ζ ( s ) = n = 1 1 n s = 1 Γ ( s ) 0 x s 1 e x 1 d x , {\displaystyle \zeta (s)=\sum _{n=1}^{\infty }{\frac {1}{n^{s}}}={\frac {1}{\Gamma (s)}}\int _{0}^{\infty }{\frac {x^{s-1}}{e^{x}-1}}\,\mathrm {d} x\,,}

where

Γ ( s ) = 0 x s 1 e x d x {\displaystyle \Gamma (s)=\int _{0}^{\infty }x^{s-1}\,e^{-x}\,\mathrm {d} x}

is the gamma function. The Riemann zeta function is defined for other complex values via analytic continuation of the function defined for σ > 1.

Leonhard Euler considered the above series in 1740 for positive integer values of s, and later Chebyshev extended the definition to Re ( s ) > 1. {\displaystyle \operatorname {Re} (s)>1.}

The above series is a prototypical Dirichlet series that converges absolutely to an analytic function for s such that σ > 1 and diverges for all other values of s. Riemann showed that the function defined by the series on the half-plane of convergence can be continued analytically to all complex values s ≠ 1. For s = 1, the series is the harmonic series which diverges to +∞, and lim s 1 ( s 1 ) ζ ( s ) = 1. {\displaystyle \lim _{s\to 1}(s-1)\zeta (s)=1.} Thus the Riemann zeta function is a meromorphic function on the whole complex plane, which is holomorphic everywhere except for a simple pole at s = 1 with residue 1.

Euler's product formula

In 1737, the connection between the zeta function and prime numbers was discovered by Euler, who proved the identity

n = 1 1 n s = p  prime 1 1 p s , {\displaystyle \sum _{n=1}^{\infty }{\frac {1}{n^{s}}}=\prod _{p{\text{ prime}}}{\frac {1}{1-p^{-s}}},}

where, by definition, the left hand side is ζ(s) and the infinite product on the right hand side extends over all prime numbers p (such expressions are called Euler products):

p  prime 1 1 p s = 1 1 2 s 1 1 3 s 1 1 5 s 1 1 7 s 1 1 11 s 1 1 p s {\displaystyle \prod _{p{\text{ prime}}}{\frac {1}{1-p^{-s}}}={\frac {1}{1-2^{-s}}}\cdot {\frac {1}{1-3^{-s}}}\cdot {\frac {1}{1-5^{-s}}}\cdot {\frac {1}{1-7^{-s}}}\cdot {\frac {1}{1-11^{-s}}}\cdots {\frac {1}{1-p^{-s}}}\cdots }

Both sides of the Euler product formula converge for Re(s) > 1. The proof of Euler's identity uses only the formula for the geometric series and the fundamental theorem of arithmetic. Since the harmonic series, obtained when s = 1, diverges, Euler's formula (which becomes Πpp/p − 1⁠) implies that there are infinitely many primes. Since the logarithm of ⁠p/p − 1⁠ is approximately ⁠1/p⁠, the formula can also be used to prove the stronger result that the sum of the reciprocals of the primes is infinite. On the other hand, combining that with the sieve of Eratosthenes shows that the density of the set of primes within the set of positive integers is zero.

The Euler product formula can be used to calculate the asymptotic probability that s randomly selected integers are set-wise coprime. Intuitively, the probability that any single number is divisible by a prime (or any integer) p is ⁠1/p⁠. Hence the probability that s numbers are all divisible by this prime is ⁠1/p⁠, and the probability that at least one of them is not is 1 − ⁠1/p⁠. Now, for distinct primes, these divisibility events are mutually independent because the candidate divisors are coprime (a number is divisible by coprime divisors n and m if and only if it is divisible by nm, an event which occurs with probability ⁠1/nm⁠). Thus the asymptotic probability that s numbers are coprime is given by a product over all primes,

p  prime ( 1 1 p s ) = ( p  prime 1 1 p s ) 1 = 1 ζ ( s ) . {\displaystyle \prod _{p{\text{ prime}}}\left(1-{\frac {1}{p^{s}}}\right)=\left(\prod _{p{\text{ prime}}}{\frac {1}{1-p^{-s}}}\right)^{-1}={\frac {1}{\zeta (s)}}.}

Riemann's functional equation

This zeta function satisfies the functional equation ζ ( s ) = 2 s π s 1   sin ( π s 2 )   Γ ( 1 s )   ζ ( 1 s )   , {\displaystyle \zeta (s)=2^{s}\pi ^{s-1}\ \sin \left({\frac {\pi s}{2}}\right)\ \Gamma (1-s)\ \zeta (1-s)\ ,} where Γ(s) is the gamma function. This is an equality of meromorphic functions valid on the whole complex plane. The equation relates values of the Riemann zeta function at the points s and 1 − s, in particular relating even positive integers with odd negative integers. Owing to the zeros of the sine function, the functional equation implies that ζ(s) has a simple zero at each even negative integer s = −2n, known as the trivial zeros of ζ(s). When s is an even positive integer, the product   sin(⁠ π s / 2 ⁠) Γ(1 − s)  on the right is non-zero because Γ(1 − s) has a simple pole, which cancels the simple zero of the sine factor.

Proof of Riemann's functional equation

A proof of the functional equation proceeds as follows: We observe that if   s > 0   , {\displaystyle \ s>0\ ,} then 0 x 1 2 s 1 e n 2 π x   d x   =     Γ ( s 2 )     n s   π s 2     . {\displaystyle \int _{0}^{\infty }x^{{\frac {1}{2}}s-1}e^{-n^{2}\pi x}\ \operatorname {d} x\ =\ {\frac {\ \Gamma \!\left({\frac {s}{2}}\right)\ }{\ n^{s}\ \pi ^{\frac {s}{2}}\ }}~.}

As a result, if   s > 1   {\displaystyle \ s>1\ } then   Γ ( s 2 )   ζ ( s )     π s 2     =   n = 1   0   x s 2 1   e n 2 π x   d x   =   0 x s 2 1 n = 1 e n 2 π x   d x   , {\displaystyle {\frac {\ \Gamma \!\left({\frac {s}{2}}\right)\ \zeta (s)\ }{\ \pi ^{\frac {s}{2}}\ }}\ =\ \sum _{n=1}^{\infty }\ \int _{0}^{\infty }\ x^{{s \over 2}-1}\ e^{-n^{2}\pi x}\ \operatorname {d} x\ =\ \int _{0}^{\infty }x^{{s \over 2}-1}\sum _{n=1}^{\infty }e^{-n^{2}\pi x}\ \operatorname {d} x\ ,} with the inversion of the limiting processes justified by absolute convergence (hence the stricter requirement on s {\displaystyle s} ).

For convenience, let ψ ( x )   :=   n = 1   e n 2 π x {\displaystyle \psi (x)\ :=\ \sum _{n=1}^{\infty }\ e^{-n^{2}\pi x}}

which is a special case of the theta function. Then ζ ( s )   =   π s 2   Γ ( s 2 )     0   x 1 2 s 1   ψ ( x )   d x   . {\displaystyle \zeta (s)\ =\ {\frac {\pi ^{s \over 2}}{\ \Gamma ({s \over 2})\ }}\ \int _{0}^{\infty }\ x^{{1 \over 2}{s}-1}\ \psi (x)\ \operatorname {d} x~.}

By the Poisson summation formula we have n =   e n 2 π   x   =   1   x       n =   e   n 2 π   x   , {\displaystyle \sum _{n=-\infty }^{\infty }\ e^{-n^{2}\pi \ x}\ =\ {\frac {1}{\ {\sqrt {x\ }}\ }}\ \sum _{n=-\infty }^{\infty }\ e^{-{\frac {\ n^{2}\pi \ }{x}}}\ ,}

so that   2   ψ ( x ) + 1   =   1   x     (   2   ψ ( 1 x ) + 1   )   . {\displaystyle \ 2\ \psi (x)+1\ =\ {\frac {1}{\ {\sqrt {x\ }}\ }}\left(\ 2\ \psi \!\left({\frac {1}{x}}\right)+1\ \right)~.}

Hence π s 2   Γ ( s 2 )   ζ ( s )   =   0 1   x s 2 1   ψ ( x )   d x + 1 x s 2 1 ψ ( x )   d x   . {\displaystyle \pi ^{-{\frac {s}{2}}}\ \Gamma \!\left({\frac {s}{2}}\right)\ \zeta (s)\ =\ \int _{0}^{1}\ x^{{\frac {s}{2}}-1}\ \psi (x)\ \operatorname {d} x+\int _{1}^{\infty }x^{{\frac {s}{2}}-1}\psi (x)\ \operatorname {d} x~.}

This is equivalent to 0 1 x s 2 1 ( 1   x       ψ ( 1 x ) + 1   2 x     1 2   )   d x + 1 x s 2 1 ψ ( x )   d x {\displaystyle \int _{0}^{1}x^{{\frac {s}{2}}-1}\left({\frac {1}{\ {\sqrt {x\ }}\ }}\ \psi \!\left({\frac {1}{x}}\right)+{\frac {1}{\ 2{\sqrt {x\ }}\ }}-{\frac {1}{2}}\ \right)\ \operatorname {d} x+\int _{1}^{\infty }x^{{s \over 2}-1}\psi (x)\ \operatorname {d} x} or 1   s 1   1   s   + 0 1   x s 2 3 2   ψ ( 1   x   )   d x + 1   x s 2 1   ψ ( x )   d x   . {\displaystyle {\frac {1}{\ s-1\ }}-{\frac {1}{\ s\ }}+\int _{0}^{1}\ x^{{\frac {s}{2}}-{\frac {3}{2}}}\ \psi \!\left({\frac {1}{\ x\ }}\right)\ \operatorname {d} x+\int _{1}^{\infty }\ x^{{\frac {s}{2}}-1}\ \psi (x)\ \operatorname {d} x~.}

So π s 2   Γ (   s   2 )   ζ ( s )   =   1   s ( s 1 )   + 1   ( x s 2 1 2 + x s 2 1 )   ψ ( x )   d x {\displaystyle \pi ^{-{\frac {s}{2}}}\ \Gamma \!\left({\frac {\ s\ }{2}}\right)\ \zeta (s)\ =\ {\frac {1}{\ s(s-1)\ }}+\int _{1}^{\infty }\ \left(x^{-{\frac {s}{2}}-{\frac {1}{2}}}+x^{{\frac {s}{2}}-1}\right)\ \psi (x)\ \operatorname {d} x}

which is convergent for all s, so holds by analytic continuation. Furthermore, note by inspection that the RHS remains the same if s is replaced by 1 − s . Hence

  Γ (   s 2   )   ζ (   s   )     π s 2       =     Γ (   1 2 s 2   )   ζ (   1 s   )     π 1 2 s 2   {\displaystyle {\frac {\ \Gamma \!\left(\ {\frac {s}{2}}\ \right)\ \zeta \!\left(\ s\ \right)\ }{\ \pi ^{{\frac {s}{2}}\ }\ }}\ =\ {\frac {\ \Gamma \!\left(\ {\frac {1}{2}}-{\frac {s}{2}}\ \right)\ \zeta \!\left(\ 1-s\ \right)\ }{\ \pi ^{{\frac {1}{2}}-{\frac {s}{2}}}\ }}}

which is the functional equation attributed to Bernhard Riemann. The functional equation above can now be obtained using the duplication formula for the gamma function.

The functional equation was established by Riemann in his 1859 paper "On the Number of Primes Less Than a Given Magnitude" and used to construct the analytic continuation in the first place.

Equivalencies

An equivalent relationship had been conjectured by Euler over a hundred years earlier, in 1749, for the Dirichlet eta function (the alternating zeta function): η ( s )   =   n = 1 ( 1 ) n + 1   n s = ( 1 2 1 s )   ζ ( s )   . {\displaystyle \eta (s)\ =\ \sum _{n=1}^{\infty }{\frac {\;(-1)^{n+1}}{\ n^{s}}}=\left(1-{2^{1-s}}\right)\ \zeta (s)~.}

Incidentally, this relation gives an equation for calculating ζ(s) in the region 0 < ℛ(s) < 1 , i.e. ζ ( s ) = 1 1 2 1 s   n = 1 ( 1 ) n + 1 n s   {\displaystyle \zeta (s)={\frac {1}{\;1-2^{1-s}\ }}\sum _{n=1}^{\infty }{\frac {(-1)^{n+1}}{\;n^{s}\ }}} where the η-series is convergent (albeit non-absolutely) in the larger half-plane s > 0 (for a more detailed survey on the history of the functional equation, see e.g. Blagouchine).

Riemann also found a symmetric version of the functional equation applying to the ξ-function: ξ ( s )   =   1 2 π s 2   s ( s 1 )   Γ ( s 2 )   ζ ( s )   , {\displaystyle \xi (s)\ =\ {\frac {1}{2}}\pi ^{-{\frac {s}{2}}}\ s(s-1)\ \Gamma \!\left({\frac {s}{2}}\right)\ \zeta (s)\ ,} which satisfies: ξ ( s ) = ξ ( 1 s )   . {\displaystyle \xi (s)=\xi (1-s)~.}

(Riemann's original ξ(t) was slightly different.)

The   π s / 2   Γ ( s / 2 )   {\displaystyle \ \pi ^{-s/2}\ \Gamma (s/2)\ } factor was not well-understood at the time of Riemann, until John Tate's (1950) thesis, in which it was shown that this so-called "Gamma factor" is in fact the local L-factor corresponding to the Archimedean place, the other factors in the Euler product expansion being the local L-factors of the non-Archimedean places.

Zeros, the critical line, and the Riemann hypothesis

Main article: Riemann hypothesis
Riemann zeta spiral along the critical line from height 999000 to a million (from red to violet)
The Riemann zeta function has no zeros to the right of σ = 1 or (apart from the trivial zeros) to the left of σ = 0 (nor can the zeros lie too close to those lines). Furthermore, the non-trivial zeros are symmetric about the real axis and the line σ = ⁠1/2⁠ and, according to the Riemann hypothesis, they all lie on the line σ = ⁠1/2⁠.
This image shows a plot of the Riemann zeta function along the critical line for real values of t running from 0 to 34. The first five zeros in the critical strip are clearly visible as the place where the spirals pass through the origin.
The real part (red) and imaginary part (blue) of the Riemann zeta function along the critical line Re(s) = 1/2. The first non-trivial zeros can be seen at Im(s) = ±14.135, ±21.022 and ±25.011.
Animation showing the Riemann zeta function along the critical line. Zeta(1/2 + I y) for y ranging from 1000 to 1005.

The functional equation shows that the Riemann zeta function has zeros at −2, −4,.... These are called the trivial zeros. They are trivial in the sense that their existence is relatively easy to prove, for example, from sin ⁠πs/2⁠ being 0 in the functional equation. The non-trivial zeros have captured far more attention because their distribution not only is far less understood but, more importantly, their study yields important results concerning prime numbers and related objects in number theory. It is known that any non-trivial zero lies in the open strip { s C : 0 < Re ( s ) < 1 } {\displaystyle \{s\in \mathbb {C} :0<\operatorname {Re} (s)<1\}} , which is called the critical strip. The set { s C : Re ( s ) = 1 / 2 } {\displaystyle \{s\in \mathbb {C} :\operatorname {Re} (s)=1/2\}} is called the critical line. The Riemann hypothesis, considered one of the greatest unsolved problems in mathematics, asserts that all non-trivial zeros are on the critical line. In 1989, Conrey proved that more than 40% of the non-trivial zeros of the Riemann zeta function are on the critical line.

For the Riemann zeta function on the critical line, see Z-function.

First few nontrivial zeros
Zero
1/2 ± 14.134725 i
1/2 ± 21.022040 i
1/2 ± 25.010858 i
1/2 ± 30.424876 i
1/2 ± 32.935062 i
1/2 ± 37.586178 i
1/2 ± 40.918719 i

Number of zeros in the critical strip

Let N ( T ) {\displaystyle N(T)} be the number of zeros of ζ ( s ) {\displaystyle \zeta (s)} in the critical strip 0 < Re ( s ) < 1 {\displaystyle 0<\operatorname {Re} (s)<1} , whose imaginary parts are in the interval 0 < Im ( s ) < T {\displaystyle 0<\operatorname {Im} (s)<T} . Trudgian proved that, if T > e {\displaystyle T>e} , then

| N ( T ) T 2 π log T 2 π e | 0.112 log T + 0.278 log log T + 3.385 + 0.2 T {\displaystyle \left|N(T)-{\frac {T}{2\pi }}\log {\frac {T}{2\pi e}}\right|\leq 0.112\log T+0.278\log \log T+3.385+{\frac {0.2}{T}}} .

The Hardy–Littlewood conjectures

In 1914, G. H. Hardy proved that ζ (⁠1/2⁠ + it) has infinitely many real zeros.

Hardy and J. E. Littlewood formulated two conjectures on the density and distance between the zeros of ζ (⁠1/2⁠ + it) on intervals of large positive real numbers. In the following, N(T) is the total number of real zeros and N0(T) the total number of zeros of odd order of the function ζ (⁠1/2⁠ + it) lying in the interval (0, T].

  1. For any ε > 0, there exists a T0(ε) > 0 such that when
    T T 0 ( ε )  and  H = T 1 4 + ε , {\displaystyle T\geq T_{0}(\varepsilon )\quad {\text{ and }}\quad H=T^{{\frac {1}{4}}+\varepsilon },}
    the interval (T, T + H] contains a zero of odd order.
  2. For any ε > 0, there exists a T0(ε) > 0 and cε > 0 such that the inequality
    N 0 ( T + H ) N 0 ( T ) c ε H {\displaystyle N_{0}(T+H)-N_{0}(T)\geq c_{\varepsilon }H}
    holds when
    T T 0 ( ε )  and  H = T 1 2 + ε . {\displaystyle T\geq T_{0}(\varepsilon )\quad {\text{ and }}\quad H=T^{{\frac {1}{2}}+\varepsilon }.}

These two conjectures opened up new directions in the investigation of the Riemann zeta function.

Zero-free region

The location of the Riemann zeta function's zeros is of great importance in number theory. The prime number theorem is equivalent to the fact that there are no zeros of the zeta function on the Re(s) = 1 line. A better result that follows from an effective form of Vinogradov's mean-value theorem is that ζ (σ + it) ≠ 0 whenever σ 1 1 57.54 ( log | t | ) 2 3 ( log log | t | ) 1 3 {\displaystyle \sigma \geq 1-{\frac {1}{57.54(\log {|t|})^{\frac {2}{3}}(\log {\log {|t|}})^{\frac {1}{3}}}}} and |t| ≥ 3.

In 2015, Mossinghoff and Trudgian proved that zeta has no zeros in the region

σ 1 1 5.573412 log | t | {\displaystyle \sigma \geq 1-{\frac {1}{5.573412\log |t|}}}

for |t| ≥ 2. This is the largest known zero-free region in the critical strip for 3.06 10 10 < | t | < exp ( 10151.5 ) 5.5 10 4408 {\displaystyle 3.06\cdot 10^{10}<|t|<\exp(10151.5)\approx 5.5\cdot 10^{4408}} .

The strongest result of this kind one can hope for is the truth of the Riemann hypothesis, which would have many profound consequences in the theory of numbers.

Other results

It is known that there are infinitely many zeros on the critical line. Littlewood showed that if the sequence (γn) contains the imaginary parts of all zeros in the upper half-plane in ascending order, then

lim n ( γ n + 1 γ n ) = 0. {\displaystyle \lim _{n\rightarrow \infty }\left(\gamma _{n+1}-\gamma _{n}\right)=0.}

The critical line theorem asserts that a positive proportion of the nontrivial zeros lies on the critical line. (The Riemann hypothesis would imply that this proportion is 1.)

In the critical strip, the zero with smallest non-negative imaginary part is ⁠1/2⁠ + 14.13472514...i (OEISA058303). The fact that

ζ ( s ) = ζ ( s ¯ ) ¯ {\displaystyle \zeta (s)={\overline {\zeta ({\overline {s}})}}}

for all complex s ≠ 1 implies that the zeros of the Riemann zeta function are symmetric about the real axis. Combining this symmetry with the functional equation, furthermore, one sees that the non-trivial zeros are symmetric about the critical line Re(s) = ⁠1/2⁠.

It is also known that no zeros lie on the line with real part 1.

Specific values

Main article: Particular values of the Riemann zeta function

For any positive even integer 2n, ζ ( 2 n ) = | B 2 n | ( 2 π ) 2 n 2 ( 2 n ) ! , {\displaystyle \zeta (2n)={\frac {|{B_{2n}}|(2\pi )^{2n}}{2(2n)!}},} where B2n is the 2n-th Bernoulli number. For odd positive integers, no such simple expression is known, although these values are thought to be related to the algebraic K-theory of the integers; see Special values of L-functions.

For nonpositive integers, one has ζ ( n ) = B n + 1 n + 1 {\displaystyle \zeta (-n)=-{\frac {B_{n+1}}{n+1}}} for n ≥ 0 (using the convention that B1 = ⁠1/2⁠). In particular, ζ vanishes at the negative even integers because Bm = 0 for all odd m other than 1. These are the so-called "trivial zeros" of the zeta function.

Via analytic continuation, one can show that ζ ( 1 ) = 1 12 {\displaystyle \zeta (-1)=-{\tfrac {1}{12}}} This gives a pretext for assigning a finite value to the divergent series 1 + 2 + 3 + 4 + ⋯, which has been used in certain contexts (Ramanujan summation) such as string theory. Analogously, the particular value ζ ( 0 ) = 1 2 {\displaystyle \zeta (0)=-{\tfrac {1}{2}}} can be viewed as assigning a finite result to the divergent series 1 + 1 + 1 + 1 + ⋯.

The value ζ ( 1 2 ) = 1.46035450880958681288 {\displaystyle \zeta {\bigl (}{\tfrac {1}{2}}{\bigr )}=-1.46035450880958681288\ldots } is employed in calculating kinetic boundary layer problems of linear kinetic equations.

Although ζ ( 1 ) = 1 + 1 2 + 1 3 + {\displaystyle \zeta (1)=1+{\tfrac {1}{2}}+{\tfrac {1}{3}}+\cdots } diverges, its Cauchy principal value lim ε 0 ζ ( 1 + ε ) + ζ ( 1 ε ) 2 {\displaystyle \lim _{\varepsilon \to 0}{\frac {\zeta (1+\varepsilon )+\zeta (1-\varepsilon )}{2}}} exists and is equal to the Euler–Mascheroni constant γ = 0.5772....

The demonstration of the particular value ζ ( 2 ) = 1 + 1 2 2 + 1 3 2 + = π 2 6 {\displaystyle \zeta (2)=1+{\frac {1}{2^{2}}}+{\frac {1}{3^{2}}}+\cdots ={\frac {\pi ^{2}}{6}}} is known as the Basel problem. The reciprocal of this sum answers the question: What is the probability that two numbers selected at random are relatively prime? The value ζ ( 3 ) = 1 + 1 2 3 + 1 3 3 + = 1.202056903159594285399... {\displaystyle \zeta (3)=1+{\frac {1}{2^{3}}}+{\frac {1}{3^{3}}}+\cdots =1.202056903159594285399...} is Apéry's constant.

Taking the limit s + {\displaystyle s\rightarrow +\infty } through the real numbers, one obtains ζ ( + ) = 1 {\displaystyle \zeta (+\infty )=1} . But at complex infinity on the Riemann sphere the zeta function has an essential singularity.

Various properties

For sums involving the zeta function at integer and half-integer values, see rational zeta series.

Reciprocal

The reciprocal of the zeta function may be expressed as a Dirichlet series over the Möbius function μ(n):

1 ζ ( s ) = n = 1 μ ( n ) n s {\displaystyle {\frac {1}{\zeta (s)}}=\sum _{n=1}^{\infty }{\frac {\mu (n)}{n^{s}}}}

for every complex number s with real part greater than 1. There are a number of similar relations involving various well-known multiplicative functions; these are given in the article on the Dirichlet series.

The Riemann hypothesis is equivalent to the claim that this expression is valid when the real part of s is greater than ⁠1/2⁠.

Universality

The critical strip of the Riemann zeta function has the remarkable property of universality. This zeta function universality states that there exists some location on the critical strip that approximates any holomorphic function arbitrarily well. Since holomorphic functions are very general, this property is quite remarkable. The first proof of universality was provided by Sergei Mikhailovitch Voronin in 1975. More recent work has included effective versions of Voronin's theorem and extending it to Dirichlet L-functions.

Estimates of the maximum of the modulus of the zeta function

Let the functions F(T;H) and G(s0;Δ) be defined by the equalities

F ( T ; H ) = max | t T | H | ζ ( 1 2 + i t ) | , G ( s 0 ; Δ ) = max | s s 0 | Δ | ζ ( s ) | . {\displaystyle F(T;H)=\max _{|t-T|\leq H}\left|\zeta \left({\tfrac {1}{2}}+it\right)\right|,\qquad G(s_{0};\Delta )=\max _{|s-s_{0}|\leq \Delta }|\zeta (s)|.}

Here T is a sufficiently large positive number, 0 < H ≪ log log T, s0 = σ0 + iT, ⁠1/2⁠ ≤ σ0 ≤ 1, 0 < Δ < ⁠1/3⁠. Estimating the values F and G from below shows, how large (in modulus) values ζ(s) can take on short intervals of the critical line or in small neighborhoods of points lying in the critical strip 0 ≤ Re(s) ≤ 1.

The case H ≫ log log T was studied by Kanakanahalli Ramachandra; the case Δ > c, where c is a sufficiently large constant, is trivial.

Anatolii Karatsuba proved, in particular, that if the values H and Δ exceed certain sufficiently small constants, then the estimates

F ( T ; H ) T c 1 , G ( s 0 ; Δ ) T c 2 , {\displaystyle F(T;H)\geq T^{-c_{1}},\qquad G(s_{0};\Delta )\geq T^{-c_{2}},}

hold, where c1 and c2 are certain absolute constants.

The argument of the Riemann zeta function

The function

S ( t ) = 1 π arg ζ ( 1 2 + i t ) {\displaystyle S(t)={\frac {1}{\pi }}\arg {\zeta \left({\tfrac {1}{2}}+it\right)}}

is called the argument of the Riemann zeta function. Here arg ζ(⁠1/2⁠ + it) is the increment of an arbitrary continuous branch of arg ζ(s) along the broken line joining the points 2, 2 + it and ⁠1/2⁠ + it.

There are some theorems on properties of the function S(t). Among those results are the mean value theorems for S(t) and its first integral

S 1 ( t ) = 0 t S ( u ) d u {\displaystyle S_{1}(t)=\int _{0}^{t}S(u)\,\mathrm {d} u}

on intervals of the real line, and also the theorem claiming that every interval (T, T + H] for

H T 27 82 + ε {\displaystyle H\geq T^{{\frac {27}{82}}+\varepsilon }}

contains at least

H ln T 3 e c ln ln T {\displaystyle H{\sqrt{\ln T}}e^{-c{\sqrt {\ln \ln T}}}}

points where the function S(t) changes sign. Earlier similar results were obtained by Atle Selberg for the case

H T 1 2 + ε . {\displaystyle H\geq T^{{\frac {1}{2}}+\varepsilon }.}

Representations

Dirichlet series

An extension of the area of convergence can be obtained by rearranging the original series. The series

ζ ( s ) = 1 s 1 n = 1 ( n ( n + 1 ) s n s n s ) {\displaystyle \zeta (s)={\frac {1}{s-1}}\sum _{n=1}^{\infty }\left({\frac {n}{(n+1)^{s}}}-{\frac {n-s}{n^{s}}}\right)}

converges for Re(s) > 0, while

ζ ( s ) = 1 s 1 n = 1 n ( n + 1 ) 2 ( 2 n + 3 + s ( n + 1 ) s + 2 2 n 1 s n s + 2 ) {\displaystyle \zeta (s)={\frac {1}{s-1}}\sum _{n=1}^{\infty }{\frac {n(n+1)}{2}}\left({\frac {2n+3+s}{(n+1)^{s+2}}}-{\frac {2n-1-s}{n^{s+2}}}\right)}

converge even for Re(s) > −1. In this way, the area of convergence can be extended to Re(s) > −k for any negative integer −k.

The recurrence connection is clearly visible from the expression valid for Re(s) > −2 enabling further expansion by integration by parts.

ζ ( s ) = 1 + 1 s 1 s 2 ! [ ζ ( s + 1 ) 1 ] s ( s + 1 ) 3 ! [ ζ ( s + 2 ) 1 ] s ( s + 1 ) ( s + 2 ) 3 ! n = 1 0 1 t 3 d t ( n + t ) s + 3 {\displaystyle {\begin{aligned}\zeta (s)=&1+{\frac {1}{s-1}}-{\frac {s}{2!}}\\-&{\frac {s(s+1)}{3!}}\\&-{\frac {s(s+1)(s+2)}{3!}}\sum _{n=1}^{\infty }\int _{0}^{1}{\frac {t^{3}dt}{(n+t)^{s+3}}}\end{aligned}}}

Mellin-type integrals

The Mellin transform of a function f(x) is defined as

0 f ( x ) x s d x x {\displaystyle \int _{0}^{\infty }f(x)x^{s}\,{\frac {\mathrm {d} x}{x}}}

in the region where the integral is defined. There are various expressions for the zeta function as Mellin transform-like integrals. If the real part of s is greater than one, we have

Γ ( s ) ζ ( s ) = 0 x s 1 e x 1 d x {\displaystyle \Gamma (s)\zeta (s)=\int _{0}^{\infty }{\frac {x^{s-1}}{e^{x}-1}}\,\mathrm {d} x\quad } and Γ ( s ) ζ ( s ) = 1 2 s 0 x s cosh ( x ) 1 d x {\displaystyle \quad \Gamma (s)\zeta (s)={\frac {1}{2s}}\int _{0}^{\infty }{\frac {x^{s}}{\cosh(x)-1}}\,\mathrm {d} x} ,

where Γ denotes the gamma function. By modifying the contour, Riemann showed that

2 sin ( π s ) Γ ( s ) ζ ( s ) = i H ( x ) s 1 e x 1 d x {\displaystyle 2\sin(\pi s)\Gamma (s)\zeta (s)=i\oint _{H}{\frac {(-x)^{s-1}}{e^{x}-1}}\,\mathrm {d} x}

for all s (where H denotes the Hankel contour).

We can also find expressions which relate to prime numbers and the prime number theorem. If π(x) is the prime-counting function, then

ln ζ ( s ) = s 0 π ( x ) x ( x s 1 ) d x , {\displaystyle \ln \zeta (s)=s\int _{0}^{\infty }{\frac {\pi (x)}{x(x^{s}-1)}}\,\mathrm {d} x,}

for values with Re(s) > 1.

A similar Mellin transform involves the Riemann function J(x), which counts prime powers p with a weight of ⁠1/n⁠, so that

J ( x ) = π ( x 1 n ) n . {\displaystyle J(x)=\sum {\frac {\pi \left(x^{\frac {1}{n}}\right)}{n}}.}

Now

ln ζ ( s ) = s 0 J ( x ) x s 1 d x . {\displaystyle \ln \zeta (s)=s\int _{0}^{\infty }J(x)x^{-s-1}\,\mathrm {d} x.}

These expressions can be used to prove the prime number theorem by means of the inverse Mellin transform. Riemann's prime-counting function is easier to work with, and π(x) can be recovered from it by Möbius inversion.

Theta functions

The Riemann zeta function can be given by a Mellin transform

2 π s 2 Γ ( s 2 ) ζ ( s ) = 0 ( θ ( i t ) 1 ) t s 2 1 d t , {\displaystyle 2\pi ^{-{\frac {s}{2}}}\Gamma \left({\frac {s}{2}}\right)\zeta (s)=\int _{0}^{\infty }{\bigl (}\theta (it)-1{\bigr )}t^{{\frac {s}{2}}-1}\,\mathrm {d} t,}

in terms of Jacobi's theta function

θ ( τ ) = n = e π i n 2 τ . {\displaystyle \theta (\tau )=\sum _{n=-\infty }^{\infty }e^{\pi in^{2}\tau }.}

However, this integral only converges if the real part of s is greater than 1, but it can be regularized. This gives the following expression for the zeta function, which is well defined for all s except 0 and 1:

π s 2 Γ ( s 2 ) ζ ( s ) = 1 s 1 1 s + 1 2 0 1 ( θ ( i t ) t 1 2 ) t s 2 1 d t + 1 2 1 ( θ ( i t ) 1 ) t s 2 1 d t . {\displaystyle \pi ^{-{\frac {s}{2}}}\Gamma \left({\frac {s}{2}}\right)\zeta (s)={\frac {1}{s-1}}-{\frac {1}{s}}+{\frac {1}{2}}\int _{0}^{1}\left(\theta (it)-t^{-{\frac {1}{2}}}\right)t^{{\frac {s}{2}}-1}\,\mathrm {d} t+{\frac {1}{2}}\int _{1}^{\infty }{\bigl (}\theta (it)-1{\bigr )}t^{{\frac {s}{2}}-1}\,\mathrm {d} t.}

Laurent series

The Riemann zeta function is meromorphic with a single pole of order one at s = 1. It can therefore be expanded as a Laurent series about s = 1; the series development is then

ζ ( s ) = 1 s 1 + n = 0 γ n n ! ( 1 s ) n . {\displaystyle \zeta (s)={\frac {1}{s-1}}+\sum _{n=0}^{\infty }{\frac {\gamma _{n}}{n!}}(1-s)^{n}.}

The constants γn here are called the Stieltjes constants and can be defined by the limit

γ n = lim m ( ( k = 1 m ( ln k ) n k ) ( ln m ) n + 1 n + 1 ) . {\displaystyle \gamma _{n}=\lim _{m\rightarrow \infty }{\left(\left(\sum _{k=1}^{m}{\frac {(\ln k)^{n}}{k}}\right)-{\frac {(\ln m)^{n+1}}{n+1}}\right)}.}

The constant term γ0 is the Euler–Mascheroni constant.

Integral

For all s ∈ ℂ, s ≠ 1, the integral relation (cf. Abel–Plana formula)

  ζ ( s )   =   1   s 1   +   1   2 + 2 0 sin (   s   arctan t   )   ( 1 + t 2 ) s / 2 ( e 2 π t 1 )     d t   {\displaystyle \ \zeta (s)\ =\ {\frac {1}{\ s-1\ }}+{\frac {\ 1\ }{2}}+2\int _{0}^{\infty }{\frac {\sin(\ s\ \arctan t\ )}{\ \left(1+t^{2}\right)^{s/2}\left(e^{2\pi t}-1\right)\ }}\ \operatorname {d} t\ }

holds true, which may be used for a numerical evaluation of the zeta function.

Rising factorial

Another series development using the rising factorial valid for the entire complex plane is

ζ ( s ) = s s 1 n = 1 ( ζ ( s + n ) 1 ) s ( s + 1 ) ( s + n 1 ) ( n + 1 ) ! . {\displaystyle \zeta (s)={\frac {s}{s-1}}-\sum _{n=1}^{\infty }{\bigl (}\zeta (s+n)-1{\bigr )}{\frac {s(s+1)\cdots (s+n-1)}{(n+1)!}}.}

This can be used recursively to extend the Dirichlet series definition to all complex numbers.

The Riemann zeta function also appears in a form similar to the Mellin transform in an integral over the Gauss–Kuzmin–Wirsing operator acting on x; that context gives rise to a series expansion in terms of the falling factorial.

Hadamard product

On the basis of Weierstrass's factorization theorem, Hadamard gave the infinite product expansion

ζ ( s ) = e ( log ( 2 π ) 1 γ 2 ) s 2 ( s 1 ) Γ ( 1 + s 2 ) ρ ( 1 s ρ ) e s ρ , {\displaystyle \zeta (s)={\frac {e^{\left(\log(2\pi )-1-{\frac {\gamma }{2}}\right)s}}{2(s-1)\Gamma \left(1+{\frac {s}{2}}\right)}}\prod _{\rho }\left(1-{\frac {s}{\rho }}\right)e^{\frac {s}{\rho }},}

where the product is over the non-trivial zeros ρ of ζ and the letter γ again denotes the Euler–Mascheroni constant. A simpler infinite product expansion is

ζ ( s ) = π s 2 ρ ( 1 s ρ ) 2 ( s 1 ) Γ ( 1 + s 2 ) . {\displaystyle \zeta (s)=\pi ^{\frac {s}{2}}{\frac {\prod _{\rho }\left(1-{\frac {s}{\rho }}\right)}{2(s-1)\Gamma \left(1+{\frac {s}{2}}\right)}}.}

This form clearly displays the simple pole at s = 1, the trivial zeros at −2, −4, ... due to the gamma function term in the denominator, and the non-trivial zeros at s = ρ. (To ensure convergence in the latter formula, the product should be taken over "matching pairs" of zeros, i.e. the factors for a pair of zeros of the form ρ and 1 − ρ should be combined.)

Globally convergent series

A globally convergent series for the zeta function, valid for all complex numbers s except s = 1 + ⁠2πi/ln 2⁠n for some integer n, was conjectured by Konrad Knopp in 1926 and proven by Helmut Hasse in 1930 (cf. Euler summation):

ζ ( s ) = 1 1 2 1 s n = 0 1 2 n + 1 k = 0 n ( n k ) ( 1 ) k ( k + 1 ) s . {\displaystyle \zeta (s)={\frac {1}{1-2^{1-s}}}\sum _{n=0}^{\infty }{\frac {1}{2^{n+1}}}\sum _{k=0}^{n}{\binom {n}{k}}{\frac {(-1)^{k}}{(k+1)^{s}}}.}

The series appeared in an appendix to Hasse's paper, and was published for the second time by Jonathan Sondow in 1994.

Hasse also proved the globally converging series

ζ ( s ) = 1 s 1 n = 0 1 n + 1 k = 0 n ( n k ) ( 1 ) k ( k + 1 ) s 1 {\displaystyle \zeta (s)={\frac {1}{s-1}}\sum _{n=0}^{\infty }{\frac {1}{n+1}}\sum _{k=0}^{n}{\binom {n}{k}}{\frac {(-1)^{k}}{(k+1)^{s-1}}}}

in the same publication. Research by Iaroslav Blagouchine has found that a similar, equivalent series was published by Joseph Ser in 1926.

In 1997 K. Maślanka gave another globally convergent (except s = 1) series for the Riemann zeta function:

ζ ( s ) = 1 s 1 k = 0 ( i = 1 k ( i s 2 ) ) A k k ! = 1 s 1 k = 0 ( 1 s 2 ) k A k k ! {\displaystyle \zeta (s)={\frac {1}{s-1}}\sum _{k=0}^{\infty }{\biggl (}\prod _{i=1}^{k}(i-{\frac {s}{2}}){\biggl )}{\frac {A_{k}}{k!}}={\frac {1}{s-1}}\sum _{k=0}^{\infty }{\biggl (}1-{\frac {s}{2}}{\biggl )}_{k}{\frac {A_{k}}{k!}}}

where real coefficients A k {\displaystyle A_{k}} are given by:

A k = j = 0 k ( 1 ) j ( k j ) ( 2 j + 1 ) ζ ( 2 j + 2 ) = j = 0 k ( k j ) B 2 j + 2 π 2 j + 2 ( 2 ) j ( 1 2 ) j {\displaystyle A_{k}=\sum _{j=0}^{k}(-1)^{j}{\binom {k}{j}}(2j+1)\zeta (2j+2)=\sum _{j=0}^{k}{\binom {k}{j}}{\frac {B_{2j+2}\pi ^{2j+2}}{\left(2\right)_{j}\left({\frac {1}{2}}\right)_{j}}}}

Here B n {\displaystyle B_{n}} are the Bernoulli numbers and ( x ) k {\displaystyle (x)_{k}} denotes the Pochhammer symbol.

Note that this representation of the zeta function is essentially an interpolation with nodes, where the nodes are points s = 2 , 4 , 6 , {\displaystyle s=2,4,6,\ldots } , i.e. exactly those where the zeta values are precisely known, as Euler showed. An elegant and very short proof of this representation of the zeta function, based on Carlson's theorem, was presented by Philippe Flajolet in 2006.

The asymptotic behavior of the coefficients A k {\displaystyle A_{k}} is rather curious: for growing k {\displaystyle k} values, we observe regular oscillations with a nearly exponentially decreasing amplitude and slowly decreasing frequency (roughly as k 2 / 3 {\displaystyle k^{-2/3}} ). Using the saddle point method, we can show that

A k 4 π 3 / 2 3 κ exp ( 3 κ 2 + π 2 4 κ ) cos ( 4 π 3 3 3 κ 2 + 3 π 2 4 κ ) {\displaystyle A_{k}\sim {\frac {4\pi ^{3/2}}{\sqrt {3\kappa }}}\exp {\biggl (}-{\frac {3\kappa }{2}}+{\frac {\pi ^{2}}{4\kappa }}{\biggl )}\cos {\biggl (}{\frac {4\pi }{3}}-{\frac {3{\sqrt {3}}\kappa }{2}}+{\frac {{\sqrt {3}}\pi ^{2}}{4\kappa }}{\biggl )}}

where κ {\displaystyle \kappa } stands for:

κ := π 2 k 3 {\displaystyle \kappa :={\sqrt{\pi ^{2}k}}}

(see for details).

On the basis of this representation, in 2003 Luis Báez-Duarte provided a new criterion for the Riemann hypothesis. Namely, if we define the coefficients c k {\displaystyle c_{k}} as

c k := j = 0 k ( 1 ) j ( k j ) 1 ζ ( 2 j + 2 ) {\displaystyle c_{k}:=\sum _{j=0}^{k}(-1)^{j}{\binom {k}{j}}{\frac {1}{\zeta (2j+2)}}}

then the Riemann hypothesis is equivalent to

c k = O ( k 3 / 4 + ε ) ( ε > 0 ) {\displaystyle c_{k}={\mathcal {O}}{\biggl (}k^{-3/4+\varepsilon }{\biggl )}\qquad (\forall \varepsilon >0)}

Rapidly convergent series

Peter Borwein developed an algorithm that applies Chebyshev polynomials to the Dirichlet eta function to produce a very rapidly convergent series suitable for high precision numerical calculations.

Series representation at positive integers via the primorial

ζ ( k ) = 2 k 2 k 1 + r = 2 ( p r 1 # ) k J k ( p r # ) k = 2 , 3 , . {\displaystyle \zeta (k)={\frac {2^{k}}{2^{k}-1}}+\sum _{r=2}^{\infty }{\frac {(p_{r-1}\#)^{k}}{J_{k}(p_{r}\#)}}\qquad k=2,3,\ldots .}

Here pn# is the primorial sequence and Jk is Jordan's totient function.

Series representation by the incomplete poly-Bernoulli numbers

The function ζ can be represented, for Re(s) > 1, by the infinite series

ζ ( s ) = n = 0 B n , 2 ( s ) ( W k ( 1 ) ) n n ! , {\displaystyle \zeta (s)=\sum _{n=0}^{\infty }B_{n,\geq 2}^{(s)}{\frac {(W_{k}(-1))^{n}}{n!}},}

where k ∈ {−1, 0}, Wk is the kth branch of the Lambert W-function, and B
n, ≥2 is an incomplete poly-Bernoulli number.

The Mellin transform of the Engel map

The function g ( x ) = x ( 1 + x 1 ) 1 {\displaystyle g(x)=x\left(1+\left\lfloor x^{-1}\right\rfloor \right)-1} is iterated to find the coefficients appearing in Engel expansions.

The Mellin transform of the map g ( x ) {\displaystyle g(x)} is related to the Riemann zeta function by the formula

0 1 g ( x ) x s 1 d x = n = 1 1 n + 1 1 n ( x ( n + 1 ) 1 ) x s 1 d x = n = 1 n s ( s 1 ) + ( n + 1 ) s 1 ( n 2 + 2 n + 1 ) + n s 1 s n 1 s ( s + 1 ) s ( n + 1 ) = ζ ( s + 1 ) s + 1 1 s ( s + 1 ) {\displaystyle {\begin{aligned}\int _{0}^{1}g(x)x^{s-1}\,dx&=\sum _{n=1}^{\infty }\int _{\frac {1}{n+1}}^{\frac {1}{n}}(x(n+1)-1)x^{s-1}\,dx\\&=\sum _{n=1}^{\infty }{\frac {n^{-s}(s-1)+(n+1)^{-s-1}(n^{2}+2n+1)+n^{-s-1}s-n^{1-s}}{(s+1)s(n+1)}}\\&={\frac {\zeta (s+1)}{s+1}}-{\frac {1}{s(s+1)}}\end{aligned}}}

Thue-Morse sequence

Certain linear combinations of Dirichlet series whose coefficients are terms of the Thue-Morse sequence give rise to identities involving the Riemann Zeta function (Tóth, 2022 ). For instance:

n 1 5 t n 1 + 3 t n n 2 = 4 ζ ( 2 ) = 2 π 2 3 , n 1 9 t n 1 + 7 t n n 3 = 8 ζ ( 3 ) , {\displaystyle {\begin{aligned}\sum _{n\geq 1}{\frac {5t_{n-1}+3t_{n}}{n^{2}}}&=4\zeta (2)={\frac {2\pi ^{2}}{3}},\\\sum _{n\geq 1}{\frac {9t_{n-1}+7t_{n}}{n^{3}}}&=8\zeta (3),\end{aligned}}}

where ( t n ) n 0 {\displaystyle (t_{n})_{n\geq 0}} is the n t h {\displaystyle n^{\rm {th}}} term of the Thue-Morse sequence. In fact, for all s {\displaystyle s} with real part greater than 1 {\displaystyle 1} , we have

( 2 s + 1 ) n 1 t n 1 n s + ( 2 s 1 ) n 1 t n n s = 2 s ζ ( s ) . {\displaystyle (2^{s}+1)\sum _{n\geq 1}{\frac {t_{n-1}}{n^{s}}}+(2^{s}-1)\sum _{n\geq 1}{\frac {t_{n}}{n^{s}}}=2^{s}\zeta (s).}

In nth dimensions

The zeta function can also be represented as an nth amount of integrals:

ζ ( n ) = [ 0 , 1 ] n d z 1 d z 2 d z n 1 z 1 z 2 z n {\displaystyle \zeta (n)=\int _{\cdots }\iint _{^{n}}{\frac {dz_{1}\,dz_{2}\ldots \,dz_{n}}{1-z_{1}\,z_{2}\ldots \,z_{n}}}}

and it only works for n N / { 1 } {\displaystyle n\in \mathbb {N} /\{1\}}

Numerical algorithms

A classical algorithm, in use prior to about 1930, proceeds by applying the Euler-Maclaurin formula to obtain, for n and m positive integers,

ζ ( s ) = j = 1 n 1 j s + 1 2 n s + n 1 s s 1 + k = 1 m T k , n ( s ) + E m , n ( s ) {\displaystyle \zeta (s)=\sum _{j=1}^{n-1}j^{-s}+{\tfrac {1}{2}}n^{-s}+{\frac {n^{1-s}}{s-1}}+\sum _{k=1}^{m}T_{k,n}(s)+E_{m,n}(s)}

where, letting B 2 k {\displaystyle B_{2k}} denote the indicated Bernoulli number,

T k , n ( s ) = B 2 k ( 2 k ) ! n 1 s 2 k j = 0 2 k 2 ( s + j ) {\displaystyle T_{k,n}(s)={\frac {B_{2k}}{(2k)!}}n^{1-s-2k}\prod _{j=0}^{2k-2}(s+j)}

and the error satisfies

| E m , n ( s ) | < | s + 2 m + 1 σ + 2 m + 1 T m + 1 , n ( s ) | , {\displaystyle |E_{m,n}(s)|<\left|{\frac {s+2m+1}{\sigma +2m+1}}T_{m+1,n}(s)\right|,}

with σ = Re(s).

A modern numerical algorithm is the Odlyzko–Schönhage algorithm.

Applications

The zeta function occurs in applied statistics including Zipf's law, Zipf–Mandelbrot law, and Lotka's law.

Zeta function regularization is used as one possible means of regularization of divergent series and divergent integrals in quantum field theory. In one notable example, the Riemann zeta function shows up explicitly in one method of calculating the Casimir effect. The zeta function is also useful for the analysis of dynamical systems.

Musical tuning

In the theory of musical tunings, the zeta function can be used to find equal divisions of the octave (EDOs) that closely approximate the intervals of the harmonic series. For increasing values of t R {\displaystyle t\in \mathbb {R} } , the value of

| ζ ( 1 2 + 2 π i ln ( 2 ) t ) | {\displaystyle \left\vert \zeta \left({\frac {1}{2}}+{\frac {2\pi {i}}{\ln {(2)}}}t\right)\right\vert }

peaks near integers that correspond to such EDOs. Examples include popular choices such as 12, 19, and 53.

Infinite series

The zeta function evaluated at equidistant positive integers appears in infinite series representations of a number of constants.

  • n = 2 ( ζ ( n ) 1 ) = 1 {\displaystyle \sum _{n=2}^{\infty }{\bigl (}\zeta (n)-1{\bigr )}=1}

In fact the even and odd terms give the two sums

  • n = 1 ( ζ ( 2 n ) 1 ) = 3 4 {\displaystyle \sum _{n=1}^{\infty }{\bigl (}\zeta (2n)-1{\bigr )}={\frac {3}{4}}}

and

  • n = 1 ( ζ ( 2 n + 1 ) 1 ) = 1 4 {\displaystyle \sum _{n=1}^{\infty }{\bigl (}\zeta (2n+1)-1{\bigr )}={\frac {1}{4}}}

Parametrized versions of the above sums are given by

  • n = 1 ( ζ ( 2 n ) 1 ) t 2 n = t 2 t 2 1 + 1 2 ( 1 π t cot ( t π ) ) {\displaystyle \sum _{n=1}^{\infty }(\zeta (2n)-1)\,t^{2n}={\frac {t^{2}}{t^{2}-1}}+{\frac {1}{2}}\left(1-\pi t\cot(t\pi )\right)}

and

  • n = 1 ( ζ ( 2 n + 1 ) 1 ) t 2 n = t 2 t 2 1 1 2 ( ψ 0 ( t ) + ψ 0 ( t ) ) γ {\displaystyle \sum _{n=1}^{\infty }(\zeta (2n+1)-1)\,t^{2n}={\frac {t^{2}}{t^{2}-1}}-{\frac {1}{2}}\left(\psi ^{0}(t)+\psi ^{0}(-t)\right)-\gamma }

with | t | < 2 {\displaystyle |t|<2} and where ψ {\displaystyle \psi } and γ {\displaystyle \gamma } are the polygamma function and Euler's constant, respectively, as well as

  • n = 1 ζ ( 2 n ) 1 n t 2 n = log ( 1 t 2 sinc ( π t ) ) {\displaystyle \sum _{n=1}^{\infty }{\frac {\zeta (2n)-1}{n}}\,t^{2n}=\log \left({\dfrac {1-t^{2}}{\operatorname {sinc} (\pi \,t)}}\right)}

all of which are continuous at t = 1 {\displaystyle t=1} . Other sums include

  • n = 2 ζ ( n ) 1 n = 1 γ {\displaystyle \sum _{n=2}^{\infty }{\frac {\zeta (n)-1}{n}}=1-\gamma }
  • n = 1 ζ ( 2 n ) 1 n = ln 2 {\displaystyle \sum _{n=1}^{\infty }{\frac {\zeta (2n)-1}{n}}=\ln 2}
  • n = 2 ζ ( n ) 1 n ( ( 3 2 ) n 1 1 ) = 1 3 ln π {\displaystyle \sum _{n=2}^{\infty }{\frac {\zeta (n)-1}{n}}\left(\left({\tfrac {3}{2}}\right)^{n-1}-1\right)={\frac {1}{3}}\ln \pi }
  • n = 1 ( ζ ( 4 n ) 1 ) = 7 8 π 4 ( e 2 π + 1 e 2 π 1 ) . {\displaystyle \sum _{n=1}^{\infty }{\bigl (}\zeta (4n)-1{\bigr )}={\frac {7}{8}}-{\frac {\pi }{4}}\left({\frac {e^{2\pi }+1}{e^{2\pi }-1}}\right).}
  • n = 2 ζ ( n ) 1 n ( ( 1 + i ) n 1 i n ) = π 4 {\displaystyle \sum _{n=2}^{\infty }{\frac {\zeta (n)-1}{n}}\Im {\bigl (}(1+i)^{n}-1-i^{n}{\bigr )}={\frac {\pi }{4}}}

where {\displaystyle \Im } denotes the imaginary part of a complex number.

Another interesting series that relates to the natural logarithm of the lemniscate constant is the following

  • n = 2 [ 2 ( 1 ) n ζ ( n ) 4 n n ( 1 ) n ζ ( n ) 2 n n ] = ln ( ϖ 2 2 ) {\displaystyle \sum _{n=2}^{\infty }\left=\ln \left({\frac {\varpi }{2{\sqrt {2}}}}\right)}


There are yet more formulas in the article Harmonic number.

Generalizations

There are a number of related zeta functions that can be considered to be generalizations of the Riemann zeta function. These include the Hurwitz zeta function

ζ ( s , q ) = k = 0 1 ( k + q ) s {\displaystyle \zeta (s,q)=\sum _{k=0}^{\infty }{\frac {1}{(k+q)^{s}}}}

(the convergent series representation was given by Helmut Hasse in 1930, cf. Hurwitz zeta function), which coincides with the Riemann zeta function when q = 1 (the lower limit of summation in the Hurwitz zeta function is 0, not 1), the Dirichlet L-functions and the Dedekind zeta function. For other related functions see the articles zeta function and L-function.

The polylogarithm is given by

Li s ( z ) = k = 1 z k k s {\displaystyle \operatorname {Li} _{s}(z)=\sum _{k=1}^{\infty }{\frac {z^{k}}{k^{s}}}}

which coincides with the Riemann zeta function when z = 1. The Clausen function Cls(θ) can be chosen as the real or imaginary part of Lis(e).

The Lerch transcendent is given by

Φ ( z , s , q ) = k = 0 z k ( k + q ) s {\displaystyle \Phi (z,s,q)=\sum _{k=0}^{\infty }{\frac {z^{k}}{(k+q)^{s}}}}

which coincides with the Riemann zeta function when z = 1 and q = 1 (the lower limit of summation in the Lerch transcendent is 0, not 1).

The multiple zeta functions are defined by

ζ ( s 1 , s 2 , , s n ) = k 1 > k 2 > > k n > 0 k 1 s 1 k 2 s 2 k n s n . {\displaystyle \zeta (s_{1},s_{2},\ldots ,s_{n})=\sum _{k_{1}>k_{2}>\cdots >k_{n}>0}{k_{1}}^{-s_{1}}{k_{2}}^{-s_{2}}\cdots {k_{n}}^{-s_{n}}.}

One can analytically continue these functions to the n-dimensional complex space. The special values taken by these functions at positive integer arguments are called multiple zeta values by number theorists and have been connected to many different branches in mathematics and physics.

See also

References

  1. "Jupyter Notebook Viewer". Nbviewer.ipython.org. Retrieved 4 January 2017.
  2. ^ Steuding, Jörn; Suriajaya, Ade Irma (1 November 2020). "Value-Distribution of the Riemann Zeta-Function Along Its Julia Lines". Computational Methods and Function Theory. 20 (3): 389–401. doi:10.1007/s40315-020-00316-x. hdl:2324/4483207. ISSN 2195-3724. S2CID 216323223. Theorem 2 implies that ζ has an essential singularity at infinity
  3. Bombieri, Enrico. "The Riemann Hypothesis – official problem description" (PDF). Clay Mathematics Institute. Archived from the original (PDF) on 22 December 2015. Retrieved 8 August 2014.
  4. Devlin, Keith (2002). The Millennium Problems: The seven greatest unsolved mathematical puzzles of our time. New York: Barnes & Noble. pp. 43–47. ISBN 978-0-7607-8659-8.
  5. Sandifer, Charles Edward (2007). How Euler Did It. Mathematical Association of America. p. 193. ISBN 978-0-88385-563-8.
  6. Titchmarsh, E.C. (1986). The Theory of the Riemann Zeta Function (2nd ed.). Oxford, UK: Oxford Science Publications. pp. 21–22. ISBN 0-19-853369-1.
  7. Blagouchine, I.V. (1 March 2018). The history of the functional equation of the zeta-function. Seminar on the History of Mathematics. St. Petersburg, RU: Steklov Institute of Mathematics; "online PDF". Archived from the original on 2 May 2018. Retrieved 2 May 2018.
  8. Blagouchine, I.V. (2014). "Rediscovery of Malmsten's integrals, their evaluation by contour integration methods and some related results". The Ramanujan Journal. 35 (1): 21–110. doi:10.1007/s11139-013-9528-5. S2CID 120943474.
    Blagouchine, I.V. (2017). "Addendum". The Ramanujan Journal. 42: 777–781. doi:10.1007/s11139-015-9763-z. S2CID 125198685. Archived from the original on 2 May 2018. Retrieved 2 May 2018.
  9. Conrey, J. B. (1989). "More than two fifths of the zeros of the Riemann zeta function are on the critical line". J. Reine Angew. Math. 1989 (399): 1–26. doi:10.1515/crll.1989.399.1. MR 1004130. S2CID 115910600.
  10. Eric Weisstein. "Riemann Zeta Function Zeros". Retrieved 24 April 2021.
  11. The L-functions and Modular Forms Database. "Zeros of ζ(s)".
  12. Trudgian, Timothy S. (2014). "An improved upper bound for the argument of the Riemann zeta function on the critical line II". J. Number Theory. 134: 280–292. arXiv:1208.5846. doi:10.1016/j.jnt.2013.07.017.
  13. Hardy, G.H. (1914). "Sur les zeros de la fonction ζ(s)". Comptes rendus de l'Académie des Sciences. 158. French Academy of Sciences: 1012–1014.
  14. Hardy, G. H.; Fekete, M.; Littlewood, J. E. (1 September 1921). "The Zeros of Riemann's Zeta-Function on the Critical Line". Journal of the London Mathematical Society. s1-1: 15–19. doi:10.1112/jlms/s1-1.1.15.
  15. Diamond, Harold G. (1982). "Elementary methods in the study of the distribution of prime numbers". Bulletin of the American Mathematical Society. 7 (3): 553–89. doi:10.1090/S0273-0979-1982-15057-1. MR 0670132.
  16. Ford, K. (2002). "Vinogradov's integral and bounds for the Riemann zeta function". Proc. London Math. Soc. 85 (3): 565–633. arXiv:1910.08209. doi:10.1112/S0024611502013655. S2CID 121144007.
  17. Mossinghoff, Michael J.; Trudgian, Timothy S. (2015). "Nonnegative trigonometric polynomials and a zero-free region for the Riemann zeta-function". J. Number Theory. 157: 329–349. arXiv:1410.3926. doi:10.1016/J.JNT.2015.05.010. S2CID 117968965.
  18. Polchinski, Joseph (1998). An Introduction to the Bosonic String. String Theory. Vol. I. Cambridge University Press. p. 22. ISBN 978-0-521-63303-1.
  19. Kainz, A. J.; Titulaer, U. M. (1992). "An accurate two-stream moment method for kinetic boundary layer problems of linear kinetic equations". J. Phys. A: Math. Gen. 25 (7): 1855–1874. Bibcode:1992JPhA...25.1855K. doi:10.1088/0305-4470/25/7/026.
  20. Further digits and references for this constant are available at OEISA059750.
  21. Sondow, Jonathan (1998). "An antisymmetric formula for Euler's constant". Mathematics Magazine. 71 (3): 219–220. doi:10.1080/0025570X.1998.11996638. Archived from the original on 4 June 2011. Retrieved 29 May 2006.
  22. Ogilvy, C. S.; Anderson, J. T. (1988). Excursions in Number Theory. Dover Publications. pp. 29–35. ISBN 0-486-25778-9.
  23. Voronin, S. M. (1975). "Theorem on the Universality of the Riemann Zeta Function". Izv. Akad. Nauk SSSR, Ser. Matem. 39: 475–486. Reprinted in Math. USSR Izv. (1975) 9: 443–445.
  24. Ramūnas Garunkštis; Antanas Laurinčikas; Kohji Matsumoto; Jörn Steuding; Rasa Steuding (2010). "Effective uniform approximation by the Riemann zeta-function". Publicacions Matemàtiques. 54 (1): 209–219. doi:10.5565/PUBLMAT_54110_12. JSTOR 43736941.
  25. Bhaskar Bagchi (1982). "A Joint Universality Theorem for Dirichlet L-Functions". Mathematische Zeitschrift. 181 (3): 319–334. doi:10.1007/bf01161980. ISSN 0025-5874. S2CID 120930513.
  26. Steuding, Jörn (2007). Value-Distribution of L-Functions. Lecture Notes in Mathematics. Vol. 1877. Berlin: Springer. p. 19. arXiv:1711.06671. doi:10.1007/978-3-540-44822-8. ISBN 978-3-540-26526-9.
  27. Karatsuba, A. A. (2001). "Lower bounds for the maximum modulus of ζ(s) in small domains of the critical strip". Mat. Zametki. 70 (5): 796–798.
  28. Karatsuba, A. A. (2004). "Lower bounds for the maximum modulus of the Riemann zeta function on short segments of the critical line". Izv. Ross. Akad. Nauk, Ser. Mat. 68 (8): 99–104. Bibcode:2004IzMat..68.1157K. doi:10.1070/IM2004v068n06ABEH000513. S2CID 250796539.
  29. Karatsuba, A. A. (1996). "Density theorem and the behavior of the argument of the Riemann zeta function". Mat. Zametki (60): 448–449.
  30. Karatsuba, A. A. (1996). "On the function S(t)". Izv. Ross. Akad. Nauk, Ser. Mat. 60 (5): 27–56.
  31. ^ Knopp, Konrad (1947). Theory of Functions, Part Two. New York, Dover publications. pp. 51–55.
  32. Riemann, Bernhard (1859). "On the number of primes less than a given magnitude". Monatsberichte der Königlich Preußischen Akademie der Wissenschaften zu Berlin. translated and reprinted in Edwards, H. M. (1974). Riemann's Zeta Function. New York: Academic Press. ISBN 0-12-232750-0. Zbl 0315.10035.
  33. Trivial exceptions of values of s that cause removable singularities are not taken into account throughout this article.
  34. Neukirch, Jürgen (1999). Algebraic number theory. Springer. p. 422. ISBN 3-540-65399-6.
  35. Hashimoto, Yasufumi; Iijima, Yasuyuki; Kurokawa, Nobushige; Wakayama, Masato (2004). "Euler's constants for the Selberg and the Dedekind zeta functions". Bulletin of the Belgian Mathematical Society, Simon Stevin. 11 (4): 493–516. doi:10.36045/bbms/1102689119. MR 2115723.
  36. "A series representation for the Riemann Zeta derived from the Gauss-Kuzmin-Wirsing Operator" (PDF). Linas.org. Retrieved 4 January 2017.
  37. ^ Blagouchine, Iaroslav V. (2018). "Three Notes on Ser's and Hasse's Representations for the Zeta-functions". INTEGERS: The Electronic Journal of Combinatorial Number Theory. 18A: 1–45. arXiv:1606.02044. Bibcode:2016arXiv160602044B.
  38. ^ Hasse, Helmut (1930). "Ein Summierungsverfahren für die Riemannsche ζ-Reihe" [A summation method for the Riemann ζ series]. Mathematische Zeitschrift (in German). 32 (1): 458–464. doi:10.1007/BF01194645. S2CID 120392534.
  39. Sondow, Jonathan (1994). "Analytic continuation of Riemann's zeta function and values at negative integers via Euler's transformation of series" (PDF). Proceedings of the American Mathematical Society. 120 (2): 421–424. doi:10.1090/S0002-9939-1994-1172954-7.
  40. Blagouchine, Iaroslav V. (2016). "Expansions of generalized Euler's constants into the series of polynomials in π and into the formal enveloping series with rational coefficients only". Journal of Number Theory. 158: 365–396. arXiv:1501.00740. doi:10.1016/j.jnt.2015.06.012.
  41. Ser, Joseph (1926). "Sur une expression de la fonction ζ(s) de Riemann" [Upon an expression for Riemann's ζ function]. Comptes rendus hebdomadaires des séances de l'Académie des Sciences (in French). 182: 1075–1077.
  42. Maślanka, Krzysztof (1997). "The Beauty of Nothingness". Acta Cosmologica. XXIII–I: 13–17.
  43. Báez-Duarte, Luis (2010). "On Maslanka's Representation for the Riemann Zeta Function". International Journal of Mathematics and Mathematical Sciences. 2010: 1–9. arXiv:math/0307214. doi:10.1155/2010/714147.
  44. Flajolet, Philippe; Vepstas, Linas (2008). "On Differences of Zeta Values". Journal of Computational and Applied Mathematics. 220 (1-2 October): 58–73. arXiv:math/0611332. Bibcode:2008JCoAM.220...58F. doi:10.1016/j.cam.2007.07.040.
  45. Maślanka, Krzysztof; Koleżyński, Andrzej (2022). "The High Precision Numerical Calculation of Stieltjes Constants. Simple and Fast Algorithm". Computational Methods in Science and Technology. 28 (2): 47–59. arXiv:2210.04609. doi:10.12921/cmst.2022.0000014. S2CID 252780397.
  46. Báez-Duarte, Luis (2003). "A New Necessary and Sufficient Condition for the Riemann Hypothesis". Number Theory. arXiv:math/0307215. Bibcode:2003math......7215B.
  47. Maślanka, Krzysztof (2006). "Báez-Duarte's Criterion for the Riemann Hypothesis and Rice's Integrals". Number Theory. arXiv:math/0603713v2. Bibcode:2006math......3713M.
  48. Wolf, Marek (2014). "Some remarks on the Báez-Duarte criterion for the Riemann Hypothesis". Computational Methods in Science and Technology. 20 (2): 39–47. doi:10.12921/cmst.2014.20.02.39-47.
  49. Borwein, Peter (2000). "An Efficient Algorithm for the Riemann Zeta Function" (PDF). In Théra, Michel A. (ed.). Constructive, Experimental, and Nonlinear Analysis. Conference Proceedings, Canadian Mathematical Society. Vol. 27. Providence, RI: American Mathematical Society, on behalf of the Canadian Mathematical Society. pp. 29–34. ISBN 978-0-8218-2167-1. Archived from the original (PDF) on 26 July 2011. Retrieved 25 November 2017.
  50. Mező, István (2013). "The primorial and the Riemann zeta function". The American Mathematical Monthly. 120 (4): 321.
  51. Komatsu, Takao; Mező, István (2016). "Incomplete poly-Bernoulli numbers associated with incomplete Stirling numbers". Publicationes Mathematicae Debrecen. 88 (3–4): 357–368. arXiv:1510.05799. doi:10.5486/pmd.2016.7361. S2CID 55741906.
  52. "A220335 - OEIS". oeis.org. Retrieved 17 April 2019.
  53. Tóth, László (2022). "Linear Combinations of Dirichlet Series Associated with the Thue-Morse Sequence". Integers. 22 (article 98). arXiv:2211.13570.
  54. Odlyzko, A. M.; Schönhage, A. (1988). "Fast algorithms for multiple evaluations of the Riemann zeta function". Trans. Amer. Math. Soc. 309 (2): 797–809. doi:10.2307/2000939. JSTOR 2000939. MR 0961614..
  55. "Work on spin-chains by A. Knauf, et. al". Empslocal.ex.ac.uk. Retrieved 4 January 2017.
  56. Gene Ward Smith. "Nearest integer to locations of increasingly large peaks of abs(zeta(0.5 + i*2*Pi/log(2)*t)) for increasing real t". The On-Line Encyclopedia of Integer Sequences. Retrieved 4 March 2022.
  57. William A. Sethares (2005). Tuning, Timbre, Spectrum, Scale (2nd ed.). Springer-Verlag London. p. 74. ...there are many different ways to evaluate the goodness, reasonableness, fitness, or quality of a scale...Under some measures, 12-tet is the winner, under others 19-tet appears best, 53-tet often appears among the victors...
  58. Most of the formulas in this section are from § 4 of J. M. Borwein et al. (2000)

Sources

External links

L-functions in number theory
Analytic examples
Algebraic examples
Theorems
Analytic conjectures
Algebraic conjectures
p-adic L-functions
Sequences and series
Integer sequences
Basic
Advanced (list)
Fibonacci spiral with square sizes up to 34.
Properties of sequences
Properties of series
Series
Convergence
Explicit series
Convergent
Divergent
Kinds of series
Hypergeometric series
Bernhard Riemann
Categories: