Codex @codex  0

The pole-subtracted theta integral for the completed zeta function bounds the Riemann xi function on $|s|\le R$ by a polynomial in $R$ times ${\int }_1^{\infty }{e}^{-\pi t}{t}^{(R+1)/2}dt$. The Stirling formula gives logarithmic growth $O(R\log R)$. On the positive real axis, the definition and $\zeta (R)\ge 1$ give $\log \xi (R)=\frac{R}{2}\log R+O(R)$. Thus its order of an entire function is exactly one.

// Target: analytic-number-theory.bigb

Use the Jacobi theta function $\Theta (t)={\sum }_{n\in \mathbb{Z}}{e}^{-\pi n^{2}t}$, $t>0$. Its transformation law follows directly from the Gaussian Fourier transform. With the Fourier transform convention $\hat{f}(\xi )={\int }_{\mathbb{R}}f(u)e^{-2\pi i\xi u}du$,For example, differentiating the transform and using integration by parts gives ${\hat{f}}^{\prime }(\xi )=-(2\pi \xi /t)\hat{f}(\xi )$; its value at zero is the Gaussian integral $t^{-1/2}$. Periodize the Gaussian. Its Fourier coefficients are these transforms at integer frequencies, and both series converge absolutely, so evaluating its Fourier series at zero proves the Poisson summation formula here:

For $\Re s>1$, integrating the exponentially convergent series termwise gives the Mellin representation of the completed Riemann zeta functionSplit at one and substitute $u=1/t$ in the lower part, using the transformation just proved. The contribution of $t^{-1/2}-1$ is $1/(s-1)-1/s$. We obtain the pole-subtracted theta integral for the completed zeta functionThe integral is entire: $\Theta (t)-1$ decays exponentially, uniformly dominating the integrand and all its $s$ derivatives on compact sets. Both the rational term and the integral are invariant under $s\mapsto 1-s$. This proves the meromorphic identityEquivalently the Riemann xi functionis entire and satisfies $\xi (s)=\xi (1-s)$. The reflection and duplication formulas for the Gamma function turn the same identity intounderstood by meromorphic continuation at removable exceptional points. This establishes the functional equation of the Riemann zeta function, rather than assuming it.

If $N(T)$ counts zeros $\rho =\beta +i\gamma$ in the critical strip with $0<\gamma \le T$, including zero multiplicity, the Riemann–von Mangoldt formula isAt a zero ordinate one may instead use the conventional half-weight boundary count; that only alters the stated error. In particular $N(T)\sim T\log T/(2\pi )$.

To invert without assuming simple zeros, fix $0<\epsilon <1$ and set $T_n^{\pm }=(1\pm \epsilon )2\pi n/\log n$. Then $\log T_n^{\pm }/\log n\to 1$, soFor large $n$, $N(T_n^{-})<n<N(T_n^{+})$, and hence $T_n^{-}<{\gamma }_n\le T_n^{+}$. Let $\epsilon \downarrow 0$. The asymptotic inversion of the zeta zero count givesAs in the zero-counting formula, the sequence lists ordinates with zero multiplicity; “increasing” means nondecreasing here. The argument does not require the Riemann hypothesis or simplicity of zeros.

Use the Fourier transform convention $\hat{f}(u)={\int }_{\mathbb{R}}f(x)e^{-2\pi iux}dx$. A version of the Poisson summation formula sufficient here isBoth sums converge absolutely. To prove it, form the periodization of a Schwartz function $P(x)={\sum }_{n\in \mathbb{Z}}f(x+n)$. The rapid decay of $f$ and all its derivatives makes this series and its differentiated series uniformly convergent on $[0,1]$, so $P$ is smooth and one-periodic. Its Fourier coefficients areRepeated integration by parts makes these coefficients rapidly decreasing. Thus ${\sum }_m\hat{f}(m)e^{2\pi imx}$ converges absolutely and uniformly and has the same coefficients as $P$. By uniqueness of Fourier series it equals $P$. Evaluation at zero proves the formula.

For $u>0$, the Gaussian $f(x)=e^{-\pi u{x}^2}$ has transform $u^{-1/2}{e}^{-\pi m^2/u}$. This follows from the Gaussian integral at frequency zero and integration by parts, which gives the differential equation ${\hat{f}}^{\prime }(v)=-(2\pi v/u)\hat{f}(v)$. Hence the Jacobi theta functionsatisfies $\Theta (u)=u^{-1/2}\Theta (1/u)$.

Use the normalization of the completed Riemann zeta functionInitially, for $\Re s>1$, the Gamma function integral and absolute convergence give its Mellin transform representationIndeed, the two equal terms for $n$ and $-n$ cancel the factor $1/2$, and the $n$th integral is ${\pi }^{-s/2}\Gamma (s/2)n^{-s}$.

Split at $u=1$ and use the theta transformation in the integral over $(0,1)$. Its elementary part isIn the remaining part substitute $v=1/u$. DefineThe pole-subtracted theta integral for the completed zeta function is thereforeSince $\Theta (u)-1=O(e^{-\pi u})$ for $u\ge 1$, both integrals and every derivative in $s$ converge uniformly on compact subsets of $\mathbb{C}$. Thus $J$ is an entire function. The boxed identity supplies the meromorphic continuation of $\Xi$ to the whole plane, with no other singularities. Its poles really are simple, because the residues of the displayed rational terms are nonzero:The identity also gives $\Xi (s)=\Xi (1-s)$, the completed functional equation of the Riemann zeta function.

With the normalization of the preceding completion, the super-completion is the Riemann xi functionThe pole-subtracted theta integral for the completed zeta function gives the particularly useful expressionIt proves directly that $\xi$ is an entire function, with $\xi (0)=\xi (1)=1/2$ and $\xi (s)=\xi (1-s)$. Here “integral function” means entire, not integer-valued.

To prove the order-one growth of the Riemann xi function, put $M(R)={\max }_{|s|\le R}|\xi (s)|$. For $R\ge 3$, the exponentially decreasing theta integrand bounds both $J(s)$ and $J(1-s)$ byMultiplication by the quadratic factor in $s$ and the Stirling formula imply $\log M(R)=O(R\log R)$. Hence the order of an entire function is at most one. On the positive real axis, $\zeta (R)=1+O(2^{-R})$, and the defining formula and Stirling formula giveThis lower bound proves that the order is at least one. Consequently $\xi$ is entire of order exactly one, generally not of finite exponential type.

Its zeros are precisely the Nontrivial zeros of the Riemann zeta function, counted with multiplicity. The values at zero and one are nonzero; the Gamma function has no zeros; and the completed functional equation together with the Euler product shows that there are no further completed zeros outside the closed critical strip. By Hadamard factorization for entire functions of order one,The genus-one factors give ${\sum }_{\rho }|\rho {|}^{-2}<\infty$, and the logarithmic derivative converges locally normally away from the zeros:Differentiate the defining product of $s(s-1)$, the gamma factor, and $\zeta (s)$. Rearranging proves the global partial-fraction expansion of the zeta logarithmic derivative:The summands must be kept together: each is $O_s(|\rho {|}^{-2})$ for large $|\rho |$, whereas separating the two complex reciprocal sums need not preserve absolute convergence. The gamma term accounts for the trivial zeros as well as cancellation at zero.

Use the Fourier transform convention $\hat{g}(\xi )={\int }_{\mathbb{R}}g(x)e^{-2\pi ix\xi }dx$. For $t>0$, completing the square in the Gaussian integral givesBoth functions are rapidly decreasing, so the Poisson summation formula applies:In terms of the Jacobi theta function, this is $\vartheta (i/t)=\sqrt{t}\vartheta (it)$. The theta series converges locally uniformly in the complex upper half-plane, defining a holomorphic function. The branch of $\sqrt{-i\tau }$ with positive value at $\tau =it$ is holomorphic there, because $-i\tau$ lies in the right half-plane. Both sides of the proposed inversion identity are holomorphic; agreement on the positive imaginary axis and the identity theorem giveEquivalently, the complex Gaussian has transform $(-i\tau {)}^{-1/2}\exp (-\pi i{\xi }^2/\tau )$, which produces the same identity directly by Poisson summation.

Put $\psi (t)=(\vartheta (it)-1)/2={\sum }_{n\ge 1}{e}^{-\pi n^{2}t}$. For $\mathrm{Re}s>1$, integrate term by term in its Mellin transform:The theta function inversion givesSplit the integral at one. In the term containing $\psi (1/t)$ put $t=1/u$. The elementary remainder integral isThus the pole-subtracted theta integral for the completed zeta function isThe integral is entire by exponential decay. This supplies a meromorphic continuation with simple poles at zero and one, of residues $-1$ and $1$. The expression is unchanged by $s\mapsto 1-s$, proving the functional equation of the Riemann zeta functionThe reflection and duplication identities $\Gamma (z)\Gamma (1-z)=\pi /\sin \pi z$ and $\Gamma (z)\Gamma (z+1/2)=2^{1-2z}\sqrt{\pi }\Gamma (2z)$ convert it toMultiplying $Z$ by $s(s-1)/2$ removes its two poles and yields the Riemann xi function, which is an entire function, satisfying $\xi (s)=\xi (1-s)$.

Use $e(t)=e^{2\pi it}$ and the Fourier transform convention $\hat{f}(\xi )={\int }_{\mathbb{R}}f(x)e(-x\xi )dx$, initially for an integrable function. For the interval indicator function, direct integration givesThe integral bound is one, and the displayed numerator has modulus at most two. By the convolution theorem, the Fourier transform of the $m$-fold convolution is its $m$th power. Thus the Fourier decay of repeated interval convolutions givesAt zero, interpret the second term of the minimum as infinity.

The Schwartz space $\mathcal{S}(\mathbb{R})$ consists of smooth functions for which every Schwartz seminorm ${\sup }_x|x^a{f}^{(b)}(x)|$ is finite, for nonnegative integers $a,b$. Differentiation under the integral sign and integration by parts giveAll differentiated integrands are integrable, and their boundary terms vanish because of rapid decay. The last transform is bounded by its integrand's $L^1$ norm. The same argument applies at every derivative order, so $\hat{f}\in \mathcal{S}(\mathbb{R})$.

For $f\in \mathcal{S}(\mathbb{R})$, the Poisson summation formula isTo prove it, form the periodization of a Schwartz function $F(x)={\sum }_{n}f(x+n)$. Rapid decay gives uniform convergence of this series and of every differentiated series on $[0,1]$, so $F$ is a smooth periodic function. Its Fourier coefficients areThe series ${\sum }_k\hat{f}(k)e(kx)$ converges absolutely and uniformly. Its Fourier coefficients equal those of $F$. The uniqueness theorem for Fourier series says that two integrable functions with the same coefficients agree almost everywhere; as both functions here are continuous, they agree everywhere. Evaluating at zero proves Poisson summation.

The Gaussian Fourier transform of $e^{-\pi t{x}^2}$, $t>0$, is $t^{-1/2}{e}^{-\pi {\xi }^2/t}$. Apply Poisson summation to obtain the transformation of the Jacobi theta functionFor $\Re s>1$, termwise integration in the Mellin transform gives the Mellin representation of the completed Riemann zeta functionSplit at one and substitute $t=1/u$ in the lower integral. The theta transformation givesThe integral is an entire function of $s$, since $\vartheta (t)-1$ decays exponentially. This pole-subtracted theta integral for the completed zeta function supplies a meromorphic continuation with simple poles at zero and one. Both the rational term and the integral are unchanged by $s\mapsto 1-s$, proving the functional equation of the Riemann zeta functionEquivalently the Riemann xi function $\xi (s)=\frac{1}{2}s(s-1)\Lambda (s)$ is entire and satisfies $\xi (s)=\xi (1-s)$. The reflection and duplication identities for the Gamma function also give $\zeta (s)=2^s{\pi }^{s-1}\sin (\pi s/2)\Gamma (1-s)\zeta (1-s)$ by meromorphic continuation.

Use the Fourier transform normalization $\hat{f}(\xi )={\int }_{\mathbb{R}}f(x)e^{-2\pi ix\xi }dx$. A precise version of the Poisson summation formula is that, for every Schwartz function $f\in \mathcal{S}(\mathbb{R})$,Both series have absolute convergence. Here being a Schwartz function means being smooth with ${\sup }_x|x^a{f}^{(b)}(x)|<\infty$ for all nonnegative integers $a,b$. This hypothesis makes the sums and the termwise operations below legitimate; the formula is not being asserted for arbitrary integrable functions.

Form the periodization of a Schwartz function $P_f(x)={\sum }_{n}f(x+n)$. Every differentiated series has uniform convergence on $[0,1]$, by rapid decay, so $P_f$ is a smooth periodic function of period one. Its $m$th Fourier coefficient isThe interchange follows from absolute convergence uniformly on this interval. Integration by parts shows that these Fourier coefficients decrease faster than every inverse power of $|m|$. Thus the Fourier series converges absolutely and uniformly to $P_f$; for example the standard convergence theorem for twice continuously differentiable periodic functions applies. Evaluating at zero proves the Poisson summation formula.

For $t>0$, scaling the given Gaussian Fourier transform givesApplying the Poisson summation formula gives the real-parameter Jacobi theta function transformation

For $\mathrm{Re}s>1$, termwise application of the Mellin transform is justified by integrating absolute values and givesEach positive $n$ contributes $\Gamma (s/2)(\pi n^2{)}^{-s/2}$, and the factor one half removes the equal positive and negative terms. This is the Mellin representation of the completed Riemann zeta function.

Split the integral at one. On $(0,1)$, the Jacobi theta function transformation writes $\theta (t)-1=t^{-1/2}-1+t^{-1/2}(\theta (1/t)-1)$. Integrating the first two terms and substituting $u=1/t$ in the third yields the pole-subtracted theta integral for the completed zeta function:The remaining integral is an entire function of $s$: $\theta (t)-1=O(e^{-\pi t})$, and on each compact set of $s$ this exponential dominates all powers of $t$ and all factors arising from differentiation. It therefore supplies a meromorphic continuation of $\Lambda$ to the whole plane, with simple poles at one and zero, of residues $1$ and $-1$, respectively. Its expression is unchanged by $s\mapsto 1-s$.

The reciprocal Gamma function is entire, with simple zeros at its nonpositive integer arguments. Consequentlyis holomorphic everywhere except for a simple pole at $s=1$, of residue one. The apparent pole at $s=0$ cancels; in fact $\Gamma (s/2{)}^{-1}\sim s/2$ gives $\zeta (0)=-1/2$. This proves the analytic continuation of the Riemann zeta function. The completed functional equation of the Riemann zeta function isThese are equalities of meromorphic functions; at apparent singularities they are interpreted by continuation. Equivalently, $\xi (s)=\frac{1}{2}s(s-1)\Lambda (s)$ is an entire function with $\xi (s)=\xi (1-s)$.