---
title: Chekhov–Eynard–Orantin Recursion
url: https://www.emergentmind.com/topics/chekhov-eynard-orantin-recursion
type: topic
---

# Chekhov–Eynard–Orantin Recursion

Chekhov–Eynard–Orantin recursion, also called Eynard–Orantin or CEO topological recursion, is a residue recursion on a spectral curve that constructs a family of symmetric meromorphic multidifferentials \(W_{g,n}\) or \(\omega_{g,n}\) for all stable pairs \(2g-2+n>0\) from the unstable data \(W_{0,1}=y\,dx\) and \(W_{0,2}=B\). In its standard form, the input is a quadruple \((\Sigma,x,y,B)\) consisting of a Riemann surface, two meromorphic functions, and a symmetric bidifferential with the diagonal double pole of the Bergman kernel; the recursion then sums residues at the simple zeros of \(dx\), using the local involution that exchanges the two sheets of the \(x\)-projection [1504.07439]. The same formalism appears in enumerative geometry, matrix and field theories, mirror symmetry, singularity theory, quantum curves, and Hamiltonian or cut-and-join descriptions, while several works also emphasize that its standard form is tied to simple ramification and must be modified in the presence of higher criticality or extra poles [1310.6022].

## 1. Spectral-curve data and the standard set-up

A spectral curve for the CEO recursion consists of a Riemann surface \(\Sigma\), meromorphic functions \(x,y:\Sigma\to\mathbb C\) or \(\mathbb P^1\), the requirement that \(dx\) have only simple zeroes, a local involution near each branch point characterized by equality of \(x\)-values on the two local sheets, and a symmetric bidifferential \(B(z_1,z_2)\) whose only singularity is a double pole on the diagonal with leading behavior \(dz_1\,dz_2/(z_1-z_2)^2\) in local coordinates [1504.07439]. In genus \(0\), one often takes \(\Sigma\cong \mathbb P^1\) and \(B(z_1,z_2)=dz_1\,dz_2/(z_1-z_2)^2\) [1411.3557].

In another standard presentation, especially for genus-zero spectral curves arising from Laplace-transform constructions, one introduces a Lagrangian immersion \(\iota:\Sigma\to T^*\mathbb C_x\) with tautological form \(\eta=y\,dx\), and defines
\[
W_{0,2}(z_1,z_2)=B_\Sigma(z_1,z_2)-\frac{dx(z_1)\,dx(z_2)}{(x(z_1)-x(z_2))^2},
\]
which is holomorphic on the diagonal \(z_1=z_2\) [1210.3006]. This difference of convention is part of the formalism rather than a contradiction: several rational-curve implementations use the Bergman kernel itself as the unstable bidifferential, while other expositions subtract the polar part in the \(x\)-plane [1210.2106].

The plane-curve formulation is not the only one. For Hitchin fibrations, the spectral curve is replaced by a smooth curve \(E\subset T^*C\) over a projective base curve \(C\), with \(W_{0,1}\) given by the restriction of the Liouville form and \(W_{0,2}\) the normalized fundamental bidifferential on \(E\times E\) [1310.6022]. Local spectral curves with simple ramification also underlie the partition functions of semi-simple cohomological field theories and their Givental–Teleman descriptions [2202.09090].

## 2. Kernel, residue formula, and low-topology structure

The unstable initial correlators are
\[
W_{0,1}(z)=y(z)\,dx(z), \qquad W_{0,2}(z_1,z_2)=B(z_1,z_2),
\]
and the standard recursion kernel is
\[
K(z_0;z)=\frac{1}{2\,\bigl(y(z)-y(\bar z)\bigr)\,dx(z)}\int_{\bar z}^{z}B(z_0,\xi),
\]
where \(\bar z\) denotes the local conjugate point with the same \(x\)-value [1504.07439]. Equivalent formulas appear throughout the literature, including formulations using \(\sigma_i(z)\) for the local involution near a branch point \(\beta_i\) or \(a_i\) [2205.12166].

For all stable pairs \(2g-2+n>0\), the CEO recursion is the universal residue formula
\[
\begin{aligned}
W_{g,n}(z_1,\dots,z_n)
=\sum_i \operatorname*{Res}_{z\to p_i} K(z_1;z)\Biggl[
&\,W_{g-1,n+1}(z,\bar z,z_2,\dots,z_n)\\
&+\sum_{\substack{g_1+g_2=g\\ I\sqcup J=\{2,\dots,n\}}}^{\prime}
W_{g_1,|I|+1}(z,z_I)\,
W_{g_2,|J|+1}(\bar z,z_J)
\Biggr],
\end{aligned}
\]
where the prime excludes unstable \((0,1)\)-terms [1504.07439]. In genus-zero and matrix-model notations the same formula is written for \(\omega_{g,n}\), with residues at the branch points \(\beta_i\) [2205.12166].

The first stable case already exhibits the structure of the recursion. For \((g,n)=(0,3)\), the only nonzero contribution is the product of two \(W_{0,2}\)'s, so \(W_{0,3}\) is a sum of branch-point residues of \(K\,W_{0,2}\,W_{0,2}\) [1504.07439]. For \((g,n)=(1,1)\), one obtains the standard genus-one residue formula involving \(B(q,\sigma(q))\) [2205.12166]. These low-topology cases are repeatedly used as base computations in ribbon-graph models, Hurwitz theories, dynamical triangulations, and quantum-curve constructions [1009.2055].

## 3. Laplace transform, unstable geometries, and rational examples

One important construction starts from unstable enumerative data. The disk and annulus contributions \(D_{0,1}\) and \(D_{0,2}\) are Laplace transformed to define \(F_{0,1}\) and \(F_{0,2}\), after which one sets
\[
x=x(t),\qquad y(t)\,dx(t)=dF_{0,1}(w(t)),\qquad
W_{0,2}(t_1,t_2)=d_{w_1}d_{w_2}F_{0,2}(w(t_1),w(t_2)).
\]
The resulting spectral curve, bidifferential, and kernel then reproduce the full CEO correlators under Laplace transform [1202.1159]. This construction is used for Grothendieck’s dessins d’enfants, intersection numbers of tautological cotangent classes, single Hurwitz numbers, and stationary Gromov–Witten invariants of \(\mathbb{CP}^1\) [1202.1159].

Several rational spectral curves became standard benchmarks. For the higher-genus Catalan numbers, the curve is
\[
x=z+\frac1z,\qquad y=-z,
\]
and the Laplace-transformed generating functions satisfy EO recursion on this curve [1210.3006]. For single Hurwitz numbers, the Lambert curve
\[
x=z\,e^{-z},\qquad y=z
\]
plays the same role, and the principal-specialized partition function satisfies both a heat-type PDE and a holonomic difference-differential equation whose semiclassical limit reproduces the curve [1210.3006].

Ribbon-graph models supply another family of examples. The Poincaré polynomials of the combinatorial moduli space of curves are the Laplace transform of counts of Grothendieck’s dessins d’enfants, and after the change of variables \(e^{-w}=(t+1)/(t-1)\), the CEO recursion appears on
\[
x(t)=\frac{(t+1)^2}{4t},\qquad y(t)=\frac{t}{t+1}
\]
with branch points at \(t=\pm1\) [1009.2135]. The Euclidean and symplectic volumes of the combinatorial moduli space of pointed smooth algebraic curves also satisfy an EO-type recursion, now on the curve \(x=y+1/y\) with the rational parametrization
\[
x(t)=2+\frac{4}{t^2-1},\qquad y(t)=\frac{t+1}{t-1},
\]
and this framework yields a new proof of Kontsevich’s constants [1009.2055].

## 4. Intersection theory, Hurwitz theory, and mirror symmetry

A central application is the reconstruction of intersection numbers on moduli spaces of curves. For the Chiodo-type spectral curve
\[
\Sigma=\mathbb{CP}^1,\qquad x(z)=-z^r+\ln z,\qquad y(z)=z^s,
\]
the coefficients in the asymptotic expansion of \(W_{g,n}\) at the puncture recover the intersection numbers
\[
\int_{\mathcal M_{g,n}} c\bigl(-R^\bullet\pi_*\mathcal S\bigr)\,\psi_1^{d_1}\cdots\psi_n^{d_n},
\]
and for \(s=r\) these coincide with the orbifold Hurwitz intersection theory, reproducing the Johnson–Pandharipande–Tseng formula [1504.07439].

Equivariant mirror symmetry for \(\mathbb P^1\) supplies another precise realization. For the equivariantly perturbed mirror Landau–Ginzburg model with superpotential
\[
W(Y)=Y+\frac{q}{Y}+w_1\ln Y+w_2\ln\!\bigl(\tfrac{q}{Y}\bigr),
\]
the EO recursion on the corresponding affine curve encodes all genus all descendants equivariant Gromov–Witten invariants of \(\mathbb P^1\); the non-equivariant limit yields the Norbury–Scott conjecture, and the large-radius limit recovers the Bouchard–Mariño conjecture on simple Hurwitz numbers [1411.3557].

The recursion also governs orbifold and monotone variants of Hurwitz theory. For monotone \(q\)-orbifold Hurwitz numbers, the spectral curve
\[
x(z)=z(1-z^q),\qquad
y(z)=\frac{z^{q-1}}{1-z^q}
\]
has simple branch points given by the zeros of \(x'(z)=1-(q+1)z^q\), and the connected generating differentials satisfy exactly the CEO recursion on this curve [1909.02302]. Deformations of spectral curves further lead to ELSV-type formulae and vanishing relations for integrals of generalized Hodge classes \(Q^{(r,s)}\); the key deformation property is
\[
\partial_t\bigl(y_t\,dx_t\bigr)=d f_t(z),
\]
which turns the variation of \(\omega_{g,n}\) into residue insertions and makes polynomiality in deformation parameters translate into vanishing of negative-power coefficients [2312.16049].

## 5. Loop equations, cut-and-join formalisms, and operator realizations

The CEO recursion admits several algebraic reformulations. For local spectral curves with simple ramification points, the partition function of topological recursion has a cubic cut-and-join operator description. In that framework one constructs a formal series \(Z_S\), proves the cut-and-join equation
\[
\frac{\partial Z}{\partial \hbar}=\widehat W\,Z,
\]
and derives \(N\) families of Virasoro constraints together with a deformed dimension constraint; these constraints imply the cut-and-join description, and for semi-simple cohomological field theories the resulting partition functions lie in the same CEO family [2202.09090].

A related formulation appears for total ancestor potentials in singularity theory. There, period integrals and phase forms replace the plane-curve data, and the EO recursion is proved equivalent to \(N\) copies of Virasoro constraints for the total ancestor potential [1211.5847]. This equivalence is one of the routes by which CEO recursion enters the Givental–Teleman description of semi-simple theories [2202.09090].

Loop equations provide a second major route. In the quartic LSZ model, Dyson–Schwinger equations imply abstract loop equations, and together with the pole structure this yields the CEO recursion for the corresponding meromorphic differentials [2205.12166]. The same paper stresses a useful limitation: in the complex LSZ model no extra poles appear away from the ramification points, so the standard CEO recursion applies directly, whereas in the hermitian Grosse–Wulkenhaar case analytic continuation produces poles on the anti-diagonals \(z=-z_j\) and at \(z=0\) in higher genus, requiring blobbed topological recursion rather than pure CEO recursion [2205.12166].

Operator and Hamiltonian realizations make the same mechanism explicit. A string-field Hamiltonian formalism introduces creation and annihilation operators, defines correlation functions as vacuum expectation values, and derives loop equations from commutators \([\,\hat H,\psi(x)\,]\); the genus-zero Schwinger–Dyson equation then recovers the classical spectral curve, while higher-order loop equations reorganize into residues at branch points, exactly matching CEO topological recursion [2512.14059]. In a non-commutative direction, higher Airy structures and \(\mathcal W\)-algebras lead to a non-commutative topological recursion computing Whittaker and Gaiotto vectors, with recovery of the classical CEO formalism in suitable self-dual limits [2104.04516].

## 6. Generalizations, limitations, and related recursions

The standard CEO formalism assumes simple ramification points. When the spectral curve has a higher-order critical point, the local involution is replaced by the full local Galois orbit of sheets, and one uses the Bouchard–Eynard generalization rather than the ordinary simple-branch recursion [1908.04147]. This is the framework used for the Bousquet–Mélou–Schaeffer numbers on the curve \(x=z/(1+z)^m\), where the unique critical point has ramification index \(m-1\) [1908.04147].

A broader geometric framework is geometric recursion. Starting from pairs of pants excision on bordered surfaces and integrating the resulting amplitudes over moduli spaces with the Weil–Petersson measure, one obtains Mirzakhani-type recursions; after Laplace transform, these become EO topological recursion on
\[
x(z)=\frac{z^2}{2},\qquad
y(z)=-\frac{\sin(2\pi z)}{2\pi},\qquad
\omega_{0,2}(z_1,z_2)=\frac{dz_1\,dz_2}{(z_1-z_2)^2}
\]
[1711.04729]. This suggests that CEO recursion can appear as the Laplace image of a more geometric excision formalism.

Quantum-curve constructions are another persistent theme. For genus-zero spectral curves such as those of generalized Catalan numbers and single Hurwitz numbers, the principal specialization of the partition function satisfies a Schrödinger equation whose classical symbol is exactly the spectral curve [1210.3006]. For Hitchin fibrations, CEO recursion on a smooth curve in \(T^*C\) produces free energies \(F_{g,n}\), a WKB wave function
\[
\Psi(z,\hbar)=\exp\!\Bigl[\sum_{m=0}^{\infty}\hbar^{m-1}S_m(z)\Bigr],
\]
and a rank-two \(D\)-module whose semiclassical limit recovers the original Hitchin spectral curve [1310.6022]. The harmonic-oscillator curve \(y^2=x^2-c^2\) gives a particularly explicit example in which CEO correlators reconstruct the WKB expansion of the corresponding Schrödinger equation [1701.08913].

The same residue technology now appears in models of two-dimensional gravity. Multicritical dynamical triangulations and causal dynamical triangulations are governed by one-cut and two-cut spectral curves on \(\mathbb P^1\), and the CEO recursion solves the Schwinger–Dyson equations order by order in the genus expansion [2512.10519]. A parallel reformulation for pure gravity on dynamical triangulations treats basic, strip, and continuum models; in the double-scaling limit both discrete models flow to the same continuum pure-gravity curve, and no modification of the recursion is needed beyond the change of spectral data [2509.18916].

In this range of applications, a common misconception is that every appearance of loop equations or every ramified spectral curve leads directly to the pure CEO formalism. The literature instead separates three cases: standard CEO recursion for simple ramification and poles only at branch points, blobbed topological recursion when extra poles survive analytic continuation, and Bouchard–Eynard recursion for higher-order critical points [2205.12166].

Source: https://www.emergentmind.com/topics/chekhov-eynard-orantin-recursion