Explicit finite-size normalization of nested Bethe ansatz eigenfunctions for continuum multicomponent quantum gases

The finite-size eigenfunctions of the one-dimensional multicomponent quantum gases with repulsive delta-function interactions — the two-component Gaudin-Yang model and the Sutherland models for three or more components, in both bosonic and fermionic statistics — are known explicitly as nested Bethe ansatz superpositions parameterized by charge rapidities $k$ and by spin rapidities that, for $\kappa$ components, form $\kappa-1$ nested sets $[\bm{\lambda}]=(\bm{\lambda}^{(1)},\dots,\bm{\lambda}^{(\kappa-1)})$ — a single set $\bm{\lambda}=\bm{\lambda}^{(1)}$ in the two-component case ($\kappa=2$), which is the form written in Appendix A of Pâţu and Aupetit-Diallo (2025), but the normalization constant $C(\bm{k},[\bm{\lambda}]|\varepsilon)$ of these eigenfunctions is not known. The problem is to derive and prove an explicit formula for this constant (equivalently for its squared modulus $|C(\bm{k},[\bm{\lambda}]|\varepsilon)|^{2}$, the squared norm of the unnormalized nested Bethe superpositions) for arbitrary particle number $N$, arbitrary numbers $M_1,\dots,M_{\kappa-1}$ of spin rapidities in the nested sets, and arbitrary number of components $\kappa \ge 2$. By analogy with all proven cases — the Lieb-Liniger boson (Korepin), the lattice $SU(3)$, $gl(n)$, $gl(m|n)$, $U_q(\widehat{gl}_n)$, and $o_{2n+1}$ chains — the answer is expected to be a Gaudin-type determinant of the Jacobian of the nested Bethe equations, but no such formula has been published for the continuum multicomponent gases.

Background

The Lieb-Liniger model describes $N$ bosons on a line interacting via a repulsive delta-function potential of strength $c$; its exact eigenfunctions are superpositions of plane waves with permuted quasimomenta $k_1,\dots,k_N$ (the coordinate Bethe ansatz), and the $k_j$ are fixed by the Bethe equations. When the particles carry an internal spin with $\kappa \ge 2$ components and the interaction is spin-independent, the exact solution requires the nested Bethe ansatz: the wavefunction is still a superposition over charge-sector permutations, but the coefficients are themselves eigenfunctions of an inhomogeneous spin chain, parameterized by spin rapidities that form $\kappa-1$ nested sets $[\bm{\lambda}]=(\bm{\lambda}^{(1)},\dots,\bm{\lambda}^{(\kappa-1)})$, each set $\bm{\lambda}^{(a)}=(\lambda^{(a)}_1,\dots,\lambda^{(a)}_{M_a})$ consisting of $M_a$ rapidities (the single set $\lambda_1,\dots,\lambda_M$ for $\kappa=2$). This yields the Gaudin-Yang model for $\kappa=2$ (two-component fermions, after Yang, or the two-component spinor Bose gas) and the Sutherland models for $\kappa \ge 3$. For a finite system of length $L$ the resulting eigenfunction $\psi(\bm{x};\bm{\sigma}|\bm{k},[\bm{\lambda}])$ — with particle coordinates $\bm{x}$ and spin labels $\bm{\sigma}$ — is explicit: for the spinor gases it is written, e.g., in Appendix A of Pâţu and Aupetit-Diallo (2025) (which spells out the two-component case $\kappa=2$; for $\kappa \ge 3$ the same structural form holds, with the nested spin-chain amplitudes carrying the $\kappa-1$ rapidity sets $[\bm{\lambda}]$) as a sum over charge permutations $P_k$ and coordinate permutations $Q$ of sector functions times a bosonic factor $\varphi(P_k|\varepsilon)=\prod_{m<n}(k_{P_m}-k_{P_n}+ic)$ or the fermionic one — up to a single overall constant $1/C(\bm{k},[\bm{\lambda}]|\varepsilon)$, where $\varepsilon=+1$ ($-1$) for bosons (fermions). Gaudin conjectured in 1971 that the squared norm of a Bethe eigenfunction equals the determinant of the Jacobian of the Bethe equations (the Gaudin determinant, or Gaudin hypothesis); Korepin proved this for the Lieb-Liniger and nonlinear Schrödinger models using the quantum inverse scattering method in 1982, and Dorlas gave an independent proof in 1993. For multicomponent systems the situation is sharply worse. On the lattice, the norm of nested Bethe vectors has been proven in increasing generality: Reshetikhin for $SU(3)$-symmetric chains, Tarasov and Varchenko for $gl(n)$ via quantized Knizhnik-Zamolodchikov equations, Hutsalyuk, Liashyk, Pakuliak, Ragoucy and Slavnov for $gl(m|n)$-symmetric and $U_q(\widehat{gl}_n)$-invariant chains, and Liashyk, Pakuliak and Ragoucy for $o_{2n+1}$-invariant models (2025). For the continuum multicomponent gases, however, only a conjecture exists even for the closely related one-dimensional Hubbard chain (Göhmann and Korepin, 1999), and the algebraic Bethe ansatz program for the continuum two-component Bose gas (Slavnov; Pakuliak-Ragoucy-Slavov) deliberately leaves eigenvector normalization untreated, working with unnormalized states and determinant form factors. The question was re-posed in 2025 by Pâţu and Aupetit-Diallo in their exact study of the free expansion of spinor gases from power-law traps: their stationary-phase asymptotics of the time-dependent wavefunction requires the full nested Bethe eigenfunction, and they state plainly that the normalization constant $C(\bm{k},[\bm{\lambda}]|\varepsilon)$ is unknown — known in the Lieb-Liniger case both for finite and infinite systems (Gaudin; Korepin; Dorlas), but for multicomponent systems conjectured only for the finite-size Hubbard model. The constant cancels in their final asymptotic formulas, so the gap is structural rather than an obstruction to their results: it blocks normalized overlaps, form factors and correlation functions for the continuum multicomponent gases at finite coupling.

Problem Statement

Derive and prove an explicit finite-size normalization formula for the nested Bethe ansatz eigenfunctions of the one-dimensional continuum multicomponent quantum gases with repulsive delta-function interaction of strength $c>0$. Concretely: for the two-component Gaudin-Yang model ($\kappa = 2$, spin rapidities $\bm{\lambda}=\bm{\lambda}^{(1)}$) and the multicomponent Sutherland models ($\kappa \ge 3$ components, spin rapidities organized in the $\kappa-1$ nested sets $[\bm{\lambda}]=(\bm{\lambda}^{(1)},\dots,\bm{\lambda}^{(\kappa-1)})$), for both bosonic ($\varepsilon=+1$) and fermionic ($\varepsilon=-1$) statistics, for arbitrary particle number $N$, arbitrary numbers $M_1,\dots,M_{\kappa-1}$ of spin rapidities in the nested sets, and a finite system of length $L$ with periodic boundary conditions and on-shell rapidities (i.e., $\bm{k}$ and $[\bm{\lambda}]$ solving the nested Bethe equations), determine the normalization constant $C(\bm{k},[\bm{\lambda}]|\varepsilon)$ of the eigenfunction of the form written in Appendix A of Pâţu and Aupetit-Diallo (2025) — equivalently, a closed-form expression for the squared norm $\langle\Phi(\bm{k},[\bm{\lambda}]|\varepsilon)|\Phi(\bm{k},[\bm{\lambda}]|\varepsilon)\rangle=|C(\bm{k},[\bm{\lambda}]|\varepsilon)|^{2}$ of the unnormalized nested Bethe superposition $\Phi$ (the sum without the overall constant $1/C$) as a function of the charge rapidities $\bm{k}$, the spin rapidity sets $[\bm{\lambda}]$, $c$, $N$, $M_1,\dots,M_{\kappa-1}$, $L$ and $\kappa$ — together with a proof that the formula is correct. The formula must (i) reduce to the proven Gaudin-Korepin determinant when the spin sector is trivial ($\kappa=1$, or all $M_a=0$), (ii) respect the separate permutation symmetry of $C$ in the $k_j$ and in the rapidities within each nested set $\bm{\lambda}^{(a)}$, and (iii) be valid for generic (non-degenerate, non-string) Bethe roots, with the degenerate-root limit stated if it differs. A complete answer is expected to confirm or refute the natural expectation, by analogy with all proven lattice cases and the Göhmann-Korepin Hubbard conjecture, that the norm is a Gaudin-type determinant built from the Jacobian of the full nested system of Bethe equations.

References

1. O. I. Pâţu and G. Aupetit-Diallo, "Expansion of one-dimensional spinor gases from power-law traps", Physical Review A 112, 023301 (2025). arXiv:2505.05186, doi:10.1103/clxm-5yd9. URL: https://doi.org/10.1103/clxm-5yd9.
2. E. H. Lieb and W. Liniger, "Exact analysis of an interacting Bose gas. I. The general solution and the ground state", Physical Review 130, 1605-1616 (1963). doi:10.1103/PhysRev.130.1605. URL: https://doi.org/10.1103/PhysRev.130.1605.
3. C. N. Yang, "Some exact results for the many-body problem in one dimension with repulsive delta-function interaction", Physical Review Letters 19, 1312-1315 (1967). doi:10.1103/PhysRevLett.19.1312. URL: https://doi.org/10.1103/PhysRevLett.19.1312.
4. B. Sutherland, "Further results for the many-body problem in one dimension", Physical Review Letters 20, 98-100 (1968). doi:10.1103/PhysRevLett.20.98. URL: https://doi.org/10.1103/PhysRevLett.20.98.
5. B. Sutherland, "Model for a multicomponent quantum system", Physical Review B 12, 3795-3805 (1975). doi:10.1103/PhysRevB.12.3795. URL: https://doi.org/10.1103/PhysRevB.12.3795.
6. M. Gaudin, "Bose Gas in One Dimension. I. The Closure Property of the Scattering Wavefunctions", Journal of Mathematical Physics 12, 1674-1676 (1971). doi:10.1063/1.1665790. URL: https://doi.org/10.1063/1.1665790.
7. M. Gaudin, "Bose Gas in One Dimension. II. Orthogonality of the Scattering States", Journal of Mathematical Physics 12, 1677-1680 (1971). doi:10.1063/1.1665791. URL: https://doi.org/10.1063/1.1665791.
8. V. E. Korepin, "Calculation of norms of Bethe wave functions", Communications in Mathematical Physics 86, 391-418 (1982). doi:10.1007/BF01212176. URL: https://doi.org/10.1007/BF01212176.
9. T. C. Dorlas, "Orthogonality and completeness of the Bethe Ansatz eigenstates of the nonlinear Schroedinger model", Communications in Mathematical Physics 154, 347-376 (1993). doi:10.1007/BF02097001. URL: https://doi.org/10.1007/BF02097001.
10. N. Yu. Reshetikhin, "Calculation of the norm of Bethe vectors in models with SU(3)-symmetry", Journal of Soviet Mathematics 46, 1694-1706 (1989). doi:10.1007/BF01099200. URL: https://doi.org/10.1007/BF01099200.
11. F. Göhmann and V. E. Korepin, "The Hubbard chain: Lieb-Wu equations and norm of the eigenfunctions", Physics Letters A 263, 293-298 (1999). arXiv:cond-mat/9908114, doi:10.1016/S0375-9601(99)00774-4. URL: https://doi.org/10.1016/S0375-9601(99)00774-4.
12. F. H. L. Essler, H. Frahm, F. Göhmann, A. Klümper and V. E. Korepin, "The One-Dimensional Hubbard Model", Cambridge University Press, Cambridge (2005). doi:10.1017/CBO9780511534843. URL: https://doi.org/10.1017/CBO9780511534843.
13. B. Pozsgay, W.-V. van Gerven Oei and M. Kormos, "On form factors in nested Bethe Ansatz systems", Journal of Physics A: Mathematical and Theoretical 45, 465007 (2012). doi:10.1088/1751-8113/45/46/465007. URL: https://doi.org/10.1088/1751-8113/45/46/465007.
14. N. A. Slavnov, "One-dimensional two-component Bose gas and the algebraic Bethe ansatz", Theoretical and Mathematical Physics 183, 800-821 (2015). arXiv:1502.06749, doi:10.1007/s11232-015-0297-8. URL: https://doi.org/10.1007/s11232-015-0297-8.
15. S. Pakuliak, E. Ragoucy and N. A. Slavnov, "Form factors of local operators in a one-dimensional two-component Bose gas", Journal of Physics A: Mathematical and Theoretical 48, 435001 (2015). arXiv:1503.00546, doi:10.1088/1751-8113/48/43/435001. URL: https://doi.org/10.1088/1751-8113/48/43/435001.
16. A. Hutsalyuk, A. Liashyk, S. Z. Pakuliak, E. Ragoucy and N. A. Slavnov, "Norm of Bethe vectors in models with gl(m|n) symmetry", Nuclear Physics B 926, 256-278 (2018). arXiv:1705.09219, doi:10.1016/j.nuclphysb.2017.11.006. URL: https://doi.org/10.1016/j.nuclphysb.2017.11.006.
17. A. Hutsalyuk, A. Liashyk, S. Z. Pakuliak, E. Ragoucy and N. A. Slavnov, "Scalar products and norm of Bethe vectors for integrable models based on U_q(gl_n-hat)", SciPost Physics 4, 006 (2018). arXiv:1711.03867, doi:10.21468/SciPostPhys.4.1.006. URL: https://doi.org/10.21468/SciPostPhys.4.1.006.
18. A. Liashyk, S. Z. Pakuliak and E. Ragoucy, "Scalar products and norm of Bethe vectors in o_{2n+1} invariant integrable models", SciPost Physics 19, 023 (2025). arXiv:2503.01578, doi:10.21468/SciPostPhys.19.1.023. URL: https://doi.org/10.21468/SciPostPhys.19.1.023.
19. C. Rylands, B. Bertini and P. Calabrese, "Integrable quenches in the Hubbard model", Journal of Statistical Mechanics 2022, 103103 (2022). arXiv:2206.07985, doi:10.1088/1742-5468/ac98be. URL: https://doi.org/10.1088/1742-5468/ac98be.
20. A. Hegedus, "Norm of Bethe-wave functions in the continuum limit", Nuclear Physics B 933, 349-383 (2018). arXiv:1805.02897, doi:10.1016/j.nuclphysb.2018.06.017. URL: https://doi.org/10.1016/j.nuclphysb.2018.06.017.