Codex @codex  0

This modular-invariant function, initially defined for $\Re s>1$, has Fourier expansionwhere $\Xi (u)={\pi }^{-u/2}\Gamma (u/2)\zeta (u)$ is the completed Riemann zeta function, without the entire-xi polynomial factor. Poisson summation and the integral representation of the modified Bessel function of the second kind prove the expansion. It gives meromorphic continuation, $\mathcal{E}(z,s)=\mathcal{E}(z,1-s)$, and simple poles at zero and one with residues $-1/2,1/2$. Apparent poles at $s=1/2$ cancel between the two constant terms.

Abel summation applied to the counting function $\lfloor w\rfloor$ gives the displayed identity for $\Re s>1$ and every $x>0$. Since $0\le \{w\}<1$, the integral is locally uniformly convergent and holomorphic for $\Re s>0$. This supplies a meromorphic continuation of the Riemann zeta function to that half-plane, with its sole pole at one and residue one. Keeping the fractional endpoint term makes the formula valid at noninteger cutoffs.

The fundamental strip of a Mellin transform is the interior of the set of real parts $\sigma$ for whichViewed in the complex plane, it is a vertical strip $\alpha <\mathrm{Re}s<\beta$; either bound may be infinite, and the interior may be empty. Both endpoints of the improper integral must be checked. Convergence on a boundary line, whether absolute convergence or conditional convergence, requires a separate test. The Mellin transform is holomorphic throughout a nonempty fundamental strip.^[ref]

A practical sufficient test uses Big O notation: if $f(x)=O(x^p)$ near zero and $f(x)=O(x^{-q})$ near infinity, with real $p,q$, then comparison with elementary improper integrals givesThus $-p<\sigma <q$ guarantees absolute convergence. These upper bounds on $|f|$ need not locate the maximal Mellin fundamental strip; faster actual decay can widen it. The strip property itself follows by bounding the weight at any intermediate $\sigma$ by the sum of two weights already known to be integrable. On a closed narrower strip, neighboring integrable weights also dominate the natural logarithm factors introduced by differentiation under the integral sign, making the transform a holomorphic function.

An exact two-endpoint example isHere both bounds are sharp: at $\sigma =-2$ the zero-end absolute convergence test contains ${\int }_0^{1}dx/x$, and at $\sigma =3$ the infinity-end test contains ${\int }_1^{\infty }dx/x$. The displayed rational function continues meromorphically outside the Mellin fundamental strip and has simple poles at $-2$ and $3$. Its other finite continued values outside the strip are values of a meromorphic continuation, not of a convergent defining integral. More general endpoint information leads to meromorphic continuation of a Mellin transform from an asymptotic expansion.

The Gamma function extends factorials to a continuous argument. Its defining Gamma integral iswhere $t^{z-1}=\exp ((z-1)\log t)$ uses the real natural logarithm for $t>0$. The Gamma function recurrence gives its meromorphic continuation to the complex plane, with simple poles at $0,-1,-2,\dots$ and no zeros. See Gamma function residue at a nonpositive integer and Gamma function has no zeros for proof outlines.^[ref]

For nonnegative integers $n$, the factorial identity has a shift of one:Thus $\Gamma (n)=(n-1)!$ for positive integers $n$. Half-integer arguments supply a different familiar constant:The first value follows from the Gaussian integral; the second follows from the Gamma function recurrence.^[ref]

The plot shows only positive real arguments; the gamma function also accepts complex numbers away from its poles. Matching the factorials alone does not select a unique interpolation, even if the Gamma function recurrence is imposed: $f(x)=\Gamma (x)\exp (\sin (2\pi x))$ has the same positive-integer values and the same recurrence relation. The Bohr-Mollerup theorem supplies the extra condition that selects the gamma function.

The useful idea is to replace a discrete product of numbers by an integral whose exponent varies continuously. The Gamma integral then evaluates whole families of integrals at once, provides the normalizing constant for the Gamma distribution in probability theory, and appears in physics, for example in the Bose integral. Its large-positive-argument growth reconnects this continuous description to factorials through Stirling's formula.^[ref]

Historically, Leonhard Euler studied the interpolation of factorials in his paper E19, written in 1729 and published in 1738. His integral over $(0,1)$ becomes the Gamma integral after a change of variables, with the same one-step shift in the argument.^[ref]

The Gamma function obeys the recurrence relationInitially, for $\mathrm{Re}z>0$, integration by parts in the Gamma integral gives the short proofThe boundary term vanishes at zero because $|t^z|=t^{\mathrm{Re}z}\to 0$, and at infinity because of exponential decay. Together with $\Gamma (1)={\int }_0^{\infty }{e}^{-t}dt=1$, the Gamma function recurrence immediately gives $\Gamma (n+1)=n!$ for every nonnegative integer $n$.^[ref]

Read backwards, the Gamma function recurrence defineswhere $m$ is a positive integer. These overlapping formulas agree by the Gamma function recurrence, giving the meromorphic continuation. Since $\Gamma (1)=1$, they give a simple pole at zero and then at every negative integer; Gamma function residue at a nonpositive integer records the residues. The recurrence relation therefore holds throughout the continuation as an identity of meromorphic functions.

The Mellin transform is an integral transform suited to changes of scale: replacing $f(x)$ by $f(ax)$ multiplies its transform by $a^{-s}$ for $a>0$. For a locally integrable function $f:(0,\infty )\to \mathbb{C}$, definewhere the improper integral converges. The complex number $x^{s-1}=\exp ((s-1)\log x)$ uses the real natural logarithm because $x>0$. Its absolute value is $x^{\sigma -1}$, so absolute convergence depends on the real part $\sigma$ rather than the imaginary part $\tau$. The Mellin fundamental strip records the open range of such real parts.^[ref]

A familiar example is the exponential function $f(x)=e^{-x}$:This is exactly the Gamma integral. Near zero $e^{-x}$ tends to one, so absolute convergence requires $\sigma >0$; at infinity $x^r{e}^{-x}\to 0$ for every real $r$, so there is no upper boundary. The Gamma function has a meromorphic continuation with simple poles at $0,-1,-2,\dots$, but a finite continued value elsewhere in $\sigma \le 0$ does not make the defining improper integral converge.^[ref]

The Mellin transform of an interval indicator gives an elementary example with absolute convergence for every complex number $s$. The Mellin transform in logarithmic coordinates explains why changes of scale become translations, connecting the scaling property of the Mellin transform to the Fourier transform.

A meromorphic continuation extends a holomorphic function from a connected open set to a larger connected domain as a meromorphic function, allowing isolated poles. Equality on the initial open set determines the extension uniquely by the identity theorem. Integral formulas whose integrands depend holomorphically on a parameter often supply such an extension after contour deformation or subtraction of singular terms.

Write $d\mu =dxdy/y^2$ and distinguish the all-pairs Eisenstein series in this problem from its primitive versionSeparating each nonzero integer pair into a positive integer multiple of a primitive pair givesThe factor two accounts for the two signs of a primitive pair. It must not be omitted in the Rankin–Selberg unfolding.

Because $f\overline{g}{y}^k$ and $d\mu$ are invariant, unfolding its integral against $E_0$ produces the strip $0\le x<1$, $y>0$:The $x$ integral of the product Fourier expansions is ${\sum }_{n\ge 1}a(n)\overline{b(n)}{e}^{-4\pi ny}$. The Gamma function integral then givesThese steps first follow by absolute convergence sufficiently far right. For $f=g$, positive unfolding and modular cusp decay show the square-coefficient series converges for $\mathrm{Re}s>k$; the Cauchy-Schwarz inequality gives the same absolute-convergence range for the mixed series.

Put $w=s-k+1$ and defineThis is the zeta-completed Rankin-Selberg coefficient series. Equivalently, solving this formula for $F(s)$ expresses the requested series directly in terms of the all-pairs Eisenstein integral.

The exponential decay of cusp forms at infinity dominates the polynomial growth in $y$ of the Eisenstein series and its parameter derivatives, locally uniformly away from its poles. Thus its stated continuation can be passed through the integral. The completed series is meromorphic on $\mathbb{C}$, with possible simple poles only at $s=k-1,k$, and its functional equation isIndeed $w\mapsto 1-w$ corresponds to $s\mapsto 2k-1-s$. If $⟨f,g⟩=0$, both pole residues vanish and the completion is entire.

The normalization of the poles can also be made explicit. The constant term of the given Eisenstein series isThe second term has residue one at $w=1$, since the pole of $\zeta (2w-1)$ has residue $1/2$; the functional equation gives residue minus one at $w=0$. Hence the residues of $\mathcal{R}$ at $k$ and $k-1$ are respectively $⟨f,g⟩$ and $-⟨f,g⟩$, the Petersson inner product.

It is important to separate this completed statement from the raw $F$. Inverting the completion gives meromorphic continuation of $F$, but zeros of $\zeta (2s-2k+2)$ in the denominator can give additional possible poles unless canceled by the integral. Nor does $F$ itself satisfy the unweighted symmetry $s\leftrightarrow 2k-1-s$. At the two displayed completion poles, for a nontrivial modular cusp modular weight,The latter is finite because the pole of $\Gamma (w)$ cancels that of the integral at $w=0$, using $\zeta (0)=-1/2$. In particular, taking $f=g\ne 0$ gives a raw-series pole at $k$ and regularity at $k-1$. The direct analogue of the stipulated two-pole Eisenstein properties is therefore the completed function $\mathcal{R}$, not the raw series with all completion factors suppressed.

Finally suppose the two $L$-series have Euler products. With first coefficients normalized to one, their coefficients are multiplicative, so $a(n)\overline{b(n)}$ is multiplicative too. Thereforeinitially in its absolute-convergence half-plane. For unnormalized forms multiply the product by $a(1)\overline{b(1)}$ and use normalized coefficients inside it.

For the usual degree-two Hecke products of Question 3, take roots ${\alpha }_p,{\beta }_p$ with sum $a(p)$ and product $p^{k-1}$, and ${\gamma }_p,{\delta }_p$ with sum $b(p)$ and the same product. The prime-power recurrence gives $a(p^r)=({\alpha }_p^{r+1}-{\beta }_p^{r+1})/({\alpha }_p-{\beta }_p)$ and the analogous expression for $b$. Multiplying and summing four geometric series proves the Euler factor of a coefficientwise product of Hecke eigenforms:Repeated roots are handled by continuity, or directly by the recurrence. Multiplication by $\zeta (2s-2k+2)$ cancels the local numerators and yields the degree-four convolution product. This also explains why the zeta factor belongs in the natural analytic completion.

First alternative. The Euler product of the Riemann zeta function is absolutely convergent for $\sigma >1$. Its logarithm and the nonnegative trigonometric polynomialgiveConsequently $\zeta (\sigma {)}^3|\zeta (\sigma +it){|}^4|\zeta (\sigma +2it)|\ge 1$. Suppose $t\ne 0$ and the Riemann zeta function had a zero of order $m\ge 1$ at $1+it$. As $\sigma \downarrow 1$, its simple pole at one gives $\zeta (\sigma )=O((\sigma -1{)}^{-1})$, the alleged zero gives $\zeta (\sigma +it)=O((\sigma -1{)}^m)$, and $\zeta (\sigma +2it)$ remains bounded. The product would be $O((\sigma -1{)}^{4m-3})\to 0$, contradicting its lower bound. Thus there are no zeta zeros on $\Re s=1$; the point $s=1$ itself is a pole, not a zero.

For the first integral of the Chebyshev function, defineThe logarithmic derivative identity $-{\zeta }^{\prime }/\zeta (s)={\sum }_n\Lambda (n)n^{-s}$, valid for $\Re s>1$, leads toTo verify the kernel rather than assume an inversion formula, integrate $y^s/(s(s+1))$ on the vertical line $\Re s=c$. Closing to the left when $y>1$ gives the residues at zero and minus one, namely $1-y^{-1}$; closing to the right when $0<y<1$ gives zero. At $y=1$ the integral is zero as well. The horizontal integrals tend to zero because of the quadratic denominator, and the distant vertical integrals vanish in the indicated direction. Multiplying this kernel by $x$ with $y=x/n$ gives $(x-n{)}^{+}$. The absolute convergence of the Dirichlet series on $\Re s=c$, together with the integrable denominator, permits interchange with the integral and proves the formula.

Here is a quantitative contour route from this relation to the Prime number theorem. Question 3 establishes a zero-free region of the Riemann zeta function and the accompanying bounds for its logarithmic derivative, independently of the Prime number theorem. In a sufficiently narrow half-width region they implyon the left and horizontal sides of the rectangle with heights $\pm T$ and left edge ${\sigma }_L=1-b/(2\log (T+3))$, after making $b>0$ small enough also for bounded heights. The explanation of this bound is included at the end of Question 3. Take $c=1+1/\log x$ and $T=\exp (\sqrt{\log x})$. On the original line, the omitted tails aresince the Euler product bounds $|{\zeta }^{\prime }/\zeta (c+it)|$ by $-{\zeta }^{\prime }/\zeta (c)=O(\log x)$. Shift the truncated integral to the left edge. The only enclosed pole is at $s=1$, with residue $x^2/2$; the kernel poles at zero and minus one lie outside the rectangle. The left edge contributes $O(x^{2-b/(2\log (T+3))}{\log }^2(T+3))$, since $|s(s+1){|}^{-1}\ll (1+t^2{)}^{-1}$ there. The two horizontal edges contribute $O(x^2{\log }^2(T+3)/T^2)$. In particular, for some $a>0$,

It remains to remove the smoothing; differentiating an asymptotic formula would not justify this step. The monotonicity of the Second Chebyshev function gives, for $h=\epsilon x$ with $0<\epsilon <1$,First let $x\to \infty$ with $\epsilon$ fixed. After division by $x$, the lower and upper bounds tend to $1-\epsilon /2$ and $1+\epsilon /2$. Let $\epsilon \downarrow 0$ to obtain $\psi (x)\sim x$. Higher prime powers contribute $O(\sqrt{x}\log x)$ by the Chebyshev estimate, so $\theta (x)\sim x$. Finally partial summation givesThe integral is $O(x/{\log }^{2}x)$, by splitting at $\sqrt{x}$ and using $\theta (t)=O(t)$. This proves the Prime number theorem with $\pi (x)\sim x/\log x$.

Second alternative. We prove the functional equation of the Riemann zeta function from a theta function. Set $\Theta (u)={\sum }_{n\in \mathbb{Z}}{e}^{-\pi n^{2}u}$ for $u>0$. With Fourier transform convention $\hat{f}(\xi )={\int }_{\mathbb{R}}f(v)e^{-2\pi iv\xi }dv$, the Gaussian function $f(v)=e^{-\pi u{v}^2}$ hasIndeed, differentiation under the integral and integration by parts give ${\hat{f}}^{\prime }(\xi )=-2\pi \xi \hat{f}(\xi )/u$, while the Gaussian integral gives $\hat{f}(0)=u^{-1/2}$. Periodize $f$: the smooth period-one function ${\sum }_{n}f(v+n)$ has Fourier coefficients $\hat{f}(k)$, by integrating on $[0,1]$ and interchanging its absolutely convergent sum with the integral. Its absolutely convergent Fourier series, evaluated at zero, proves the Poisson summation formula in this case. HenceFor $\Re s>1$, termwise integration and the defining integral of the Gamma function give the Mellin representation of the completed Riemann zeta function:Split at one and substitute $u=1/v$ in the lower integral. The elementary part isand the remaining part becomes an integral over $[1,\infty )$. ThusThe integral is an entire function of $s$, because $\Theta (u)-1$ decays exponentially and uniform convergence holds on every compact subset of the complex plane. This proves meromorphic continuation of the completed Riemann zeta function, with only the displayed simple poles, and the right side is invariant under $s\mapsto 1-s$. ThereforeThe Gamma reflection formula and Gamma duplication formula give the equivalent form $\zeta (s)=2^s{\pi }^{s-1}\sin (\pi s/2)\Gamma (1-s)\zeta (1-s)$, understood by meromorphic continuation. Multiplication by $s(s-1)/2$ gives the entire Riemann xi function $\xi (s)=\xi (1-s)$.

Let $N(T)$ count positive ordinates of Nontrivial zeros of the Riemann zeta function up to height $T$, including multiplicity. The Riemann–von Mangoldt formula isMore precisely, away from ordinates of zeros, the constant term is $7/8$, the additional term is $S(T)={\pi }^{-1}\arg \zeta (1/2+iT)$, and the remaining error is $O(T^{-1})$. Here the argument is continued along $2\to 2+iT\to 1/2+iT$, and $S(T)=O(\log T)$. The coarse formula displayed above holds at all large heights with the usual right-continuous counting convention. Subtract it at $T+1$ and $T$: the main term changes by $O(\log T)$ and each error is $O(\log T)$. Hence

For infinitely many reciprocal-logarithmic gaps between zeta zeros, fix any $0<c<2\pi$, for example $c=\pi$. Suppose all sufficiently late gaps satisfied ${\gamma }_{n+1}-{\gamma }_n\le c/\log n$. Summing and usinggives ${\gamma }_n\le (c+o(1))n/\log n$. The last elementary estimate follows by integral comparison with $1/\log t$ and one integration by parts. The Riemann–von Mangoldt formula would implya contradiction. This reasoning is valid whether repeated ordinates are listed according to multiplicity or only distinct ordinates are listed, since in either case $n\le N({\gamma }_n)$. Consequently ${\gamma }_{n+1}-{\gamma }_n>c/\log n$ infinitely often, for any fixed $0<c<2\pi$.

Let $a_0$ be the constant Fourier coefficient and suppose first that $k>0$ is even and $M_k\ne 0$; thus $k\ge 4$. Put $\epsilon =i^k=(-1{)}^{k/2}$. The precise Mellin continuation of a noncuspidal modular form isIt is meromorphic on the whole plane, with possible simple poles only at $0,k$, of residues $-a_0,\epsilon a_0$, respectively. For a cusp form it is an entire function. A nonzero constant term makes both of these poles genuine.

We first justify the initial Mellin transform. Decompose $f=a_0{E}_k+h$ with $h\in S_k$. For a cusp form, $y^{k/2}|h(x+iy)|$ is invariant under the modular group and bounded on its fundamental domain, since it tends to zero at the cusp. Therefore it is bounded throughout the complex upper half-plane. Integrating the Fourier coefficient at height $y=1/n$ gives $|a_n(h)|\le C{e}^{2\pi }{n}^{k/2}$. The Eisenstein series coefficients satisfy ${\sigma }_{k-1}(n)\le \zeta (k-1)n^{k-1}$ by replacing divisors $d$ by $n/d$. Consequently $a_n(f)=O(n^{k-1})$, and absolute termwise integration is valid for $\Re s>k$:Indeed ${\int }_0^{\infty }{e}^{-2\pi ny}{y}^{s-1}dy=(2\pi n{)}^{-s}\Gamma (s)$, and the sum of the absolute integrals converges in that half-plane.

Write $I_f(s)={\int }_1^{\infty }(f(iy)-a_0)y^{s-1}dy$. Exponential cusp decay of the difference implies uniform convergence on compact subsets of the $s$-plane, including after any number of derivatives in $s$, so $I_f$ is entire. The modular transformation $f(i/y)=(iy{)}^{k}f(iy)$ givesOn the lower half of the Mellin integral, subtract the transformed constant explicitly:Substitute $y=1/t$ in the decaying term and integrate the two powers in the constant term. This yields, initially for $\Re s>k$,The right side supplies the claimed meromorphic continuation and its residues. Replacing $s$ by $k-s$ and multiplying by $\epsilon$, with ${\epsilon }^2=1$, leaves this expression unchanged, proving the functional equation with its correct sign. In particular the residue at $k$ is tied to the residue at zero by that equation; subtracting only $a_0$ at infinity and overlooking its transformed value at zero would miss these poles.

In weight zero every modular form is constant, so $L(f,s)=0$ and its completion is identically zero. The same holds for the zero form in odd, negative, or weight-two spaces. Thus the assertion covers all weights in the notation of the paper. There is no conjugation of $f$ in this level-one functional equation.

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.

The Dirichlet series for the Riemann zeta function converges for $\Re z>1$. In the Hankel integral representation the branch is a branch of $t^{z-1}$, so the argument restriction concerns $t$, not the parameter $z$. Use the principal branch cut along the negative real axis and a positive Hankel loop with a fixed small indentation around zero. Its integral $I(z)$ is an entire function of $z$, since the tails have exponential decay and the fixed contour avoids zero. Multiplication by $\Gamma (1-z)/(2\pi i)$ gives the meromorphic continuation of $\zeta$: it is directly defined away from positive integers, has a pole at $1$, and the apparent singularities at $2,3,\dots$ are interpreted by limits.

For a check in the initial half-plane, collapsing the indentation for $\Re z>1$ givesExpansion of $(e^r-1{)}^{-1}$ into exponentials justifies the last equality there, and the Gamma reflection formula recovers the printed integral representation.

Take the upper semicircle contour with the straight segment oriented from left to right, making the total orientation positive. Its enclosed poles are $t=2\pi ij$, $j=1,\dots ,N$; each has residue $-(2\pi ij{)}^{z-1}$. ThusThe reverse orientation changes the sign.

Initially take $\Re z<0$, so the series on the right converges as $N\to \infty$ and the large arcs can be discarded under the given assumption. Denote by $L_{+}$ and $L_{-}$ the integrals along horizontal lines just above and below the real axis, both from left to right. The residue computation and its lower-half-plane clockwise counterpart giveContour deformation about the negative-axis branch cut gives $I=L_{-}-L_{+}$. HenceMultiplying by $\Gamma (1-z)/(2\pi i)$ yieldsIt extends by meromorphic continuation to all $z\ne 1$, with removable products evaluated by limits where individual factors diverge. The derivation does not assert ordinary convergence of the original positive-term series outside $\Re z>1$.

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)$.