ODE/IM Correspondence

• ODE/IM Correspondence is an integrable framework that equates distinguished differential equation solutions and associated spectral determinants with analytic objects like Baxter T and Q functions.
• It employs techniques such as WKB approximations, Stokes multipliers, and fusion identities to map the spectral data to Bethe ansatz equations and quantum integrable structures.
• The framework extends to diverse settings—including conformal field theories, affine Toda and supersymmetric models—and connects with gauge theories and geometric minimal surfaces.

The ODE/IM correspondence is the exact equivalence between spectral determinants associated with certain generalized Sturm–Liouville problems and the Baxter $T$ and $Q$ functions emerging within the Bethe Ansatz framework, and more generally the identification of spectral data of a differential equation with eigenvalues of commuting operators in an integrable quantum field theory (Dorey et al., 2019). In the conformal setting it links ordinary differential equations, Stokes multipliers, connection coefficients, and WKB periods to the integrable structure of quantum KdV and $W$-algebra conformal field theories; in its off-critical and generalized forms it extends to modified affine Toda equations, thermodynamic Bethe ansatz systems, minimal surfaces in $AdS$, supersymmetric models, line defects, and quantum spectral curves (Ito et al., 2013, Dorey et al., 2019, Degano, 2024).

1. Conceptual dictionary

At the level of its basic dictionary, the correspondence identifies distinguished solutions of a differential equation with the analytic objects of an integrable model. On the ODE side one works with local solutions near regular singular points, subdominant solutions near irregular singularities, Wronskians, Stokes multipliers, and spectral determinants. On the IM side one encounters Baxter $\mathcal Q$-operators, transfer matrices, $T$- and $Y$-systems, Bethe ansatz equations, and commuting local or non-local integrals of motion (Dorey et al., 2019, Degano, 2024).

A standard example is the anharmonic oscillator

$$\frac{d^2\psi}{dx^2}=\left(x^{2\alpha}+\frac{\ell(\ell+1)}{x^2}-E\right)\psi,$$

viewed on the universal cover of $\mathbb C^\*$. In the quantum-KdV realization, one introduces Frobenius solutions near $x=0$, denoted $Q$0, and Sibuya solutions near $Q$1, denoted $Q$2, and defines spectral determinants by the Wronskians

$Q$3

These $Q$4 are exactly the objects whose zeros encode the central connection problem, and under the ODE/IM dictionary they correspond to the ground-state eigenvalues of the quantum $Q$5-operators, while the Stokes multiplier corresponds to the transfer-matrix eigenvalue (Degano, 2024).

The same pattern recurs in higher-rank and higher-spin settings. The distinguished ODE solutions are expanded in local bases, the expansion coefficients become $Q$6-functions, and Wronskian or Plücker identities among these solutions become the functional relations of the integrable model. A concise summary, already explicit in the modified affine Toda framework, is

$Q$7

(Ito et al., 2013).

2. Differential-equation structures and analytic mechanisms

The ODE side is organized by Stokes geometry. For the anharmonic oscillator, the irregular singularity at infinity produces Stokes sectors, sectorial subdominant solutions, and Stokes phenomena under analytic continuation. In the complex WKB treatment, the relevant reduced potential is the Langer-modified expression

$Q$8

and the basic WKB approximant is

$Q$9

The fundamental theorem of the WKB approximation on curves is then formulated through Volterra-type integral equations, with admissible curves selected by the quadratic differential $W$0 (Degano et al., 10 Jan 2025).

This analytic framework does more than produce formal asymptotics. It supplies basis-independent global coordinates on the space of connection data. The lectures on complex WKB introduce asymptotic values, cross-ratio coordinates of Fock–Goncharov type, and a distinguished coordinate $W$1 encoding the spectral condition; for the anharmonic oscillator,

$W$2

A WKB approximation of such coordinates is given by exponential period integrals of $W$3 along suitable cycles, making the spectral problem into a question about periods of the quadratic differential (Degano et al., 10 Jan 2025).

In higher-rank systems the same role is played by $W$4-systems, Wronskians, and Symanzik rotations. In the modified $W$5 affine Toda case, wedge-product relations among subdominant solutions in the vector and spinor representations yield algebraic identities among the associated $W$6-functions, and evaluating these identities at zeros gives Bethe ansatz equations. The resulting $W$7- and $W$8-systems are of $W$9-type, with monodromy around the origin controlling the boundary conditions (Ito et al., 2016). Numerical work on untwisted affine Lie algebras makes this mechanism explicit by representing $AdS$0-functions as inner products of a dual local solution with the subdominant solution and computing their zeros by Cheng’s algorithm (Ito et al., 2020).

A central symmetry on the ODE side is discrete rotational covariance. In the large-degree study of the quantum-KdV anharmonic oscillator, the quantum monodromy or Dorey–Tateo symmetry

$AdS$1

underlies the functional equations for $AdS$2-functions and transfer matrices (Degano, 2024). In broader ODE/IM settings, Symanzik rotation plays the same structural role by relating sectorial solutions and enforcing the functional relations that become Bethe equations or $AdS$3-systems (Ito et al., 2013, Ito et al., 2017).

3. Modified affine Toda origins and algebraic generalizations

A major structural formulation begins from affine Toda field theory rather than from a guessed ODE. For an affine Lie algebra $AdS$4, the modified affine Toda equation is obtained by a conformal transformation and a shift by the co-Weyl vector, introducing a prescribed function $AdS$5 into the affine-root term. The nonlinear equation is equivalent to the flatness condition of a Lax connection $AdS$6, and in the conformal limit this linear problem collapses to an ODE or pseudo-ODE whose spectral data reproduce the Bethe equations of the associated integrable model (Ito et al., 2013).

In the simply-laced $AdS$7 case this reduction produces a scalar differential operator of generalized Gel'fand–Dikii type. In the $AdS$8 case it yields a pseudo-differential equation, and for non-simply-laced algebras the relevant operator is naturally associated with the Langlands dual affine algebra. The identifications

$AdS$9

are explicit examples of this rule (Ito et al., 2013). This Langlands-dual appearance is not an auxiliary feature; it is the algebra naturally associated with the integrable-model Bethe equations.

The massive generalization replaces the conformal ODE by the linear problem of a modified affine Toda field equation at finite mass scale. For $\mathcal Q$0, the $\mathcal Q$1-system of subdominant solutions yields Bethe ansatz equations, from which one derives nonlinear integral equations and an ultraviolet effective central charge. With trivial monodromy around the origin, the ultraviolet limit reproduces the effective central charge of non-unitary $\mathcal Q$2 minimal models (Ito et al., 2018). The modified $\mathcal Q$3 theory extends the massive correspondence to a non-simply-laced affine algebra, with a $\mathcal Q$4-/$\mathcal Q$5-system related to $\mathcal Q$6 and monodromy-dependent boundary data (Ito et al., 2016).

Supersymmetric extensions replace ordinary affine Lie algebras by affine Lie superalgebras with purely odd simple root systems. After superconformal reduction and bosonic projection, the linear problem typically splits into a couple of ODEs, one for each bosonic subalgebra. In the special case $\mathcal Q$7, the resulting second-order equation has squared potential,

$\mathcal Q$8

which is related to $\mathcal Q$9 supersymmetric minimal models (Ito et al., 2022). A more recent treatment diagonalizes the full $T$0-type linear problem, computes WKB periods and non-local conserved quantities, and verifies that the local periods match Neveu–Schwarz-sector local integrals of motion up to sixth order (Tanabe, 16 Apr 2026).

Exceptional Lie algebras also fit into this pattern. For $T$1, the ODE side is a first-order linear system in the $T$2-dimensional fundamental representation. Its diagonalized WKB periods agree, up to sixth order, with the eigenvalues of local integrals of motion in the $T$3 conformal field theory, after the same kind of parameter identification familiar from the $T$4 and $T$5 cases (Ide et al., 9 Apr 2026).

4. Semiclassical regimes, period integrals, and special-function limits

The semiclassical side of the correspondence is governed by WKB period integrals, turning points, and asymptotic zero distributions. In the large-degree analysis of

$T$6

the limit $T$7 turns the potential $T$8 into a sharply localized barrier near $T$9, and the spectral determinants admit explicit limiting forms. When $Y$0 and $Y$1 remain in bounded domains, the central spectral determinants converge uniformly to expressions written in terms of Bessel functions of order $Y$2, and their zeros converge to the corresponding Bessel zeros (Degano, 2024).

In the scaled regime

$Y$3

with $Y$4, the determinant becomes highly oscillatory and admits a WKB cosine asymptotic with error $Y$5. The real zeros in an interval accumulate with continuous density

$Y$6

Near the turning point $Y$7, under the scaling

$Y$8

the spectral determinant converges instead to an Airy-function expression, and the zeros converge to Airy zeros (Degano, 2024). These Bessel, oscillatory, and Airy regimes provide a concrete description of how an ODE-side $Y$9-function degenerates into classical special-function asymptotics in a singular semiclassical limit.

The same WKB-period logic persists in higher rank. For $$\frac{d^2\psi}{dx^2}=\left(x^{2\alpha}+\frac{\ell(\ell+1)}{x^2}-E\right)\psi,$$0 and $$\frac{d^2\psi}{dx^2}=\left(x^{2\alpha}+\frac{\ell(\ell+1)}{x^2}-E\right)\psi,$$1 conformal field theories, the higher-order ODEs obtained from $$\frac{d^2\psi}{dx^2}=\left(x^{2\alpha}+\frac{\ell(\ell+1)}{x^2}-E\right)\psi,$$2 and $$\frac{d^2\psi}{dx^2}=\left(x^{2\alpha}+\frac{\ell(\ell+1)}{x^2}-E\right)\psi,$$3 modified affine Toda theories admit WKB expansions along a Pochhammer contour, and the coefficients $$\frac{d^2\psi}{dx^2}=\left(x^{2\alpha}+\frac{\ell(\ell+1)}{x^2}-E\right)\psi,$$4 are interpreted as classical conserved densities of the Drinfeld–Sokolov hierarchy. The period integrals $$\frac{d^2\psi}{dx^2}=\left(x^{2\alpha}+\frac{\ell(\ell+1)}{x^2}-E\right)\psi,$$5 agree with vacuum eigenvalues of local integrals of motion up to sixth order, with $$\frac{d^2\psi}{dx^2}=\left(x^{2\alpha}+\frac{\ell(\ell+1)}{x^2}-E\right)\psi,$$6 data further computed up to eighth order (Ito et al., 2024). For the $$\frac{d^2\psi}{dx^2}=\left(x^{2\alpha}+\frac{\ell(\ell+1)}{x^2}-E\right)\psi,$$7 case, the same strategy yields vacuum and first excited-state eigenvalues of the quantum Boussinesq charges (Ashok et al., 2024).

For Virasoro and $$\frac{d^2\psi}{dx^2}=\left(x^{2\alpha}+\frac{\ell(\ell+1)}{x^2}-E\right)\psi,$$8 systems, a systematic WKB algorithm makes the relation operational. Starting from the oper

$$\frac{d^2\psi}{dx^2}=\left(x^{2\alpha}+\frac{\ell(\ell+1)}{x^2}-E\right)\psi,$$9

one expands

$\mathbb C^\*$0

and interprets the coefficients $\mathbb C^\*$1 as meromorphic one-forms on the WKB curve. Period integrals of these one-forms reproduce eigenvalues of local integrals of motion in Virasoro, $\mathbb C^\*$2, and $\mathbb C^\*$3 theories, with explicit checks in Argyres–Douglas minimal models (Kudrna et al., 28 Aug 2025).

5. Conformal field theory, gauge theory, and geometry

On the conformal-field-theory side, the ODE/IM correspondence equates WKB periods or spectral determinants with eigenvalues of local integrals of motion. In the $\mathbb C^\*$4 theory, the spins $\mathbb C^\*$5 of the $\mathbb C^\*$6-currents mirror the Casimir data appearing in the $\mathbb C^\*$7 WKB expansion, and the quantities $\mathbb C^\*$8, $\mathbb C^\*$9, and $x=0$0 match the corresponding integrals of motion $x=0$1, $x=0$2, and $x=0$3 up to overall factors (Ide et al., 9 Apr 2026). For $x=0$4 and $x=0$5, the vacuum eigenvalues on the cylinder match higher-order WKB periods of the associated higher-order ODEs through sixth order (Ito et al., 2024). For $x=0$6, vacuum and level-one excited-state data extracted from the ODE coincide with quantum Boussinesq charges reconstructed from torus correlators and Zhu recursion (Ashok et al., 2024).

The correspondence is equally natural in gauge-theoretic settings. For Argyres–Douglas theories in the Nekrasov–Shatashvili limit, the quantum Seiberg–Witten curve becomes an ODE such as

$x=0$7

in the $x=0$8 case, or higher-order analogues for $x=0$9 theories. The Stokes data of these ODEs yield $Q$00-type or more general $Q$01- and $Q$02-systems, and the $Q$03-functions are exponentials of quantum-corrected Seiberg–Witten periods. The exact Bohr–Sommerfeld condition is written directly in terms of $Q$04-functions, and for $Q$05-type theories the associated integrable model is identified with a non-unitary coset model

$Q$06

(Ito et al., 2017).

A distinct but closely related off-critical branch of the subject is geometric. In the review on geometric aspects, the modified sinh-Gordon equation

$Q$07

is presented as the Gauss–Mainardi–Codazzi equation for minimal surfaces in $Q$08. The associated linear problem produces Stokes data whose $Q$09- and $Q$10-systems lead to TBA equations, with WKB periods giving mass parameters. For polygonal boundaries determined by a polynomial $Q$11, the regularized area of the minimal surface equals the free energy of the TBA system, linking ODE/IM to minimal surfaces, Wilson loops, and strong-coupling amplitude problems (Dorey et al., 2019).

The same off-critical logic survives irrelevant deformations. For the sine-Gordon/modified sinh-Gordon correspondence, a $Q$12 deformation acts quantum mechanically by the standard spectral flow of finite-volume energies and classically by a dynamical change of coordinates. Both deformations satisfy the same Burgers-type equation

$Q$13

with momentum fixed, so the ODE/IM identification of conserved quantities persists along the flow (Aramini et al., 2022).

6. Exact identities, completeness, and computational directions

Recent work has extended the correspondence from spectral matching to exact algebraic control of global spectral data. For $Q$14-symmetric and Hermitian oscillators unified by the ODE/IM framework, the fusion hierarchy of the $Q$15 $Q$16-system generates exact sum rules and zeta generating formulas for sector-dependent spectral zeta functions

$Q$17

The exact sum rules are polynomial relations among zeta values within a fixed sector, while the zeta generating formulas give explicit sector-to-sector maps between entire towers $Q$18. Their structure is governed by a selection rule controlled by the $Q$19 Symanzik symmetry together with Chebyshev-polynomial reductions of the fusion identities (Kamata, 8 Aug 2025).

Representation-theoretic refinements also appear. Guided by ODE/IM, conjectural polynomial relations among $Q$20-characters of $Q$21-modules have been proposed for quantum affine algebras of types $Q$22, $Q$23, $Q$24, and $Q$25. On the ODE side these relations mirror $Q$26-system identities for canonical solutions of a Langlands-dual connection; on the representation-theoretic side they take the form of discrete Wronskian or Casorati determinant identities among formal series built from asymptotic Kirillov–Reshetikhin characters (Sun, 2012).

A particularly strong exact result is available at the free-fermion point $Q$27. There the ODE/IM correspondence for quantum KdV is proved to be complete: every admissible solution of the Bethe or $Q$28 system is realized by a monodromy-free rational extension of the harmonic oscillator, and conversely every such ODE solution comes from the $Q$29 system. The ODE side is classified by Laguerre Wronskians and pairs of partitions, while on the IM side the first three Hamiltonians are diagonalized explicitly, with eigenstates given by Schur functions and eigenvalues by shifted symmetric functions on partitions (Masoero et al., 23 May 2026).

Computationally, the subject now combines exact analysis, WKB recursion, and numerical methods. For untwisted affine Lie algebras, Cheng’s algorithm together with a dual-linear-problem formula

$Q$30

allows direct numerical determination of Bethe roots, with agreement against nonlinear integral equations for simply-laced types including $Q$31, $Q$32, and $Q$33, and with non-simply-laced cases obtained by folding Dynkin diagrams (Ito et al., 2020). In higher-spin and $Q$34 theories, formal WKB algorithms now extract local-integral-of-motion eigenvalues directly from Bethe roots and period integrals, making the ODE/IM dictionary substantially more explicit (Kudrna et al., 28 Aug 2025).

Taken together, these developments show that the ODE/IM correspondence is not a single relation between one Schrödinger equation and one Bethe ansatz system. It is a broad integrable framework in which Stokes data, Wronskians, period integrals, fusion identities, and spectral zeta data organize the spectra of conformal, massive, supersymmetric, and defect integrable models into a common analytic structure (Dorey et al., 2019, Degano, 2024).