Conjectured lattice-to-continuum operator identity and spectral classification for the $\mathbb{Z}_2$-invariant inhomogeneous six-vertex model (Bazhanov–Kotousov–Koval–Lukyanov, 2021)

Bazhanov, Kotousov, Koval and Lukyanov (Nucl. Phys. B 965 (2021) 115337) analyzed the scaling limit of the critical $\mathbb{Z}_2$-invariant inhomogeneous (staggered) six-vertex model with twisted boundary conditions and proposed that its continuum limit is a non-compact conformal field theory, describable as a gauged $\mathrm{SL}(2,\mathbb{R})$ WZW model (the quantum Lorentzian black-hole nonlinear sigma model). The analytic backbone of this proposal is a conjectural instance of the ODE/IQFT correspondence: the properly normalized scaling limit of the eigenvalues of the lattice Baxter $Q$-operators $\mathbb{A}_\pm(\zeta)$ on any low-energy Bethe state is required to coincide with the connection coefficients $D_\pm(\mu)$ of a second-order differential equation carrying multiple apparent singularities whose positions are fixed by a nonlinear algebraic system, together with an operator-level version of this identity and a conjectured classification of the discrete part of the spectrum. Despite extensive numerical verification by the authors, a series of follow-up works generalizing the framework to $\mathbb{Z}_r$ symmetry (2023–2026), and rigorous proofs of the analogous correspondence in the homogeneous case (monster potentials of the quantum KdV model), all of these conjectures remain open; the most recent paper in the series (2025/2026) states explicitly that its conjectures were not proven. The problem is to prove or refute these conjectures for all low-energy states and all admissible parameter values.

Background

The six-vertex model is an exactly solvable two-dimensional lattice model of statistical mechanics whose local Boltzmann weights are parameterized by an anisotropy $q$ and a spectral parameter $\zeta$; its row-to-row transfer matrix $\mathbb{T}(\zeta)$ generates a commuting family of operators, which is what makes the model solvable by the Bethe ansatz. Baxter's 1971 generalized (inhomogeneous) six-vertex model allows the weights to carry site-dependent inhomogeneity parameters $\{\eta_J\}_{J=1}^{N}$ while preserving this integrability. Besides $\mathbb{T}(\zeta)$, the commuting family contains two Baxter $Q$-operators $\mathbb{A}_\pm(\zeta)$: their eigenvalues $A_\pm(\zeta)$ are polynomials whose zeros $\{\zeta_m\}$ (the Bethe roots) satisfy a closed system of algebraic Bethe ansatz equations, from which the transfer-matrix spectrum is reconstructed. The present problem concerns the two-site periodic (staggered) choice $\eta_J={\rm i}(-1)^{J-1}$, which endows the model with a global $\mathbb{Z}_2$ symmetry, studied with twisted boundary conditions $\omega={\rm e}^{{\rm i}\pi{\tt k}}$ on a lattice of $N$ columns. The model is critical whenever $q$ and $\omega$ are unimodular; for the $\mathbb{Z}_2$-invariant case two critical domains occur, and the source paper works in $\arg(q^2)\in(0,\pi)$, parameterized as $q={\rm e}^{{\rm i}\pi/(n+2)}$ with $n>0$. Writing $\beta^2=2/(n+2)$, so that $n=2/\beta^2-2$ and $\arg(q^2)=\pi\beta^2$, this is the domain $0<\beta^2<1$.

In the scaling limit one sends $N\to\infty$ along sequences of Bethe states (renormalization-group, or RG, trajectories) that keep the scaled energy finite. The trajectories are characterized by RG invariants: the chiral parameters $p,\bar p$ (fixed by the spin $S^z$, the twist ${\tt k}$ and the winding number ${\tt w}$), a parameter $s=\mathrm{s}\!\!\!\lim_{N\to\infty} b(N)$ built from an eigenvalue of the quasi-shift operator (the lattice translation by two sites), and, at excitation level $\mathtt{L}$, a non-ordered set $\bm w=\{w_a\}_{a=1}^{\mathtt{L}}$ constrained by a coupled nonlinear algebraic system. The ODE/IQFT (ordinary-differential-equation / integrable-quantum-field-theory) correspondence (Dorey–Tateo; Bazhanov–Lukyanov–Zamolodchikov) asserts, broadly, that scaling limits of $Q$-operator eigenvalues of critical integrable lattice models coincide with spectral data of differential equations; for the homogeneous XXZ chain the scaled Bethe roots reproduce the spectrum of a quantum anharmonic oscillator. For the $\mathbb{Z}_2$-inhomogeneous model the proposed counterpart is the second-order equation

$\Bigg[-\frac{{\rm d}^2}{{\rm d}z^2}+\frac{p^2-\tfrac14}{z^2}+\frac{2{\rm i}s}{z}+1+\sum_{a=1}^{\mathtt{L}}\bigg(\frac{2}{(z-w_a)^2}+\frac{n}{z(z-w_a)}\bigg)+\mu^{-2-n}z^n\Bigg]\,\Phi=0\,,$

which has a regular singularity at $z=0$ with exponent $p$, an irregular singularity at $z=\infty$ of degree $n$, a term carrying the continuous parameter $s$, and $\mathtt{L}$ apparent singularities (regular points at which solutions can nevertheless pick up monodromy) at the positions $w_a$. The connection coefficients $D_\pm(\mu\,|\,\bm w,p,s)$, defined as Wronskians between the two basis solutions with fixed $z\to0$ asymptotics and a distinguished subdominant solution as $z\to+\infty$, play the role of the continuum $Q$-functions. The central conjecture (eq. (11.1) of the source) is the eigenvalue identity

$\mathrm{s}\!\!\!\lim_{N\to\infty\atop b(N)\to s}G^{(N/2)}\big(-\mu^2\,\big|\,\tfrac{2}{n+2}\big)\,A_\pm\big({\rm i}\,\big(N/(2N_0)\big)^{-\frac{n}{n+2}}\mu\big)={D}_\pm(\mu\,|\,\bm w,p,s)\,,$

where $G^{(N)}$ is a state-independent normalizing factor ensuring the limit exists and $N_0$ is a known constant. Accepting (11.1), the limit becomes an operator identity (13.13): continuum operators $\mathbf a_\pm(\lambda)$, constructed explicitly as traces over a $q$-oscillator algebra of path-ordered exponentials of vertex operators in two-copy Heisenberg-algebra Fock spaces, should act diagonally on vectors $\bm\psi_{p,s}(\bm w)$ with eigenvalues $D_\pm$. Beyond the identity itself, the source conjectures a complete spectral classification (sec. 10): every low-energy state flows either to the continuous sector (real limiting $s$, forming the space $\mathcal H^{(\rm cont)}$) or to the discrete sector $\mathcal H^{(\rm disc)}$ with $s=\pm{\rm i}\mathfrak q_a$ or $\pm{\rm i}\bar{\mathfrak q}_a$, where $\mathfrak q_a=-p-\tfrac12-a$ and $\bar{\mathfrak q}_a=-\bar p-\tfrac12-a$ with $a$ an integer in an explicitly bounded range, and the number of RG trajectories at level $(\mathtt{L},\bar{\mathtt{L}})$ with a given limiting value is $d_{S^z+a}(\bar{\mathtt{L}})\,d_a(\mathtt{L})$ (resp. its barred version), where the integers $d_a(\mathtt{L})$ are defined by the generating function $\chi_a({\tt q})=({\tt q},{\tt q})_\infty^{-2}\sum_{m\ge0}(-1)^m{\tt q}^{ma+m(m+1)/2}=\sum_{\mathtt{L}\ge0}d_a(\mathtt{L})\,{\tt q}^{\mathtt{L}}$.

The physical stakes are the identity of the continuum limit. Ikhlef–Jacobsen–Saleur (2012) conjectured that the staggered $\mathbb{Z}_2$ chain flows to the Euclidean black-hole sigma model $\mathrm{SL}(2,\mathbb{R})/U(1)$, the exact conformal field theory of Witten's two-dimensional black hole; its partition function indeed matches half of the lattice scaling-limit partition function. Bazhanov–Kotousov–Koval–Lukyanov rejected this identification on unitarity grounds (the Euclidean black hole is unitary, the lattice model is not) and proposed instead the Lorentzian black-hole nonlinear sigma model, a gauged $\mathrm{SL}(2,\mathbb{R})$ Wess–Zumino–Witten (WZW) model — a two-dimensional conformal field theory built from group-valued fields with a topological coupling — with specific boundary and reality conditions, whose space of states carries an extended conformal algebra $\bar W_\infty\otimes W_\infty$ with central charge $c<2$ and a continuous component in its spectrum of conformal dimensions. The source paper is explicit that this identification is tentative ('The local CFT, if it exists, ...') precisely because the correspondence (11.1)/(13.13) underlying it is unproven, resting on analytic continuation arguments, sum-rule checks at low orders, and extensive numerics. The source additionally records further unresolved issues: closed analytic forms for the non-universal norm constants $C(\beta)$, $C_0(\beta)$ (introduced earlier for the homogeneous chain by Kotousov–Lukyanov) and $C_0^{(\rm alt)}$; the exceptional loci $n=\tfrac{2}{2k-1}$ (equivalently $\beta^2=1-\tfrac1{2k}$) where the continuum operators require renormalization; and the fact that the continuous and discrete sectors inherit two inequivalent conjugations from the lattice, so that whether they are sectors of one local quantum field theory remains unresolved.

Problem Statement

Prove or refute the conjectured ODE/IQFT correspondence formulated by Bazhanov–Kotousov–Koval–Lukyanov for the scaling limit of the critical $\mathbb{Z}_2$-invariant inhomogeneous six-vertex model with twisted boundary conditions, in the parameter domain $0<\beta^2<1$ (i.e. $q={\rm e}^{{\rm i}\pi/(n+2)}$ with $\arg(q^2)=\pi\beta^2\in(0,\pi)$, $n=2/\beta^2-2>0$) with generic twist ${\tt k}$ and spin sector $S^z\ge0$. The correspondence comprises three linked claims, to be established for every low-energy Bethe state (RG trajectory) characterized by the scaling data $p,\bar p,s,\bm w,\bar{\bm w}$:

(1) Eigenvalue identity (eq. (11.1) of the source). The normalized scaling limit

$\mathrm{s}\!\!\!\lim_{N\to\infty\atop b(N)\to s}G^{(N/2)}\big(-\mu^2\,\big|\,\tfrac{2}{n+2}\big)\,A_\pm\big({\rm i}\,\big(N/(2N_0)\big)^{-\frac{n}{n+2}}\mu\big)$

exists and equals the connection coefficient $D_\pm(\mu\,|\,\bm w,p,s)$ of the ODE (11.2) with the normalizations (11.3)–(11.5) of the source, for all levels $\mathtt{L},\bar{\mathtt{L}}\ge0$, all admissible $s$, and all $n>0$, including a consistent treatment of the exceptional loci $n=\tfrac{2}{2k-1}$ where the normalizing factor is modified and the continuum operators require renormalization.

(2) Operator identity (eq. (13.13)). The continuum operators $\mathbf a_\pm(\lambda)$, defined through the $q$-oscillator traces of path-ordered vertex-operator exponentials over the Heisenberg Fock spaces (eqs. (13.9)–(13.10) of the source), admit eigenvectors $\bm\psi_{p,s}(\bm w)$ with eigenvalues exactly equal to $D_\pm(\mu\,|\,\bm w,p,s)$ (with the $\lambda$–$\mu$ relation (11.10) of the source), so that the lattice scaling limit realizes a diagonalizable commuting family on the conformal towers.

(3) Spectral classification (Conjectures (I), (II), (A), (B), (III) of sec. 10). Every low-energy Bethe state flows either to the continuous sector (real limiting $s$) or to the discrete sector with $s=\pm{\rm i}\mathfrak q_a$ or $\pm{\rm i}\bar{\mathfrak q}_a$, where $\mathfrak q_a=-p-\tfrac12-a$ and $\bar{\mathfrak q}_a=-\bar p-\tfrac12-a$ with $a$ an integer in the bounded range of eq. (10.27); and the number of RG trajectories at level $(\mathtt{L},\bar{\mathtt{L}})$ with a given limiting value is exactly $d_{S^z+a}(\bar{\mathtt{L}})\,d_a(\mathtt{L})$ (resp. $d_a(\bar{\mathtt{L}})\,d_{S^z+a}(\mathtt{L})$), with $d_a$ defined by the generating function (10.25) — establishing in particular existence and uniqueness of the RG trajectory for each admissible data set.

These three components form one scientific objective — establishing the conjectured continuum description of the scaling limit — since (2) is the operator-level content of (1) and (3) is the completeness statement of the same correspondence. An answer may be a proof of (1)–(3) in full (a precisely delimited partial proof is acceptable only if it states exactly which states, levels and parameters remain unsettled), or a rigorous refutation: an explicitly identified state and parameter values inside the conjecture's domain for which the scaling limit in (1) fails to exist or differs from $D_\pm$, or for which the classification (3) fails (a low-energy state whose $b(N)$ has a limit outside the predicted set, or degeneracies violating the counting formula). The exact analytic determination of the norm constants $C(\beta)$, $C_0(\beta)$, $C_0^{(\rm alt)}$ and the construction of the Lorentzian black-hole QFT as a complete local theory are related open issues recorded in the source, but lie outside the completion boundary of this problem.

References

1. V. V. Bazhanov, G. A. Kotousov, S. M. Koval, S. L. Lukyanov, Scaling limit of the $\mathbb{Z}_2$ invariant inhomogeneous six-vertex model, Nucl. Phys. B 965 (2021) 115337, doi:10.1016/j.nuclphysb.2021.115337, https://arxiv.org/abs/2010.10613
2. V. V. Bazhanov, G. A. Kotousov, S. M. Koval, S. L. Lukyanov, On the scaling behaviour of the alternating spin chain, JHEP 08 (2019) 087, doi:10.1007/JHEP08(2019)087, https://arxiv.org/abs/1903.05033
3. V. V. Bazhanov, G. A. Kotousov, S. L. Lukyanov, Equilibrium density matrices for the 2D black hole sigma models from an integrable spin chain, JHEP 03 (2021) 169, doi:10.1007/JHEP03(2021)169, https://arxiv.org/abs/2010.10603
4. V. V. Bazhanov, G. A. Kotousov, S. M. Koval, S. L. Lukyanov, Some Algebraic Aspects of the Inhomogeneous Six-Vertex Model, SIGMA 17 (2021) 025, doi:10.3842/SIGMA.2021.025, https://arxiv.org/abs/2010.10615
5. R. J. Baxter, Generalized ferroelectric model on a square lattice, Stud. Appl. Math. 50 (1971) 51-69, doi:10.1002/sapm197150151, https://doi.org/10.1002/sapm197150151
6. P. Dorey, R. Tateo, Anharmonic oscillators, the thermodynamic Bethe ansatz, and nonlinear integral equations, J. Phys. A: Math. Gen. 32 (1999) L419-L425, doi:10.1088/0305-4470/32/38/102, arXiv:hep-th/9812211, https://arxiv.org/abs/hep-th/9812211
7. V. V. Bazhanov, S. L. Lukyanov, A. B. Zamolodchikov, Spectral determinants for Schrödinger equation and Q operators of conformal field theory, J. Stat. Phys. 102 (2001) 567-576, doi:10.1023/A:1004838616921, arXiv:hep-th/9812247, https://arxiv.org/abs/hep-th/9812247
8. Y. Ikhlef, J. Lykke Jacobsen, H. Saleur, An Integrable spin chain for the SL(2,R)/U(1) black hole sigma model, Phys. Rev. Lett. 108 (2012) 081601, doi:10.1103/PhysRevLett.108.081601, https://arxiv.org/abs/1109.1119
9. E. Witten, String theory and black holes, Phys. Rev. D 44 (1991) 314, doi:10.1103/PhysRevD.44.314, https://doi.org/10.1103/PhysRevD.44.314
10. G. A. Kotousov, S. L. Lukyanov, Bethe state norms for the Heisenberg spin chain in the scaling limit, Nucl. Phys. B 947 (2019) 114748, doi:10.1016/j.nuclphysb.2019.114748, https://arxiv.org/abs/1906.07081
11. G. A. Kotousov, S. L. Lukyanov, On the scaling behaviour of an integrable spin chain with $\mathbb{Z}_r$ symmetry, Nucl. Phys. B 993 (2023) 116269, doi:10.1016/j.nuclphysb.2023.116269, https://arxiv.org/abs/2305.03620
12. S. Gehrmann, G. A. Kotousov, S. L. Lukyanov, Scaling limit of the ground state Bethe roots for the inhomogeneous XXZ spin-$\frac12$ chain, Nucl. Phys. B 1006 (2024) 116624, doi:10.1016/j.nuclphysb.2024.116624, https://arxiv.org/abs/2406.12102
13. G. A. Kotousov, S. L. Lukyanov, D. A. Shabetnik, Scaling limit of Q-functions for the $\mathbb{Z}_r$ invariant inhomogeneous XXZ spin-$\frac12$ chain near free fermion point, Nucl. Phys. B 1024 (2026) 117315, doi:10.1016/j.nuclphysb.2026.117315, https://arxiv.org/abs/2511.20530
14. M. Azhari, A. Klümper, Managing singular kernels and logarithmic corrections in the staggered six-vertex model, JHEP 12 (2024) 040, doi:10.1007/JHEP12(2024)040, https://arxiv.org/abs/2406.09889
15. H. Frahm, S. Gehrmann, G. A. Kotousov, Scaling limit of the staggered six-vertex model with $U_q(\mathfrak{sl}(2))$ invariant boundary conditions, SciPost Phys. 16 (2024) 149, doi:10.21468/SciPostPhys.16.6.149, https://arxiv.org/abs/2312.11238
16. H. Frahm, S. Gehrmann, Finite size spectrum of the staggered six-vertex model with $U_q(\mathfrak{sl}(2))$-invariant boundary conditions, JHEP 01 (2022) 070, doi:10.1007/JHEP01(2022)070, https://arxiv.org/abs/2111.00850
17. R. Conti, D. Masoero, Counting monster potentials, JHEP 02 (2021) 059, doi:10.1007/JHEP02(2021)059, https://arxiv.org/abs/2009.14638
18. R. Conti, D. Masoero, On Solutions of the Bethe Ansatz for the Quantum KdV Model, Commun. Math. Phys. 402 (2023) 335-390, doi:10.1007/s00220-023-04728-3, https://arxiv.org/abs/2112.14625
19. D. Masoero, G. Ruzza, ODE/IM Correspondence at the Free-Fermion Point. Laguerre Wronskians, Shifted Symmetric Functions, and Quantum KdV, arXiv:2605.24563 (2026), https://arxiv.org/abs/2605.24563