Integrable Many-body Systems Research Articles

Commutative families in DIM algebra, integrable many-body systems and q, t matrix models

We extend our consideration of commutative subalgebras (rays) in different representations of the W1+∞ algebra to the elliptic Hall algebra (or, equivalently, to the Ding-Iohara-Miki (DIM) algebra Uq,tgl̂̂1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {U}_{q,t}\\left({\\hat{\\hat{\\mathfrak{gl}}}}_1\\right) $$\\end{document}). Its advantage is that it possesses the Miki automorphism, which makes all commutative rays equivalent. Integrable systems associated with these rays become finite-difference and, apart from the trigonometric Ruijsenaars system not too much familiar. We concentrate on the simplest many-body and Fock representations, and derive explicit formulas for all generators of the elliptic Hall algebra en,m. In the one-body representation, they differ just by normalization from znqmD̂\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {z}^n{q}^{m\\hat{D}} $$\\end{document} of the W1+∞ Lie algebra, and, in the N -body case, they are non-trivially generalized to monomials of the Cherednik operators with action restricted to symmetric polynomials. In the Fock representation, the resulting operators are expressed through auxiliary polynomials of n variables, which define weights in the residues formulas. We also discuss q, t-deformation of matrix models associated with constructed commutative subalgebras.

Particle scattering and fusion for the Ablowitz–Ladik chain

The Ablowitz–Ladik (AL) chain is an integrable discretized version of the nonlinear Schrödinger equation. We report on a novel underlying Hamiltonian particle system with properties similar to the ones known for the classical Toda chain and Calogero fluid with 1/sinh2 pair interaction. Boundary conditions are imposed such that, both in the distant past and future, particles have a constant velocity. We establish the many-particle scattering for the AL chain and obtain properties known for generic integrable many-body systems. For a specific choice of the chain, real initial data remain real in the course of time. Then, asymptotically, particles move in pairs with a velocity-dependent size and scattering shifts are governed by the fusion rule.

New Classical Integrable Systems from Generalized TT[over ¯]-Deformations.

We introduce and study a novel class of classical integrable many-body systems obtained by generalized TT[over ¯] deformations of free particles. Deformation terms are bilinears in densities and currents for the continuum of charges counting asymptotic particles of different momenta. In these models, which we dub "semiclassical Bethe systems" for their link with the dynamics of Bethe ansatz wave packets, many-body scattering processes are factorized, and two-body scattering shifts can be set to an almost arbitrary function of momenta. The dynamics is local but inherently different from that of known classical integrable systems. At short scales, the geometry of the deformation is dynamically resolved: either particles are slowed down (more space available), or accelerated via a novel classical particle-pair creation and annihilation process (less space available). The thermodynamics both at finite and infinite volumes is described by the equations of (or akin to) the thermodynamic Bethe ansatz, and at large scales generalized hydrodynamics emerge.

Three-stage thermalization of a quasi-integrable system

We consider a system of classical hard rods or billiard balls in one dimension, initially prepared in a Bragg-pulse state at a given temperature and subjected to external periodic fields. We show that at late times the system thermalizes in the thermodynamic limit via a three-stages process characterized by: an early phase where the dynamics is well described by Euler hydrodynamics, a subsequent phase where a (weak) turbulent phase is observed and where hydrodynamic gradient expansion can be broken, and a final one where the gas thermalizes according to viscous hydrodynamics. As the hard rod gas shares the same large-scale hydrodynamics as other quantum and classical integrable systems, we expect these features to universally characterize all many-body integrable systems in generic external potentials. Published by the American Physical Society 2024

Quantum quenches in driven-dissipative quadratic fermionic systems with parity-time symmetry

We study the quench dynamics of noninteracting fermionic quantum many-body systems that are subjected to Markovian drive and dissipation and are described by a quadratic Liouvillian which has parity-time (PT) symmetry. In recent work, we have shown that such systems relax locally to a maximum entropy ensemble that we have dubbed the PT-symmetric generalized Gibbs ensemble (PTGGE), in analogy to the generalized Gibbs ensemble that describes the steady state of isolated integrable quantum many-body systems after a quench. Here, using driven-dissipative versions of the Su-Schrieffer-Heeger (SSH) model and the Kitaev chain as paradigmatic model systems, we corroborate and substantially expand upon our previous results. In particular, we confirm the validity of a dissipative quasiparticle picture at finite dissipation by demonstrating light cone spreading of correlations and the linear growth and saturation to the PTGGE prediction of the quasiparticle-pair contribution to the subsystem entropy in the PT-symmetric phase. Further, we introduce the concept of directional pumping phases, which is related to the non-Hermitian topology of the Liouvillian and based upon qualitatively different dynamics of the dual string order parameter and the subsystem fermion parity in the SSH model and the Kitaev chain, respectively: depending on the postquench parameters, there can be pumping of string order and fermion parity through both ends of a subsystem corresponding to a finite segment of the one-dimensional lattice, through only one end, or there can be no pumping at all. We show that transitions between dynamical pumping phases give rise to a new and independent type of dynamical critical behavior of the rates of directional pumping, which are determined by the soft modes of the PTGGE. Published by the American Physical Society 2024

Commutative families in W∞, integrable many-body systems and hypergeometric τ-functions

We explain that the set of new integrable systems, generalizing the Calogero family and implied by the study of WLZZ models, which was described in arXiv:2303.05273, is only the tip of the iceberg. We provide its wide generalization and explain that it is related to commutative subalgebras (Hamiltonians) of the W1+∞ algebra. We construct many such subalgebras and explain how they look in various representations. We start from the even simpler w∞ contraction, then proceed to the one-body representation in terms of differential operators on a circle, further generalizing to matrices and in their eigenvalues, in finally to the bosonic representation in terms of time-variables. Moreover, we explain that some of the subalgebras survive the β-deformation, an intermediate step from W1+∞ to the affine Yangian. The very explicit formulas for the corresponding Hamiltonians in these cases are provided. Integrable many-body systems generalizing the rational Calogero model arise in the representation in terms of eigenvalues. Each element of W1+∞ algebra gives rise to KP/Toda τ-functions. The hidden symmetry given by the families of commuting Hamiltonians is in charge of the special, (skew) hypergeometric τ-functions among these.

Dynamical transition from localized to uniform scrambling in locally hyperbolic systems.

Fast scrambling of quantum correlations, reflected by the exponential growth of out-of-time-order correlators (OTOCs) on short pre-Ehrenfest time scales, is commonly considered as a major quantum signature of unstable dynamics in quantum systems with a classical limit. In two recent works [Phys.Rev. Lett. 123, 160401 (2019)0031-900710.1103/PhysRevLett.123.160401] and [Phys. Rev.Lett. 124, 140602 (2020)10.1103/PhysRevLett.124.140602], a significant difference in the scrambling rate of integrable (many-body) systems was observed, depending on the initial state being semiclassically localized around unstable fixed points or fully delocalized (infinite temperature). Specifically, the quantum Lyapunov exponent λ_{q} quantifying the OTOC growth is given, respectively, by λ_{q}=2λ_{s} or λ_{q}=λ_{s} in terms of the stability exponent λ_{s} of the hyperbolic fixed point. Here we show that a wave packet, initially localized around this fixed point, features a distinct dynamical transition between these two regions. We present an analytical semiclassical approach providing a physical picture of this phenomenon, and support our findings by extensive numerical simulations in the whole parameter range of locally unstable dynamics of a Bose-Hubbard dimer. Our results suggest that the existence of this crossover is a hallmark of unstable separatrix dynamics in integrable systems, thus opening the possibility to distinguish the latter, on the basis of this particular observable, from genuine chaotic dynamics generally featuring uniform exponential growth of the OTOC.

Integrable crosscap states: from spin chains to 1D Bose gas

The notion of a crosscap state, a special conformal boundary state first defined in 2d CFT, was recently generalized to 2d massive integrable quantum field theories and integrable spin chains. It has been shown that the crosscap states preserve integrability. In this work, we first generalize this notion to the Lieb-Liniger model, which is a prototype of integrable non-relativistic many-body systems. We then show that the defined crosscap state preserves integrability. We derive the exact overlap formula of the crosscap state and the on-shell Bethe states. As a byproduct, we prove the conjectured overlap formula for integrable spin chains rigorously by coordinate Bethe ansatz. It turns out that the overlap formula for both models take the same form as a ratio of Gaudin-like determinants with a trivial prefactor. Finally we study quench dynamics of the crosscap state, which turns out to be surprisingly simple. The stationary density distribution is simply a constant. We also derive the analytic formula for dynamical correlation functions in the Tonks-Girardeau limit.

Many-body integrable systems implied by WLZZ models

We provide some details about the recently discovered integrable systems implied by commutativity of W operators along the rays on the plane of roots of w∞-algebra. The simplest system of this type is the rational Calogero model, other systems escaped attention in the past. Existence of these systems is intimately tied to the very interesting WLZZ matrix models, which are now under intensive study.

Periodic orbit theory of Bethe-integrable quantum systems: an N-particle Berry–Tabor trace formula

One of the fundamental results of semiclassical theory is the existence of trace formulae showing how spectra of quantum mechanical systems emerge from massive interference among amplitudes related with time-periodic structures of the corresponding classical limit. If it displays the properties of Hamiltonian integrability, this connection is given by the celebrated Berry–Tabor trace formula, and the periodic structures it is built on are tori supporting closed trajectories in phase space. Here we show how to extend this connection into the domain of quantum many-body systems displaying integrability in the sense of the Bethe ansatz, where a classical limit cannot be rigorously defined due to the presence of singular potentials. Formally following the original derivation of Berry and Tabor (1976 Proc. R. Soc. A 349 101), but applied to the Bethe equations without underlying classical structure, we obtain a many-particle trace formula for the density of states of N interacting bosons on a ring, the Lieb–Liniger model. Our semiclassical expressions are in excellent agreement with quantum mechanical results for and 4 particles. For N = 2 we relate our results to the quantization of billiards with mixed boundary conditions. Our work paves the way towards the treatment of the important class of integrable many-body systems by means of semiclassical trace formulae pioneered by Michael Berry in the single-particle context.

Relaxation to a Parity-Time Symmetric Generalized Gibbs Ensemble after a Quantum Quench in a Driven-Dissipative Kitaev Chain.

The construction of the generalized Gibbs ensemble, to which isolated integrable quantum many-body systems relax after a quantum quench, is based upon the principle of maximum entropy. In contrast, there are no universal and model-independent laws that govern the relaxation dynamics and stationary states of open quantum systems, which are subjected to Markovian drive and dissipation. Yet, as we show, relaxation of driven-dissipative systems after a quantum quench can, in fact, be determined by a maximum entropy ensemble, if the Liouvillian that generates the dynamics of the system has parity-time symmetry. Focusing on the specific example of a driven-dissipative Kitaev chain, we show that, similar to isolated integrable systems, the approach to a parity-time symmetric generalized Gibbs ensemble becomes manifest in the relaxation of local observables and the dynamics of subsystem entropies. In contrast, the directional pumping of fermion parity, which is induced by nontrivial non-Hermitian topology of the Kitaev chain, represents a phenomenon that is unique to relaxation dynamics in driven-dissipative systems. Upon increasing the strength of dissipation, parity-time symmetry is broken at a finite critical value, which thus constitutes a sharp dynamical transition that delimits the applicability of the principle of maximum entropy. We show that these results, which we obtain for the specific example of the Kitaev chain, apply to broad classes of noninteracting fermionic models, and we discuss their generalization to a noninteracting bosonic model and an interacting spin chain.

One-Dimensional Quantum Systems with Ground State of Jastrow Form Are Integrable.

Exchange operator formalism describes many-body integrable systems using phase-space variables involving an exchange operator that acts on any pair of particles. We establish an equivalence between models described by exchange operator formalism and the complete infinite family of parent Hamiltonians describing quantum many-body models with ground states of Jastrow form. This makes it possible to identify the invariants of motion for any model in the family and establish its integrability, even in the presence of an external potential. Using this construction we establish the integrability of the long-range Lieb-Liniger model, describing bosons in a harmonic trap and subject to contact and Coulomb interactions in one dimension. We further identify a variety of models exemplifying the integrability of Hamiltonians in this family.

Emergent tracer dynamics in constrained quantum systems

We show how the tracer motion of tagged, distinguishable particles can effectively describe transport in various homogeneous quantum many-body systems with constraints. We consider systems of spinful particles on a one-dimensional lattice subjected to constrained spin interactions, such that some or even all multipole moments of the effective spin pattern formed by the particles are conserved. On the one hand, when all moments - and thus the entire spin pattern - are conserved, dynamical spin correlations reduce to tracer motion identically, generically yielding a subdiffusive dynamical exponent $z=4$. This provides a common framework to understand the dynamics of several constrained lattice models, including models with XNOR or $tJ_z$ - constraints. We consider random unitary circuit dynamics with such a conserved spin pattern and use the tracer picture to obtain exact expressions for their late-time dynamical correlations. Our results can also be extended to integrable quantum many-body systems that feature a conserved spin pattern but whose dynamics is insensitive to the pattern, which includes for example the folded XXZ spin chain. On the other hand, when only a finite number of moments of the pattern are conserved, the dynamics is described by a convolution of the internal hydrodynamics of the spin pattern with a tracer distribution function. As a consequence, we find that the tracer universality is robust in generic systems if at least the quadrupole moment of the pattern remains conserved. In cases where only total magnetization and dipole moment of the pattern are constant, we uncover an intriguing coexistence of two processes with equal dynamical exponent but different scaling functions, which we relate to phase coexistence at a first order transition.

Integrable quantum many-body sensors for AC field sensing

Quantum sensing is inevitably an elegant example of the supremacy of quantum technologies over their classical counterparts. One of the desired endeavors of quantum metrology is AC field sensing. Here, by means of analytical and numerical analysis, we show that integrable many-body systems can be exploited efficiently for detecting the amplitude of an AC field. Unlike the conventional strategies in using the ground states in critical many-body probes for parameter estimation, we only consider partial access to a subsystem. Due to the periodicity of the dynamics, any local block of the system saturates to a steady state which allows achieving sensing precision well beyond the classical limit, almost reaching the Heisenberg bound. We associate the enhanced quantum precision to closing of the Floquet gap, resembling the features of quantum sensing in the ground state of critical systems. We show that the proposed protocol can also be realized in near-term quantum simulators, e.g. ion-traps, with a limited number of qubits. We show that in such systems a simple block magnetization measurement and a Bayesian inference estimator can achieve very high precision AC field sensing.

Dualities in quantum integrable many-body systems and integrable probabilities. Part I

In this study we map the dualities observed in the framework of integrable probabilities into the dualities familiar in a realm of integrable many-body systems. The dualities between the pairs of stochastic processes involve one representative from Macdonald-Schur family, while the second representative is from stochastic higher spin six-vertex model of TASEP family. We argue that these dualities are counterparts and generalizations of the familiar quantum-quantum (QQ) dualities between pairs of integrable systems. One integrable system from QQ dual pair belongs to the family of inhomogeneous XXZ spin chains, while the second to the Calogero-Moser-Ruijsenaars-Schneider (CM-RS) family. The wave functions of the Hamiltonian system from CM-RS family are known to be related to solutions to (q)KZ equations at the inhomogeneous spin chain side. When the wave function gets substituted by the measure, bilinear in wave functions, a similar correspondence holds true. As an example, we have elaborated in some details a new duality between the discrete-time inhomogeneous multispecies TASEP model on the circle and the quantum Goldfish model from the RS family. We present the precise map of the inhomogeneous multispecies TASEP and 5-vertex model to the trigonometric and rational Goldfish models respectively, where the TASEP local jump rates get identified as the coordinates in the Goldfish model. Some comments concerning the relation of dualities in the stochastic processes with the dualities in SUSY gauge models with surface operators included are made.

Hydrodynamic equations for the Ablowitz–Ladik discretization of the nonlinear Schrödinger equation

Ablowitz and Ladik discovered a discretization that preserves the integrability of the nonlinear Schrödinger equation in one dimension. We compute the generalized free energy of this model and determine the generalized Gibbs ensemble averaged fields and their currents. They are linked to the mean-field circular unitary matrix ensemble. The resulting hydrodynamic equations follow the pattern already known from other integrable many-body systems. The discretized modified Korteweg–de-Vries equation is also studied, which turns out to be related to the beta Jacobi log gas.

Introduction to the Special Issue on Emergent Hydrodynamics in Integrable Many-Body Systems

Introduction to the Special Issue on Emergent Hydrodynamics in Integrable Many-Body Systems

On Cherednik and Nazarov-Sklyanin large N limit construction for integrable many-body systems with elliptic dependence on momenta

The infinite number of particles limit in the dual to elliptic Ruijsenaars model (coordinate trigonometric degeneration of quantum double elliptic model) is proposed using the Nazarov-Sklyanin approach. For this purpose we describe double-elliptization of the Cherednik construction. Namely, we derive explicit expression in terms of the Cherednik operators, which reduces to the generating function of Dell commuting Hamiltonians on the space of symmetric functions. Although the double elliptic Cherednik operators do not commute, they can be used for construction of the N → ∞ limit.

Generalized-hydrodynamic approach to inhomogeneous quenches: correlations, entanglement and quantum effects

We give a pedagogical introduction to the generalized hydrodynamic approach to inhomogeneous quenches in integrable many-body quantum systems. We review recent applications of the theory, focusing in particular on two classes of problems: bipartitioning protocols and trap quenches, which represent two prototypical examples of broken translational symmetry in either the system initial state or post-quench Hamiltonian. We report on exact results that have been obtained for generic time-dependent correlation functions and entanglement evolution, and discuss in detail the range of applicability of the theory. Finally, we present some open questions and suggest perspectives on possible future directions.

Neural Network Approach to Construction of Classical Integrable Systems

Integrable systems have provided various insights into physical phenomena and mathematics. The way of constructing many-body integrable systems is limited to few ansatzes for the Lax pair, except for highly inventive findings of conserved quantities. Machine learning techniques have recently been applied to broad physics fields and proven powerful for building non-trivial transformations and potential functions. We here propose a machine learning approach to a systematic construction of classical integrable systems. Given the Hamiltonian or samples in latent space, our neural network simultaneously learns the corresponding natural Hamiltonian in real space and the canonical transformation between the latent space and the real space variables. We also propose a loss function for building integrable systems and demonstrate successful unsupervised learning for the Toda lattice. Our approach enables exploring new integrable systems without any prior knowledge about the canonical transformation or any ansatz for the Lax pair.