This thesis explores various approaches of studying the long-range colour order of antiferromagnetic SU(N) Heisenberg models with the linear flavour-wave theory (LFWT). The LFWT is an extension of the well-known SU(2) spin-wave theory to SU(N), and this semi-classical method has been used to study SU(N) models with colour-ordered ground states in the fundamental irreducible representation (irrep). The aim of this thesis is to study various SU(N) Heisenberg models in different irreps and in various colour configurations, in part by extending the LFWT method to different irreps of SU(N). This will be achieved by using three different methods all using different bosonic representations. First, the LFWT willbe performed using the multiboson method that introduces a boson for each state of a given irrep. Then, a different SU(3) bosonic representation first introduced by Mathur & Sen will be presented and used to perform the LFWT calculations. Finally, another way of applying the LFWT will be shown using the bosons first used by Read & Sachdev for SU(N) irreps with rectangular Young tableaux. The specific models that will be treated by these methods are the fully antisymmetric SU(N) models on the square, honeycomb, and triangular lattices with m>1 particles per site, as well as the bipartite SU(3) Heisenberg chain in the adjoint irrep to show how the LFWT can be used for mixed symmetries. The last chapter of the thesis will be dedicated to another problem related to the accidental line of zero modes in the harmonic spectrum of the antiferromagnetic SU(3) Heisenberg model on the square lattice with one particle per site in a three-sublattice order. The objective will be to try to lift the accidental zero modes in the spectrum.
Using linear flavor-wave theory (LFWT) and auxiliary field quantum Monte Carlo (QMC), we investigate the properties of the SU(4) Heisenberg model on the anisotropic square lattice in the fully antisymmetric six-dimensional irreducible representation, a model that describes interacting fermions with four flavors at half-filling. Thanks to the calculations on very large systems, we have been able to convincingly demonstrate that QMC results are consistent with a small but finite antiferromagnetic moment at the isotropic point, in qualitative agreement with LFWT obtained earlier [F. H. Kim et al., Phys. Rev. B 96, 205142 (2017)], and in quantitative agreement with results obtained previously on the Hubbard model [D. Wang et al., Phys. Rev. Lett. 112, 156403 (2014)] after extrapolation to infinite
U/t. The presence of a long-range antiferromagnetic order has been further confirmed by showing that a phase transition takes place into a valence-bond solid (VBS) phase not too far from the isotropic point when reducing the coupling constant along one direction on the way to decoupled chains.
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 first show that, for the SU(N) Heisenberg model with isotropic two-site interactions on an arbitrary lattice of Ns sites, the full spectrum in its singlet sector (when N divides Ns) is the negative of the full spectrum in its singlet sector of the SU(Ns/N) Heisenberg model defined on the same lattice. This theorem, valid for models with the fundamental irreducible representation (irrep) at each lattice site, follows naturally from a transposition operation in the basis of Standard Young Tableaux (SYTs). Second, this basis -- previously recognized as particularly convenient for exact diagonalization (ED) of SU(N) Heisenberg models -- admits an efficient, memory-saving implementation, thereby enabling large-scale diagonalizations for unprecedented system sizes Ns and numbers of colors N>3. We apply these ideas to an ED study of the SU(4) and SU(7) Heisenberg models on the triangular lattice (TL), extending the calculations up to Ns=28 sites. In particular, starting from the SU(4) case, whose nearest-neighbor ground state remains debated, we confirm that next-nearest-neighbor couplings stabilize four-sublattice long-range order. We further find that increasing the number of colors from N=4 to N=7 provides another route toward a Néel-like ground state, even without such additional couplings. N = 4 to N = 7 provides another route toward a Néel-like ground state, even without such additional couplings.
We show how to perform exact diagonalizations of SU(N) Fermi-Hubbard models on L-site clusters separately in each irreducible representation (irrep) of SU(N). Using the representation theory of the unitary group U(L), we demonstrate that a convenient orthonormal basis, on which matrix elements of the Hamiltonian are very simple, is given by the set of semistandard Young tableaux (or, equivalently, the Gelfand-Tsetlin patterns) corresponding to the targeted irrep. As an application of this color factorization, we study the robustness of some SU(N) phases predicted in the Heisenberg limit upon decreasing the on-site interaction U on various lattices of size L≤12 and for 2≤N≤6. In particular, we show that a long-range color ordered phase emerges for intermediate U for N=4 at filling 1/4 on the triangular lattice.
We show how to perform exact diagonalizations of
SU(N) Fermi-Hubbard models on
L-site clusters separately in each irreducible representation ({irrep}) of
SU(N). Using the representation theory of the unitary group
U(L), we demonstrate that a convenient orthonormal basis, on which matrix elements of the Hamiltonian are very simple, is given by the set of {\it semistandard Young tableaux} (or, equivalently the Gelfand-Tsetlin patterns) corresponding to the targeted irrep. As an application of this color factorization, we study the robustness of some
SU(N) phases predicted in the Heisenberg limit upon decreasing the on-site interaction
U on various lattices of size
L≤12 and for
2≤N≤6. In particular, we show that a long-range color ordered phase emerges for intermediate
U for
N=4 at filling
1/4 on the triangular lattice.