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:1102.1819 · Journal of Physics Condensed Matter

Variational study of the quantum phase transition in the bilayer Heisenberg model with bosonic RVB wavefunction

Haijun Liao, Tao Li

Published Nov 10, 2011• 6 citations
arXivHTMLPDF

Abstract

We study the ground state phase diagram of the bilayer Heisenberg model on a square lattice with a bosonic resonating valence bond (RVB) wavefunction. The wavefunction has the form of a Gutzwiller projected Schwinger boson mean-field ground state and involves two variational parameters. We find the wavefunction provides an accurate description of the system on both sides of the quantum phase transition. In particular, through the analysis of the spin structure factor, ground state fidelity susceptibility and the Binder moment ratio Q(2), a continuous transition from the antiferromagnetic ordered state to the quantum disordered state is found at the critical coupling of α(c) = J(⊥)/J(∥) ≈ 2.62, in good agreement with the result of quantum Monte Carlo simulation. The critical exponent estimated from the finite size scaling analysis (1/ν ≈ 1.4) is consistent with that of the classical 3D Heisenberg universality class.

DOI
10.1088/0953-8984/23/47/475602

Related papers

Gapped spin-liquid phase in the J1−J2 Heisenberg model by a bosonic resonating valence-bond ansatz

Tao Li, Federico Becca, Wen-Jun Hu +1 more

We study the ground-state phase diagram of the spin-$1/2$${J}_{1}\ensuremath{-}{J}_{2}$ Heisenberg model on the square lattice with an accurate bosonic resonating-valence-bond (RVB) wave function. In contrast to the RVB ansatz based on Schwinger fermions, the representation based on Schwinger bosons, supplemented by a variational Monte Carlo technique enforcing the exact projection onto the physical subspace, is able to describe a fully gapped spin liquid in the strongly frustrated regime. In particular, a fully symmetric ${Z}_{2}$ spin liquid is stable between two antiferromagnetic phases; a continuous transition at ${J}_{2}=0.4{J}_{1}$, when the Marshall sign rule begins to be essentially violated, and a first-order transition around ${J}_{2}=0.6{J}_{1}$ are present. Most importantly, the triplet gap is found to have a nonmonotonic behavior, reaching a maximum around ${J}_{2}=0.51{J}_{1}$, when the lowest spinon excitation moves from the $\ensuremath{\Gamma}$ to the $M$ point, i.e., $k=(\ensuremath{\pi},0)$.

arXiv:1205.3838•Physical Review B•Aug 7, 2012•62 citations

SU(N) quantum spin models: a variational wavefunction study

Arun Paramekanti, John Bradley Marston

SU ( N ) quantum spin systems may be realized in a variety of physical systems including ultracold atoms in optical lattices. The study of such models also leads to insights into possible novel quantum phases and phase transitions of SU (2) spin models. Here we use Gutzwiller projected fermionic variational wavefunctions to explore the phase diagram and correlation functions of SU ( N ) quantum spin models in the self-conjugate representation. In one dimension, the ground state of the SU (4) spin chain with Heisenberg bilinear and biquadratic interactions is studied by examining instabilities of the Gutzwiller projected free fermion ground state to various broken symmetries. The variational phase diagram so obtained agrees well with exact results. The spin–spin and dimer–dimer correlation functions of the Gutzwiller projected free fermion state with N flavours of fermions are in good agreement with exact and 1/ N calculations for the critical points of SU ( N ) spin chains. In two dimensions, the phase diagram of the antiferromagnetic Heisenberg model on the square lattice is obtained by finding instabilities of the Gutzwiller projected π-flux state. In the absence of biquadratic interactions, the model exhibits long-range Néel order for N = 2 and 4, and spin Peierls (columnar dimer) order for N >4. Upon including biquadratic interactions in the SU (4) model (with a sign appropriate to a fermionic Hubbard model), the Néel order diminishes and eventually disappears, giving way to an extended valence bond crystal. In the case of the SU (6) model, the dimerized ground state melts at sufficiently large biquadratic interaction, yielding a projected π-flux spin liquid phase which in turn undergoes a transition into an extended valence bond crystal at even larger biquadratic interaction. The spin correlations of the projected π-flux state at N = 4 are in good agreement with 1/ N calculations. We find that the state shows strongly enhanced dimer correlations, in qualitative agreement with recent theoretical predictions. We also compare our results with a recent quantum Monte Carlo study of the SU (4) Heisenberg model.

arXiv:cond-mat/0608691•Journal of Physics Condensed Matter•Mar 6, 2007•54 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

High-precision finite-size scaling analysis of the quantum-critical point of S=1∕2 Heisenberg antiferromagnetic bilayers

Ling Wang, K. S. D. Beach, Anders W. Sandvik

We use quantum Monte Carlo (stochastic series expansion) and finite-size scaling to study the quantum critical points of two $S=1∕2$ Heisenberg antiferromagnets in two dimensions: a bilayer and a Kondo-lattice-like system (incomplete bilayer), each with intraplane and interplane couplings $J$ and ${J}_{\ensuremath{\perp}}$. We discuss the ground-state finite-size scaling properties of three different quantities---the Binder moment ratio, the spin stiffness, and the long-wavelength magnetic susceptibility---which we use to extract the critical value of the coupling ratio $g={J}_{\ensuremath{\perp}}∕J$. The individual estimates of ${g}_{c}$ are consistent provided that subleading finite-size corrections are properly taken into account. For both models, we find that the spin stiffness has the smallest subleading finite-size corrections; in the case of the incomplete bilayer we find that the first subleading correction vanishes or is extremely small. In agreement with predictions, we find that at the critical point the Binder ratio has a universal value and the product of the spin stiffness and the long-wavelength susceptibility scales as $1∕{L}^{2}$ with a universal prefactor. Our results for the critical coupling ratios are ${g}_{c}=2.5220(1)$ (full bilayer) and ${g}_{c}=1.3888(1)$ (incomplete bilayer), which represent improvements of more than an order of magnitude over the previous best estimates. For the correlation length exponent we obtain $\ensuremath{\nu}=0.7106(9)$, consistent with the expected three-dimensional Heisenberg universality class.

arXiv:cond-mat/0509747•Physical Review B•Jan 27, 2006•175 citations

Quantum phases in the frustrated Heisenberg model on the bilayer honeycomb lattice

Hao Zhang, Marcelo Arlego, Carlos Alberto Lamas

We use a combination of analytical and numerical techniques to study the phase diagram of the frustrated Heisenberg model on the bilayer honeycomb lattice. Using the Schwinger-boson description of the spin operators followed by a mean-field decoupling, the magnetic phase diagram is studied as a function of the frustration coupling ${J}_{2}$ and the interlayer coupling ${J}_{\ensuremath{\perp}}$. The presence of both magnetically ordered and disordered phases is investigated by means of the evaluation of ground-state energy, spin gap, local magnetization, and spin-spin correlations. We observe a phase with a spin gap and short-range N\'eel correlations that survives for nonzero next-nearest-neighbor interaction and interlayer coupling. Furthermore, we detect signatures of a reentrant behavior in the melting of the N\'eel phase and symmetry restoring when the system undergoes a transition from an on-layer nematic valence-bond crystal phase to an interlayer valence-bond crystal phase. We complement our work with exact diagonalization on small clusters and dimer-series expansion calculations, together with a linear spin-wave approach to study the phase diagram as a function of the spin $S$, the frustration, and the interlayer couplings.

arXiv:1312.6324•Physical Review B•Jan 6, 2014•31 citations

Quantum Monte Carlo simulation of S=1/2 Heisenberg model with four spin interaction

M. Tsukamoto, Kenji Harada, Naoki Kawashima

The spin S = 1/2 Heisenberg model with four-spin interaction on the square lattice is studied by using quantum Monte Carlo method. When the four-spin interaction is dominant, the system has a VBS ground state. In this case, we find a finite-temperature second-order phase transition to the VBS state. The universality class of the transition is investigated. We estimate the critical exponents ν and η from the finite size scaling analysis and find ν = 0.68(1) and η = 0.55(2).

Journal of Physics Conference Series•Mar 1, 2009•3 citations
We study the ground state phase diagram of the bilayer Heisenberg model on a square lattice with a bosonic resonating valence bond (RVB) wavefunction. The wavefunction has the form of a Gutzwiller projected Schwinger boson mean-field ground state and involves two variational parameters. We find the wavefunction provides an accurate description of the system on both sides of the quantum phase transition. In particular, through the analysis of the spin structure factor, ground state fidelity susceptibility and the Binder moment ratio Q(2), a continuous transition from the antiferromagnetic ordered state to the quantum disordered state is found at the critical coupling of α(c) = J(⊥)/J(∥) ≈ 2.62, in good agreement with the result of quantum Monte Carlo simulation. The critical exponent estimated from the finite size scaling analysis (1/ν ≈ 1.4) is consistent with that of the classical 3D Heisenberg universality class.
We study the ground-state phase diagram of the spin-1/21/21/2${J}_{1}\ensuremath{-}{J}_{2}$ Heisenberg model on the square lattice with an accurate bosonic resonating-valence-bond (RVB) wave function. In contrast to the RVB ansatz based on Schwinger fermions, the representation based on Schwinger bosons, supplemented by a variational Monte Carlo technique enforcing the exact projection onto the physical subspace, is able to describe a fully gapped spin liquid in the strongly frustrated regime. In particular, a fully symmetric Z2{Z}_{2}Z2​ spin liquid is stable between two antiferromagnetic phases; a continuous transition at J2=0.4J1{J}_{2}=0.4{J}_{1}J2​=0.4J1​, when the Marshall sign rule begins to be essentially violated, and a first-order transition around J2=0.6J1{J}_{2}=0.6{J}_{1}J2​=0.6J1​ are present. Most importantly, the triplet gap is found to have a nonmonotonic behavior, reaching a maximum around J2=0.51J1{J}_{2}=0.51{J}_{1}J2​=0.51J1​, when the lowest spinon excitation moves from the $\ensuremath{\Gamma}$ to the MMM point, i.e., $k=(\ensuremath{\pi},0)$.
SU ( N ) quantum spin systems may be realized in a variety of physical systems including ultracold atoms in optical lattices. The study of such models also leads to insights into possible novel quantum phases and phase transitions of SU (2) spin models. Here we use Gutzwiller projected fermionic variational wavefunctions to explore the phase diagram and correlation functions of SU ( N ) quantum spin models in the self-conjugate representation. In one dimension, the ground state of the SU (4) spin chain with Heisenberg bilinear and biquadratic interactions is studied by examining instabilities of the Gutzwiller projected free fermion ground state to various broken symmetries. The variational phase diagram so obtained agrees well with exact results. The spin–spin and dimer–dimer correlation functions of the Gutzwiller projected free fermion state with N flavours of fermions are in good agreement with exact and 1/ N calculations for the critical points of SU ( N ) spin chains. In two dimensions, the phase diagram of the antiferromagnetic Heisenberg model on the square lattice is obtained by finding instabilities of the Gutzwiller projected π-flux state. In the absence of biquadratic interactions, the model exhibits long-range Néel order for N = 2 and 4, and spin Peierls (columnar dimer) order for N >4. Upon including biquadratic interactions in the SU (4) model (with a sign appropriate to a fermionic Hubbard model), the Néel order diminishes and eventually disappears, giving way to an extended valence bond crystal. In the case of the SU (6) model, the dimerized ground state melts at sufficiently large biquadratic interaction, yielding a projected π-flux spin liquid phase which in turn undergoes a transition into an extended valence bond crystal at even larger biquadratic interaction. The spin correlations of the projected π-flux state at N = 4 are in good agreement with 1/ N calculations. We find that the state shows strongly enhanced dimer correlations, in qualitative agreement with recent theoretical predictions. We also compare our results with a recent quantum Monte Carlo study of the SU (4) Heisenberg model.
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 use quantum Monte Carlo (stochastic series expansion) and finite-size scaling to study the quantum critical points of two S=1∕2S=1∕2S=1∕2 Heisenberg antiferromagnets in two dimensions: a bilayer and a Kondo-lattice-like system (incomplete bilayer), each with intraplane and interplane couplings JJJ and ${J}_{\ensuremath{\perp}}$. We discuss the ground-state finite-size scaling properties of three different quantities---the Binder moment ratio, the spin stiffness, and the long-wavelength magnetic susceptibility---which we use to extract the critical value of the coupling ratio $g={J}_{\ensuremath{\perp}}∕J$. The individual estimates of gc{g}_{c}gc​ are consistent provided that subleading finite-size corrections are properly taken into account. For both models, we find that the spin stiffness has the smallest subleading finite-size corrections; in the case of the incomplete bilayer we find that the first subleading correction vanishes or is extremely small. In agreement with predictions, we find that at the critical point the Binder ratio has a universal value and the product of the spin stiffness and the long-wavelength susceptibility scales as 1∕L21∕{L}^{2}1∕L2 with a universal prefactor. Our results for the critical coupling ratios are gc=2.5220(1){g}_{c}=2.5220(1)gc​=2.5220(1) (full bilayer) and gc=1.3888(1){g}_{c}=1.3888(1)gc​=1.3888(1) (incomplete bilayer), which represent improvements of more than an order of magnitude over the previous best estimates. For the correlation length exponent we obtain $\ensuremath{\nu}=0.7106(9)$, consistent with the expected three-dimensional Heisenberg universality class.
We use a combination of analytical and numerical techniques to study the phase diagram of the frustrated Heisenberg model on the bilayer honeycomb lattice. Using the Schwinger-boson description of the spin operators followed by a mean-field decoupling, the magnetic phase diagram is studied as a function of the frustration coupling J2{J}_{2}J2​ and the interlayer coupling ${J}_{\ensuremath{\perp}}$. The presence of both magnetically ordered and disordered phases is investigated by means of the evaluation of ground-state energy, spin gap, local magnetization, and spin-spin correlations. We observe a phase with a spin gap and short-range N\'eel correlations that survives for nonzero next-nearest-neighbor interaction and interlayer coupling. Furthermore, we detect signatures of a reentrant behavior in the melting of the N\'eel phase and symmetry restoring when the system undergoes a transition from an on-layer nematic valence-bond crystal phase to an interlayer valence-bond crystal phase. We complement our work with exact diagonalization on small clusters and dimer-series expansion calculations, together with a linear spin-wave approach to study the phase diagram as a function of the spin SSS, the frustration, and the interlayer couplings.