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:0712.1806 · Physical Review B

Quantum transverse-field Ising model on an infinite tree from matrix product states

Daniel Nagaj, Edward Farhi, Jeffrey Goldstone, Peter W. Shor, Igor Andrade Sylvester

Published Jun 25, 2008• 69 citations
arXivHTMLPDF

Abstract

We give a generalization to an infinite tree geometry of Vidal's infinite time-evolving block decimation (iTEBD) algorithm [G. Vidal, Phys. Rev. Lett. 98, 070201 (2007)] for simulating an infinite line of quantum spins. We numerically investigate the quantum Ising model in a transverse field on the Bethe lattice using the matrix product state ansatz. We observe a second order phase transition, with certain key differences from the transverse field Ising model on an infinite spin chain. We also investigate a transverse field Ising model with a specific longitudinal field. When the transverse field is turned off, this model has a highly degenerate ground state as opposed to the pure Ising model whose ground state is only doubly degenerate.

DOI
10.1103/physrevb.77.214431

Related papers

Efficient simulation of infinite tree tensor network states on the Bethe lattice

Wei Li, Jan von Delft, Tao Xiang

We show that the simple update approach proposed by Jiang et al. [H. C. Jiang, Z. Y. Weng, and T. Xiang, Phys. Rev. Lett. 101, 090603 (2008)] is an efficient and accurate method for determining the infinite tree tensor network states on the Bethe lattice. Ground-state properties of the quantum transverse Ising model and the Heisenberg XXZ model on the Bethe lattice are studied. The transverse Ising model is found to undergo a second-order quantum phase transition with a diverging magnetic susceptibility but a finite correlation length which is upper bounded by $1/\mathrm{ln}(q\ensuremath{-}1)$ even at the transition point ($q$ is the coordinate number of the Bethe lattice). An intuitive explanation on this peculiar ``critical'' phenomenon is given. The XXZ model on the Bethe lattice undergoes a first-order quantum phase transition at the isotropic point. Furthermore, the simple update scheme is found to be related with the Bethe approximation. Finally, by applying the simple update to various tree tensor clusters, we can obtain rather nice and scalable approximations for two-dimensional lattices.

arXiv:1209.2387•Physical Review B•Nov 26, 2012•84 citations

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

Characterizing quantum criticality in the transverse field Ising model with staggered Γ interaction

Ai Min Chen, Dongchang Liu, Xia Jin +2 more

Abstract This study investigates the quantum phase transitions and critical behaviors in the transverse Ising model with staggered Γ interaction by utilizing the infinite-time evolving block decimation method in the infinite matrix product states representation. The ground-state phase diagram shows that the system exhibits a phase transition from the transverse Néel phase to the paramagnetic phase as the transverse magnetic field increases, which are characterized by the local magnetization and non-local transverse Néel order. The chiral order and nematic order are calculated to explore the helical and nematic magnetic structure caused by the Γ interaction. From the spin-spin correlation, a modulated commensurate power-law oscillating decay behavior provides insight in analyzing the ground state characteristics of each phase. In order to study the quantum criticality, the critical exponents β , ν and central charge c are extracted and numerically consistent with the Ising class transition, demonstrating that the critical behavior in our model occurs is independent of the emergent Gamma interaction.

Physica Scripta•Oct 1, 2025•1 citation

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

Xiangyue Zhang, Dizhou Xie, Yongqiang Li

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.

Entropy•Apr 26, 2026•0 citations

Quantum Imaginary Time Evolution on an Infinite 1D Chain

Hao-Ti Hung, Tung Tsao, Ying-Jer Kao

We introduce a quantum-circuit algorithm for performing imaginary-time evolution on infinite one-dimensional lattice systems. The method uses a parameterized quantum circuit to represent a uniform matrix product state ansatz. We derive the ITE algorithm using the time-dependent variational principle and employ the quantum Lanczos algorithm to improve the ground-state energy estimate. As a benchmark, we simulate the transverse-field Ising model using both classical simulators and IBM Quantum devices. Our analysis includes a statistical study of the distributions of the cost function and energy density obtained from quantum measurements, illustrating the effects of finite-sampling noise on convergence.

quant-pharXiv:2608.30363•Sep 1, 2026•0 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
We give a generalization to an infinite tree geometry of Vidal's infinite time-evolving block decimation (iTEBD) algorithm [G. Vidal, Phys. Rev. Lett. 98, 070201 (2007)] for simulating an infinite line of quantum spins. We numerically investigate the quantum Ising model in a transverse field on the Bethe lattice using the matrix product state ansatz. We observe a second order phase transition, with certain key differences from the transverse field Ising model on an infinite spin chain. We also investigate a transverse field Ising model with a specific longitudinal field. When the transverse field is turned off, this model has a highly degenerate ground state as opposed to the pure Ising model whose ground state is only doubly degenerate.
We show that the simple update approach proposed by Jiang et al. [H. C. Jiang, Z. Y. Weng, and T. Xiang, Phys. Rev. Lett. 101, 090603 (2008)] is an efficient and accurate method for determining the infinite tree tensor network states on the Bethe lattice. Ground-state properties of the quantum transverse Ising model and the Heisenberg XXZ model on the Bethe lattice are studied. The transverse Ising model is found to undergo a second-order quantum phase transition with a diverging magnetic susceptibility but a finite correlation length which is upper bounded by $1/\mathrm{ln}(q\ensuremath{-}1)$ even at the transition point (qqq is the coordinate number of the Bethe lattice). An intuitive explanation on this peculiar ``critical'' phenomenon is given. The XXZ model on the Bethe lattice undergoes a first-order quantum phase transition at the isotropic point. Furthermore, the simple update scheme is found to be related with the Bethe approximation. Finally, by applying the simple update to various tree tensor clusters, we can obtain rather nice and scalable approximations for two-dimensional lattices.
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.
Abstract This study investigates the quantum phase transitions and critical behaviors in the transverse Ising model with staggered Γ interaction by utilizing the infinite-time evolving block decimation method in the infinite matrix product states representation. The ground-state phase diagram shows that the system exhibits a phase transition from the transverse Néel phase to the paramagnetic phase as the transverse magnetic field increases, which are characterized by the local magnetization and non-local transverse Néel order. The chiral order and nematic order are calculated to explore the helical and nematic magnetic structure caused by the Γ interaction. From the spin-spin correlation, a modulated commensurate power-law oscillating decay behavior provides insight in analyzing the ground state characteristics of each phase. In order to study the quantum criticality, the critical exponents β , ν and central charge c are extracted and numerically consistent with the Ising class transition, demonstrating that the critical behavior in our model occurs is independent of the emergent Gamma interaction.
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.
We introduce a quantum-circuit algorithm for performing imaginary-time evolution on infinite one-dimensional lattice systems. The method uses a parameterized quantum circuit to represent a uniform matrix product state ansatz. We derive the ITE algorithm using the time-dependent variational principle and employ the quantum Lanczos algorithm to improve the ground-state energy estimate. As a benchmark, we simulate the transverse-field Ising model using both classical simulators and IBM Quantum devices. Our analysis includes a statistical study of the distributions of the cost function and energy density obtained from quantum measurements, illustrating the effects of finite-sampling noise on convergence.
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​.