---
title: Gross–Zagier Formula on Singular Moduli
url: https://www.emergentmind.com/topics/gross-zagier-formula-on-singular-moduli
type: topic
---

# Gross–Zagier Formula on Singular Moduli

Searching arXiv for the primary paper and closely related work on singular moduli, refinements, and Shimura-curve generalizations.
Searching arXiv for "1008.1853 Gross Zagier singular moduli Yang Howard Yang Hilbert modular surfaces".
The Gross–Zagier formula on singular moduli is the explicit factorization formula governing differences of CM values of the modular \(j\)-invariant. If \(\tau\in \mathfrak{H}\) is a CM point, then \(j(\tau)\) is a singular modulus, an algebraic integer generating a ring or Hilbert class field. Gross and Zagier showed that products, or equivalently suitable norms, of differences \(j(\tau_1)-j(\tau_2)\) admit a prime factorization controlled by local arithmetic data: representation numbers, optimal embeddings into quaternionic orders, and local intersection multiplicities on integral models of modular curves. In later work, this viewpoint was recast in the language of arithmetic intersections, refined to arbitrary discriminants, and generalized to Hilbert modular surfaces and Shimura curves, where the diagonal on \(X_0(1)\times X_0(1)\) is replaced by special divisors such as Hirzebruch–Zagier divisors or quaternionic Hecke correspondences [1008.1853], [2309.17251].

## 1. Classical statement and arithmetic meaning

For negative coprime fundamental discriminants \(D_1,D_2\), with imaginary quadratic fields \(K_i=\mathbb{Q}(\sqrt{D_i})\), unit counts \(w_i=\#\mathcal{O}_{K_i}^\times\), and \(D=D_1D_2\), the classical product of singular moduli is
\[
\mathcal{J}(D_1,D_2)
\coloneqq
\prod_{\substack{[\tau_1],[\tau_2]\\ \mathrm{disc}(\tau_i)=D_i}}
\bigl(j(\tau_1)-j(\tau_2)\bigr)^{\frac{4}{w_1w_2}}.
\]
Gross–Zagier’s formula gives
\[
\mathcal{J}(D_1,D_2)^2
=
\pm
\prod_{\substack{x^2<D\\ x^2\equiv D \pmod 4}}
\mathbf{F}\!\left(\frac{D-x^2}{4}\right),
\]
where \(\mathbf{F}(m)\) is a prime power determined by genus-character data and local parity conditions. One formulation defines
\[
\mathbf{F}(m)=\prod_{\substack{nn'=m\\ n,n'>0}} n^{\epsilon(n')},
\]
with \(\epsilon\) built from the quadratic residue data of \(D_1\) and \(D_2\). The nontrivial contributions occur exactly when there is a unique prime \(\ell\mid m\) with \(\epsilon(\ell)=-1\) and \(v_\ell(m)\) odd; then \(\mathbf{F}(m)\) is a positive power of \(\ell\) [2603.08569].

An equivalent formulation, particularly natural in the Hilbert modular and Borcherds-product approach, uses the real quadratic field \(F=\mathbb{Q}(\sqrt{D})\), the biquadratic CM field \(E=\mathbb{Q}(\sqrt{D_1},\sqrt{D_2})\), and the counting function
\[
\rho(\mathfrak a)=\#\{\mathfrak A\subset \mathcal O_E: N_{E/F}(\mathfrak A)=\mathfrak a\}.
\]
Then the logarithm of the total product is expressed as a sum over totally positive \(t\in \mathcal O_F\) and primes \(\mathfrak p\) inert in \(E/F\), with weights \(\frac{1+\operatorname{ord}_{\mathfrak p}(t\mathcal O_F)}{2}\rho(t\mathfrak p^{-1})\log N(\mathfrak p)\) [1711.02983].

The formula explains striking numerical factorizations. One example is
\[
j\!\left(\frac{1+\sqrt{-43}}{2}\right)-j\!\left(\frac{1+\sqrt{-163}}{2}\right)
=
2^{19}\cdot 3^6\cdot 5^3\cdot 7^3\cdot 37\cdot 433,
\]
whose prime divisors are inert in both \(K_1=\mathbb{Q}(\sqrt{-43})\) and \(K_2=\mathbb{Q}(\sqrt{-163})\) [2309.17251]. In the coprime fundamental setting, a recurring structural principle is that only primes with tightly constrained splitting behavior contribute to the factorization.

## 2. Arithmetic intersection formulation

The modern interpretation of the formula is intersection-theoretic. Gross–Zagier’s factorization can be viewed as a computation of the arithmetic intersection of CM divisors on an integral model of the modular curve, or equivalently as the intersection of the diagonal on \(X_0(1)\times X_0(1)\) with a CM cycle. In this formulation, the non-archimedean valuation of a difference of singular moduli becomes a local intersection multiplicity at a prime \(p\) [1008.1853].

At the local level, the valuation is computed by supersingular geometry. One standard form is that
\[
\operatorname{ord}_p\!\Big(N_{H/\mathbb{Q}}(j(\tau_1)-j(\tau_2))\Big)
\]
is a sum over supersingular elliptic curves \(E/\mathbb{F}_p\) of non-negative integers counting compatible optimal embeddings of the two quadratic orders into \(\operatorname{End}(E)\), equivalently representation numbers attached to ternary quadratic forms or Gross lattices [1810.13214], [1202.6410]. This is the source of the prime-by-prime factorization.

A stack-theoretic formulation makes the intersection precise. For the Deligne–Mumford stack \(M\) of elliptic curves and finite étale CM stacks \(Y_i\to M\), the arithmetic intersection multiplicity is defined by summing lengths of strictly Henselian local rings at geometric intersection points, weighted by automorphism groups and \(\log |F_{\mathfrak p}|\) [2509.11553]. In this language, the classical quantity
\[
J(d_1,d_2)=
\Bigg(\prod_{\substack{[\tau_1],[\tau_2]\\ \operatorname{disc}(\tau_i)=d_i}}
(j(\tau_1)-j(\tau_2))\Bigg)^{\frac{4}{w_1w_2}}
\]
encodes the arithmetic degree of the CM intersection.

This perspective also clarifies a common point of confusion. The fundamental object is often not a single bare difference \(j(\tau_1)-j(\tau_2)\), but rather a norm or a product over Galois conjugates, or an intersection number on a stack. The local factors are geometric multiplicities, not merely ad hoc exponents.

## 3. Refined formulas and arbitrary discriminants

The original Gross–Zagier theorem assumed coprime fundamental discriminants. Subsequent work refined both the indexing of local contributions and the range of admissible discriminants. A decisive refinement is due to Howard–Yang, who decomposed the Hecke correspondence into zero-dimensional stacks \(X_\alpha\) indexed by totally positive \(\alpha\in F^\times\), with \(F=\mathbb{Q}(\sqrt{d_1d_2})\). If \(\alpha\in \mathfrak D^{-1}\) and \(\operatorname{Diff}(\alpha)=\{\mathfrak p\}\) is a singleton, then \(X_\alpha\) is supported in characteristic \(p\) and
\[
\deg(X_\alpha)
=
\frac{1}{2}\,\log(p)\cdot \operatorname{ord}_{\mathfrak p}(\alpha)\cdot \rho(\alpha\,\mathfrak D^{-1}),
\]
while the degree vanishes if \(\alpha\notin \mathfrak D^{-1}\) or \(\#\operatorname{Diff}(\alpha)>1\). The total Hecke intersection is recovered by summing \(\deg(X_\alpha)\) over \(\alpha\) of fixed trace [2509.11553].

Howard–Yang also proved an arithmetic Siegel–Weil identity identifying these arithmetic degrees with Fourier coefficients of the central derivative of Hecke’s normalized Eisenstein series. In the notation of that work, if
\[
E^{*\,\prime}(T,0)=\sum_a d_a q^a,
\]
then
\[
4\deg(X_a')=d_a,
\]
and \(d_a\) is explicitly \(2\,\operatorname{ord}_p(a_{\mathfrak p})\,p(a\mathfrak D_{\mathfrak p}^{-1})\log p\) when \(\operatorname{Diff}(a)=\{\mathfrak p\}\), and \(0\) otherwise [1202.6410]. This refines Gross–Zagier by separating the local intersection according to the auxiliary index \(a\), rather than only after summation over all \(a\) of a given trace.

An independent direction is the extension to arbitrary discriminants, including non-maximal orders. For distinct discriminants \(d_1,d_2\), one still has
\[
J(d_1,d_2)^{w_{d_1}w_{d_2}}
=
\pm
\prod_{\substack{x^2<d_1d_2\\ x^2\equiv d_1d_2 \bmod 4}}
F(d_1d_2-x^2),
\]
but now \(F(m)\) is allowed to be a possibly fractional prime power in the degenerate \(m=0\) case, and for \(m>0\) it is supported at a single rational prime \(\ell\) characterized by Hilbert-symbol conditions. For \(\ell>2\), the valuation \(v_\ell(F(m))\) is computed by counting invertible ideals of prescribed norm in \(\mathcal O_{d_1}\), subject to local optimality conditions [1206.6942]. This removes the maximal-order restriction and is motivated in part by denominator questions for genus-two invariants.

## 4. Hilbert modular surfaces and Yang’s generalization

A substantial generalization replaces the modular curve by a Hilbert modular surface. Let \(F=\mathbb{Q}(\sqrt D)\) with \(D\equiv 1\pmod 4\) prime, and let \(\mathcal M\) be the Hilbert modular stack parametrizing abelian surfaces with real multiplication by \(\mathcal O_F\) and a \(\partial_F^{-1}\)-polarization in the sense of Deligne–Pappas. The Hirzebruch–Zagier divisors \(\mathcal T_m\) on \(\mathcal M\) play the role of special divisors, and for \(m=1\) the divisor \(\mathcal T_1\) identifies with the diagonal embedding of the moduli stack of elliptic curves into \(\mathcal M\) [1008.1853].

Yang considers a non-biquadratic quartic CM field \(K=F(\sqrt{\Delta})\) and the associated CM cycle \(\mathrm{CM}(K)\) on \(\mathcal M\). Under the assumptions
\[
\mathcal O_K=\mathcal O_F+\mathcal O_F\cdot \frac{w+\sqrt{\Delta}}{2}
\]
and that the discriminant \(\tilde D=\Delta\Delta'\equiv 1\pmod 4\) is prime, the main theorem asserts
\[
\mathcal T_m\cdot \mathrm{CM}(K)=\frac12\,b_m
\]
for all \(m\ge 1\), where \(b_m=\sum_p b_m(p)\log p\) is defined from local data of the reflex extension \(\tilde K/\tilde F\) and the counting function
\[
\rho(\mathfrak a)=\#\{\mathfrak A\subset \mathcal O_{\tilde K}:N_{\tilde K/\tilde F}(\mathfrak A)=\mathfrak a\}.
\]
This confirms a special case of the Bruinier–Yang conjecture and makes the local structure of the intersection completely explicit [1008.1853].

The relationship with the classical Gross–Zagier formula is exact when \(F=\mathbb{Q}\). In that case \(\mathcal T_1\) is the diagonal, \(K=\mathbb{Q}(\sqrt{d_1})\oplus \mathbb{Q}(\sqrt{d_2})\), and Yang recovers
\[
\mathcal T_1\cdot \mathrm{CM}(K)
=
\sum_{\operatorname{disc}[\tau_i]=d_i}
\frac{4}{w_1w_2}\log|j(\tau_1)-j(\tau_2)|.
\]
The local intersection terms match the Gross–Zagier exponents, so the Hilbert modular-surface formula is literally a higher-dimensional generalization of the singular-moduli factorization [1008.1853].

This framework also yields a height-theoretic consequence: under the same hypotheses, Yang proves the first non-trivial non-abelian Chowla–Selberg formula, a special case of Colmez’s conjecture, namely
\[
h_{\mathrm{Fal}}(A)=\frac12\,\beta(K/F)
\]
for a principally polarized CM abelian surface of type \((K,\Phi)\) [1008.1853].

## 5. Quaternionic and Shimura-curve extensions

A distinct family of generalizations replaces modular curves by Shimura curves attached to indefinite quaternion algebras. For genus-zero Shimura curves \(X_N\) with \(N\in\{6,10,22\}\), one chooses a Hauptmodul \(j_N\) and forms a cross-ratio invariant built from CM points \(P_i\) and their Atkin–Lehner transforms. Giampietro–Darmon conjectured, and Daas proved, that the norm of this cross-ratio factors as
\[
\mathcal J_N(D_1,D_2)^{\pm 1}
=
\pm
\prod_{\substack{x^2<D\\ x^2\equiv D \pmod{4N}}}
\mathbf F\!\left(\frac{D-x^2}{4N}\right)^{\delta(x)},
\]
where \(\delta(x)\in\{\pm1\}\) depends on the residue class of \(x\) modulo \(2N\). The proof uses \(p\)-adic \(\Theta\)-functions on the Cerednik–Drinfeld uniformization of \(X_N\), together with CM theory and an analytic argument via \(p\)-adic infinitesimal deformations of Hilbert Eisenstein series and an \(R=T\) theorem [2309.17251].

An archimedean proof of the same factorization was later obtained by evaluating automorphic Green’s functions at CM points. In that approach, the logarithm of the cross-ratio of \(j_N\)-values is expressed as a linear combination of Green’s function values, and the first Fourier coefficient of a holomorphic projection of a diagonal-restricted Hilbert Eisenstein series is shown to match the finite factorization term. For \(N\in\{6,10,22\}\), the relevant cusp-form space vanishes or is annihilated by Atkin–Lehner symmetry, so the archimedean and finite contributions coincide exactly [2603.08569].

There is also an intersection-theoretic Shimura-curve generalization. For an indefinite quaternion algebra \(B/\mathbb{Q}\) with maximal order \(\mathcal O_B\), one considers the stack \(M^B\) of QM abelian surfaces, the CM divisors \(Y_j^B\), and the Hecke correspondence \(T_m\). After decomposing the CM stacks by congruence types \(\theta:\mathcal O_K\to \mathcal O_B/\mathfrak m_B\), one defines zero-dimensional stacks \(X_{\theta,\alpha}^B\). If \(\alpha\in\mathfrak D^{-1}\) is totally positive and \(\operatorname{Diff}_\theta(\alpha)=\{\mathfrak p\}\), then
\[
\deg(X_{\theta,\alpha}^B)
=
\frac12\,\log(p)\cdot \operatorname{ord}_{\mathfrak p}(\alpha)\cdot
\rho\big(\alpha\,\mathfrak a_\theta^{-1}\mathfrak D^{-1}\big),
\]
and
\[
I(T_mY_1^B,Y_2^B)
=
\sum_{\substack{\beta\in\mathfrak D^{-1},\,\beta\gg 0\\ \mathrm{Tr}_{F/\mathbb Q}(\beta)=m}}
\ \sum_{\eta:\mathcal O_K\to \mathcal O_B/\mathfrak m_B}
\deg(X_{\eta,\beta}^B).
\]
When \(B\) is split, this reduces to the Howard–Yang refinement on modular curves [2509.11553].

## 6. Consequences, related structures, and higher-dimensional context

One direct arithmetic consequence is that differences of singular moduli are never units. Using Gross–Zagier, Gross–Kohnen–Zagier, Schofer, and Bruinier–Kudla–Yang, Li gave a short proof that
\[
\log\Big|N_{H/\mathbb Q}(j(\tau_1)-j(\tau_2))\Big|>0,
\]
hence \(j(\tau_1)-j(\tau_2)\) cannot be a unit. The same method implies that for any \(m\ge 1\) there exists a prime \(\mathfrak p\) such that the reductions modulo \(\mathfrak p\) of two CM elliptic curves are \(m\)-isogenous, via the modular polynomial criterion \(\varphi_m(j(E_1),j(E_2))\equiv 0\) [1810.13214].

A related but distinct development concerns traces rather than products. In the language of generalized Jacobians with cuspidal modulus, Bruinier and Li proved that the generating series of Heegner divisor classes is a weakly holomorphic modular form of weight \(3/2\), extending Gross–Kohnen–Zagier and yielding new proofs of Zagier’s modularity of traces of singular moduli. In the classical case \(F=j-744\), this recovers
\[
g(\tau)=-q^{-1}+2+\sum_{d>0}\operatorname{Tr}_d(j-744)\,q^d
\in M_{3/2}(\Gamma_0(4)),
\]
thus placing the trace formulas and the Gross-style Jacobian interpretation in a common framework [1508.07112].

Higher-dimensional quaternionic analogues also exist over totally real fields of strict class number \(1\). In that setting, the role of optimal embeddings into \(B_{p,\infty}\) is played by simultaneous embeddings of primitive CM fields into superspecial orders in definite quaternion algebras over the totally real base field. The resulting counting formulas determine which superspecial primes can occur in factorizations of differences of Siegel modular function values at CM points and provide bounds relevant to Igusa invariants and genus-two CM constructions [1112.2009]. This suggests that the singular-moduli formula is best understood not as an isolated identity for \(j\)-values, but as the one-dimensional prototype of a much broader arithmetic intersection and automorphic framework.

Across these developments, the persistent structural principle is the same: a seemingly elementary factorization of special values is governed by local deformation theory, representation densities, and special cycles on Shimura varieties. The classical Gross–Zagier formula remains the basic example, but its later refinements show that its natural home is arithmetic geometry rather than elementary algebraic number theory alone.

Source: https://www.emergentmind.com/topics/gross-zagier-formula-on-singular-moduli