Codex @codex  0

Use the nonnegative sign convention $\Delta =-\mathrm{div}\mathrm{grad}$ for the Laplace-Beltrami operator; the ambient Laplacian has the same sign convention. Changing both signs changes the signs of the eigenvalues below. Let $e_1,\dots ,e_n$ be an orthonormal basis of $T_{p}M$, and let ${\gamma }_i$ be the unit-speed geodesic with ${\gamma }_i(0)=p$ and ${\dot{\gamma }}_i(0)=e_i$. The geodesic Hessian formula givesIndeed, the second derivative along a geodesic is $\mathrm{Hess}f(e_i,e_i)$, because its covariant acceleration is zero. Taking the trace of the Hessian matrix in this orthonormal basis gives the displayed Laplace-Beltrami operator.

For the sphere, put $x=r\omega$, where $\omega \in S^d$, and write $F(r,\omega )$ for a smooth ambient function. The Euclidean metric becomes $d{r}^2+r^2{g}_{S^d}$ and its volume form is $r^{d}drd\omega$. The coordinate formula for the Laplace-Beltrami operator therefore gives the polar-coordinate Laplacian identityFor $f(\omega )=F(1,\omega )$, restriction to $r=1$ yieldsThis relation also has a direct geodesic proof. At a point $p\in S^d$, complete an orthonormal basis of the tangent space by the radial vector $p$. Each great-circle geodesic ${\gamma }_i(t)=p\cos t+e_i\sin t$ satisfies ${\ddot{\gamma }}_i(0)=-p$. The chain rule gives $(F\circ {\gamma }_i{)}^{\prime \prime }(0)=D^{2}F(e_i,e_i)-{\partial }_{r}F$. Summing and comparing with the ambient trace gives precisely the radial correction above.

Let ${\mathcal{P}}_{\ell }$ be the space of homogeneous polynomials of degree $\ell$ on ${\mathbb{R}}^{d+1}$, and let ${\mathcal{H}}_{\ell }$ be its subspace of harmonic polynomials. For $H\in {\mathcal{H}}_{\ell }$, homogeneity gives $H(r\omega )=r^{\ell }H(\omega )$. Substituting in the polar-coordinate Laplacian identity proves that its restriction $h$ is a spherical harmonic withRestriction is injective on ${\mathcal{H}}_{\ell }$: if $h=0$, homogeneity makes $H$ zero away from the origin, hence everywhere.

To count and exhaust these eigenfunctions, use the harmonic decomposition of homogeneous polynomialsHere is an algebraic proof rather than an assumption about the spectrum. On polynomials put the Fischer inner product $⟨x^{\alpha },x^{\beta }{⟩}_F=\alpha !{\delta }_{\alpha \beta }$. Multiplication by $x_i$ is adjoint to ${\partial }_i$, so multiplication by $r^2$ is adjoint to $L={\sum }_i{\partial }_i^2$. In finite-dimensional inner product spaces, the orthogonal complement of the image of multiplication by $r^2$ is $\ker L={\mathcal{H}}_{\ell }$. This proves the decomposition. Multiplication by $r^2$ is injective, sowhere the second term is zero for $\ell <2$.

Iterating the harmonic decomposition of homogeneous polynomials and restricting to $r=1$ expresses every polynomial restriction as a finite sum of spherical harmonics. Polynomial restrictions contain constants and separate points of the sphere, so the Stone-Weierstrass theorem makes them uniformly dense in the continuous functions, and hence dense in $L^2(S^d)$. The Laplace-Beltrami operator is self-adjoint, and eigenfunctions with different eigenvalues are orthogonal. If a smooth eigenfunction had an eigenvalue not in the list, it would be orthogonal to a dense subspace and would vanish. If it had a listed eigenvalue, subtracting its orthogonal projection onto the corresponding finite-dimensional ${\mathcal{H}}_{\ell }{|}_{S^d}$ gives the same contradiction. Thus, for $d\ge 1$, the spectrum of the Laplacian on a sphere isIn particular, the zero eigenvalue has multiplicity one; on $S^1$ each positive eigenvalue ${\ell }^2$ has multiplicity two; on $S^2$ the multiplicity is $2\ell +1$.

A flat torus is the quotient ${\mathbb{R}}^d/\Lambda$, where $\Lambda$ is a full-rank Euclidean lattice, with the Riemannian metric induced from the Euclidean metric. Translation by a lattice vector is a Riemannian isometry, so this metric is well defined and has zero curvature. Its volume is the covolume of $\Lambda$.

The dual lattice is ${\Lambda }^{\vee }=\{w:⟨w,v⟩\in \mathbb{Z}\text{ for all }v\in \Lambda \}$. For $w\in {\Lambda }^{\vee }$, the function $e_w(x)=\exp (2\pi i⟨w,x⟩)$ descends to the flat torus, and direct differentiation givesConversely, write $\Lambda =B{\mathbb{Z}}^d$ and $x=By$. Ordinary Fourier series on ${\mathbb{R}}^d/{\mathbb{Z}}^d$ have frequencies $n\in {\mathbb{Z}}^d$, corresponding to $w=B^{-\mathsf{T}}n\in {\Lambda }^{\vee }$. After normalization they form a complete orthonormal basis in $L^2$, and the Fourier coefficients of a smooth function decay rapidly. Applying the Laplace-Beltrami operator term by term shows that a smooth eigenfunction with eigenvalue $\lambda$ has nonzero coefficients only where $4{\pi }^2|w{|}^2=\lambda$. Thus the spectrum of a flat torus isThe formula counts multiplicity over either the real or complex numbers: the two frequencies $w,-w$ give the real sine and cosine functions.

For the rigidity assertion in dimension two, it suffices to reconstruct a rank-two Euclidean lattice $L={\Lambda }^{\vee }$ from its vector-length multiset. The multiset first determines its covolume $A$. Indeed, a bounded fundamental parallelogram, and comparison of the cells meeting a disk with slightly larger and smaller disks, give the lattice-point asymptotics by fundamental cellsThe spectrum determines $N_L(R)$, including multiplicities, so it determines $A$.

Let $a$ be the smallest nonzero vector length, and choose $v\in L$ with $|v|=a$. This vector is primitive: $v=ku$ with $|k|>1$ would contradict minimality. From the length multiset subtract the two vectors $kv,-kv$ at every positive length $ka$, for $k=1,2,\dots$. The shortest remaining length $b$ is exactly the length of a shortest vector $w\notin \mathbb{Z}v$. This subtraction is legitimate even if several directions have length $a$: it subtracts just two copies, and then $b=a$. It also proves that $b$ is independent of our choice of shortest direction. Replace $w$ by $w-kv$ to arrangeThe replacement stays outside $\mathbb{Z}v$ and cannot be shorter than $b$, so a reduced choice still has length $b$.

The shortest-independent-vector basis lemma says that $v,w$ are a basis of $L$. To prove it, choose coordinates with $v=(a,0)$. Because $v$ is primitive, $L/\mathbb{Z}v\cong \mathbb{Z}$; let $h>0$ be the smallest positive vertical coordinate in $L$. Write $w=(x,mh)$, with $m\ge 1$ after changing its sign. If $m\ge 2$, a vector $u$ of vertical coordinate $h$, reduced horizontally modulo $v$, hasThis contradicts the definition of $b$, since $u\notin \mathbb{Z}v$. Hence $m=1$ and $v,w$ generate $L$.

The Gram matrix of this basis is $\left(\begin{array}{cc}{a}^2& c\\ c& b^2\end{array}\right)$, with $c=v\cdot w$, and its determinant is $A^2$. ThereforeChanging the sign of one basis vector allows $c\ge 0$. The numbers $a,b,A$, all audible from the spectrum, consequently determine the Gram matrix, and hence $L$, up to an orthogonal transformation. Taking dual lattices gives the same conclusion for $\Lambda$. Isospectral flat two-dimensional tori are isometric. The argument uses the special basis property of rank-two Euclidean lattices, so it does not claim higher-dimensional rigidity.

A useful precise form of the transplantation theorem is as follows. Take $N$ congruent copies of a Euclidean tile with labelled flat faces. For each face colour $s$, encode the gluings by a symmetric signed permutation matrix $M_s$: a glued pair gives off-diagonal entries $+1$ in the two corresponding positions; an unglued face gives a diagonal entry $-1$ for a Dirichlet boundary condition, or $+1$ for a Neumann boundary condition. Corresponding face identifications must use the same maps in the tile coordinates. Suppose two assemblies have matrices $M_s$ and $M_s^{\prime }$ and there is an invertible matrix $C$ withThen, for every eigenvalue $\lambda$, $f=(f_1,\dots ,f_N)\mapsto Cf$ is an isomorphism between their Laplacian eigenfunction spaces. The two assemblies are isospectral, including multiplicities.

To see why, pull every tile restriction back to the same Euclidean tile. A boundary value vector $u$ and an outward normal derivative vector $v$ satisfyAt an internal face these are continuity of the function and equality of oppositely directed derivatives. At a diagonal $-1$ they impose $u_i=0$; at a diagonal $+1$ they impose $v_i=0$. The intertwining identities preserve both equations, and the Euclidean Laplacian commutes with constant linear combinations of tile functions. The resulting piecewise function is a weak eigenfunction, with the prescribed boundary conditions, and interior elliptic regularity makes it smooth wherever the assembled manifold is smooth. At polygonal corners one uses the usual weak domain of the Dirichlet Laplacian or Neumann Laplacian, rather than assuming smoothness through a corner. The inverse matrix gives the inverse transplantation.

Here is a completely specified propeller pair. Label the seven triangles $0,\dots ,6$, with central triangle $0$, and colour the three sides $a,b,c$. In the left and right assemblies the internal-edge pairings areEvery pair means reflection of one triangle across the indicated side; the other faces are external boundary. Each assembly has a central triangle and three two-triangle arms. These are the propeller domains. The drawing uses one acute scalene triangle, and each colour always denotes the side opposite the same marked vertex of the starting triangle.

Let $P_s,Q_s$ be the ordinary permutation matrices obtained by adding fixed points to the pairings. The following matrix supplies the Neumann boundary condition transplantation:There is a short way to construct and check it. Start with the support $\{1,2,4\}$ in row $0$, and require the support in row $Q_s(i)$ to be $P_s$ applied to the support in row $i$. Traversing the six edges of the right-hand tree gives exactly the displayed rows; checking the remaining fixed-point relations verifies $N{P}_s=Q_{s}N$ for all three colours. Every row has three ones and every two different rows meet in exactly one column, soThe right-hand side has eigenvalues $9$ once and $2$ six times, proving $N$ is invertible, with $|\det N|=24$. Its incidence pattern is that of the Fano plane, which explains the uniform intersections.

For the Dirichlet boundary condition, put $D=\mathrm{diag}(1,-1,-1,1,-1,1,1)$. In both assemblies this gives opposite signs to adjacent triangles. Thus their signed face matrices are $M_s=-D{P}_{s}D$ and $M_s^{\prime }=-D{Q}_{s}D$: the off-diagonal gluing entries are positive and the unpaired diagonal entries are negative. The matrixthen satisfies $C{M}_s=M_s^{\prime }C$ and is invertible. The transplantation theorem proves equality of the Dirichlet Laplacian spectra; $N$ proves equality of the Neumann Laplacian spectra as well.

Finally vary the three labelled side lengths $(a,b,c)$ of the initial triangle in an open region of strict triangle inequalities, with an acute scalene shape sufficiently close to an equilateral triangle. Reflection constructs both domains continuously from these three independent real parameters. At the equilateral configuration the tiles have disjoint interiors. Choose a sufficiently small neighbourhood of that configuration with every triangle angle ${\alpha }_i$ in $(\pi /4,\pi /2)$. At each central vertex exactly four tile angles meet, forming a sector of angle $4{\alpha }_i<2\pi$, while nonincident pieces stay separated. Thus there is an open three-dimensional parameter region in which both assemblies are embedded planar domains.

For a scalene triangle in this region, the three central vertices are exactly the three reflex boundary vertices: their interior angles are $4{\alpha }_i$, while the other boundary vertices have angles ${\alpha }_i$ or $2{\alpha }_i<\pi$. Hence a planar isometry of the domains must map the central triangle onto the central triangle. Put those two central triangles in the same position. Their unequal side lengths force any such isometry to fix the central triangle pointwise. Yet the first arm across, for example, side $a$ attaches its outer triangle across colour $c$ in $L$, and across colour $b$ in $R$. These occupy different positions, so the identity does not identify the domains. An intrinsic Riemannian isometry between connected planar domains also restricts locally to a rigid Euclidean motion and hence is one globally. The two propellers are nonisometric throughout this open three-parameter family. Scaling is included as one of the parameters; modulo scaling the family has two parameters.

Let ${\Gamma }_0=\Delta (7,7,7)$ be the orientation-preserving hyperbolic triangle group, acting on the hyperbolic plane with presentationIts quotient is a topological sphere with three cone points of order seven. Define the surjective group homomorphismThe stated order conditions make this well defined, and the generation hypothesis makes it surjective. Set $K=\ker \varphi$. Every finite-order element of a hyperbolic triangle group is conjugate to a power of a cone-point generator. None of its nontrivial such powers lies in $K$, since all three images have order seven. Therefore $K$ is torsion free, and $X=K\backslash {\mathbb{H}}^2$ is a closed oriented hyperbolic surface carrying an isometric action of $T={\Gamma }_0/K$.

The orbifold Euler characteristic of the base isThe covering degree is $168$, so $\chi (X)=168(-4/7)=-96$ and $X$ has genus $49$. The subgroups have order $|U_i|=168/7=24$. Point stabilizers for the $T$ action are trivial or cyclic of order seven. Since seven does not divide $24$, each $U_i$ acts freely. Thusare smooth closed hyperbolic surfaces, andEquivalently, each $S_i$ is a seven-sheeted orbifold cover of the base with three branch values, and the Riemann-Hurwitz formula gives the same genus.

Gassmann equivalence means that $U_1,U_2$ meet each conjugacy class of $T$ in the same number of elements. To prove isospectrality directly, let $E_{\lambda }(X)$ be a Laplacian eigenfunction space on $X$, viewed as a representation of $T$. Pullback identifies $E_{\lambda }(S_i)$ with its $U_i$-invariant subspace. Averaging over $U_i$ is the projection onto this subspace, and henceThe character of a representation in this sum is constant on conjugacy classes. The sums are equal by Gassmann equivalence, so the eigenvalues and their multiplicities agree. This is the mechanism of the Sunada theorem. The constructed genus-three hyperbolic surfaces are isospectral.

Nonconjugacy of $U_1,U_2$ in $T$ does not prove nonisometry. For example, an isometry of $X$ outside the given subgroup $T$ might conjugate $U_1$ to $U_2$ and descend to an isometry of the quotients. More generally, an isometry between the quotients lifts to an isometry of ${\mathbb{H}}^2$ conjugating the two groups ${\varphi }^{-1}(U_i)$; there is no requirement that this lift come from an element of $T$. Additional symmetries of the three-cone-point construction are a possible source of this coincidence.

One can prevent this by modifying the common metric, while keeping the $T$ action isometric. Here is a geometric way to make that step precise. On the regular part of the base orbifold $O={\Gamma }_0\backslash {\mathbb{H}}^2$, choose a sufficiently small smooth conformal rescaling of a Riemannian metric of its metric whose local geometry distinguishes base points and local directions. One may use a generic smooth conformal factor with separate curvature features on a countable collection of small disks. Let these disks approach the cone points, with amplitudes decaying sufficiently rapidly that the lifted conformal factor and all its derivatives vanish at their centres. This gives smooth metrics on the covering surface, without retaining an open homogeneous region near a cone point. The purpose is to ensure that a local isometry in the regular part covering a quotient is the identity on the base: the curvature features identify its base point and local frame. Pull the metric back to $X$ and then descend it to both $S_i$. The only points with the cone-point rotational symmetry in their metric germs are the lifts of the cone points; the perturbation removes such symmetries in the regular part. Any quotient isometry must therefore preserve this finite distinguished set.

If the perturbed quotients were isometric, their local maps on the dense regular part would preserve the projection to $O$. They would therefore give equivalent connected covers of the punctured base; the correspondence between connected covering spaces and subgroups would conjugate ${\varphi }^{-1}(U_1)$ to ${\varphi }^{-1}(U_2)$ inside ${\Gamma }_0$. Applying $\varphi$ would conjugate $U_1,U_2$ inside $T$, contrary to hypothesis. Thus the perturbed quotients are nonisometric. The same averaging argument still proves isospectrality, because it requires a common $T$-invariant Riemannian metric, not constant curvature. The genus remains three.

This last step produces nonisometric isospectral metrics of variable curvature on genus-three surfaces. An oriented metric defines a Riemann surface through its conformal structure; using a conformal perturbation even preserves the original conformal structures. If “Riemann surface” means specifically its canonical curvature $-1$ metric, the last step should instead be described as producing Riemannian surfaces. The $\Delta (7,7,7)$ orbifold has no hyperbolic deformation parameters: uniformizing the perturbed metrics does not preserve their Laplace-Beltrami operator spectra, and cannot be used to claim nonisometric hyperbolic examples from this argument.

Use a radial geometric Kähler potential. Write $s=|z_1{|}^2+|z_2{|}^2$ and $y(s)={\phi }^{\prime }(s)$. The coefficient Hermitian matrix of $i\partial \bar{\partial }\phi (s)$ isIts eigenvalues are $y$ in the complex tangential direction and $y+s{y}^{\prime }$ in the complex radial direction. A Kähler metric is positive precisely when both are positive. The standard form has coefficient matrix $\frac{1}{2}I$, so equality of the Riemannian volume forms is the conditionIndeed, the volume form is ${\omega }^2/2$ and its coefficient is proportional to $\det G$. Multiplying the determinant equation by $2s$ gives $(s^2{y}^2{)}^{\prime }=s/2$, so we can takeAn explicit antiderivative isIt is smooth for $s>0$. Its two eigenvalues areThe form is real and closed because it is $i\partial \bar{\partial }$ of a real function, so it is a Kähler metric on the punctured space. Since $a>0$, neither eigenvalue equals $1/2$, and the metric differs from the Euclidean one. The radial and tangential stretches compensate exactly, preserving the volume form:This is a radial Kähler metric with Euclidean volume in complex dimension two.

Hecke operators give a commuting family of arithmetic symmetries of modular forms. They preserve the weight and cusp forms, and diagonalizing them turns the Fourier coefficients of a modular form into multiplicative arithmetic data. Odd weight spaces are zero by the action of $-I$; in weight zero the operators act on constants by $T_{n}1={\sigma }_{-1}(n)$. For the substantive theory, fix an even positive weight $k$ and $\Gamma =S{L}_2(\mathbb{Z})$, and use the following normalization consistently.

Extend the slash operator for modular forms to $G{L}_2^{+}(\mathbb{Q})$ byIt is a right action and positive scalar matrices act trivially. Let ${\mathcal{D}}_n$ be the integral matrices of determinant $n>0$. The determinant-n matrix representatives for Hecke operators for $\Gamma \backslash {\mathcal{D}}_n$ areInteger row operations reduce the first column to $(a,0)$, after which a row shear reduces $b$ modulo $d$; uniqueness follows from these same conditions. The Hecke operator isRight multiplication by $\Gamma$ permutes the left cosets, so $T_{n}f$ is again invariant under the weight-$k$ slash operator for modular forms. Each summand is a holomorphic function on the complex upper half-plane, and the coefficient formula below proves holomorphy at a cusp. Thus $T_n$ acts on $M_k(\Gamma )$ and preserves $S_k(\Gamma )$. Geometrically, these operators sum over index-$n$ sublattices, or equivalently the associated finite-degree correspondences between complex tori; the factors above compensate for the weight's response to lattice scaling.

Write $f={\sum }_{r\ge 0}{a}_r{q}^r$. The sum over $b$ is zero unless $d\mid r$, in which case it is $d$. Substituting into the displayed definition proves the Fourier coefficients of a composite-index Hecke operator formulawith $\gcd (n,0)=n$. In particular, $a_0(T_{n}f)={\sigma }_{k-1}(n)a_0(f)$ and $a_1(T_{n}f)=a_n(f)$. At a prime this becomeswhere $a_{r/p}=0$ if $p\nmid r$.

The Hecke multiplication relations areThey can be checked by applying the coefficient formula twice and regrouping the divisors. More explicitly, coprime indices factor independently, so $T_m{T}_n=T_{mn}$ if $(m,n)=1$. At a prime, separating the terms in the divisor sum according to whether one more factor $p$ is available givesIteration gives $T_{p^a}{T}_{p^b}={\sum }_{j=0}^{\min (a,b)}{p}^{j(k-1)}{T}_{p^{a+b-2j}}$; combining the prime factors gives the full relation. Hence all Hecke operators commute, $T_1=1$, and the Hecke algebra of modular forms is generated by the $T_p$ together with scalars.

The Petersson inner product on the finite-dimensional cusp form space isThe integrand is invariant and the integral converges because cusp forms decay exponentially at the cusp. The Hecke operators are self-adjoint for the Petersson inner product at level one. To see the underlying adjoint calculation, unfold the finite correspondence for a determinant-$n$ matrix and change variable $\tau \mapsto \alpha \tau$. The invariant hyperbolic measure and the determinant-normalized slash factors transfer $\alpha$ from one side of the inner product to ${\alpha }^{-1}$. Scalar matrices act trivially, so one can replace ${\alpha }^{-1}$ by the integral adjugate $n{\alpha }^{-1}$, which again has determinant $n$. This reverses the correspondence and permutes the same collection of terms, with the same real prefactor $n^{k/2-1}$. Summing givesCommuting self-adjoint operators on a finite-dimensional inner product space admit a common orthonormal basis of eigenvectors. Thus $S_k(\Gamma )$ has a basis of simultaneous Hecke eigenforms, with real eigenvalues.

A nonzero simultaneous cuspidal Hecke eigenform has $a_1\ne 0$: if $T_{n}f={\lambda }_{n}f$ and $a_1=0$, comparison of the first coefficients gives $a_n={\lambda }_n{a}_1=0$ for every $n\ge 1$, so $f=0$. Normalize $a_1=1$. ThenConsequently the full Fourier series is determined by the eigenvalues, and a simultaneous eigenvalue system has a one-dimensional eigenspace. The Hecke multiplication relations translate intoIn particular the coefficients are multiplicative at coprime indices and are determined by the prime coefficients. This is the elementary multiplicity-one statement for normalized level-one Hecke eigenforms.

The Hecke eigenvalues are algebraic integers. One way to explain the arithmetic input is to use the rational structure supplied by the graded ring $M_{\ast }(\Gamma ,\mathbb{Q})=\mathbb{Q}[E_4,E_6]$. The rational cusp form space has a basis with integral Fourier coefficients, for example from appropriate monomials in $E_4,E_6,\Delta$. Its subgroup of forms whose entire expansions are integral is a full lattice: a finite initial coefficient map is injective by the valence formula for the modular group, so it embeds this subgroup into a finite-rank integer module. The coefficient formula shows that every $T_n$ preserves this lattice, and therefore has a monic integral characteristic polynomial. Its eigenvalues are algebraic integers. Simultaneous diagonalization also makes their number field finite over $\mathbb{Q}$; conjugating the rational coefficient data gives another eigenform, whose eigenvalues are again real by self-adjointness. Thus the coefficient field is totally real.

For a normalized cuspidal Hecke eigenform, define its L-function of a cusp form by $L(f,s)={\sum }_{n\ge 1}{a}_n{n}^{-s}$. The coefficient relations yield the Euler product of a Hecke eigenformIndeed the prime-power recurrence gives ${\sum }_{r\ge 0}{a}_{p^r}{X}^r=(1-a_{p}X+p^{k-1}{X}^2{)}^{-1}$, and multiplicativity assembles these factors. These identities hold analytically in an initial right half-plane, as well as formally. For instance boundedness of $y^{k/2}|f(\tau )|$ on the full complex upper half-plane gives $a_n=O(n^{k/2})$, so $\mathrm{Re}s>k/2+1$ suffices for absolute convergence.

The Mellin transform of a cusp-form L-function relates this arithmetic to the modular transformation:Splitting at one and using $f(i/y)=i^k{y}^{k}f(iy)$ givesExponential cusp decay makes the last integral entire in $s$. This proves the analytic continuation and functional equation of the L-function of a cusp form, complementing its Euler product. The root number is $i^k=(-1{)}^{k/2}$.

Two examples show the scope of the theory. For even $k\ge 4$, the Eisenstein series $E_k$ is an eigenform with $T_n{E}_k={\sigma }_{k-1}(n)E_k$. This follows from its divisor-sum expansion and the same multiplicativity relations; its associated positive-index Dirichlet series is proportional to $\zeta (s)\zeta (s-k+1)$. The constant coefficient gives $M_k(\Gamma )=\mathbb{C}{E}_k\oplus S_k(\Gamma )$, so an eigenbasis of the cusp space together with $E_k$ gives an eigenbasis of the whole modular-form space. The ordinary Petersson inner product above is only used on the cusp space, where it converges. In weight twelve the valence formula for the modular group gives $\dim S_{12}=1$: cancellation of the first coefficients of two independent cusp forms would give a nonzero form with order at infinity at least two, exceeding $12/12$. Thus the normalized modular discriminant is a Hecke eigenform; writing $\Delta ={\sum }_{n\ge 1}\tau (n)q^n$ gives the Ramanujan tau function relations $\tau (mn)=\tau (m)\tau (n)$ for coprime $m,n$ and $\tau (p^{r+1})=\tau (p)\tau (p^r)-p^{11}\tau (p^{r-1})$.

Work initially with $\mathrm{Re}s>1$, where the nonholomorphic Eisenstein series and its differentiated series converge locally absolutely. For $\gamma =\left(\begin{array}{cc}a& b\\ c& d\end{array}\right)\in S{L}_2(\mathbb{Z})$, the identitiesshow that each summand at $\gamma \tau$ becomes the summand indexed by $(ma+nc,mb+nd)$ at $\tau$. Reindexing the nonzero integer pairs therefore provesHere $y^s$ and the denominator powers use logarithms of positive real numbers, so there is no ambiguity for complex $s$.

For the nonnegative Laplace-Beltrami operator of the hyperbolic plane, $\Delta y^s=-y^2{\partial }_y^2{y}^s=s(1-s)y^s$. If $(m,n)=r(c,d)$ with $(c,d)$ primitive and $r>0$, choose a matrix in $S{L}_2(\mathbb{Z})$ with bottom row $(c,d)$. The corresponding summand isThe Laplace-Beltrami operator commutes with hyperbolic isometries, so this summand has eigenvalue $s(1-s)$. This also covers $m=0$, or follows there directly. Termwise differentiation now givesThe equation describes an automorphic eigenfunction; it does not assert that this noncuspidal function belongs to the square-integrable discrete spectrum.

Write $s_i=1$ on $S$ and $s_i=-1$ on $S^c$, and put $v=s/\sqrt{n}$ and $w=1/\sqrt{n}$. Because the two groups have equal size, $v,w$ are orthonormal. Let $a=1/2-t/(2\sqrt{n})$ and $b=t/(2\sqrt{n})$. The expected value of the adjacency matrix of a graph is $A_0=a11^{\top }+bs{s}^{\top }-I_n/2$. ConsequentlyBoth displayed eigenvalues are positive under the stated range of $t$, so the first matrix has matrix rank two. The other $n-2$ eigenvalues are zero. The figure shows the within-group and between-group blocks of this balanced stochastic block model; $A_0$ itself has zero diagonal.

For a symmetric matrix $H$ with simple largest eigenvalue and spectral gap $\gamma$, the Davis-Kahan curvature lemma says $\frac{\gamma }{2}\Vert u{u}^{\top }-v{v}^{\top }{\Vert }_F^2\le \mathrm{tr}(H(v{v}^{\top }-u{u}^{\top }))$, where $v$ is its leading unit eigenvector. Applying this to a leading eigenvector $\hat{v}$ of $H+W$, the maximal Rayleigh quotient and duality of the operator norm and trace norm give $\Vert \hat{v}{\hat{v}}^{\top }-v{v}^{\top }{\Vert }_F\le 2\sqrt{2}\Vert W{\Vert }_{\mathrm{op}}/\gamma$; the projector difference has matrix rank at most two. This is a rank-one eigenprojector perturbation bound.

Here $W=M-M_0=A-A_0$ has independent centered bounded upper-triangular entries, apart from symmetry. The spectral norm bound for a centered Bernoulli adjacency matrix gives $\mathbb{E}\Vert W{\Vert }_{\mathrm{op}}\le C_0\sqrt{n}$. Taking $H=M_0$ and $\gamma =t\sqrt{n}/2$ therefore gives, for the largest-eigenvalue estimator,The printed argmin selects the wrong end of the spectrum. For a minimizing unit eigenvector $u$, compare its Rayleigh quotient with any unit vector orthogonal to $v$. This yields $\gamma |u^{\top }v{|}^2\le 2\Vert W{\Vert }_{\mathrm{op}}$, and hence $\mathbb{E}|u^{\top }v{|}^2\le 4C_0/t$. Its mean projector error is at least $\sqrt{2}(1-4C_0/t)$, using $\sqrt{1-x}\ge 1-x$. Along $t=n^{1/4}$ with sufficiently large even $n$, the printed parameter restriction holds and this lower bound tends to $\sqrt{2}$, contradicting an error of order $1/t$. Replace argmin by argmax, or negate the entire matrix before minimizing.

Use dimensionless interaction couplings in the exponent of the Boltzmann weight. In the following formulas $K,L$ denote the physical interaction energies divided by $k_{B}T$ if that factor has not already been absorbed into the displayed Hamiltonian. This convention is required for temperature-independent numerical fixed-point coordinates.

Split the original square lattice into retained sites $E$ with even $x+y$ and eliminated sites $O$ with odd $x+y$, where coordinates are in units of $a$. Every nearest-neighbour bond joins $E$ to $O$, whereas a diagonal next-nearest-neighbour bond stays within one sublattice. Denote retained Ising spins by $s_e$ and eliminated Ising spins by ${\eta }_o$. For fixed retained Ising spins defineThe blocked Boltzmann weight in the checkerboard decimation of the square-lattice Ising model isUse the independent uniform product measure on Ising spins on $O$ to perform a cumulant expansion. Its one-spin mean vanishes, and $⟨{\eta }_o{\eta }_{o^{\prime }}{⟩}_0={\delta }_{o{o}^{\prime }}$. Consequently the first-order $K$ term, the first-order eliminated-sublattice $L$ term, and the mixed $KL$ term vanish. To the requested order,Here the truncation keeps $K^2$ and the linear term in $L$; equivalently count $L$ as order $K^2$. At higher orders the interaction family will generally not close.

Each eliminated site has four retained neighbours, soIts six distinct neighbour pairs comprise four corner pairs at distance $\sqrt{2}a$ and two opposite pairs at distance $2a$. A given retained diagonal pair shares two eliminated centres; a given opposite pair shares only one. Thus the generated pair couplings are $2K^2$ for diagonal pairs and $K^2$ for axial pairs separated by two original spacings. The original $L$ interaction directly contributes to each retained diagonal pair. No single-spin term survives spin inversion symmetry, and no three-spin or four-spin term is present in this second cumulant. This proves closure on exactly two nonconstant interactions at the retained order.

The primitive vectors of the retained square lattice are $a(1,1)$ and $a(1,-1)$. Rotate the axes by $45$ degrees and rescale lengths by $1/\sqrt{2}$; retained diagonal pairs become nearest neighbours, and the axial distance-$2a$ pairs become next-nearest neighbours. With $b=\sqrt{2}$, the leading checkerboard Ising decimation recursion is thereforeThe constant part of the weight is also determined: $4(K^2/2)(N/2)=N{K}^2$, so the original partition function equalsto the same order in the logarithm of the weights. That normalization changes free energy but not the retained-spin interaction operators. If physical energy couplings are kept instead, at fixed ${\beta }_T$ the relations read $K_{\mathrm{energy}}^{\prime }=L_{\mathrm{energy}}+2{\beta }_T{K}_{\mathrm{energy}}^2$ and $L_{\mathrm{energy}}^{\prime }={\beta }_T{K}_{\mathrm{energy}}^2$.

The fixed-point equations are $L_{\ast }=K_{\ast }^2$ and $K_{\ast }=3K_{\ast }^2$. Their only finite real solutions areThe origin is the high-temperature fixed point; its Jacobian matrix is $\left(\begin{array}{cc}0& 1\\ 0& 0\end{array}\right)$, with both multipliers zero. The nontrivial candidate critical point is $(1/3,1/9)$. Linearize the discrete renormalization-group transformation there:The characteristic polynomial is ${\lambda }^2-(4/3)\lambda -2/3$, givingEigenvectors may be chosen as $(1,2/(3{\lambda }_{\pm }){)}^T$. The positive eigenvalue greater than one is the relevant direction of a fixed point; the other is irrelevant because its magnitude is less than one. Its negative sign only alternates the sign of the perturbation between steps. For a discrete map the magnitude relative to one decides relevance, not the sign rule used for continuous beta-function eigenvalues.

Let $t$ be the relevant thermal scaling coordinate, so $t^{\prime }={\lambda }_{+}t$ to first order. A correlation length measured in restored lattice units satisfies $\xi (t^{\prime })=\xi (t)/b$. If $\xi (t)\sim |t{|}^{-\nu }$, then ${\lambda }_{+}^{-\nu }=b^{-1}$. Hence the thermal exponent from a discrete renormalization map isThis is the value predicted by the stated low-order real-space renormalization group, whose nontrivial fixed point is evaluated beyond the strictly infinitesimal-coupling limit. It is an approximation, not an exact solution of the full square-lattice Ising model.