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:2602.10008

Ground-state phases of $S = 1/2$ Heisenberg models on the body-centered cubic lattice

Rajah P. Nutakki, Filippo Vicentini

Published Feb 10, 2026• 0 citations
arXivHTMLPDF

Abstract

Simulating low-temperature properties of three-dimensional frustrated quantum magnets is challenging due to the sign problem and the system sizes required to mitigate substantial finite-size effects. However, there are many experimental examples of three-dimensional crystals that could host exotic low-temperature states of matter, such as quantum spin liquids. We calculate the ground-state phase diagrams of frustrated quantum spin models on the body-centered cubic lattice using neural quantum states. First, we study the antiferromagnetic $J_1-J_2$ model where we find a direct first-order phase transition between Néel and collinear long-range-ordered phases at $(J_2/J_1)_c = 0.705$, consistent with previous studies. Then, in a tetragonally-distorted variant, proposed as a minimal model of NaCa$_2$Cu$_2$(VO$_4$)$_3$, we find no evidence of a quantum paramagnetic ground state, with a first-order phase transition between Néel and chain phases at $(J_{2ab}/J_1)_c = 1.0375$. Therefore, the ground state of the tetragonally-distorted model does not reproduce the low-temperature magnetic properties of NaCa$_2$Cu$_2$(VO$_4$)$_3$, and the inclusion of other effects is necessary to rationalize experimental observations.

DOI
10.48550/arxiv.2602.10008

Related papers

Ground-state ordering of the J1−J2 model on the simple cubic and body-centered cubic lattices

DAMIAN J. J. FARNELL, Oliver Götze, Johannes Richter

The ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model is a ``canonical'' model in the field of quantum magnetism in order to study the interplay between frustration and quantum fluctuations as well as quantum phase transitions driven by frustration. Here we apply the coupled cluster method (CCM) to study the spin-half ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ model with antiferromagnetic nearest-neighbor bonds ${J}_{1}>0$ and next-nearest-neighbor bonds ${J}_{2}>0$ for the simple cubic (sc) and body-centered cubic (bcc) lattices. In particular, we wish to study the ground-state ordering of these systems as a function of the frustration parameter $p={z}_{2}{J}_{2}/{z}_{1}{J}_{1}$, where ${z}_{1} ({z}_{2})$ is the number of nearest (next-nearest) neighbors. We wish to determine the positions of the phase transitions using the CCM and we aim to resolve the nature of the phase transition points. We consider the ground-state energy, order parameters, spin-spin correlation functions, as well as the spin stiffness in order to determine the ground-state phase diagrams of these models. We find a direct first-order phase transition at a value of $p=0.528$ from a state of nearest-neighbor N\'eel order to next-nearest-neighbor N\'eel order for the bcc lattice. For the sc lattice the situation is more subtle. CCM results for the energy, the order parameter, the spin-spin correlation functions, and the spin stiffness indicate that there is no direct first-order transition between ground-state phases with magnetic long-range order, rather it is more likely that two phases with antiferromagnetic long range are separated by a narrow region of a spin-liquid-like quantum phase around $p=0.55$. Thus the strong frustration present in the ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model on the sc lattice may open a window for an unconventional quantum ground state in this three-dimensional spin model.

arXiv:1606.00388•Physical review. B./Physical review. B•Jun 14, 2016•39 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}$.

arXiv:1109.6229•Physical Review B•Feb 21, 2012•31 citations

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

Quantum phase transition in the three-dimensional anisotropic frustrated Heisenberg antiferromagnetic model

J. Roberto Viana, J. Ricardo de Sousa, M. A. Contínentino

The phase diagram of the anisotropic quantum spin-1/2 Heisenberg antiferromagnet in the presence of competitive nearest $({J}_{1})$- and next-nearest $({J}_{2}=\ensuremath{\alpha}{J}_{1})$-neighbor interactions on a simple cubic lattice is studied within the framework of the differential operator technique. The Hamiltonian is solved by using an effective-field theory in a two-spin cluster. We propose a function for the free energy, which is similar to the Landau expansion, in order to find the ground state $(T=0)$ and the finite-$T$ phase diagram. At zero temperature, we find a first-order quantum phase transition in the point ${\ensuremath{\alpha}}_{c}(\ensuremath{\Delta})$ between the antiferromagnetic (AF) and superantiferromagnetic (SAF) phases. This is dependent on the anisotropy parameter $\ensuremath{\Delta}$ ($\ensuremath{\Delta}=0$ and 1 correspond to the isotropic Heisenberg and Ising models, respectively). The quantum critical values of ${\ensuremath{\alpha}}_{c}(\ensuremath{\Delta}\ensuremath{\ne}1)={\ensuremath{\alpha}}_{c}^{q}(\ensuremath{\Delta})$ are smaller than those of the classical model ${\ensuremath{\alpha}}_{c}(\ensuremath{\Delta}=1)={\ensuremath{\alpha}}_{c}^{c}=1/4$ (for example, for $\ensuremath{\Delta}=0$, we find ${\ensuremath{\alpha}}_{c}^{q}(0)=0.21$). The former monotonically increases with the anisotropy parameter $\ensuremath{\Delta}$. The phase diagram in $T\text{\ensuremath{-}}\ensuremath{\alpha}$ plane presents a second-order transition at ${T}_{2}(\ensuremath{\alpha})$ between the AF and paramagnetic (P) phases. The phase transition between the SAF and P phases is of first-order ${T}_{1}(\ensuremath{\alpha})$ for all values of $\ensuremath{\Delta}∊[0,1]$. The phase transition lines ${T}_{1}(\ensuremath{\alpha})$ and ${T}_{2}(\ensuremath{\alpha})$ merge in a critical end point (CEP) that is dependent on $\ensuremath{\Delta}$ [i.e., ${T}_{\text{CEP}}(\ensuremath{\Delta})$ and ${\ensuremath{\alpha}}_{\text{CEP}}(\ensuremath{\Delta})$ are the temperature and frustration parameters at the CEP]. The ordered phases (AF and SAF) are separated by first-order transitions and exhibit a nontrivial re-entrant behavior at low temperatures.

Physical Review B•May 27, 2008•24 citations

Ground-state phases of the spin- 1 J 1 − J 2 Heisenberg antiferromagnet on the honeycomb lattice

P. H. Y. Li, R. F. Bishop

We study the zero-temperature quantum phase diagram of a spin-1 Heisenberg antiferromagnet on the honeycomb lattice with both nearest-neighbor exchange coupling ${J}_{1}>0$ and frustrating next-nearest-neighbor coupling ${J}_{2}\ensuremath{\equiv}\ensuremath{\kappa}{J}_{1}>0$, using the coupled cluster method implemented to high orders of approximation, and based on model states with different forms of classical magnetic order. For each we calculate directly in the bulk thermodynamic limit both ground-state low-energy parameters (including the energy per spin, magnetic order parameter, spin stiffness coefficient, and zero-field uniform transverse magnetic susceptibility) and their generalized susceptibilities to various forms of valence-bond crystalline (VBC) order, as well as the energy gap to the lowest-lying spin-triplet excitation. In the range $0<\ensuremath{\kappa}<1$ we find evidence for four distinct phases. Two of these are quasiclassical phases with antiferromagnetic long-range order, one with two-sublattice N\'eel order for $\ensuremath{\kappa}<{\ensuremath{\kappa}}_{{c}_{1}}=0.250(5)$, and another with four-sublattice N\'eel-II order for $\ensuremath{\kappa}>{\ensuremath{\kappa}}_{{c}_{2}}=0.340(5)$. Two different paramagnetic phases are found to exist in the intermediate region. Over the range ${\ensuremath{\kappa}}_{{c}_{1}}<\ensuremath{\kappa}<{\ensuremath{\kappa}}_{c}^{i}=0.305(5)$ we find a gapless phase with no discernible magnetic order, which is a strong candidate for being a quantum spin liquid, while over the range ${\ensuremath{\kappa}}_{c}^{i}<\ensuremath{\kappa}<{\ensuremath{\kappa}}_{{c}_{2}}$ we find a gapped phase, which is most likely a lattice nematic with staggered dimer VBC order that breaks the lattice rotational symmetry.

arXiv:1602.08915•Physical review. B./Physical review. B•Jun 28, 2016•21 citations

Quantum phase transitions in the two-dimensional J1−J2 model

Oleg P. Sushkov, J. Oitmaa, Zheng Weihong

We analyze the phase diagram of the frustrated Heisenberg antiferromagnet, the ${J}_{1}{\ensuremath{-}J}_{2}$ model, in two dimensions. Two quantum phase transitions in the model are already known: the second-order transition from the N\'eel state to the spin liquid state at ${(J}_{2}{/J}_{1}{)}_{c2}\ensuremath{\approx}0.38,$ and the first-order transition from the spin liquid state to the collinear state at ${(J}_{2}{/J}_{1}{)}_{c4}\ensuremath{\approx}0.60.$ We have found evidence for two second-order phase transitions: the transition from the spin columnar dimerized state to the state with plaquette-type modulation at ${(J}_{2}{/J}_{1}{)}_{c3}=0.50\ifmmode\pm\else\textpm\fi{}0.02,$ and the transition from the simple N\'eel state to the N\'eel state with spin columnar dimerization at ${(J}_{2}{/J}_{1}{)}_{c1}=0.34\ifmmode\pm\else\textpm\fi{}0.04.$ We also present an independent calculation of ${(J}_{2}{/J}_{1}{)}_{c2}\ensuremath{\approx}0.38$ using another approach.

arXiv:cond-mat/0007329•Physical review. B, Condensed matter•Feb 20, 2001•151 citations
Ground-state phases of S=1/2S = 1/2S=1/2 Heisenberg models on the body-centered cubic lattice
Simulating low-temperature properties of three-dimensional frustrated quantum magnets is challenging due to the sign problem and the system sizes required to mitigate substantial finite-size effects. However, there are many experimental examples of three-dimensional crystals that could host exotic low-temperature states of matter, such as quantum spin liquids. We calculate the ground-state phase diagrams of frustrated quantum spin models on the body-centered cubic lattice using neural quantum states. First, we study the antiferromagnetic J1−J2J_1-J_2J1​−J2​ model where we find a direct first-order phase transition between Néel and collinear long-range-ordered phases at (J2/J1)c=0.705(J_2/J_1)_c = 0.705(J2​/J1​)c​=0.705, consistent with previous studies. Then, in a tetragonally-distorted variant, proposed as a minimal model of NaCa2_22​Cu2_22​(VO4_44​)3_33​, we find no evidence of a quantum paramagnetic ground state, with a first-order phase transition between Néel and chain phases at (J2ab/J1)c=1.0375(J_{2ab}/J_1)_c = 1.0375(J2ab​/J1​)c​=1.0375. Therefore, the ground state of the tetragonally-distorted model does not reproduce the low-temperature magnetic properties of NaCa2_22​Cu2_22​(VO4_44​)3_33​, and the inclusion of other effects is necessary to rationalize experimental observations.
The ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model is a ``canonical'' model in the field of quantum magnetism in order to study the interplay between frustration and quantum fluctuations as well as quantum phase transitions driven by frustration. Here we apply the coupled cluster method (CCM) to study the spin-half ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ model with antiferromagnetic nearest-neighbor bonds J1>0{J}_{1}>0J1​>0 and next-nearest-neighbor bonds J2>0{J}_{2}>0J2​>0 for the simple cubic (sc) and body-centered cubic (bcc) lattices. In particular, we wish to study the ground-state ordering of these systems as a function of the frustration parameter p=z2J2/z1J1p={z}_{2}{J}_{2}/{z}_{1}{J}_{1}p=z2​J2​/z1​J1​, where z1(z2){z}_{1} ({z}_{2})z1​(z2​) is the number of nearest (next-nearest) neighbors. We wish to determine the positions of the phase transitions using the CCM and we aim to resolve the nature of the phase transition points. We consider the ground-state energy, order parameters, spin-spin correlation functions, as well as the spin stiffness in order to determine the ground-state phase diagrams of these models. We find a direct first-order phase transition at a value of p=0.528p=0.528p=0.528 from a state of nearest-neighbor N\'eel order to next-nearest-neighbor N\'eel order for the bcc lattice. For the sc lattice the situation is more subtle. CCM results for the energy, the order parameter, the spin-spin correlation functions, and the spin stiffness indicate that there is no direct first-order transition between ground-state phases with magnetic long-range order, rather it is more likely that two phases with antiferromagnetic long range are separated by a narrow region of a spin-liquid-like quantum phase around p=0.55p=0.55p=0.55. Thus the strong frustration present in the ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model on the sc lattice may open a window for an unconventional quantum ground state in this three-dimensional spin model.
We study the ground-state (GS) properties of the frustrated spin-12\frac{1}{2}21​ J1{J}_{1}J1​--J2{J}_{2}J2​--J3{J}_{3}J3​ Heisenberg model on the two-dimensional honeycomb lattice with ferromagnetic nearest-neighbor (${J}_{1}=\ensuremath{-}1$) exchange and frustrating antiferromagnetic next-nearest-neighbor (J2>0{J}_{2}>0J2​>0) and next-next-nearest-neighbor (J3>0{J}_{3}>0J3​>0) exchanges, for the case J3=J2{J}_{3}={J}_{2}J3​=J2​. 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 J2cl=17{J}_{2}^{\mathrm{cl}}=\frac{1}{7}J2cl​=71​ ($\ensuremath{\approx}0.143$) for the classical version of the model, in which we momentarily ignore the intervening noncollinear spiral phase in the region 110<J2<15\frac{1}{10}<{J}_{2}<\frac{1}{5}101​<J2​<51​. 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 J1=+1{J}_{1}=+1J1​=+1). The magnetic order parameter is found to behave similarly for the ferromagnetic and the antiferromagnetic models for large values of the frustration parameter J2{J}_{2}J2​. 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 J1=+1{J}_{1}=+1J1​=+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 J1=+1{J}_{1}=+1J1​=+1), where a plaquette valence-bond crystal phase intrudes between the two corresponding quasiclassical antiferromagnetic phases (with N\'eel and striped order) for J2c1<J2<J2c2{J}_{2}^{{c}_{1}}<{J}_{2}<{J}_{2}^{{c}_{2}}J2c1​​<J2​<J2c2​​, 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 110<J2<15\frac{1}{10}<{J}_{2}<\frac{1}{5}101​<J2​<51​.
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 phase diagram of the anisotropic quantum spin-1/2 Heisenberg antiferromagnet in the presence of competitive nearest (J1)({J}_{1})(J1​)- and next-nearest $({J}_{2}=\ensuremath{\alpha}{J}_{1})$-neighbor interactions on a simple cubic lattice is studied within the framework of the differential operator technique. The Hamiltonian is solved by using an effective-field theory in a two-spin cluster. We propose a function for the free energy, which is similar to the Landau expansion, in order to find the ground state (T=0)(T=0)(T=0) and the finite-TTT phase diagram. At zero temperature, we find a first-order quantum phase transition in the point ${\ensuremath{\alpha}}_{c}(\ensuremath{\Delta})$ between the antiferromagnetic (AF) and superantiferromagnetic (SAF) phases. This is dependent on the anisotropy parameter $\ensuremath{\Delta}$ ($\ensuremath{\Delta}=0$ and 1 correspond to the isotropic Heisenberg and Ising models, respectively). The quantum critical values of ${\ensuremath{\alpha}}_{c}(\ensuremath{\Delta}\ensuremath{\ne}1)={\ensuremath{\alpha}}_{c}^{q}(\ensuremath{\Delta})$ are smaller than those of the classical model ${\ensuremath{\alpha}}_{c}(\ensuremath{\Delta}=1)={\ensuremath{\alpha}}_{c}^{c}=1/4$ (for example, for $\ensuremath{\Delta}=0$, we find ${\ensuremath{\alpha}}_{c}^{q}(0)=0.21$). The former monotonically increases with the anisotropy parameter $\ensuremath{\Delta}$. The phase diagram in $T\text{\ensuremath{-}}\ensuremath{\alpha}$ plane presents a second-order transition at ${T}_{2}(\ensuremath{\alpha})$ between the AF and paramagnetic (P) phases. The phase transition between the SAF and P phases is of first-order ${T}_{1}(\ensuremath{\alpha})$ for all values of $\ensuremath{\Delta}∊[0,1]$. The phase transition lines ${T}_{1}(\ensuremath{\alpha})$ and ${T}_{2}(\ensuremath{\alpha})$ merge in a critical end point (CEP) that is dependent on $\ensuremath{\Delta}$ [i.e., ${T}_{\text{CEP}}(\ensuremath{\Delta})$ and ${\ensuremath{\alpha}}_{\text{CEP}}(\ensuremath{\Delta})$ are the temperature and frustration parameters at the CEP]. The ordered phases (AF and SAF) are separated by first-order transitions and exhibit a nontrivial re-entrant behavior at low temperatures.
We study the zero-temperature quantum phase diagram of a spin-1 Heisenberg antiferromagnet on the honeycomb lattice with both nearest-neighbor exchange coupling J1>0{J}_{1}>0J1​>0 and frustrating next-nearest-neighbor coupling ${J}_{2}\ensuremath{\equiv}\ensuremath{\kappa}{J}_{1}>0$, using the coupled cluster method implemented to high orders of approximation, and based on model states with different forms of classical magnetic order. For each we calculate directly in the bulk thermodynamic limit both ground-state low-energy parameters (including the energy per spin, magnetic order parameter, spin stiffness coefficient, and zero-field uniform transverse magnetic susceptibility) and their generalized susceptibilities to various forms of valence-bond crystalline (VBC) order, as well as the energy gap to the lowest-lying spin-triplet excitation. In the range $0<\ensuremath{\kappa}<1$ we find evidence for four distinct phases. Two of these are quasiclassical phases with antiferromagnetic long-range order, one with two-sublattice N\'eel order for $\ensuremath{\kappa}<{\ensuremath{\kappa}}_{{c}_{1}}=0.250(5)$, and another with four-sublattice N\'eel-II order for $\ensuremath{\kappa}>{\ensuremath{\kappa}}_{{c}_{2}}=0.340(5)$. Two different paramagnetic phases are found to exist in the intermediate region. Over the range ${\ensuremath{\kappa}}_{{c}_{1}}<\ensuremath{\kappa}<{\ensuremath{\kappa}}_{c}^{i}=0.305(5)$ we find a gapless phase with no discernible magnetic order, which is a strong candidate for being a quantum spin liquid, while over the range ${\ensuremath{\kappa}}_{c}^{i}<\ensuremath{\kappa}<{\ensuremath{\kappa}}_{{c}_{2}}$ we find a gapped phase, which is most likely a lattice nematic with staggered dimer VBC order that breaks the lattice rotational symmetry.
We analyze the phase diagram of the frustrated Heisenberg antiferromagnet, the ${J}_{1}{\ensuremath{-}J}_{2}$ model, in two dimensions. Two quantum phase transitions in the model are already known: the second-order transition from the N\'eel state to the spin liquid state at ${(J}_{2}{/J}_{1}{)}_{c2}\ensuremath{\approx}0.38,$ and the first-order transition from the spin liquid state to the collinear state at ${(J}_{2}{/J}_{1}{)}_{c4}\ensuremath{\approx}0.60.$ We have found evidence for two second-order phase transitions: the transition from the spin columnar dimerized state to the state with plaquette-type modulation at ${(J}_{2}{/J}_{1}{)}_{c3}=0.50\ifmmode\pm\else\textpm\fi{}0.02,$ and the transition from the simple N\'eel state to the N\'eel state with spin columnar dimerization at ${(J}_{2}{/J}_{1}{)}_{c1}=0.34\ifmmode\pm\else\textpm\fi{}0.04.$ We also present an independent calculation of ${(J}_{2}{/J}_{1}{)}_{c2}\ensuremath{\approx}0.38$ using another approach.