We show that the kinetic theory of quantum and classical Calogero particles reduces to the free-particle Boltzmann equation. We reconcile this simple emergent behaviour with the strongly interacting character of the model by developing a Bethe-Lax co
rrespondence in the classical case. This demonstrates explicitly that the freely propagating degrees of freedom are not bare particles, but rather quasiparticles corresponding to eigenvectors of the Lax matrix. We apply the resulting kinetic theory to classical Calogero particles in external trapping potentials and find excellent agreement with numerical simulations in all cases, both for harmonic traps that preserve integrability and exhibit perfect revivals, and for anharmonic traps that break microscopic integrability. Our framework also yields a simple description of multi-soliton solutions in a harmonic trap, with solitons corresponding to sharp peaks in the quasiparticle density. Extensions to quantum systems of Calogero particles are discussed.
In this paper, we compute exactly the average density of a harmonically confined Riesz gas of $N$ particles for large $N$ in the presence of a hard wall. In this Riesz gas, the particles repel each other via a pairwise interaction that behaves as $|x
_i - x_j|^{-k}$ for $k>-2$, with $x_i$ denoting the position of the $i^{rm th}$ particle. This density can be classified into three different regimes of $k$. For $k geq 1$, where the interactions are effectively short-ranged, the appropriately scaled density has a finite support over $[-l_k(w),w]$ where $w$ is the scaled position of the wall. While the density vanishes at the left edge of the support, it approaches a nonzero constant at the right edge $w$. For $-1<k<1$, where the interactions are weakly long-ranged, we find that the scaled density is again supported over $[-l_k(w),w]$. While it still vanishes at the left edge of the support, it diverges at the right edge $w$ algebraically with an exponent $(k-1)/2$. For $-2<k< -1$, the interactions are strongly long-ranged that leads to a rather exotic density profile with an extended bulk part and a delta-peak at the wall, separated by a hole in between. Exactly at $k=-1$ the hole disappears. For $-2<k< -1$, we find an interesting first-order phase transition when the scaled position of the wall decreases through a critical value $w=w^*(k)$. For $w<w^*(k)$, the density is a pure delta-peak located at the wall. The amplitude of the delta-peak plays the role of an order parameter which jumps to the value $1$ as $w$ is decreased through $w^*(k)$. Our analytical results are in very good agreement with our Monte-Carlo simulations.
We show that the one dimensional discrete nonlinear Schrodinger chain (DNLS) at finite temperature has three different dynamical regimes (ultra-low, low and high temperature regimes). This has been established via (i) one point macroscopic thermodyna
mic observables (temperature $T$ , energy density $epsilon$ and the relationship between them), (ii) emergence and disappearance of an additional almost conserved quantity (total phase difference) and (iii) classical out-of-time-ordered correlators (OTOC) and related quantities (butterfly speed and Lyapunov exponents). The crossover temperatures $T_{textit{l-ul}}$ (between low and ultra-low temperature regimes) and $T_{textit{h-l}}$ (between high and low temperature regimes) extracted from these three different approaches are consistent with each other. The analysis presented here is an important step forward towards the understanding of DNLS which is ubiquitous in many fields and has a non-separable Hamiltonian form. Our work also shows that the different methods used here can serve as important tools to identify dynamical regimes in other interacting many body systems.
It is very common in the literature to write down a Markovian quantum master equation in Lindblad form to describe a system with multiple degrees of freedom and weakly connected to multiple thermal baths which can, in general, be at different tempera
tures and chemical potentials. However, the microscopically derived quantum master equation up to leading order in system-bath coupling is of the so-called Redfield form which is known to not preserve complete positivity in most cases. We analytically show that, in such cases, enforcing complete positivity by imposing any Lindblad form, via any further approximation, necessarily leads to either violation of thermalization, or inaccurate coherences in the energy eigenbasis which then cause a violation of local conservation laws in the non-equilibrium steady state (NESS). In other words, a weak system-bath coupling quantum master equation that is completely positive, shows thermalization and preserves local conservation laws in NESS is fundamentally impossible in generic situations. On the other hand, the Redfield equation, although generically not completely positive, shows thermalization, always preserves local conservation laws and gives correct coherences to leading order. We exemplify our analytical results numerically in an interacting open quantum spin system.
We present exact results for the classical version of the Out-of-Time-Order Commutator (OTOC) for a family of power-law models consisting of $N$ particles in one dimension and confined by an external harmonic potential. These particles are interactin
g via power-law interaction of the form $propto sum_{substack{i, j=1 (i eq j)}}^N|x_i-x_j|^{-k}$ $forall$ $k>1$ where $x_i$ is the position of the $i^text{th}$ particle. We present numerical results for the OTOC for finite $N$ at low temperatures and short enough times so that the system is well approximated by the linearized dynamics around the many body ground state. In the large-$N$ limit, we compute the ground-state dispersion relation in the absence of external harmonic potential exactly and use it to arrive at analytical results for OTOC. We find excellent agreement between our analytical results and the numerics. We further obtain analytical results in the limit where only linear and leading nonlinear (in momentum) terms in the dispersion relation are included. The resulting OTOC is in agreement with numerics in the vicinity of the edge of the light cone. We find remarkably distinct features in OTOC below and above $k=3$ in terms of going from non-Airy behaviour ($1<k<3$) to an Airy universality class ($k>3$). We present certain additional rich features for the case $k=2$ that stem from the underlying integrability of the Calogero-Moser model. We present a field theory approach that also assists in understanding certain aspects of OTOC such as the sound speed. Our findings are a step forward towards a more general understanding of the spatio-temporal spread of perturbations in long-range interacting systems.
In this paper, we look at four generalizations of the one dimensional Aubry-Andre-Harper (AAH) model which possess mobility edges. We map out a phase diagram in terms of population imbalance, and look at the system size dependence of the steady state
imbalance. We find non-monotonic behaviour of imbalance with system parameters, which contradicts the idea that the relaxation of an initial imbalance is fixed only by the ratio of number of extended states to number of localized states. We propose that there exists dimensionless parameters, which depend on the fraction of single particle localized states, single particle extended states and the mean participation ratio of these states. These ingredients fully control the imbalance in the long time limit and we present numerical evidence of this claim. Among the four models considered, three of them have interesting duality relations and their location of mobility edges are known. One of the models (next nearest neighbour coupling) has no known duality but mobility edge exists and the model has been experimentally realized. Our findings are an important step forward to understanding non-equilibrium phenomena in a family of interesting models with incommensurate potentials.
We compute exactly the average spatial density for $N$ spinless noninteracting fermions in a $2d$ harmonic trap rotating with a constant frequency $Omega$ in the presence of an additional repulsive central potential $gamma/r^2$. We find that, in the
large $N$ limit, the bulk density has a rich and nontrivial profile -- with a hole at the center of the trap and surrounded by a multi-layered wedding cake structure. The number of layers depends on $N$ and on the two parameters $Omega$ and $gamma$ leading to a rich phase diagram. Zooming in on the edge of the $k^{rm th}$ layer, we find that the edge density profile exhibits $k$ kinks located at the zeroes of the $k^{rm th}$ Hermite polynomial. Interestingly, in the large $k$ limit, we show that the edge density profile approaches a limiting form, which resembles the shape of a propagating front, found in the unitary evolution of certain quantum spin chains. We also study how a newly formed droplet grows in size on top of the last layer as one changes the parameters.
We theoretically study correlations present deep in the spectrum of many-body-localized systems. An exact analytical expression for the spectral form factor of Poisson spectra can be obtained and is shown to agree well with numerical results on two m
odels exhibiting many-body-localization: a disordered quantum spin chain and a phenomenological $l$-bit model based on the existence of local integrals of motion. We also identify a universal regime that is insensitive to the global density of states as well as spectral edge effects.
Many-body non-equilibrium steady states can be described by a Landau-Ginzburg theory if one allows non-analytic terms in the potential. We substantiate this claim by working out the case of the Ising magnet in contact with a thermal bath and undergoi
ng stochastic reheating: It is reset to a paramagnet at random times. By a combination of stochastic field theory and Monte Carlo simulations, we unveil how the usual $varphi^4$ potential is deformed by non-analytic operators of intrinsic non-equilibrium nature. We demonstrate their infrared relevance at low temperatures by a renormalization-group analysis of the non-equilibrium steady state. The equilibrium ferromagnetic fixed point is thus destabilized by stochastic reheating, and we identify the new non-equilibrium fixed point.
Out-of-time-ordered correlators (OTOC) have been extensively used as a major tool for exploring quantum chaos and also recently, there has been a classical analogue. Studies have been limited to closed systems. In this work, we probe an open classica
l many-body system, more specifically, a spatially extended driven dissipative chain of coupled Duffing oscillators using the classical OTOC to investigate the spread and growth (decay) of an initially localized perturbation in the chain. Correspondingly, we find three distinct types of dynamical behavior, namely the sustained chaos, transient chaos and non-chaotic region, as clearly exhibited by different geometrical shapes in the heat map of OTOC. To quantify such differences, we look at instantaneous speed (IS), finite time Lyapunov exponents (FTLE) and velocity dependent Lyapunov exponents (VDLE) extracted from OTOC. Introduction of these quantities turn out to be instrumental in diagnosing and demarcating different regimes of dynamical behavior. To gain control over open nonlinear systems, it is important to look at the variation of these quantities with respect to parameters. As we tune drive, dissipation and coupling, FTLE and IS exhibit transition between sustained chaos and non-chaotic regimeswith intermediate transient chaos regimes and highly intermittent sustained chaos points. In the limit of zero nonlinearity, we present exact analytical results for the driven dissipative harmonic system and we find that our analytical results can very well describe the non-chaotic regime as well as the late time behavior in the transient regime of the Duffing chain. We believe, this analysis is an important step forward towards understanding nonlinear dynamics, chaos and spatio-temporal spread of perturbations in many-particle open systems.
