The Lieb-Robinson correlation function captures propagation of quantum correlations in a many-body system. We calculate the value of the leading order of the correlation function, not its bound, for a system of interacting qubits at early times. The general analytical result is compared to numerical calculations and is applied to regular qubit lattices in one, two, and three dimensions. The Lieb-Robinson velocity and the approximately exponential leading edge of correlations emerge in the limit of large arrays.
The Lieb-Robinson correlation function is the norm of a commutator between local operators acting on separate subsystems at different times. This provides a useful state-independent measure for characterizing the specifically quantum interaction between spatially separated qubits. The finite propagation velocity for this correlator defines a “light cone” of quantum influence. We calculate the Lieb-Robinson correlation function for one-dimensional qubit arrays described by the transverse-field Ising model. Direct calculations of this correlation function have been limited by the exponential increase in the size of the state space with the number of qubits. We introduce a technique that avoids this barrier by transforming the calculation to a sum over Pauli walks which results in linear scaling with system size. We can then explore propagation in arrays of hundreds of qubits and observe the effects of the quantum phase transition in the system. We observe the emergence of two distinct velocities of propagation: a correlation front velocity, which is affected by the phase transition, and the Lieb-Robinson velocity which is not. The correlation front velocity is equal to the maximum group velocity of single quasiparticle excitations. The Lieb-Robinson velocity describes the extreme leading edge of correlations when the value of the correlation function itself is still very small. For the semi-infinite chain of qubits at the quantum critical point, we derive an analytical result for the correlation function. Published by the American Physical Society 2024
The Lieb-Robinson correlation function can be used to characterize the propagation of quantum information in Ising chains. Considerable work has been done to establish bounds on this correlation function in various circumstances. To actually calculate the value of the correlation function directly typically requires a state space which grows exponentially with system size, and so is intractable for all but relatively small systems. We employ a recently-developed method that enables direct calculation of the Lieb--Robinson correlation function and which scales linearly with system size. This enables the computation for systems with more than a thousand qubits, revealing the propagation of quantum information down the chain. We extend this technique to the problem of Ising chains with randomly disordered coupling strengths. Increasing disorder causes localization of the quantum correlations and halts propagation of quantum information down the chain. We see light-cone isocontours therefore becoming vertical. The decay of correlations down the chain is characterized by a correlation length which in the paramagnetic phase increases quadratically with mean coupling strength and inverse quadratically with disorder level. The correlation length peaks at the phase transition. We map our general expression for the correlation function in the Ising system onto the probability density in the single particle tight-binding Su-Schrieffer-Heeger model for a specific initial condition.
Lieb-Robinson bounds demonstrate the emergence of locality in many-body quantum systems. Intuitively, Lieb-Robinson bounds state that, with local or exponentially decaying interactions, the correlation that can be built up between two sites separated by distance r after a time t decays as exp ( v t − r ) , where v is the emergent Lieb-Robinson velocity. In many problems, it is important to also capture how much of an operator grows to act on r d sites in d spatial dimensions. Perturbation theory and cluster expansion methods suggest that, at short times, these volume-filling operators are suppressed as exp ( − r d ) . We confirm this intuition, showing that, for r > v t , the volume-filling operator is suppressed by exp ( − ( r − v t ) d / ( v t ) d − 1 ) . This closes a conceptual and practical gap between the cluster expansion and the Lieb-Robinson bound. We then present two very different applications of this new bound. Firstly, we obtain improved bounds on the classical computational resources necessary to simulate many-body dynamics with error tolerance ε for any finite time t : as ε becomes sufficiently small, only ε − O ( t d − 1 ) resources are needed. A protocol that likely saturates this bound is given. Secondly, we prove that disorder operators have volume-law suppression near the “solvable (Ising) point” in quantum phases with spontaneous symmetry breaking, which implies a new diagnostic for distinguishing many-body phases of quantum matter.
We derive an extension of the Lieb--Robinson bound that accounts for memory effects in quantum lattice systems evolving under discrete-time non-Markovian dynamics generated by local interactions with finite-dimensional memory subsystems associated with each bond. The global evolution is unitary on an extended Hilbert space
Hext=⨂i=1NHi⊗⨂bondsHM(i), while non-Markovian behavior, arising from repeated coupling to the same bond memory at successive time steps, leads to temporal correlations in the system. In this work, “finite memory” refers specifically to memory subsystems of fixed finite Hilbert space dimension whose reduced dynamics satisfy an exponential mixing condition in operator norm, meaning that correlation functions decay exponentially over time, with a characteristic timescale
τ. More precisely, the memory channels satisfy
∥TMs(X)−Tr(X)σM∥≤ce−s/τ∥X∥, for all memory observables
X, where
σM is a stationary state, and
c>0 depends only on the memory dynamics. For a one-dimensional lattice with nearest-neighbor interactions, we prove that the connected correlation function, defined as
Ci,j(t):=supAi=Ai†,Bj=Bj†∥Ai∥≤1,∥Bj∥≤1∣Tr[(Ai⊗Bj)ρ(t)]−Tr[Aiρi(t)]Tr[Bjρj(t)]∣, satisfies an exponential bound of the form
Ci,j(t)≤C0exp[−ξd(i,j)−veff(τ)t], where
d(i,j)=∣i−j∣ and
ξ>0 is independent of system size and time. The effective velocity
veff(τ) is explicitly bounded in terms of the local interaction strength and the memory mixing parameters, ensuring the system exhibits a finite light cone, meaning that information propagation is confined to a region growing linearly with time, with a well-defined maximum speed. In the limit
τ→∞, the bound no longer guarantees a finite propagation speed, reflecting the transition to a persistent memory regime where correlations do not decay exponentially. As an operational consequence, the entanglement fidelity between two-site reduced states
ρi,j(t) and a maximally entangled state satisfies
Fi,j(t)=⟨Φ+∣ρi,j(t)∣Φ+⟩≤21[1+Ci,j(t)]. These results extend Lieb--Robinson theory to a class of discrete-time quantum dynamics with memory subsystems associated with each bond, enabling resummation of multi-step memory effects into an effective correlation propagation kernel that governs the decay of temporal correlations. This establishes a quantitative relationship between the timescale of memory mixing and the effective speed of information propagation in the system.
A quantum quench is the simplest protocol to investigate nonequilibrium many-body quantum dynamics. Previous studies on the entanglement properties of quenched quantum many-body systems mainly focus on the growth of entanglement entropy. Several rigorous results and phenomenological guiding principles have been established, such as the no-faster-than-linear entanglement growth generated by generic local Hamiltonians and the peculiar logarithmic growth for many-body localized systems. However, little is known about the dynamical behavior of the full entanglement spectrum, which is a refined character closely related to the topological nature of the wave function. Here, we establish a rigorous and universal result for the entanglement spectra of 1D SPT systems evolving out of equilibrium. Our result is derived both for free-fermion SPT systems and interacting ones. For free-fermion systems with AZ symmetries, we prove that the single-particle entanglement gap after quenches obeys essentially the same Lieb-Robinson bound as that on the equal-time correlation, provided that there is no dynamical symmetry breaking. As a notable byproduct, we obtain a new type of Lieb-Robinson velocity which is related to the band dispersion with a complex wave number and reaches the minimum as the maximal (relative) group velocity. Within the framework of tensor networks, i.e., for SPT MPSs evolved by symmetric and trivial MPUs, we also identify a Lieb-Robinson bound on the many-body entanglement gap for general quenched interacting SPT systems. This result suggests high potential of tensor-network approaches for exploring rigorous results on long-time quantum dynamics. Influence of partial symmetry breaking, effects of disorder, and the relaxation property in the long-time limit are also discussed. Our work establishes a paradigm for exploring rigorous results of SPT systems out of equilibrium.