arXiv:1711.10121

Revisiting a variational study of the quantum $J_{1}-J_{1}^{'}-J_{2}$ spin-1/2 Heisenberg antiferromagnet

Kevin J. U. Vidarte, Luigui A. Chero Nunura, Octavio D. Rodrı́guez Salmon, José Ricardo de Sousa

Published Nov 28, 2017• 0 citations

Abstract

The variational study of the ground state of the spin$-1/2$ anisotropic Heisenberg antiferromagnet has been revisited on a square lattice by improving and correcting past numerical results found in Sol. State. Comm. 165, 33 (2013). The Hamiltonian has been implemented on a square lattice with antiferromagnetic interactions between nearest- and next-nearest neighbors. The nearest-neighbor couplings have different strengths, namely, $J_{1}$ and $J_{1}^{'}$, for the x and y directions, respectively. These couplings compete with the next-nearest ones denoted by $J_{2}$. We obtained a new phase diagram in the $λ- α$ plane, where $λ= J_{1}^{'}/J_{1}$ and $α= J_{2}/J_{1}$, whose topology is slightly different of that previously found. There is no direct frontier dividing the collinear (CAF) and the antiferromagnetic order (AF), rather, the quantum paramagnetic phase (QP) separates these two phases for all positive values of $λ$ and $α$. The true nature of the frontiers has been obtained by scanning rigorously the relevant points of the $λ-α$ plane.

DOI: 10.48550/arxiv.1711.10121

Related papers

Phase diagram of a frustrated Heisenberg antiferromagnet on the honeycomb lattice: The J1 - J2 - J3 model

P. H. Y. Li, R. F. Bishop, DAMIAN J. J. FARNELL +1 more

arXiv:1204.2390Physical Review BOct 9, 201299 citations

Ground-state phase diagram of the spin- 12 square-lattice J1 - J2 model with plaquette structure

Oliver Götze, Sven E. Krüger, F. Fleck +2 more

arXiv:1205.0681Physical Review BJun 21, 201216 citations

Spin- 12 J1 - J2 Heisenberg antiferromagnet on a square lattice: A plaquette renormalized tensor network study

Ji-Feng Yu, Ying-Jer Kao

arXiv:1108.3393Physical Review BMar 5, 201270 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

arXiv:1410.7557Journal of the Physical Society of JapanJan 29, 2015122 citations

Quantum phase diagram of spin-1 J1−J2 Heisenberg model on the square lattice: An infinite projected entangled-pair state and density matrix renormalization group study

Reza Haghshenas, Wang-Wei Lan, Shou-Shu Gong +1 more

arXiv:1802.00874Physical review. B./Physical review. BMay 31, 201821 citations

Ground-state phases of the frustrated spin- 12 J1 – J2 – J3 Heisenberg ferromagnet ( J1<0 ) on the honeycomb lattice with J3=J2>0

P. H. Y. Li, R. F. Bishop, DAMIAN J. J. FARNELL +2 more

arXiv:1109.6229Physical Review BFeb 21, 201231 citations

arXiv:1711.10121

Revisiting a variational study of the quantum $J_{1}-J_{1}^{'}-J_{2}$ spin-1/2 Heisenberg antiferromagnetRevisiting a variational study of the quantum $J_{1}-J_{1}^{'}-J_{2}$ spin-1/2 Heisenberg antiferromagnet

Kevin J. U. Vidarte, Luigui A. Chero Nunura, Octavio D. Rodrı́guez Salmon, José Ricardo de Sousa

Published Nov 28, 2017• 0 citations

arXiv HTML PDF

Abstract

The variational study of the ground state of the spin

-1/2

anisotropic Heisenberg antiferromagnet has been revisited on a square lattice by improving and correcting past numerical results found in Sol. State. Comm. 165, 33 (2013). The Hamiltonian has been implemented on a square lattice with antiferromagnetic interactions between nearest- and next-nearest neighbors. The nearest-neighbor couplings have different strengths, namely,

J_{1}

and

J_{1}^{'}

, for the x and y directions, respectively. These couplings compete with the next-nearest ones denoted by

J_{2}

. We obtained a new phase diagram in the

λ- α

plane, where

λ= J_{1}^{'}/J_{1}

and

α= J_{2}/J_{1}

, whose topology is slightly different of that previously found. There is no direct frontier dividing the collinear (CAF) and the antiferromagnetic order (AF), rather, the quantum paramagnetic phase (QP) separates these two phases for all positive values of

λ

and

α

. The true nature of the frontiers has been obtained by scanning rigorously the relevant points of the

λ-α

plane.

DOI: 10.48550/arxiv.1711.10121

Related papers

Phase diagram of a frustrated Heisenberg antiferromagnet on the honeycomb lattice: The J1 - J2 - J3 model

P. H. Y. Li, R. F. Bishop, DAMIAN J. J. FARNELL +1 more

We use the coupled-cluster method in high orders of approximation to make a comprehensive study of the ground-state (GS) phase diagram of the spin-1/2 ${J}_{1}$ - ${J}_{2}$ - ${J}_{3}$ model on a two-dimensional honeycomb lattice with antiferromagnetic (AFM) interactions up to third-nearest neighbors. Results are presented for the GS energy and the average local onsite magnetization. With the nearest-neighbor coupling strength ${J}_{1}\ensuremath{\equiv}1$, we find four magnetically ordered phases in the parameter window ${J}_{2},{J}_{3}\ensuremath{\in}[0,1]$, namely, the N\'eel, striped, and N\'eel-II collinear AFM phases, plus a spiral phase. The N\'eel-II phase appears as a stable GS phase in the classical version of the model only for values ${J}_{3}<0$ . Each of these four ordered phases shares a boundary with a disordered quantum paramagnetic (QP) phase, and at several widely separated points on the phase boundaries the QP phase has an infinite susceptibility to plaquette valence-bond crystalline order. We identify all of the phase boundaries with good precision in the parameter window studied, and we find three tricritical quantum critical points therein at (a) $({J}_{2}^{{c}_{1}},{J}_{3}^{{c}_{1}})=(0.51\ifmmode\pm\else\textpm\fi{}0.01,0.69\ifmmode\pm\else\textpm\fi{}0.01)$ between the N\'eel, striped, and QP phases; (b) $({J}_{2}^{{c}_{2}},{J}_{3}^{{c}_{2}})=(0.65\ifmmode\pm\else\textpm\fi{}0.02,0.55\ifmmode\pm\else\textpm\fi{}0.01)$ between the striped, spiral, and QP phases; and (c) $({J}_{2}^{{c}_{3}},{J}_{3}^{{c}_{3}})=(0.69\ifmmode\pm\else\textpm\fi{}0.01,0.12\ifmmode\pm\else\textpm\fi{}0.01)$ between the spiral, N\'eel-II, and QP phases.We use the coupled-cluster method in high orders of approximation to make a comprehensive study of the ground-state (GS) phase diagram of the spin-1/2 ${J}_{1}$-${J}_{2}$-${J}_{3}$ model on a two-dimensional honeycomb lattice with antiferromagnetic (AFM) interactions up to third-nearest neighbors. Results are presented for the GS energy and the average local onsite magnetization. With the nearest-neighbor coupling strength ${J}_{1}\ensuremath{\equiv}1$, we find four magnetically ordered phases in the parameter window ${J}_{2},{J}_{3}\ensuremath{\in}[0,1]$, namely, the N\'eel, striped, and N\'eel-II collinear AFM phases, plus a spiral phase. The N\'eel-II phase appears as a stable GS phase in the classical version of the model only for values ${J}_{3}<0$. Each of these four ordered phases shares a boundary with a disordered quantum paramagnetic (QP) phase, and at several widely separated points on the phase boundaries the QP phase has an infinite susceptibility to plaquette valence-bond crystalline order. We identify all of the phase boundaries with good precision in the parameter window studied, and we find three tricritical quantum critical points therein at (a) $({J}_{2}^{{c}_{1}},{J}_{3}^{{c}_{1}})=(0.51\ifmmode\pm\else\textpm\fi{}0.01,0.69\ifmmode\pm\else\textpm\fi{}0.01)$ between the N\'eel, striped, and QP phases; (b) $({J}_{2}^{{c}_{2}},{J}_{3}^{{c}_{2}})=(0.65\ifmmode\pm\else\textpm\fi{}0.02,0.55\ifmmode\pm\else\textpm\fi{}0.01)$ between the striped, spiral, and QP phases; and (c) $({J}_{2}^{{c}_{3}},{J}_{3}^{{c}_{3}})=(0.69\ifmmode\pm\else\textpm\fi{}0.01,0.12\ifmmode\pm\else\textpm\fi{}0.01)$ between the spiral, N\'eel-II, and QP phases.

arXiv:1204.2390Physical Review BOct 9, 201299 citations

Ground-state phase diagram of the spin- 12 square-lattice J1 - J2 model with plaquette structure

Oliver Götze, Sven E. Krüger, F. Fleck +2 more

Using the coupled-cluster method for high orders of approximation and Lanczos exact diagonalization we study the ground-state phase diagram of a quantum spin- $\frac{1}{2}$ ${J}_{1}$ - ${J}_{2}$ model on the square lattice with plaquette structure. We consider antiferromagnetic ( ${J}_{1}>0$ ) as well as ferromagnetic ( ${J}_{1}<0$ ) nearest-neighbor interactions together with frustrating antiferromagnetic next-nearest-neighbor interaction ${J}_{2}>0$ . The strength of interplaquette interaction $\ensuremath{\lambda}$ varies between $\ensuremath{\lambda}=1$ (which corresponds to the uniform ${J}_{1}$ - ${J}_{2}$ model) and $\ensuremath{\lambda}=0$ (which corresponds to isolated frustrated 4-spin plaquettes). While on the classical level ($s\ensuremath{\rightarrow}\ensuremath{\infty}$) both versions of models (i.e., with ferro- and antiferromagnetic ${J}_{1}$ ) exhibit the same ground-state behavior, the ground-state phase diagram differs basically for the quantum case $s=1/2$ . For the antiferromagnetic case ( ${J}_{1}>0$ ) N\'eel antiferromagnetic long-range order at small ${J}_{2}/{J}_{1}$ and $\ensuremath{\lambda}\ensuremath{\gtrsim}0.47$ as well as collinear striped antiferromagnetic long-range order at large ${J}_{2}/{J}_{1}$ and $\ensuremath{\lambda}\ensuremath{\gtrsim}0.30$ appear which correspond to their classical counterparts. Both semiclassical magnetic phases are separated by a nonmagnetic quantum paramagnetic phase. The parameter region, where this nonmagnetic phase exists, increases with decreasing $\ensuremath{\lambda}$. For the ferromagnetic case ( ${J}_{1}<0$ ) we have the trivial ferromagnetic ground state at small ${J}_{2}/|{J}_{1}|$ . By increasing ${J}_{2}$ this classical phase gives way for a semiclassical plaquette phase, where the plaquette block spins of length $s=2$ are antiferromagnetically long-range ordered. Further increasing of ${J}_{2}$ then yields collinear striped antiferromagnetic long-range order for $\ensuremath{\lambda}\ensuremath{\gtrsim}0.38$, but a nonmagnetic quantum paramagnetic phase $\ensuremath{\lambda}\ensuremath{\lesssim}0.38$.Using the coupled-cluster method for high orders of approximation and Lanczos exact diagonalization we study the ground-state phase diagram of a quantum spin-$\frac{1}{2}$ ${J}_{1}$-${J}_{2}$ model on the square lattice with plaquette structure. We consider antiferromagnetic (${J}_{1}>0$) as well as ferromagnetic (${J}_{1}<0$) nearest-neighbor interactions together with frustrating antiferromagnetic next-nearest-neighbor interaction ${J}_{2}>0$. The strength of interplaquette interaction $\ensuremath{\lambda}$ varies between $\ensuremath{\lambda}=1$ (which corresponds to the uniform ${J}_{1}$-${J}_{2}$ model) and $\ensuremath{\lambda}=0$ (which corresponds to isolated frustrated 4-spin plaquettes). While on the classical level ($s\ensuremath{\rightarrow}\ensuremath{\infty}$) both versions of models (i.e., with ferro- and antiferromagnetic ${J}_{1}$) exhibit the same ground-state behavior, the ground-state phase diagram differs basically for the quantum case $s=1/2$. For the antiferromagnetic case (${J}_{1}>0$) N\'eel antiferromagnetic long-range order at small ${J}_{2}/{J}_{1}$ and $\ensuremath{\lambda}\ensuremath{\gtrsim}0.47$ as well as collinear striped antiferromagnetic long-range order at large ${J}_{2}/{J}_{1}$ and $\ensuremath{\lambda}\ensuremath{\gtrsim}0.30$ appear which correspond to their classical counterparts. Both semiclassical magnetic phases are separated by a nonmagnetic quantum paramagnetic phase. The parameter region, where this nonmagnetic phase exists, increases with decreasing $\ensuremath{\lambda}$. For the ferromagnetic case (${J}_{1}<0$) we have the trivial ferromagnetic ground state at small ${J}_{2}/|{J}_{1}|$. By increasing ${J}_{2}$ this classical phase gives way for a semiclassical plaquette phase, where the plaquette block spins of length $s=2$ are antiferromagnetically long-range ordered. Further increasing of ${J}_{2}$ then yields collinear striped antiferromagnetic long-range order for $\ensuremath{\lambda}\ensuremath{\gtrsim}0.38$, but a nonmagnetic quantum paramagnetic phase $\ensuremath{\lambda}\ensuremath{\lesssim}0.38$.

arXiv:1205.0681Physical Review BJun 21, 201216 citations

Spin- 12 J1 - J2 Heisenberg antiferromagnet on a square lattice: A plaquette renormalized tensor network study

Ji-Feng Yu, Ying-Jer Kao

We apply the plaquette renormalization scheme of tensor network states [Phys. Rev. E 83, 056703 (2011)] to study the spin-1/2 frustrated Heisenberg ${J}_{1}$ - ${J}_{2}$ model on an $L\ifmmode\times\else\texttimes\fi{}L$ square lattice with $L$ $=$ 8, 16, and 32. By treating tensor elements as variational parameters, we obtain the ground states for different ${J}_{2}/{J}_{1}$ values, and investigate staggered magnetizations, nearest-neighbor spin-spin correlations, and plaquette order parameters. In addition to the well-known N\'eel order and collinear order at low and high ${J}_{2}/{J}_{1}$ , we observe a plaquettelike order at ${J}_{2}/{J}_{1}\ensuremath{\approx}0.5$. A continuous transition between the N\'eel order and the plaquettelike order near ${J}_{2}^{{c}_{1}}\ensuremath{\approx}0.40{J}_{1}$ is observed. The collinear order emerges at ${J}_{2}^{{c}_{2}}\ensuremath{\approx}0.62{J}_{1}$ through a first-order phase transition.We apply the plaquette renormalization scheme of tensor network states [Phys. Rev. E 83, 056703 (2011)] to study the spin-1/2 frustrated Heisenberg ${J}_{1}$-${J}_{2}$ model on an $L\ifmmode\times\else\texttimes\fi{}L$ square lattice with $L$ $=$ 8, 16, and 32. By treating tensor elements as variational parameters, we obtain the ground states for different ${J}_{2}/{J}_{1}$ values, and investigate staggered magnetizations, nearest-neighbor spin-spin correlations, and plaquette order parameters. In addition to the well-known N\'eel order and collinear order at low and high ${J}_{2}/{J}_{1}$, we observe a plaquettelike order at ${J}_{2}/{J}_{1}\ensuremath{\approx}0.5$. A continuous transition between the N\'eel order and the plaquettelike order near ${J}_{2}^{{c}_{1}}\ensuremath{\approx}0.40{J}_{1}$ is observed. The collinear order emerges at ${J}_{2}^{{c}_{2}}\ensuremath{\approx}0.62{J}_{1}$ through a first-order phase transition.

arXiv:1108.3393Physical Review BMar 5, 201270 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.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.7557Journal of the Physical Society of JapanJan 29, 2015122 citations

Quantum phase diagram of spin-1 J1−J2 Heisenberg model on the square lattice: An infinite projected entangled-pair state and density matrix renormalization group study

Reza Haghshenas, Wang-Wei Lan, Shou-Shu Gong +1 more

We study the spin-1 Heisenberg model on the square lattice with the antiferromagnetic nearest-neighbor ${J}_{1}$ and the next-nearest-neighbor ${J}_{2}$ couplings by using the infinite projected entangled-pair state (iPEPS) ansatz and the density matrix renormalization group (DMRG) calculation. The iPEPS simulation, which studies the model directly in the thermodynamic limit, finds a crossing of the ground state from the N\'eel magnetic state to the stripe magnetic state at ${J}_{2}/{J}_{1}\ensuremath{\simeq}0.549$, showing a direct phase transition. In the finite-size DMRG calculation on the cylinder geometry up to the cylinder width ${L}_{y}=10$ , we find that around the same critical point the N\'eel and the stripe orders are strongly suppressed, which implies the absence of an intermediate phase. Both calculations identify that the stripe order comes with a first-order transition at ${J}_{2}/{J}_{1}\ensuremath{\simeq}0.549$. Our results indicate that unlike the spin-1/2 ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ square model, quantum fluctuations in the spin-1 model may not be strong enough to stabilize an intermediate nonmagnetic phase.We study the spin-1 Heisenberg model on the square lattice with the antiferromagnetic nearest-neighbor ${J}_{1}$ and the next-nearest-neighbor ${J}_{2}$ couplings by using the infinite projected entangled-pair state (iPEPS) ansatz and the density matrix renormalization group (DMRG) calculation. The iPEPS simulation, which studies the model directly in the thermodynamic limit, finds a crossing of the ground state from the N\'eel magnetic state to the stripe magnetic state at ${J}_{2}/{J}_{1}\ensuremath{\simeq}0.549$, showing a direct phase transition. In the finite-size DMRG calculation on the cylinder geometry up to the cylinder width ${L}_{y}=10$, we find that around the same critical point the N\'eel and the stripe orders are strongly suppressed, which implies the absence of an intermediate phase. Both calculations identify that the stripe order comes with a first-order transition at ${J}_{2}/{J}_{1}\ensuremath{\simeq}0.549$. Our results indicate that unlike the spin-1/2 ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ square model, quantum fluctuations in the spin-1 model may not be strong enough to stabilize an intermediate nonmagnetic phase.

arXiv:1802.00874Physical review. B./Physical review. BMay 31, 201821 citations

Ground-state phases of the frustrated spin- 12 J1 – J2 – J3 Heisenberg ferromagnet ( J1<0 ) on the honeycomb lattice with J3=J2>0

P. H. Y. Li, R. F. Bishop, DAMIAN J. J. FARNELL +2 more

We study the ground-state (GS) properties of the frustrated spin- $\frac{1}{2}$ ${J}_{1}$ -- ${J}_{2}$ -- ${J}_{3}$ Heisenberg model on the two-dimensional honeycomb lattice with ferromagnetic nearest-neighbor (${J}_{1}=\ensuremath{-}1$) exchange and frustrating antiferromagnetic next-nearest-neighbor ( ${J}_{2}>0$ ) and next-next-nearest-neighbor ( ${J}_{3}>0$ ) exchanges, for the case ${J}_{3}={J}_{2}$ . We use the coupled-cluster method implemented to high orders of approximation, complemented by the Lanczos exact diagonalization of a large finite lattice with 32 sites, in order to calculate the GS energy, magnetic order parameter, and spin-spin correlation functions. In one scenario we find a quantum phase transition point between regions characterized by ferromagnetic order and a form of antiferromagnetic (``striped'') collinear order at ${J}_{2}^{c}\ensuremath{\approx}0.1095\ifmmode\pm\else\textpm\fi{}0.0005$, which is below the corresponding hypothetical transition point at ${J}_{2}^{\mathrm{cl}}=\frac{1}{7}$ ($\ensuremath{\approx}0.143$) for the classical version of the model, in which we momentarily ignore the intervening noncollinear spiral phase in the region $\frac{1}{10}<{J}_{2}<\frac{1}{5}$ . Hence we see that quantum fluctuations appear to stabilize somewhat the collinear antiferromagnetic order in preference to the ferromagnetic order in this model. We compare results for the present ferromagnetic case (with ${J}_{1}=\ensuremath{-}1$) to previous results for the corresponding antiferromagnetic case (with ${J}_{1}=+1$ ). The magnetic order parameter is found to behave similarly for the ferromagnetic and the antiferromagnetic models for large values of the frustration parameter ${J}_{2}$ . However, there are considerable differences in the behavior of the order parameters for the two models for ${J}_{2}/|{J}_{1}|\ensuremath{\lesssim}0.6$. For example, the quasiclassical collinear magnetic long-range order for the antiferromagnetic model (with ${J}_{1}=+1$ ) breaks down at ${J}_{2}^{{c}_{2}}\ensuremath{\approx}0.60$, whereas the ``equivalent'' point for the ferromagnetic model (with ${J}_{1}=\ensuremath{-}1$) occurs at ${J}_{2}^{c}\ensuremath{\approx}0.11$. Unlike in the antiferromagnetic model (with ${J}_{1}=+1$ ), where a plaquette valence-bond crystal phase intrudes between the two corresponding quasiclassical antiferromagnetic phases (with N\'eel and striped order) for ${J}_{2}^{{c}_{1}}<{J}_{2}<{J}_{2}^{{c}_{2}}$ , with ${J}_{2}^{{c}_{1}}\ensuremath{\approx}0.47$, we find no clear indications at all in the ferromagnetic model for an intermediate magnetically disordered phase between the corresponding phases exhibiting ferromagnetic and striped order. Instead the evidence for the ferromagnetic model (with ${J}_{1}=\ensuremath{-}1$) points to one of two scenarios: either there is a direct first-order transition between the two magnetically ordered phases, as mentioned above; or there exists an intervening phase between them in the very narrow range $0.10\ensuremath{\lesssim}{J}_{2}\ensuremath{\lesssim}0.12$, which is probably a remnant of the spiral phase that exists in the classical counterpart of the model over the larger range $\frac{1}{10}<{J}_{2}<\frac{1}{5}$ .We study the ground-state (GS) properties of the frustrated spin-$\frac{1}{2}$ ${J}_{1}$--${J}_{2}$--${J}_{3}$ Heisenberg model on the two-dimensional honeycomb lattice with ferromagnetic nearest-neighbor (${J}_{1}=\ensuremath{-}1$) exchange and frustrating antiferromagnetic next-nearest-neighbor (${J}_{2}>0$) and next-next-nearest-neighbor (${J}_{3}>0$) exchanges, for the case ${J}_{3}={J}_{2}$. We use the coupled-cluster method implemented to high orders of approximation, complemented by the Lanczos exact diagonalization of a large finite lattice with 32 sites, in order to calculate the GS energy, magnetic order parameter, and spin-spin correlation functions. In one scenario we find a quantum phase transition point between regions characterized by ferromagnetic order and a form of antiferromagnetic (``striped'') collinear order at ${J}_{2}^{c}\ensuremath{\approx}0.1095\ifmmode\pm\else\textpm\fi{}0.0005$, which is below the corresponding hypothetical transition point at ${J}_{2}^{\mathrm{cl}}=\frac{1}{7}$ ($\ensuremath{\approx}0.143$) for the classical version of the model, in which we momentarily ignore the intervening noncollinear spiral phase in the region $\frac{1}{10}<{J}_{2}<\frac{1}{5}$. Hence we see that quantum fluctuations appear to stabilize somewhat the collinear antiferromagnetic order in preference to the ferromagnetic order in this model. We compare results for the present ferromagnetic case (with ${J}_{1}=\ensuremath{-}1$) to previous results for the corresponding antiferromagnetic case (with ${J}_{1}=+1$). The magnetic order parameter is found to behave similarly for the ferromagnetic and the antiferromagnetic models for large values of the frustration parameter ${J}_{2}$. However, there are considerable differences in the behavior of the order parameters for the two models for ${J}_{2}/|{J}_{1}|\ensuremath{\lesssim}0.6$. For example, the quasiclassical collinear magnetic long-range order for the antiferromagnetic model (with ${J}_{1}=+1$) breaks down at ${J}_{2}^{{c}_{2}}\ensuremath{\approx}0.60$, whereas the ``equivalent'' point for the ferromagnetic model (with ${J}_{1}=\ensuremath{-}1$) occurs at ${J}_{2}^{c}\ensuremath{\approx}0.11$. Unlike in the antiferromagnetic model (with ${J}_{1}=+1$), where a plaquette valence-bond crystal phase intrudes between the two corresponding quasiclassical antiferromagnetic phases (with N\'eel and striped order) for ${J}_{2}^{{c}_{1}}<{J}_{2}<{J}_{2}^{{c}_{2}}$, with ${J}_{2}^{{c}_{1}}\ensuremath{\approx}0.47$, we find no clear indications at all in the ferromagnetic model for an intermediate magnetically disordered phase between the corresponding phases exhibiting ferromagnetic and striped order. Instead the evidence for the ferromagnetic model (with ${J}_{1}=\ensuremath{-}1$) points to one of two scenarios: either there is a direct first-order transition between the two magnetically ordered phases, as mentioned above; or there exists an intervening phase between them in the very narrow range $0.10\ensuremath{\lesssim}{J}_{2}\ensuremath{\lesssim}0.12$, which is probably a remnant of the spiral phase that exists in the classical counterpart of the model over the larger range $\frac{1}{10}<{J}_{2}<\frac{1}{5}$.

arXiv:1109.6229Physical Review BFeb 21, 201231 citations

Revisiting a variational study of the quantum

J_{1}-J_{1}^{'}-J_{2}

spin-1/2 Heisenberg antiferromagnet

The variational study of the ground state of the spin

-1/2

J_{1}

and

J_{1}^{'}

, for the x and y directions, respectively. These couplings compete with the next-nearest ones denoted by

J_{2}

. We obtained a new phase diagram in the

λ- α

plane, where

λ= J_{1}^{'}/J_{1}

and

α= J_{2}/J_{1}

λ

and

α

. The true nature of the frontiers has been obtained by scanning rigorously the relevant points of the

λ-α

plane.

We use the coupled-cluster method in high orders of approximation to make a comprehensive study of the ground-state (GS) phase diagram of the spin-1/2

{J}_{1}

{J}_{2}

{J}_{3}

model on a two-dimensional honeycomb lattice with antiferromagnetic (AFM) interactions up to third-nearest neighbors. Results are presented for the GS energy and the average local onsite magnetization. With the nearest-neighbor coupling strength ${J}_{1}\ensuremath{\equiv}1$, we find four magnetically ordered phases in the parameter window ${J}_{2},{J}_{3}\ensuremath{\in}[0,1]$, namely, the N\'eel, striped, and N\'eel-II collinear AFM phases, plus a spiral phase. The N\'eel-II phase appears as a stable GS phase in the classical version of the model only for values

{J}_{3}<0

. Each of these four ordered phases shares a boundary with a disordered quantum paramagnetic (QP) phase, and at several widely separated points on the phase boundaries the QP phase has an infinite susceptibility to plaquette valence-bond crystalline order. We identify all of the phase boundaries with good precision in the parameter window studied, and we find three tricritical quantum critical points therein at (a) $({J}_{2}^{{c}_{1}},{J}_{3}^{{c}_{1}})=(0.51\ifmmode\pm\else\textpm\fi{}0.01,0.69\ifmmode\pm\else\textpm\fi{}0.01)$ between the N\'eel, striped, and QP phases; (b) $({J}_{2}^{{c}_{2}},{J}_{3}^{{c}_{2}})=(0.65\ifmmode\pm\else\textpm\fi{}0.02,0.55\ifmmode\pm\else\textpm\fi{}0.01)$ between the striped, spiral, and QP phases; and (c) $({J}_{2}^{{c}_{3}},{J}_{3}^{{c}_{3}})=(0.69\ifmmode\pm\else\textpm\fi{}0.01,0.12\ifmmode\pm\else\textpm\fi{}0.01)$ between the spiral, N\'eel-II, and QP phases.

Using the coupled-cluster method for high orders of approximation and Lanczos exact diagonalization we study the ground-state phase diagram of a quantum spin-

\frac{1}{2}

{J}_{1}

{J}_{2}

model on the square lattice with plaquette structure. We consider antiferromagnetic (

{J}_{1}>0

) as well as ferromagnetic (

{J}_{1}<0

) nearest-neighbor interactions together with frustrating antiferromagnetic next-nearest-neighbor interaction

{J}_{2}>0

. The strength of interplaquette interaction $\ensuremath{\lambda}$ varies between $\ensuremath{\lambda}=1$ (which corresponds to the uniform

{J}_{1}

{J}_{2}

model) and $\ensuremath{\lambda}=0$ (which corresponds to isolated frustrated 4-spin plaquettes). While on the classical level ($s\ensuremath{\rightarrow}\ensuremath{\infty}$) both versions of models (i.e., with ferro- and antiferromagnetic

{J}_{1}

) exhibit the same ground-state behavior, the ground-state phase diagram differs basically for the quantum case

s=1/2

. For the antiferromagnetic case (

{J}_{1}>0

) N\'eel antiferromagnetic long-range order at small

{J}_{2}/{J}_{1}

and $\ensuremath{\lambda}\ensuremath{\gtrsim}0.47$ as well as collinear striped antiferromagnetic long-range order at large

{J}_{2}/{J}_{1}

and $\ensuremath{\lambda}\ensuremath{\gtrsim}0.30$ appear which correspond to their classical counterparts. Both semiclassical magnetic phases are separated by a nonmagnetic quantum paramagnetic phase. The parameter region, where this nonmagnetic phase exists, increases with decreasing $\ensuremath{\lambda}$. For the ferromagnetic case (

{J}_{1}<0

) we have the trivial ferromagnetic ground state at small

{J}_{2}/|{J}_{1}|

. By increasing

{J}_{2}

this classical phase gives way for a semiclassical plaquette phase, where the plaquette block spins of length

s=2

are antiferromagnetically long-range ordered. Further increasing of

{J}_{2}

then yields collinear striped antiferromagnetic long-range order for $\ensuremath{\lambda}\ensuremath{\gtrsim}0.38$, but a nonmagnetic quantum paramagnetic phase $\ensuremath{\lambda}\ensuremath{\lesssim}0.38$.

We apply the plaquette renormalization scheme of tensor network states [Phys. Rev. E 83, 056703 (2011)] to study the spin-1/2 frustrated Heisenberg

{J}_{1}

{J}_{2}

model on an $L\ifmmode\times\else\texttimes\fi{}L$ square lattice with

L

=

8, 16, and 32. By treating tensor elements as variational parameters, we obtain the ground states for different

{J}_{2}/{J}_{1}

values, and investigate staggered magnetizations, nearest-neighbor spin-spin correlations, and plaquette order parameters. In addition to the well-known N\'eel order and collinear order at low and high

{J}_{2}/{J}_{1}

, we observe a plaquettelike order at ${J}_{2}/{J}_{1}\ensuremath{\approx}0.5$. A continuous transition between the N\'eel order and the plaquettelike order near ${J}_{2}^{{c}_{1}}\ensuremath{\approx}0.40{J}_{1}$ is observed. The collinear order emerges at ${J}_{2}^{{c}_{2}}\ensuremath{\approx}0.62{J}_{1}$ through a first-order phase transition.

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

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.

We study the spin-1 Heisenberg model on the square lattice with the antiferromagnetic nearest-neighbor

{J}_{1}

and the next-nearest-neighbor

{J}_{2}

couplings by using the infinite projected entangled-pair state (iPEPS) ansatz and the density matrix renormalization group (DMRG) calculation. The iPEPS simulation, which studies the model directly in the thermodynamic limit, finds a crossing of the ground state from the N\'eel magnetic state to the stripe magnetic state at ${J}_{2}/{J}_{1}\ensuremath{\simeq}0.549$, showing a direct phase transition. In the finite-size DMRG calculation on the cylinder geometry up to the cylinder width

{L}_{y}=10

, we find that around the same critical point the N\'eel and the stripe orders are strongly suppressed, which implies the absence of an intermediate phase. Both calculations identify that the stripe order comes with a first-order transition at ${J}_{2}/{J}_{1}\ensuremath{\simeq}0.549$. Our results indicate that unlike the spin-1/2 ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ square model, quantum fluctuations in the spin-1 model may not be strong enough to stabilize an intermediate nonmagnetic phase.

We study the ground-state (GS) properties of the frustrated spin-

\frac{1}{2}

{J}_{1}

{J}_{2}

{J}_{3}

Heisenberg model on the two-dimensional honeycomb lattice with ferromagnetic nearest-neighbor (${J}_{1}=\ensuremath{-}1$) exchange and frustrating antiferromagnetic next-nearest-neighbor (

{J}_{2}>0

) and next-next-nearest-neighbor (

{J}_{3}>0

) exchanges, for the case

{J}_{3}={J}_{2}

. We use the coupled-cluster method implemented to high orders of approximation, complemented by the Lanczos exact diagonalization of a large finite lattice with 32 sites, in order to calculate the GS energy, magnetic order parameter, and spin-spin correlation functions. In one scenario we find a quantum phase transition point between regions characterized by ferromagnetic order and a form of antiferromagnetic (``striped'') collinear order at ${J}_{2}^{c}\ensuremath{\approx}0.1095\ifmmode\pm\else\textpm\fi{}0.0005$, which is below the corresponding hypothetical transition point at

{J}_{2}^{\mathrm{cl}}=\frac{1}{7}

($\ensuremath{\approx}0.143$) for the classical version of the model, in which we momentarily ignore the intervening noncollinear spiral phase in the region

\frac{1}{10}<{J}_{2}<\frac{1}{5}

. Hence we see that quantum fluctuations appear to stabilize somewhat the collinear antiferromagnetic order in preference to the ferromagnetic order in this model. We compare results for the present ferromagnetic case (with ${J}_{1}=\ensuremath{-}1$) to previous results for the corresponding antiferromagnetic case (with

{J}_{1}=+1

). The magnetic order parameter is found to behave similarly for the ferromagnetic and the antiferromagnetic models for large values of the frustration parameter

{J}_{2}

. However, there are considerable differences in the behavior of the order parameters for the two models for ${J}_{2}/|{J}_{1}|\ensuremath{\lesssim}0.6$. For example, the quasiclassical collinear magnetic long-range order for the antiferromagnetic model (with

{J}_{1}=+1

) breaks down at ${J}_{2}^{{c}_{2}}\ensuremath{\approx}0.60$, whereas the ``equivalent'' point for the ferromagnetic model (with ${J}_{1}=\ensuremath{-}1$) occurs at ${J}_{2}^{c}\ensuremath{\approx}0.11$. Unlike in the antiferromagnetic model (with

{J}_{1}=+1

), where a plaquette valence-bond crystal phase intrudes between the two corresponding quasiclassical antiferromagnetic phases (with N\'eel and striped order) for

{J}_{2}^{{c}_{1}}<{J}_{2}<{J}_{2}^{{c}_{2}}

, with ${J}_{2}^{{c}_{1}}\ensuremath{\approx}0.47$, we find no clear indications at all in the ferromagnetic model for an intermediate magnetically disordered phase between the corresponding phases exhibiting ferromagnetic and striped order. Instead the evidence for the ferromagnetic model (with ${J}_{1}=\ensuremath{-}1$) points to one of two scenarios: either there is a direct first-order transition between the two magnetically ordered phases, as mentioned above; or there exists an intervening phase between them in the very narrow range $0.10\ensuremath{\lesssim}{J}_{2}\ensuremath{\lesssim}0.12$, which is probably a remnant of the spiral phase that exists in the classical counterpart of the model over the larger range

\frac{1}{10}<{J}_{2}<\frac{1}{5}