Quantum spin liquids (QSLs) are long-range entangled phases of frustrated magnets exhibiting fractionalised spin excitations. In two dimensions, there is limited analytical understanding of their excitation spectra beyond parton mean-field theories, which fail to capture many features of the finite frequency dynamical response. We use a self-consistent mean-field theory combined with the random phase approximation (RPA) for the $J_1$-$J_2$ Heisenberg model on the triangular lattice to describe the strong spinon-spinon interactions of the U(1) Dirac QSL. We obtain quantitative results for the dynamical spin structure factor and phase diagram compatible with comprehensive numerical efforts. We show the continuum response at the Brillouin zone corners is described by a spinon-exciton hybridised with monopole gauge excitations. We further show that the transition from the QSL to 120-degree coplanar order is in the QED$_3$ O(3)-GN universality class. We extend the method to chiral QSLs and XXZ anisotropy, and discuss its broad range of applicability to other models and for describing inelastic neutron scattering experiments in KYbSe$_2$.
Dynamics and stability of U(1) spin liquids beyond mean-field theory: Triangular-lattice J1-J2 Heisenberg model
Quantum spin liquids (QSLs) are long-range entangled phases of frustrated magnets exhibiting fractionalised spin excitations. In two dimensions, there is limited analytical understanding of their excitation spectra beyond parton mean-field theories, which fail to capture many features of the finite frequency dynamical response. We use a self-consistent mean-field theory combined with the random phase approximation (RPA) for the J1-J2 Heisenberg model on the triangular lattice to describe the strong spinon-spinon interactions of the U(1) Dirac QSL. We obtain quantitative results for the dynamical spin structure factor and phase diagram compatible with comprehensive numerical efforts. We show the continuum response at the Brillouin zone corners is described by a spinon-exciton hybridised with monopole gauge excitations. We further show that the transition from the QSL to 120-degree coplanar order is in the QED3 O(3)-GN universality class. We extend the method to chiral QSLs and XXZ anisotropy, and discuss its broad range of applicability to other models and for describing inelastic neutron scattering experiments in KYbSe2.
Frustrated magnetism provides a rich platform for realizing exotic quantum phases in condensed matter systems. Among these, quantum spin liquids (QSLs) are a particularly intriguing class of states that exhibit no magnetic ordering and preserve all symmetries down to zero temperature. Just as conventional magnets can be thought of as spin crystals,QSLs represent their liquid analogs. This thesis explores two broad types of QSLs in two-dimensional magnets, with a focus on their dynamic properties and the nature of transitions from conventional magnetic states into QSL phases. QSLs are topological phases that host fractionalized excitations known as spinons. In these states, mobile spinons are coupled to an emergent gauge field, forming a highly entangled phase of matter. The first part of this thesis focuses on Dirac spin liquids (DSLs), a class of QSLs where spinons exhibit a gapless Dirac dispersion. To study the dynamicsof these spin liquids beyond leading-order mean-field theory, we develop a novel large-N fermionic limit of gapless spinons. In this limit, all gauge excitations are suppressed, and the leading-order corrections instead arise from order parameter fluctuations, which may condense and drive transitions into symmetry-broken phases from the parent DSL. Incorporating random phase approximation (RPA) corrections to the spin susceptibility, we find that sharp, spin-wave–like modes are recovered in the excitation spectrum. This methodology successfully captures characteristic features of the dynamics near phase transitions and in weakly symmetry-broken phases proximate to a DSL. We apply this approach to microscopic models, including the J1-J2 XXZ model on the triangular lattice and the J1-J3 XXZ model on the honeycomb lattice. These models closely represent a wide range of frustrated materials, such as triangular and honeycomb cobaltates and triangular ytterbium compounds. Our theoretical predictions show excellent agreement with experimentally measured inelastic neutron scattering spectra from these materials. Gapless spinons can also sometimes become gapped through chiral symmetry breaking, leading to the formation of a chiral spin liquid (CSL). The CSL is an archetypal bosonic realization of a fractional quantum Hall (FQH) state at half-filling, and has been proposed as a time-reversal–breaking ground state in a variety of frustrated spin models. The second part of this thesis addresses the intriguing possibility of realizing a CSL through the "melting" of a chiral spin crystal—a conventional magnet with non-coplanar ordering. Such a transition offers a physically appealing analogue to the melting of a Wigner crystal into a fermionic FQH liquid. Focusing on the kagome lattice, where chiral interactions are knownto stabilize a CSL ground state, we study the effects of additional long-range Heisenberg interactions and demonstrate a phase transition into a non-coplanar octahedral magnetic order. To describe the universal class of transition between a CSL and a non-coplanar anti-ferromagnetic order, we then develop a novel Chern-Simons–matter field theory framework, constrained by symmetry principles and symmetry fractionalization patterns. This thesis thus provides new insights into the dynamical signatures and phase transitions of quantum spin liquids, offering a bridge between theoretical models and experimental observations in frustrated magnetic systems.
Investigation of the J1-J2 Heisenberg model on the triangular lattice: A study with projected entangled-pair states
The nature of the quantum spin liquid (QSL) phase in the frustrated J1-J2 Heisenberg model on the triangular lattice remains an open and actively debated problem. In this work, we employ the infinite projected entangled-pair state (PEPS) to systematically investigate the model under different symmetry constraints. Our simulations reveal a direct transition from the 120∘ Néel state to a putative QSL at J2/J1≈0.08, signaled by the collapse of magnetic order. We further show that, through either an appropriate unitary rotation or spontaneous spin long-range order, the stripe antiferromagnetic phase can also be accurately captured within the infinite PEPS framework. A central focus of our study is the role played by the PEPS symmetry in approximating the QSL ground-state sandwiched between the two magnetic phases. We first found that a fully-symmetric topological Z2 Resonating Valence Bond state, which can be written as a simple PEPS with bond dimension D=3, exhibits a reasonably good variational energy. Motivated by this finding, we have further constructed generic Z2-symmetric PEPS of larger bond dimension (up to D=7). We found that, under wavefunction optimization, spinons condense and, simultaneously, topological vison excitations get confined, hence precluding Z2 topological order. This strongly indicates the gapless (or critical) nature of the QSL phase, which is most naturally consistent with a U(1) Dirac spin liquid scenario.
The U(1) Dirac spin liquid might realize an exotic phase of matter whose low-energy properties are described by quantum electrodynamics in 2+1 dimensions, where gapless modes exist but spinons and gauge fields are strongly coupled. Its existence was proposed in frustrated Heisenberg models in the presence of frustrating superexchange interactions by the (Abrikosov) fermionic representation of the spin operators [Wen, Phys. Rev. B 65, 165113 (2002)], supplemented by the Gutzwiller projection. Here, we construct charge-Q monopole excitations in the Heisenberg model on the triangular lattice with nearest-neighbor (J1) and next-nearest-neighbor (J2) couplings. In the highly frustrated regime, singlet and triplet monopoles with Q=1 become gapless in the thermodynamic limit; in addition, the energies for generic Q agree with field-theoretical predictions, obtained for a large number of gapless fermion modes. Finally, we consider localized gauge excitations, in which magnetic $\ensuremath{\pi}$ fluxes are concentrated in the triangular plaquettes (in analogy with Z2 visons), showing that these kinds of states do not play a relevant role at low energies. All our findings lend support to a stable U(1) Dirac spin liquid in the ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model on the triangular lattice.
Quantum spin liquids (QSLs) host exotic fractionalized magnetic and gauge-field excitations whose microscopic origins and experimental verification remain frustratingly elusive. In the absence of static magnetic order, the spin excitation spectrum constitutes the crucial probe of QSL behavior, but its theoretical computation is a serious challenge. Here we employ state-of-the-art tensor-network methods to obtain the full dynamical spectral function of the J1-J2 kagome Heisenberg model and benchmark our results by tracking their evolution across the magnetically ordered and QSL phases. Reducing ∣J2∣/J1 causes increasingly strong spin-wave renormalization, flattening these modes then merging them into a continuum characteristic of deconfined spinons at all finite energies in the QSL. The low-energy continuum and the occurrence of gap closure at multiple high-symmetry points identify this gapless QSL as the U(1) Dirac spin liquid. These results establish a unified understanding of spin excitations in highly frustrated quantum magnets and provide clear spectral fingerprints for experimental detection in candidate kagome QSL materials.
Understanding the nature of the excitation spectrum in quantum spin liquids is of fundamental importance, in particular for the experimental detection of candidate materials.However, current theoretical and numerical techniques have limited capabilities, especially in obtaining the dynamical structure factor, which gives a crucial characterization of the ultimate nature of the quantum state and may be directly assessed by inelastic neutron scattering.In this work, we investigate the low-energy properties of the S ¼ 1=2 Heisenberg model on the triangular lattice, including both nearest-neighbor J 1 and next-nearestneighbor J 2 superexchanges, by a dynamical variational Monte Carlo approach that allows accurate results on spin models.For J 2 ¼ 0, our calculations are compatible with the existence of a well-defined magnon in the whole Brillouin zone, with gapless excitations at K points (i.e., at the corners of the Brillouin zone).The strong renormalization of the magnon branch (also including rotonlike minima around the M points, i.e., midpoints of the border zone) is described by our Gutzwiller-projected state, where Abrikosov fermions are subject to a nontrivial magnetic π flux threading half of the triangular plaquettes.When increasing the frustrating ratio J 2 =J 1 , we detect a progressive softening of the magnon branch at M, which eventually becomes gapless within the spin-liquid phase.This feature is captured by the band structure of the unprojected wave function (with two Dirac points for each spin component).In addition, we observe an intense signal at low energies around the K points, which cannot be understood within the unprojected picture and emerges only when the Gutzwiller projection is considered, suggesting the relevance of gauge fields for the low-energy physics of spin liquids.
Recent work has shown that the triangular lattice spin-$\frac{1}{2} {J}_{1}\ensuremath{-}{J}_{2}$ Heisenberg and XXZ antiferromagnets may exhibit coplanar or supersolid orders proximate to a gapless Dirac spin-liquid phase. We explore a distinct $\text{SU}(2N)\ifmmode\times\else\texttimes\fi{}\text{SU}(M)$ fermionic parton approach, complemented by variational Monte Carlo calculations for the spin-21 model, to study the phase diagram of these models. We also calculate their dynamical spin response including parton interactions within a random-phase approximation, and discuss implications for neutron scattering on triangular lattice cobaltates Ba3CoSb2O9,Na2BaCo(PO4)2,K2Co(SeO3)2,Rb2Co(SeO3)2, and Yb-based magnet KYbSe2.