Search

We present a mathematical theory of metastable pure states in closed many-body quantum systems with finite-dimensional Hilbert space. Given a Hamiltonian, a pure state is defined to be metastable when all sufficiently local operators either stabilize the state, or raise its average energy. We prove that short-range entangled metastable states are necessarily eigenstates (scars) of a perturbatively close Hamiltonian. Given any metastable eigenstate of a Hamiltonian, in the presence of perturbations, we prove the presence of prethermal behavior: local correlation functions decay at a rate bounded by a time scale nonperturbatively long in the inverse metastability radius, rather than Fermi's Golden Rule. Inspired by this general theory, we prove that the lifetime of the false vacuum in certain $d$-dimensional quantum models grows at least as fast as $\exp(\epsilon^{-d} \log^{-3} \epsilon^{-1} )$, where $\epsilon\rightarrow 0$ is the relative energy density of the false vacuum; up to logarithms, this lower bound matches, for the first time, explicit calculations using quantum field theory. We identify metastable states at finite energy density in the PXP model, along with exponentially many metastable states in "helical" spin chains and the two-dimensional Ising model. Our inherently quantum formalism reveals precise connections between many problems, including prethermalization, robust quantum scars, and quantum nucleation theory, and applies to systems without known semiclassical and/or field theoretic limits., Comment: 26+34 pages; 6+3 figures

We prove the existence of extensive many-body Hamiltonians with few-body interactions and a many-body mobility edge: all eigenstates below a nonzero energy density are localized in an exponentially small fraction of "energetically allowed configurations" within Hilbert space. Our construction is based on quantum perturbations to a classical low-density parity check code. In principle, it is possible to detect this eigenstate localization by measuring few-body correlation functions in efficiently preparable mixed states., Comment: 5+8 pages, 1+0 figures. Published version

We bound the time necessary to implement the shift unitary on a one-dimensional ring, both using local Hamiltonians and those with power-law interactions. This time is constrained by the Frobenius light cone; hence we prove that (for certain power law exponents) shift unitaries cannot be implemented in the same amount of time needed to prepare long-range Bell pairs. We note an intriguing similarity between the proof of our results, and the hardness of preparing symmetry-protected topological states with symmetry-preserving Hamiltonians., Comment: 17 pages, 1 figure. Published version

We discuss the generic slowing down of quantum dynamics in low energy density states of spatially local Hamiltonians. Beginning with quantum walks of a single particle, we prove that for certain classes of Hamiltonians (deformations of lattice-regularized $H\propto p^{2k}$), the ``butterfly velocity" of particle motion at low energies has an upper bound that must scale as $E^{(2k-1)/2k}$, as expected from dimensional analysis. We generalize these results to obtain bounds on the typical velocities of particles in many-body systems with repulsive interactions, where for certain families of Hubbard-like models we obtain similar scaling., Comment: 12 pages, 0 figures

We show that it is possible to perform Heisenberg-limited metrology on GHZ-like states, in the presence of generic spatially local, possibly strong interactions during the measurement process. An explicit protocol, which relies on single-qubit measurements and feedback based on polynomial-time classical computation, achieves the Heisenberg limit. In one dimension, matrix product state methods can be used to perform this classical calculation, while in higher dimensions the cluster expansion underlies the efficient calculations. The latter approach is based on an efficient classical sampling algorithm for short-time quantum dynamics, which may be of independent interest., Comment: 30+7 pages, 3+0 figures. Presented at the 27th Annual Conference on Quantum Information Processing (QIP 2024) as a talk. Accepted for publication in Quantum

The computational complexity of simulating quantum many-body systems generally scales exponentially with the number of particles. This enormous computational cost prohibits first principles simulations of many important problems throughout science, ranging from simulating quantum chemistry to discovering the thermodynamic phase diagram of quantum materials or high-density neutron stars. We present a classical algorithm that samples from a high-temperature quantum Gibbs state in a computational (product state) basis. The runtime grows polynomially with the number of particles, while error vanishes polynomially. This algorithm provides an alternative strategy to existing quantum Monte Carlo methods for overcoming the sign problem. Our result implies that measurement-based quantum computation on a Gibbs state can provide exponential speed up only at sufficiently low temperature, and further constrains what tasks can be exponentially faster on quantum computers., Comment: 4+4 pages; 0+1 figure

Circuit quantization is an extraordinarily successful theory that describes the behavior of quantum circuits with high precision. The most widely used approach of circuit quantization relies on introducing a classical Lagrangian whose degrees of freedom are either magnetic fluxes or electric charges in the circuit. By combining nonlinear circuit elements (such as Josephson junctions or quantum phase slips), it is possible to build circuits where a standard Lagrangian description (and thus the standard quantization method) does not exist. Inspired by the mathematics of symplectic geometry and graph theory, we address this challenge, and present a Hamiltonian formulation of non-dissipative electrodynamic circuits. The resulting procedure for circuit quantization is independent of whether circuit elements are linear or nonlinear, or if the circuit is driven by external biases. We explain how to re-derive known results from our formalism, and provide an efficient algorithm for quantizing circuits, including those that cannot be quantized using existing methods., Comment: 30 pages, 8 figures

We review the mathematical speed limits on quantum information processing in many-body systems. After the proof of the Lieb-Robinson Theorem in 1972, the past two decades have seen substantial developments in its application to other questions, such as the simulatability of quantum systems on classical or quantum computers, the generation of entanglement, and even the properties of ground states of gapped systems. Moreover, Lieb-Robinson bounds have been extended in non-trivial ways, to demonstrate speed limits in systems with power-law interactions or interacting bosons, and even to prove notions of locality that arise in cartoon models for quantum gravity with all-to-all interactions. We overview the progress which has occurred, highlight the most promising results and techniques, and discuss some central outstanding questions which remain open. To help bring newcomers to the field up to speed, we provide self-contained proofs of the field's most essential results., Comment: review article. 93 pages, 10 figures, 1 table. v2: minor changes

We prove that prethermalization is a generic property of gapped local many-body quantum systems, subjected to small perturbations, in any spatial dimension. More precisely, let $H_0$ be a Hamiltonian, spatially local in $d$ spatial dimensions, with a gap $\Delta$ in the many-body spectrum; let $V$ be a spatially local Hamiltonian consisting of a sum of local terms, each of which is bounded by $\epsilon \ll \Delta$. Then, the approximation that quantum dynamics is restricted to the low-energy subspace of $H_0$ is accurate, in the correlation functions of local operators, for stretched exponential time scale $\tau \sim \exp[(\Delta/\epsilon)^a]$ for any $a<1/(2d-1)$. This result does not depend on whether the perturbation closes the gap. It significantly extends previous rigorous results on prethermalization in models where $H_0$ was frustration-free. We infer the robustness of quantum simulation in low-energy subspaces, the existence of athermal ``scarred" correlation functions in gapped systems subject to generic perturbations, the long lifetime of false vacua in symmetry broken systems, and the robustness of quantum information in non-frustration-free gapped phases with topological order., Comment: 5+35 pages; 1+4 figures, 1+0 tables. v2: added some results, improved presentation

The speed of light $c$ sets a strict upper bound on the speed of information transfer in both classical and quantum systems. In nonrelativistic quantum systems, the Lieb-Robinson Theorem imposes an emergent speed limit $v \hspace{-0.2mm} \ll \hspace{-0.2mm} c$, establishing locality under unitary evolution and constraining the time needed to perform useful quantum tasks. We extend the Lieb-Robinson Theorem to quantum dynamics with measurements. In contrast to the expectation that measurements can arbitrarily violate spatial locality, we find at most an $(M \hspace{-0.5mm} +\hspace{-0.5mm} 1)$-fold enhancement to the speed $v$ of quantum information, provided the outcomes of measurements in $M$ local regions are known. This holds even when classical communication is instantaneous, and extends beyond projective measurements to weak measurements and other nonunitary channels. Our bound is asymptotically optimal, and saturated by existing measurement-based protocols. We tightly constrain the resource requirements for quantum computation, error correction, teleportation, and generating entangled resource states (Bell, GHZ, quantum-critical, Dicke, W, and spin-squeezed states) from short-range-entangled initial states. Our results impose limits on the use of measurements and active feedback to speed up quantum information processing, resolve fundamental questions about the nature of measurements in quantum dynamics, and constrain the scalability of a wide range of proposed quantum technologies., Comment: 24 pages main text + 61 pages SM; main text significantly expanded, includes new Lieb-Robinson bounds, additional discussion, more clarity on main results

We prove that quantum information propagates with a finite velocity in any model of interacting bosons whose (possibly time-dependent) Hamiltonian contains spatially local single-boson hopping terms along with arbitrary local density-dependent interactions. More precisely, with density matrix $\rho \propto \exp[-\mu N]$ (with $N$ the total boson number), ensemble averaged correlators of the form $\langle [A_0,B_r(t)]\rangle $, along with out-of-time-ordered correlators, must vanish as the distance $r$ between two local operators grows, unless $t \ge r/v$ for some finite speed $v$. In one dimensional models, we give a useful extension of this result that demonstrates the smallness of all matrix elements of the commutator $[A_0,B_r(t)]$ between finite density states if $t/r$ is sufficiently small. Our bounds are relevant for physically realistic initial conditions in experimentally realized models of interacting bosons. In particular, we prove that $v$ can scale no faster than linear in number density in the Bose-Hubbard model: this scaling matches previous results in the high density limit. The quantum walk formalism underlying our proof provides an alternative method for bounding quantum dynamics in models with unbounded operators and infinite-dimensional Hilbert spaces, where Lieb-Robinson bounds have been notoriously challenging to prove., Comment: 35 pages, 1 figure. v2: substantially extended and published version

In many-body quantum systems with spatially local interactions, quantum information propagates with a finite velocity, reminiscent of the ``light cone" of relativity. In systems with long-range interactions which decay with distance $r$ as $1/r^\alpha$, however, there are multiple light cones which control different information theoretic tasks. We show an optimal (up to logarithms) ``Frobenius light cone" obeying $t\sim r^{\min(\alpha-1,1)}$ for $\alpha>1$ in one-dimensional power-law interacting systems with finite local dimension: this controls, among other physical properties, the butterfly velocity characterizing many-body chaos and operator growth. We construct an explicit random Hamiltonian protocol that saturates the bound and settles the optimal Frobenius light cone in one dimension. We partially extend our constraints on the Frobenius light cone to a several operator $p$-norms, and show that Lieb-Robinson bounds can be saturated in at most an exponentially small $e^{-\Omega(r)}$ fraction of the many-body Hilbert space., Comment: 17+20 pages, 3+1 figures. v2: expanded and published version

We present a framework for understanding the dynamics of operator size, and bounding the growth of out-of-time-ordered correlators, in models of large-$S$ spins. Focusing on the dynamics of a single spin, we show the finiteness of the Lyapunov exponent in the large-$S$ limit; our bounds are tighter than the best known Lieb-Robinson-type bounds on these systems. We numerically find our upper bound on Lyapunov exponents is within an order of magnitude of numerically computed values in classical and quantum kicked top models. Generalizing our results to coupled large-$S$ spins on lattices, we show that the butterfly velocity, which characterizes the spatial speed of quantum information scrambling, is finite as $S\rightarrow\infty$. We emphasize qualitative differences between operator growth in semiclassical large-spin models, and quantum holographic systems including the Sachdev-Ye-Kitaev model., Comment: 7+3 pages, 2 figures. v2: modified/improved Lieb-Robinson and quantum walk bounds. v3: published version

We study operator growth in a model of $N(N-1)/2$ interacting Majorana fermions, which live on the edges of a complete graph of $N$ vertices. Terms in the Hamiltonian are proportional to the product of $q$ fermions which live on the edges of cycles of length $q$. This model is a cartoon "matrix model": the interaction graph mimics that of a single-trace matrix model, which can be holographically dual to quantum gravity. We prove (non-perturbatively in $1/N$, and without averaging over any ensemble) that the scrambling time of this model is at least of order $\log N$, consistent with the fast scrambling conjecture. We comment on apparent similarities and differences between operator growth in our "matrix model" and in the melonic models., Comment: 17 pages, 1 figure

Large-scale electrical and thermal currents in ordinary metals are well approximated by effective medium theory: global transport properties are governed by the solution to homogenized coupled diffusion equations. In some metals, including the Dirac fluid of nearly charge neutral graphene, microscopic transport is not governed by diffusion, but by a more complicated set of linearized hydrodynamic equations, which form a system of degenerate elliptic equations coupled with the Stokes equation for fluid velocity. In sufficiently inhomogeneous media, these hydrodynamic equations reduce to homogenized diffusion equations. We re-cast the hydrodynamic transport equations as the infimum of a functional over conserved currents, and present a functional framework to model and compute the homogenized diffusion tensor relating electrical and thermal currents to charge and temperature gradients. We generalize to this system two well-known results in homogenization theory: Tartar's proof of local convergence to the homogenized theory in periodic and highly oscillatory media, and sub-additivity of the above functional in random media with highly oscillatory, stationary and ergodic coefficients., Comment: 27 pages

In quantum many-body systems with local interactions, quantum information and entanglement cannot spread outside of a linear light cone, which expands at an emergent velocity analogous to the speed of light. Local operations at sufficiently separated spacetime points approximately commute -- given a many-body state, $\mathcal{O}_x(t) \mathcal{O}_y |\psi\rangle \approx \mathcal{O}_y\mathcal{O}_x(t) |\psi\rangle$ with arbitrarily small errors -- so long as $|x-y|\gtrsim vt$, where $v$ is finite. Yet most non-relativistic physical systems realized in nature have long-range interactions: two degrees of freedom separated by a distance $r$ interact with potential energy $V(r) \propto 1/r^{\alpha}$. In systems with long-range interactions, we rigorously establish a hierarchy of linear light cones: at the same $\alpha$, some quantum information processing tasks are constrained by a linear light cone while others are not. In one spatial dimension, this linear light cone exists for every many-body state when $\alpha>3$ (Lieb-Robinson light cone); for a typical state chosen uniformly at random from the Hilbert space when $\alpha>\frac{5}{2}$ (Frobenius light cone); for every state of a non-interacting system when $\alpha>2$ (free light cone). These bounds apply to time-dependent systems and are optimal up to subalgebraic improvements. Our theorems regarding the Lieb-Robinson and free light cones -- and their tightness -- also generalize to arbitrary dimensions. We discuss the implications of our bounds on the growth of connected correlators and of topological order, the clustering of correlations in gapped systems, and the digital simulation of systems with long-range interactions. In addition, we show that universal quantum state transfer, as well as many-body quantum chaos, are bounded by the Frobenius light cone, and therefore are poorly constrained by all Lieb-Robinson bounds., Comment: 36 pages; 6 figures; v2: revised and expanded introduction, a few extra results; v3: minor revisions. v4: corrections in Section 7

We prove non-perturbative bounds on the time evolution of the probability distribution of operator size in the $q$-local Sachdev-Ye-Kitaev model with $N$ fermions, for any even integer $q>2$ and any positive even integer $N>2q$. If the couplings in the Hamiltonian are independent and identically distributed Rademacher random variables, the infinite temperature many-body Lyapunov exponent is almost surely finite as $N\rightarrow\infty$. In the limit $q \rightarrow \infty$, $N\rightarrow \infty$, $q^{6+\delta}/N \rightarrow 0$, the shape of the size distribution of a growing fermion, obtained by leading order perturbation calculations in $1/N$ and $1/q$, is similar to a distribution that locally saturates our constraints. Our proof is not based on Feynman diagram resummation; instead, we note that the operator size distribution obeys a continuous time quantum walk with bounded transition rates, to which we apply concentration bounds from classical probability theory., Comment: 23 pages. v2: published version

In a locally interacting many-body system, two isolated qubits, separated by a large distance $r$, become correlated and entangled with each other at a time $t \ge r/v$. This finite speed $v$ of quantum information scrambling limits quantum information processing, thermalization and even equilibrium correlations. Yet most experimental systems contain long range power law interactions -- qubits separated by $r$ have potential energy $V(r)\propto r^{-\alpha}$. Examples include the long range Coulomb interactions in plasma ($\alpha=1$) and dipolar interactions between spins ($\alpha=3$). In one spatial dimension, we prove that the speed of quantum scrambling remains finite for sufficiently large $\alpha$. This result parametrically improves previous bounds, compares favorably with recent numerical simulations, and can be realized in quantum simulators with dipolar interactions. Our new mathematical methods lead to improved algorithms for classically simulating quantum systems, and improve bounds on environmental decoherence in experimental quantum information processors., Comment: 4+10 pages; 1+1 figures; v2: corrected an error; v3: published version

Let $A$ and $B$ be local operators in Hamiltonian quantum systems with $N $ degrees of freedom and finite-dimensional Hilbert space. We prove that the commutator norm $\lVert [A(t),B]\rVert$ is upper bounded by a topological combinatorial problem: counting irreducible weighted paths between two points on the Hamiltonian's factor graph. Our bounds sharpen existing Lieb-Robinson bounds by removing extraneous growth. In quantum systems drawn from zero-mean random ensembles with few-body interactions, we prove stronger bounds on the ensemble-averaged out-of-time-ordered correlator $\mathbb{E}\left[ \lVert [A(t),B]\rVert_F^2\right]$. In such quantum systems on Erd\"os-R\'enyi factor graphs, we prove that the scrambling time $t_{\mathrm{s}}$, at which $\lvert [A(t),B]\rVert_F=\mathrm{\Theta}(1)$, is almost surely $t_{\mathrm{s}}=\mathrm{\Omega}(\sqrt{\log N})$; we further prove $t_{\mathrm{s}}=\mathrm{\Omega}(\log N) $ to high order in perturbation theory in $1/N$. We constrain infinite temperature quantum chaos in the $q$-local Sachdev-Ye-Kitaev model at any order in $1/N$; at leading order, our upper bound on the Lyapunov exponent is within a factor of 2 of the known result at any $q>2$. We also speculate on the implications of our theorems for conjectured holographic descriptions of quantum gravity., Comment: 49 pages, 14 figures. v2: published version with errors fixed

We construct non-equilibrium steady states in the Klein-Gordon theory in arbitrary space dimension $d$ following a local quench. We consider the approach where two independently thermalized semi-infinite systems, with temperatures $T_{\rm L}$ and $T_{\rm R}$, are connected along a $d-1$-dimensional hypersurface. A current-carrying steady state, described by thermally distributed modes with temperatures $T_{\rm L}$ and $T_{\rm R}$ for left and right-moving modes, respectively, emerges at late times. The non-equilibrium density matrix is the exponential of a non-local conserved charge. We obtain exact results for the average energy current and the complete distribution of energy current fluctuations. The latter shows that the long-time energy transfer can be described by a continuum of independent Poisson processes, for which we provide the exact weights. We further describe the full time evolution of local observables following the quench. Averages of generic local observables, including the stress-energy tensor, approach the steady state with a power-law in time, where the exponent depends on the initial conditions at the connection hypersurface. We describe boundary conditions and special operators for which the steady state is reached instantaneously on the connection hypersurface. A semiclassical analysis of freely propagating modes yields the average energy current at large distances and late times. We conclude by comparing and contrasting our findings with results for interacting theories and provide an estimate for the timescale governing the crossover to hydrodynamics. As a modification of our Klein-Gordon analysis we also include exact results for free Dirac fermions., Comment: 42 pages, 7 figures