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

Ground states and dynamical properties of the S>1/2 quantum Heisenberg model on the 1/5-depleted square lattice

Jun-Han Huang, Zenan Liu, Han-Qing Wu, Dao‐Xin Yao

Published Aug 1, 2022• 5 citations
arXivHTMLPDF

Abstract

We study the $S>1/2$ antiferromagnetic Heisenberg model on the 1/5-depleted square lattice as a function of the ratio of the intraplaquette coupling to the interplaquette coupling. Using stochastic series expansion quantum Monte Carlo simulations, we numerically identify three quantum phases, including the dimer phase, N\'eel phase, and plaquette valence-bond-solid phase. We also obtain the accurate quantum critical points that belong to the $O$(3) universality class using the large-scale finite-size scaling. Most importantly, we study the dynamic spin structure factors of different phases, which can be measured by inelastic neutron scattering experiments. The low-energy excitations can be explained as triplons in the dimer phase and plaquette valence-bond-solid phase, while in the N\'eel phase, the more prominent magnon mode can be found as the spin magnitude increases. Furthermore, we find a broader continuum at smaller $S$, which may be the dynamical signature of nearly deconfined spinon excitations.

DOI
10.1103/physrevb.106.085101

Related papers

Spectral signatures of fractionalization in the frustrated Heisenberg model on the square lattice

Francesco Ferrari, Federico Becca

We employ a variational Monte Carlo approach to efficiently obtain the dynamical structure factor for the spin-$1/2\phantom{\rule{0.16em}{0ex}}{J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model on the square lattice. Upon increasing the frustrating ratio ${J}_{2}/{J}_{1}$, the ground state undergoes a continuous transition from a N\'eel antiferromagnet to a ${\mathbb{Z}}_{2}$ gapless spin liquid. We identify the characteristic spectral features in both phases and highlight the existence of a broad continuum of excitations in the proximity of the spin-liquid phase. The magnon branch, which dominates the spectrum of the unfrustrated Heisenberg model, gradually loses its spectral weight, thus releasing nearly deconfined spinons, whose signatures are visible even in the magnetically ordered state. Our results provide an important example on how magnons fractionalize into deconfined spinons across a quantum critical point.

arXiv:1805.09287•Physical review. B./Physical review. B•Sep 17, 2018•82 citations

Plaquette Ordered Phase and Quantum Phase Diagram in the Spin- 12 J1−J2 Square Heisenberg Model

Shou-Shu Gong, Wei Kang Zhu, D. N. Sheng +2 more

We study the spin-1/2 Heisenberg model on the square lattice with first- and second-neighbor antiferromagnetic interactions J(1) and J(2), which possesses a nonmagnetic region that has been debated for many years and might realize the interesting Z(2) spin liquid. We use the density matrix renormalization group approach with explicit implementation of SU(2) spin rotation symmetry and study the model accurately on open cylinders with different boundary conditions. With increasing J(2), we find a Néel phase and a plaquette valence-bond (PVB) phase with a finite spin gap. From the finite-size scaling of the magnetic order parameter, we estimate that the Néel order vanishes at J(2)/J(1)≃0.44. For 0.5<J(2)/J(1)<0.61, we find dimer correlations and PVB textures whose decay lengths grow strongly with increasing system width, consistent with a long-range PVB order in the two-dimensional limit. The dimer-dimer correlations reveal the s-wave character of the PVB order. For 0.44<J(2)/J(1)<0.5, spin order, dimer order, and spin gap are small on finite-size systems, which is consistent with a near-critical behavior. The critical exponents obtained from the finite-size spin and dimer correlations could be compatible with the deconfined criticality in this small region. We compare and contrast our results with earlier numerical studies.

arXiv:1311.5962•Physical Review Letters•Jul 7, 2014•275 citations

Phase diagram and magnetic excitations of J1–J3 Heisenberg model on the square lattice

Muwei Wu, Shou-Shu Gong, Dao‐Xin Yao +1 more

We study the phase diagram and the dynamical spin structure factor of the spin-1/2 ${J}_{1}--{J}_{3}$ Heisenberg model on the square lattice using density matrix renormalization group, exact diagonalization (ED), and cluster perturbation theory (CPT). By extrapolating the order parameters and studying the level crossings of the low-lying energy and entanglement spectra, we obtain the phase diagram of this model and identify a narrow region of quantum spin liquid (QSL) phase followed by a plaquette valence-bond solid (PVBS) state in the intermediate region, whose nature has been controversial for many years. More importantly, we use CPT and ED to study the dynamical spin structure factor in the QSL and the PVBS phase. In the QSL phase, the high-energy magnon mode completely turns into some dispersive weak excitations around the $X$ and $M$ points. For the PVBS phase, the low-energy spectrum is characterized by a gapped triplet excitation, and at the high energy, we find another branch of dispersive excitation with broad continua, which is unlike the plaquette phase in the $2\ifmmode\times\else\texttimes\fi{}2$ checkerboard model. In the latter case, the second branch of excitation is nearly flat due to the weak effective interactions between the local excitations of the plaquettes. And in the ${J}_{1}--{J}_{3}$ Heisenberg model, the uniform interactions and the spontaneously translational symmetry breaking of the PVBS phase make the difference in the excitation spectra.

arXiv:2205.07910•Physical review. B./Physical review. B•Sep 16, 2022•7 citations

Magnetization Plateaus in the Two-dimensional S = 1/2 Heisenberg Model with a 3$\times$3 Checkerboard Structure

Liang, Xuyang, Dao‐Xin Yao

We investigate the S=1/2 antiferromagnetic Heisenberg model with a 3$\times$3 checkerboard lattice structure in a longitudinal magnetic field. By using the stochastic series expansion quantum Monte Carlo (SSE-QMC) method, we obtain the properties of the non-plateau XY phase, 1/9, 3/9, 5/9, 7/9 magnetization plateau phases, and fully polarized phase. Then, we determine the precise phase transition critical points belonging to the 3D XY universality class through finite-size scaling. Moreover, we study the longitudinal and transverse dynamic spin structure factors of this model in different phases. For the non-plateau XY phase, the energy spectra present a gap between the low-energy gapless branch and the high-energy part under the competition of magnetic field and interaction. The gapless branch can be described by the spin wave theory in the canted antiferromagnetic phase of the effective "block spin" model. In the magnetization plateau phase, we identify that the excitation arises from localized disturbances within the sublattice, which are capable of spreading in momentum space. This study offers theoretical insights and interpretations for the characteristics of the ground state and the inelastic neutron scattering spectrum in two-dimensional quantum magnetic materials. Specifically, it focuses on materials with a checkerboard unit cell structure and odd-spin configurations under the influence of a longitudinal magnetic field.

arXiv:2408.11899•Aug 21, 2024•0 citations

Tensor-product state approach to spin- 12 square J1−J2 antiferromagnetic Heisenberg model: Evidence for deconfined quantum criticality

Ling Wang, Zheng‐Cheng Gu, Frank Verstraete +1 more

The ground state phase of a spin-$\frac{1}{2} {J}_{1}\text{\ensuremath{-}}{J}_{2}$ antiferromagnetic Heisenberg model on a square lattice around the maximally frustrated regime $({J}_{2}\ensuremath{\sim}0.5{J}_{1})$ has been debated for decades. Here we study this model using the cluster update algorithm for tensor-product states (TPSs). The ground state energies at finite sizes and in the thermodynamic limit (with finite size scaling) are in good agreement with exact diagonalization study. Through finite size scaling of the spin correlation function, we find the critical point ${J}_{2}^{{c}_{1}}=0.572(5){J}_{1}$ and critical exponents $\ensuremath{\nu}=0.50(8),\phantom{\rule{0.16em}{0ex}}{\ensuremath{\eta}}_{s}=0.28(6)$. In the range of $0.572<{J}_{2}/{J}_{1}\ensuremath{\le}0.6$ we find a paramagnetic ground state with an exponentially decaying spin-spin correlation. Up to a $24\ifmmode\times\else\texttimes\fi{}24$ system size, we observe power law decaying dimer-dimer and plaquette-plaquette correlations with an anomalous plaquette scaling exponent ${\ensuremath{\eta}}_{p}=0.24(1)$ and an anomalous columnar scaling exponent ${\ensuremath{\eta}}_{c}=0.28(1)$ at ${J}_{2}/{J}_{1}=0.6$. These results are consistent with a potential gapless $U(1)$ spin-liquid phase. However, since the $U(1)$ spin liquid is unstable due to the instanton effect, a valence bond solid order with very small amplitude might develop in the thermodynamic limit. Thus, our numerical results strongly indicate a deconfined quantum critical point at ${J}_{2}^{{c}_{1}}$. Remarkably, all the observed critical exponents are consistent with the $J\text{\ensuremath{-}}Q$ model.

arXiv:1112.3331•Physical review. B./Physical review. B•Aug 19, 2016•102 citations

Magnetization plateaus in the two-dimensional S = 1 2 Heisenberg model with a 3 × 3 checkerboard structure

Xuyang Liang, Dao‐Xin Yao

We investigate the $S=1/2$ antiferromagnetic Heisenberg model with a $3\ifmmode\times\else\texttimes\fi{}3$ checkerboard lattice structure in a longitudinal magnetic field. By using the stochastic series expansion quantum Monte Carlo method, we obtain the properties of the nonplateau XY phase; 1/9, 3/9, 5/9, and 7/9 magnetization plateau phases; and fully polarized phase. Then, we determine the precise phase transition critical points belonging to the three-dimensional XY universality class through finite-size scaling. Moreover, we study the longitudinal and transverse dynamic spin structure factors of this model in different phases. For the nonplateau XY phase, the energy spectra present a gap between the low-energy gapless branch and the high-energy part due to the competition of magnetic field and interaction. The gapless branch can be described by the spin wave theory of the effective ``block spin'' model in the canted antiferromagnetic phase. In the magnetization plateau phase, we identify that the low-energy excitations arise from the localized disturbances within the sublattice and are capable of spreading in momentum space. The high-energy quasiparticles could be the multiplet excitations similar to the ``triplon'' in the dimerized system and ``doublon'' and ``quarton'' in the trimerized system. This paper offers theoretical insights and interpretations for the characteristics of the ground state and the inelastic neutron scattering spectra in two-dimensional quantum magnetic materials. Specifically, it focuses on materials with odd-spin configurations in the checkerboard unit cells under the influence of a longitudinal magnetic field.

Physical review. B./Physical review. B•Mar 19, 2025•3 citations
We study the S>1/2S>1/2S>1/2 antiferromagnetic Heisenberg model on the 1/5-depleted square lattice as a function of the ratio of the intraplaquette coupling to the interplaquette coupling. Using stochastic series expansion quantum Monte Carlo simulations, we numerically identify three quantum phases, including the dimer phase, N'eel phase, and plaquette valence-bond-solid phase. We also obtain the accurate quantum critical points that belong to the OOO(3) universality class using the large-scale finite-size scaling. Most importantly, we study the dynamic spin structure factors of different phases, which can be measured by inelastic neutron scattering experiments. The low-energy excitations can be explained as triplons in the dimer phase and plaquette valence-bond-solid phase, while in the N'eel phase, the more prominent magnon mode can be found as the spin magnitude increases. Furthermore, we find a broader continuum at smaller SSS, which may be the dynamical signature of nearly deconfined spinon excitations.
We employ a variational Monte Carlo approach to efficiently obtain the dynamical structure factor for the spin-$1/2\phantom{\rule{0.16em}{0ex}}{J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model on the square lattice. Upon increasing the frustrating ratio J2/J1{J}_{2}/{J}_{1}J2​/J1​, the ground state undergoes a continuous transition from a N\'eel antiferromagnet to a Z2{\mathbb{Z}}_{2}Z2​ gapless spin liquid. We identify the characteristic spectral features in both phases and highlight the existence of a broad continuum of excitations in the proximity of the spin-liquid phase. The magnon branch, which dominates the spectrum of the unfrustrated Heisenberg model, gradually loses its spectral weight, thus releasing nearly deconfined spinons, whose signatures are visible even in the magnetically ordered state. Our results provide an important example on how magnons fractionalize into deconfined spinons across a quantum critical point.
We study the spin-1/2 Heisenberg model on the square lattice with first- and second-neighbor antiferromagnetic interactions J(1) and J(2), which possesses a nonmagnetic region that has been debated for many years and might realize the interesting Z(2) spin liquid. We use the density matrix renormalization group approach with explicit implementation of SU(2) spin rotation symmetry and study the model accurately on open cylinders with different boundary conditions. With increasing J(2), we find a Néel phase and a plaquette valence-bond (PVB) phase with a finite spin gap. From the finite-size scaling of the magnetic order parameter, we estimate that the Néel order vanishes at J(2)/J(1)≃0.44. For 0.5<J(2)/J(1)<0.61, we find dimer correlations and PVB textures whose decay lengths grow strongly with increasing system width, consistent with a long-range PVB order in the two-dimensional limit. The dimer-dimer correlations reveal the s-wave character of the PVB order. For 0.44<J(2)/J(1)<0.5, spin order, dimer order, and spin gap are small on finite-size systems, which is consistent with a near-critical behavior. The critical exponents obtained from the finite-size spin and dimer correlations could be compatible with the deconfined criticality in this small region. We compare and contrast our results with earlier numerical studies.
We study the phase diagram and the dynamical spin structure factor of the spin-1/2 J1−−J3{J}_{1}--{J}_{3}J1​−−J3​ Heisenberg model on the square lattice using density matrix renormalization group, exact diagonalization (ED), and cluster perturbation theory (CPT). By extrapolating the order parameters and studying the level crossings of the low-lying energy and entanglement spectra, we obtain the phase diagram of this model and identify a narrow region of quantum spin liquid (QSL) phase followed by a plaquette valence-bond solid (PVBS) state in the intermediate region, whose nature has been controversial for many years. More importantly, we use CPT and ED to study the dynamical spin structure factor in the QSL and the PVBS phase. In the QSL phase, the high-energy magnon mode completely turns into some dispersive weak excitations around the XXX and MMM points. For the PVBS phase, the low-energy spectrum is characterized by a gapped triplet excitation, and at the high energy, we find another branch of dispersive excitation with broad continua, which is unlike the plaquette phase in the $2\ifmmode\times\else\texttimes\fi{}2$ checkerboard model. In the latter case, the second branch of excitation is nearly flat due to the weak effective interactions between the local excitations of the plaquettes. And in the J1−−J3{J}_{1}--{J}_{3}J1​−−J3​ Heisenberg model, the uniform interactions and the spontaneously translational symmetry breaking of the PVBS phase make the difference in the excitation spectra.
We investigate the S=1/2 antiferromagnetic Heisenberg model with a 3$\times$3 checkerboard lattice structure in a longitudinal magnetic field. By using the stochastic series expansion quantum Monte Carlo (SSE-QMC) method, we obtain the properties of the non-plateau XY phase, 1/9, 3/9, 5/9, 7/9 magnetization plateau phases, and fully polarized phase. Then, we determine the precise phase transition critical points belonging to the 3D XY universality class through finite-size scaling. Moreover, we study the longitudinal and transverse dynamic spin structure factors of this model in different phases. For the non-plateau XY phase, the energy spectra present a gap between the low-energy gapless branch and the high-energy part under the competition of magnetic field and interaction. The gapless branch can be described by the spin wave theory in the canted antiferromagnetic phase of the effective "block spin" model. In the magnetization plateau phase, we identify that the excitation arises from localized disturbances within the sublattice, which are capable of spreading in momentum space. This study offers theoretical insights and interpretations for the characteristics of the ground state and the inelastic neutron scattering spectrum in two-dimensional quantum magnetic materials. Specifically, it focuses on materials with a checkerboard unit cell structure and odd-spin configurations under the influence of a longitudinal magnetic field.
The ground state phase of a spin-$\frac{1}{2} {J}_{1}\text{\ensuremath{-}}{J}_{2}$ antiferromagnetic Heisenberg model on a square lattice around the maximally frustrated regime $({J}_{2}\ensuremath{\sim}0.5{J}_{1})$ has been debated for decades. Here we study this model using the cluster update algorithm for tensor-product states (TPSs). The ground state energies at finite sizes and in the thermodynamic limit (with finite size scaling) are in good agreement with exact diagonalization study. Through finite size scaling of the spin correlation function, we find the critical point J2c1=0.572(5)J1{J}_{2}^{{c}_{1}}=0.572(5){J}_{1}J2c1​​=0.572(5)J1​ and critical exponents $\ensuremath{\nu}=0.50(8),\phantom{\rule{0.16em}{0ex}}{\ensuremath{\eta}}_{s}=0.28(6)$. In the range of $0.572<{J}_{2}/{J}_{1}\ensuremath{\le}0.6$ we find a paramagnetic ground state with an exponentially decaying spin-spin correlation. Up to a $24\ifmmode\times\else\texttimes\fi{}24$ system size, we observe power law decaying dimer-dimer and plaquette-plaquette correlations with an anomalous plaquette scaling exponent ${\ensuremath{\eta}}_{p}=0.24(1)$ and an anomalous columnar scaling exponent ${\ensuremath{\eta}}_{c}=0.28(1)$ at J2/J1=0.6{J}_{2}/{J}_{1}=0.6J2​/J1​=0.6. These results are consistent with a potential gapless U(1)U(1)U(1) spin-liquid phase. However, since the U(1)U(1)U(1) spin liquid is unstable due to the instanton effect, a valence bond solid order with very small amplitude might develop in the thermodynamic limit. Thus, our numerical results strongly indicate a deconfined quantum critical point at J2c1{J}_{2}^{{c}_{1}}J2c1​​. Remarkably, all the observed critical exponents are consistent with the $J\text{\ensuremath{-}}Q$ model.
We investigate the S=1/2S=1/2S=1/2 antiferromagnetic Heisenberg model with a $3\ifmmode\times\else\texttimes\fi{}3$ checkerboard lattice structure in a longitudinal magnetic field. By using the stochastic series expansion quantum Monte Carlo method, we obtain the properties of the nonplateau XY phase; 1/9, 3/9, 5/9, and 7/9 magnetization plateau phases; and fully polarized phase. Then, we determine the precise phase transition critical points belonging to the three-dimensional XY universality class through finite-size scaling. Moreover, we study the longitudinal and transverse dynamic spin structure factors of this model in different phases. For the nonplateau XY phase, the energy spectra present a gap between the low-energy gapless branch and the high-energy part due to the competition of magnetic field and interaction. The gapless branch can be described by the spin wave theory of the effective ``block spin'' model in the canted antiferromagnetic phase. In the magnetization plateau phase, we identify that the low-energy excitations arise from the localized disturbances within the sublattice and are capable of spreading in momentum space. The high-energy quasiparticles could be the multiplet excitations similar to the ``triplon'' in the dimerized system and ``doublon'' and ``quarton'' in the trimerized system. This paper offers theoretical insights and interpretations for the characteristics of the ground state and the inelastic neutron scattering spectra in two-dimensional quantum magnetic materials. Specifically, it focuses on materials with odd-spin configurations in the checkerboard unit cells under the influence of a longitudinal magnetic field.