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:1511.07895 · Journal of Physics Conference Series

Mott glass phase in a diluted bilayer Heisenberg quantum antiferromagnet

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

Published Sep 28, 2015• 1 citation
arXivHTMLPDF

Abstract

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.

DOI
10.1088/1742-6596/640/1/012045

Related papers

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

Criticality and Mott glass phase in a disordered two-dimensional quantum spin system

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

We use quantum Monte Carlo simulations to study a disordered $S=1/2$ Heisenberg quantum spin model with three different nearest-neighbor interactions, ${J}_{1}\ensuremath{\le}{J}_{2}\ensuremath{\le}{J}_{3}$, on the square lattice. We consider the regime in which ${J}_{1}$ represents weak bonds, and ${J}_{2}$ and ${J}_{3}$ correspond to two kinds of stronger bonds (dimers) which are randomly distributed on columns forming coupled two-leg ladders. When increasing the average intradimer coupling $({J}_{2}+{J}_{3})/2$, the system undergoes a N\'eel to quantum glass transition of the ground state and later a second transition into a quantum paramagnet. The quantum glass phase is of the gapless Mott glass type (i.e., in boson language it is incompressible at temperature $T=0$), and we find that the temperature dependence of the uniform magnetic susceptibility follows the stretched exponential form $\ensuremath{\chi}\ensuremath{\sim}exp(\ensuremath{-}b/{T}^{\ensuremath{\alpha}})$, with $0<\ensuremath{\alpha}<1$. At the N\'eel-glass transition, we observe the standard O(3) critical exponents, which implies that the Harris criterion for the relevance of the disorder is violated in this system.

arXiv:1407.1113•Physical Review B•Sep 26, 2014•13 citations

Quantum Criticality and Percolation in Dimer-Diluted Two-Dimensional Antiferromagnets

Anders W. Sandvik

The S = 1/2 Heisenberg model is considered on bilayer and single-layer square lattices with couplings J1, J2, with each spin belonging to one J2-coupled dimer. A transition from a Néel to disordered ground state occurs at a critical value of g = J2/J1. The systems are here studied at their dimer-dilution percolation points p*. The multicritical point (g*,p*) previously found for the bilayer is not reproduced for the single layer. Instead, there is a line of critical points (g < g*, p*) with continuously varying exponents. The uniform magnetic susceptibility diverges as T(-alpha) with alpha element of [1/2,1]. This unusual behavior is attributed to an effective free-moment density approximately T(1-alpha). The susceptibility of the bilayer is not divergent but exhibits remarkably robust quantum-critical scaling.

arXiv:cond-mat/0603179•Physical Review Letters•May 22, 2006•31 citations

Quantum critical behavior of the superfluid-Mott glass transition

Thomas Vojta, Jack Crewse, Martin Puschmann +2 more

We investigate the zero-temperature superfluid to insulator transitions in a diluted two-dimensional quantum rotor model with particle-hole symmetry. We map the Hamiltonian onto a classical $(2+1)$-dimensional $XY$ model with columnar disorder which we analyze by means of large-scale Monte Carlo simulations. For dilutions below the lattice percolation threshold, the system undergoes a generic superfluid-Mott glass transition. In contrast to other quantum phase transitions in disordered systems, its critical behavior is of conventional power-law type with universal (dilution-independent) critical exponents $z=1.52(3), \ensuremath{\nu}=1.16(5), \ensuremath{\beta}/\ensuremath{\nu}=0.48(2), \ensuremath{\gamma}/\ensuremath{\nu}=2.52(4)$, and $\ensuremath{\eta}=\ensuremath{-}0.52(4)$. These values agree with and improve upon earlier Monte Carlo results [Phys. Rev. Lett. 92, 015703 (2004)] while (partially) excluding other findings in the literature. As a further test of universality, we also consider a soft-spin version of the classical Hamiltonian. In addition, we study the percolation quantum phase transition across the lattice percolation threshold; its critical behavior is governed by the lattice percolation exponents in agreement with recent theoretical predictions. We relate our results to a general classification of phase transitions in disordered systems, and we briefly discuss experiments.

arXiv:1607.01860•Physical review. B./Physical review. B•Oct 3, 2016•28 citations

Quantum critical behavior of a three-dimensional superfluid-Mott glass transition

Jack Crewse, Cameron Lerch, Thomas Vojta

The superfluid to insulator quantum phase transition of a three-dimensional particle-hole symmetric system of disordered bosons is studied. To this end, a site-diluted quantum rotor Hamiltonian is mapped onto a classical (3+1)-dimensional $XY$ model with columnar disorder and analyzed by means of large-scale Monte Carlo simulations. The superfluid-Mott insulator transition of the clean, undiluted system is in the four-dimensional $XY$ universality class and shows mean-field critical behavior with logarithmic corrections. The clean correlation length exponent $\ensuremath{\nu}=1/2$ violates the Harris criterion, indicating that disorder must be a relevant perturbation. For nonzero dilutions below the lattice percolation threshold of ${p}_{c}=0.688\phantom{\rule{0.16em}{0ex}}392$, our simulations yield conventional power-law critical behavior with dilution-independent critical exponents $z=1.67(6)$, $\ensuremath{\nu}=0.90(5)$, $\ensuremath{\beta}/\ensuremath{\nu}=1.09(3)$, and $\ensuremath{\gamma}/\ensuremath{\nu}=2.50(3)$. The critical behavior of the transition across the lattice percolation threshold is controlled by the classical percolation exponents. Our results are discussed in the context of a classification of disordered quantum phase transitions, as well as experiments in superfluids, superconductors, and magnetic systems.

arXiv:1805.07466•Physical review. B./Physical review. B•Aug 22, 2018•5 citations

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

Siying Wu, Xiaoxue Ran, Binbin Yin +4 more

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.

arXiv:2210.06764•Physical review. B./Physical review. B•Apr 10, 2023•15 citations
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.
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.
We use quantum Monte Carlo simulations to study a disordered S=1/2S=1/2S=1/2 Heisenberg quantum spin model with three different nearest-neighbor interactions, ${J}_{1}\ensuremath{\le}{J}_{2}\ensuremath{\le}{J}_{3}$, on the square lattice. We consider the regime in which J1{J}_{1}J1​ represents weak bonds, and J2{J}_{2}J2​ and J3{J}_{3}J3​ correspond to two kinds of stronger bonds (dimers) which are randomly distributed on columns forming coupled two-leg ladders. When increasing the average intradimer coupling (J2+J3)/2({J}_{2}+{J}_{3})/2(J2​+J3​)/2, the system undergoes a N\'eel to quantum glass transition of the ground state and later a second transition into a quantum paramagnet. The quantum glass phase is of the gapless Mott glass type (i.e., in boson language it is incompressible at temperature T=0T=0T=0), and we find that the temperature dependence of the uniform magnetic susceptibility follows the stretched exponential form $\ensuremath{\chi}\ensuremath{\sim}exp(\ensuremath{-}b/{T}^{\ensuremath{\alpha}})$, with $0<\ensuremath{\alpha}<1$. At the N\'eel-glass transition, we observe the standard O(3) critical exponents, which implies that the Harris criterion for the relevance of the disorder is violated in this system.
The S = 1/2 Heisenberg model is considered on bilayer and single-layer square lattices with couplings J1, J2, with each spin belonging to one J2-coupled dimer. A transition from a Néel to disordered ground state occurs at a critical value of g = J2/J1. The systems are here studied at their dimer-dilution percolation points p*. The multicritical point (g*,p*) previously found for the bilayer is not reproduced for the single layer. Instead, there is a line of critical points (g < g*, p*) with continuously varying exponents. The uniform magnetic susceptibility diverges as T(-alpha) with alpha element of [1/2,1]. This unusual behavior is attributed to an effective free-moment density approximately T(1-alpha). The susceptibility of the bilayer is not divergent but exhibits remarkably robust quantum-critical scaling.
We investigate the zero-temperature superfluid to insulator transitions in a diluted two-dimensional quantum rotor model with particle-hole symmetry. We map the Hamiltonian onto a classical (2+1)(2+1)(2+1)-dimensional XYXYXY model with columnar disorder which we analyze by means of large-scale Monte Carlo simulations. For dilutions below the lattice percolation threshold, the system undergoes a generic superfluid-Mott glass transition. In contrast to other quantum phase transitions in disordered systems, its critical behavior is of conventional power-law type with universal (dilution-independent) critical exponents $z=1.52(3), \ensuremath{\nu}=1.16(5), \ensuremath{\beta}/\ensuremath{\nu}=0.48(2), \ensuremath{\gamma}/\ensuremath{\nu}=2.52(4)$, and $\ensuremath{\eta}=\ensuremath{-}0.52(4)$. These values agree with and improve upon earlier Monte Carlo results [Phys. Rev. Lett. 92, 015703 (2004)] while (partially) excluding other findings in the literature. As a further test of universality, we also consider a soft-spin version of the classical Hamiltonian. In addition, we study the percolation quantum phase transition across the lattice percolation threshold; its critical behavior is governed by the lattice percolation exponents in agreement with recent theoretical predictions. We relate our results to a general classification of phase transitions in disordered systems, and we briefly discuss experiments.
The superfluid to insulator quantum phase transition of a three-dimensional particle-hole symmetric system of disordered bosons is studied. To this end, a site-diluted quantum rotor Hamiltonian is mapped onto a classical (3+1)-dimensional XYXYXY model with columnar disorder and analyzed by means of large-scale Monte Carlo simulations. The superfluid-Mott insulator transition of the clean, undiluted system is in the four-dimensional XYXYXY universality class and shows mean-field critical behavior with logarithmic corrections. The clean correlation length exponent $\ensuremath{\nu}=1/2$ violates the Harris criterion, indicating that disorder must be a relevant perturbation. For nonzero dilutions below the lattice percolation threshold of pc=0.688392{p}_{c}=0.688\phantom{\rule{0.16em}{0ex}}392pc​=0.688392, our simulations yield conventional power-law critical behavior with dilution-independent critical exponents z=1.67(6)z=1.67(6)z=1.67(6), $\ensuremath{\nu}=0.90(5)$, $\ensuremath{\beta}/\ensuremath{\nu}=1.09(3)$, and $\ensuremath{\gamma}/\ensuremath{\nu}=2.50(3)$. The critical behavior of the transition across the lattice percolation threshold is controlled by the classical percolation exponents. Our results are discussed in the context of a classification of disordered quantum phase transitions, as well as experiments in superfluids, superconductors, and magnetic systems.
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.