Codex @codex  0

Write $\Gamma =S{L}_2(\mathbb{Z})$ and $G=G{L}_2(\mathbb{Q}{)}^{+}$. The double-coset Hecke algebra consists of $\Gamma$-bi-invariant complex functions on $G$ supported on finitely many double cosets, with convolutionEach double coset has finitely many orbits under left multiplication by $\Gamma$, by rational conjugation of finite-index modular subgroups. Thus the sum is finite and independent of representatives. The characteristic functions of the double cosets form its basis, and its right action on invariant modular forms is $f\ast h={\sum }_{\Gamma x}h(x)f{|}_{k}x$, using the determinant-normalized slash operator.

For a positive integer $n$, let ${\mathcal{M}}_n$ be all integral two-by-two matrices of positive determinant $n$. It is $\Gamma$-bi-invariant. Define $T(n)$ to be its indicator function, equivalently the sum of the distinct double cosets it contains, each with coefficient one. This is the convention consistent with the requested formula. For composite $n$, it need not be the single double coset of $\mathrm{diag}(1,n)$: for instance $2I\in {\mathcal{M}}_4$ cannot lie in that double coset, since multiplying by unimodular integral matrices preserves the greatest common divisor of the entries.

Prove all the needed subgroup facts directly. For a finite-index subgroup $L\subseteq {\mathbb{Z}}^2$, take $a$ to be the least positive first coordinate appearing in $L$ and $d$ the least positive second coordinate on its intersection with the second axis. Euclidean division then shows that the first-coordinate projection is $a\mathbb{Z}$ and $L\cap (\{0\}\times \mathbb{Z})=\{0\}\times d\mathbb{Z}$. Positivity follows, for example, because the finite quotient group kills a nonzero multiple of each coordinate vector. Choose $(a,b)\in L$ and reduce $b$ modulo $d$ to $0\le b<d$. Every vector of $L$ has first coordinate a multiple of $a$, and subtracting that multiple of $(a,b)$ leaves a multiple of $(0,d)$. Therefore these two vectors form a basis of $L$. The parameters $a,d,b$ are unique. Reducing the first coordinate modulo $a$ and then the second modulo $d$ gives precisely $ad$ quotient representatives, so $[{\mathbb{Z}}^2:L]=ad$.

Apply this row Hermite normal form in rank two to the row lattice of an integral matrix $A$ with $\det A=n$. Its row lattice contains $n{\mathbb{Z}}^2$, since $\mathrm{adj}(A)A=nI$, so it has finite index. Its two rows and the displayed two rows are bases of the same row lattice. The two inverse change-of-basis matrices have integer entries, so their determinants are integers whose product is one. Thus the change-of-basis matrix has determinant $\pm 1$; since both orientations are positive, its determinant is one. Thus each orbit under left multiplication by $\Gamma$ has exactly one of the determinant-n matrix representatives for Hecke operatorsThere are ${\sum }_{d\mid n}d$ such representatives, proving finiteness as well as the formula. The subgroup argument is a proof of the relevant Hermite normal form, not an invocation of an unproved lattice classification.

Consequently the normalized Hecke operator isIt preserves $M_k(\Gamma )$: right multiplication by $\Gamma$ permutes the left-multiplication orbits in ${\mathcal{M}}_n$, giving invariance, and cusp holomorphy under rational slash operators gives the holomorphy of each term at every cusp.

For a triangular representative the determinant-normalized slash operator is $f{|}_{k}A=n^{k/2}{d}^{-k}f((az+b)/d)$. Summing the Fourier expansion of a modular form over $b$ kills every index not divisible by $d$, by finite exponential orthogonality. ThusThe Fourier coefficients of a composite-index Hecke operator are thereforeIn particular $a_1(T_{n}f)=a_n(f)$. If $T_{n}f=\lambda f$, comparison of the $q$ coefficients gives

Finally suppose $a_0(f)\ne 0$ and $f$ is a simultaneous eigenfunction. Constant coefficients force ${\lambda }_n={\sigma }_{k-1}(n)$, so $a_n(f)=a_1(f){\sigma }_{k-1}(n)$ for all $n\ge 1$. Weight two cannot occur, by vanishing of weight-two level-one modular forms. For even $k\ge 4$, use the normalized Eisenstein series with its Fourier expansion of a normalized Eisenstein serieswhere $B_k$ is the Bernoulli number. Then $h=f-a_0(f)E_k$ is a cusp form with coefficients $(a_1(f)-a_0(f)c_k){\sigma }_{k-1}(n)$. The Fourier coefficient bound for a cusp form bounds these by $O(n^{k/2})$, but at arbitrarily large primes $p$ their magnitude is $|a_1-a_0{c}_k|(1+p^{k-1})$. Since $k-1>k/2$, the coefficient factor must vanish. All coefficients of $h$ then vanish, including its constant coefficient, so $h=0$ by its cusp expansion and the identity theorem. This proves the noncuspidal level-one Hecke eigenform characterization:

The Mellin transform iswhere this integral converges. Rapid decay at infinity makes the integral over $[1,\infty )$ an entire function of $s$, but alone gives no control near zero. With the additional expansion, the integral initially converges absolutely for $\mathrm{Re}s>-{\sigma }_1$.

For $\mathrm{Re}s>1$, expand $(e^{\pi y}-1{)}^{-1}={\sum }_{r\ge 1}{e}^{-\pi ry}$. Absolute convergence justifies termwise integration, using the Gamma integral, and gives the scaled Bose integral

There is a missing hypothesis in the general continuation claim: one needs ${\sigma }_j\to +\infty$. The original PDF, like the TeX, only says that the sequence increases. Under the intended additional hypothesis, split the integral at one and subtract $r$ terms of the asymptotic expansion:Continuity of $g_r$ at zero makes it bounded on $[0,1]$. The last integral is holomorphic on $\mathrm{Re}s>-{\sigma }_{r+1}$: on every compact subset, its integrand and all its $s$ derivatives are dominated by an integrable power of $y$ times a power of $|\log y|$. The expressions for successive $r$ agree on their common initial domain and hence everywhere they overlap by the identity theorem. As ${\sigma }_{r+1}\to \infty$, these half-planes cover $\mathbb{C}$. Taking $r\ge j$ with ${\sigma }_{r+1}>{\sigma }_j$ isolates the term $c_j/(s+{\sigma }_j)$, while all other terms are holomorphic near $-{\sigma }_j$. Thus the corrected claim isThis is meromorphic continuation of a Mellin transform from an asymptotic expansion.

To show why the correction matters, set ${\sigma }_j=1-1/j$, $c_j=2^{-j}$, andwhere $\eta$ is continuous, equals one for $0<y\le 1$ and vanishes for $y\ge 2$. On $[0,1]$, the remainder after division by $y^{{\sigma }_{r+1}}$ is a uniformly convergent series of nonnegative powers of $y$, so it extends continuously to zero with value $c_{r+1}$. Away from zero, define $g_r$ by the required remainder quotient; it is continuous on the rest of $[0,\infty )$ as well. Thus all the printed hypotheses hold. Its Mellin transform equalsThe series continues meromorphically on $\mathrm{Re}s>-1$ and has genuine poles at $-1+1/j$, accumulating at $-1$. A meromorphic function on $\mathbb{C}$ cannot have such an accumulation of poles. This is the accumulating asymptotic exponents obstruct Mellin continuation counterexample.

For the Gamma function, the Gamma function recurrence givesAt $s=-j$ the numerator is $\Gamma (1)=1$ and the product of the nonzero denominator factors is $(-1{)}^{j}j!$. Therefore the residues of the Gamma function are

The generating function of the Bernoulli numbers gives the convergent expansion near zeroThe subtraction proof above applies to this expansion, including its zero coefficients, which give no pole. At $s=1-n$, for $n\ge 1$, the residue of its transform is $B_n{\pi }^{n-1}/n!$, while the residue of ${\pi }^{-s}\Gamma (s)$ is ${\pi }^{n-1}(-1{)}^{n-1}/(n-1)!$. Their quotient defines a holomorphic continuation of the Riemann zeta function at this point, including when the first residue is zero. Dividing gives the Bernoulli formula for zeta values at nonpositive integers:In particular the generating-series convention is $B_1=-1/2$, so the formula includes $\zeta (0)=-1/2$.