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:2402.11080 · Physical Review Research

Lieb-Robinson correlation function for the quantum transverse-field Ising model

Brendan J. Mahoney, Craig S. Lent

Published Jun 17, 2024• 3 citations
arXivHTMLPDF

Abstract

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

DOI
10.1103/physrevresearch.6.023286

Related papers

The Value of the Early-Time Lieb-Robinson Correlation Function for Qubit Arrays

Brendan J. Mahoney, Craig S. Lent

The Lieb-Robinson correlation function is one way to capture the propagation of quantum entanglement and correlations in many-body systems. We consider arrays of qubits described by the tranverse-field Ising model and examine correlations as the expanding front of entanglement first reaches a particular qubit. Rather than a new bound for the correlation function, we calculate its value, both numerically and analytically. A general analytical result is obtained that enables us to analyze very large arrays of qubits. The velocity of the entanglement front saturates to a constant value, for which an analytic expression is derived. At the leading edge of entanglement, the correlation function is well-described by an exponential reduced by the square root of the distance. This analysis is extended to arbitrary arrays with general coupling and topologies. For regular two and three dimensional qubit arrays with near-neighbor coupling we find the saturation values for the direction-dependent Lieb-Robinson velocity. The symmetry of the underlying 2D or 3D lattice is evident in the shape of surfaces of constant entanglement, even as the correlations front expands over hundreds of qubits.

Symmetry•Oct 26, 2022•2 citations

The Lieb-Robinson correlation function for long disordered transverse-field Ising chains

Brendan J. Mahoney, Craig S. Lent

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.

quant-pharXiv:2512.22395•Physical Review Research•Dec 30, 2025•0 citations

The early-time Lieb-Robinson correlation function for qubit arrays

Brendan J. Mahoney, Craig S. Lent

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.

arXiv:2203.07251•Mar 14, 2022•0 citations

Dynamics of correlation spreading in low-dimensional transverse-field Ising models

Ryui Kaneko, Ippei Danshita

We investigate the dynamical spreading of spatial correlations after a quantum quench starting from a magnetically disordered state in the transverse-field Ising model at one (1D) and two spatial dimensions (2D). We analyze specifically the longitudinal and transverse spin-spin correlation functions at equal time with use of several methods. From the comparison of the results in 1D obtained by the linear spin-wave approximation (LSWA) and those obtained by the rigorous analytical approach, we show that the LSWA can asymptotically reproduce the exact group velocity in the limit of strong transverse fields while it fails to capture the detailed time dependence of the correlation functions. By applying the LSWA to the 2D case, in which the rigorous analytical approach is unavailable, we estimate the propagation velocity to be $Ja/(2\ensuremath{\hbar})$ at the strong-field limit, where $J$ is the Ising interaction and $a$ is the lattice spacing. We also utilize the tensor-network method based on the projected-entangled pair states for 2D and quantitatively compute the time evolution of the correlation functions for a relatively short time. Our findings provide useful benchmarks for quantum simulation experiments of correlation spreading and theoretical refinement of the Lieb-Robinson bound in the future.

arXiv:2301.01407•Physical Review A•Aug 1, 2023•13 citations

Limits to Velocity of Signal Propagation in Many-Body Systems: A Quantum Information Perspective

Piotr Wysocki, Jan Chwedeńczuk

The Lieb-Robinson bound (LRB) states that the range and strength of interactions between the constituents of a complex many-body system impose upper limits to how fast the signal can propagate. It manifests in a light conelike growth of correlation function connecting two distant subsystems. Here, we employ the techniques of quantum information to demonstrate that the LRB can be determined from local measurements performed on a single qubit that is connected to a many-body system. This formulation provides an operational recipe for estimating the LRB in complex systems, replacing the measurement of the correlation function with simple single-particle manipulations. We demonstrate the potency of this approach by deriving the upper limit to the speed of signal propagation in the XY spin chain.

arXiv:2405.03751•Physical Review Letters•Jan 17, 2025•4 citations

Spread of Correlations in Long-Range Interacting Quantum Systems

Philipp Hauke, Luca Tagliacozzo

The nonequilibrium response of a quantum many-body system defines its fundamental transport properties and how initially localized quantum information spreads. However, for long-range-interacting quantum systems little is known. We address this issue by analyzing a local quantum quench in the long-range Ising model in a transverse field, where interactions decay as a variable power law with distance ∝r(-α), α>0. Using complementary numerical and analytical techniques, we identify three dynamical regimes: short-range-like with an emerging light cone for α>2, weakly long range for 1<α<2 without a clear light cone but with a finite propagation speed of almost all excitations, and fully nonlocal for α<1 with instantaneous transmission of correlations. This last regime breaks generalized Lieb-Robinson bounds and thus locality. Numerical calculation of the entanglement spectrum demonstrates that the usual picture of propagating quasiparticles remains valid, allowing an intuitive interpretation of our findings via divergences of quasiparticle velocities. Our results may be tested in state-of-the-art trapped-ion experiments.

arXiv:1304.7725•Physical Review Letters•Nov 12, 2013•335 citations
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 is one way to capture the propagation of quantum entanglement and correlations in many-body systems. We consider arrays of qubits described by the tranverse-field Ising model and examine correlations as the expanding front of entanglement first reaches a particular qubit. Rather than a new bound for the correlation function, we calculate its value, both numerically and analytically. A general analytical result is obtained that enables us to analyze very large arrays of qubits. The velocity of the entanglement front saturates to a constant value, for which an analytic expression is derived. At the leading edge of entanglement, the correlation function is well-described by an exponential reduced by the square root of the distance. This analysis is extended to arbitrary arrays with general coupling and topologies. For regular two and three dimensional qubit arrays with near-neighbor coupling we find the saturation values for the direction-dependent Lieb-Robinson velocity. The symmetry of the underlying 2D or 3D lattice is evident in the shape of surfaces of constant entanglement, even as the correlations front expands over hundreds of qubits.
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.
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.
We investigate the dynamical spreading of spatial correlations after a quantum quench starting from a magnetically disordered state in the transverse-field Ising model at one (1D) and two spatial dimensions (2D). We analyze specifically the longitudinal and transverse spin-spin correlation functions at equal time with use of several methods. From the comparison of the results in 1D obtained by the linear spin-wave approximation (LSWA) and those obtained by the rigorous analytical approach, we show that the LSWA can asymptotically reproduce the exact group velocity in the limit of strong transverse fields while it fails to capture the detailed time dependence of the correlation functions. By applying the LSWA to the 2D case, in which the rigorous analytical approach is unavailable, we estimate the propagation velocity to be $Ja/(2\ensuremath{\hbar})$ at the strong-field limit, where JJJ is the Ising interaction and aaa is the lattice spacing. We also utilize the tensor-network method based on the projected-entangled pair states for 2D and quantitatively compute the time evolution of the correlation functions for a relatively short time. Our findings provide useful benchmarks for quantum simulation experiments of correlation spreading and theoretical refinement of the Lieb-Robinson bound in the future.
The Lieb-Robinson bound (LRB) states that the range and strength of interactions between the constituents of a complex many-body system impose upper limits to how fast the signal can propagate. It manifests in a light conelike growth of correlation function connecting two distant subsystems. Here, we employ the techniques of quantum information to demonstrate that the LRB can be determined from local measurements performed on a single qubit that is connected to a many-body system. This formulation provides an operational recipe for estimating the LRB in complex systems, replacing the measurement of the correlation function with simple single-particle manipulations. We demonstrate the potency of this approach by deriving the upper limit to the speed of signal propagation in the XY spin chain.
The nonequilibrium response of a quantum many-body system defines its fundamental transport properties and how initially localized quantum information spreads. However, for long-range-interacting quantum systems little is known. We address this issue by analyzing a local quantum quench in the long-range Ising model in a transverse field, where interactions decay as a variable power law with distance ∝r(-α), α>0. Using complementary numerical and analytical techniques, we identify three dynamical regimes: short-range-like with an emerging light cone for α>2, weakly long range for 1<α<2 without a clear light cone but with a finite propagation speed of almost all excitations, and fully nonlocal for α<1 with instantaneous transmission of correlations. This last regime breaks generalized Lieb-Robinson bounds and thus locality. Numerical calculation of the entanglement spectrum demonstrates that the usual picture of propagating quasiparticles remains valid, allowing an intuitive interpretation of our findings via divergences of quasiparticle velocities. Our results may be tested in state-of-the-art trapped-ion experiments.