Abstract We propose an effective Bethe ansatz (EBA) for solving quantum many-body systems near an integrable point. Our approach retains the functional form of the Bethe wave function while renormalizing the Bethe roots to account for integrability-breaking interactions. These effective roots are determined by minimizing physically motivated cost functions. The resulting off-shell Bethe states serve as approximate eigenstates of the non-integrable models. We assess the quality of the approximation using various physical observables, including the energy eigenvalue, state fidelity, and bipartite entanglement entropy. Our tests show that for models with weak integrability-breaking, the EBA provides a high-quality approximation to the exact eigenstates over a wide range of deformation parameters. In contrast, for models with strong integrability-breaking interactions, the efficacy of EBA degrades relatively quickly as the deformation parameter increases. The efficacy of the method thus offers a useful probe for characterizing the strength of integrability breaking. Within its regime of accuracy, it also provides a new representation of the eigenstates of nearly integrable models, enabling one to exploit the algebraic structure inherited from integrability.
Abstract We propose an effective Bethe ansatz (EBA) for solving quantum many-body systems near an integrable point. Our approach retains the functional form of the Bethe wave function while renormalizing the Bethe roots to account for integrability-breaking interactions. These effective roots are determined by minimizing physically motivated cost functions. The resulting off-shell Bethe states serve as approximate eigenstates of the non-integrable models. We assess the quality of the approximation using various physical observables, including the energy eigenvalue, state fidelity, and bipartite entanglement entropy. Our tests show that for models with weak integrability-breaking, the EBA provides a high-quality approximation to the exact eigenstates over a wide range of deformation parameters. In contrast, for models with strong integrability-breaking interactions, the efficacy of EBA degrades relatively quickly as the deformation parameter increases. The efficacy of the method thus offers a useful probe for characterizing the strength of integrability breaking. Within its regime of accuracy, it also provides a new representation of the eigenstates of nearly integrable models, enabling one to exploit the algebraic structure inherited from integrability.
Abstract This work presents a comprehensive benchmark and validation of a recently proposed method called Effective Bethe Ansatz (EBA). It is a variational method that deforms the exact Bethe wavefunctions of one-dimensional spin chains at integrable points to approximate non-integrable systems. We apply this method to the non-integrable regime of the spin-1 bilinear–biquadratic chain. By performing EBA method starting from the two integrable endpoints, the Takhtajan-Babujian point and the Lai-Sutherland point, we systematically evaluate the accuracy of the EBA for the ground state and first excited state. Our validation is based on a direct comparison with exact diagonalization, assessing energy, fidelity, and entanglement entropy. The results confirm that the EBA provides a quantitatively accurate description in a finite window around the integrable points, while its fidelity and entanglement properties degrade in a controlled way as the perturbation increases. The method successfully captures key finite-size effects, such as level crossings, manifested as sharp drops in fidelity, and provides a probe to potential phase transitions. This study establishes the EBA as a reliable and efficient semi-analytical tool, clarifying its scope and limitations for studying low-energy physics in non-integrable quantum spin chains.
We investigate the quantum entanglement of Bethe states across a family of integrable spin chains, including the XXX21 model, its higher-spin generalizations (XXXs), and the non-compact SL(2,R) chain. For on-shell eigenstates, we perform a comprehensive scan of the bipartite entanglement entropy across the entire spectrum of finite chains with periodic boundary conditions, and identify the Bethe solutions that minimize and maximize the entanglement. These extremal solutions follow systematic, spin-dependent patterns in the Bethe quantum numbers. In the XXX21 spin chain, for the antiferromagnetic chain, the state with minimal entropy always coincides with the lowest-energy state (the ground state) within a given fixed-magnon sector. For the higher-spin XXXs model, however, the lowest-entropy state is not always identical to the ground state, and can even be the state of highest energy. By contrast, the Bethe roots that maximize entropy exhibit considerably more intricate structure. Our analysis further reveals how special Bethe root configurations, such as singular and strange solutions, affect entanglement, and it uncovers characteristic entanglement features in the non-compact SL(2,R) chain that are absent from compact spin chains. For off-shell Bethe states, we develop an optimization algorithm that extremizes the entanglement entropy over rapidity distributions, enabling us to explore the maximum entanglement achievable by a Bethe state without imposing the Bethe ansatz equations.
The Bethe Ansatz is a powerful method in the theory of quantum integrable models, essential for determining the energy spectrum of dynamical systems - from spin chains in magnetism to models in high-energy physics. This book provides a comprehensive introduction to the Bethe ansatz, from its historical roots to modern developments. First introduced by Hans Bethe in 1931, the method has evolved into a universal framework encompassing algebraic, analytic, thermodynamic, and functional forms. The book explores various Bethe ansatz techniques and their interrelations, covering both coordinate and algebraic versions, with particular attention to nested structures and functional relations involving transfer matrices. Advanced tools such as the separation of variables method are presented in detail. With a wealth of worked examples and precise calculations, this volume serves as an accessible and rigorous reference for graduate students and researchers in mathematical physics and integrable systems.
This thesis presents an introduction to the class of Richardson-Gaudin integrable models, with special focus on the Bethe ansatz wave function, and investigates ways of applying the properties of Richardson-Gaudin models both in and out of integrability. A framework is outlined for the numerical and theoretical treatment of these systems, exposing a duality allowing the Bethe equations to be solved numerically. This is extended to the calculation of inner products and correlation functions. Using this framework, the influence of particle exchange on the Bethe ansatz is discussed, after which it is shown how the Bethe ansatz is able to accurately model wave functions of non-integrable models in two different settings. First, a variational approach is outlined for stationary models where integrability-breaking perturbations are explicitly introduced. Second, an alternative way of breaking integrability is through the introduction of dynamics and periodic driving, where it is shown how integrability can be used to model the resulting Floquet many-body resonances. Throughout this work, it is shown how the clear-cut structure and relatively large freedom in Richardson-Gaudin models makes them ideal for an investigation of the general principles of integrability, as well as being a perfect testing ground for the development of new quantum many-body techniques beyond integrability.
Integrable systems take a special niche in quantum many-body systems. Due to the underlying algebraic structure, it is possible to capture all the quantum correlations exactly within a Bethe Ansatz wavefunction at a favourably scaling computational cost. Richardson [1] and Gaudin [2] have shown that the reduced Bardeen-Cooper-Schrieffer (BCS) Hamiltonian of superconductivity [3] is integrable and that the Bethe Ansatz state corresponds to a product wave function of generalized Cooper pairs [4]. In this formulation, the Bethe ansatz state is equivalent to a non-orthogonal geminal, which is widely used in quantum chemistry to describe electron pair correlations in atoms and molecules [5]. Inspired by Richardson-Gaudin integrable systems, it is possible to construct classes of non-orthogonal geminals, allowing for an accurate description of realistic and non-integrable strongly correlated atoms and molecules, while retaining the favorable computational scaling [6, 7]. In this contribution, I will discuss Richardson-Gaudin integrability and why it is useful for strongly correlated systems.
We present a one-dimensional multicomponent model, known to be partially integrable when restricted to the subspaces made of only two components. By constructing fully antisymmetrized bases, we find integrable excited eigenstates corresponding to the totally antisymmetric irreducible representation of the permutation operator in the otherwise nonintegrable subspaces. We establish rigorously the breakdown of integrability in those subspaces by showing explicitly the violation of the Yang-Baxter equation. We further solve the constraints from the Yang-Baxter equation to find exceptional momenta that allows Bethe ansatz solutions of solitonic bound states. These integrable eigenstates have distinct dynamical consequence from the embedded integrable subspaces previously known, as they do not span their separate Krylov subspaces, and a generic initial state can partly overlap with them and therefore have slow thermalization. However, this novel form of weak ergodicity breaking contrasts with that of quantum many-body scars in that the integrable eigenstates involved do not have necessarily low entanglement. Our approach provides a complementary route to arrive at exact excited states in nonintegrable models: instead of solving towers of single-mode excited states based on a solvable ground state in a nonintegrable model, we identify the integrable eigenstates that survive in a deformation of the Hamiltonian away from its integrable point.