First-principles ODE/IQFT derivation and the quantization-condition conjecture for the staggered six-vertex model with $U_q(\mathfrak{sl}(2))$-invariant open boundary conditions

The staggered (inhomogeneous) six-vertex model in its critical regime, equipped with integrable open boundary conditions that leave it invariant under the quantum group $U_q(\mathfrak{sl}(2))$, has a scaling limit whose spectrum of conformal dimensions contains a continuous component, making it a leading lattice realization of the non-compact conformal field theories associated with two-dimensional Euclidean and Lorentzian black hole backgrounds. Frahm, Gehrmann and Kotousov (SciPost Phys. 16, 149 (2024)) showed that the entire low-energy scaling limit — the density of the continuous spectrum, a discrete set of levels, and the scaling partition function — is organized by an ODE/IQFT-type correspondence in which a Schrödinger-type equation with extra singularities, constrained by an algebraic system, has a spectral determinant that reproduces the lattice Baxter Q-function data along renormalization-group trajectories. The central statement, a one-to-one correspondence between lattice RG trajectories and solutions of a quantization condition, is a conjecture supported by numerical data at lattice size $L=10$; the authors state they know of no way to rigorously prove the underlying scaling relation, the counting of admissible singularity configurations rests on a non-rigorous limiting argument, and the auxiliary assumptions on the pole/zero structure of the phase shift were verified numerically only for descendant level $d\le 3$. The open problem is to establish this framework from first principles: derive the scaling relation with its exact prefactors, prove or refute the quantization-condition conjecture, and establish the analytic structure of the auxiliary algebraic system for general level $d$. No resolution is known in the literature to date: follow-up work by the same research groups remains at the level of conjecture plus numerics, and rigorous results exist only for adjacent, simpler systems.

Background

The six-vertex model assigns arrows (equivalently spin-$\frac{1}{2}$ degrees of freedom) to the edges of a square lattice with Boltzmann weights that conserve the arrow flux through each vertex; its anisotropy is parametrized by a parameter $q$. In the staggered (inhomogeneous) variant the weights alternate between the two sublattices, breaking translation symmetry to a $\mathbb{Z}_2$ sublattice symmetry. Throughout this problem the staggered inhomogeneities are fixed at the model's self-dual point, the standard choice for this family. Staggered six-vertex models remain integrable — they possess an infinite family of commuting transfer matrices — and in the critical regime $q=\mathrm{e}^{\mathrm{i}\pi/(n+2)}$ with $n>0$ (so $|q|=1$ and $0<\arg q<\pi/2$; the analysis excludes the exceptional values $n=\frac{2}{2k-1}$, $k=1,2,\ldots$, which require separate treatment) they exhibit a remarkable phenomenon: although each site carries a finite-dimensional degree of freedom, the scaling limit develops a non-compact degree of freedom and the spectrum of conformal dimensions (the scaling exponents that label the states of the continuum-limit theory) acquires a continuous component. For (quasi-)periodic boundary conditions the model was conjectured by Ikhlef, Jacobsen and Saleur to realize the two-dimensional Euclidean black hole sigma model, i.e. the $SL(2,\mathbb{R})/U(1)$ Wess–Zumino–Witten coset ("cigar") conformal field theory arising in two-dimensional string theory; Bazhanov, Kotousov, Koval and Lukyanov subsequently developed a detailed ODE/IQFT description of its scaling limit, solved the spectral problem of that CFT, and refined the conjecture to include a Lorentzian-signature sector.

The ODE/IM correspondence (between ordinary differential equations and integrable models; in its quantum-field-theoretic form, ODE/IQFT), discovered by Dorey and Tateo and by Bazhanov, Lukyanov and Zamolodchikov, states that spectral determinants of certain Schrödinger-type ordinary differential equations — generalizations of the anharmonic oscillator — satisfy the same functional ($TQ$) relations as the Baxter $Q$-operators of integrable quantum field theories and lattice models. Here the Baxter $Q$-operator is a commuting operator whose eigenvalues encode the Bethe roots, the auxiliary variables parametrizing eigenstates in the Bethe-ansatz solution of the transfer matrix, and the spectral determinant is a function built from the subdominant solutions of the ODE that vanishes exactly at its discrete eigenvalues. The correspondence turns the lattice scaling limit into a problem of one-dimensional spectral theory.

The present problem concerns the open-boundary version studied by Robertson, Jacobsen and Saleur, and by Frahm and Gehrmann. Open integrable boundary conditions for vertex models are introduced through Sklyanin's reflection formalism, in which the boundary weights are encoded in $K$-matrices satisfying reflection equations. A particular one-parameter family of diagonal $K$-matrices depends on a boundary parameter $\epsilon$; its $U_q(\mathfrak{sl}(2))$-invariant member at $\epsilon=0$ leaves the staggered chain invariant under the quantum group $U_q(\mathfrak{sl}(2))$ — the $q$-deformation of the $\mathfrak{sl}(2)$ symmetry algebra, whose finite-dimensional irreducible multiplets $\mathcal{M}_{\mathcal{S}}$ of spin $\mathcal{S}=0,1,2,\ldots$ organize the spectrum; the ground state has non-zero spin (its value depends on the anisotropy) and so spontaneously breaks this symmetry. Working in the Hamiltonian picture with a chain of $2L$ sites, an RG trajectory $\{|{}_L\rangle\}$ is a sequence of $U_q(\mathfrak{sl}(2))$ multiplets, one per (even) lattice size, tracked to large $L$ by means of the $Q$-operator, for which the source paper found a new finite matrix-element formula (valid in the sector $S^{z}=0$ of zero total spin projection, for the one-parameter diagonal $K$ family, with no infinite sums and for generic complex $q$ and $\epsilon$). Along a trajectory the energy obeys

$E \asymp L\,e_{\infty}+f_{\infty}+\frac{\pi v_{\mathrm{F}}}{L}\Big(\frac{p^{2}}{n+2}+\frac{b^{2}}{n}-\frac{1}{12}+d\Big)+o(L^{-1-\varepsilon})\,,$

where $2p=2\mathcal{S}+1-(n+2)$, $d=0,1,2,\ldots$ is a descendant level, and the nontrivial quantity is $b(L)=\frac{n}{2\pi}\log\sqrt{B}$, with $B$ the eigenvalue of the quasi-shift operator — a conserved charge commuting with the transfer matrix (related to the lowest transfer matrix of the hierarchy) that plays the role of lattice quasi-momentum. In the scaling limit $b(L)\to s$ with $s$ real for the continuous part of the spectrum and $s$ pure imaginary for a discrete set of levels; the conformal dimensions are $\Delta-\frac{c}{24}=\frac{p^{2}}{n+2}+\frac{s^{2}}{n}-\frac{1}{12}+d$ with central charge $c=\frac{2(n-1)}{n+2}<2$ (the number characterizing the conformal field theory).

The ODE/IQFT description identifies the normalized product $A(\zeta)=\prod_{j=1}^{L-\mathcal{S}}(1-\zeta/\zeta_{j})(1-\zeta\zeta_{j})$ over the Bethe roots with the spectral determinant $D(\mu)$ of

$\Big[-\frac{\mathrm{d}^{2}}{\mathrm{d}z^{2}}+\frac{p^{2}-\frac{1}{4}}{z^{2}}+\frac{2\mathrm{i}s}{z}+1+\sum_{a=1}^{d}\Big(\frac{2}{(z-w_{a})^{2}}+\frac{n}{z(z-w_{a})}\Big)+\mu^{-2-n}z^{n}\Big]\psi=0\,,$

a Schrödinger equation with inverse-square and Coulomb-type terms, an exponential potential, and — for $d>0$ — extra inverse-square singularities at complex points $z=w_{a}$. Single-valuedness of solutions near $z=w_{a}$ forces the $w_{a}$ to obey a coupled algebraic system of $d$ equations with cubic singularities; for generic parameters the admissible solution sets $\boldsymbol{w}=(w_{1},\ldots,w_{d})$ are finite in number, predicted to be $N(d\,|\,\mathcal{S})=\mathrm{par}_{2}(d)-\mathrm{par}_{2}(d-2\mathcal{S}-1)$, where $\mathrm{par}_{2}(d)$ is the number of bipartitions of the integer $d$ (partitions of $d$ into parts of two colours). The identification itself, $\lim_{L\to\infty,\,b(L)\to s}G^{(L)}(-\mu^{2})\,A((L/L_{0})^{-\frac{n}{n+2}}\mathrm{i}\mu)=D(\mu)$ with an explicit divergent-sum renormalization factor $G^{(L)}$ and scale $L_{0}$, is asserted without proof; the authors write that they "do not know of a way of rigorously proving" it, its support being numerical sum-rule checks inherited from the quasi-periodic analysis.

From the large-$L$ Bethe-root asymptotics the authors derive the quantization condition

$\Big(\frac{L}{L_{0}}\Big)^{2\mathrm{i}s}\,\mathrm{e}^{\frac{\mathrm{i}}{2}\delta(\boldsymbol{w},p,s)}\Big|_{s=b(L)}=\sigma+O(L^{-\varepsilon})\,,\qquad \sigma=(-1)^{L-\mathcal{S}}\,,$

where the phase shift $\delta$ is the logarithm of a ratio of coefficients $\mathfrak{C}^{(\pm)}_{p,s}(\boldsymbol{w})$ given by exact determinant formulae of Kotousov and Lukyanov depending on $p$, $s$ and the solution set $\boldsymbol{w}$. The paper's central result, stated as a Conjecture, is a one-to-one correspondence: every RG trajectory labeled by $(\mathcal{S},d,\boldsymbol{w})$ yields a $b(L)$ obeying the quantization condition, and conversely every solution $b_{*}(L)$ of the quantization condition (correction terms ignored) in the strip $-\frac{n}{2}<\operatorname{Im}b_{*}\le\frac{n}{2}$ corresponds to a unique $U_q(\mathfrak{sl}(2))$ multiplet, with $\sqrt{B}-\exp(\frac{2\pi}{n}b_{*}(L))$ decaying faster than any power of $\log L$. Accepting the conjecture, the paper derives the density of states of the continuous spectrum, the discrete spectrum (imaginary $s=\pm s_{a}$, $s_{a}=\mathrm{i}(\frac{n}{2}-\mathcal{S}-a)$ with $0\le a+\mathcal{S}<\frac{n}{2}$), the decomposition of the scaling-limit space of states into irreducible representations of the $\mathcal{W}_{\infty}$ algebra (the algebra of extended conformal symmetry generated by infinitely many higher-spin currents), and the scaling partition function $Z^{(\mathrm{scl})}=Z^{(\mathrm{cont})}+Z^{(\mathrm{disc})}$.

What is not established: the conjecture is supported only by comparison with exact-diagonalization data at $L=10$; the scaling relation between Bethe-root products and the spectral determinant lacks proof; and the assumptions on the phase shift's poles and zeros (that they are simple, and that poles and zeros coming from different solution sets do not coincide) — on which the multiplicities $N(d\,|\,\mathcal{S})$ and $N_{a}^{\pm}$ in the density-of-states and character formulae depend — were verified numerically only for $d\le 3$. In addition, the computed $Z^{(\mathrm{scl})}$ does not match the published brane (boundary-state) partition functions of the 2D black hole CFTs (in particular those of Ribault and Schomerus), a mismatch the authors flag as requiring a separate investigation.

Problem Statement

Establish, from first principles, the ODE/IQFT description of the scaling limit of the staggered six-vertex model with $U_q(\mathfrak{sl}(2))$-invariant open boundary conditions at the self-dual point, in the critical regime $q=\mathrm{e}^{\mathrm{i}\pi/(n+2)}$ with $n>0$ (excluding the exceptional values $n=\frac{2}{2k-1}$, $k=1,2,\ldots$). Concretely, the following linked results must be proved — or the central one rigorously refuted:

1. Scaling relation. Starting from the lattice Hamiltonian and its Bethe-ansatz / Baxter $Q$-operator solution (diagonal one-parameter $K$-matrix family), prove the scaling relation identifying the normalized Bethe-root product $A(\zeta)$ along an arbitrary RG trajectory $\{|{}_L\rangle\}$ with the spectral determinant $D(\mu)$ of the ODE $\Big[-\frac{\mathrm{d}^{2}}{\mathrm{d}z^{2}}+\frac{p^{2}-\frac{1}{4}}{z^{2}}+\frac{2\mathrm{i}s}{z}+1+\sum_{a=1}^{d}\big(\frac{2}{(z-w_{a})^{2}}+\frac{n}{z(z-w_{a})}\big)+\mu^{-2-n}z^{n}\Big]\psi=0$, with $2p=2\mathcal{S}+1-(n+2)$ and the singularity positions $\boldsymbol{w}$ constrained by the associated algebraic system, including the exact normalization factors $G^{(L)}$ and $L_{0}$ and the branch conventions for $b(L)=\frac{n}{2\pi}\log\sqrt{B}$.

2. Quantization-condition conjecture. Prove or refute the conjecture of Frahm, Gehrmann and Kotousov: (a) for any RG trajectory labeled by the $U_q(\mathfrak{sl}(2))$ spin $\mathcal{S}$, the descendant level $d$ and a solution set $\boldsymbol{w}$ of the algebraic system, the quantity $b(L)$ computed from the quasi-shift eigenvalue obeys the quantization condition $\big(\frac{L}{L_{0}}\big)^{2\mathrm{i}s}\mathrm{e}^{\frac{\mathrm{i}}{2}\delta(\boldsymbol{w},p,s)}\big|_{s=b(L)}=\sigma+O(L^{-\varepsilon})$, $\sigma=(-1)^{L-\mathcal{S}}$; and (b) conversely, every solution $b_{*}(L)$ of this condition (correction terms ignored) with $-\frac{n}{2}<\operatorname{Im}b_{*}\le\frac{n}{2}$ corresponds to a unique $U_q(\mathfrak{sl}(2))$ multiplet whose quasi-shift eigenvalue satisfies $\sqrt{B}-\exp\big(\frac{2\pi}{n}b_{*}(L)\big)\to 0$ faster than any power of $\log L$.

3. Analytic structure of the auxiliary system. Establish, for general descendant level $d$: the counting formula $N(d\,|\,\mathcal{S})=\mathrm{par}_{2}(d)-\mathrm{par}_{2}(d-2\mathcal{S}-1)$ for the admissible solution sets of the algebraic system, so far supported only by a non-rigorous limiting argument; and the simplicity of the poles and zeros of the phase-shift function $\mathrm{e}^{\frac{\mathrm{i}}{2}\delta(\boldsymbol{w},p,s)}$ together with their non-coincidence across different solution sets, so far verified numerically only for $d\le 3$ — the properties on which the multiplicities in the density-of-states and discrete-spectrum formulae rest.

The problem is answered by a complete derivation of items 1–3, or by a rigorous refutation of the conjecture in item 2 (for example, an explicit RG trajectory violating the quantization condition beyond its stated corrections, or a quantization-condition solution in the strip with no corresponding multiplet). Two further open directions recorded in the source paper — extending the finite Q-operator matrix-element formula beyond the $S^{z}=0$ sector and the one-parameter diagonal $K$ family, and identifying the scaling partition function $Z^{(\mathrm{scl})}$ with boundary partition functions of the 2D black hole CFTs — are separable follow-up questions and do not belong to the completion boundary of this problem.

References

1. H. Frahm, S. Gehrmann and G. A. Kotousov, Scaling limit of the staggered six-vertex model with $U_q(\mathfrak{sl}(2))$ invariant boundary conditions, SciPost Phys. 16, 149 (2024), DOI 10.21468/SciPostPhys.16.6.149, arXiv:2312.11238, https://arxiv.org/abs/2312.11238
2. V. V. Bazhanov, G. A. Kotousov, S. M. Koval and S. L. Lukyanov, Scaling limit of the ${\cal Z}_2$ invariant inhomogeneous six-vertex model, Nucl. Phys. B 965 (2021) 115337, DOI 10.1016/j.nuclphysb.2021.115337, arXiv:2010.10613, https://arxiv.org/abs/2010.10613
3. H. Frahm and S. Gehrmann, Finite size spectrum of the staggered six-vertex model with $U_q(\mathfrak{sl}(2))$-invariant boundary conditions, J. High Energy Phys. 01 (2022) 070, DOI 10.1007/JHEP01(2022)070, arXiv:2111.00850, https://arxiv.org/abs/2111.00850
4. N. F. Robertson, J. L. Jacobsen and H. Saleur, Lattice regularisation of a non-compact boundary conformal field theory, J. High Energy Phys. 02 (2021) 180, DOI 10.1007/JHEP02(2021)180, arXiv:2012.07757, https://arxiv.org/abs/2012.07757
5. Y. Ikhlef, J. L. Jacobsen and 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, arXiv:1109.1119, https://arxiv.org/abs/1109.1119
6. G. A. Kotousov and S. L. Lukyanov, Spectrum of the reflection operators in different integrable structures, J. High Energy Phys. 02 (2020) 029, DOI 10.1007/JHEP02(2020)029, arXiv:1910.05947, https://arxiv.org/abs/1910.05947
7. S. Ribault and V. Schomerus, Branes in the 2D black hole, J. High Energy Phys. 02 (2004) 019, DOI 10.1088/1126-6708/2004/02/019, arXiv:hep-th/0310024, https://arxiv.org/abs/hep-th/0310024
8. P. Dorey and R. Tateo, Anharmonic oscillators, the thermodynamic Bethe ansatz, and nonlinear integral equations, J. Phys. A 32 (1999) L419, DOI 10.1088/0305-4470/32/38/102, arXiv:hep-th/9812211, https://arxiv.org/abs/hep-th/9812211
9. V. V. Bazhanov, S. L. Lukyanov and A. B. Zamolodchikov, Spectral determinants for Schrödinger equation and Q-operators of conformal field theory, J. Stat. Phys. 102 (2001) 567, DOI 10.1023/A:1004838616921, arXiv:hep-th/9812247, https://arxiv.org/abs/hep-th/9812247
10. S. Gehrmann, G. A. Kotousov and S. L. Lukyanov, Scaling limit of the ground state Bethe roots for the inhomogeneous XXZ spin-$\frac{1}{2}$ chain, Nucl. Phys. B 1006 (2024) 116624, DOI 10.1016/j.nuclphysb.2024.116624, arXiv:2406.12102, https://arxiv.org/abs/2406.12102
11. G. A. Kotousov, S. L. Lukyanov and D. A. Shabetnik, Scaling limit of $Q$-functions for the ${\cal Z}_r$ invariant inhomogeneous XXZ spin-$\frac{1}{2}$ chain near free fermion point, Nucl. Phys. B 1024 (2026) 117315, DOI 10.1016/j.nuclphysb.2026.117315, arXiv:2511.20530, https://arxiv.org/abs/2511.20530
12. H. Frahm and S. Gehrmann, Finite-size spectrum of the staggered six-vertex model with antidiagonal boundary conditions, Nucl. Phys. B 1006 (2024) 116655, DOI 10.1016/j.nuclphysb.2024.116655, arXiv:2405.20919, https://arxiv.org/abs/2405.20919
13. H. Frahm, S. Gehrmann, R. I. Nepomechie and A. L. Retore, Quantum-group-invariant $D^{(2)}_{n+1}$ models: Bethe ansatz and finite-size spectrum, J. High Energy Phys. 12 (2025) 117, DOI 10.1007/JHEP12(2025)117, arXiv:2509.00610, https://arxiv.org/abs/2509.00610
14. P. Lu, J. Cao, W.-L. Yang, I. Marquette and Y.-Z. Zhang, Exact physical quantities of the $D^{(2)}_{2}$ spin chain model with generic open boundary conditions, J. High Energy Phys. 05 (2025) 137, DOI 10.1007/JHEP05(2025)137, arXiv:2502.09262, https://arxiv.org/abs/2502.09262
15. M. Azhari and A. Klümper, Managing singular kernels and logarithmic corrections in the staggered six-vertex model, J. High Energy Phys. 12 (2024) 040, DOI 10.1007/JHEP12(2024)040, arXiv:2406.09889, https://arxiv.org/abs/2406.09889
16. R. Conti and D. Masoero, On solutions of the Bethe Ansatz for the Quantum KdV model, Commun. Math. Phys. 402 (2023) 335, DOI 10.1007/s00220-023-04728-3, arXiv:2112.14625, https://arxiv.org/abs/2112.14625
17. G. Degano, ODE/IM correspondence in the semiclassical limit: large degree asymptotics of the spectral determinants for the ground state potential, arXiv:2409.07866, https://arxiv.org/abs/2409.07866
18. S. Gehrmann, The influence of boundary conditions on the emergence of non-compact degrees of freedom in the scaling limit of integrable lattice models, PhD thesis, Leibniz Universität Hannover (2024), DOI 10.15488/17537, https://inspirehep.net/literature/2805471
19. Y. Liu, H. Shimizu, D. Liu and K. Kawabata, Extracting Boundary Conformal Data from Periodic Non-Hermitian Critical Chains, arXiv:2606.16785, https://arxiv.org/abs/2606.16785