We construct a family of short-range resonating-valence-bond wave functions on a layered cubic lattice, allowing for a tunable anisotropy in the amplitudes assigned to nearest-neighbour valence bonds along one axis. Monte Carlo simulations reveal that four phases are stabilized over the full range of the anisotropy parameter. They are separated from one another by a sequence of continuous quantum phase transitions. An antiferromagnetic phase, centred on the perfect isotropy point, intervenes between two distinct quantum spin liquid states. One of them is continuously deformable to the two-dimensional U(1) spin liquid, which is known to exhibit critical bond correlations. The other has both spin and bond correlations that decay exponentially. The existence of this second phase is proof that, contrary to expectations, neither a bipartite lattice structure nor a conventional Marshall sign rule is an impediment to realizing a fully gapped quantum spin liquid.
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
Z2 spin liquid is stable between two antiferromagnetic phases; a continuous transition at
J2=0.4J1, when the Marshall sign rule begins to be essentially violated, and a first-order transition around
J2=0.6J1 are present. Most importantly, the triplet gap is found to have a nonmonotonic behavior, reaching a maximum around
J2=0.51J1, when the lowest spinon excitation moves from the $\ensuremath{\Gamma}$ to the
M point, i.e., $k=(\ensuremath{\pi},0)$.
We investigate the arguably simplest
SU(2)-invariant wave functions capable of accounting for spin-liquid behavior, expressed in terms of nearest-neighbor valence-bond states on the square lattice and characterized by different topological invariants. While such wave functions are known to exhibit short-range spin correlations, we perform Monte Carlo simulations and show that four-point correlations decay algebraically with an exponent 1.16(4). This is reminiscent of the classical dimer problem, albeit with a slower decay. Furthermore, these correlators are found to be spatially modulated according to a wave vector related to the topological invariants. We conclude that a recently proposed spin Hamiltonian that stabilizes the here considered wave function(s) as its (degenerate) ground state(s) should exhibit gapped spin and gapless nonmagnetic excitations.
Classification of possible quantum spin liquid (QSL) states of interacting spin-1/2's in two dimensions has been a fascinating topic of condensed matter for decades, resulting in enormous progress in our understanding of low-dimensional quantum matter. By contrast, relatively little work exists on the identification, let alone classification, of QSL phases for spin-1 systems in dimensions higher than one. Employing the powerful ideas of tensor network theory and its classification, we develop general methods for writing QSL wave functions of spin-1 respecting all the lattice symmetries, spin rotation, and time reversal with trivial gauge structure on the square lattice. We find
25 distinct classes characterized by five binary quantum numbers. Several explicit constructions of such wave functions are given for bond dimensions
D ranging from two to four, along with thorough numerical analyses to identify their physical characters. Both gapless and gapped states are found. The topological entanglement entropy of the gapped states is close to zero, indicative of topologically trivial states. In
D=4, several different tensors can be linearly combined to produce a family of states within the same symmetry class. A rich ``phase diagram'' can be worked out among the phases of these tensors, as well as the phase transitions among them. Among the states we identified in this putative phase diagram is the plaquette-ordered phase, gapped resonating valence bond phase, and a critical phase. A continuous transition separates the plaquette-ordered phase from the resonating valence bond phase.
Abstract We analyze a Higgs transition from a U(1) Dirac spin liquid to a gapless ℤ2 spin liquid. This ℤ2 spin liquid is of relevance to the spin S = 1/2 square lattice antiferromagnet, where recent numerical studies have given evidence for such a phase existing in the regime of high frustration between nearest neighbor and next-nearest neighbor antiferromagnetic interactions (the J1-J2 model), appearing in a parameter regime between the vanishing of Néel order and the onset of valence bond solid ordering. The proximate Dirac spin liquid is unstable to monopole proliferation on the square lattice, ultimately leading to Néel or valence bond solid ordering. As such, we conjecture that this Higgs transition describes the critical theory separating the gapless ℤ2 spin liquid of the J1-J2 model from one of the two proximate ordered phases. The transition into the other ordered phase can be described in a unified manner via a transition into an unstable SU(2) spin liquid, which we have analyzed in prior work. By studying the deconfined critical theory separating the U(1) Dirac spin liquid from the gapless ℤ2 spin liquid in a 1/Nf expansion, with Nf proportional to the number of fermions, we find a stable fixed point with an anisotropic spinon dispersion and a dynamical critical exponent z ≠ 1. We analyze the consequences of this anisotropic dispersion by calculating the angular profiles of the equal-time Néel and valence bond solid correlation functions, and we find them to be distinct. We also note the influence of the anisotropy on the scaling dimension of monopoles.
We present results for the phase diagram of an SU(
N) generalization of the Heisenberg antiferromagnet on a bipartite three-dimensional anisotropic cubic (tetragonal) lattice as a function of
N and the lattice anisotropy $\ensuremath{\gamma}$. In the ``isotropic'' $\ensuremath{\gamma}=1$ cubic limit, we find a transition from N\'eel to valence bond solid (VBS) between
N=9 and
N=10. We follow the N\'eel-VBS transition to the limiting cases of $\ensuremath{\gamma}\ensuremath{\ll}1$ (weakly coupled layers) and $\ensuremath{\gamma}\ensuremath{\gg}1$ (weakly coupled chains). Throughout the phase diagram we find a direct first-order transition from N\'eel at small
N to VBS at large
N. In the three-dimensional models studied here, we find no evidence for either an intervening spin-liquid ``photon'' phase or a continuous transition, even close to the limit $\ensuremath{\gamma}\ensuremath{\ll}1$ where the isolated layers undergo continuous N\'eel-VBS transitions.
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,
J1, and those between next-nearest-neighbor sites,
J2. 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.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.5, while the gapped valence-bond-crystal (VBC) phase is stabilized for
0.5<J2/J1≤0.6. The VBC order is likely to have the columnar symmetry with a spontaneous symmetry breaking of the
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+η at a long distance
r appears to increase from
z+η∼1 at
J2/J1∼0.4 toward the continuous transition to the VBC phase at
J1/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.