The large-N limit of a nearest-neighbor SU(N) antiferromagnet on a bipartite lattice exhibits in dimensions d\ensuremath{\ge}2 a zero-temperature phase transition between a N\'eel ordered state and a resonating-valence-bond state. Here it is shown in d=1,2 that topological effects produce spin-Peierls or valence-bond-solid order in the non-N\'eel phase with a ground-state degeneracy which varies periodically with ``spin'' for fixed N with periodicity given by the coordination number of the lattice. Thus a non-N\'eel phase of the spin-(1/2 Heisenberg model on a square lattice would be a spin-Peierls state with a fourfold degeneracy due to broken lattice rotational symmetry.
We examine a recently developed sine-Gordon model description of the non-N\'eel phase of quantum antiferromagnets on two-dimensional bipartite lattices. Using a spatial dimensionality d=1+\ensuremath{\epsilon} expansion we argue that the model always scales to its strong-coupling limit and displays spin Peierls or valence-bond-solid order. The structure of the theory in this strong-coupling limit bears a remarkable resemblance to a fermionic large-N limit of the nearest-neighbor SU(N) antiferromagnet.
We investigate the ground state phase diagram of an SU(
N)-symmetric antiferromagnetic spin model on a square lattice where each site hosts an irreducible representation of SU(
N) described by a square Young tableau of
N/2 rows and
2S columns. We show that negative sign free fermion Monte Carlo simulations can be carried out for this class of quantum magnets at any
S and even values of
N. In the large-
N limit, the saddle point approximation favors a four-fold degenerate valence bond solid phase. In the large
S-limit, the semi-classical approximation points to Néel state. On a line set by
N=8S+2 in the
S versus
N phase diagram, we observe a variety of phases proximate to the Néel state. At
S=1/2 and
3/2 we observe the aforementioned four fold degenerate valence bond solid state. At
S=1 a two fold degenerate spin nematic state in which the C
4 lattice symmetry is broken down to C
2 emerges. Finally at
S=2 we observe a unique ground state that pertains to a two-dimensional version of the Affleck-Kennedy-Lieb-Tasaki state. For our specific realization, this symmetry protected topological state is characterized by an SU(18),
S=1/2 boundary state, that has a dimerized ground state. These phases that are proximate to the Néel state are consistent with the notion of monopole condensation of the antiferromagnetic order parameter. In particular one expects spin disordered states with degeneracy set by
mod(4,2S).
We generalize the SU(N=2) S=1/2 square-lattice quantum magnet with nearest-neighbor antiferromagnetic coupling (J(1)) and next-nearest-neighbor ferromagnetic coupling (J(2)) to arbitrary N. For all N>4, the ground state has valence-bond-solid order for J(2)=0 and Néel order for J(2)/J(1)≫1, allowing us access to the transition between these types of states for large N. Using quantum Monte Carlo simulations, we show that both order parameters vanish at a single quantum-critical point, whose universal exponents for large enough N (here up to N=12) approach the values obtained in a 1/N expansion of the noncompact CP(N-1) field theory. These results lend strong support to the deconfined quantum-criticality theory of the Néel-valence-bond-solid transition.
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 use stochastic series expansion (SSE) quantum Monte Carlo (QMC) methods to study the phases and transitions displayed by a class of sign-free designer Hamiltonians for
SU(N) analogs of spin
S=1 quantum antiferromagnets on the square lattice. The
SU(N) spins are generators of the single-row two-column representation (complex conjugate of single-row two-column representation) on
A (
B) sublattices of the square lattice, and the Hamiltonian is designed to explore the competition between the nearest-neighbor antiferromagnetic exchange couplings
J and four-spin interactions
Q that favor a plaquette-ordered valence bond solid (p-VBS) ground state. We find that this state is indeed established at large
Q/J for all
N>3. For $3<N\ensuremath{\le}9$, the ground state exhibits a direct first-order quantum phase transition from a small-
Q/J N\'eel ordered antiferromagnetic state to this large-
Q/J p-VBS state. The ground state at
Q/J=0 for $N\ensuremath{\ge}10$ has been previously reported to be a valence bond nematic state, dubbed the Haldane nematic in recent literature. For small nonzero
Q/J and $N\ensuremath{\ge}10$, we additionally find an unusual intermediate state in which the four-site plaquette energy has a Bragg peak at wave vector $(\ensuremath{\pi},\ensuremath{\pi})$ with no accompanying Bragg peaks at wave vectors $(\ensuremath{\pi},0)$ and $(0,\ensuremath{\pi})$. This $(\ensuremath{\pi},\ensuremath{\pi})$ state also appears to be metastable at
Q/J=0, as evidenced by low-temperature histograms of the Haldane nematic and $(\ensuremath{\pi},\ensuremath{\pi})$ order parameters. Deep in the p-VBS phase of the ground state phase diagram, we find the temperature-driven melting of the p-VBS order is in the Ashkin-Teller universality class. In this regime, we identify an interesting signature of Ashkin-Teller criticality in four-site plaquette singlet correlations at wave vector $(\ensuremath{\pi},\ensuremath{\pi})$; this is in addition to the expected critical fluctuations of the conventional p-VBS order parameter at wave vectors $(\ensuremath{\pi},0)$ and $(0,\ensuremath{\pi})$.