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:2210.06764 · Physical review. B./Physical review. B

Classical model emerges in quantum entanglement: Quantum Monte Carlo study for an Ising-Heisenberg bilayer

Siying Wu, Xiaoxue Ran, Binbin Yin, Qifang Li, Bin-Bin Mao, Yan-Cheng Wang, Yan Qing Zheng

Published Apr 10, 2023• 15 citations
arXivHTMLPDF

Abstract

By developing a cluster sampling of the stochastic series expansion quantum Monte Carlo method, we investigate a spin-$\frac{1}{2}$ model on a bilayer square lattice with intralayer ferromagnetic (FM) Ising coupling and interlayer antiferromagnetic Heisenberg interaction. The continuous quantum phase transition which occurs at ${g}_{c}=3.045(2)$ between the FM Ising phase and the dimerized phase is studied via large-scale simulations. From analysis of the critical exponents we show that this phase transition belongs to the ($2+1$)-dimensional Ising universality class. In addition, the quantum entanglement is strong between the two layers, especially in the dimerized phase. The effective Hamiltonian of a single layer seems like a transverse-field Ising model. However, we found that the quantum entanglement Hamiltonian is a pure classical Ising model without any quantum fluctuations. Furthermore, we give a more general explanation about how a classical entanglement Hamiltonian emerges.

DOI
10.1103/physrevb.107.155121

Related papers

Monte Carlo simulation of quantum Potts model

Chengxiang Ding, Yangcheng Wang, Youjin Deng +1 more

Using Monte Carlo simulations in the frame of stochastic series expansion (SSE), we study the three-state quantum Potts model. The cluster algorithm we used is a direct generalization of that for the quantum Ising model. The simulations include the one dimensional and two dimensional ferromagnetic three-state quantum Potts model and the two dimensional antiferromagnetic three-state quantum Potts model. Our results show that the phase transition of the one dimensional ferromagnetic quantum Potts model belongs to the same universality class of the two dimensional classical Potts model, the two dimensional ferromagnetic quantum Potts model undergoes a first order transition, which is also in analogy to its classical correspondence. The phase transition of the antiferromagnetic quantum Potts model is continuous, whose universality class belongs to the three-dimensional classical XY model, owing to an `emergent' O(2) symmetry at the critical point, although its ordered phase breaks the Z_6 symmetry.

arXiv:1702.02675•Feb 9, 2017•3 citations

Mott glass phase in a diluted bilayer Heisenberg quantum antiferromagnet

Nvsen Ma, Anders W. Sandvik, Dao‐Xin Yao

We use quantum Monte Carlo simulations to study a dimer-diluted S = 1/2 Heisenberg model on a bilayer square lattice with intralayer interaction J 1 and interlayer interaction J 2 . Below the classical percolation threshold p c , the system has three phases reachable by tuning the interaction ratio g = J 2 / J 1 : a Néel ordered phase, a gapless quantum glass phase, and a gapped quantum paramagnetic phase. We present the ground-state phase diagram in the plane of dilution p and interaction ratio g . The quantum glass phase is certified to be of the gapless Mott glass type, having a uniform susceptibility vanishing at zero temperature T and following a stretched exponential form at T > 0; χ u exp(- b / T α ) with α < 1. At the phase transition point from Neel ordered to Mott glass, we find that the critical exponents are different from those of the clean system described by the standard O (3) universality class in 2+1 dimensions.

arXiv:1511.07895•Journal of Physics Conference Series•Sep 28, 2015•1 citation

Ising-type critical exponents of the fully frustrated spin-1/2 Heisenberg FM/AF square bilayer at a critical magnetic field

Katarína Karľová, Jozef Strečka

The fully frustrated spin-1/2 Heisenberg FM/AF square bilayer in a magnetic field with the ferromagnetic inter-dimer interaction and the antiferromagnetic intra-dimer interaction is explored by the use of localized many-magnon approach, which allows to connect the original purely quantum Heisenberg spin model on a square bilayer with the effective ferromagnetic Ising model on a simple square lattice. Magnetization and specific heat are investigated exactly at a field-driven phase transition from the singlet-dimer phase towards the fully saturated ferromagnetic phase, which changes from a discontinuous phase transition to a continuous one at a certain critical temperature. The mapping correspondence between the spin-1/2 Heisenberg FM/AF square bilayer and the ferromagnetic Ising square lattice suggests for this special critical point of the spin-1/2 Heisenberg FM/AF square bilayer critical exponents from the standard two-dimensional Ising universality class.

Phase Transitions•Feb 7, 2019•1 citation

Mott glass and criticality in an S = 1 2 bilayer Heisenberg model with interlayer bond dilution

K. J. Li, Han-Qing Wu, Dao‐Xin Yao

We employ the stochastic series expansion quantum Monte Carlo method to investigate the $S=1/2$ antiferromagnetic Heisenberg model on a bilayer square lattice with diluted interlayer couplings. Both regular and random dilution patterns are considered. In systems with regular dilution, tuning the interlayer interaction drives a quantum phase transition from a N\'eel-ordered phase to a quantum disordered phase, consistent with the $O(3)$ universality class. In contrast, random dilution gives rise to a two-step transition: from the N\'eel phase to an intermediate Mott glass (MG) phase, followed by a transition to the quantum disordered phase. Within the MG phase, the uniform magnetic susceptibility exhibits a stretched-exponential temperature dependence ${\ensuremath{\chi}}_{u}\ensuremath{\sim}exp(\ensuremath{-}b/{T}^{\ensuremath{\alpha}}), 0<\ensuremath{\alpha}<1$. At the N\'eel-to-glass transition, the quenched disorder modifies the critical exponents in a manner consistent with the Harris criterion. These findings provide new insights into disorder-driven quantum phase transitions and the emergence of glassy phases in diluted bilayer quantum magnets.

Physical review. B./Physical review. B•Nov 8, 2025•0 citations

Quantum-critical properties of the long-range transverse-field Ising model from quantum Monte Carlo simulations

Jan Alexander Koziol, Anja Langheld, Sebastian C. Kapfer +1 more

The quantum-critical properties of the transverse-field Ising model with algebraically decaying interactions are investigated by means of stochastic series expansion quantum Monte Carlo, on both the one-dimensional linear chain and the two-dimensional square lattice. We extract the critical exponents $\ensuremath{\nu}$ and $\ensuremath{\beta}$ as a function of the decay exponent of the long-range interactions. For ferromagnetic Ising interactions, we resolve the limiting regimes known from field theory, ranging from the nearest-neighbor Ising to the long-range Gaussian universality classes, as well as the intermediate regime with continuously varying critical exponents. In the long-range Gaussian regime, we treat the effect of dangerous irrelevant variables on finite-size scaling forms. For antiferromagnetic and therefore competing Ising interactions, the stochastic series expansion algorithm displays growing autocorrelation times leading to a reduced performance. Nevertheless, our results are consistent with the nearest-neighbor Ising universality for all investigated interaction ranges both on the linear chain and the square lattice.

arXiv:2103.09469•Physical review. B./Physical review. B•Jun 23, 2021•47 citations

Evidence for an Unconventional Universality Class from a Two-Dimensional Dimerized Quantum Heisenberg Model

S. C. Wenzel, L. Bogacz, Wolfhard Janke E

The two-dimensional $J\mathrm{\text{\ensuremath{-}}}{J}^{\ensuremath{'}}$ dimerized quantum Heisenberg model is studied on the square lattice by means of (stochastic series expansion) quantum Monte Carlo simulations as a function of the coupling ratio $\ensuremath{\alpha}={J}^{\ensuremath{'}}/J$. The critical point of the order-disorder quantum phase transition in the $J\mathrm{\text{\ensuremath{-}}}{J}^{\ensuremath{'}}$ model is determined as ${\ensuremath{\alpha}}_{c}=2.5196(2)$ by finite-size scaling for up to approximately 10 000 quantum spins. By comparing six dimerized models we show, contrary to the current belief, that the critical exponents of the $J\mathrm{\text{\ensuremath{-}}}{J}^{\ensuremath{'}}$ model are not in agreement with the three-dimensional classical Heisenberg universality class. This lends support to the notion of nontrivial critical excitations at the quantum critical point.

arXiv:0805.2500•Physical Review Letters•Sep 19, 2008•88 citations
By developing a cluster sampling of the stochastic series expansion quantum Monte Carlo method, we investigate a spin-12\frac{1}{2}21​ model on a bilayer square lattice with intralayer ferromagnetic (FM) Ising coupling and interlayer antiferromagnetic Heisenberg interaction. The continuous quantum phase transition which occurs at gc=3.045(2){g}_{c}=3.045(2)gc​=3.045(2) between the FM Ising phase and the dimerized phase is studied via large-scale simulations. From analysis of the critical exponents we show that this phase transition belongs to the (2+12+12+1)-dimensional Ising universality class. In addition, the quantum entanglement is strong between the two layers, especially in the dimerized phase. The effective Hamiltonian of a single layer seems like a transverse-field Ising model. However, we found that the quantum entanglement Hamiltonian is a pure classical Ising model without any quantum fluctuations. Furthermore, we give a more general explanation about how a classical entanglement Hamiltonian emerges.
Using Monte Carlo simulations in the frame of stochastic series expansion (SSE), we study the three-state quantum Potts model. The cluster algorithm we used is a direct generalization of that for the quantum Ising model. The simulations include the one dimensional and two dimensional ferromagnetic three-state quantum Potts model and the two dimensional antiferromagnetic three-state quantum Potts model. Our results show that the phase transition of the one dimensional ferromagnetic quantum Potts model belongs to the same universality class of the two dimensional classical Potts model, the two dimensional ferromagnetic quantum Potts model undergoes a first order transition, which is also in analogy to its classical correspondence. The phase transition of the antiferromagnetic quantum Potts model is continuous, whose universality class belongs to the three-dimensional classical XY model, owing to an `emergent' O(2) symmetry at the critical point, although its ordered phase breaks the Z_6 symmetry.
We use quantum Monte Carlo simulations to study a dimer-diluted S = 1/2 Heisenberg model on a bilayer square lattice with intralayer interaction J 1 and interlayer interaction J 2 . Below the classical percolation threshold p c , the system has three phases reachable by tuning the interaction ratio g = J 2 / J 1 : a Néel ordered phase, a gapless quantum glass phase, and a gapped quantum paramagnetic phase. We present the ground-state phase diagram in the plane of dilution p and interaction ratio g . The quantum glass phase is certified to be of the gapless Mott glass type, having a uniform susceptibility vanishing at zero temperature T and following a stretched exponential form at T > 0; χ u exp(- b / T α ) with α < 1. At the phase transition point from Neel ordered to Mott glass, we find that the critical exponents are different from those of the clean system described by the standard O (3) universality class in 2+1 dimensions.
The fully frustrated spin-1/2 Heisenberg FM/AF square bilayer in a magnetic field with the ferromagnetic inter-dimer interaction and the antiferromagnetic intra-dimer interaction is explored by the use of localized many-magnon approach, which allows to connect the original purely quantum Heisenberg spin model on a square bilayer with the effective ferromagnetic Ising model on a simple square lattice. Magnetization and specific heat are investigated exactly at a field-driven phase transition from the singlet-dimer phase towards the fully saturated ferromagnetic phase, which changes from a discontinuous phase transition to a continuous one at a certain critical temperature. The mapping correspondence between the spin-1/2 Heisenberg FM/AF square bilayer and the ferromagnetic Ising square lattice suggests for this special critical point of the spin-1/2 Heisenberg FM/AF square bilayer critical exponents from the standard two-dimensional Ising universality class.
We employ the stochastic series expansion quantum Monte Carlo method to investigate the S=1/2S=1/2S=1/2 antiferromagnetic Heisenberg model on a bilayer square lattice with diluted interlayer couplings. Both regular and random dilution patterns are considered. In systems with regular dilution, tuning the interlayer interaction drives a quantum phase transition from a N\'eel-ordered phase to a quantum disordered phase, consistent with the O(3)O(3)O(3) universality class. In contrast, random dilution gives rise to a two-step transition: from the N\'eel phase to an intermediate Mott glass (MG) phase, followed by a transition to the quantum disordered phase. Within the MG phase, the uniform magnetic susceptibility exhibits a stretched-exponential temperature dependence ${\ensuremath{\chi}}_{u}\ensuremath{\sim}exp(\ensuremath{-}b/{T}^{\ensuremath{\alpha}}), 0<\ensuremath{\alpha}<1$. At the N\'eel-to-glass transition, the quenched disorder modifies the critical exponents in a manner consistent with the Harris criterion. These findings provide new insights into disorder-driven quantum phase transitions and the emergence of glassy phases in diluted bilayer quantum magnets.
The quantum-critical properties of the transverse-field Ising model with algebraically decaying interactions are investigated by means of stochastic series expansion quantum Monte Carlo, on both the one-dimensional linear chain and the two-dimensional square lattice. We extract the critical exponents $\ensuremath{\nu}$ and $\ensuremath{\beta}$ as a function of the decay exponent of the long-range interactions. For ferromagnetic Ising interactions, we resolve the limiting regimes known from field theory, ranging from the nearest-neighbor Ising to the long-range Gaussian universality classes, as well as the intermediate regime with continuously varying critical exponents. In the long-range Gaussian regime, we treat the effect of dangerous irrelevant variables on finite-size scaling forms. For antiferromagnetic and therefore competing Ising interactions, the stochastic series expansion algorithm displays growing autocorrelation times leading to a reduced performance. Nevertheless, our results are consistent with the nearest-neighbor Ising universality for all investigated interaction ranges both on the linear chain and the square lattice.
The two-dimensional $J\mathrm{\text{\ensuremath{-}}}{J}^{\ensuremath{'}}$ dimerized quantum Heisenberg model is studied on the square lattice by means of (stochastic series expansion) quantum Monte Carlo simulations as a function of the coupling ratio $\ensuremath{\alpha}={J}^{\ensuremath{'}}/J$. The critical point of the order-disorder quantum phase transition in the $J\mathrm{\text{\ensuremath{-}}}{J}^{\ensuremath{'}}$ model is determined as ${\ensuremath{\alpha}}_{c}=2.5196(2)$ by finite-size scaling for up to approximately 10 000 quantum spins. By comparing six dimerized models we show, contrary to the current belief, that the critical exponents of the $J\mathrm{\text{\ensuremath{-}}}{J}^{\ensuremath{'}}$ model are not in agreement with the three-dimensional classical Heisenberg universality class. This lends support to the notion of nontrivial critical excitations at the quantum critical point.