Chekhov–Eynard–Orantin Recursion

• Chekhov–Eynard–Orantin Recursion is a residue-based method that constructs symmetric meromorphic multidifferentials from spectral curve data.
• The technique systematically processes unstable inputs using a kernel and residue sums to produce correlators applicable in enumerative geometry, Hurwitz theory, and quantum curves.
• Its versatility in handling simple and higher-order branch points underpins applications in matrix models, topological recursion, and operator formalisms within mathematical physics.

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 (Lewanski et al., 2015). 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 (Dumitrescu et al., 2013).

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 $\omega_{g,n}$0, the requirement that $\omega_{g,n}$1 have only simple zeroes, a local involution near each branch point characterized by equality of $\omega_{g,n}$2-values on the two local sheets, and a symmetric bidifferential $\omega_{g,n}$3 whose only singularity is a double pole on the diagonal with leading behavior $\omega_{g,n}$4 in local coordinates (Lewanski et al., 2015). In genus $\omega_{g,n}$5, one often takes $\omega_{g,n}$6 and $\omega_{g,n}$7 (Fang et al., 2014).

In another standard presentation, especially for genus-zero spectral curves arising from Laplace-transform constructions, one introduces a Lagrangian immersion $\omega_{g,n}$8 with tautological form $\omega_{g,n}$9, and defines

$2g-2+n>0$0

which is holomorphic on the diagonal $2g-2+n>0$1 (Mulase et al., 2012). 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 $2g-2+n>0$2-plane (Mulase, 2012).

The plane-curve formulation is not the only one. For Hitchin fibrations, the spectral curve is replaced by a smooth curve $2g-2+n>0$3 over a projective base curve $2g-2+n>0$4, with $2g-2+n>0$5 given by the restriction of the Liouville form and $2g-2+n>0$6 the normalized fundamental bidifferential on $2g-2+n>0$7 (Dumitrescu et al., 2013). Local spectral curves with simple ramification also underlie the partition functions of semi-simple cohomological field theories and their Givental–Teleman descriptions (Alexandrov, 2022).

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

The unstable initial correlators are

$2g-2+n>0$8

and the standard recursion kernel is

$2g-2+n>0$9

where $W_{0,1}=y\,dx$0 denotes the local conjugate point with the same $W_{0,1}=y\,dx$1-value (Lewanski et al., 2015). Equivalent formulas appear throughout the literature, including formulations using $W_{0,1}=y\,dx$2 for the local involution near a branch point $W_{0,1}=y\,dx$3 or $W_{0,1}=y\,dx$4 (Branahl et al., 2022).

For all stable pairs $W_{0,1}=y\,dx$5, the CEO recursion is the universal residue formula

$W_{0,1}=y\,dx$6

where the prime excludes unstable $W_{0,1}=y\,dx$7-terms (Lewanski et al., 2015). In genus-zero and matrix-model notations the same formula is written for $W_{0,1}=y\,dx$8, with residues at the branch points $W_{0,1}=y\,dx$9 (Branahl et al., 2022).

The first stable case already exhibits the structure of the recursion. For $W_{0,2}=B$0, the only nonzero contribution is the product of two $W_{0,2}=B$1's, so $W_{0,2}=B$2 is a sum of branch-point residues of $W_{0,2}=B$3 (Lewanski et al., 2015). For $W_{0,2}=B$4, one obtains the standard genus-one residue formula involving $W_{0,2}=B$5 (Branahl et al., 2022). These low-topology cases are repeatedly used as base computations in ribbon-graph models, Hurwitz theories, dynamical triangulations, and quantum-curve constructions (Chapman et al., 2010).

3. Laplace transform, unstable geometries, and rational examples

One important construction starts from unstable enumerative data. The disk and annulus contributions $W_{0,2}=B$6 and $W_{0,2}=B$7 are Laplace transformed to define $W_{0,2}=B$8 and $W_{0,2}=B$9, after which one sets

$(\Sigma,x,y,B)$0

The resulting spectral curve, bidifferential, and kernel then reproduce the full CEO correlators under Laplace transform (Dumitrescu et al., 2012). 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 $(\Sigma,x,y,B)$1 (Dumitrescu et al., 2012).

Several rational spectral curves became standard benchmarks. For the higher-genus Catalan numbers, the curve is

$(\Sigma,x,y,B)$2

and the Laplace-transformed generating functions satisfy EO recursion on this curve (Mulase et al., 2012). For single Hurwitz numbers, the Lambert curve

$(\Sigma,x,y,B)$3

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 (Mulase et al., 2012).

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 $(\Sigma,x,y,B)$4, the CEO recursion appears on

$(\Sigma,x,y,B)$5

with branch points at $(\Sigma,x,y,B)$6 (Mulase et al., 2010). 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 $(\Sigma,x,y,B)$7 with the rational parametrization

$(\Sigma,x,y,B)$8

and this framework yields a new proof of Kontsevich’s constants (Chapman et al., 2010).

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,x,y,B)$9

the coefficients in the asymptotic expansion of $dx$0 at the puncture recover the intersection numbers

$dx$1

and for $dx$2 these coincide with the orbifold Hurwitz intersection theory, reproducing the Johnson–Pandharipande–Tseng formula (Lewanski et al., 2015).

Equivariant mirror symmetry for $dx$3 supplies another precise realization. For the equivariantly perturbed mirror Landau–Ginzburg model with superpotential

$dx$4

the EO recursion on the corresponding affine curve encodes all genus all descendants equivariant Gromov–Witten invariants of $dx$5; the non-equivariant limit yields the Norbury–Scott conjecture, and the large-radius limit recovers the Bouchard–Mariño conjecture on simple Hurwitz numbers (Fang et al., 2014).

The recursion also governs orbifold and monotone variants of Hurwitz theory. For monotone $dx$6-orbifold Hurwitz numbers, the spectral curve

$dx$7

has simple branch points given by the zeros of $dx$8, and the connected generating differentials satisfy exactly the CEO recursion on this curve (Kramer et al., 2019). Deformations of spectral curves further lead to ELSV-type formulae and vanishing relations for integrals of generalized Hodge classes $dx$9; the key deformation property is

$x$0

which turns the variation of $x$1 into residue insertions and makes polynomiality in deformation parameters translate into vanishing of negative-power coefficients (Borot et al., 2023).

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 $x$2, proves the cut-and-join equation

$x$3

and derives $x$4 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 (Alexandrov, 2022).

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 $x$5 copies of Virasoro constraints for the total ancestor potential (Milanov, 2012). This equivalence is one of the routes by which CEO recursion enters the Givental–Teleman description of semi-simple theories (Alexandrov, 2022).

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 (Branahl et al., 2022). 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 $x$6 and at $x$7 in higher genus, requiring blobbed topological recursion rather than pure CEO recursion (Branahl et al., 2022).

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 $x$8; 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 (Fuji et al., 16 Dec 2025). In a non-commutative direction, higher Airy structures and $x$9-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 (Borot et al., 2021).

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 (Bychkov et al., 2019). This is the framework used for the Bousquet–Mélou–Schaeffer numbers on the curve $\Sigma$0, where the unique critical point has ramification index $\Sigma$1 (Bychkov et al., 2019).

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

$\Sigma$2

(Andersen et al., 2017). 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 (Mulase et al., 2012). For Hitchin fibrations, CEO recursion on a smooth curve in $\Sigma$3 produces free energies $\Sigma$4, a WKB wave function

$\Sigma$5

and a rank-two $\Sigma$6-module whose semiclassical limit recovers the original Hitchin spectral curve (Dumitrescu et al., 2013). The harmonic-oscillator curve $\Sigma$7 gives a particularly explicit example in which CEO correlators reconstruct the WKB expansion of the corresponding Schrödinger equation (Cutimanco et al., 2017).

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 $\Sigma$8, and the CEO recursion solves the Schwinger–Dyson equations order by order in the genus expansion (Fuji et al., 11 Dec 2025). 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 (Fuji et al., 23 Sep 2025).

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 (Branahl et al., 2022).