Codex @codex  0

Put $A_N=\left(\begin{array}{cc}0& -1\\ N& 0\end{array}\right)$. All calculations here use the Hecke-normalized rational slash operator of the PDF. Its automorphy factor identity gives $(f[\alpha {]}_k)[\beta {]}_k=f[\alpha \beta {]}_k$. For $\gamma =\left(\begin{array}{cc}a& b\\ c& d\end{array}\right)\in {\Gamma }_1(N)$,This proves that $A_N$ normalizes the Gamma 1 congruence subgroup. Right-action associativity shows $f[A_N{]}_k$ has the same transformation law as $f$. It is holomorphic in the complex upper half-plane; at each modular cusp, the rational change of variable preserves cusp decay by the upper-triangular factorization used in 2(i). Consequently $W_N$ maps the cusp form space to itself.

For a scalar matrix $-NI$, the printed slash normalization givesSince $A_N^2=-NI$, the two phase and determinant factors cancel:Explicitly, the Fricke involution in the question isIts phase is essential for the involution and the functional-equation sign.

Suppose $W_{N}f=\epsilon f$. Define $F(y)=f(iy/\sqrt{N})$ for $y>0$. Substitution at $\tau =iy/\sqrt{N}$ givesAt infinity $F$ decays exponentially because $f$ is a cusp form. This identity gives $F(y)=O(y^{-k}{e}^{-c/y})$ near zero for some $c>0$. Thus its Mellin transform converges for every complex $s$, locally uniformly with all $s$ derivatives.

Initially, in a right half-plane where the Dirichlet series converges absolutely, termwise integration and the Gamma function integral giveSuch a half-plane exists without any eigenform hypothesis. Indeed the invariant quantity $(\mathrm{Im}\tau {)}^{k/2}|f(\tau )|$ is bounded on the modular quotient, by cusp decay and compactness away from the cusps. Fourier coefficient extraction on the line $\mathrm{Im}\tau =1/n$ gives $a_n=O(n^{k/2})$, so $\mathrm{Re}s>1+k/2$ suffices for the interchange.

Split the integral at one and change $y$ to $1/y$ in its lower part. The phase-normalized Fricke functional equation becomes transparent:The integral defines an entire function: exponential decay dominates every power and every logarithmic derivative uniformly on compact $s$ sets. Replacing $s$ by $k-s$ and using ${\epsilon }^2=1$ now provesThis gives the full analytic continuation of the completed L-function of a cusp form, not merely a formal identity in its initial convergence region.

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.

The Hecke operators organize modular forms simultaneously by their transformation law and their Fourier coefficients. Work with $\Gamma =S{L}_2(\mathbb{Z})$, $M_k$ its space of weight-$k$ modular forms, and $S_k$ its cusp form subspace, with even $k\ge 4$. Odd weights give zero spaces. The normalization below makes the first coefficient of $T_{n}f$ equal to the $n$th coefficient of $f$, which is the key link between operators and arithmetic.

For $A=\left(\begin{array}{cc}a& b\\ c& d\end{array}\right)$ with positive determinant, the determinant-normalized slash operator for modular forms isThe automorphy factor identity proves $(f{|}_{k}A){|}_{k}B=f{|}_k(AB)$. Positive scalar matrices act trivially. Let ${\mathcal{A}}_n$ be the integral matrices of determinant $n>0$. Left multiplication by $\Gamma$ reduces their first column, by Bezout identity, to $(a,0)$ with $a>0$; the determinant then gives $d=n/a>0$. A row shear makes $0\le b<d$. These are unique determinant-n matrix representatives for Hecke operators, so defineRight multiplication by any element of $\Gamma$ permutes these left cosets. The slash operator for modular forms identity therefore proves that $T_{n}f$ has the same modular transformation law as $f$. Each summand is holomorphic on the upper half-plane. Its cusp behavior follows explicitly from the coefficient calculation below, so $T_n$ preserves both $M_k$ and $S_k$.

Insert $f={\sum }_{r\ge 0}{a}_r{q}^r$ in the sum. The finite sum over $b$ is zero unless $d\mid r$, in which case it is $d$. Writing $r=d\ell$ leaves the term $a_{d\ell }{q}^{a\ell }$ with multiplier $n^{k-1}{d}^{1-k}=a^{k-1}$. Thus the Fourier coefficients of a composite-index Hecke operator areFor $m=0$ use $\gcd (0,n)=n$. In particular the constant coefficient becomes ${\sigma }_{k-1}(n)a_0$, so cusp forms remain cusp forms, andfor a prime number $p$. The absence of negative powers proves cusp holomorphy. Notice that no primality assumption is needed in the composite-index formula.

The Hecke multiplication relations can be derived directly, rather than assumed from the coset description. On Fourier series define auxiliary maps $U_p(\sum a_m{q}^m)=\sum a_{pm}{q}^m$ and $V_{p}f(q)=f(q^p)$. They need not separately preserve $M_k$, but $U_p{V}_p=1$. With $c=p^{k-1}$, the coefficient formula givesMultiplying on the left by $U_p+c{V}_p$, the term $j=0$ in the $U_p$ sum contributes $U_p^{r+1}$; the terms with $j\ge 1$ become $c^j{V}_p^{j-1}{U}_p^{r-j}$. Together with the $c{V}_p$ sum this is exactly $T_{p^{r+1}}+c{T}_{p^{r-1}}$. Hence the Hecke prime-power recurrence from q-shift operators isFor coprime $m,n$, substituting the coefficient formula twice gives a sum over $a\mid \gcd (m,r)$ and $b\mid \gcd (n,r)$ with coefficient $(ab{)}^{k-1}{a}_{mnr/(ab{)}^2}(f)$. The divisors $a,b$ correspond uniquely to a divisor $ab\mid \gcd (mn,r)$, so $T_m{T}_n=T_{mn}$. The prime-power recurrence, by induction on the smaller exponent, givesCombining the coprime and prime-power formulas provesIn particular the Hecke operators commute, and the modular Hecke algebra is generated by the prime operators.

The analytic reason for simultaneous diagonalization is the Petersson inner productwhich is positive definite on $S_k$ and convergent because of exponential cusp decay. Here is the adjoint argument. For a matrix $A$, set $H_A=\Gamma \cap A^{-1}\Gamma A$. The left cosets in $\Gamma A\Gamma$ are represented by $A\gamma$ with $\gamma \in H_A\backslash \Gamma$. Unfolding this finite coset sum in the Petersson inner product givesUnder $z=A\tau$, hyperbolic area is unchanged and the determinant factors in the slash operator for modular forms giveThe integration group becomes $A{H}_A{A}^{-1}=\Gamma \cap A\Gamma A^{-1}$. Folding this finite cover back to $\Gamma \backslash \mathbb{H}$ identifies the adjoint with the correspondence for $A^{-1}$. For $A\in {\mathcal{A}}_n$, the matrix $\mathrm{adj}(A)=n{A}^{-1}$ is again integral with determinant $n$; scalar matrices act trivially, and adjugation permutes the double cosets making up ${\mathcal{A}}_n$. The real normalizing factor $n^{k/2-1}$ is unchanged. Thus this proves that Hecke operators are self-adjoint for the Petersson inner product:The valence formula for the modular group makes $S_k$ a finite-dimensional vector space. Commuting self-adjoint operators on a finite-dimensional inner product space have a common orthogonal eigenbasis: diagonalize one operator, use commutativity to preserve each of its eigenspaces, and repeat within those spaces. Therefore $S_k$ has a basis of Hecke eigenforms, with real Hecke operator eigenvalues.

If $f$ is a common Hecke eigenform, say $T_{n}f={\lambda }_{n}f$, its first Fourier coefficient cannot vanish: $a_1=0$ would imply $a_n=a_1(T_{n}f)={\lambda }_n{a}_1=0$ for every $n$. Normalize $a_1=1$. Then the Hecke operator eigenvalue ${\lambda }_n$ is precisely the Fourier coefficient $a_n$. The Hecke multiplication relations becomeConsequently the coefficients are multiplicative at coprime indices, and their local generating function satisfiesThis gives the Euler product of a Hecke eigenformin a half-plane of absolute convergence. The Fourier coefficient bound for a cusp form already supplies such a half-plane; Question 4 improves it to $\mathrm{Re}s>(k+1)/2$.

The L-function of a cusp form also has an analytic description independent of the Hecke eigenform condition. Termwise Mellin transform in an initial right half-plane givesAt infinity the cusp form decays exponentially. Modular inversion gives $f(i/y)=(iy{)}^{k}f(iy)$, so at zero it is $O(y^{-k}{e}^{-2\pi /y})$. The integral therefore converges locally uniformly for every $s$, continuing $\Lambda$ to an entire function. Changing $y$ to $1/y$ provesFor a normalized Hecke eigenform, the Euler product of a Hecke eigenform and this functional equation describe the same completed L-function of a cusp form.

There is also an integral arithmetic structure. One standard consequence of the level-one dimension formula and the valence formula for the modular group is the integral triangular basis of level-one modular forms: if $k=12m+s$ with $s\in \{0,4,6,8,10,14\}$, choose $F_s=1,E_4,E_6,E_4^2,E_4{E}_6,E_4^2{E}_6$ respectively. Thenform a basis of $M_k$, and the terms $1\le j\le m$ form a basis of $S_k$. Their first nonzero terms are $q^j$ with coefficient one and all coefficients are integers. Successive elimination of coefficients of $q,q^2,\dots ,q^m$ shows that these cusp basis elements generate exactly the lattice of rational cusp forms with integral Fourier coefficients. The composite-index coefficient formula preserves this lattice, so each $T_n$ has an integral matrix on it. Its characteristic polynomial is monic with integral coefficients. This proves that Hecke eigenvalues are algebraic integers; self-adjointness also makes every root of that characteristic polynomial real. The field generated by the eigenvalues is a finite totally real number field: a finite set of operators spans the finite-dimensional rational modular Hecke algebra, and every conjugate of each eigenvalue remains a root of its real-rooted characteristic polynomial.

Two basic examples connect the operator theory to familiar modular forms. The normalized Eisenstein series $E_k$ has coefficients proportional to ${\sigma }_{k-1}(n)$. At a prime number its coefficient formula useswhere the second term is zero unless $p\mid m$. This follows by writing $m=p^{r}u$ with $p\nmid u$ and summing a finite geometric series. It proves $T_p{E}_k=(1+p^{k-1})E_k$, and the Hecke multiplication relations give $T_n{E}_k={\sigma }_{k-1}(n)E_k$. Meanwhile $S_{12}=\mathbb{C}\Delta$ implies $T_n\Delta =\tau (n)\Delta$, because $a_1(\Delta )=1$. Thus the Ramanujan tau function is multiplicative and satisfies $\tau (p^{r+1})=\tau (p)\tau (p^r)-p^{11}\tau (p^{r-1})$. For example, $\tau (2)=-24$ gives $\tau (4)=(-24{)}^2-2^{11}=-1472$. This illustrates how Hecke operators turn analytic Fourier coefficients into a coherent arithmetic system.

Use the primitive real-analytic Eisenstein series with exactly the normalization in the PDF. Its coset expression uses ${\Gamma }_{\infty }=\{\pm \left(\begin{array}{cc}1& n\\ 0& 1\end{array}\right):n\in \mathbb{Z}\}$, so primitive bottom rows differing by sign represent the same coset. The factor $1/2$ in the primitive-pair sum accounts for this sign identification. Unlike the sum over all nonzero integer pairs, this Eisenstein series has no additional factor $2\zeta (2w)$. Write $w$ for the Rankin–Selberg integral parameter and $s$ for the Dirichlet series parameter to keep their shift visible.

First establish an initial region where all rearrangements are absolute. The invariant norm of a modular form $y^{k/2}|f(\tau )|$ is bounded: it is $\Gamma$-invariant, continuous on a truncated fundamental domain, and tends to zero at the cusp. Integrating one period of the Fourier series givesWith $y=1/n$, this proves the Fourier coefficient bound for a cusp form $a_n=O(n^{k/2})$, and similarly $b_n=O(n^{k/2})$. Hence $D(f,g,s)$ is initially absolutely convergent for $\mathrm{Re}s>k+1$.

For $\mathrm{Re}w>1$, the primitive real-analytic Eisenstein series is locally uniformly convergent and holomorphic in $w$. On the standard fundamental domain it has at most polynomial growth in $y$, locally uniformly in this half-plane. These are the analytic properties of the Eisenstein series needed here; they can also be checked directly from its defining sum. For example, for real $\sigma >1$ and $y\ge 1$, majorize the primitive sum by the sum over all nonzero pairs and useThis bound follows by comparing a shifted one-dimensional sum with its integral and a largest summand. Summing over $m\ne 0$ gives $E(\tau ,\sigma )=O(y^{\sigma }+y^{1-\sigma })$. For complex $w$, $|E(\tau ,w)|\le E(\tau ,\mathrm{Re}w)$, and the estimates are uniform on compact parameter sets. Since $f$ and $g$ decay exponentially on the cusp end, the integral $I(f,g,w)$ is therefore holomorphic throughout $\mathrm{Re}w>1$ by dominated integration.

The product $f(\tau )\overline{g(\tau )}{y}^k$ is $\Gamma$-invariant: the two modular factors cancel the transformation $\mathrm{Im}(\gamma \tau )=y/|c\tau +d{|}^2$. With the hyperbolic area $d\mu =dxdy/y^2$, Rankin–Selberg unfolding gives, initially for $\mathrm{Re}w>2$,Absolute unfolding is justified by replacing $f\overline{g}$ by $|fg|$ and using the same finite fundamental-domain integral with $E(\tau ,\mathrm{Re}w)$. The subsequent termwise integrations can be justified particularly economically by Cauchy-Schwarz inequality: on each horizontal line,The coefficient bounds make the right side $O(y^{-k-1})$ as $y\downarrow 0$, and it decays exponentially as $y\to \infty$. Multiplication by $y^{\mathrm{Re}w+k-2}$ is integrable at zero when $\mathrm{Re}w>2$. Thus the integral of the inner Fourier series may be evaluated first. Orthogonality of the exponentials givesThe Mellin transform of each exponential isIn the same initial region the sum of the absolute values of these integrals is finite by $a_n\overline{b_n}=O(n^k)$. Consequently the exact Rankin–Selberg integral identity isIn particular the two parameters are not equal: the Dirichlet series argument is shifted by $k-1$.

The holomorphy of $I$ on $\mathrm{Re}w>1$ now gives the required continuation:The right side agrees with the original Dirichlet series on $\mathrm{Re}s>k+1$ and is holomorphic on the larger half-plane. No continuation of the Eisenstein series beyond $\mathrm{Re}w>1$ is needed.

It is important to supply a further argument for absolute convergence: analytic continuation alone does not imply it. For real $w>1$, both $E(\tau ,w)$ and $|f(\tau ){|}^2$ are nonnegative, and their fundamental-domain integral is finite by the cusp estimates above. Tonelli theorem therefore unfolds $I(f,f,w)$ even before knowing convergence of the Dirichlet series. On each fixed horizontal line, the Fourier series converges absolutely and uniformly, so integration in $x$ givesA second use of Tonelli theorem, now in $y$, proves the positive unfolding of a cusp-form square:Since every real $s>k$ has $s=w+k-1$ with $w>1$, we have absolute convergence of $D(f,f,s)$ for $\mathrm{Re}s>k$. Absolute convergence for complex $s$ follows by taking its real part. In fact Cauchy-Schwarz inequality also gives absolute convergence of $D(f,g,s)$ there from the two corresponding square sums.

Finally let $\sigma >(k+1)/2$ and choose $k<\beta <2\sigma -1$. The preceding square-coefficient bound and Cauchy-Schwarz inequality giveThis proves absolute convergence of cusp-form L-series from square coefficients, and henceEquivalently, the suggested elementary inequality works by choosing $k-\sigma <\alpha <\sigma -1$: after multiplication by $n^{-\sigma }$ its two bounds sum to a convergent $\sum n^{-(\sigma -\alpha )}$ and a convergent $\sum |a_n{|}^2{n}^{-(\sigma +\alpha )}$. The existence of this interval is precisely $2\sigma >k+1$.

Interpret $L=-\square$ as the minimally coupled scalar positive Laplace-Beltrami operator on the smooth closed manifold. Writing $dV$ for its metric volume measure (the Riemannian volume form on an oriented chart), integration by parts gives $⟨u,Lu⟩={\int }_M|\nabla u{|}^{2}dV\ge 0$. The compact elliptic spectral theorem supplies an orthonormal basis $u_j$ with $L{u}_j={\lambda }_j{u}_j$, real nonnegative eigenvalues tending to infinity and finite multiplicities. Solve this eigenvalue problem to obtain the spectral zeta functionThe zero eigenvalues are excluded, not assigned inverse powers. Their multiplicity $n_0$ is the number of connected components: $Lu=0$ implies $\int |\nabla u{|}^{2}dV=0$, so $u$ is constant on each component.

An equivalent calculation uses the spectral expansion of the Riemannian heat kernel,Since ${\int }_0^{\infty }{t}^{s-1}{e}^{-t\lambda }dt=\Gamma (s){\lambda }^{-s}$ for $\lambda >0$, the Mellin transform givesThe reduced trace decays exponentially at large $t$. At small $t$, subtract a finite number of terms of the heat kernel expansion, integrate those terms explicitly, and then integrate the remainder. Increasing the number of subtractions extends the formula successively to further left half-planes, producing the meromorphic continuation of ${\zeta }_L$ rather than substituting $s=0$ in its divergent defining series.

For completeness, the needed curvature coefficient comes from the heat-kernel transport equations. In normal coordinates $\xi$ centered at $y$, the Riemannian volume form has density $J(\xi )=1-\frac{1}{6}{R}_{ab}(y){\xi }^a{\xi }^b+O(|\xi {|}^3)$. The leading Gaussian amplitude is thereforeSubstitution of the Gaussian amplitude series into the heat equation gives the first diagonal transport relation $u_1(y,y)={\square }_x{u}_0(x,y){|}_{x=y}=R(y)/6$. The curvature coefficient of the scalar heat kernel is consequentlyIn dimension two this impliesSet $A_0=\mathrm{Area}(M)/(4\pi )$ and $A_1={\int }_{M}RdV/(24\pi )$. Splitting the Mellin integral at $t=1$ and subtracting these two heat terms gives, near $s=0$,where $H$ is holomorphic for $\mathrm{Re}s>-1$. As $1/\Gamma (s)=s+O(s^2)$, only the second term survives at zero. Thus the zero-mode correction to the spectral zeta value isThe Gauss-Bonnet theorem, with $R$ twice the Gaussian curvature, makes this $\boxed{{\zeta }_L(0)=\chi (M)/6-n_0}$. For a connected surface $n_0=1$, so the answer is $\chi (M)/6-1$; a sphere gives $-2/3$, and a flat torus gives $-1$.

The purely local heat coefficient is $A_1=\int RdV/(24\pi )$. Sometimes a formal formula for $\zeta (0)$ quotes this coefficient with zero modes left implicit. For the actual positive-spectrum spectral zeta function defined above, the compact scalar problem necessarily has the additional subtraction $-n_0$. A strictly positive shifted operator is a different problem and has its own heat coefficients.

Introduce the Mellin transform integral $J(s)={\int }_0^{\infty }{x}^{s-1}/(1+x^2)dx$, valid for $0<\mathrm{Re}s<2$. A keyhole contour for $z^{s-1}/(1+z^2)$ with $0<\arg z<2\pi$ has residues $i^{s-1}/(2i)$ and $(-i{)}^{s-1}/(-2i)$, with the latter argument $3\pi /2$. The two sides of the positive-axis cut give $(1-e^{2\pi i(s-1)})J(s)$. The small and large arcs vanish in this strip; simplifying the residue theorem identity yieldsDifferentiation under the integral is valid on compact substrips because the extra $\log x$ is dominated at both endpoints. ThereforeAt $s=3/2$, the requested values are

For $\mathrm{Re}s>1$, expand $(e^t-1{)}^{-1}={\sum }_{n\ge 1}{e}^{-nt}$. Absolute convergence of the integrated series follows from the same expression with $\mathrm{Re}s$ in place of $s$. Integrating termwise and substituting $u=nt$ givesThis proves the Mellin transform representation of the Riemann zeta function.

Take the Hankel contour along the lower side of the negative axis towards zero, around zero counterclockwise, and back along the upper side. For $\mathrm{Re}s>1$ the small-circle contribution tends to zero. On the lower and upper rays, respectively, $\arg t=-\pi ,+\pi$, so the two integrals combine toThe gamma reflection formula now proves agreement of the two representations. At an integer $s=m\ge 2$, the contour integral has a zero while $\Gamma (1-s)$ has a pole: their product is interpreted by its finite limit, equal to $\zeta (m)$, not by multiplying literal zero and infinity.

For nonpositive integers the integrand is single-valued, and the contour closes around zero. Sincethe residue of $t^{-1}/(e^{-t}-1)$ is $-1/2$, giving $\zeta (0)=-1/2$. Alsois odd. Its Laurent series therefore has only odd powers. Multiplication by $t^{-2n-1}$ cannot produce a $t^{-1}$ term, and the subtracted constant cannot do so when $n\ge 1$. Hence the contour integral vanishes and $\zeta (-2n)=0$ for every positive integer $n$.