
We provide numerical evidence in favor of spontaneous chiral symmetry
breaking in the Mott phase of threecomponent fermions on the triangular
lattice, that we describe by an SU(3) symmetric Hubbard model with hopping t
and onsite interactions U. Our approach relies on effective models derived in
the strongcoupling limit in powers of t/U for general SU(N) and arbitrary
flux, which are studied using exact diagonalization and variational Monte Carlo
simulations for N=3. Up to third order in t/U, there is a large chiral phase
that encompasses a range of real ring exchange where timereversal is
spontaneously broken. For the Hubbard model, this phase is stabilized below
U/t~13, as revealed by including higherorder effects. How to realize it with
ultra cold atoms without artificial gauge fields is briefly discussed.

Using determinant quantum Monte Carlo (DQMC) simulations, we systematically
study the doping dependence of the crossover from one to two dimensions and its
impact on the magnetic properties of the Hubbard model. A square lattice of
chains is used, in which the dimensionality can be tuned by varying the
interchain coupling $t_\perp$. The dynamical spin structure factor and static
quantities, such as the static spin susceptibility and nearestneighbor spin
correlation function, are characterized in the one and twodimensional limits
as a benchmark. When the dimensionality is tuned between these limits, the
magnetic properties, while evolving smoothly from one to two dimensions,
drastically change regardless of the doping level. This suggests that the spin
excitations in the twodimensional Hubbard model, even in the heavily doped
case, cannot be explained using the spinon picture known from one dimension.
The DQMC calculations are complemented by cluster perturbation theory studies
to form a more complete picture of how the crossover occurs as a function of
doping and how doped holes impact magnetic order.

We investigate manifestations of topological order in monoaxial helimagnet
Cr$_{1/3}$NbS$_2$ by performing transport measurements on ultrathin crystals.
Upon sweeping the magnetic field perpendicularly to the helical axis, crystals
thicker than one helix pitch (48 nm) but much thinner than the magnetic domain
size ($\sim$1 $\mu$m) are found to exhibit sharp and hysteretic resistance
jumps. We show that these phenomena originate from transitions between
topological sectors with different number of magnetic solitons. This is
confirmed by measurements on crystals thinner than 48 nm in which the
topological sector cannot change that do not exhibit any jump or hysteresis.
Our results show the ability to deterministically control the topological
sector of finitesize Cr$_{1/3}$NbS$_2$ and to detect intersector transitions
by transport measurements.

Quantum antiferromagnets have proven to be some of the cleanest realizations
available for theoretical, numerical, and experimental studies of quantum
fluctuation effects. At finite temperatures, however, the additional effects of
thermal fluctuations in the restricted phase space of a lowdimensional system
have received much less attention, particularly the situation in frustrated
quantum magnets, where the excitations may be complex collective (bound or even
fractionalized) modes. We investigate this problem by studying the
thermodynamic properties of the frustrated twoleg S=1/2 spin ladder, with
particular emphasis on the fully frustrated case. We present numerical results
for the magnetic specific heat and susceptibility, obtained from exact
diagonalization and quantum Monte Carlo studies, which we show can be rendered
free of the sign problem even in a strongly frustrated system and which allow
us to reach unprecedented sizes of L=200 ladder rungs. We find that frustration
effects cause an unconventional evolution of the thermodynamic response across
the full parameter regime of the model. However, close to the firstorder
transition they cause a highly anomalous reduction in temperature scales with
no concomitant changes in the gap; the specific heat shows a very narrow peak
at very low energies and the susceptibility rises abruptly at extremely low
temperatures. Unusually, the two quantities have different gaps over an
extended region of the parameter space. We demonstrate that these results
reflect the presence of large numbers of multiparticle boundstate
excitations, whose energies fall below the onetriplon gap in the transition
region.

We show that the critical fields of the magnetization plateaus of the
ShastrySutherland model decrease significantly upon increasing the ratio of
inter to intradimer coupling, and accordingly that the magnetization plateaus
of SrCu_2(BO_3)_2 shift to lower field under pressure, making the first two
plateaus at 1/8 and 2/15 potentially accessible to neutron scattering
experiments. These conclusions are based on the derivation of an effective
classical model of interacting pinwheelshaped spin2 bound states using a
combination of perturbative and graphbased continuous unitary transformations,
showing that pinwheel crystals are indeed the lowestenergy plateau structures
at low magnetization, and that a simple model of intermediaterange twobody
repulsion between pinwheels is able to account quantitatively for the plateau
sequence.

Using quantum Monte Carlo simulations along with higherorder spinwave
theory, bondoperator and strongcoupling expansions, we analyse the dynamical
spin structure factor of the spinhalf Heisenberg model on the squarelattice
bilayer. We identify distinct contributions from the lowenergy Goldstone modes
in the magnetically ordered phase and the gapped triplon modes in the quantum
disordered phase. In the antisymmetric (with respect to layer inversion)
channel, the dynamical spin structure factor exhibits a continuous evolution of
spectral features across the quantum phase transition, connecting the two types
of modes. Instead, in the symmetric channel we find a depletion of the spectral
weight when moving from the ordered to the disordered phase. While the
dynamical spin structure factor does not exhibit a welldefined distinct
contribution from the amplitude (or Higgs) mode in the ordered phase, we
identify an only marginallydamped amplitude mode in the dynamical singlet
structure factor, obtained from interlayer bond correlations, in the vicinity
of the quantum critical point. These findings provide quantitative information
in direct relation to possible neutron or light scattering experiments in a
fundamental twodimensional quantumcritical spin system.

The magnetic excitation spectrum in the bilayer iridate Sr$_3$Ir$_2$O$_7$ has
been investigated using highresolution resonant inelastic xray scattering
(RIXS) performed at the iridium L$_3$ edge and theoretical techniques. A study
of the systematic dependence of the RIXS spectrum on the orientation of the
wavevector transfer, $\mathbf{Q}$, with respect to the iridiumoxide bilayer
has revealed that the magnon dispersion is comprised of two branches well
separated in energy and gapped across the entire Brillouin zone. Our results
contrast with those of an earlier study which reported the existence of a
single dominant branch. While these earlier results were interpreted as two
overlapping modes within a spinwave model of weakly coupled iridiumoxide
planes, our results are more reminiscent of those expected for a system of
weakly coupled dimers. In this latter approach the lower and higher energy
modes find a natural explanation as those corresponding to transverse and
longitudinal fluctuations, respectively. We have therefore developed a
bondoperator theory which describes the magnetic dispersion in
Sr$_3$Ir$_2$O$_7$ in terms of quantum dimer excitations. In our model
dimerisation is produced by the leading Heisenberg exchange, $J_c$, which
couples iridium ions in adjacent planes of the bilayer. The Hamiltonian also
includes in plane exchange, $J$, as well as further neighbour couplings and
relevant anisotropies. The bondoperator theory provides an excellent account
of the dispersion of both modes, while the measured $\mathbf{Q}$ dependence of
the RIXS intensities is in reasonable qualitative accord with the spinspin
correlation function calculated from the theory. We discuss our results in the
context of the quantum criticality of bilayer dimer systems in the presence of
anisotropic interactions derived from strong spinorbit coupling.

We report on the real space profile of spin polarons in the quasi
twodimensional frustrated dimer spin system SrCu2(BO3)2 doped with 0.16% of
Zn. The 11B nuclear magnetic resonance spectrum exhibits 15 additional boron
sites near nonmagnetic Zn impurities. With the help of exact diagonalizations
of finite clusters, we have deduced from the boron spectrum the distribution of
local magnetizations at the Cu sites with fine spatial resolution, providing
direct evidence for an extended spin polaron. The results are confronted with
those of other experiments performed on doped and undoped samples of
SrCu2(BO3)2.

We investigate the structure of the spectrum of antiferromagnetically coupled
spin1 bosons on a square lattice using degenerate perturbation theory and
exact diagonalizations of finite clusters. We show that the superfluid phase
develops an Anderson tower of states typical of nematic longrange order with
broken SU(2) symmetry.We further show that this order persists into the Mott
insulating phase down to zero hopping for one boson per site, and down to a
critical hopping for two bosons per site, in agreement with meanfield and
Quantum Monte Carlo results. The connection with the transition between a
fragmented condensate and a polar one in a single trap is briefly discussed.

The magnetization process of the orthogonaldimer antiferromagnet SrCu2(BO3)2
is investigated in high magnetic fields of up to 118 T. A 1/2 plateau is
clearly observed in the field range 84 to 108 T in addition to 1/8, 1/4 and 1/3
plateaux at lower fields. Using a combination of stateoftheart numerical
simulations, the main features of the highfield magnetization, a 1/2 plateau
of width 24 T, a 1/3 plateau of width 34 T, and no 2/5 plateau, are shown to
agree quantitatively with the ShastrySutherland model if the ratio of inter
to intradimer exchange interactions J'/J=0.63. It is further predicted that
the intermediate phase between the 1/3 and 1/2 plateau is not uniform but
consists of a 1/3 supersolid followed by a 2/5 supersolid and possibly a
domainwall phase, with a reentrance into the 1/3 supersolid above the 1/2
plateau.

We study the spin S=1 antiferromagnetic Heisenberg model on the square
lattice with, in addition to the nearestneighbor interaction, a threesite
interaction of the form (S_i S_j)*(S_j S_k) + h.c.. This interaction appears
naturally in a strong coupling exansion of the twoorbital, halffilled Hubbard
model. For spin 1/2, this model reduces to a Heisenberg model with bilinear
interactions up to third neighbors, with a secondneighbor interaction twice as
large the thirdneighbor one, a very frustrated model with an infinite family
of helical classical ground states in a large parameter range. Using a variety
of analytical and numerical methods, we show that the spin1 case is also very
frustrated, and that its phase diagram is even richer, with possibly the
succession of seven different phases as a function of the ratio of the
threesite interaction to the bilinear one. The phases are either purely
magnetic phases with collinear order, or of mixed magnetic and quadrupolar
character with helical order.

In order to understand the nature of the twodimensional BoseEinstein
condensed (BEC) phase in BaCuSi2O6, we performed detailed 63Cu and 29Si NMR
above the critical magnetic field, Hc1= 23.4 T. The two different alternating
layers present in the system have very different local magnetizations close to
Hc1; one is very weak, and its size and field dependence are highly sensitive
to the nature of interlayer coupling. Its precise value could only be
determined by "onsite" 63Cu NMR, and the data are fully reproduced by a model
of interacting hardcore bosons in which the perfect frustration associated to
tetragonal symmetry is slightly lifted, leading to the conclusion that the
population of the less populated layers is not fully incoherent but must be
partially condensed.

Motivated by the intriguing report, in some frustrated quantum
antiferromagnets, of magnetization plateaus whose simple collinear structure is
{\it not} stabilized by an external magnetic field in the classical limit, we
develop a semiclassical method to estimate the zeropoint energy of collinear
configurations even when they do not correspond to a local minimum of the
classical energy. For the spin1/2 frustrated squarelattice antiferromagnet,
this approach leads to the stabilization of a large 1/2 plateau with
"upupupdown" structure for J_2/J_1>1/2, in agreement with exact
diagonalization results, while for the spin1/2 anisotropic triangular
antiferromagnet, it predicts that the 1/3 plateau with "upupdown" structure
is stable far from the isotropic point, in agreement with the properties of
Cs_2CuBr_4.

In addition to lowenergy spin fluctuations, which distinguish them from band
insulators, Mott insulators often possess orbital degrees of freedom when
crystalfield levels are partially filled. While in most situations spins and
orbitals develop longrange order, the possibility for the ground state to be a
quantum liquid opens new perspectives. In this paper, we provide clear evidence
that the SU(4) symmetric KugelKhomskii model on the honeycomb lattice is a
quantum spinorbital liquid. The absence of any form of symmetry breaking 
lattice or SU(N)  is supported by a combination of semiclassical and numerical
approaches: flavorwave theory, tensor network algorithm, and exact
diagonalizations. In addition, all properties revealed by these methods are
very accurately accounted for by a projected variational wavefunction based on
the \piflux state of fermions on the honeycomb lattice at 1/4filling. In that
state, correlations are algebraic because of the presence of a Dirac point at
the Fermi level, suggesting that the symmetric KugelKhomskii model on the
honeycomb lattice is an algebraic quantum spinorbital liquid. This model
provides a good starting point to understand the recently discovered
spinorbital liquid behavior of Ba_3CuSb_2O_9. The present results also suggest
to choose optical lattices with honeycomb geometry in the search for quantum
liquids in ultracold fourcolor fermionic atoms.

We investigate the zerotemperature behavior of the classical Heisenberg
model on the triangular lattice in which the competition between exchange
interactions of different orders favors a relative angle between neighboring
spins in the interval (0,2pi/3). In this situation, the ground states are
noncoplanar and have an infinite discrete degeneracy. In the generic case, the
set of the ground states is in one to one correspondence (up to a global
rotation) with the noncrossing loop coverings of the three equivalent
honeycomb sublattices into which the bonds of the triangular lattice can be
partitioned. This allows one to identify the order parameter space as an
infinite Cayley tree with coordination number 3. Building on the duality
between a similar loop model and the ferromagnetic O(3) model on the honeycomb
lattice, we argue that a typical ground state should have longrange order in
terms of spin orientation. This conclusion is further supported by the
comparison with the fourstate antiferromagnetic Potts model [describing the
case when the angle between neighboring spins is equal to arccos(1/3)], which
at zero temperature is critical and in terms of the solidonsolid
representation is located exactly at the point of roughening transition. At
other values of the angle between neighboring spins an additional constraint
appears, whose presence drives the system into an ordered phase (unless this
angle is equal to pi/2, when another constraint is removed and the model
becomes trivially exactly solvable).

A small value of the spin gap in quantum antiferromagnets with strong
frustration makes them susceptible to nominally small deviations from the ideal
Heisenberg model. One of such perturbations, the anisotropic
DzyaloshinskiiMoriya interaction, is an important perturbation for the $S=1/2$
kagome antiferromagnet, one of the current candidates for a quantumdisordered
ground state. We study the influence of the DM term in a related
onedimensional system, the sawtooth chain that has valencebond order in its
ground state. Through a combination of analytical and numerical methods, we
show that a relatively weak DM coupling, $0.115J$, is sufficient to destroy the
valencebond order, close the spin gap, and turn the system into a Luttinger
liquid with algebraic spin correlations. A similar mechanism may be at work in
the kagome antiferromagnet.

Motivated by the ongoing investigation of SrCu$_2$(BO$_3$)$_2$ under
pressure, we study a variant of the twodimensional ShastrySutherland (SS)
spin1/2 model with two types of dimers. Combined with the frustration of the
SS model, this modification induces, in a large parameter range, a dimensional
reduction at low energies, with nearly decoupled effective S=1 Haldane chains
forming along one of the diagonals of the lattice. We also present evidence
that the intermediate plaquette solid phase of the undistorted SS model remains
stable in a finite region of the phase diagram.

We show that the spin liquid phase of the halffilled Hubbard model on the
triangular lattice can be described by a pure spin model. This is based on a
highorder strong coupling expansion (up to order 12) using perturbative
continuous unitary transformations. The resulting spin model is consistent with
a transition from threesublattice longrange magnetic order to an insulating
spin liquid phase, and with a jump of the double occupancy at the transition.
Exact diagonalizations of both models show that the effective spin model is
quantitatively accurate well into the spin liquid phase, and a comparison with
the Gutzwiller projected Fermi sea suggests a gapless spectrum and a spinon
Fermi surface.

Motivated by the current interest in the quantum dimer model on the
triangular lattice, we investigate the phase diagram of the closely related
fullyfrustrated transverse field Ising model on the honeycomb lattice using
classical and semiclassical approximations. We show that, in addition to the
fully polarized phase at large field, the classical model possesses a multitude
of phases that break the translational symmetry which in the dimer language,
correspond to a plaquette phase and a columnar phase separated by an infinite
cascade of mixed phases. The modification of the phase diagram by quantum
fluctuations has been investigated in the context of linear spinwave theory.
The extrapolation of the semiclassical energies suggests that the plaquette
phase extends down to zero field for spin 1/2, in agreement with the
$\sqrt{12}\times\sqrt{12}$ phase of the quantum dimer model on the triangular
lattice with only kinetic energy.

Magnetic excitations in copper pyrimidine dinitrate, a spin1/2
antiferromagnetic chain with alternating $g$tensor and DzyaloshinskiiMoriya
interactions that exhibits a fieldinduced spin gap, are probed by means of
pulsedfield electron spin resonance spectroscopy. In particular, we report on
a minimum of the gap in the vicinity of the saturation field $H_{sat}=48.5$ T
associated with a transition from the sineGordon region (with solitonbreather
elementary excitations) to a spinpolarized state (with magnon excitations).
This interpretation is fully confirmed by the quantitative agreement over the
entire field range of the experimental data with the DMRG investigation of the
spin1/2 Heisenberg chain with a staggered transverse field.

We present a thorough $^{77}$Se NMR study of a single crystal of the
magnetoelectric compound Cu$_2$OSeO$_3$. The temperature dependence of the
local electronic moments extracted from the NMR data is fully consistent with a
magnetic phase transition from the highT paramagnetic phase to a lowT
ferrimagnetic state with 3/4 of the Cu$^{2+}$ ions aligned parallel and 1/4
aligned antiparallel to the applied field of 14.09 T. The transition to this
3up1down magnetic state is not accompanied by any splitting of the NMR lines
or any abrupt modification in their broadening, hence there is no observable
reduction of the crystalline symmetry from its highT cubic \textit{P}2$_1$3
space group. These results are in agreement with high resolution xray
diffraction and magnetization data on powder samples reported previously by Bos
{\it et al.} [Phys. Rev. B, {\bf 78}, 094416 (2008)]. We also develop a mean
field theory description of the problem based on a microscopic spin Hamiltonian
with one antiferromagnetic ($J_\text{afm}\simeq 68$ K) and one ferromagnetic
($J_\text{fm}\simeq 50$ K) nearestneighbor exchange interaction.

When looking for analytical approaches to treat frustrated quantum magnets,
it is often very useful to start from a limit where the ground state is highly
degenerate. This chapter discusses several ways of deriving {effective
Hamiltonians} around such limits, starting from standard {degenerate
perturbation theory} and proceeding to modern approaches more appropriate for
the derivation of highorder effective Hamiltonians, such as the perturbative
continuous unitary transformations or contractor renormalization. In the course
of this exposition, a number of examples taken from the recent literature are
discussed, including frustrated ladders and other dimerbased Heisenberg models
in a field, as well as the mapping between frustrated Ising models in a
transverse field and quantum dimer models.

We study the T = 0 magnetization of frustrated twoleg spin ladders with
arbitrary value of the spin S. In the strong rung limit, we use degenerate
perturbation theory to prove that frustration leads to magnetization plateaux
at fractional values of the magnetization for all spins S, and to determine the
critical ratios of parallel to diagonal interrung couplings for the appearance
of these plateaux. These ratios depend both on the plateau and on the spin. To
confirm these results, and to investigate the properties of these ladders away
from the strong coupling limit, we have performed extensive density matrix
renormalization group (DMRG) calculations for S <= 2. For large enough
interrung couplings, all plateaux simply disappear, leading to a magnetization
curve typical of integerspin chains in a magnetic field. The intermediate
region turns out to be surprisingly rich however, with, upon increasing the
interrung couplings, the development of magnetization jumps and, in some
cases, the appearance of one or more phase transitions inside a given plateau.

Controlled doping of SrCu$_2$(BO$_3$)$_2$, a faithful realization of the
Heisenberg spin1/2 antiferromagnet on the ShastrySutherland lattice, with
nonmagnetic impurities generates boundstates below the spin gap. These
boundstates and their symmetry properties are investigated by exact
diagonalisation of small clusters and within a simple effective model
describing a spinon submitted to an attractive extended potential. It is shown
that Raman spectroscopy is a unique technique to probe these boundstates.
Quantitative theoretical Raman spectra are numerically obtained.

We discuss the appearance of supersolid phases for interacting hardcore
bosons on the square lattice when, in addition to the standard nearest neighbor
hopping and repulsion, correlated or nextnearest neighbor hopping is present.
Having in mind dimerbased quantum magnets in a field described by effective
bosonic models of this kind, we put special emphasis on a comparison between
the different cases of relative signs of the kinetic processes, which
correspond to unfrustrated or frustrated magnetic models. In the unfrustrated
case, we compare Quantum Monte Carlo simulations with a meanfield (classical)
approach, which is shown to give qualitatively correct results. Using this
classical approach for the frustrated case, we find that the phase diagram is
generically richer than in the unfrustrated case. We also investigate in detail
the differences between standard nextnearest neighbour and correlated hopping
over the diagonal, with the conclusion that both cases are similar if
checkerboard order is present at halffilling, while a supersolid phase can be
stabilized without any adjacent solid phase only in the case of correlated
hopping.