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:2005.03104 · SciPost Physics

Quantum quench dynamics in the transverse-field Ising model: A numerical expansion in linked rectangular clusters

Jonas Richter, Tjark Heitmann, Robin Steinigeweg

Published Sep 3, 2020• 28 citations
arXivHTMLPDF

Abstract

We study quantum quenches in the transverse-field Ising model defined on different lattice geometries such as chains, two- and three-leg ladders, and two-dimensional square lattices. Starting from fully polarized initial states, we consider the dynamics of the transverse and the longitudinal magnetization for quenches to weak, strong, and critical values of the transverse field. To this end, we rely on an efficient combination of numerical linked cluster expansions (NLCEs) and a forward propagation of pure states in real time. As a main result, we demonstrate that NLCEs comprising solely rectangular clusters provide a promising approach to study the real-time dynamics of two-dimensional quantum many-body systems directly in the thermodynamic limit. By comparing to existing data from the literature, we unveil that NLCEs yield converged results on time scales which are competitive to other state-of-the-art numerical methods.

DOI
10.21468/scipostphys.9.3.031

Related papers

Hybrid quantum-classical algorithm for the transverse-field Ising model in the thermodynamic limit

Sumeet, Max Hörmann, Kai Phillip Schmidt

We describe a hybrid quantum-classical approach to treat quantum many-body systems in the thermodynamic limit. This is done by combining numerical linked-cluster expansions (NLCE) with the variational quantum eigensolver (VQE). Here, the VQE algorithm is used as a cluster solver within the NLCE. We test our hybrid quantum-classical algorithm ($\mathrm{NLCE}+\mathrm{VQE}$) for the ferromagnetic transverse-field Ising model on the one-dimensional chain and the two-dimensional square lattice. The calculation of ground-state energies on each open cluster demands a modified Hamiltonian variational ansatz for the VQE. One major finding is convergence of $\mathrm{NLCE}+\mathrm{VQE}$ to the conventional NLCE result in the thermodynamic limit when at least $N/2$ circuit repetitions (known as layers) are used in the VQE ansatz for each cluster with $N$ sites. Our approach demonstrates the fruitful connection of techniques known from correlated quantum many-body systems with hybrid algorithms explored on existing quantum-computing devices.

arXiv:2310.07600•Physical review. B./Physical review. B•Oct 15, 2024•9 citations

Computational dynamics of fermion and spin systems

Ian G. White

The study of quantum dynamics is an increasingly important aspect of atomic, molecular and optical physics, as well as condensed-matter physics and other disciplines. It is essential in order to understand non-equilibrium phenomena which do not fit within the framework of equilibrium thermodynamics, such as many-body localization, dynamical quantum phase transitions, and non-equilibrium steady states in driven-dissipative systems. However, the exact calculation of dynamical observables is extremely difficult in general for several reasons, which makes it challenging for theoretical predictions and analyses to keep pace with rapidly advancing experiments. In this thesis I use and develop several computational methods to investigate the dynamics of quantum many-body systems, with the two goals of better understanding the evolution of local observables and correlation functions in spin and Fermion systems as well as creating and testing more effective numerical methods for solving quantum dynamics problems. I explore quantum dynamics by using Hamiltonian parameter quenches in several commonly-used models in AMO physics and elsewhere: the Fermi-Hubbard model, the XXZ model, and the transverse-Ising model. The first represents a system of short-range interacting Fermions, while the latter two represent spin systems with different spin-spin couplings and magnetic field terms. All of these models can be realized in current experiments. First, I derive an exact solution for an extreme interaction quench in the Fermi-Hubbard model, from a strongly-interacting initial state to a noninteracting final Hamiltonian. From this result I show that there are nontrivial transient connected correlations, despite the fact that the system has very high temperature and a noninteracting Hamiltonian. Those features would suppress such correlations in equilibrium, but the dynamical correlations are still significant. Then, I develop and examine two numerical methods for solving dynamics of the spin expectation values and correlation functions in the transverse-Ising and XXZ models. Both methods rely on approximating a thermodynamically large system with a much smaller and more computationally tractable cluster of lattice sites, or a series of clusters. The first method is a dynamical numerical linked cluster expansion (NLCE). NLCEs have already been used successfully on systems in thermodynamic equilibrium, and I extended the technique to directly compute the dynamics of an observable after a quantum quench. I review the general NLCE method, describe my dynamical variant of it, and compute the quench dynamics of spin and correlation observables in the transverse-Ising and XXZ models. The second method, which I term the adaptive-boundary cluster method, relies on recently refined bounds on information propagation in quantum systems to choose a cluster of lattice sites that minimizes the finite size error bound of a dynamical computation. I briefly review these bounds, and describe a new algorithm which uses them to construct a cluster whose boundary is chosen according to a related error estimate. Such an adaptive-boundary cluster has the property that there exists a computable upper bound on the finite size error of observables measured at the cluster center, and this error bound is nearly minimal when compared to that of other clusters of the same size. I then compute the quench dynamics of spin and correlation observables in the transverse-Ising model using the adaptive boundary clusters. I apply this method to both uniform and strongly disordered systems. To evaluate the effectiveness of both dynamical NLCE and adaptive-boundary clusters, I compare the results obtained using them to results from a standard and general purpose numerical method: exact solvers applied a rectangular cluster with periodic boundary conditions. Both of the new methods usually converge faster and are more accurate than the standard method for comparable cluster sizes. This thesis illustrates that choosing numerical methods which take advantage of the geometric and entanglement structure of the initial state and Hamiltonian can substantially enhance our ability to study quantum dynamics, by reducing the computational cost of numerically exact solutions. I conclude with a discussion of various avenues for future research to build on this work.

•Apr 23, 2020•0 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

Quantum dynamics from a numerical linked cluster expansion

Ian G. White, Bhuvanesh Sundar, Kaden R. A. Hazzard

We demonstrate that a numerical linked cluster expansion method is a powerful tool to calculate quantum dynamics. We calculate the dynamics of the magnetization and spin correlations in the two-dimensional transverse field Ising and XXZ models evolved from a product state. Such dynamics are directly probed in ongoing experiments in ultracold atoms, molecules, and ions. We show that a numerical linked cluster expansion gives dramatically more accurate results at short-to-moderate times than exact diagonalization, and simultaneously requires fewer computational resources. More specifically, the cluster expansion frequently produces more accurate results than an exact diagonalization calculation that would require $10^{5}$--$10^{10}$ more computational operations and memory.

arXiv:1710.07696•Oct 20, 2017•15 citations

Efficient quantum cluster algorithms for frustrated transverse field Ising antiferromagnets and Ising gauge theories

Biswas, Sounak, Kedar Damle

Working within the Stochastic Series Expansion (SSE) framework, we construct efficient quantum cluster algorithms for transverse field Ising antiferromagnets on the pyrochlore lattice and the planar pyrochlore lattice, for the fully frustrated square lattice Ising model in a transverse field (dual to the 2+1 dimensional odd Ising gauge theory), and for a transverse field Ising model with multi-spin interactions on the square lattice, which is dual to a 2+1 dimensional even Ising gauge theory (and reduces to the two dimensional quantum loop model in a certain limit). Our cluster algorithms use a microcanonical update procedure that generalizes and exploits the notion of "pre-marked motifs" introduced earlier in the context of a quantum cluster algorithm for triangular lattice transverse field Ising antiferromagnets. We demonstrate that the resulting algorithms are significantly more efficient than the standard link percolation based quantum cluster approach. We also introduce a new canonical update scheme that leads to a further improvement in measurement of some observables arising from its ability to make one-dimensional clusters in the "imaginary time" direction. Finally, we demonstrate that refinements in the choice of premarking strategies can lead to additional improvements in the efficiency of the microcanonical updates. As a first example of the physics that can be studied using these algorithmic developments, we obtain evidence for a power-law ordered intermediate-temperature phase associated with the two-step melting of long-range order in the fully frustrated square lattice transverse field Ising model.

arXiv:1812.05326•Dec 13, 2018•4 citations

Griffiths-McCoy singularities in the dilute transverse-field Ising model: A numerical linked cluster expansion study

Foster Thompson, Rajiv R. P. Singh

We use numerical linked cluster expansions (NLCEs) to study the site-diluted transverse-field Ising model on the square lattice at T=0. NLCE with a self-consistent mean field on the boundary of the clusters is used to obtain the ground-state magnetization, susceptibility, and structure factor as a function of transverse field h and exchange constant J. Adding site dilution to the model turns NLCE into a series expansion in the dilution parameter p. Studying the divergence of the structure factor allows us to establish the phase diagram in the h/J and p plane. By studying the magnetization of the system in a longitudinal field, we investigate the Griffiths-McCoy singularities. We find that the magnetization develops nonlinearities in the Griffiths phase with exponents that vary continuously with h. Additionally, the probability distribution of the local susceptibility develops long tails in the Griffiths phase, which is studied in terms of its moments.

arXiv:1811.05504•Physical review. E•Mar 25, 2019•6 citations
We study quantum quenches in the transverse-field Ising model defined on different lattice geometries such as chains, two- and three-leg ladders, and two-dimensional square lattices. Starting from fully polarized initial states, we consider the dynamics of the transverse and the longitudinal magnetization for quenches to weak, strong, and critical values of the transverse field. To this end, we rely on an efficient combination of numerical linked cluster expansions (NLCEs) and a forward propagation of pure states in real time. As a main result, we demonstrate that NLCEs comprising solely rectangular clusters provide a promising approach to study the real-time dynamics of two-dimensional quantum many-body systems directly in the thermodynamic limit. By comparing to existing data from the literature, we unveil that NLCEs yield converged results on time scales which are competitive to other state-of-the-art numerical methods.
We describe a hybrid quantum-classical approach to treat quantum many-body systems in the thermodynamic limit. This is done by combining numerical linked-cluster expansions (NLCE) with the variational quantum eigensolver (VQE). Here, the VQE algorithm is used as a cluster solver within the NLCE. We test our hybrid quantum-classical algorithm (NLCE+VQE\mathrm{NLCE}+\mathrm{VQE}NLCE+VQE) for the ferromagnetic transverse-field Ising model on the one-dimensional chain and the two-dimensional square lattice. The calculation of ground-state energies on each open cluster demands a modified Hamiltonian variational ansatz for the VQE. One major finding is convergence of NLCE+VQE\mathrm{NLCE}+\mathrm{VQE}NLCE+VQE to the conventional NLCE result in the thermodynamic limit when at least N/2N/2N/2 circuit repetitions (known as layers) are used in the VQE ansatz for each cluster with NNN sites. Our approach demonstrates the fruitful connection of techniques known from correlated quantum many-body systems with hybrid algorithms explored on existing quantum-computing devices.
The study of quantum dynamics is an increasingly important aspect of atomic, molecular and optical physics, as well as condensed-matter physics and other disciplines. It is essential in order to understand non-equilibrium phenomena which do not fit within the framework of equilibrium thermodynamics, such as many-body localization, dynamical quantum phase transitions, and non-equilibrium steady states in driven-dissipative systems. However, the exact calculation of dynamical observables is extremely difficult in general for several reasons, which makes it challenging for theoretical predictions and analyses to keep pace with rapidly advancing experiments. In this thesis I use and develop several computational methods to investigate the dynamics of quantum many-body systems, with the two goals of better understanding the evolution of local observables and correlation functions in spin and Fermion systems as well as creating and testing more effective numerical methods for solving quantum dynamics problems. I explore quantum dynamics by using Hamiltonian parameter quenches in several commonly-used models in AMO physics and elsewhere: the Fermi-Hubbard model, the XXZ model, and the transverse-Ising model. The first represents a system of short-range interacting Fermions, while the latter two represent spin systems with different spin-spin couplings and magnetic field terms. All of these models can be realized in current experiments. First, I derive an exact solution for an extreme interaction quench in the Fermi-Hubbard model, from a strongly-interacting initial state to a noninteracting final Hamiltonian. From this result I show that there are nontrivial transient connected correlations, despite the fact that the system has very high temperature and a noninteracting Hamiltonian. Those features would suppress such correlations in equilibrium, but the dynamical correlations are still significant. Then, I develop and examine two numerical methods for solving dynamics of the spin expectation values and correlation functions in the transverse-Ising and XXZ models. Both methods rely on approximating a thermodynamically large system with a much smaller and more computationally tractable cluster of lattice sites, or a series of clusters. The first method is a dynamical numerical linked cluster expansion (NLCE). NLCEs have already been used successfully on systems in thermodynamic equilibrium, and I extended the technique to directly compute the dynamics of an observable after a quantum quench. I review the general NLCE method, describe my dynamical variant of it, and compute the quench dynamics of spin and correlation observables in the transverse-Ising and XXZ models. The second method, which I term the adaptive-boundary cluster method, relies on recently refined bounds on information propagation in quantum systems to choose a cluster of lattice sites that minimizes the finite size error bound of a dynamical computation. I briefly review these bounds, and describe a new algorithm which uses them to construct a cluster whose boundary is chosen according to a related error estimate. Such an adaptive-boundary cluster has the property that there exists a computable upper bound on the finite size error of observables measured at the cluster center, and this error bound is nearly minimal when compared to that of other clusters of the same size. I then compute the quench dynamics of spin and correlation observables in the transverse-Ising model using the adaptive boundary clusters. I apply this method to both uniform and strongly disordered systems. To evaluate the effectiveness of both dynamical NLCE and adaptive-boundary clusters, I compare the results obtained using them to results from a standard and general purpose numerical method: exact solvers applied a rectangular cluster with periodic boundary conditions. Both of the new methods usually converge faster and are more accurate than the standard method for comparable cluster sizes. This thesis illustrates that choosing numerical methods which take advantage of the geometric and entanglement structure of the initial state and Hamiltonian can substantially enhance our ability to study quantum dynamics, by reducing the computational cost of numerically exact solutions. I conclude with a discussion of various avenues for future research to build on this work.
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.
We demonstrate that a numerical linked cluster expansion method is a powerful tool to calculate quantum dynamics. We calculate the dynamics of the magnetization and spin correlations in the two-dimensional transverse field Ising and XXZ models evolved from a product state. Such dynamics are directly probed in ongoing experiments in ultracold atoms, molecules, and ions. We show that a numerical linked cluster expansion gives dramatically more accurate results at short-to-moderate times than exact diagonalization, and simultaneously requires fewer computational resources. More specifically, the cluster expansion frequently produces more accurate results than an exact diagonalization calculation that would require 10510^{5}105--101010^{10}1010 more computational operations and memory.
Working within the Stochastic Series Expansion (SSE) framework, we construct efficient quantum cluster algorithms for transverse field Ising antiferromagnets on the pyrochlore lattice and the planar pyrochlore lattice, for the fully frustrated square lattice Ising model in a transverse field (dual to the 2+1 dimensional odd Ising gauge theory), and for a transverse field Ising model with multi-spin interactions on the square lattice, which is dual to a 2+1 dimensional even Ising gauge theory (and reduces to the two dimensional quantum loop model in a certain limit). Our cluster algorithms use a microcanonical update procedure that generalizes and exploits the notion of "pre-marked motifs" introduced earlier in the context of a quantum cluster algorithm for triangular lattice transverse field Ising antiferromagnets. We demonstrate that the resulting algorithms are significantly more efficient than the standard link percolation based quantum cluster approach. We also introduce a new canonical update scheme that leads to a further improvement in measurement of some observables arising from its ability to make one-dimensional clusters in the "imaginary time" direction. Finally, we demonstrate that refinements in the choice of premarking strategies can lead to additional improvements in the efficiency of the microcanonical updates. As a first example of the physics that can be studied using these algorithmic developments, we obtain evidence for a power-law ordered intermediate-temperature phase associated with the two-step melting of long-range order in the fully frustrated square lattice transverse field Ising model.
We use numerical linked cluster expansions (NLCEs) to study the site-diluted transverse-field Ising model on the square lattice at T=0. NLCE with a self-consistent mean field on the boundary of the clusters is used to obtain the ground-state magnetization, susceptibility, and structure factor as a function of transverse field h and exchange constant J. Adding site dilution to the model turns NLCE into a series expansion in the dilution parameter p. Studying the divergence of the structure factor allows us to establish the phase diagram in the h/J and p plane. By studying the magnetization of the system in a longitudinal field, we investigate the Griffiths-McCoy singularities. We find that the magnetization develops nonlinearities in the Griffiths phase with exponents that vary continuously with h. Additionally, the probability distribution of the local susceptibility develops long tails in the Griffiths phase, which is studied in terms of its moments.