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.
Motivated by recent experimental progress in the context of ultracold multicolor fermionic atoms in optical lattices, we have developed a method to exactly diagonalize the Heisenberg
SU(N) Hamiltonian with several particles per site living in a fully symmetric or antisymmetric representation of
SU(N). The method, based on the use of standard Young tableaux, takes advantage of the full
SU(N) symmetry, allowing one to work directly in each irreducible representation of the global
SU(N) group. Since the
SU(N) singlet sector is often much smaller than the full Hilbert space, this enables one to reach much larger system sizes than with conventional exact diagonalizations. The method is applied to the study of Heisenberg chains in the symmetric representation with two and three particles per site up to
N=10 and up to 20 sites. For the length scales accessible to this approach, all systems except the Haldane chain [
SU(2) with two particles per site] appear to be gapless, and the central charge and scaling dimensions extracted from the results are consistent with a critical behavior in the
SU(N) level
k Wess-Zumino-Witten universality class, where
k is the number of particles per site. These results point to the existence of a crossover between this universality class and the asymptotic low-energy behavior with a gapped spectrum or a critical behavior in the
SU(N) level 1 WZW universality class.
SU ( N ) quantum spin systems may be realized in a variety of physical systems including ultracold atoms in optical lattices. The study of such models also leads to insights into possible novel quantum phases and phase transitions of SU (2) spin models. Here we use Gutzwiller projected fermionic variational wavefunctions to explore the phase diagram and correlation functions of SU ( N ) quantum spin models in the self-conjugate representation. In one dimension, the ground state of the SU (4) spin chain with Heisenberg bilinear and biquadratic interactions is studied by examining instabilities of the Gutzwiller projected free fermion ground state to various broken symmetries. The variational phase diagram so obtained agrees well with exact results. The spin–spin and dimer–dimer correlation functions of the Gutzwiller projected free fermion state with N flavours of fermions are in good agreement with exact and 1/ N calculations for the critical points of SU ( N ) spin chains. In two dimensions, the phase diagram of the antiferromagnetic Heisenberg model on the square lattice is obtained by finding instabilities of the Gutzwiller projected π-flux state. In the absence of biquadratic interactions, the model exhibits long-range Néel order for N = 2 and 4, and spin Peierls (columnar dimer) order for N >4. Upon including biquadratic interactions in the SU (4) model (with a sign appropriate to a fermionic Hubbard model), the Néel order diminishes and eventually disappears, giving way to an extended valence bond crystal. In the case of the SU (6) model, the dimerized ground state melts at sufficiently large biquadratic interaction, yielding a projected π-flux spin liquid phase which in turn undergoes a transition into an extended valence bond crystal at even larger biquadratic interaction. The spin correlations of the projected π-flux state at N = 4 are in good agreement with 1/ N calculations. We find that the state shows strongly enhanced dimer correlations, in qualitative agreement with recent theoretical predictions. We also compare our results with a recent quantum Monte Carlo study of the SU (4) Heisenberg model.
We compute the dynamical spin structure factor $S(k,\ensuremath{\omega})$ of the SU(3) Heisenberg chain variationally using a truncated Hilbert space spanned by the Gutzwiller projected particle-hole excitations of the Fermi sea, introduced by Dalla Piazza et al. [Nat. Phys. 11, 62 (2015)], with a modified importance sampling. We check the reliability of the method by comparing the $S(k,\ensuremath{\omega})$ to exact diagonalization results for 18 sites and to the two-soliton continuum of the Bethe ansatz for 72 sites. We get an excellent agreement in both cases. Detailed analysis of the finite-size effects shows that the method captures the critical Wess-Zumino-Witten SU(3)
1 behavior and reproduces the correct exponent, with the exception of the size dependence of the weight of the bottom of the conformal tower. We also calculate the single-mode approximation for the
SU(N) Heisenberg model and determine the velocity of excitations. Finally, we apply the method to the SU(3) Haldane-Shastry model and find that the variational method gives the exact wave function for the lowest excitation at $k=\ifmmode\pm\else\textpm\fi{}2\ensuremath{\pi}/3$.
By using Gutzwiller projected fermionic wave functions and variational Monte Carlo technique, we study the spin-
1/2 Heisenberg model with the first-neighbor (
J1), second-neighbor (
J2), and additional scalar chiral interaction ${J}_{\ensuremath{\chi}}{\mathbf{S}}_{i}\ifmmode\cdot\else\textperiodcentered\fi{}({\mathbf{S}}_{j}\ifmmode\times\else\texttimes\fi{}{\mathbf{S}}_{k})$ on the triangular lattice. In the nonmagnetic phase of the ${J}_{1}\ensuremath{-}{J}_{2}$ triangular model with $0.08\ensuremath{\lesssim}{J}_{2}/{J}_{1}\ensuremath{\lesssim}0.16$, recent density-matrix renormalization group (DMRG) studies [Zhu and White, Phys. Rev. B 92, 041105(R) (2015) and Hu, Gong, Zhu, and Sheng, Phys. Rev. B 92, 140403(R) (2015)] find a possible gapped spin liquid with the signal of a competition between a chiral and a
Z2 spin liquid. Motivated by the DMRG results, we consider the chiral interaction ${J}_{\ensuremath{\chi}}{\mathbf{S}}_{i}\ifmmode\cdot\else\textperiodcentered\fi{}({\mathbf{S}}_{j}\ifmmode\times\else\texttimes\fi{}{\mathbf{S}}_{k})$ as a perturbation for this nonmagnetic phase. We find that with growing ${J}_{\ensuremath{\chi}}$, the gapless U(1) Dirac spin liquid, which has the best variational energy for ${J}_{\ensuremath{\chi}}=0$, exhibits the energy instability towards a gapped spin liquid with nontrivial magnetic fluxes and nonzero chiral order. We calculate topological Chern number and ground-state degeneracy, both of which identify this flux state as the chiral spin liquid with fractionalized Chern number
C=1/2 and twofold topological degeneracy. Our results indicate a positive direction to stabilize a chiral spin liquid near the nonmagnetic phase of the ${J}_{1}\ensuremath{-}{J}_{2}$ triangular model.
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.