We combine the pseudofermion functional renormalization group (PFFRG) method with a self-consistent Fock-like mean-field scheme to calculate low-energy effective theories for emergent spinon excitations in spin-1/2 quantum spin liquids. Using effective spin interactions from PFFRG as an input for the Fock equation and allowing for the most general types of free spinon Ans\"atze as classified by the projective symmetry group (PSG) method, we are able to systematically determine spinon band structures for spin-liquid candidate systems beyond mean-field theory. We apply this approach to the antiferromagnetic ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model on the square lattice and to the antiferromagnetic nearest-neighbor Heisenberg model on the kagome lattice. For the ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ model, we find that in the regime of maximal frustration a SU(2) $\ensuremath{\pi}$-flux state with Dirac spinons yields the largest mean-field amplitudes. For the kagome model, we identify a gapless ${\mathbb{Z}}_{2}$ spin liquid with a small circular spinon Fermi surface and approximate Dirac cones at low but finite energies.
We combine the pseudofermion functional renormalization group (PFFRG) method with a self-consistent Fock-like mean-field scheme to calculate low-energy effective theories for emergent spinon excitations in spin-1/2 quantum spin liquids. Using effective spin interactions from PFFRG as an input for the Fock equation and allowing for the most general types of free spinon Ans"atze as classified by the projective symmetry group (PSG) method, we are able to systematically determine spinon band structures for spin-liquid candidate systems beyond mean-field theory. We apply this approach to the antiferromagnetic J1\ensuremath-J2 Heisenberg model on the square lattice and to the antiferromagnetic nearest-neighbor Heisenberg model on the kagome lattice. For the J1\ensuremath-J2 model, we find that in the regime of maximal frustration a SU(2) \ensuremathπ-flux state with Dirac spinons yields the largest mean-field amplitudes. For the kagome model, we identify a gapless Z2 spin liquid with a small circular spinon Fermi surface and approximate Dirac cones at low but finite energies.
The ruby lattice has been the subject of much interest recently due its realization in Rydberg atom arrays, where its rich variety of frustrated interactions gives rise to topologically ordered quantum spin liquids. Similarly, numerical studies of ruby-lattice spin models, with both isotropic and anisotropic interactions, have provided evidence of gapped and gapless spin-liquid ground states with different low-energy gauge structures. Motivated by these findings, we perform a projective symmetry group (PSG) classification of U(1) and Z2 fermionic spinon mean-field theories$\unicode{x2014}$respecting space-group and time-reversal symmetries$\unicode{x2014}$for S=1/2 spins. We obtain a total of 50 U(1) and 64 Z2 PSGs, and upon restricting their realization via mean-field Ansa¨tze with up to second-nearest-neighbor singlet amplitudes (relevant to the models studied here), only 8 U(1) and 18 Z2 distinct states are obtained. We present the singlet fields for all Ansa¨tze up to third-nearest-neighboring bonds and discuss their spinon dispersions as well as their dynamical spin structure factors. Building on this information, we also obtain the phase diagram of the Heisenberg model in the presence of first (J1), second (J1′), and third (J2) neighbor antiferromagnetic couplings within a self-consistent mean-field approximation.
Multiloop pseudofermion functional renormalization for quantum spin systems: Application to the spin-21 kagome Heisenberg model
We present a multiloop pseudofermion functional renormalization group (pffRG) approach to quantum spin systems. As a test case, we study the spin-21 Heisenberg model on the kagome lattice, a prime example of a geometrically frustrated magnet believed to host a quantum spin liquid. Our main physical result is that, at pure nearest-neighbor coupling, the system shows indications for an algebraic spin liquid through slower-than-exponential decay with distance for the static spin susceptibility, while the pseudofermion self-energy develops intriguing low-energy features. Methodologically, the pseudofermion representation of spin models inherently yields a strongly interacting system, and the quantitative reliability of a truncated fRG flow is \textit{a priori} unclear. Our main technical result is the demonstration of convergence in loop number within multiloop pffRG. Through correspondence with the self-consistent parquet equations, this provides further evidence for the internal consistency of the approach. The loop order required for convergence of the pseudofermion vertices is rather large, but the spin susceptibility is more benign, appearing almost fully converged for loop orders ℓ≥5. The multiloop flow remains stable as the infrared cutoff Λ is reduced relative to the microscopic exchange interaction J, allowing us to reach values of Λ/J on the subpercent level in the spin-liquid phase. By contrast, solving the parquet equations via direct fixed-point iteration becomes increasingly difficult for low Λ/J. We also scrutinize the pseudofermion constraint of single occupation per site, which is only fulfilled on average in pffRG, by explicitly computing fermion-number fluctuations. Although the latter are not entirely suppressed, we find that they do not affect the qualitative conclusions drawn from the spin susceptibility.
We enhance the pseudo-fermion functional-renormalisation-group (PFFRG) method for the investigation of highly frustrated quantum spin systems in two significant ways. On one hand, we generalise the method towards off-diagonal Dzyaloshinsky-Moriya (DM) interactions and single-ion anisotropies. On the other hand, we develop a self-consistent Fock-like mean-field technique which is based on vertex functions from a PFFRG analysis and calculates effective free spinon models for the considered systems at low energies. In Chap. 1, we introduce the fundamental mathematical concepts on which our studies rely. The so-called flow equations for the PFFRG are then derived in Chap. 3. The first part of our original work is presented in Chap. 4 where we, for the first time, extend the PFFRG scheme towards off-diagonal spin interactions and therefore confirm the quite general applicability of this method. There are two important results from our investigations of the spin-1/2 Heisenberg-DM model on the kagome lattice in this chapter. One is that a finite DM interaction enlarges magnetically ordered phases and thereby shrinks the paramagnetic phases which possibly house spin liquid states. The other important result is that the mineral herbertsmithite which inspired this project lies in close proximity to a quantum critical point. We obtain a phase diagram for the considered Hamiltonian and qualitatively good results in comparison to experimental data for herbertsmithite [68,69]. Our collaboration with a research group from Cologne is summarised in Chap. 5. Here, we investigate two different XXZ models with tetragonally split-up second-neighbour Heisenberg interactions and single-ion anisotropies on the diamond lattice in order to study the spin-1 A-site spinel NiRh2O4. We are able to confirm the existence of a paramagnetic phase for one of the considered models. However, such a paramagnetic behaviour cannot be verified for the Hamiltonians suggested in Refs. [29] and [30]. The second main part of this thesis is contained in Chap. 6. The developed PFFRG plus mean-field method enables us to characterise the effective low-energy spinon band structure of the spin-1/2 J1-J2 (nearest-neighbour) Heisenberg model on the square (kagome) lattice. This is achieved by comparing self-consistently determined overall amplitudes for different mean-field ansätze which are taken from corresponding PSG classifications in Refs. [182] and [111]. Our studies imply that the square lattice model realises a SU(2) spin liquid with Pi-fluxes through elementary plaquettes and that the spin-1/2 kagome antiferromagnet is most likely described by a Z2 spin liquid called Z2[0,Pi]Alpha. We argue that our analysed models have current experimental relevance and that the developed techniques show plenty of possibilities for future applications. A detailed summary and discussion of all methods and results is provided in Chap. 7.
Quantum spin liquids are novel phases of matter whose physical properties cannot be captured within the conventional framework of phase transitions. These quantum states are characterized by excitations, called spinons, carrying a fractional value of the integer spin. A description based on a mean-field approximation leads to effective gauge theory. Due to this gauge freedom symmetry group operations act via projective symmetry group (PSG) representations, which can be used to classify quantum states beyond the conventional method. This introducing discussion is found in Chapters 2-6. A preceding PSG classification on the square lattice, which took spin rotational symmetry breaking terms into account, showed that the resulting mean-field states could exhibit topologically non-trivial spinon bands. In Chapter 7, only those el- ements of this Classification are studied, which lead to short-ranged mean-field models. The mutual effect of spinons and visons is scrutinized by adding static configurations of a background gauge field. Spinons and visons form composite ob- jects, which in the case of topological spinon bands, are described by Majorana zero modes. An exact mapping corroborates the relation to topological supercon- ductors and vortices. The effect of multiple visons is investigated, and a possible experimental signature is exposed. Chapter 8 explains how spinon mean-field theory can be used to derive a phe- nomenological effective model for the compound Ca 10 Cr 7 O 28 . A material that was recently proposed as a quantum spin liquid. This low-energy theory consists of spinon hopping processes. The underlying bilayer kagome structure can be simpli- fied to an effective honeycomb lattice. The dynamical spin structure factor is calcu- lated and compared to neutron scattering data. It is discussed how f -wave pairing of spinons can explain small deviations from a linear temperature dependency observed in heat capacity measurements at the low temperatures. This pairing mechanism breaks the U (1) invariance group of the pure hopping model and renders a coherent picture of Ca 10 Cr 7 O 28 being a Z 2 quantum spin liquid. Chapter 9 carries out a PSG classification for the simple-, body centered-, and face centered cubic lattices. Despite a large number of possible PSG representations, only very few describe short-ranged mean-field states. These mean-field models are constructed up to third neighbor couplings and calculated in a self-consistent manner. It is found that the energy is lower for non-trivial sign patterns induced by non-trivial representations of the symmetry generators. Higher neighbor extensions, including pairing terms, yield lower mean-field ground state energies compared to extensions with only hopping terms. For the fcc lattice, a first neighbor state could be identified, which exhibits symmetry protected zero-energy modes.
We develop a generalized pseudofermion functional renormalization group (PFFRG) approach that can be applied to arbitrary Heisenberg models with spins ranging from the quantum case S=1/2 to the classical limit $S\ensuremath{\rightarrow}\ensuremath{\infty}$. Within this framework, spins of magnitude S are realized by implementing M=2S copies of spin-1/2 degrees of freedom on each lattice site. We confirm that even without explicitly projecting onto the highest spin sector of the Hilbert space, ground states tend to select the largest possible local spin magnitude. This justifies the average treatment of the pseudofermion constraint in previous spin-1/2 PFFRG studies. We apply this method to the antiferromagnetic ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ honeycomb Heisenberg model with nearest-neighbor J1>0 and second-neighbor J2>0 interactions. Mapping out the phase diagram in the ${J}_{2}/{J}_{1}\text{\ensuremath{-}}S$ plane, we find that upon increasing S, quantum fluctuations are rapidly decreasing. In particular, already at S=1 we find no indication for a magnetically disordered phase. In the limit $S\ensuremath{\rightarrow}\ensuremath{\infty}$, the known phase diagram of the classical system is exactly reproduced. More generally, we prove that for $S\ensuremath{\rightarrow}\ensuremath{\infty}$ the PFFRG approach is identical to the Luttinger-Tisza method.
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.