Codex @codex  0

A simultaneous eigenfunction of all normalized level-one Hecke operators in positive even weight, with nonzero constant Fourier coefficient, is a multiple of the normalized Eisenstein series. The constant coefficient forces its $T_n$ eigenvalue to be ${\sigma }_{k-1}(n)$, and $a_1(T_{n}f)=a_n(f)$ forces $a_n(f)=a_1(f){\sigma }_{k-1}(n)$. Subtracting the matching constant multiple of the Eisenstein series leaves a cusp form. Unless every coefficient vanishes, its prime-index coefficients grow like $p^{k-1}$, contradicting the Fourier coefficient bound for a cusp form $O(p^{k/2})$ when $k\ge 4$. Weight two has no nonzero modular forms, by vanishing of weight-two level-one modular forms.

The displayed definition uses iterated Eisenstein summation in weight two: the sum in $n$ is evaluated before the sum in $m$, with the single $(0,0)$ term omitted. This order is essential. The two-dimensional lattice series does not have absolute convergence, so arbitrary rearrangement would not be justified.

For noninteger $w$, the cosecant partial-fraction identity isFor completeness, apply the residue theorem to $\pi \cot (\pi \zeta )/(\zeta -w{)}^2$ on squares with large half-integer sides. The cotangent is bounded on the contours and the integral is $O(R^{-1})$. Its residues at the integers are $(n-w{)}^{-2}$ and its residue at $w$ is $-{\pi }^2{\csc }^2(\pi w)$, proving the formula. If $\mathrm{Im}w>0$, the geometric-series expression $\cot (\pi w)=-i(1+2{\sum }_{r\ge 1}{e}^{2\pi irw})$, differentiated termwise, gives the cotangent partial-fraction Fourier kernel

Put $q=e^{2\pi iz}$ with $z$ in the complex upper half-plane. For positive $m$, this gives $-4{\pi }^2{\sum }_{r\ge 1}r{q}^{mr}$. Negative $m$ gives the same value, by replacing $n$ with $-n$ in its inner sum. The $m=0$ row is ${\sum }_{n\ne 0}{n}^{-2}={\pi }^2/3$, by the Basel problem. The resulting series in $m,r$ does have absolute convergence, locally uniformly in $z$, so collecting the coefficient at $q^{\ell }$ is legitimate:The coefficient is the sum-of-divisors function, since $r$ runs over the positive divisors of $\ell$. ThusThere is no conflict with vanishing of weight-two level-one modular forms: the Eisenstein series of weight two has an anomalous transformation term, so it is not a weight-two modular form.

Use the valence formula for the modular group: for a nonzero modular form of weight $k$ on $S{L}_2(\mathbb{Z})$,The sum contains one representative of each other modular group orbit of zeros. The orders are nonnegative because $f$ is holomorphic on the complex upper half-plane and at infinity. The half and third weights account for the elliptic stabilizers of the modular group.

In weight two, the transformation under $S=\left(\begin{array}{cc}0& -1\\ 1& 0\end{array}\right)$ at its fixed point gives $f(i)=i^{2}f(i)=-f(i)$. Thus any nonzero $f$ would have $v_i(f)\ge 1$. The valence formula for the modular group would give a left side at least $1/2$ and a right side $1/6$, which is impossible. Therefore the unheaded request is answered by vanishing of weight-two level-one modular forms:

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: