Out-of-time-ordered correlator in the one-dimensional Kuramoto-Sivashinsky and Kardar-Parisi-Zhang equations.

The out-of-time-ordered correlator (OTOC) has emerged as an interesting object in both classical and quantum systems for probing the spatial spread and temporal growth of initially local perturbations in spatially extended chaotic systems. Here, we study the (classical) OTOC and its "light cone" in the nonlinear Kuramoto-Sivashinsky (KS) equation, using extensive numerical simulations. We also show that the linearized KS equation exhibits a qualitatively similar OTOC and light cone, which can be understood via a saddle-point analysis of the linearly unstable modes. Given the deep connection between the KS (deterministic) and the Kardar-Parisi-Zhang (KPZ, which is stochastic) equations, we also explore the OTOC in the KPZ equation. While our numerical results in the KS case are expected to hold in the continuum limit, for the KPZ case it is valid in a discretized version of the KPZ equation. More broadly, our work unravels the intrinsic interplay between noise/instability, nonlinearity, and dissipation in partial differential equations (deterministic or stochastic) through the lens of OTOC.

Many-Body Resonances in the Avalanche Instability of Many-Body Localization.

Many-body localized (MBL) systems fail to reach thermal equilibrium under their own dynamics, even though they are interacting, nonintegrable, and in an extensively excited state. One instability toward thermalization of MBL systems is the so-called "avalanche," where a locally thermalizing rare region is able to spread thermalization through the full system. The spreading of the avalanche may be modeled and numerically studied in finite one-dimensional MBL systems by weakly coupling an infinite-temperature bath to one end of the system. We find that the avalanche spreads primarily via strong many-body resonances between rare near-resonant eigenstates of the closed system. Thus we find and explore a detailed connection between many-body resonances and avalanches in MBL systems.

Many-Body Quantum Chaos and Emergence of Ginibre Ensemble.

We show that non-Hermitian Ginibre random matrix behaviors emerge in spatially extended many-body quantum chaotic systems in the space direction, just as Hermitian random matrix behaviors emerge in chaotic systems in the time direction. Starting with translational invariant models, which can be associated with dual transfer matrices with complex-valued spectra, we show that the linear ramp of the spectral form factor necessitates that the dual spectra have nontrivial correlations, which in fact fall under the universality class of the Ginibre ensemble, demonstrated by computing the level spacing distribution and the dissipative spectral form factor. As a result of this connection, the exact spectral form factor for the Ginibre ensemble can be used to universally describe the spectral form factor for translational invariant many-body quantum chaotic systems in the scaling limit where t and L are large, while the ratio between L and L_{Th}, the many-body Thouless length is fixed. With appropriate variations of Ginibre models, we analytically demonstrate that our claim generalizes to models without translational invariance as well. The emergence of the Ginibre ensemble is a genuine consequence of the strongly interacting and spatially extended nature of the quantum chaotic systems we consider, unlike the traditional emergence of Hermitian random matrix ensembles.

Self-dual quasiperiodic percolation.

How does the percolation transition behave in the absence of quenched randomness? To address this question, we study two nonrandom self-dual quasiperiodic models of square-lattice bond percolation. In both models, the critical point has emergent discrete scale invariance, but none of the additional emergent conformal symmetry of critical random percolation. From the discrete sequences of critical clusters, we find fractal dimensions of D_{f}=1.911943(1) and D_{f}=1.707234(40) for the two models, significantly different from D_{f}=91/48=1.89583... of random percolation. The critical exponents ν, determined through a numerical study of cluster sizes and wrapping probabilities on a torus, are also well below the ν=4/3 of random percolation. While these new models do not appear to belong to a universality class, they demonstrate how the removal of randomness can fundamentally change the critical behavior.

Probing site-resolved correlations in a spin system of ultracold molecules.

Synthetic quantum systems with interacting constituents play an important role in quantum information processing and in explaining fundamental phenomena in many-body physics. Following impressive advances in cooling and trapping techniques, ensembles of ultracold polar molecules have emerged as a promising platform that combines several advantageous properties1-11. These include a large set of internal states with long coherence times12-17 and long-range, anisotropic interactions. These features could enable the exploration of intriguing phases of correlated quantum matter, such as topological superfluids18, quantum spin liquids19, fractional Chern insulators20 and quantum magnets21,22. Probing correlations in these phases is crucial to understanding their properties, necessitating the development of new experimental techniques. Here we use quantum gas microscopy23 to measure the site-resolved dynamics of quantum correlations of polar 23Na87Rb molecules confined in a two-dimensional optical lattice. By using two rotational states of the molecules, we realize a spin-1/2 system with dipolar interactions between particles, producing a quantum spin-exchange model21,22,24,25. We study the evolution of correlations during the thermalization process of an out-of-equilibrium spin system for both spatially isotropic and anisotropic interactions. Furthermore, we examine the correlation dynamics of a spin-anisotropic Heisenberg model engineered from the native spin-exchange model by using periodic microwave pulses26-28. These experiments push the frontier of probing and controlling interacting systems of ultracold molecules, with prospects for exploring new regimes of quantum matter and characterizing entangled states that are useful for quantum computation29,30 and metrology31.

Unusual ergodic and chaotic properties of trapped hard rods.

We investigate ergodicity, chaos, and thermalization for a one-dimensional classical gas of hard rods confined to an external quadratic or quartic trap, which breaks microscopic integrability. To quantify the strength of chaos in this system, we compute its maximal Lyapunov exponent numerically. The approach to thermal equilibrium is studied by considering the time evolution of particle position and velocity distributions and comparing the late-time profiles with the Gibbs state. Remarkably, we find that quadratically trapped hard rods are highly nonergodic and do not resemble a Gibbs state even at extremely long times, despite compelling evidence of chaos for four or more rods. On the other hand, our numerical results reveal that hard rods in a quartic trap exhibit both chaos and thermalization, and equilibrate to a Gibbs state as expected for a nonintegrable many-body system.

Many-body quantum chaos and space-time translational invariance.

We study the consequences of having translational invariance in space and time in many-body quantum chaotic systems. We consider ensembles of random quantum circuits as minimal models of translational invariant many-body quantum chaotic systems. We evaluate the spectral form factor as a sum over many-body Feynman diagrams in the limit of large local Hilbert space dimension q. At sufficiently large t, diagrams corresponding to rigid translations dominate, reproducing the random matrix theory (RMT) behaviour. At finite t, we show that translational invariance introduces additional mechanisms via two novel Feynman diagrams which delay the emergence of RMT. Our analytics suggests the existence of exact scaling forms which describe the approach to RMT behavior in the scaling limit where both t and L are large while the ratio between L and LTh(t), the many-body Thouless length, is fixed. We numerically demonstrate, with simulations of two distinct circuit models, that the resulting scaling functions are universal in the scaling limit.

Emergence of Fermi's Golden Rule.

Fermi's golden rule applies in the limit where an initial quantum state is weakly coupled to a continuum of other final states overlapping its energy. Here we investigate what happens away from this limit, where the set of final states is discrete, with a nonzero mean level spacing; this question arises in a number of recently investigated many-body systems. For different symmetry classes, we analytically and/or numerically calculate the universal crossovers in the average decay of the initial state as the level spacing is varied, with the golden rule emerging in the limit of a continuum. Among the corrections to the exponential decay of the initial state given by Fermi's golden rule is the appearance of the spectral form factor in the longtime regime for small but nonzero level spacing.

Field Theory of Charge Sharpening in Symmetric Monitored Quantum Circuits.

Monitored quantum circuits (MRCs) exhibit a measurement-induced phase transition between area-law and volume-law entanglement scaling. MRCs with a conserved charge additionally exhibit two distinct volume-law entangled phases that cannot be characterized by equilibrium notions of symmetry-breaking or topological order, but rather by the nonequilibrium dynamics and steady-state distribution of charge fluctuations. These include a charge-fuzzy phase in which charge information is rapidly scrambled leading to slowly decaying spatial fluctuations of charge in the steady state, and a charge-sharp phase in which measurements collapse quantum fluctuations of charge without destroying the volume-law entanglement of neutral degrees of freedom. By taking a continuous-time, weak-measurement limit, we construct a controlled replica field theory description of these phases and their intervening charge-sharpening transition in one spatial dimension. We find that the charge fuzzy phase is a critical phase with continuously evolving critical exponents that terminates in a modified Kosterlitz-Thouless transition to the short-range correlated charge-sharp phase. We numerically corroborate these scaling predictions also hold for discrete-time projective-measurement circuit models using large-scale matrix-product state simulations, and discuss generalizations to higher dimensions.

Operator Scaling Dimensions and Multifractality at Measurement-Induced Transitions.

Repeated local measurements of quantum many-body systems can induce a phase transition in their entanglement structure. These measurement-induced phase transitions (MIPTs) have been studied for various types of dynamics, yet most cases yield quantitatively similar critical exponents, making it unclear how many distinct universality classes are present. Here, we probe the properties of the conformal field theories governing these MIPTs using a numerical transfer-matrix method, which allows us to extract the effective central charge, as well as the first few low-lying scaling dimensions of operators at these critical points for (1+1)-dimensional systems. Our results provide convincing evidence that the generic and Clifford MIPTs for qubits lie in different universality classes and that both are distinct from the percolation transition for qudits in the limit of large on-site Hilbert space dimension. For the generic case, we find strong evidence of multifractal scaling of correlation functions at the critical point, reflected in a continuous spectrum of scaling dimensions.

Spatiotemporal spread of perturbations in power-law models at low temperatures: Exact results for classical out-of-time-order correlators.

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 interacting via power-law interaction of the form ∝∑_{i,j=1(i≠j)}^{N}|x_{i}-x_{j}|^{-k}∀k>1 where x_{i} is the position of the ith 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 behavior (13). 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 spatiotemporal spread of perturbations in long-range interacting systems.

Quantum Ergodicity in the Many-Body Localization Problem.

We generalize Page's result on the entanglement entropy of random pure states to the many-body eigenstates of realistic disordered many-body systems subject to long-range interactions. This extension leads to two principal conclusions: first, for increasing disorder the "shells" of constant energy supporting a system's eigenstates fill only a fraction of its full Fock space and are subject to intrinsic correlations absent in synthetic high-dimensional random lattice systems. Second, in all regimes preceding the many-body localization transition individual eigenstates are thermally distributed over these shells. These results, corroborated by comparison to exact diagonalization for an SYK model, are at variance with the concept of "nonergodic extended states" in many-body systems discussed in the recent literature.

Transport, correlations, and chaos in a classical disordered anharmonic chain.

We explore transport properties in a disordered nonlinear chain of classical harmonic oscillators, and thereby identify a regime exhibiting behavior analogous to that seen in quantum many-body-localized systems. Through extensive numerical simulations of this system connected at its ends to heat baths at different temperatures, we computed the heat current and the temperature profile in the nonequilibrium steady state as a function of system size N, disorder strength Δ, and temperature T. The conductivity κ_{N}, obtained for finite length (N), saturates to a value κ_{∞}>0 in the large N limit, for all values of disorder strength Δ and temperature T>0. We show evidence that for any Δ>0 the conductivity goes to zero faster than any power of T in the (T/Δ)→0 limit, and find that the form κ_{∞}∼e^{-B|ln(CΔ/T)|^{3}} fits our data. This form has earlier been suggested by a theory based on the dynamics of multioscillator chaotic islands. The finite-size effect can be κ_{N}<κ_{∞} due to boundary resistance when the bulk conductivity is high (the weak disorder case), or κ_{N}>κ_{∞} due to direct bath-to-bath coupling through bulk localized modes when the bulk is weakly conducting (the strong disorder case). We also present results on equilibrium dynamical correlation functions and on the role of chaos on transport properties. Finally, we explore the differences in the growth and propagation of chaos in the weak and strong chaos regimes by studying the classical version of the out-of-time-ordered commutator.

Scalable Probes of Measurement-Induced Criticality.

We uncover a local order parameter for measurement-induced phase transitions: the average entropy of a single reference qubit initially entangled with the system. Using this order parameter, we identify scalable probes of measurement-induced criticality that are immediately applicable to advanced quantum computing platforms. We test our proposal on a 1+1 dimensional stabilizer circuit model that can be classically simulated in polynomial time. We introduce the concept of a "decoding light cone" to establish the local and efficiently measurable nature of this probe. We also estimate bulk and surface critical exponents for the transition. Developing scalable probes of measurement-induced criticality in more general models may be a useful application of noisy intermediate scale quantum devices, as well as point to more efficient realizations of fault-tolerant quantum computation.

Localization as an Entanglement Phase Transition in Boundary-Driven Anderson Models.

The Anderson localization transition is one of the most well studied examples of a zero temperature quantum phase transition. On the other hand, many open questions remain about the phenomenology of disordered systems driven far out of equilibrium. Here we study the localization transition in the prototypical three-dimensional, noninteracting Anderson model when the system is driven at its boundaries to induce a current carrying nonequilibrium steady state. Recently we showed that the diffusive phase of this model exhibits extensive mutual information of its nonequilibrium steady-state density matrix. We show that this extensive scaling persists in the entanglement and at the localization critical point, before crossing over to a short-range (area-law) scaling in the localized phase. We introduce an entanglement witness for fermionic states that we name the mutual coherence, which, for fermionic Gaussian states, is also a lower bound on the mutual information. Through a combination of analytical arguments and numerics, we determine the finite-size scaling of the mutual coherence across the transition. These results further develop the notion of entanglement phase transitions in open systems, with direct implications for driven many-body localized systems, as well as experimental studies of driven-disordered systems.

Bad metallic transport in a cold atom Fermi-Hubbard system.

Strong interactions in many-body quantum systems complicate the interpretation of charge transport in such materials. To shed light on this problem, we study transport in a clean quantum system: ultracold lithium-6 in a two-dimensional optical lattice, a testing ground for strong interaction physics in the Fermi-Hubbard model. We determine the diffusion constant by measuring the relaxation of an imposed density modulation and modeling its decay hydrodynamically. The diffusion constant is converted to a resistivity by using the Nernst-Einstein relation. That resistivity exhibits a linear temperature dependence and shows no evidence of saturation, two characteristic signatures of a bad metal. The techniques we developed in this study may be applied to measurements of other transport quantities, including the optical conductivity and thermopower.

Light-Cone Spreading of Perturbations and the Butterfly Effect in a Classical Spin Chain.

We find that the effects of a localized perturbation in a chaotic classical many-body system-the classical Heisenberg chain at infinite temperature-spread ballistically with a finite speed even when the local spin dynamics is diffusive. We study two complementary aspects of this butterfly effect: the rapid growth of the perturbation, and its simultaneous ballistic (light-cone) spread, as characterized by the Lyapunov exponents and the butterfly speed, respectively. We connect this to recent studies of the out-of-time-ordered commutators (OTOC), which have been proposed as an indicator of chaos in a quantum system. We provide a straightforward identification of the OTOC with a natural correlator in our system and demonstrate that many of its interesting qualitative features are present in the classical system. Finally, by analyzing the scaling forms, we relate the growth, spread, and propagation of the perturbation with the growth of one-dimensional interfaces described by the Kardar-Parisi-Zhang equation.

Weyl Semimetal to Metal Phase Transitions Driven by Quasiperiodic Potentials.

We explore the stability of three-dimensional Weyl and Dirac semimetals subject to quasiperiodic potentials. We present numerical evidence that the semimetal is stable for weak quasiperiodic potentials, despite being unstable for weak random potentials. As the quasiperiodic potential strength increases, the semimetal transitions to a metal, then to an "inverted" semimetal, and then finally to a metal again. The semimetal and metal are distinguished by the density of states at the Weyl point, as well as by level statistics, transport, and the momentum-space structure of eigenstates near the Weyl point. The critical properties of the transitions in quasiperiodic systems differ from those in random systems: we do not find a clear critical scaling regime in energy; instead, at the quasiperiodic transitions, the density of states appears to jump abruptly (and discontinuously to within our resolution).

Obtaining highly excited eigenstates of the localized XX chain via DMRG-X.

We benchmark a variant of the recently introduced density matrix renormalization group (DMRG)-X algorithm against exact results for the localized random field XX chain. We find that the eigenstates obtained via DMRG-X exhibit a highly accurate l-bit description for system sizes much bigger than the direct, many-body, exact diagonalization in the spin variables is able to access. We take advantage of the underlying free fermion description of the XX model to accurately test the strengths and limitations of this algorithm for large system sizes. We discuss the theoretical constraints on the performance of the algorithm from the entanglement properties of the eigenstates, and its actual performance at different values of disorder. A small but significant improvement to the algorithm is also presented, which helps significantly with convergence. We find that, at high entanglement, DMRG-X shows a bias towards eigenstates with low entanglement, but can be improved with increased bond dimension. This result suggests that one must be careful when applying the algorithm for interacting many-body localized spin models near a transition.This article is part of the themed issue 'Breakdown of ergodicity in quantum systems: from solids to synthetic matter'.

Thermal inclusions: how one spin can destroy a many-body localized phase.

Many-body localized (MBL) systems lie outside the framework of statistical mechanics, as they fail to equilibrate under their own quantum dynamics. Even basic features of MBL systems, such as their stability to thermal inclusions and the nature of the dynamical transition to thermalizing behaviour, remain poorly understood. We study a simple central spin model to address these questions: a two-level system interacting with strength J with N≫1 localized bits subject to random fields. On increasing J, the system transitions from an MBL to a delocalized phase on the vanishing scale Jc(N)∼1/N, up to logarithmic corrections. In the transition region, the single-site eigenstate entanglement entropies exhibit bimodal distributions, so that localized bits are either 'on' (strongly entangled) or 'off' (weakly entangled) in eigenstates. The clusters of 'on' bits vary significantly between eigenstates of the same sample, which provides evidence for a heterogeneous discontinuous transition out of the localized phase in single-site observables. We obtain these results by perturbative mapping to bond percolation on the hypercube at small J and by numerical exact diagonalization of the full many-body system. Our results support the arguments that the MBL phase is unstable in systems with short-range interactions and quenched randomness in dimensions d that are high but finite.This article is part of the themed issue 'Breakdown of ergodicity in quantum systems: from solids to synthetic matter'.