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:2301.01407 · Physical Review A

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

Ryui Kaneko, Ippei Danshita

Published Aug 1, 2023• 13 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physreva.108.023301

Related papers

Tensor-network study of correlation-spreading dynamics in the two-dimensional Bose-Hubbard model

Ryui Kaneko, Ippei Danshita

Abstract Recent developments in analog quantum simulators based on cold atoms and trapped ions call for cross-validating the accuracy of quantum-simulation experiments with use of quantitative numerical methods; however, it is particularly challenging for dynamics of systems with more than one spatial dimension. Here we demonstrate that a tensor-network method running on classical computers is useful for this purpose. We specifically analyze real-time dynamics of the two-dimensional Bose-Hubbard model after a sudden quench starting from the Mott insulator by means of the infinite projected entangled pair state algorithm. Calculated single-particle correlation functions are found to be in good agreement with a recent experiment [Y. Takasu et al., Sci. Adv. 6, eaba9255 (2020)]. By estimating the phase and group velocities from the single-particle and density-density correlation functions, we predict how these velocities vary in the moderate interaction region, which serves as a quantitative benchmark for future experiments.

arXiv:2108.11051•Communications Physics•Mar 21, 2022•30 citations

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

Brendan J. Mahoney, Craig S. Lent

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

arXiv:2402.11080•Physical Review Research•Jun 17, 2024•3 citations

Time evolution of the quantum Ising model in two dimensions using tree tensor networks

Wladislaw Krinitsin, Niklas Tausendpfund, Markus Heyl +2 more

The numerical simulation of two-dimensional quantum many-body systems away from equilibrium constitutes a major challenge for all known computational methods. We investigate the utility of tree tensor network (TTN) states to solve the dynamics of the quantum Ising model in two dimensions. Within the perturbative regime of small transverse fields, TTNs faithfully reproduce analytically known, but nontrivial and physically interesting dynamics of magnetic domains, for lattices up to 16 × 16 sites. Limitations of the method related to the rapid growth of entanglement entropy are explored within more general, paradigmatic quench settings. We provide and discuss comprehensive benchmarks regarding the benefit of GPU acceleration and the impact of using local operator sums on the performance.

arXiv:2505.07612•Physical review. B./Physical review. B•Sep 29, 2025•3 citations

Simulating dynamics of the two-dimensional transverse-field Ising model: a comparative study of large-scale classical numerics

Joseph Vovrosh, Sergi Julià-Farré, Wladislaw Krinitsin +11 more

The quantum dynamics of many-qubit systems is an outstanding problem that has recently driven significant advances in both numerical methods and programmable quantum processing units. In this work, we employ a comprehensive toolbox of state-of-the-art numerical approaches to classically simulate the dynamics of the two-dimensional transverse field Ising model. Our methods include three different tensor network techniques -- matrix product states, tree-tensor networks, and two-dimensional tensor-networks under the belief propagation approximation -- as well as time-dependent variational Monte Carlo with Neural Quantum States. We focus on two paradigmatic dynamical protocols: (i) quantum annealing through a critical point and (ii) post-quench dynamics. Our extensive results show the quantitative predictions of various state-of-the-art numerical methods providing a benchmark for future numerical investigations and experimental studies with the aim to push the limitations on classical and QPUs. In particular, our work connects classical simulability to different regimes associated with quantum dynamics in Rydberg arrays - namely, quasi-adiabatic dynamics, the Kibble-Zurek mechanism, and quantum quenches.

quant-pharXiv:2511.19340•Physical Review Research•Nov 25, 2025•1 citation

Probing the Quench Dynamics of Antiferromagnetic Correlations in a 2D Quantum Ising Spin System

Elmer Guardado-Sanchez, Peter T. Brown, Debayan Mitra +4 more

Simulating the real-time evolution of quantum spin systems far out of equilibrium poses a major theoretical challenge, especially in more than one dimension.We experimentally explore quench dynamics in a two-dimensional Ising spin system with transverse and longitudinal fields.We realize the system with a near unit-occupancy atomic array of over 200 atoms obtained by loading a spin-polarized band insulator of fermionic lithium into an optical lattice and induce short-range interactions by direct excitation to a lowlying Rydberg state.Using site-resolved microscopy, we probe antiferromagnetic correlations in the system after a sudden quench from a paramagnetic state and compare our measurements to numerical calculations using state-of-the-art techniques.We achieve many-body states with longer-range antiferromagnetic correlations by implementing a near-adiabatic quench of the longitudinal field and study the buildup of correlations as we vary the rate with which we change the field.

arXiv:1711.00887•Physical Review X•Jun 18, 2018•195 citations

Spread of correlations in long-range interacting systems

Philipp Hauke, Luca Tagliacozzo

Understanding the dynamics of many-body systems is crucial for understanding, e.g., thermalization or transmission of information. Nevertheless, little is known in the case of quantum systems with long-range interactions. Here, we analyze the long-range Ising model in a transverse field, where interactions decay as a power-law with distance $\propto r^{-\alpha}$, $\alpha>0$. Using complementary numerical and analytical techniques, we identify three dynamical regimes: short-range-like with an emerging light cone for $\alpha>2$; weakly long-range for $1<\alpha<2$ without a clear light cone but with a finite propagation speed of excitations; and fully non-local for $\alpha<1$ with instantaneous transmission of correlations. This last regime breaks generalized Lieb--Robinson bounds. Numerical calculation of the entanglement spectrum demonstrates that the usual picture of propagating quasi-particles remains valid for long-range interactions. This allows an intuitive interpretation in terms of qualitative changes to the spin-wave dispersion, leading to diverging quasi-particle velocities in the long-range regime. Our results may be tested in state-of-the-art trapped-ion experiments.

•Apr 29, 2013•3 citations
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ℏ)Ja/(2\ensuremath{\hbar})Ja/(2\ensuremathℏ) 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.
Abstract Recent developments in analog quantum simulators based on cold atoms and trapped ions call for cross-validating the accuracy of quantum-simulation experiments with use of quantitative numerical methods; however, it is particularly challenging for dynamics of systems with more than one spatial dimension. Here we demonstrate that a tensor-network method running on classical computers is useful for this purpose. We specifically analyze real-time dynamics of the two-dimensional Bose-Hubbard model after a sudden quench starting from the Mott insulator by means of the infinite projected entangled pair state algorithm. Calculated single-particle correlation functions are found to be in good agreement with a recent experiment [Y. Takasu et al., Sci. Adv. 6, eaba9255 (2020)]. By estimating the phase and group velocities from the single-particle and density-density correlation functions, we predict how these velocities vary in the moderate interaction region, which serves as a quantitative benchmark for future experiments.
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 numerical simulation of two-dimensional quantum many-body systems away from equilibrium constitutes a major challenge for all known computational methods. We investigate the utility of tree tensor network (TTN) states to solve the dynamics of the quantum Ising model in two dimensions. Within the perturbative regime of small transverse fields, TTNs faithfully reproduce analytically known, but nontrivial and physically interesting dynamics of magnetic domains, for lattices up to 16 × 16 sites. Limitations of the method related to the rapid growth of entanglement entropy are explored within more general, paradigmatic quench settings. We provide and discuss comprehensive benchmarks regarding the benefit of GPU acceleration and the impact of using local operator sums on the performance.
The quantum dynamics of many-qubit systems is an outstanding problem that has recently driven significant advances in both numerical methods and programmable quantum processing units. In this work, we employ a comprehensive toolbox of state-of-the-art numerical approaches to classically simulate the dynamics of the two-dimensional transverse field Ising model. Our methods include three different tensor network techniques -- matrix product states, tree-tensor networks, and two-dimensional tensor-networks under the belief propagation approximation -- as well as time-dependent variational Monte Carlo with Neural Quantum States. We focus on two paradigmatic dynamical protocols: (i) quantum annealing through a critical point and (ii) post-quench dynamics. Our extensive results show the quantitative predictions of various state-of-the-art numerical methods providing a benchmark for future numerical investigations and experimental studies with the aim to push the limitations on classical and QPUs. In particular, our work connects classical simulability to different regimes associated with quantum dynamics in Rydberg arrays - namely, quasi-adiabatic dynamics, the Kibble-Zurek mechanism, and quantum quenches.
Simulating the real-time evolution of quantum spin systems far out of equilibrium poses a major theoretical challenge, especially in more than one dimension.We experimentally explore quench dynamics in a two-dimensional Ising spin system with transverse and longitudinal fields.We realize the system with a near unit-occupancy atomic array of over 200 atoms obtained by loading a spin-polarized band insulator of fermionic lithium into an optical lattice and induce short-range interactions by direct excitation to a lowlying Rydberg state.Using site-resolved microscopy, we probe antiferromagnetic correlations in the system after a sudden quench from a paramagnetic state and compare our measurements to numerical calculations using state-of-the-art techniques.We achieve many-body states with longer-range antiferromagnetic correlations by implementing a near-adiabatic quench of the longitudinal field and study the buildup of correlations as we vary the rate with which we change the field.
Understanding the dynamics of many-body systems is crucial for understanding, e.g., thermalization or transmission of information. Nevertheless, little is known in the case of quantum systems with long-range interactions. Here, we analyze the long-range Ising model in a transverse field, where interactions decay as a power-law with distance ∝r−α\propto r^{-\alpha}∝r−α, α>0\alpha>0α>0. Using complementary numerical and analytical techniques, we identify three dynamical regimes: short-range-like with an emerging light cone for α>2\alpha>2α>2; weakly long-range for 1<α<21<\alpha<21<α<2 without a clear light cone but with a finite propagation speed of excitations; and fully non-local for α<1\alpha<1α<1 with instantaneous transmission of correlations. This last regime breaks generalized Lieb--Robinson bounds. Numerical calculation of the entanglement spectrum demonstrates that the usual picture of propagating quasi-particles remains valid for long-range interactions. This allows an intuitive interpretation in terms of qualitative changes to the spin-wave dispersion, leading to diverging quasi-particle velocities in the long-range regime. Our results may be tested in state-of-the-art trapped-ion experiments.