Bethe Ansatz Research Articles

Separation of spherically and translationally covariant finite quantum spaces within the XXX model

A method of separation of subspaces with definite values of the total spin, its z-projection and quasimomentum, within the state space of the XXX model, i.e. the N-th tensor power of a single qubit, is presented. We construct a complete system of orthogonal projectors onto such subspaces by use of the basis of wavelets and chosen measuring operators based on the total spin. We also demonstrate, on examples of N=6 and 8 magnetic spin rings, that the constructed projectors can be further decomposed into the orthogonal sum of density matrices of eigenstates of the XXX model in the corresponding subspaces. This decomposition exhibits a Galois symmetry, stemming from roots of indecomposable factors of characteristic polynomial of the Heisenberg Hamiltonian. In this way, we propose an alternative for Bethe Ansatz for moderate lengths of spin chains, by a return to the original spectral problem, avoiding the highly nonlinear system of Bethe Ansatz equations. Such results are easy to adapt in quantum information processing.

$q$-opers, $QQ$-systems, and Bethe Ansatz

We introduce the notions of $(G,q)$-opers and Miura $(G,q)$-opers, where $G$ is a simply connected simple complex Lie group, and prove some general results about their structure. We then establish a one-to-one correspondence between the set of $(G,q)$-opers of a certain kind and the set of nondegenerate solutions of a system of Bethe Ansatz equations. This may be viewed as a $q$DE/IM correspondence between the spectra of a quantum integrable model (IM) and classical geometric objects ($q$-differential equations). If $\mathfrak{g}$ is simply laced, the Bethe Ansatz equations we obtain coincide with the equations that appear in the quantum integrable model of XXZ-type associated to the quantum affine algebra $U\_q \widehat{\mathfrak{g}}$. However, if $\mathfrak{g}$ is non-simply-laced, then these equations correspond to a different integrable model, associated to $U\_q {}^L\widehat{\mathfrak{g}}$ where $^L\widehat {\mathfrak{g}}$ is the Langlands dual (twisted) affine algebra. A key element in this $q$DE/IM correspondence is the $QQ$-system that has appeared previously in the study of the ODE/IM correspondence and the Grothendieck ring of the category $\mathcal{O}$ of the relevant quantum affine algebra.

Asymptotics in an asymptotic CFT

In this work we illustrate the resurgent structure of the λ-deformation; a two-dimensional integrable quantum field theory that has an RG flow with an SU(N)k Wess-Zumino-Witten conformal fixed point in the UV. To do so we use modern matched asymptotic techniques applied to the thermodynamic Bethe ansatz formulation to compute the free energy to 38 perturbative orders in an expansion of large applied chemical potential. We find numerical evidence for factorial asymptotic behaviour with both alternating and non-alternating character which we match to an analytic expression. A curiosity of the system is that the leading non-alternating factorial growth vanishing when k divides N. The ambiguities associated to Borel resummation of this series are suggestive of non-perturbative contributions. This is verified with an analytic study of the TBA system demonstrating a cancellation between perturbative and non-perturbative ambiguities.

The Bethe ansatz for sticky Brownian motions

We consider a multi-dimensional diffusion whose coordinates behave as one-dimensional Brownian motions, evolving independently when apart, but with a sticky interaction when they coincide. We derive the Kolmogorov backwards equation and show that for a specific choice of interaction it can be solved exactly with the Bethe ansatz. The diffusion in Rn can be viewed as the n-point motions of a stochastic flow of kernels. We use our formulae to study the flow of kernels and show that atoms in the flow are asymptotically exponentially distributed in size at large times.

Elliptic BCS-Richardson model and the modified algebraic Bethe ansatz

We consider the elliptic Gaudin-type model in an external magnetic field (Skrypnyk T 2005 Phys. Lett. A 334 390–9; Skrypnyk T 2005 Phys. Lett. A 347 266–7; Skrypnyk T 2006 J. Geom. Phys. 57 53–67; Skrypnyk T 2006 J. Math. Phys. 47; Skrypnyk T 2007 J. Phys. A 40 1611–23; Skrypnyk T 2019 Nucl. Phys. B 941 225–48) associated with non-skew-symmetric elliptic r-matrix (Skrypnyk T 2005 Phys. Lett. A 334 390–9; Skrypnyk T 2005 Phys. Lett. A 347 266–7; Skrypnyk T 2006 J. Geom. Phys. 57 53–67; Skrypnyk T 2006 J. Math. Phys. 47). Using them we construct a new integrable fermion Hamiltonian of the Richardson type. We use the modified algebraic Bethe ansatz obtained for integrable models with the considered elliptic r-matrix in (Skrypnyk T 2023 Nucl. Phys. B 988 116102) and find the spectrum of the obtained Richardson-type Hamiltonian in terms of solutions of the modified Bethe equations. The obtained results generalize our previous results on Richardson-type models (Skrypnyk T 2022 Nucl. Phys. B 975 115679).

Manifestation of spin nematic ordering in the spin-1 chain system

Spin-1 chain material is considered. Recently, it was predicted that the spin nematic ordered phase can be realized due to the spin-elastic coupling. In this study we show that the appearance of the spin nematic ordering and phase transitions to such ordered states can be studied using magneto-acoustic experiments. One can observe such phases and phase transitions considering relative changes of the velocity of sound (related to the changes of the elastic modulus of the crystal) as a function of temperature and the external magnetic field. It is shown that those relative changes are proportional to the spin quadrupole susceptibility of the system. The results are obtained using the exact Bethe ansatz solution and Quantum Monte Carlo simulations for several typical values of the biquadratic spin-spin exchange interaction and for various values of the spin-elastic coupling and elastic modules.

Critical site percolation on the triangular lattice: from integrability to conformal partition functions

Critical site percolation on the triangular lattice is described by the Yang–Baxter solvable dilute loop model with crossing parameter specialized to , corresponding to the contractible loop fugacity . We study the functional relations satisfied by the commuting transfer matrices of this model and the associated Bethe ansatz equations. The single and double row transfer matrices are respectively endowed with strip and periodic boundary conditions, and are elements of the ordinary and periodic dilute Temperley–Lieb algebras. The standard modules for these algebras are labelled by the number of defects d and, in the latter case, also by the twist . Nonlinear integral equation techniques are used to analytically solve the Bethe ansatz functional equations in the scaling limit for the central charge c = 0 and conformal weights . For the ground states, we find for strip boundary conditions and for periodic boundary conditions, where . We give explicit conjectures for the scaling limit of the trace of the transfer matrix in each standard module. For , these conjectures are supported by numerical solutions of the logarithmic form of the Bethe ansatz equations for the leading 20 or more conformal eigenenergies. With these conjectures, we apply the Markov traces to obtain the conformal partition functions on the cylinder and torus. These precisely coincide with our previous results for critical bond percolation on the square lattice, described by the dense loop model with . The concurrence of all this conformal data provides compelling evidence supporting a strong form of universality between these two stochastic models as logarithmic conformal field theories.

Correlation properties of a one-dimensional repulsive Bose gas at finite temperature

We present a comprehensive study shedding light on how thermal fluctuations affect correlations in a Bose gas with contact repulsive interactions in one spatial dimension. The pair correlation function, the static structure factor, and the one-body density matrix are calculated as a function of the interaction strength and temperature with the exact ab-initio Path Integral Monte Carlo method. We explore all possible gas regimes from weak to strong interactions and from low to high temperatures. We provide a detailed comparison with a number of theories, such as perturbative (Bogoliubov and decoherent classical), effective (Luttinger liquid) and exact (ground-state and thermal Bethe Ansatz) ones. Our Monte Carlo results exhibit an excellent agreement with the tractable limits and provide a fundamental benchmark for future observations which can be achieved in atomic gases, cavity quantum-electrodynamic and superconducting-circuit platforms.

Exact solution of the q-deformed D3(1) vertex model with open boundaries

In this paper, we study the exact solution of the $q$-deformed $D^{(1)}_3$ quantum lattice model with non-diagonal open boundary condition. We demonstrate the crossing symmetry of the transfer matrix and obtain the quantum determinant. We construct the independent transfer matrix fusion identities and show that the fusion processes can be closed. Based on the fusion hierarchies and polynomial analysis, we obtain the inhomogeneous $T-Q$ relations, exact energy spectrum and Bethe ansatz equations of the system.

Dynamics of quantum double dark-solitons and an exact finite-size scaling of Bose–Einstein condensation

We show several novel aspects in the exact non-equilibrium dynamics of quantum double dark-soliton states in the Lieb–Liniger model for the one-dimensional Bose gas with repulsive interactions. We also show an exact finite-size scaling of the fraction of the quasi-Bose–Einstein condensation (BEC) in the ground state, which should characterize the quasi-BEC in quantum double dark-soliton states that we assume to occur in the weak coupling regime. First, we show the exact time evolution of the density profile in the quantum state associated with a quantum double dark-soliton by the Bethe ansatz. Secondly, we derive a kind of macroscopic quantum wave-function effectively by exactly evaluating the square amplitude and phase profiles of the matrix element of the field operator between the quantum double dark-soliton states. The profiles are close to those of dark-solitons particularly in the weak-coupling regime. Then, the scattering of two notches in the quantum double dark-soliton state is exactly demonstrated. It is suggested from the above observations that the quasi-BEC should play a significant role in the dynamics of quantum double dark-soliton states. If the condensate fraction is close to 1, the quantum state should be well approximated by the quasi-BEC state where the mean-field picture is valid.

Integrability of two-species partially asymmetric exclusion processes

We work towards the classification of all one-dimensional exclusion processes with two species of particles that can be solved by a nested coordinate Bethe ansatz (BA). Using the Yang–Baxter equations, we obtain conditions on the model parameters that ensure that the underlying system is integrable. Three classes of integrable models are thus found. Of these, two classes are well known in literature, but the third has not been studied until recently, and never in the context of the BA. The Bethe equations are derived for the latter model as well as for the associated dynamics encoding the large deviation of the currents.

The cumulant Green’s functions method for the Hubbard model

We use the cumulant Green’s functions method (CGFM) to study the single-band Hubbard model. The starting point of the method is to diagonalize a cluster (‘seed’) containing N correlated sites and employ the cumulants calculated from the cluster solution to obtain the full Green’s functions for the lattice. All calculations are done directly; no variational or self-consistent process is needed. We benchmark the one-dimensional results for the gap, the double occupancy, and the ground-state energy as functions of the electronic correlation at half-filling and the occupation numbers as functions of the chemical potential obtained from the CGFM against the corresponding results of the thermodynamic Bethe ansatz and the quantum transfer matrix methods. The particle-hole symmetry of the density of states is fulfilled, and the gap, occupation numbers, and ground-state energy tend systematically to the known results as the cluster size increases. We include a straightforward application of the CGFM to simulate the singles occupation of an optical lattice experiment with lithium-6 atoms in an eight-site Fermi-Hubbard chain near half-filling. The method can be applied to any parameter space for one, two, or three-dimensional Hubbard Hamiltonians and extended to other strongly correlated models, like the Anderson Hamiltonian, the t − J, Kondo, and Coqblin-Schrieffer models.

All-Orders Quadratic-Logarithmic Behavior for Amplitudes.

We classify origin limits of maximally helicity violating multigluon scattering amplitudes in planar N=4 super-Yang-Mills theory, where a large number of cross ratios approach zero, with the help of cluster algebras. By analyzing existing perturbative data and bootstrapping new data, we provide evidence that the amplitudes become the exponential of a quadratic polynomial in the large logarithms. With additional input from the thermodynamic Bethe ansatz at strong coupling, we conjecture exact expressions for amplitudes with up to eight gluons in all origin limits. Our expressions are governed by the tilted cusp anomalous dimension evaluated at various values of the tilt angle.

Rational $Q$-systems, Higgsing and mirror symmetry

The rational QQ-system is an efficient method to solve Bethe ansatz equations for quantum integrable spin chains. We construct the rational QQ-systems for generic Bethe ansatz equations described by an A_{\ell-1}Aℓ−1 quiver, which include models with multiple momentum carrying nodes, generic inhomogeneities, generic diagonal twists and qq-deformation. The rational QQ-system thus constructed is specified by two partitions. Under Bethe/Gauge correspondence, the rational QQ-system is in a one-to-one correspondence with a 3d \mathcal{N}=4𝒩=4 quiver gauge theory of the type {T}_{{\rho}}^{{\sigma}}[\mathrm{SU}(n)]Tρσ[SU(n)], which is also specified by the same partitions. This shows that the rational QQ-system is a natural language for the Bethe/Gauge correspondence, because known features of the {T}_{{\rho}}^{{\sigma}}[\mathrm{SU}(n)]Tρσ[SU(n)] theories readily translate. For instance, we show that the Higgs and Coulomb branch Higgsing correspond to modifying one of the partitions in the rational QQ-system while keeping the other untouched. Similarly, mirror symmetry is realized in terms of the rational QQ-system by simply swapping the two partitions - exactly as for {T}_{{\rho}}^{{\sigma}}[\mathrm{SU}(n)]Tρσ[SU(n)]. We exemplify the computational efficiency of the rational QQ-system by evaluating topologically twisted indices for 3d \mathcal{N}=4𝒩=4\mathrm{U}(n)U(n) SQCD theories with n=1,\ldots,5n=1,…,5.

Exact quantization and analytic continuation

In this paper we give a streamlined derivation of the exact quantization condition (EQC) on the quantum periods of the Schrödinger problem in one dimension with a general polynomial potential, based on Wronskian relations. We further generalize the EQC to potentials with a regular singularity, describing spherical symmetric quantum mechanical systems in a given angular momentum sector. We show that the thermodynamic Bethe ansatz (TBA) equations that govern the quantum periods undergo nontrivial monodromies as the angular momentum is analytically continued between integer values in the complex plane. The TBA equations together with the EQC are checked numerically against Hamiltonian truncation at real angular momenta and couplings, and are used to explore the analytic continuation of the spectrum on the complex angular momentum plane in examples.

Open-system spin transport and operator weight dissipation in spin chains

We use non-equilibrium steady states to study the effect of dissipation-assisted operator evolution (DAOE) on the scaling behavior of transport in one-dimensional spin chains. We consider three models in the XXZ family: the XXZ model with staggered anisotropy, which is chaotic; XXZ model with no external field and tunable interaction, which is Bethe ansatz integrable and (in the zero interaction limit) free fermion integrable; and the disordered XY model, which is free-fermion integrable and Anderson localized. We find evidence that DAOE's effect on transport is controlled by its effect on the system's conserved quantities. To the extent that DAOE preserves those symmetries, it preserves the scaling of the system's transport properties; to the extent it breaks those conserved quantities, it pushes the system towards diffusive scaling of transport.

Effect of Transverse Confinement on a Quasi-One-Dimensional Dipolar Bose Gas

We study a gas of bosonic dipolar atoms in the presence of a transverse harmonic trapping potential by using an improved variational Bethe ansatz, which includes the transverse width of the atomic cloud as a variational parameter. Our calculations show that the system behavior evolves from quasi-one dimensional to a strictly one-dimensional one by changing the atom–atom interaction, or the axial density, or the frequency of the transverse confinement. Quite remarkably, in the droplet phase induced by the attractive dipolar interaction the system becomes sub-one dimensional when the transverse width is smaller than the characteristic length of the transverse harmonic confinement.

Exact densities of loops in O(1) dense loop model and of clusters in critical percolation on a cylinder: II. Rotated lattice

This work continues the study started in Povolotsky (2021 J. Phys. A: Math. Theor. 54 22LT01), where the exact densities of loops in the O(1) dense loop model on an infinite strip of the square lattice with periodic boundary conditions were obtained. These densities are also equal to the densities of critical percolation clusters on the 45∘ rotated square lattice rolled into a cylinder. Here, we extend those results to the square lattice with a tilt. This in particular allows us to obtain the densities of critical percolation clusters on the cylinder of the square lattice of standard orientation extensively studied before. We obtain exact densities of contractible and non-contractible loops or equivalently the densities of critical percolation clusters, which do not and do wrap around the cylinder, respectively. The solution uses the mapping of O(1) dense loop model to the six-vertex model in the Razumov–Stroganov point, while the effective tilt is introduced via the inhomogeneous transfer matrix proposed by Fujimoto. The further solution is based on the Bethe ansatz and Fridkin–Stroganov–Zagier’s solution of Baxter’s T–Q equation. The results are represented in terms of the solution of two explicit systems of linear algebraic equations, which can be performed either analytically for small circumferences of the cylinder or numerically for larger ones. We present exact rational values of the densities on the cylinders of small circumferences and several lattice orientations and use the results of high precision numerical calculations to study the finite-size corrections to the densities, in particular their dependence on the tilt of the lattice.

An exact solution of the homogenous trimer Bose-Hubbard model

It is shown that the homogenous three-site (trimer) Bose–Hubbard model with a periodic boundary condition can be solved exactly by using an extended Bethe ansatz based on the solution of the dimer model and the site-permutation symmetry. A solution of the model within the S 3 symmetric or non-symmetric subspace is presented. Coupled differential equations of a series of extended Heine–Stieltjes polynomials, the zeros of which are related to the solution of the model, are derived. Numerical examples of the solution of the model with bosons, including the related Heine–Stieltjes polynomials and the Van Vleck polynomials, are presented, which serve to validate the procedure and illustrate the completeness of the solutions they render.

Power-law decay of correlations after a global quench in the massive XXZ chain

We investigate the relaxation dynamics of equal-time correlations in the antiferromagnetic phase of the XXZ spin-1/2 chain following a global quantum quench of the anisotropy parameter. We focus, in particular, on the relaxation dynamics starting from an initial N\'eel state. Using state-of-the-art density-matrix renormalization group simulations, the exact solution of an effective free-fermion model, and the quench-action approach within the thermodynamic Bethe ansatz, we show that the late-time relaxation is characterized by a power-law decay $\sim t^{-3/2}$ independent of anisotropy. This is in contrast to the previously studied exponential decay of the antiferromagnetic order parameter. Remarkably, the effective model describes the numerical data extremely well even on a quantitative level if higher order corrections to the leading asymptotic behavior are taken into account.