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:

We use left multiplication by $\Gamma =S{L}_2(\mathbb{Z})$. For an integral matrix $A$ of determinant $n>0$, the Bezout identity gives a determinant-one row operation sending its first column to $(a,0{)}^t$, where $a$ is the positive gcd of that column. The resulting matrix is $\left(\begin{array}{cc}a& b\\ 0& d\end{array}\right)$ with $d=n/a>0$. Adding a multiple of the second row to the first makes $0\le b<d$.

These representatives are unique. Left multiplication by a unimodular matrix preserves the gcd of the first column, so two representatives in one orbit have the same $a$ and hence $d$. A matrix taking $(a,0{)}^t$ to itself has the form $\left(\begin{array}{cc}1& t\\ 0& 1\end{array}\right)$, so it changes $b$ by $td$; the prescribed range makes $b$ unique. Thus the determinant-n matrix representatives for Hecke operators are exactly ${\Pi }_n$. In the printed set, $0\le b<d$ already forces $d>0$, and $ad=n>0$ forces $a>0$.

For positive determinant extend the slash operator for modular forms byThe determinant factor makes this a right action: $(f{|}_{k}A){|}_{k}B=f{|}_k(AB)$. Right multiplication by $\delta \in \Gamma$ permutes the left cosets represented by ${\Pi }_n$, since it preserves the set of determinant-$n$ integral matrices. Hence $T_{n}f$ satisfies the modular transformation law. Its summands are holomorphic on the upper half-plane. For an upper-triangular representative,Substitute the Fourier expansion of a modular form. The sum over $b$ is zero unless the original index is divisible by $d$, in which case it equals $d$. Writing that index as $dr$ givesOnly nonnegative powers occur, so $T_{n}f$ is holomorphic at infinity, and all level-one cusps are equivalent to infinity. We have proved $T_{n}f\in M_k(\Gamma )$ and the Fourier coefficients of a composite-index Hecke operator formulaHere $\gcd (n,0)=n$, giving $a_0(T_{n}f)={\sum }_{a\mid n}{a}^{k-1}{a}_0(f)$. In particular cusp forms stay cusp forms. For nonzero level-one cusp forms the weight is an even integer at least twelve, so the powers $a^{k-1}$ are integers; the zero cusp space causes no exception. Thus every $T_n$ preserves integral Fourier coefficients, and closure under addition, multiplication and integer scalars provesThis is the integral Hecke algebra of level-one cusp forms.

We will also need commutativity. The coefficient formula directly givesFor the first identity, the divisors of $u$ and $v$ in the iterated coefficient sum combine uniquely into a divisor of $uv$. For the second, write a coefficient index as $p^t{r}_0$ with $p\nmid r_0$; the formula isApplying $T_p$ splits this into the two sums with shifts $t+1$ and $t-1$; their overlapping terms give exactly the stated recurrence. Since $T_1$ is the identity, induction expresses each $T_{p^r}$ as a polynomial in $T_p$, and coprime multiplicativity expresses every $T_n$ as a product of these polynomials. Distinct prime operators commute by the coprime identity. Hence all Hecke operators commute, and so do all elements of $\mathbb{T}$.

Let $\Lambda =S_k(\Gamma ,\mathbb{Z})$ and $m=\dim S_k(\Gamma )$. In the permitted basis $F_j={\Delta }^j{E}_4^{a_j}{E}_6^{b_j}$, the modular discriminant begins $q+O(q^2)$ and the normalized Eisenstein series begin $1+O(q)$, soThe matrix of $a_1,\dots ,a_m$ on this integral basis is unitriangular and has determinant one. Integer elimination therefore gives an integral basis $f_1,\dots ,f_m$ with $a_i(f_j)={\delta }_{ij}$ for $1\le i,j\le m$. Thus $a_1,\dots ,a_m$ are a $\mathbb{Z}$-basis of the dual module ${\mathrm{Hom}}_{\mathbb{Z}}(\Lambda ,\mathbb{Z})$.

The map $\alpha$ is well-defined by lattice preservation and is $\mathbb{Z}$-linear. Because $\alpha (T_j)=a_j$, it is surjective. For injectivity suppose $\alpha (T)=0$. For every $f\in \Lambda$ and every $n\ge 1$, commutativity givesThe constant term is also zero because $Tf$ is a cusp form. The identity theorem applied to its Fourier expansion of a modular form gives $Tf=0$. The integral basis spans the complex cusp space, so $T=0$ as an endomorphism. This proves the perfect integral Hecke pairing andIf $m=0$, both modules are zero and the asserted basis is empty.

For a finite-index subgroup $L\le {\mathbb{Z}}^2$, its projection to the first coordinate is $a\mathbb{Z}$ and its intersection with the second axis is $\{0\}\times d\mathbb{Z}$, with $a,d>0$. Choose $(a,b)\in L$, reducing $b$ modulo $d$. Together with $(0,d)$ it is a basis of this row lattice: subtracting a multiple of $(a,b)$ from any vector leaves a vector on the second axis. These parameters are unique, and reduction of the two coordinates shows $[{\mathbb{Z}}^2:L]=ad$. Consequently an integral matrix of positive determinant $n$ has a unique representative of this form under left multiplication by $S{L}_2(\mathbb{Z})$, with $ad=n$. This proves the lattice facts underlying determinant-n matrix representatives for Hecke operators.