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

Unconventional Phases in Doped or Frustrated Quantum Antiferromagnets: A Systematic QMC Study

Wen-Jun Hu

Published Dec 16, 2013• 0 citations
OpenAlex

Abstract

We investigate the ground-state properties of two strongly-correlated systems: the two-dimensional $t-J$ model that can be used to study the high-temperature superconducting phase in the doped antiferromagnet, and the frustrated N\'eel antiferromagnet described by the $J_{1}-J_{2}$ Heisenberg model on the square lattice, which is widely considered as the prototype model for spin frustration. In this thesis, we apply to these two systems state-of-the-art quantum Monte Carlo techniques, including the variational and the Green's function Monte Carlo with the fixed-node approximation. Few Lanczos steps are used to systematically improve the accuracy of the trial wave functions. By introducing a suitable regularization scheme for the variational calculations with few Lanczos steps, stable and controllable simulations can be performed up to very large cluster sizes with very good accuracy. In the two-dimensional $t{-}J$ model at $J/t=0.4$, we show that the accuracy of the Gutzwiller-projected variational state (containing $d_{x^2-y^2}$ pairing) can be improved much by few Lanczos steps; in addition, the fixed-node Monte Carlo with these systematically improvable trial wave functions gives results that are comparable with the best accurate DMRG ones. Our main outcome is that the ground state is homogeneous and no evidence of stripes is detected around the doping $\delta=1/8$. Indeed, our best approximation to the ground state does not show any tendency towards charge inhomogeneity. Furthermore, our results show that a uniform state containing superconductivity and antiferromagnetism is stabilized at low hole doping, i.e., $\delta \lesssim 0.1$. In the $J_{1}-J_{2}$ Heisenberg model on the square lattice, we use the projected mean-field state that is built from Abrikosov fermions having a $Z_2$ gauge structure and four Dirac points in the spinon spectrum. No spin or dimer order is found in the strongly frustrated regime and our calculations imply that a spin liquid phase may faithfully represent the exact ground state around $J_{2}/J_{1}=0.5$. The few Lanczos step technique is used to systematically improve the accuracy of the variational states both for the ground state and for few relevant low-energy excitations. This procedure allows us to estimate, in a valuable and systemic way, the spin gaps within thermodynamical limit and to show a solid evidence of an unconventional gapless excitation spectrum in the strongly frustrated regime, i.e., $J_{2}/J_{1}\simeq 0.5$. In particular, we found gapless triplet excitations at momenta $(\pi,0)$ and $(0,\pi)$, which are compatible with the presence of four Dirac points at momenta $(\pm \frac{\pi}{2},\pm \frac{\pi}{2})$ in the spinon spectrum.

Related papers

Ground state phases of the spin-1/2 J1–J2 Heisenberg antiferromagnet on the square lattice: A high-order coupled cluster treatment

R. Darradi, Oleg Derzhko, Ronald Zinke +3 more

Using the coupled cluster method for high orders of approximation and complementary exact diagonalization studies we investigate the ground state properties of the spin-1/2 ${J}_{1}--{J}_{2}$ frustrated Heisenberg antiferromagnet on the square lattice. We have calculated the ground state energy, the magnetic order parameter, the spin stiffness, and several generalized susceptibilities to probe magnetically disordered quantum valence-bond phases. We have found that the quantum critical points for both the N\'eel and collinear orders are ${J}_{2}^{c1}\ensuremath{\approx}(0.44\ifmmode\pm\else\textpm\fi{}0.01){J}_{1}$ and ${J}_{2}^{c2}\ensuremath{\approx}(0.59\ifmmode\pm\else\textpm\fi{}0.01){J}_{1}$, respectively, which are in good agreement with the results obtained by other approximations. In contrast to the recent study by [Sirker et al. Phys. Rev. B 73, 184420 (2006)], our data do not provide evidence for the transition from the N\'eel to the valence-bond solid state to be first order. Moreover, our results are in favor of the deconfinement scenario for that phase transition. We also discuss the nature of the magnetically disordered quantum phase.

arXiv:0806.3825•Physical Review B•Dec 11, 2008•183 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

Simulating a two-dimensional frustrated spin system with fermionic resonating-valence-bond states

Chung-Pin Chou, Hong‐Yi Chen

The frustrated Heisenberg ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ model on a square lattice is numerically investigated by variational Monte Carlo simulations. We propose an antiferromagnetic fermion resonating-valence-bond (AF-fRVB) state that has the ability to examine the entire phase diagram in the ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ model. Two phase transition points, the second order around ${J}_{2}/{J}_{1}=0.45$ and the first order around ${J}_{2}/{J}_{1}=0.6$, can be extracted more clearly than the conventional bosonic RVB state. At the maximally frustrated point $({J}_{2}/{J}_{1}=0.5)$, the AF-fRVB state shows the variational ground-state energy in the thermodynamic limit very close to the one estimated by the projected entangled pair state at the largest bond dimension available. On the other hand, in the frustrated regime $0.4\ensuremath{\lesssim}{J}_{2}/{J}_{1}\ensuremath{\le}0.5$, AF-fRVB states with $\text{ext}\phantom{\rule{0.16em}{0ex}}{s}_{2}$ (using the terminology in the field of iron-based superconductors) and ${d}_{xy}$ pairing symmetries are degenerate in the thermodynamic limit, implying the existence of gapless Dirac excitations in the spinon spectrum.

arXiv:1403.6879•Physical Review B•Jul 22, 2014•15 citations

Quantum Spin Liquid in Spin 1/2J1–J2Heisenberg Model on Square Lattice: Many-Variable Variational Monte Carlo Study Combined with Quantum-Number Projections

Satoshi Morita, Ryui Kaneko, Masatoshi Imada

The nature of quantum spin liquids is studied for the spin-$1/2$ antiferromagnetic Heisenberg model on a square lattice containing exchange interactions between nearest-neighbor sites, $J_1$, and those between next-nearest-neighbor sites, $J_2$. We perform variational Monte Carlo simulations together with the quantum-number-projection technique and clarify the phase diagram in the ground state together with its excitation spectra. We obtain the nonmagnetic phase in the region $0.4< J_2/J_1\le 0.6$ sandwiched by the staggered and stripe antiferromagnetic (AF) phases. Our direct calculations of the spin gap support the notion that the triplet excitation from the singlet ground state is gapless in the range of $0.4 < J_2/J_1 \le 0.5$, while the gapped valence-bond-crystal (VBC) phase is stabilized for $0.5 < J_2/J_1 \le 0.6$. The VBC order is likely to have the columnar symmetry with a spontaneous symmetry breaking of the $C_{4v}$ symmetry. The power-law behaviors of the spin-spin and dimer-dimer correlation functions in the gapless region are consistent with the emergence of the algebraic quantum-spin-liquid phase (critical phase). The exponent of the spin correlation $\langle S(0)S(r)\rangle \propto 1/r^{z+η}$ at a long distance $r$ appears to increase from $z+η\sim 1$ at $J_2/J_1\sim0.4$ toward the continuous transition to the VBC phase at $J_1/J_1\sim0.5$. Our results, however, do not fully exclude the possibility of a direct quantum transition between the staggered AF and VBC phases with a wide critical region and deconfined criticality.

arXiv:1410.7557•Journal of the Physical Society of Japan•Jan 29, 2015•122 citations

Near-Term Quantum Spin Simulation of the Spin-$\frac{1}{2}$ Square $J_{1}-J_{2}$ Heisenberg Model

Dylan Sheils, Trevor David Rhone

Simulating complex spin systems, known for high frustration and entanglement, presents significant challenges due to their intricate energy landscapes. This study focuses on the $J_{1}-J_{2}$ Heisenberg model, renowned for its rich phase behavior on the square lattice, to investigate strongly correlated spin systems. We conducted the first experimental quantum computing study of this model using the 127-qubit IBM Rensselear Eagle processor and the Variational Quantum Eigensolver (VQE) algorithm. By employing classical warm-starting ($+40\%$ ground state energy approximation) and a newly developed ansatz ($+9.31\%$ improvement compared to prior best), we improved ground state approximation accuracy on the 16-site variant, achieving usable results with approximately $10^{3}$ iterations, significantly fewer than the $10^{4}-10^{5}$ steps proposed by previous theoretical studies. We utilized existing error mitigation strategies and introduced a novel Classically-Reinforced VQE error mitigation scheme, achieving $93\%$ ground state accuracy, compared to $83.7\%$ with the Quantum Moments algorithm and $60\%$ with standard error mitigation. These strategies reduced the average error of observable prediction from $\approx 20\%$ to $5\%$, enhancing phase prediction from qualitative to quantitative alignment. Additionally, we explored an experimental implementation of the Quantum Lanczos (QLanczos) algorithm using Variational-Fast Forwarding (VFF) on a 4-qubit site, achieving $\approx 97\%$ ground state approximation. Theoretical simulations indicated that Krylov-based methods outperform VQE, with the Lanczos basis converging faster than the real-time basis. Our study demonstrates that near-term quantum devices can predict phase-relevant observables for the $J_1-J_2$ Heisenberg model, transitioning focus from theoretical to experimental, and suggesting general improvements to VQE-based methods.

arXiv:2406.18474•Jun 26, 2024•1 citation

J1−J2 frustrated two-dimensional Heisenberg model: Random phase approximation and functional renormalization group

Johannes Reuther, P. Wölfle

We study the ground-state properties of the two-dimensional spin-1/2 ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model on a square lattice, within diagrammatic approximations using an auxiliary-fermion formulation with exact projection. In a first approximation, we assume a phenomenological width of the pseudofermion spectral function to calculate the magnetization, susceptibilities, and the spin-correlation length within random-phase approximation, demonstrating the appearance of a paramagnetic phase between the N\'eel-ordered and Collinear-ordered phases, at sufficiently large pseudofermion damping. Second we use a functional renormalization-group formulation. We find that the conventional truncation scheme omitting three-particle and higher-order vertices is not sufficient. We therefore include self-energy renormalizations in the single-scale propagator as recently proposed by Katanin, to preserve Ward identities in a better way. We find N\'eel order at $g={J}_{2}/{J}_{1}\ensuremath{\lesssim}{g}_{\text{c}1}\ensuremath{\approx}0.4\dots{}0.45$ and Collinear order at $g\ensuremath{\gtrsim}{g}_{\text{c}2}\ensuremath{\approx}0.66\dots{}0.68$, which is in good agreement with results obtained by numerical studies. In the intervening quantum paramagnetic phase, we find enhanced columnar dimer and plaquette fluctuations of equal strength.

arXiv:0912.0860•Physical Review B•Apr 8, 2010•161 citations
We investigate the ground-state properties of two strongly-correlated systems: the two-dimensional t−Jt-Jt−J model that can be used to study the high-temperature superconducting phase in the doped antiferromagnet, and the frustrated N'eel antiferromagnet described by the J1−J2J_{1}-J_{2}J1​−J2​ Heisenberg model on the square lattice, which is widely considered as the prototype model for spin frustration. In this thesis, we apply to these two systems state-of-the-art quantum Monte Carlo techniques, including the variational and the Green's function Monte Carlo with the fixed-node approximation. Few Lanczos steps are used to systematically improve the accuracy of the trial wave functions. By introducing a suitable regularization scheme for the variational calculations with few Lanczos steps, stable and controllable simulations can be performed up to very large cluster sizes with very good accuracy. In the two-dimensional t−Jt{-}Jt−J model at J/t=0.4J/t=0.4J/t=0.4, we show that the accuracy of the Gutzwiller-projected variational state (containing dx2−y2d_{x^2-y^2}dx2−y2​ pairing) can be improved much by few Lanczos steps; in addition, the fixed-node Monte Carlo with these systematically improvable trial wave functions gives results that are comparable with the best accurate DMRG ones. Our main outcome is that the ground state is homogeneous and no evidence of stripes is detected around the doping δ=1/8\delta=1/8δ=1/8. Indeed, our best approximation to the ground state does not show any tendency towards charge inhomogeneity. Furthermore, our results show that a uniform state containing superconductivity and antiferromagnetism is stabilized at low hole doping, i.e., δ≲0.1\delta \lesssim 0.1δ≲0.1. In the J1−J2J_{1}-J_{2}J1​−J2​ Heisenberg model on the square lattice, we use the projected mean-field state that is built from Abrikosov fermions having a Z2Z_2Z2​ gauge structure and four Dirac points in the spinon spectrum. No spin or dimer order is found in the strongly frustrated regime and our calculations imply that a spin liquid phase may faithfully represent the exact ground state around J2/J1=0.5J_{2}/J_{1}=0.5J2​/J1​=0.5. The few Lanczos step technique is used to systematically improve the accuracy of the variational states both for the ground state and for few relevant low-energy excitations. This procedure allows us to estimate, in a valuable and systemic way, the spin gaps within thermodynamical limit and to show a solid evidence of an unconventional gapless excitation spectrum in the strongly frustrated regime, i.e., J2/J1≃0.5J_{2}/J_{1}\simeq 0.5J2​/J1​≃0.5. In particular, we found gapless triplet excitations at momenta (π,0)(\pi,0)(π,0) and (0,π)(0,\pi)(0,π), which are compatible with the presence of four Dirac points at momenta (±π2,±π2)(\pm \frac{\pi}{2},\pm \frac{\pi}{2})(±2π​,±2π​) in the spinon spectrum.
Using the coupled cluster method for high orders of approximation and complementary exact diagonalization studies we investigate the ground state properties of the spin-1/2 J1−−J2{J}_{1}--{J}_{2}J1​−−J2​ frustrated Heisenberg antiferromagnet on the square lattice. We have calculated the ground state energy, the magnetic order parameter, the spin stiffness, and several generalized susceptibilities to probe magnetically disordered quantum valence-bond phases. We have found that the quantum critical points for both the N\'eel and collinear orders are ${J}_{2}^{c1}\ensuremath{\approx}(0.44\ifmmode\pm\else\textpm\fi{}0.01){J}_{1}$ and ${J}_{2}^{c2}\ensuremath{\approx}(0.59\ifmmode\pm\else\textpm\fi{}0.01){J}_{1}$, respectively, which are in good agreement with the results obtained by other approximations. In contrast to the recent study by [Sirker et al. Phys. Rev. B 73, 184420 (2006)], our data do not provide evidence for the transition from the N\'eel to the valence-bond solid state to be first order. Moreover, our results are in favor of the deconfinement scenario for that phase transition. We also discuss the nature of the magnetically disordered quantum phase.
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.
The frustrated Heisenberg ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ model on a square lattice is numerically investigated by variational Monte Carlo simulations. We propose an antiferromagnetic fermion resonating-valence-bond (AF-fRVB) state that has the ability to examine the entire phase diagram in the ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ model. Two phase transition points, the second order around J2/J1=0.45{J}_{2}/{J}_{1}=0.45J2​/J1​=0.45 and the first order around J2/J1=0.6{J}_{2}/{J}_{1}=0.6J2​/J1​=0.6, can be extracted more clearly than the conventional bosonic RVB state. At the maximally frustrated point (J2/J1=0.5)({J}_{2}/{J}_{1}=0.5)(J2​/J1​=0.5), the AF-fRVB state shows the variational ground-state energy in the thermodynamic limit very close to the one estimated by the projected entangled pair state at the largest bond dimension available. On the other hand, in the frustrated regime $0.4\ensuremath{\lesssim}{J}_{2}/{J}_{1}\ensuremath{\le}0.5$, AF-fRVB states with exts2\text{ext}\phantom{\rule{0.16em}{0ex}}{s}_{2}exts2​ (using the terminology in the field of iron-based superconductors) and dxy{d}_{xy}dxy​ pairing symmetries are degenerate in the thermodynamic limit, implying the existence of gapless Dirac excitations in the spinon spectrum.
The nature of quantum spin liquids is studied for the spin-1/21/21/2 antiferromagnetic Heisenberg model on a square lattice containing exchange interactions between nearest-neighbor sites, J1J_1J1​, and those between next-nearest-neighbor sites, J2J_2J2​. We perform variational Monte Carlo simulations together with the quantum-number-projection technique and clarify the phase diagram in the ground state together with its excitation spectra. We obtain the nonmagnetic phase in the region 0.4<J2/J1≤0.60.4< J_2/J_1\le 0.60.4<J2​/J1​≤0.6 sandwiched by the staggered and stripe antiferromagnetic (AF) phases. Our direct calculations of the spin gap support the notion that the triplet excitation from the singlet ground state is gapless in the range of 0.4<J2/J1≤0.50.4 < J_2/J_1 \le 0.50.4<J2​/J1​≤0.5, while the gapped valence-bond-crystal (VBC) phase is stabilized for 0.5<J2/J1≤0.60.5 < J_2/J_1 \le 0.60.5<J2​/J1​≤0.6. The VBC order is likely to have the columnar symmetry with a spontaneous symmetry breaking of the C4vC_{4v}C4v​ symmetry. The power-law behaviors of the spin-spin and dimer-dimer correlation functions in the gapless region are consistent with the emergence of the algebraic quantum-spin-liquid phase (critical phase). The exponent of the spin correlation ⟨S(0)S(r)⟩∝1/rz+η\langle S(0)S(r)\rangle \propto 1/r^{z+η}⟨S(0)S(r)⟩∝1/rz+η at a long distance rrr appears to increase from z+η∼1z+η\sim 1z+η∼1 at J2/J1∼0.4J_2/J_1\sim0.4J2​/J1​∼0.4 toward the continuous transition to the VBC phase at J1/J1∼0.5J_1/J_1\sim0.5J1​/J1​∼0.5. Our results, however, do not fully exclude the possibility of a direct quantum transition between the staggered AF and VBC phases with a wide critical region and deconfined criticality.
Near-Term Quantum Spin Simulation of the Spin-12\frac{1}{2}21​ Square J1−J2J_{1}-J_{2}J1​−J2​ Heisenberg Model
Simulating complex spin systems, known for high frustration and entanglement, presents significant challenges due to their intricate energy landscapes. This study focuses on the J1−J2J_{1}-J_{2}J1​−J2​ Heisenberg model, renowned for its rich phase behavior on the square lattice, to investigate strongly correlated spin systems. We conducted the first experimental quantum computing study of this model using the 127-qubit IBM Rensselear Eagle processor and the Variational Quantum Eigensolver (VQE) algorithm. By employing classical warm-starting (+40%+40\%+40% ground state energy approximation) and a newly developed ansatz (+9.31%+9.31\%+9.31% improvement compared to prior best), we improved ground state approximation accuracy on the 16-site variant, achieving usable results with approximately 10310^{3}103 iterations, significantly fewer than the 104−10510^{4}-10^{5}104−105 steps proposed by previous theoretical studies. We utilized existing error mitigation strategies and introduced a novel Classically-Reinforced VQE error mitigation scheme, achieving 93%93\%93% ground state accuracy, compared to 83.7%83.7\%83.7% with the Quantum Moments algorithm and 60%60\%60% with standard error mitigation. These strategies reduced the average error of observable prediction from ≈20%\approx 20\%≈20% to 5%5\%5%, enhancing phase prediction from qualitative to quantitative alignment. Additionally, we explored an experimental implementation of the Quantum Lanczos (QLanczos) algorithm using Variational-Fast Forwarding (VFF) on a 4-qubit site, achieving ≈97%\approx 97\%≈97% ground state approximation. Theoretical simulations indicated that Krylov-based methods outperform VQE, with the Lanczos basis converging faster than the real-time basis. Our study demonstrates that near-term quantum devices can predict phase-relevant observables for the J1−J2J_1-J_2J1​−J2​ Heisenberg model, transitioning focus from theoretical to experimental, and suggesting general improvements to VQE-based methods.
We study the ground-state properties of the two-dimensional spin-1/2 ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model on a square lattice, within diagrammatic approximations using an auxiliary-fermion formulation with exact projection. In a first approximation, we assume a phenomenological width of the pseudofermion spectral function to calculate the magnetization, susceptibilities, and the spin-correlation length within random-phase approximation, demonstrating the appearance of a paramagnetic phase between the N\'eel-ordered and Collinear-ordered phases, at sufficiently large pseudofermion damping. Second we use a functional renormalization-group formulation. We find that the conventional truncation scheme omitting three-particle and higher-order vertices is not sufficient. We therefore include self-energy renormalizations in the single-scale propagator as recently proposed by Katanin, to preserve Ward identities in a better way. We find N\'eel order at $g={J}_{2}/{J}_{1}\ensuremath{\lesssim}{g}_{\text{c}1}\ensuremath{\approx}0.4\dots{}0.45$ and Collinear order at $g\ensuremath{\gtrsim}{g}_{\text{c}2}\ensuremath{\approx}0.66\dots{}0.68$, which is in good agreement with results obtained by numerical studies. In the intervening quantum paramagnetic phase, we find enhanced columnar dimer and plaquette fluctuations of equal strength.