Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2205.11544 · Physical review. B./Physical review. B

Real-time correlators in chaotic quantum many-body systems

Adam Nahum, Sthitadhi Roy, Sagar Vijay, Tianci Zhou

Published Dec 26, 2022• 17 citations
arXivHTMLPDF

Abstract

We study real-time correlators $\ensuremath{\langle}\mathcal{O}(\mathbf{x},t)\mathcal{O}(0,0)\ensuremath{\rangle}$ of local operators in chaotic quantum many-body systems. These correlators show universal structure at late times, determined by the geometry of the dominant operator-space Feynman trajectories for the evolving operator $\mathcal{O}(\mathbit{x},t)$. These trajectories, which involve the operator contracting to a point at both the initial and final time, are qualitatively different from those that dominate the out-of-time-order correlation function. In the absence of conservation laws, local correlations decay exponentially: $\ensuremath{\langle}\mathcal{O}(\mathbf{x},t)\mathcal{O}(0,0)\ensuremath{\rangle}\ensuremath{\sim}exp(\ensuremath{-}{s}_{\text{eq}}\phantom{\rule{0.16em}{0ex}}r(\mathbit{v})\phantom{\rule{0.16em}{0ex}}t)$, where $\mathbit{v}=\mathbit{x}/t$ defines a ray in spacetime, and $r(\mathbit{v})$ is a rate function associated with this ray. We express $r(\mathbit{v})$ in terms of cost functions for various spacetime structures. In $1+1\mathrm{D}$ the operator histories can exhibit a phase transition at a critical value ${v}_{c}$ of the ray velocity, which leads to a singular behavior in $r(\mathbit{v})$. At low velocities, the dominant Feynman histories are ``fat'': The operator grows to a size of order ${t}^{\ensuremath{\alpha}}$ (with $\ensuremath{\alpha}=1/2$ in the simplest case) before contracting to a point again. At high velocities the trajectories are ``thin'': The operator always remains of order-one size. In a Haar-random quantum circuit, this transition maps to a simple binding transition for a pair of random walks, which represent the left and right spatial boundaries of the operator. In higher-dimensional systems, thin trajectories always dominate. We discuss the circumstances in which the butterfly velocity ${v}_{B}$ can be deduced from a time-ordered two-point function, rather than the out-of-time ordered correlator. In the random circuit, correlators may also be computed in the framework of an effective Ising-like statistical mechanics model: we describe this calculation, as well as a special feature of the weights in the case of a $1+1\mathrm{D}$ Haar-random brickwork circuit. The present paper addresses lattice models, but also suggests the possibility of morphological phase transitions for real-time Feynman diagrams in continuum quantum field theories.

DOI
10.1103/physrevb.106.224310

Related papers

Many-body quantum chaos and time reversal symmetry

Weijun Wu, Saumya Shivam, Amos Chan

We investigate universal signatures of quantum chaos in the presence of time reversal symmetry (TRS) in generic many body quantum chaotic systems (gMBQCs). We study three classes of minimal models of gMBQCs with TRS, realized through random quantum circuits with (i) local TRS, (ii) global TRS, and (iii) TRS combined with discrete time-translation symmetry. In large local Hilbert space dimension $q$, we derive the emergence of random matrix theory (RMT) universality in the spectral form factor (SFF) at times larger than the Thouless time $t_{\mathrm{Th}}$, which diverges with system sizes in gMBQCs. At times before $t_{\mathrm{Th}}$, we identify universal behaviour beyond RMT by deriving explicit scaling functions of SFF in the thermodynamic limit. In particular, in the simplest non-trivial setting - preserving global TRS while breaking time translation symmetry - we show that the SFF is mapped to the partition function of an emergent classical ferromagnetic Ising model, where the Ising spins correspond to the time-parallel and time-reversed pairings of Feynman paths, and external magnetic fields are induced by TRS-breaking mechanisms. Without relying on the large-$q$ limit, we develop a second independent derivation of the Ising scaling behaviour of SFF using space-time duality and parity symmetric non-Hermitian Ginibre ensembles. Moreover, we show that many body effects originating from time-reversed pairings of Feynman paths manifest in the two-point autocorrelation function (2PAF), the out-of-time-ordered correlator (OTOC), and the partial spectral form factor - quantities sensitive to both eigenvalue and eigenstate correlations. We establish that the fluctuations of 2PAF are governed by an emergent three-state Potts model, leading to an exponential scaling with the operator support size, at a rate set by the three-state Potts model. [See full abstract in the paper]

arXiv:2503.23475•Mar 30, 2025•0 citations

Eigenstate Correlations, the Eigenstate Thermalization Hypothesis, and Quantum Information Dynamics in Chaotic Many-Body Quantum Systems

Dominik Hahn, David J. Luitz, J. T. Chalker

We consider the statistical properties of eigenstates of the time-evolution operator in chaotic many-body quantum systems. Our focus is on correlations between eigenstates that are specific to spatially extended systems and that characterize entanglement dynamics and operator spreading. In order to isolate these aspects of dynamics from those arising as a result of local conservation laws, we consider Floquet systems in which there are no conserved densities. The correlations associated with scrambling of quantum information lie outside the standard framework established by the eigenstate thermalization hypothesis (ETH). In particular, ETH provides a statistical description of matrix elements of local operators between pairs of eigenstates, whereas the aspects of dynamics we are concerned with arise from correlations among sets of four or more eigenstates. We establish the simplest correlation function that captures these correlations and discuss features of its behavior that are expected to be universal at long distances and low energies. We also propose a maximum-entropy ansatz for the joint distribution of a small number n of eigenstates. In the case n = 2 , this ansatz reproduces ETH. For n = 4 it captures both the growth with time of entanglement between subsystems, as characterized by the purity of the time-evolution operator, and also operator spreading, as characterized by the behavior of the out-of-time-order correlator. We test these ideas by comparing results from Monte Carlo sampling of our ansatz with exact diagonalization studies of Floquet quantum circuits. Published by the American Physical Society 2024

arXiv:2309.12982•Physical Review X•Aug 16, 2024•18 citations

Dynamical scaling laws of out-of-time-ordered correlators

Bo-Bo Wei, Gaoyong Sun, Myung-Joong Hwang

The out-of-time-ordered correlator (OTOC) is central to the understanding of information scrambling in quantum many-body systems. In this work, we show that the OTOC in a quantum many-body system close to its critical point obeys dynamical scaling laws which are specified by a few universal critical exponents of the quantum critical point. Such scaling laws of the OTOC imply a universal form for the butterfly velocity of a chaotic system in the quantum critical region and allow one to locate the quantum critical point and extract all universal critical exponents of the quantum phase transitions. We numerically confirm the universality of the butterfly velocity in a chaotic model, namely, the transverse axial next-nearest-neighbor Ising model, and show the feasibility of extracting the critical properties of quantum phase transitions from OTOC using the Lipkin-Meshkov-Glick model.

arXiv:1906.00533•Physical review. B./Physical review. B•Nov 6, 2019•21 citations

Onset of many-body chaos in the O(N) model

Debanjan Chowdhury, Brian G. Swingle

The growth of commutators of initially commuting local operators diagnoses the onset of chaos in quantum many-body systems. We compute such commutators of local field operators with $N$ components in the ($2+1$)-dimensional $O(N)$ nonlinear sigma model to leading order in $1/N$. The system is taken to be in thermal equilibrium at a temperature $T$ above the zero temperature quantum critical point separating the symmetry broken and unbroken phases. The commutator grows exponentially in time with a rate denoted ${\ensuremath{\lambda}}_{L}$. At large $N$ the growth of chaos as measured by ${\ensuremath{\lambda}}_{L}$ is slow because the model is weakly interacting, and we find ${\ensuremath{\lambda}}_{L}\ensuremath{\approx}3.2T/N$. The scaling with temperature is dictated by conformal invariance of the underlying quantum critical point. We also show that operators grow ballistically in space with a ``butterfly velocity'' given by ${v}_{B}/c\ensuremath{\approx}1$ where $c$ is the Lorentz-invariant speed of particle excitations in the system. We briefly comment on the behavior of ${\ensuremath{\lambda}}_{L}$ and ${v}_{B}$ in the neighboring symmetry broken and unbroken phases.

arXiv:1703.02545•Physical review. D/Physical review. D.•Sep 12, 2017•139 citations

Many-Body Quantum Chaos: Analytic Connection to Random Matrix Theory

Pavel Kos, Marko Ljubotina, Tomaž Prosen

A key goal of quantum chaos is to establish a relationship between widely observed universal spectral fluctuations of clean quantum systems and random matrix theory (RMT).Most prominent features of such RMT behavior with respect to a random spectrum, both encompassed in spectral pair correlation function, are statistical suppression of small level spacings (correlation hole) and enhanced stiffness of the spectrum at large spectral ranges.For single particle systems with fully chaotic classical counterparts, the problem has been partly solved by Berry [Proc.R. Soc.London, A400, 229 (1985)] within the so-called diagonal approximation of semiclassical periodic-orbit sums, while the derivation of the full RMT spectral form factor K(t) (Fourier transform of spectral pair correlation function), from semiclassics has been completed by Müller et al. [Phys.Rev. Lett.93, 014103 (2004)].In recent years, the questions of long-time dynamics at high energies, for which the full many-body energy spectrum becomes relevant, are coming at the forefront even for simple many-body quantum systems, such as locally interacting spin chains.Such systems display two universal types of behaviour which are termed as the 'many-body localized phase' and 'ergodic phase'.In the ergodic phase, the spectral fluctuations are excellently described by RMT, even for very simple interactions and in the absence of any external source of disorder.Here we provide a clear theoretical explanation for these observations.We compute K(t) in the leading two orders in t and show its agreement with RMT for non-integrable, time-reversal invariant many-body systems without classical counterparts, a generic example of which are Ising spin 1/2 models in a periodically kicking transverse field.In particular, we relate K(t) to partition functions of a class of twisted classical Ising models on a ring of size t, hence the leading order RMT behavior K(t)2t is a consequence of translation and reflection symmetry of the Ising partition function.

arXiv:1712.02665•Physical Review X•Jun 8, 2018•183 citations

Operator Spreading in Random Unitary Circuits

Adam Nahum, Sagar Vijay, Jeongwan Haah

Random quantum circuits yield minimally structured models for chaotic quantum dynamics, which are able to capture, for example, universal properties of entanglement growth.We provide exact results and coarse-grained models for the spreading of operators by quantum circuits made of Haar-random unitaries.We study both 1 þ 1D and higher dimensions and argue that the coarse-grained pictures carry over to operator spreading in generic many-body systems.In 1 þ 1D, we demonstrate that the out-of-time-order correlator (OTOC) satisfies a biased diffusion equation, which gives exact results for the spatial profile of the OTOC and determines the butterfly speed v B .We find that in 1 þ 1D, the "front" of the OTOC broadens diffusively, with a width scaling in time as t 1=2 .We address fluctuations in the OTOC between different realizations of the random circuit, arguing that they are negligible in comparison to the broadening of the front within a realization.Turning to higher dimensions, we show that the averaged OTOC can be understood exactly via a remarkable correspondence with a purely classical droplet growth problem.This implies that the width of the front of the averaged OTOC scales as t 1=3 in 2 þ 1D and as t 0.240 in 3 þ 1D (exponents of the Kardar-Parisi-Zhang universality class).We support our analytic argument with simulations in 2 þ 1D.We point out that, in two or higher spatial dimensions, the shape of the spreading operator at late times is affected by underlying lattice symmetries and, in general, is not spherical.However, when full spatial rotational symmetry is present in 2 þ 1D, our mapping implies an exact asymptotic form for the OTOC, in terms of the Tracy-Widom distribution.For an alternative perspective on the OTOC in 1 þ 1D, we map it to the partition function of an Ising-like statistical mechanics model.As a result of special structure arising from unitarity, this partition function reduces to a random walk calculation which can be performed exactly.We also use this mapping to give exact results for entanglement growth in 1 þ 1D circuits.

arXiv:1705.08975•Physical Review X•Apr 11, 2018•791 citations
We study real-time correlators \ensuremath⟨O(x,t)O(0,0)\ensuremath⟩\ensuremath{\langle}\mathcal{O}(\mathbf{x},t)\mathcal{O}(0,0)\ensuremath{\rangle}\ensuremath⟨O(x,t)O(0,0)\ensuremath⟩ of local operators in chaotic quantum many-body systems. These correlators show universal structure at late times, determined by the geometry of the dominant operator-space Feynman trajectories for the evolving operator O(\mathbitx,t)\mathcal{O}(\mathbit{x},t)O(\mathbitx,t). These trajectories, which involve the operator contracting to a point at both the initial and final time, are qualitatively different from those that dominate the out-of-time-order correlation function. In the absence of conservation laws, local correlations decay exponentially: \ensuremath⟨O(x,t)O(0,0)\ensuremath⟩\ensuremath∼exp(\ensuremath−seqr(\mathbitv)t)\ensuremath{\langle}\mathcal{O}(\mathbf{x},t)\mathcal{O}(0,0)\ensuremath{\rangle}\ensuremath{\sim}exp(\ensuremath{-}{s}_{\text{eq}}\phantom{\rule{0.16em}{0ex}}r(\mathbit{v})\phantom{\rule{0.16em}{0ex}}t)\ensuremath⟨O(x,t)O(0,0)\ensuremath⟩\ensuremath∼exp(\ensuremath−seq​r(\mathbitv)t), where \mathbitv=\mathbitx/t\mathbit{v}=\mathbit{x}/t\mathbitv=\mathbitx/t defines a ray in spacetime, and r(\mathbitv)r(\mathbit{v})r(\mathbitv) is a rate function associated with this ray. We express r(\mathbitv)r(\mathbit{v})r(\mathbitv) in terms of cost functions for various spacetime structures. In 1+1D1+1\mathrm{D}1+1D the operator histories can exhibit a phase transition at a critical value vc{v}_{c}vc​ of the ray velocity, which leads to a singular behavior in r(\mathbitv)r(\mathbit{v})r(\mathbitv). At low velocities, the dominant Feynman histories are fat'': The operator grows to a size of order ${t}^{\ensuremath{\alpha}}$ (with $\ensuremath{\alpha}=1/2$ in the simplest case) before contracting to a point again. At high velocities the trajectories are thin'': The operator always remains of order-one size. In a Haar-random quantum circuit, this transition maps to a simple binding transition for a pair of random walks, which represent the left and right spatial boundaries of the operator. In higher-dimensional systems, thin trajectories always dominate. We discuss the circumstances in which the butterfly velocity vB{v}_{B}vB​ can be deduced from a time-ordered two-point function, rather than the out-of-time ordered correlator. In the random circuit, correlators may also be computed in the framework of an effective Ising-like statistical mechanics model: we describe this calculation, as well as a special feature of the weights in the case of a 1+1D1+1\mathrm{D}1+1D Haar-random brickwork circuit. The present paper addresses lattice models, but also suggests the possibility of morphological phase transitions for real-time Feynman diagrams in continuum quantum field theories.
We investigate universal signatures of quantum chaos in the presence of time reversal symmetry (TRS) in generic many body quantum chaotic systems (gMBQCs). We study three classes of minimal models of gMBQCs with TRS, realized through random quantum circuits with (i) local TRS, (ii) global TRS, and (iii) TRS combined with discrete time-translation symmetry. In large local Hilbert space dimension qqq, we derive the emergence of random matrix theory (RMT) universality in the spectral form factor (SFF) at times larger than the Thouless time tTht_{\mathrm{Th}}tTh​, which diverges with system sizes in gMBQCs. At times before tTht_{\mathrm{Th}}tTh​, we identify universal behaviour beyond RMT by deriving explicit scaling functions of SFF in the thermodynamic limit. In particular, in the simplest non-trivial setting - preserving global TRS while breaking time translation symmetry - we show that the SFF is mapped to the partition function of an emergent classical ferromagnetic Ising model, where the Ising spins correspond to the time-parallel and time-reversed pairings of Feynman paths, and external magnetic fields are induced by TRS-breaking mechanisms. Without relying on the large-qqq limit, we develop a second independent derivation of the Ising scaling behaviour of SFF using space-time duality and parity symmetric non-Hermitian Ginibre ensembles. Moreover, we show that many body effects originating from time-reversed pairings of Feynman paths manifest in the two-point autocorrelation function (2PAF), the out-of-time-ordered correlator (OTOC), and the partial spectral form factor - quantities sensitive to both eigenvalue and eigenstate correlations. We establish that the fluctuations of 2PAF are governed by an emergent three-state Potts model, leading to an exponential scaling with the operator support size, at a rate set by the three-state Potts model. [See full abstract in the paper]
We consider the statistical properties of eigenstates of the time-evolution operator in chaotic many-body quantum systems. Our focus is on correlations between eigenstates that are specific to spatially extended systems and that characterize entanglement dynamics and operator spreading. In order to isolate these aspects of dynamics from those arising as a result of local conservation laws, we consider Floquet systems in which there are no conserved densities. The correlations associated with scrambling of quantum information lie outside the standard framework established by the eigenstate thermalization hypothesis (ETH). In particular, ETH provides a statistical description of matrix elements of local operators between pairs of eigenstates, whereas the aspects of dynamics we are concerned with arise from correlations among sets of four or more eigenstates. We establish the simplest correlation function that captures these correlations and discuss features of its behavior that are expected to be universal at long distances and low energies. We also propose a maximum-entropy ansatz for the joint distribution of a small number n of eigenstates. In the case n = 2 , this ansatz reproduces ETH. For n = 4 it captures both the growth with time of entanglement between subsystems, as characterized by the purity of the time-evolution operator, and also operator spreading, as characterized by the behavior of the out-of-time-order correlator. We test these ideas by comparing results from Monte Carlo sampling of our ansatz with exact diagonalization studies of Floquet quantum circuits. Published by the American Physical Society 2024
The out-of-time-ordered correlator (OTOC) is central to the understanding of information scrambling in quantum many-body systems. In this work, we show that the OTOC in a quantum many-body system close to its critical point obeys dynamical scaling laws which are specified by a few universal critical exponents of the quantum critical point. Such scaling laws of the OTOC imply a universal form for the butterfly velocity of a chaotic system in the quantum critical region and allow one to locate the quantum critical point and extract all universal critical exponents of the quantum phase transitions. We numerically confirm the universality of the butterfly velocity in a chaotic model, namely, the transverse axial next-nearest-neighbor Ising model, and show the feasibility of extracting the critical properties of quantum phase transitions from OTOC using the Lipkin-Meshkov-Glick model.
The growth of commutators of initially commuting local operators diagnoses the onset of chaos in quantum many-body systems. We compute such commutators of local field operators with NNN components in the (2+12+12+1)-dimensional O(N)O(N)O(N) nonlinear sigma model to leading order in 1/N1/N1/N. The system is taken to be in thermal equilibrium at a temperature TTT above the zero temperature quantum critical point separating the symmetry broken and unbroken phases. The commutator grows exponentially in time with a rate denoted ${\ensuremath{\lambda}}_{L}$. At large NNN the growth of chaos as measured by ${\ensuremath{\lambda}}_{L}$ is slow because the model is weakly interacting, and we find ${\ensuremath{\lambda}}_{L}\ensuremath{\approx}3.2T/N$. The scaling with temperature is dictated by conformal invariance of the underlying quantum critical point. We also show that operators grow ballistically in space with a ``butterfly velocity'' given by ${v}_{B}/c\ensuremath{\approx}1$ where ccc is the Lorentz-invariant speed of particle excitations in the system. We briefly comment on the behavior of ${\ensuremath{\lambda}}_{L}$ and vB{v}_{B}vB​ in the neighboring symmetry broken and unbroken phases.
A key goal of quantum chaos is to establish a relationship between widely observed universal spectral fluctuations of clean quantum systems and random matrix theory (RMT).Most prominent features of such RMT behavior with respect to a random spectrum, both encompassed in spectral pair correlation function, are statistical suppression of small level spacings (correlation hole) and enhanced stiffness of the spectrum at large spectral ranges.For single particle systems with fully chaotic classical counterparts, the problem has been partly solved by Berry [Proc.R. Soc.London, A400, 229 (1985)] within the so-called diagonal approximation of semiclassical periodic-orbit sums, while the derivation of the full RMT spectral form factor K(t) (Fourier transform of spectral pair correlation function), from semiclassics has been completed by Müller et al. [Phys.Rev. Lett.93, 014103 (2004)].In recent years, the questions of long-time dynamics at high energies, for which the full many-body energy spectrum becomes relevant, are coming at the forefront even for simple many-body quantum systems, such as locally interacting spin chains.Such systems display two universal types of behaviour which are termed as the 'many-body localized phase' and 'ergodic phase'.In the ergodic phase, the spectral fluctuations are excellently described by RMT, even for very simple interactions and in the absence of any external source of disorder.Here we provide a clear theoretical explanation for these observations.We compute K(t) in the leading two orders in t and show its agreement with RMT for non-integrable, time-reversal invariant many-body systems without classical counterparts, a generic example of which are Ising spin 1/2 models in a periodically kicking transverse field.In particular, we relate K(t) to partition functions of a class of twisted classical Ising models on a ring of size t, hence the leading order RMT behavior K(t)2t is a consequence of translation and reflection symmetry of the Ising partition function.
Random quantum circuits yield minimally structured models for chaotic quantum dynamics, which are able to capture, for example, universal properties of entanglement growth.We provide exact results and coarse-grained models for the spreading of operators by quantum circuits made of Haar-random unitaries.We study both 1 þ 1D and higher dimensions and argue that the coarse-grained pictures carry over to operator spreading in generic many-body systems.In 1 þ 1D, we demonstrate that the out-of-time-order correlator (OTOC) satisfies a biased diffusion equation, which gives exact results for the spatial profile of the OTOC and determines the butterfly speed v B .We find that in 1 þ 1D, the "front" of the OTOC broadens diffusively, with a width scaling in time as t 1=2 .We address fluctuations in the OTOC between different realizations of the random circuit, arguing that they are negligible in comparison to the broadening of the front within a realization.Turning to higher dimensions, we show that the averaged OTOC can be understood exactly via a remarkable correspondence with a purely classical droplet growth problem.This implies that the width of the front of the averaged OTOC scales as t 1=3 in 2 þ 1D and as t 0.240 in 3 þ 1D (exponents of the Kardar-Parisi-Zhang universality class).We support our analytic argument with simulations in 2 þ 1D.We point out that, in two or higher spatial dimensions, the shape of the spreading operator at late times is affected by underlying lattice symmetries and, in general, is not spherical.However, when full spatial rotational symmetry is present in 2 þ 1D, our mapping implies an exact asymptotic form for the OTOC, in terms of the Tracy-Widom distribution.For an alternative perspective on the OTOC in 1 þ 1D, we map it to the partition function of an Ising-like statistical mechanics model.As a result of special structure arising from unitarity, this partition function reduces to a random walk calculation which can be performed exactly.We also use this mapping to give exact results for entanglement growth in 1 þ 1D circuits.