We provide numerical evidence in favor of spontaneous chiral symmetry breaking in the Mott phase of three-component fermions on the triangular lattice, that we describe by an SU(3) symmetric Hubbard model with hopping t and on-site interactions U. Our approach relies on effective models derived in the strong-coupling 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 time-reversal is spontaneously broken. For the Hubbard model, this phase is stabilized below U/t~13, as revealed by including higher-order 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 nearest-neighbor spin correlation function, are characterized in the one- and two-dimensional 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 two-dimensional 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 ultra-thin 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 finite-size Cr$_{1/3}$NbS$_2$ and to detect inter-sector 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 low-dimensional 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 two-leg 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 first-order 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 multi-particle bound-state excitations, whose energies fall below the one-triplon gap in the transition region.

We show that the critical fields of the magnetization plateaus of the Shastry-Sutherland model decrease significantly upon increasing the ratio of inter- to intra-dimer 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 pinwheel-shaped spin-2 bound states using a combination of perturbative and graph-based continuous unitary transformations, showing that pinwheel crystals are indeed the lowest-energy plateau structures at low magnetization, and that a simple model of intermediate-range two-body repulsion between pinwheels is able to account quantitatively for the plateau sequence.

Using quantum Monte Carlo simulations along with higher-order spin-wave theory, bond-operator and strong-coupling expansions, we analyse the dynamical spin structure factor of the spin-half Heisenberg model on the square-lattice bilayer. We identify distinct contributions from the low-energy 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 well-defined distinct contribution from the amplitude (or Higgs) mode in the ordered phase, we identify an only marginally-damped 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 two-dimensional quantum-critical spin system.

The magnetic excitation spectrum in the bilayer iridate Sr$_3$Ir$_2$O$_7$ has been investigated using high-resolution resonant inelastic x-ray 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 iridium-oxide 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 spin-wave model of weakly coupled iridium-oxide 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 bond-operator 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 bond-operator 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 spin-spin 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 spin-orbit coupling.

We report on the real space profile of spin polarons in the quasi two-dimensional 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 non-magnetic 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 spin-1 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 long-range 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 mean-field 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 orthogonal-dimer 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 state-of-the-art numerical simulations, the main features of the high-field 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 Shastry-Sutherland model if the ratio of inter- to intra-dimer 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 domain-wall 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 nearest-neighbor interaction, a three-site 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 two-orbital, half-filled Hubbard model. For spin 1/2, this model reduces to a Heisenberg model with bilinear interactions up to third neighbors, with a second-neighbor interaction twice as large the third-neighbor 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 spin-1 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 three-site 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 two-dimensional Bose-Einstein 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 inter-layer coupling. Its precise value could only be determined by "on-site" 63Cu NMR, and the data are fully reproduced by a model of interacting hard-core 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 zero-point energy of collinear configurations even when they do not correspond to a local minimum of the classical energy. For the spin-1/2 frustrated square-lattice antiferromagnet, this approach leads to the stabilization of a large 1/2 plateau with "up-up-up-down" structure for J_2/J_1>1/2, in agreement with exact diagonalization results, while for the spin-1/2 anisotropic triangular antiferromagnet, it predicts that the 1/3 plateau with "up-up-down" structure is stable far from the isotropic point, in agreement with the properties of Cs_2CuBr_4.

In addition to low-energy spin fluctuations, which distinguish them from band insulators, Mott insulators often possess orbital degrees of freedom when crystal-field levels are partially filled. While in most situations spins and orbitals develop long-range 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 Kugel-Khomskii model on the honeycomb lattice is a quantum spin-orbital liquid. The absence of any form of symmetry breaking - lattice or SU(N) - is supported by a combination of semiclassical and numerical approaches: flavor-wave theory, tensor network algorithm, and exact diagonalizations. In addition, all properties revealed by these methods are very accurately accounted for by a projected variational wave-function based on the \pi-flux state of fermions on the honeycomb lattice at 1/4-filling. In that state, correlations are algebraic because of the presence of a Dirac point at the Fermi level, suggesting that the symmetric Kugel-Khomskii model on the honeycomb lattice is an algebraic quantum spin-orbital liquid. This model provides a good starting point to understand the recently discovered spin-orbital 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 ultra-cold four-color fermionic atoms.

We investigate the zero-temperature 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 non-crossing 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 long-range order in terms of spin orientation. This conclusion is further supported by the comparison with the four-state 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 solid-on-solid 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 Dzyaloshinskii-Moriya interaction, is an important perturbation for the $S=1/2$ kagome antiferromagnet, one of the current candidates for a quantum-disordered ground state. We study the influence of the DM term in a related one-dimensional system, the sawtooth chain that has valence-bond 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 valence-bond 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 on-going investigation of SrCu$_2$(BO$_3$)$_2$ under pressure, we study a variant of the two-dimensional Shastry-Sutherland (SS) spin-1/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 half-filled Hubbard model on the triangular lattice can be described by a pure spin model. This is based on a high-order strong coupling expansion (up to order 12) using perturbative continuous unitary transformations. The resulting spin model is consistent with a transition from three-sublattice long-range 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 fully-frustrated transverse field Ising model on the honeycomb lattice using classical and semi-classical 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 spin-wave 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 spin-1/2 antiferromagnetic chain with alternating $g$-tensor and Dzyaloshinskii-Moriya interactions that exhibits a field-induced spin gap, are probed by means of pulsed-field 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 sine-Gordon region (with soliton-breather elementary excitations) to a spin-polarized 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 spin-1/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 high-T paramagnetic phase to a low-T 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 3up-1down 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 high-T cubic \textit{P}2$_1$3 space group. These results are in agreement with high resolution x-ray 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) nearest-neighbor 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 high-order 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 dimer-based 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 two-leg 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 inter-rung 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 inter-rung couplings, all plateaux simply disappear, leading to a magnetization curve typical of integer-spin chains in a magnetic field. The intermediate region turns out to be surprisingly rich however, with, upon increasing the inter-rung 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 spin-1/2 antiferromagnet on the Shastry-Sutherland lattice, with non-magnetic impurities generates bound-states below the spin gap. These bound-states 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 bound-states. 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 next-nearest neighbor hopping is present. Having in mind dimer-based 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 mean-field (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 next-nearest neighbour and correlated hopping over the diagonal, with the conclusion that both cases are similar if checkerboard order is present at half-filling, while a supersolid phase can be stabilized without any adjacent solid phase only in the case of correlated hopping.