Integrable Many-body Systems Research Articles

Exact Thermalization Dynamics in the "Rule 54" Quantum Cellular Automaton.

We study the out-of-equilibrium dynamics of the quantum cellular automaton known as "Rule 54." For a class of low-entangled initial states, we provide an analytic description of the effect of the global evolution on finite subsystems in terms of simple quantum channels, which gives access to the full thermalization dynamics at the microscopic level. As an example, we provide analytic formulas for the evolution of local observables and Rényi entropies. We show that, in contrast to other known examples of exactly solvable quantum circuits, Rule 54 does not behave as a simple Markovian bath on its own parts, and displays typical nonequilibrium features of interacting integrable many-body quantum systems such as finite relaxation rate and interaction-induced dressing effects. Our study provides a rare example where the full thermalization dynamics can be solved exactly at the microscopic level.

Euler-scale dynamical correlations in integrable systems with fluid motion

We devise an iterative scheme for numerically calculating dynamical two-point correlation functions in integrable many-body systems, in the Eulerian scaling limit. Expressions for these were originally derived in Ref. [1] by combining the fluctuation-dissipation principle with generalized hydrodynamics. Crucially, the scheme is able to address non-stationary, inhomogeneous situations, when motion occurs at the Euler-scale of hydrodynamics. In such situations, in interacting systems, the simple correlations due to fluid modes propagating with the flow receive subtle corrections, which we test. Using our scheme, we study the spreading of correlations in several integrable models from inhomogeneous initial states. For the classical hard-rod model we compare our results with Monte-Carlo simulations and observe excellent agreement at long time-scales, thus providing the first demonstration of validity for the expressions derived in Ref. [1]. We also observe the onset of the Euler-scale limit for the dynamical correlations.

Relaxation to Generalized Gibbs Ensembles in Quadratic Quantum Open Systems

Integrable quantum many-body systems are considered to equilibrate to generalized Gibbs ensembles (GGEs). A general mechanism of such relaxation to GGE has been shown for unitary evolutions of exac...

Weak Integrability Breaking: Chaos with Integrability Signature in Coherent Diffusion.

We study how perturbations affect dynamics of integrable many-body quantum systems, causing transition from integrability to chaos. Looking at spin transport in the Heisenberg chain with impurities we find that in the thermodynamic limit transport gets diffusive already at an infinitesimal perturbation. Small extensive perturbations therefore cause an immediate transition from integrability to chaos. Nevertheless, there is a remnant of integrability encoded in the dependence of the diffusion constant on the impurity density, namely, at small densities it is proportional to the square root of the inverse density, instead of to the inverse density as would follow from Matthiessen's rule. We show that Matthiessen's rule has to be modified in nonballistic systems. Results also highlight a nontrivial role of interacting scattering on a single impurity, and that there is a regime where adding more impurities can actually increase transport.

On the hydrodynamics of unstable excitations

The generalized hydrodynamic (GHD) approach has been extremely successful in describing the out-of-equilibrium properties of a great variety of integrable many-body quantum systems. It naturally extracts the large-scale dynamical degrees of freedom of the system, and is thus a particularly good probe for emergent phenomena. One such phenomenon is the presence of unstable particles, traditionally seen via special analytic structures of the scattering matrix. Because of their finite lifetime and energy threshold, these are especially hard to study. In this paper we apply the GHD approach to a model possessing both unstable excitations and quantum integrability. The largest family of relativistic integrable quantum field theories known to have these features are the homogeneous sine-Gordon models. We consider the simplest non-trivial example of such theories and investigate the effect of an unstable excitation on various physical quantities, both at equilibrium and in the non-equilibrium state arising from the partitioning protocol. The hydrodynamic approach sheds new light onto the physics of the unstable particle, going much beyond its definition via the analytic structure of the scattering matrix, and clarifies its effects both on the equilibrium and out-of-equilibrium properties of the theory. Crucially, within this dynamical perspective, we identify unstable particles as finitely-lived bound states of co-propagating stable particles of different types, and observe how stable populations of unstable particles emerge in large-temperature thermal baths.

Lecture notes on Generalised Hydrodynamics

These are lecture notes for a series of lectures given at the Les Houches Summer School on Integrability in Atomic and Condensed Matter Physics, 30 July to 24 August 2018. The same series of lectures has also been given at the Tokyo Institute of Technology, October 2019. I overview in a pedagogical fashion the main aspects of the theory of generalised hydrodynamics, a hydrodynamic theory for quantum and classical many-body integrable systems. Only very basic knowledge of hydrodynamics and integrable systems is assumed.

Collision rate ansatz for the classical Toda lattice.

We consider a generalized Gibbs ensemble of the classical Toda lattice. We establish that the collision rate ansatz follows because (i) the charge-current susceptibility matrix is symmetric and (ii) the stretch current is proportional to the momentum, hence conserved. The method applies also to other integrable many-body systems, either classical or quantum, provided there is a self-conserved current.

Observation of dynamical fermionization.

The wave function of a Tonks-Girardeau (T-G) gas of strongly interacting bosons in one dimension maps onto the absolute value of the wave function of a noninteracting Fermi gas. Although this fermionization makes many aspects of the two gases identical, their equilibrium momentum distributions are quite different. We observed dynamical fermionization, where the momentum distribution of a T-G gas evolves from bosonic to fermionic after its axial confinement is removed. The asymptotic momentum distribution after expansion in one dimension is the distribution of rapidities, which are the conserved quantities associated with many-body integrable systems. Our measurements agree well with T-G gas theory. We also studied momentum evolution after the trap depth is suddenly changed to a new nonzero value, and we observed the theoretically predicted bosonic-fermionic oscillations.

Quantum-classical duality for Gaudin magnets with boundary

We establish a remarkable relationship between the quantum Gaudin models with boundary and the classical many-body integrable systems of Calogero-Moser type associated with the root systems of classical Lie algebras (B, C and D). We show that under identification of spectra of the Gaudin Hamiltonians HjG with particles velocities q˙j of the classical model all integrals of motion of the latter take zero values. This is the generalization of the quantum-classical duality observed earlier for Gaudin models with periodic boundary conditions and Calogero-Moser models associated with the root system of the type A.

A Remark on the Notion of Independence of Quantum Integrals of Motion in the Thermodynamic Limit

Studies of integrable quantum many-body systems have a long history with an impressive record of success. However, surprisingly enough, an unambiguous definition of quantum integrability remains a matter of an ongoing debate. We contribute to this debate by dwelling upon an important aspect of quantum integrability -- the notion of independence of quantum integrals of motion (QIMs). We point out that a widely accepted definition of functional independence of QIMs is flawed, and suggest a new definition. Our study is motivated by the PXP model -- a model of $N$ spins $1/2$ possessing an extensive number of binary QIMs. The number of QIMs which are independent according to the common definition turns out to be equal to the number of spins, $N$. A common wisdom would then suggest that the system is completely integrable, which is not the case. We discuss the origin of this conundrum and demonstrate how it is resolved when a new definition of independence of QIMs is employed.

Dynamical glass in weakly nonintegrable Klein-Gordon chains.

Integrable many-body systems are characterized by a complete set of preserved actions. Close to an integrable limit, a nonintegrable perturbation creates a coupling network in action space which can be short or long ranged. We analyze the dynamics of observables which become the conserved actions in the integrable limit. We compute distributions of their finite time averages and obtain the ergodization time scale T_{E} on which these distributions converge to δ distributions. We relate T_{E} to the statistics of fluctuation times of the observables, which acquire fat-tailed distributions with standard deviations σ_{τ}^{+} dominating the means μ_{τ}^{+} and establish that T_{E}∼(σ_{τ}^{+})^{2}/μ_{τ}^{+}. The Lyapunov time T_{Λ} (the inverse of the largest Lyapunov exponent) is then compared to the above time scales. We use a simple Klein-Gordon chain to emulate long- and short-range coupling networks by tuning its energy density. For long-range coupling networks T_{Λ}≈σ_{τ}^{+}, which indicates that the Lyapunov time sets the ergodization time, with chaos quickly diffusing through the coupling network. For short-range coupling networks we observe a dynamical glass, where T_{E} grows dramatically by many orders of magnitude and greatly exceeds the Lyapunov time, which satisfies T_{Λ}≲μ_{τ}^{+}. This effect arises from the formation of highly fragmented inhomogeneous distributions of chaotic groups of actions, separated by growing volumes of nonchaotic regions. These structures persist up to the ergodization time T_{E}.

Elliptic solutions to integrable nonlinear equations and many-body systems

We review elliptic solutions to integrable nonlinear partial differential and difference equations (KP, mKP, BKP, Toda) and derive equations of motion for poles of the solutions. The pole dynamics is given by an integrable many-body system (Calogero–Moser, Ruijsenaars–Schneider). The basic tool is the auxiliary linear problems for the wave function which yield equations of motion together with their Lax representation. We also discuss integrals of motion and properties of the spectral curves.

Relaxation to Gaussian generalized Gibbs ensembles in quadratic bosonic systems in the thermodynamic limit.

Integrable quantum many-body systems are considered to equilibrate to generalized Gibbs ensembles (GGEs) characterized by the expectation values of integrals of motion. We study the dynamics of exactly solvable quadratic bosonic systems in the thermodynamic limit, and show a general mechanism for the relaxation to GGEs, in terms of the diagonal singularity. We show analytically and explicitly that a free bosonic system relaxes from a general (not necessarily Gaussian) initial state under certain physical conditions to a Gaussian GGE. We also show the relaxation to a Gaussian GGE in an exactly solvable coupled system, a harmonic oscillator linearly interacting with bosonic reservoirs.

On factorized Lax pairs for classical many-body integrable systems

In this paper, we study the factorization formulae for the Lax matrices of the classical Ruijsenaars–Schneider and Calogero–Moser models. We review the already known results and discuss their possible origins. The first origin comes from the IRF-Vertex relations and the properties of the intertwining matrices. The second origin is based on the Schlesinger transformations generated by modifications of underlying vector bundles. We show that both approaches provide explicit formulae for [Formula: see text]-matrices of the integrable systems in terms of the intertwining matrices (and/or modification matrices). In the end, we discuss the Calogero–Moser models related to classical root systems. The factorization formulae are proposed for a number of special cases.

Global Description of Action-Angle Duality for a Poisson\u2013Lie Deformation of the Trigonometric varvec{mathrm {BC}_n} Sutherland System

Integrable many-body systems of Ruijsenaars–Schneider–van Diejen type displaying action-angle duality are derived by Hamiltonian reduction of the Heisenberg double of the Poisson–Lie group mathrm{SU}(2n). New global models of the reduced phase space are described, revealing non-trivial features of the two systems in duality with one another. For example, after establishing that the symplectic vector space mathbb {C}^nsimeq mathbb {R}^{2n} underlies both global models, it is seen that for both systems the action variables generate the standard torus action on mathbb {C}^n, and the fixed point of this action corresponds to the unique equilibrium positions of the pertinent systems. The systems in duality are found to be non-degenerate in the sense that the functional dimension of the Poisson algebra of their conserved quantities is equal to half the dimension of the phase space. The dual of the deformed Sutherland system is shown to be a limiting case of a van Diejen system.

Typicality of Prethermalization.

Prethermalization refers to the remarkable relaxation behavior which an integrable many-body system in the presence of a weak integrability-breaking perturbation may exhibit: After initial transients have died out, it stays for a long time close to some nonthermal steady state, but on even much larger time scales, it ultimately switches over to the proper thermal equilibrium behavior. By extending Deutsch's conceptual framework from Phys. Rev. A 43, 2046 (1991)PLRAAN1050-294710.1103/PhysRevA.43.2046, we analytically predict that prethermalization is a typical feature for a very general class of such weakly perturbed systems.

Entangling power of time-evolution operators in integrable and nonintegrable many-body systems

The entangling power and operator entanglement entropy are state independent measures of entanglement. Their growth and saturation is examined in the time-evolution operator of quantum many-body systems that can range from the integrable to the fully chaotic. An analytically solvable integrable model of the kicked transverse field Ising chain is shown to have ballistic growth of operator von Neumann entanglement entropy and exponentially fast saturation of the linear entropy with time. Surprisingly a fully chaotic model with longitudinal fields turned on shares the same growth phase, and is consistent with a random matrix model that is also exactly solvable for the linear entropy entanglements. However an examination of the entangling power shows that its largest value is significantly less than the nearly maximal value attained by the nonintegrable one. The importance of long-range spectral correlations, and not just the nearest neighbor spacing, is pointed out in determing the growth of entanglement in nonintegrable systems. Finally an interesting case that displays some features peculiar to both integrable and nonintegrable systems is briefly discussed.

Statistical properties of eigenstate amplitudes in complex quantum systems.

We study the eigenstates of quantum systems with large Hilbert spaces, via their distribution of wave-function amplitudes in a real-space basis. For single-particle "quantum billiards," these real-space amplitudes are known to have Gaussian distribution for chaotic systems. In this work, we formulate and address the corresponding question for many-body lattice quantum systems. For integrable many-body systems, we examine the deviation from Gaussianity and provide evidence that the distribution generically tends toward power-law behavior in the limit of large sizes. We relate the deviation from Gaussianity to the entanglement content of many-body eigenstates. For integrable billiards, we find several cases where the distribution has power-law tails.

Multi-component Ermakov and non-autonomous many-body system connections

It is established that classes of autonomous 3-body and 4-body systems previously set down in the literature may be embedded in solvable, non-autonomous, multi-component Ermakov systems as originally derived in a multi-layer hydrodynamic context. The Ermakov connection is exploited, in particular, to construct integrable non-autonomous many-body type systems parametrised in terms of an arbitrary function J(y / x) where x, y are Jacobi variables.

A geometric viewpoint on generalized hydrodynamics

Generalized hydrodynamics (GHD) is a large-scale theory for the dynamics of many-body integrable systems. It consists of an infinite set of conservation laws for quasi-particles traveling with effective (“dressed”) velocities that depend on the local state. We show that these equations can be recast into a geometric dynamical problem. They are conservation equations with state-independent quasi-particle velocities, in a space equipped with a family of metrics, parametrized by the quasi-particles' type and speed, that depend on the local state. In the classical hard rod or soliton gas picture, these metrics measure the free length of space as perceived by quasi-particles; in the quantum picture, they weigh space with the density of states available to them. Using this geometric construction, we find a general solution to the initial value problem of GHD, in terms of a set of integral equations where time appears explicitly. These integral equations are solvable by iteration and provide an extremely efficient solution algorithm for GHD.