We study the
T=0 phase diagrams of models of bilayers of
S=1/2 square lattices antiferromagnets with SU(2) Heisenberg symmetry that have 2, 4, and 6 spin exchanges. We study two families of bilayer models with distinct internal symmetries and, hence, different phase diagram topologies. A traditional bilayer model in which the interlayer interaction is Heisenberg so that the two layers can exchange spin (and energy) with each other, making it possible to achieve a simple dimerized valence bond liquid-like state. The resulting phase diagram is rich with Néel, valence bond solid and simple dimer phases, and both first-order and continuous transitions, which we demonstrate are consistent with the conventional Landau theory of order parameters. In the second family of models in which the layers can exchange only energy but no spin (reminiscent of the Ashkin-Teller coupling), the simple dimer state cannot occur. The phase diagrams reveal a number of phase transitions that are accessed for the first time. We find that the phase transition between Néel and VBS is first order in both the spin-spin and energy-energy coupled models, although they have strikingly distinct finite-size scaling behavior and that the transition from VBS to dimer in the spin-spin coupling model deviates from the expected scenario of an XY model with dangerously irrelevant four-fold anisotropy.
We present a detailed study of the destruction of SU(
N) magnetic order in square lattice bilayer antiferromagnets using unbiased quantum Monte Carlo numerical simulations and field theoretic techniques. We study phase transitions from an SU(
N) N\'eel state into two distinct quantum disordered ``valence-bond'' phases: a valence-bond liquid (VBL) with no broken symmetries and a lattice-symmetry-breaking valence-bond solid (VBS) state. For finite interlayer coupling, the cancellation of Berry phases between the layers has dramatic consequences on the two phase transitions: the N\'eel-VBS transition is first order for all $N\ensuremath{\ge}5$ accesible in our model, whereas the N\'eel-VBL transition is continuous for
N=2 and first order for $N\ensuremath{\ge}4$; for
N=3 the N\'eel-VBL transition show no signs of first-order behavior.
We study the N\'eel to fourfold columnar valence bond solid (cVBS) quantum phase transition in a sign-free
S=1 square-lattice model. This is the same kind of transition that for
S=1/2 has been argued to realize the prototypical deconfined critical point. Extensive numerical simulations of the square-lattice
S=1/2 N\'eel-VBS transition have found consistency with the deconfined critical point scenario with no direct evidence for first-order behavior. In contrast to the
S=1/2 case, in our quantum Monte Carlo simulations for the
S=1 model, we present unambiguous evidence for a direct conventional first-order quantum phase transition. Classic signs of a first-order transition demonstrating coexistence including double-peaked histograms and switching behavior are observed. The sharp contrast from the
S=1/2 case is remarkable; we hypothesize that this is a striking demonstration of the role of the size of the quantum spin in the phase diagram of two-dimensional lattice models.
We use quantum Monte Carlo methods to study the ground-state phase diagram of a
S=1/2 honeycomb lattice magnet in which a nearest-neighbor antiferromagnetic exchange
J (favoring N\'eel order) competes with two different multispin interaction terms: a six-spin interaction
Q3 that favors columnar valence-bond solid (VBS) order, and a four-spin interaction
Q2 that favors staggered VBS order. For ${Q}_{3}\ensuremath{\sim}{Q}_{2}\ensuremath{\gg}J$, we establish that the competition between the two different VBS orders stabilizes N\'eel order in a large swath of the phase diagram even when
J is the smallest energy scale in the Hamiltonian. When ${Q}_{3}\ensuremath{\gg}({Q}_{2},J)$ $\mathrm{[}{Q}_{2}\ensuremath{\gg}({Q}_{3},J)\text{]}$, this model exhibits at zero temperature phase transition from the N\'eel state to a columnar (staggered) VBS state. We establish that the N\'eel-columnar VBS transition is continuous for all values of
Q2, and that critical properties along the entire phase boundary are well characterized by critical exponents and amplitudes of the noncompact
CP1 (NCCP1) theory of deconfined criticality, similar to what is observed on a square lattice. However, a surprising threefold anisotropy of the phase of the VBS order parameter at criticality, whose presence was recently noted at the
Q2=0 deconfined critical point, is seen to persist all along this phase boundary. We use a classical analogy to explore this by studying the critical point of a three-dimensional
XY model with a fourfold anisotropy field which is known to be weakly irrelevant at the three-dimensional
XY critical point. In this case, we again find that the critical anisotropy appears to saturate to a nonzero value over the range of sizes accessible to our simulations.
We report on a valence bond projector Monte Carlo simulation of the cubic lattice quantum Heisenberg model with additional higher-order exchange interactions in each unit cell. The model supports two different valence bond solid (VBS) ground states. In one of these states, the dimer pattern is a three-dimensional analogue of the columnar pattern familiar from two dimensions. In the other, the dimers are regularly arranged along the four main diagonals in 1/8 of the unit cells. The phases are separated from one another and from a Néel phase by strongly first-order boundaries. Our results strengthen the case for exotic transitions in two dimensions, where no discontinuities have been detected at the Heisenberg Néel-VBS transition driven by four-spin plaquette interactions.
We study a spin-1/2 SU(2) model on the honeycomb lattice with nearest-neighbor antiferromagnetic exchange J that favors Néel order and competing six-spin interactions Q that favor a valence-bond-solid (VBS) state in which the bond energies order at the "columnar" wave vector K=(2π/3,-2π/3). We present quantum Monte Carlo evidence for a direct continuous quantum phase transition between Néel and VBS states, with exponents and logarithmic violations of scaling consistent with those at analogous deconfined critical points on the square lattice. Although this strongly suggests a description in terms of deconfined criticality, the measured threefold anisotropy of the phase of the VBS order parameter shows unusual near-marginal behavior at the critical point.