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

Entropy

Tree Tensor Network Simulation of Dynamical Quantum Phase Transitions in the 2D Transverse-Field Ising Model

Xiangyue Zhang, Dizhou Xie, Yongqiang Li

Published Apr 26, 2026• 0 citations
Publisher

Abstract

The discovery of dynamical quantum phase transitions (DQPTs) has fundamentally challenged the traditional view that phase transitions only occur in thermal equilibrium. Experimental platforms and 1D numerical methods, like matrix product states (MPS), have made great progress. However, exploring true 2D DQPTs remains difficult due to finite-size limitations and the geometric biases of quasi-1D cylinder mappings. Here, we bypass these limitations by deploying a tree tensor network (TTN) approach. This allows us to directly compute the quench dynamics of the transverse-field Ising model (TFIM) on an open 2D square lattice. Because the TTN architecture naturally mirrors 2D lattice connectivity, we can extract the global Loschmidt echo. Our simulations reveal that while deep quenches yield standard DQPTs, quenching within the ferromagnetic phase produces an anomalous dynamical response. In this regime, the rate function exhibits sharp non-analytic peaks even as the macroscopic order parameter maintains its initial sign. This decoupled behavior strongly indicates that local spin excitations drive 2D DQPTs, rather than the macroscopic domain-wall motions seen in 1D chains. These results provide a quantitative numerical baseline for understanding non-equilibrium quantum matter in higher dimensions.

DOI
10.3390/e28050495

Related papers

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

Dynamical quantum phase transitions on random networks

Tomohiro Hashizume, Felix Herbort, Joseph Tindall +1 more

Abstract We investigate two types of dynamical quantum phase transitions (DQPTs) in the transverse-field Ising model on ensembles of Erdős–Rényi networks of size N. These networks consist of vertices connected randomly with probability p ( 0 < p ⩽ 1 ). Using analytical derivations and numerical techniques, we compare the characteristics of the transitions for p < 1 against the fully connected network (p = 1). We analytically show that the overlap between the wave function after a quench and the wave function of the fully connected network after the same quench deviates by at most O ( N − 1 / 2 ) . For a DQPT defined by an order parameter, the critical point remains unchanged for all p. For a DQPT defined by the rate function of the Loschmidt echo, we find that the rate function deviates from the p = 1 limit near vanishing points of the overlap with the initial state, while the critical point remains independent for all p. Our analysis suggests that this divergence arises from persistent non-trivial global many-body correlations absent in the p = 1 limit.

arXiv:2503.04891•New Journal of Physics•Jun 1, 2025•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

Dynamical phase transitions in the two-dimensional transverse-field Ising model

Tomohiro Hashizume, Ian P. McCulloch, Jad C. Halimeh

We investigate two separate notions of dynamical phase transitions in the two-dimensional nearest-neighbor transverse-field Ising model on a square lattice using matrix product states and a hybrid infinite time-evolving block decimation algorithm, where the model is implemented on an infinitely long cylinder with a finite diameter along which periodic boundary conditions are employed. Starting in an ordered initial state, our numerical results suggest that quenches below the dynamical critical point give rise to a ferromagnetic long-time steady state with the Loschmidt return rate exhibiting anomalous cusps even when the order parameter never crosses zero. Within the accessible timescales of our numerics, quenches above the dynamical critical point suggest a paramagnetic long-time steady state with the return rate exhibiting regular cusps connected to zero crossings of the order parameter. Additionally, our simulations indicate that quenching slightly above the dynamical critical point leads to a coexistence region where both anomalous and regular cusps appear in the return rate. Quenches from the disordered phase further confirm our main conclusions. Our work supports the recent finding that anomalous cusps arise only when local spin excitations are the energetically dominant quasiparticles. Our results are accessible in modern Rydberg experiments.

arXiv:1811.09275•Physical Review Research•Mar 31, 2022•63 citations

Quantum phase diagram of the two-dimensional transverse-field Ising model: Unconstrained tree tensor network and mapping analysis

M. Sadrzadeh, Reza Haghshenas, A. Langari

We investigate the ground-state phase diagram of the frustrated transverse-field Ising (TFI) model on the checkerboard lattice (CL), which consists of N\'eel, collinear, quantum paramagnet, and plaquette-valence bond solid (-VBS) phases. We implement a numerical simulation that is based on the recently developed unconstrained tree tensor network ansatz, which systematically improves the accuracy over the conventional methods as it exploits the internal gauge selections. At the highly frustrated region $({J}_{2}={J}_{1})$, we observe a second-order phase transition from the plaquette-VBS state to the paramagnet phase at the critical magnetic-field ${\mathrm{\ensuremath{\Gamma}}}_{c}=0.28$ with the associated critical exponents $\ensuremath{\nu}=1$ and $\ensuremath{\gamma}\ensuremath{\simeq}0.4$, which are obtained within the finite-size scaling analysis on different lattice sizes $N=4\ifmmode\times\else\texttimes\fi{}4,6\ifmmode\times\else\texttimes\fi{}6,8\ifmmode\times\else\texttimes\fi{}8$. The stability of the plaquette-VBS phase at low magnetic fields is examined by a spin-spin correlation function, which verifies the presence of plaquette-VBS at ${J}_{2}={J}_{1}$ and rules out the existence of a N\'eel phase. In addition, our numerical results suggest that the transition from the N\'eel (for ${J}_{2}<{J}_{1})$ to the plaquette-VBS phase is a deconfined phase transition. Moreover, we introduce a mapping, which renders the low-energy effective theory of the TFI on the CL to be the same model on the ${J}_{1}\ensuremath{-}{J}_{2}$ square lattice (SL). We show that the plaquette-VBS phase of the highly frustrated point ${J}_{2}={J}_{1}$ on the CL is mapped to the emergent string-VBS phase on the SL at ${J}_{2}=0.5{J}_{1}$.

arXiv:1901.00278•Physical review. B./Physical review. B•Apr 16, 2019•11 citations

Finite-Temperature Dynamical Phase Diagram of the $2+1$D Quantum Ising Model

Lucas Katschke, Roland C. Farrell, Umberto Borla +2 more

Mapping finite-temperature dynamical phase diagrams of quantum many-body models is a necessary step towards establishing a framework of far-from-equilibrium quantum many-body universality. However, this is quite difficult due, in part, to the severe challenges in representing the volume-law entanglement that is generated under nonequilibrium dynamics at finite temperatures. Here, we address these challenges with an efficient equilibrium quantum Monte Carlo (QMC) framework for computing the finite-temperature dynamical phase diagram. Our method uses energy conservation and the self-thermalizing properties of ergodic quantum systems to determine observables at late times after a quantum quench. We use this technique to chart the dynamical phase diagram of the $2+1$D quantum Ising model generated by quenches of the transverse field in initial thermal states. Our approach allows us to track the evolution of dynamical phases as a function of both the initial temperature and transverse field. Surprisingly, we identify quenches in the ordered phase that cool the system as well as an interval of initial temperatures where it is possible to quench from the paramagnetic (PM) to ferromagnetic (FM) phases. Our method gives access to dynamical properties without explicitly simulating unitary time evolution, and is immediately applicable to other lattice geometries and interacting many-body systems. Finally, we propose a quantum simulation experiment on state-of-the-art digital quantum hardware to directly probe the predicted dynamical phases and their real-time formation.

quant-pharXiv:2602.16772•Feb 20, 2026•0 citations
The discovery of dynamical quantum phase transitions (DQPTs) has fundamentally challenged the traditional view that phase transitions only occur in thermal equilibrium. Experimental platforms and 1D numerical methods, like matrix product states (MPS), have made great progress. However, exploring true 2D DQPTs remains difficult due to finite-size limitations and the geometric biases of quasi-1D cylinder mappings. Here, we bypass these limitations by deploying a tree tensor network (TTN) approach. This allows us to directly compute the quench dynamics of the transverse-field Ising model (TFIM) on an open 2D square lattice. Because the TTN architecture naturally mirrors 2D lattice connectivity, we can extract the global Loschmidt echo. Our simulations reveal that while deep quenches yield standard DQPTs, quenching within the ferromagnetic phase produces an anomalous dynamical response. In this regime, the rate function exhibits sharp non-analytic peaks even as the macroscopic order parameter maintains its initial sign. This decoupled behavior strongly indicates that local spin excitations drive 2D DQPTs, rather than the macroscopic domain-wall motions seen in 1D chains. These results provide a quantitative numerical baseline for understanding non-equilibrium quantum matter in higher dimensions.
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.
Abstract We investigate two types of dynamical quantum phase transitions (DQPTs) in the transverse-field Ising model on ensembles of Erdős–Rényi networks of size N. These networks consist of vertices connected randomly with probability p ( 0 < p ⩽ 1 ). Using analytical derivations and numerical techniques, we compare the characteristics of the transitions for p < 1 against the fully connected network (p = 1). We analytically show that the overlap between the wave function after a quench and the wave function of the fully connected network after the same quench deviates by at most O ( N − 1 / 2 ) . For a DQPT defined by an order parameter, the critical point remains unchanged for all p. For a DQPT defined by the rate function of the Loschmidt echo, we find that the rate function deviates from the p = 1 limit near vanishing points of the overlap with the initial state, while the critical point remains independent for all p. Our analysis suggests that this divergence arises from persistent non-trivial global many-body correlations absent in the p = 1 limit.
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 investigate two separate notions of dynamical phase transitions in the two-dimensional nearest-neighbor transverse-field Ising model on a square lattice using matrix product states and a hybrid infinite time-evolving block decimation algorithm, where the model is implemented on an infinitely long cylinder with a finite diameter along which periodic boundary conditions are employed. Starting in an ordered initial state, our numerical results suggest that quenches below the dynamical critical point give rise to a ferromagnetic long-time steady state with the Loschmidt return rate exhibiting anomalous cusps even when the order parameter never crosses zero. Within the accessible timescales of our numerics, quenches above the dynamical critical point suggest a paramagnetic long-time steady state with the return rate exhibiting regular cusps connected to zero crossings of the order parameter. Additionally, our simulations indicate that quenching slightly above the dynamical critical point leads to a coexistence region where both anomalous and regular cusps appear in the return rate. Quenches from the disordered phase further confirm our main conclusions. Our work supports the recent finding that anomalous cusps arise only when local spin excitations are the energetically dominant quasiparticles. Our results are accessible in modern Rydberg experiments.
We investigate the ground-state phase diagram of the frustrated transverse-field Ising (TFI) model on the checkerboard lattice (CL), which consists of N\'eel, collinear, quantum paramagnet, and plaquette-valence bond solid (-VBS) phases. We implement a numerical simulation that is based on the recently developed unconstrained tree tensor network ansatz, which systematically improves the accuracy over the conventional methods as it exploits the internal gauge selections. At the highly frustrated region (J2=J1)({J}_{2}={J}_{1})(J2​=J1​), we observe a second-order phase transition from the plaquette-VBS state to the paramagnet phase at the critical magnetic-field ${\mathrm{\ensuremath{\Gamma}}}_{c}=0.28$ with the associated critical exponents $\ensuremath{\nu}=1$ and $\ensuremath{\gamma}\ensuremath{\simeq}0.4$, which are obtained within the finite-size scaling analysis on different lattice sizes $N=4\ifmmode\times\else\texttimes\fi{}4,6\ifmmode\times\else\texttimes\fi{}6,8\ifmmode\times\else\texttimes\fi{}8$. The stability of the plaquette-VBS phase at low magnetic fields is examined by a spin-spin correlation function, which verifies the presence of plaquette-VBS at J2=J1{J}_{2}={J}_{1}J2​=J1​ and rules out the existence of a N\'eel phase. In addition, our numerical results suggest that the transition from the N\'eel (for J2<J1){J}_{2}<{J}_{1})J2​<J1​) to the plaquette-VBS phase is a deconfined phase transition. Moreover, we introduce a mapping, which renders the low-energy effective theory of the TFI on the CL to be the same model on the ${J}_{1}\ensuremath{-}{J}_{2}$ square lattice (SL). We show that the plaquette-VBS phase of the highly frustrated point J2=J1{J}_{2}={J}_{1}J2​=J1​ on the CL is mapped to the emergent string-VBS phase on the SL at J2=0.5J1{J}_{2}=0.5{J}_{1}J2​=0.5J1​.
Finite-Temperature Dynamical Phase Diagram of the 2+12+12+1D Quantum Ising Model
Mapping finite-temperature dynamical phase diagrams of quantum many-body models is a necessary step towards establishing a framework of far-from-equilibrium quantum many-body universality. However, this is quite difficult due, in part, to the severe challenges in representing the volume-law entanglement that is generated under nonequilibrium dynamics at finite temperatures. Here, we address these challenges with an efficient equilibrium quantum Monte Carlo (QMC) framework for computing the finite-temperature dynamical phase diagram. Our method uses energy conservation and the self-thermalizing properties of ergodic quantum systems to determine observables at late times after a quantum quench. We use this technique to chart the dynamical phase diagram of the 2+12+12+1D quantum Ising model generated by quenches of the transverse field in initial thermal states. Our approach allows us to track the evolution of dynamical phases as a function of both the initial temperature and transverse field. Surprisingly, we identify quenches in the ordered phase that cool the system as well as an interval of initial temperatures where it is possible to quench from the paramagnetic (PM) to ferromagnetic (FM) phases. Our method gives access to dynamical properties without explicitly simulating unitary time evolution, and is immediately applicable to other lattice geometries and interacting many-body systems. Finally, we propose a quantum simulation experiment on state-of-the-art digital quantum hardware to directly probe the predicted dynamical phases and their real-time formation.