The quantum inverse problem — reconstructing the local operators of an integrable model in terms of algebraic Bethe ansatz monodromy-matrix entries, the step that makes form-factor expansions of correlation functions rigorous — is solved for fundamental lattice models and has been demonstrated for the one-component (gl(2)) Bose gas, but it remains open for continuum quantum gases solvable by the nested Bethe ansatz, such as the one-dimensional spin-1/2 Fermi gas (Gaudin–Yang model, of gl(2|1) symmetry in the algebraic Bethe ansatz). As a consequence, the finite-volume expansion theorem and the LeClair–Mussardo series for nested Bethe ansatz systems derived by Hutsalyuk, Pozsgay and Pristyák (2021) rest on an unproven assumption: that form factors of generic local operators in these gases satisfy the same residue relations as monodromy-matrix entries, an analogy supported only by a scaling-limit argument and the gl(2) case. This problem asks for a solution of the quantum inverse problem for local operators of nested continuum gases — an explicit reconstruction of the fields, density bilinears and their products via monodromy-matrix entries — or an equivalent rigorous proof (or refutation) of the assumed residue relations, thereby putting the nested LeClair–Mussardo expansion on firm ground.

Background

Many one-dimensional quantum many-body models are exactly solvable by the Bethe ansatz: their energy eigenstates are labeled by sets of rapidities satisfying algebraic Bethe equations. Solvability of the spectrum does not by itself give correlation functions. These are assembled from form factors — matrix elements of local operators (operators acting near a point xx, such as a density or a field) between eigenstates. In the algebraic Bethe ansatz (ABA) one organizes the Hilbert space through a monodromy matrix T(u)T(u), a matrix-valued operator whose entries Tij(u)T_{ij}(u) obey the RTT commutation relations (the quadratic exchange relations fixed by the model's R-matrix); Bethe eigenstates are created from a reference state by the off-diagonal entries. For multi-component models (several species or spin) the ABA is nested: one level of rapidities carries the momenta and further levels carry the internal degrees of freedom. The Gaudin–Yang model — the one-dimensional gas of spin-1/2 fermions with repulsive contact interaction, solved by Yang (1967) and by Gaudin — is the paradigmatic nested model; in the ABA it is described by the graded (supersymmetric) algebra gl(2|1) with two nesting levels, momentum rapidities and spin rapidities.

A pivotal structural fact for lattice models is the solution of the quantum inverse problem: in fundamental models, where every lattice site carries the fundamental representation of the R-matrix, every local operator can be written as a polynomial in the monodromy-matrix entries evaluated at shifted spectral parameters (Maillet–Terras 2000; Göhmann–Korepin 2000). This converts statements about local operators into statements about monodromy entries, for which the ABA gives complete control. No analogous construction is known for continuum quantum gases: there is no lattice-site structure, and for multi-component (nested) gases such as the Gaudin–Yang model, generic local operators — the fields, the spin density, spin-flip bilinears, pairing operators, and their products — have not been reconstructed in terms of monodromy entries or Bethe-ansatz data. For the one-component (gl(2)) Lieb–Liniger Bose gas, the corresponding analytic properties of local-operator form factors were demonstrated by taking a scaling limit of the XXZ spin chain (the anisotropic spin-1/2 Heisenberg chain; Golzer–Holz 1987; Seel et al. 2007; Pozsgay 2011); this gl(2) demonstration is the only concrete precedent.

The stakes of this gap were made explicit by Hutsalyuk, Pozsgay and Pristyák (2021), who generalized the LeClair–Mussardo (LM) series to nested Bethe ansatz systems. The LM series is an integral series for one- and two-point functions of local operators in finite-density Bethe states, and underlies exact finite-temperature and generalized-Gibbs-ensemble (statistical ensembles that fix every conserved charge, describing post-quench steady states) correlation functions. Its derivation rests on a finite-volume expansion theorem: the mean value of an operator in an on-shell Bethe state equals a sum over bipartite partitions of the rapidity sets, involving Gaudin matrices (the Jacobians of the Bethe equations) and symmetric form factors (a diagonal limit of off-shell form factors). The proof is an induction on particle numbers and requires residue relations for the operator's form factors: precise pole and residue behavior when a bra-state rapidity approaches a ket-state rapidity, in each nesting level (equations (3.24)–(3.27) of that work). The authors prove these relations for monodromy-matrix entries in gl(2|1) and gl(3) models from multiple-action formulas. For the Gaudin–Yang model they then assume that form factors of physical local operators obey the same relations, stating that in models of quantum gases the solution of the quantum inverse problem is not known and that the assumed analogy is supported only by a scaling-limit argument and the gl(2) demonstration. For the specific density operators of their explicit application they bypass the assumption through an exact generating function built from partial (spatial-subsegment) monodromy matrices and their zero modes (the conserved particle numbers of a segment), but the theorem as stated for generic local operators retains the unproven assumption.

Adjacent partial results exist but do not close the gap. For non-fundamental lattice models, form factors of local operators were computed in a weak sense via composite/partial monodromy matrices and zero modes, and in the supersymmetric t-J chain (a lattice model of correlated electrons with constrained hopping and spin exchange) via universal form factors that explicitly bypass the missing solution of the quantum inverse problem (Pakuliak–Ragoucy–Slavnov; Fuksa–Slavnov). For the continuum two-component Bose gas, determinant form factors are known for the elementary fields and density bilinears, obtained from partial zero modes (Pakuliak–Ragoucy–Slavnov 2015) — the most advanced continuum-nested result — but generic composite local operators are not covered, and no proof of the residue relations or of the expansion theorem exists for that model. An alternative axiomatic route, the thermodynamic form-factor bootstrap program (Cortés Cubero–Panfil 2019; Panfil–Konik 2023), targets form factors of local operators on finite-energy-density states but relies on partly conjectural axioms and its explicit computations concern single-component models.

Problem Statement

Solve the quantum inverse problem for local operators in continuum multi-component quantum gases solvable by the nested Bethe ansatz — the class exemplified by the Gaudin–Yang model (the one-dimensional spin-1/2 delta-interacting Fermi gas, of gl(2|1) symmetry in the algebraic Bethe ansatz) and including the gl(3) two-component Bose gas. The class is the target: a reconstruction must hold for the whole class of nested continuum gases, not for a single member such as the Gaudin–Yang model alone. Concretely: express generic local operators of these models — at a minimum the elementary fields and their conjugates, the density bilinears for all component pairs (including spin-flip and pairing combinations), and their normal-ordered products and derivatives — as operator identities on the physical Hilbert space, each written as a polynomial combination of the total or partial monodromy-matrix entries evaluated at the appropriate shifted spectral parameters, in the way available for fundamental lattice models since Maillet–Terras and Göhmann–Korepin.

The reconstruction immediately yields, and an alternative acceptable resolution directly establishes, the load-bearing application: prove (or disprove) that form factors of such generic local operators between twisted on-shell Bethe states (states satisfying the Bethe equations, taken with twisted boundary conditions that shift the rapidities and lift degeneracies) satisfy the same residue relations as the monodromy-matrix entries — that is, the prescribed pole and residue behavior when a bra-state rapidity approaches a ket-state rapidity, in each nesting level (equations (3.24)–(3.27) of Hutsalyuk–Pozsgay–Pristyák 2021) — thereby justifying — or refuting — the induction step of the finite-volume expansion theorem (Theorem 3.1 of that work) and the LeClair–Mussardo series for nested continuum systems. A valid solution may proceed by a direct operator reconstruction, by a controlled scaling limit from a fundamental lattice model that provably preserves the residue structure through both nesting levels and the graded gl(2|1) structure, or by any other rigorous route; a valid refutation is an explicit local operator of one of these models together with explicit finite-volume Bethe states whose form factors violate the assumed residue relations.

References

  1. A. Hutsalyuk, B. Pozsgay, L. Pristyák, 'The LeClair-Mussardo series and nested Bethe Ansatz', Nuclear Physics B 964 (2021) 115306. DOI: 10.1016/j.nuclphysb.2021.115306; arXiv:2009.13203; https://arxiv.org/abs/2009.13203
  2. J. M. Maillet, V. Terras, 'On the quantum inverse scattering problem', Nuclear Physics B 575 (2000) 627-644. DOI: 10.1016/S0550-3213(00)00097-3; https://doi.org/10.1016/S0550-3213(00)00097-3
  3. F. Göhmann, V. E. Korepin, 'Solution of the quantum inverse problem', Journal of Physics A: Mathematical and General 33 (2000) 1199-1220. DOI: 10.1088/0305-4470/33/6/308; https://doi.org/10.1088/0305-4470/33/6/308
  4. B. Pozsgay, 'Local correlations in the 1D Bose gas from a scaling limit of the XXZ chain', Journal of Statistical Mechanics (2011) P11017. DOI: 10.1088/1742-5468/2011/11/P11017; arXiv:1108.6224; https://arxiv.org/abs/1108.6224
  5. A. Seel, T. Bhattacharyya, F. Göhmann, A. Klümper, 'A note on the spin-1/2 XXZ chain concerning its relation to the Bose gas', Journal of Statistical Mechanics (2007) P08030. DOI: 10.1088/1742-5468/2007/08/P08030; arXiv:0705.3569; https://arxiv.org/abs/0705.3569
  6. B. Golzer, A. Holz, 'The nonlinear Schrödinger model as a special continuum limit of the anisotropic Heisenberg model', Journal of Physics A: Mathematical and General 20 (1987) 3327-3338. DOI: 10.1088/0305-4470/20/11/035; https://doi.org/10.1088/0305-4470/20/11/035
  7. S. Pakuliak, E. Ragoucy, N. A. Slavnov, 'Form factors of local operators in a one-dimensional two-component Bose gas', Journal of Physics A: Mathematical and Theoretical 48 (2015) 435001. DOI: 10.1088/1751-8113/48/43/435001; arXiv:1503.00546; https://arxiv.org/abs/1503.00546
  8. J. Fuksa, N. A. Slavnov, 'Form factors of local operators in supersymmetric quantum integrable models', Journal of Statistical Mechanics (2017) 043106. DOI: 10.1088/1742-5468/aa6686; arXiv:1701.05866; https://arxiv.org/abs/1701.05866
  9. S. Pakuliak, E. Ragoucy, N. A. Slavnov, 'GL(3)-based quantum integrable composite models. II. Form factors of local operators', Symmetry, Integrability and Geometry: Methods and Applications 11 (2015) 064. DOI: 10.3842/SIGMA.2015.064; arXiv:1502.01966; https://arxiv.org/abs/1502.01966
  10. M. Borsi, B. Pozsgay, L. Pristyák, 'Current operators in integrable models: a review', Journal of Statistical Mechanics (2021) 094001. DOI: 10.1088/1742-5468/ac0f6b; arXiv:2103.12160; https://arxiv.org/abs/2103.12160
  11. A. Hutsalyuk, A. Liashyk, 'Master equation for correlation functions in algebra symmetry gl(2|1) related models', arXiv:2102.05017 (2021, unpublished SciPost submission); https://arxiv.org/abs/2102.05017
  12. M. Panfil, R. M. Konik, 'Extending the thermodynamic form factor bootstrap program: multiple particle-hole excitations, crossing symmetry, and reparameterization invariance', Journal of High Energy Physics 2023, article 072. DOI: 10.1007/JHEP10(2023)072; arXiv:2306.10748; https://arxiv.org/abs/2306.10748
  13. A. Cortés Cubero, M. Panfil, 'Thermodynamic bootstrap program for integrable QFT's: form factors and correlation functions at finite energy density', Journal of High Energy Physics 2019, article 104. DOI: 10.1007/JHEP01(2019)104; arXiv:1809.02044; https://arxiv.org/abs/1809.02044
  14. C. N. Yang, 'Some Exact Results for the Many-Body Problem in one Dimension with Repulsive Delta-Function Interaction', Physical Review Letters 19 (1967) 1312-1315. DOI: 10.1103/PhysRevLett.19.1312; https://doi.org/10.1103/PhysRevLett.19.1312