Bethe Ansatz Research Articles

New integrable RG flows with parafermions

Abstract We consider irrelevant deformations of massless RSOS scattering theories by an infinite number of higher [$$ T\overline{T} $$ T T ¯ ]s operators which introduce extra non-trivial CDD factors between left-movers and right-movers. It is shown that the resulting theories can be UV complete after bypassing typical Hagedorn-like singularities if the coefficients of the deformations are fine-tuned. By classifying all integrable cases, we have discovered a new UV complete QFT associated to the $$ {\mathcal{M}}_p $$ M p (p ≥ 3) minimal CFT based on the integrable structure of the RSOS scattering theory. This new theory is the massless ℤp−1 parafermionic sinh-Gordon (PShG) model with a self-dual coupling constant. This correspondence is confirmed by showing that the scale-dependent vacuum energies computed by the thermodynamic Bethe ansatz derived from the S-matrices match those from the quantization conditions for the PShG models using the reflection amplitudes.

Lifted TASEP: A Solvable Paradigm for Speeding up Many-Particle Markov Chains

Virtually all Markov-chain Monte Carlo algorithms used for sampling a given distribution are reversible, and they satisfy the detailed-balance condition. For local chains, this leads to a slow, diffusive exploration of sample space. Significant speedups can be achieved through nonreversible algorithms with the given distribution as a targeted steady state. However, nonreversible algorithms for sampling are difficult to set up and to analyze, and exact speedup results for interacting many-particle systems are very rare. Here, we introduce the “lifted” totally asymmetric simple exclusion process (TASEP) as an exactly solvable paradigm for nonreversible many-particle Markov chains. It samples the same hard-sphere distribution as the Metropolis algorithm for symmetrically diffusing hard-core particles on a one-dimensional lattice. We solve the lifted TASEP by an unusual kind of coordinate Bethe ansatz and show that it exhibits polynomial (in particle number) speedups in the relaxation time for the asymptotic approach of the steady state, as well as the nonasymptotic mixing time, compared to both Metropolis and Kardar-Parisi-Zhang-based dynamics. The lifted TASEP is the reduction onto the one-dimensional lattice of the successful hard-sphere event-chain Monte Carlo algorithm, and we discuss that it can likewise be generalized to soft interaction potentials. Published by the American Physical Society 2024

Bethe vectors and recurrence relations for twisted Yangian based models

We study Olshanski twisted Yangian based models, known as one-dimensional “soliton non-preserving” open spin chains, by means of algebraic Bethe Ansatz. The even case, when the bulk symmetry is \mathfrak{gl}_{2n}𝔤𝔩2n and the boundary symmetry is \mathfrak{sp}_{2n}𝔰𝔭2n or \mathfrak{so}_{2n}𝔰𝔬2n, was studied in [Ann. Henri Poincaré 20, 339 (2018)]. In the present work, we focus on the odd case, when the bulk symmetry is \mathfrak{gl}_{2n+1}𝔤𝔩2n+1 and the boundary symmetry is \mathfrak{so}_{2n+1}𝔰𝔬2n+1. We explicitly construct Bethe vectors and present a more symmetric form of the trace formula. We use the composite model approach and Y(\mathfrak{gl}_n)Y(𝔤𝔩n)-type recurrence relations to obtain recurrence relations for twisted Yangian based Bethe vectors, for both even and odd cases.

Bosonization and correlators of the one-dimensional Hubbard model

Abstract We present a simple derivation of the Luttinger liquid relation for the one-dimensional Hubbard model, considering a finite U and in the U = ∞ limit. We describe a simple solution of the Hubbard model in the infinite repulsion limit and use it to calculate the correlators of the model in this limit in a simple and physical way using the bosonization technique. We then calculate the asymptotics of the correlators of the model at arbitrary U through a single parameter, which can be calculated from the Bethe Ansatz solution. Our derivation of the critical exponents is simple and allows one to express different physical operators of the Hubbard model through the charge and spin Bose fields in a direct and physically transparent way.

Infinite temperature spin dc conductivity of the spin-1/2 XXZ chain

Abstract Using the Bethe ansatz method and the thermodynamic Bethe ansatz equations for the higher-spin integrable XXZ chain, the regular zero frequency contribution to the spin current correlation (spin dc conductivity) is analyzed for the spin-1/2 XXZ chain with the anisotropy 0 ⩽ Δ < 1 . In the high temperature limit, we write down the dressed scattering kernels by one quasi-particle bare energies and exactly evaluate the spin dc conductivity at zero magnetic field, which is proportional to the inverse temperature β at the β → 0 limit. We find that the high temperature proportionality constant σ 0 reg of the spin dc conductivity is discontinuous at all rational numbers of the anisotropy parameter p 0 = π / cos − 1 ⁡ Δ ⩾ 2 with the gap increasing larger than the second power of growing magnetization on one quasi-particle. The isotropic Δ = 1 point is exceptional. Close to this point, σ 0 reg slowly increases in proportion to the first power of the one quasi-particle magnetization. On the other hand, σ 0 reg is proportional to the second power of the same when p 0 approaches irrational numbers.

Superconformal indices of 3d N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 2 SCFTs and holography

We study the superconformal index of 3d N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 2 superconformal field theories on S1 ×ωS2 in the Cardy-like limit where the radius of the S1 is much smaller than that of the S2. We show that the first two leading terms in this Cardy-like expansion are dictated by the Bethe Ansatz formulation of the topologically twisted index of the same theory. We apply this relation to 3d N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 2 holographic superconformal field theories describing the low-energy dynamics of N M2-branes and derive closed form expressions, valid to all orders in the 1/N expansion, for the two leading terms in the Cardy-like expansion of the superconformal index. We also discuss the implications of our results for the entropy of supersymmetric Kerr-Newman black holes in AdS4 and the four-derivative corrections to 4d gauged supergravity.

Higher-rank sectors in the hexagon formalism and marginal deformations

Abstract The hexagon approach provides an integrability framework for the computation of structure constants in N = 4 super Yang-Mills theory in four dimensions. Three-point functions are cut into two hexagonal patches, on which the excitations of the long-range Bethe ansatz of the spectrum problem scatter. To this end, the Bethe states representing the operators also need to be cut into two parts. In rank-one sectors such entangled states are fairly straightforward to construct so that most applications of the method have so far been restricted to this simplest case. In this article we construct entangled states for operators in psu(1, 1|2) sectors, importing a minimum of information from the nested Bethe ansatz. The idea is successfully tested against free field theory for a sample set of correlators with up to three higher-rank operators.
Further, we take a look at the same correlators in the presence of marginal deformations of the theory. While a systematic modification of the hexagon procedure remains out of reach for now, in practical applications the undeformed amplitudes are surprisingly efficient, especially for a certain one-parameter deformation.

Exact results of the one-dimensional repulsive Hubbard model

We present analytical results of the fundamental properties of the one-dimensional (1D) Hubbard model with a repulsive interaction. The new model results with arbitrary external fields include: (I) using the exact solutions of the Bethe ansatz equations of the Hubbard model, we first rigorously calculate the gapless spin and charge excitations, exhibiting exotic features of fractionalized spinons and holons. We then investigate the gapped excitations in terms of the spin string and thek-Λstring bound states at arbitrary driving fields, showing subtle differences in spin magnons and chargeη-pair excitations. (II) For a high-density and high spin magnetization region, i.e. near the quadruple critical point, we further analytically obtain the thermodynamic properties, dimensionless ratios and scaling functions near quantum phase transitions. (III) Importantly, we give the general scaling functions at quantum criticality for arbitrary filling and interaction strength. These can directly apply to other integrable models. (IV) Based on the fractional excitations and the scaling laws, the spin-incoherent Luttinger liquid (SILL) with only the charge propagation mode is elucidated by the asymptotic of the two-point correlation functions with the help of conformal field theory. We also, for the first time, obtain the analytical results of the thermodynamics for the SILL. (V) Finally, to capture deeper insights into the Mott insulator and interaction-driven criticality, we further study the double occupancy and propose its associated contact and contact susceptibilities, through which an adiabatic cooling scheme based upon quantum criticality is proposed. In this scenario, we build up general relations among arbitrary external- and internal-potential-driven quantum phase transitions, providing a comprehensive understanding of quantum criticality. Our methods offer rich perspectives of quantum integrability and offer promising guidance for future experiments with interacting electrons and ultracold atoms, both with and without a lattice.

Approach to stationarity for the KPZ fixed point with boundaries

Current fluctuations for the one-dimensional totally asymmetric exclusion process (TASEP) connected to reservoirs of particles, and their large scale limit to the KPZ fixed point in finite volume, are studied using exact methods. Focusing on the maximal current phase for TASEP, corresponding to infinite boundary slopes for the KPZ height field, we obtain for general initial condition an exact expression for the late time correction to stationarity, involving extreme value statistics of Brownian paths. In the special cases of stationary and narrow wedge initial conditions, a combination of Bethe ansatz and numerical conjectures alternatively provide fully explicit exact expressions.

On Bethe equations of 2d conformal field theory

We study the higher spin algebras of two-dimensional conformal field theory from the perspective of quantum integrability. Starting from Maulik-Okounkov instanton R-matrix and applying the procedure of algebraic Bethe ansatz, we obtain infinite commuting families of Hamiltonians of quantum ILW hierarchy parametrized by the shape of the auxiliary torus. We calculate explicitly the first five of these Hamiltonians. Then, we numerically verify that their joint spectrum can be parametriezed by solutions of Litvinov’s Bethe ansatz equations and we conjecture a general formula for the joint spectrum of all ILW Hamiltonians, based on results of Feigin, Jimbo, Miwa and Mukhin.There are two interesting degeneration limits, the infinitely thick and the infinitely thin auxiliary torus. In one of these limits, the ILW hierarchy degenerates to Yangian or Benjamin-Ono hiearchy and the Bethe equations can be easily solved. The other limit is singular but we can nevertheless extract local Hamiltonians corresponding to quantum KdV or KP hierarchy. Litvinov’s Bethe equations in this local limit provide an alternative to Bethe ansatz equations of Bazhanov, Lukyanov and Zamolodchikov, but are more transparent, work at any rank and are manifestly symmetric under triality symmetry of W\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{W} $$\\end{document}1+∞. Finally, we illustrate analytic properties of the solutions of Bethe equations in minimal models, in particular for Lee-Yang CFT. The analytic structure of Bethe roots is very rich as it captures the representation theory of W\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{W} $$\\end{document}N minimal models.

Peierls Transition in Gross-Neveu Model from Bethe Ansatz

The two-dimensional Gross-Neveu model is anticipated to undergo a crystalline phase transition at high baryon charge densities. This conclusion is drawn from the mean-field approximation, which closely resembles models of Peierls instability. We demonstrate that this transition indeed occurs when both the rank of the symmetry group and the dimension of the particle representation contributing to the baryon density are large (the large N limit). We derive this result through the exact solution of the model, developing the large N limit of the Bethe ansatz. Our analytical construction of the large-N solution of the Bethe ansatz equations aligns perfectly with the periodic (finite-gap) solution of the Korteweg–de Vries (KdV) of the mean-field analysis. Published by the American Physical Society 2024

Solving the Yang-Baxter, tetrahedron and higher simplex equations using Clifford algebras

Bethe Ansatz was discovered in 1932. Half a century later its algebraic structure was unearthed: Yang-Baxter equation was discovered, as well as its multidimensional generalizations [tetrahedron equation and d-simplex equations]. Here we describe a universal method to solve these equations using Clifford algebras. The Yang-Baxter equation (d=2), Zamolodchikov's tetrahedron equation (d=3) and the Bazhanov-Stroganov equation (d=4) are special cases. Our solutions form a linear space. This helps us to include spectral parameters. Potential applications are discussed.

Novel ASEP-inspired solutions of the Yang-Baxter equation

Abstract We explore the algebraic structure of a particular ansatz of the Yang-Baxter equation (YBE), which is inspired by the Bethe Ansatz treatment of the asymmetric simple exclusion process spin-model. Various classes of Hamiltonian density arriving from the two types of R-matrices are found, which also appear as solutions of the constant YBE. We identify the idempotent and nilpotent categories of such constant R-matrices and perform a rank-1 numerical search for the lowest dimension. A summary of the final results reveals general non-Hermitian spin-1/2 chain models.

Generalised hydrodynamics of TT¯\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \ extrm{T}\\overline{\ extrm{T}} $$\\end{document}-deformed integrable quantum field theories

In this paper we evaluate the averages of conserved densities and currents associated to charges of generic spin in (1+1)-dimensional massive integrable Quantum Field Theories perturbed by the irrelevant TT¯\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \ extrm{T}\\overline{\ extrm{T}} $$\\end{document} operator. By making use of the Thermodynamic Bethe Ansatz approach and of the theory of Generalised Hydrodynamics, we study the non-equilibrium steady state averages of conserved densities and currents in a partitioning protocol. We show that, in particular limits, averages can be evaluated exactly in terms of quantities known from the unperturbed theory. In the massless limit we recover known results for the energy and momentum currents and generalise those to any higher spin conserved quantities. We extend some of our results to perturbations of the generalised TT¯\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \ extrm{T}\\overline{\ extrm{T}} $$\\end{document} type. For the massive free fermion theory, we find an analytic expression for the effective inverse temperature after at TT¯\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \ extrm{T}\\overline{\ extrm{T}} $$\\end{document} perturbation in terms of the bare inverse temperature by making use of Lambert’s W function.

Thermodynamic formulation of vacuum energy density in flat spacetime and potential implications for the cosmological constant

We propose a thermodynamical definition of the vacuum energy density ρvac, defined as 〈vac|Tμν|vac〉 = − ρvacgμν, in quantum field theory in flat Minkowski space in D spacetime dimensions, which can be computed in the limit of high temperature, namely in the limit β = 1/T → 0. It takes the form ρvac = const ∙ mD where m is a fundamental mass scale and “const” is a computable constant which can be positive or negative depending on interaction couplings. Due to modular invariance ρvac can also be computed in a different non-thermodynamic channel where one spatial dimension is compactifed on a circle of circumference β and we confirm this modularity for free massive theories for both bosons and fermions for D = 2, 3, 4. We list various properties of ρvac that are generally required, for instance ρvac = 0 for conformal field theories, and others, such as the constraint that ρvac has opposite signs for free bosons verses fermions of the same mass, which is related to constraints from supersymmetry. Using the Thermodynamic Bethe Ansatz we compute ρvac exactly for 2 classes of integrable QFT’s in 2D and interpreting some previously known results. We apply our definition of ρvac to Lattice QCD data with two light quarks (up and down) and one additional massive flavor (the strange quark), and find it is negative, ρvac ≈ − (200 MeV)4. Finally we make some remarks on the Cosmological Constant Problem since ρvac is central to any discussion of it.

Analytic thermodynamic properties of the Lieb-Liniger gas

We present a comprehensive review on the state-of-the-art of the approximate analytic approaches describing the finite-temperature thermodynamic quantities of the Lieb-Liniger model of the one-dimensional (1D) Bose gas with contact repulsive interactions. This paradigmatic model of quantum many-body-theory plays an important role in many areas of physics-thanks to its integrability and possible experimental realization using, e.g., ensembles of ultracold bosonic atoms confined to quasi-1D geometries. The thermodynamics of the uniform Lieb-Liniger gas can be obtained numerically using the exact thermal Bethe ansatz (TBA) method, first derived in 1969 by Yang and Yang. However, the TBA numerical calculations do not allow for the in-depth understanding of the underlying physical mechanisms that govern the thermodynamic behavior of the Lieb-Liniger gas at finite temperature. Our work is then motivated by the insights that emerge naturally from the transparency of closed-form analytic results, which are derived here in six different regimes of the gas and which exhibit an excellent agreement with the TBA numerics. Our findings can be further adopted for characterising the equilibrium properties of inhomogeneous (e.g., harmonically trapped) 1D Bose gases within the local density approximation and for the development of improved hydrodynamic theories, allowing for the calculation of breathing mode frequencies which depend on the underlying thermodynamic equation of state. Our analytic approaches can be applied to other systems including impurities in a quantum bath, liquid helium-4, and ultracold Bose gas mixtures.

Functional Bethe Ansatz for a sinh-Gordon Model with Real q

Recently, Bazhanov and Sergeev have described an Ising-type integrable model which can be identified as a sinh-Gordon-type model with an infinite number of states but with a real parameter q. This model is the subject of Sklyanin’s Functional Bethe Ansatz. We develop in this paper the whole technique of the FBA which includes: (1) Construction of eigenstates of an off-diagonal element of a monodromy matrix. The most important ingredients of these eigenstates are the Clebsh-Gordan coefficients of the corresponding representation. (2) Separately, we discuss the Clebsh-Gordan coefficients, as well as the Wigner’s 6j symbols, in details. The later are rather well known in the theory of 3D indices. Thus, the Sklyanin basis of the quantum separation of variables is constructed. The matrix elements of an eigenstate of the auxiliary transfer matrix in this basis are products of functions satisfying the Baxter equation. Such functions are called usually the Q-operators. We investigate the Baxter equation and Q-operators from two points of view. (3) In the model considered the most convenient Bethe-type variables are the zeros of a Wronskian of two well defined particular solutions of the Baxter equation. This approach works perfectly in the thermodynamic limit. We calculate the distribution of these roots in the thermodynamic limit, and so we reproduce in this way the partition function of the model. (4) The real parameter q, which is the standard quantum group parameter, plays the role of the absolute temperature in the model considered. Expansion with respect to q (tropical expansion) gives an alternative way to establish the structure of the eigenstates. In this way we classify the elementary excitations over the ground state.

Quantum phase transitions in one-dimensional nanostructures: a comparison between DFT and DMRG methodologies.

In the realm of quantum chemistry, the accurate prediction of electronic structure and properties of nanostructures remains a formidable challenge. Density functional theory (DFT) and density matrix renormalization group (DMRG) have emerged as two powerful computational methods for addressing electronic correlation effects in diverse molecular systems. We compare ground-state energies ( ), density profiles ( ), and average entanglement entropies ( ) in metals, insulators and at the transition from metal to insulator, in homogeneous, superlattices, and harmonically confined chains described by the fermionic one-dimensional Hubbard model. While for the homogeneous systems, there is a clear hierarchy between the deviations, , and all the deviations decrease with the chain size; for superlattices and harmonic confinement, the relation among the deviations is less trivial and strongly dependent on the superlattice structure and the confinement strength considered. For the superlattices, in general, increasing the number of impurities in the unit cell represents lower precision in the DFT calculations. For the confined chains, DFT performs better for metallic phases, while the highest deviations appear for the Mott and band-insulator phases. This work provides a comprehensive comparative analysis of these methodologies, shedding light on their respective strengths, limitations, and applications. The DFT calculations were performed using the standard Kohn-Sham scheme within the BALDA approach. It integrated the numerical Bethe-Ansatz (BA) solution of the Hubbard model as the homogeneous density functional within a local-density approximation (LDA) for the exchange-correlation energy. The DMRG algorithms were implemented using the ITensor library, which is based on the matrix product states (MPS) ansatz. The calculations were performed until the energy reaches convergence of at least .

Root patterns and exact surface energy of the spin-1 Heisenberg model with generic open boundaries

We investigate the thermodynamic limit and exact surface energy of the isotropic spin-1 Heisenberg chain with integrable generic open boundary conditions by a novel Bethe ansatz method. We obtain the homogeneous (or two-term) Bethe ansatz like equations for the zero roots of the transfer matrix. Based on the patterns of the zero roots, we analytical calculate the densities of zero roots and the surface energies of the model in all regimes of the boundary parameters.

Scaling limit of the ground state Bethe roots for the inhomogeneous XXZ spin - [formula omitted] chain

It is known that for the Heisenberg XXZ spin - 12 chain in the critical regime, the scaling limit of the vacuum Bethe roots yields an infinite set of numbers that coincide with the energy spectrum of the quantum mechanical 3D anharmonic oscillator. The discovery of this curious relation, among others, gave rise to the approach referred to as the ODE/IQFT (or ODE/IM) correspondence. Here we consider a multiparametric generalization of the Heisenberg spin chain, which is associated with the inhomogeneous six-vertex model. When quasi-periodic boundary conditions are imposed the lattice system may be explored within the Bethe Ansatz technique. We argue that for the critical spin chain, with a properly formulated scaling limit, the scaled Bethe roots for the ground state are described by second order differential equations, which are multi-parametric generalizations of the Schrödinger equation for the anharmonic oscillator.