The
SU(N) symmetric antiferromagnetic Heisenberg model with multicolumn representations on the two-dimensional square lattice is investigated by quantum Monte Carlo simulations. For the representation of a Young diagram with two columns, we confirm that a valence-bond solid (VBS) order appears as soon as the N\'eel order disappears at
N=10, indicating no intermediate phase. In the case of the representation with three columns, there is no evidence for either the N\'eel or the VBS ordering for $N\ensuremath{\ge}15$. This is actually consistent with the large-
N theory, which predicts that the VBS state immediately follows the N\'eel state, because the expected spontaneous order is too weak to be detected.
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 map out the ground-state phase diagram of a one-dimensional
SU(N) Kondo-Heisenberg lattice model at half-filling and in the fully antisymmetric self-adjoint representation as a function of
N1 and Kondo coupling
Jk/t. On the basis of auxiliary field quantum Monte Carlo (QMC) simulations with even
N up to 8, we show that the enlarged $\mathrm{SU}(N\ensuremath{\ge}4$) symmetry realizes a quantum phase transition separating a valence bond solid (VBS) phase occupying a weak-coupling part of the phase diagram and the Kondo insulator (KI) state dominating in the strong
Jk/t limit. Along the phase boundary, we always observe critical exponents that belong to a two-dimensional classical Ising universality class. We next trace the evolution of the composite fermion and spin spectral functions across the phase boundary and conclude that VBS order triggers a bond order wave state of conduction electrons, both coexisting with Kondo screening. Upon further reducing
Jk/t, we observe that composite quasiparticles lose their spectral weight indicating that VBS order gradually liberates localized
f spins from forming Kondo singlets with conduction electron spins. We contrast the QMC results with a static large-
N approximation, and we show that in the limit of infinite degeneracy $N\ensuremath{\rightarrow}\ensuremath{\infty}$ the order-disorder transition becomes of first order and is accompanied by a full decoupling of conduction electrons and localized
f spins. We complete our analysis by considering the limit of a vanishing Heisenberg coupling
Jh=0, and we provide evidence that the enlarged
SU(N) symmetry and low dimensionality of the Ruderman-Kittel-Kasuya-Yosida exchange interaction are not sufficient to induce VBS order in the conventional
SU(N) Kondo chain, which hosts solely a KI phase. Finally, we quantify the role of increasing
N in reducing the separation of charge and spin timescales specific to the KI state, and we show that the latter is adiabatically connected to a trivial band insulator captured by the large-
N approximation.
To examine the validity of the scenario of the deconfined critical phenomena, we carry out a quantum Monte Carlo simulation for the SU(
N) generalization of the Heisenberg model with four-body and six-body interactions. The quantum phase transition between the SU(
N) N\'eel and valence-bond solid phases is characterized for
N=2, 3, and 4 on the square and honeycomb lattices. While finite-size scaling analysis works well up to the maximum lattice size (
L=256) and indicates the continuous nature of the phase transition, a clear systematic change towards the first-order transition is observed in the estimates of the critical exponent $y\ensuremath{\equiv}1/\ensuremath{\nu}$ as the system size increases. We also confirm the relevance of a squared valence-bond solid field ${\ensuremath{\Psi}}^{2}$ for the SU(3) model.
Motivated by recent experimental progress in the context of ultra-cold multi-colour fermionic atoms in optical lattices, this thesis investigates the properties of the antiferromagnetic SU(N) Heisenberg models with fully antisymmetric irreducible representations labelled by a Young tableau with
m boxes on each site of various lattices. Thanks to an intensive use of variational Monte Carlo method based on Gutzwiller projected fermionic wave functions, we aim to determine the ground state and its related properties of systems which are beyond the capabilities of other numerical methods. This typically means that we will study large lattices with many particles per site for
N>2. The first part of this thesis will introduce the basic theoretical tools connected to the SU(N) Heisenberg model. A very short review of group theory introducing Young tableaux, will clarify the concept of fully antisymmetric irreducible representation of the SU(N) Lie algebra. Some properties of the Heisenberg model will be discussed and a mean-field solution will be given. The numerical methods used during this thesis and their limitations wil also be introduced. We will investigate the properties of 1-dimensional systems in the second part of this work. In particular, the chains have been studied for arbitrary
N and
m with non-abelian bosonisation leading to predictions about the nature of the ground state (gapped or critical) in most but not all cases. We have been able to verify these predictions for a representative number of cases with
N≤10 and
m≤N/2, and we have shown that the opening of a gap is associated with a spontaneous dimerisation or trimerisation depending on the value of
m and
N. We have also investigated the marginal cases, where non-abelian bosonisation did not lead to any prediction. In these cases, variational Monte Carlo predicts that the ground state is critical with exponents that are consistent with conformal field theory. The other 1-dimensional system that has been studied is the ladder for which field-theory gives some predictions. For instance, the ground state of the SU(4) case in the fundamental irreducible representation is believed to consist of a 4-site plaquette state for any non-zero inter-leg coupling. We have confirmed the presence of this state, even if in the weak coupling regime, the optimal variational state is gapless. The study has been extended to various
N and
m values providing interesting results. The last part will present the results for 2-dimensional lattices. The square lattice will first be investigated. A mean-field analysis predicts a plaquette ground state for
3≤N/m<5 and of a chiral ground state for
N/m≥5. By systematically studying the SU(3m) and SU(6m) Heisenberg model, we determined that the mean-field solution is valid even at small
m (
m>1 for SU(3m) and
m>2 for SU(6m)). Then, for the SU(6m) Heisenberg model with
m particles per site on a honeycomb lattice, we make the connection between the
m=1 case, for which the ground state has recently been shown to be in a plaquette singlet state, and the
m→∞ limit, where a mean-field approach has established that the ground state has chiral order. We were able to show that this system has a clear tendency toward the chiral order as soon as
m≥2. This demonstrates that the chiral phase can be stabilised for not too large values of
m, opening the way to its experimental realisations in other lattices. We will finally present results for the SU(3m) on the triangular lattice for which we find a plaquette state, when
m>1.
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).