Ever since the introduction of resonating valence bond states by P.\\ W.\\ Anderson in 1970s, the search for exotic many--body quantum behavior in peculiar states of matter known as quantum spin liquids has been a central challenge in the field of condensed matter physics. An important searching ground for such \\emph{quantum--ness} is the low energy physics of frustrated magnets. In a quantum system, frustration---the inability of a physical system to achieve a global state that consistently minimizes its energy locally---enables the unavoidable quantum fluctuations to hybridize an extensive number of classical states into highly entangled quantum states at low temperatures. It can be argued that this is the most practical approach to achieve genuine many--body quantum behavior at large scales. The fundamental source of frustration in many physical systems, is the existence of local non--commuting terms. In the case of antiferromagnet spin theories the simplest examples are triangles of nearest--neighboring sites. Consequently, lattices with triangles as building blocks have long been studied in hopes of discovering spin liquids physics. In particular the kagome lattice, which is composed of corner--sharing triangles, has been the center of at least four decades of spin liquid research. This is due to the underconstrained nature of system made out of corner--sharing triangles. After the first three chapters, introduction in chapter~\\ref{chap:intro}, preliminary physics in chapter~\\ref{chap:prelim}, and a brief introduction to spin--liquids in chapter~\\ref{chap:QSLs}, in chapter~\\ref{chap:kagome} I tackle a simpler version of the long-standing problem of spin--$1/2$ Heisenberg antiferromagnets by considering a quasi--1D lattice consisting entirely of corner--sharing triangles, \\emph{kagome strip}. I illustrate that the standard Heisenberg antiferromagnet Hamiltonian over kagome strip is an extended gapless quantum phase, that is well characterized by two fermionic/bosonic gapless modes and power--law decaying spin and bond--energy correlations. I also demonstrate that the correlation functions oscillate at tunably incommensurate wave vectors. It turns out that this phase can be identified by a particular marginal instability of a two-band spinon Fermi surface coupled to an emergent U(1) gauge field. This interpretation is supported by analytic Abelian bosonization and with extensive numerical large--scale density matrix renormalization group study as well as variational Monte Carlo calculations on Gutzwiller ans\\"atz wave functions. This intriguing result is the first numerical demonstration of emergent fermionic spinons in a simple SU(2) invariant nearest-neighbor Heisenberg model beyond the strictly 1D (Bethe chain) limit. The unexpected success of the fermionic spinons in describing the low energy physics of kagome antiferromagnet, reintroduces the questions about the validity of spinon physics as an effective theory of quantum spin liquids. However, numerical methods are still far behind the theoretical advances. And still even today, an accurate description of the quantum spin liquid states using tensor network methods is notoriously challenging. It is known that for large quasi-1D systems, the density matrix renormalization group and related methods usually require significant computational resources and sometimes fail to converge to a satisfactory state. On the other hand, variational wavefunctions acquired from the Gutzwiller projection of gaussian fermionic theories have long served as both a theoretical starting point for the construction of such spin liquid states and as an inspiration for numerical variational Monte Carlo (VMC) to calculate observables of interest. Noting this observation I examine a different method by exploring the possibility of constructing a matrix product state (MPS) representation for a Gutzwiller--projected state from two given MPS representations of gaussian fermionic theories in the \\ref{chap:gutz} chapter. I investigate the complexity of different approaches to achieve Gutzwiller projection for MPSs and introduce the novel algorithm which we call the Guzwiller zipper method. The performance of the algorithm is tested against two copies of a single half-filled band of spin-1/2 fermionic spinons. In a successful attempt to describe the nature of spin liquid states on quasi-1D strips of triangular and kagome-like lattices, I apply this method to two MPS of multi-band fermionic spinon theories and compare with the complexity of the traditional VMC approach. In particular, we methodically disprove the conjecture of a spinon fermi surface spin liquid for the triangular lattice. Finally, I conclude this thesis by laying out a bird's eye view of the current and future of spin liquid physics in kagome and triangular lattices and pointing out numerous possible applications of the novel Gutzwiller zipper method. In addition, I would also discuss its possible extensions to more complicated tensor networks as well as versions for higher spatial dimensions.
Ever since the introduction of resonating valence bond states by P.\ W.\ Anderson in 1970s, the search for exotic many--body quantum behavior in peculiar states of matter known as quantum spin liquids has been a central challenge in the field of condensed matter physics. An important searching ground for such \emph{quantum--ness} is the low energy physics of frustrated magnets. In a quantum system, frustration---the inability of a physical system to achieve a global state that consistently minimizes its energy locally---enables the unavoidable quantum fluctuations to hybridize an extensive number of classical states into highly entangled quantum states at low temperatures. It can be argued that this is the most practical approach to achieve genuine many--body quantum behavior at large scales. The fundamental source of frustration in many physical systems, is the existence of local non--commuting terms. In the case of antiferromagnet spin theories the simplest examples are triangles of nearest--neighboring sites. Consequently, lattices with triangles as building blocks have long been studied in hopes of discovering spin liquids physics. In particular the kagome lattice, which is composed of corner--sharing triangles, has been the center of at least four decades of spin liquid research. This is due to the underconstrained nature of system made out of corner--sharing triangles. After the first three chapters, introduction in chapter~\ref{chap:intro}, preliminary physics in chapter~\ref{chap:prelim}, and a brief introduction to spin--liquids in chapter~\ref{chap:QSLs}, in chapter~\ref{chap:kagome} I tackle a simpler version of the long-standing problem of spin--1/2 Heisenberg antiferromagnets by considering a quasi--1D lattice consisting entirely of corner--sharing triangles, \emph{kagome strip}. I illustrate that the standard Heisenberg antiferromagnet Hamiltonian over kagome strip is an extended gapless quantum phase, that is well characterized by two fermionic/bosonic gapless modes and power--law decaying spin and bond--energy correlations. I also demonstrate that the correlation functions oscillate at tunably incommensurate wave vectors. It turns out that this phase can be identified by a particular marginal instability of a two-band spinon Fermi surface coupled to an emergent U(1) gauge field. This interpretation is supported by analytic Abelian bosonization and with extensive numerical large--scale density matrix renormalization group study as well as variational Monte Carlo calculations on Gutzwiller ans\"atz wave functions. This intriguing result is the first numerical demonstration of emergent fermionic spinons in a simple SU(2) invariant nearest-neighbor Heisenberg model beyond the strictly 1D (Bethe chain) limit. The unexpected success of the fermionic spinons in describing the low energy physics of kagome antiferromagnet, reintroduces the questions about the validity of spinon physics as an effective theory of quantum spin liquids. However, numerical methods are still far behind the theoretical advances. And still even today, an accurate description of the quantum spin liquid states using tensor network methods is notoriously challenging. It is known that for large quasi-1D systems, the density matrix renormalization group and related methods usually require significant computational resources and sometimes fail to converge to a satisfactory state. On the other hand, variational wavefunctions acquired from the Gutzwiller projection of gaussian fermionic theories have long served as both a theoretical starting point for the construction of such spin liquid states and as an inspiration for numerical variational Monte Carlo (VMC) to calculate observables of interest. Noting this observation I examine a different method by exploring the possibility of constructing a matrix product state (MPS) representation for a Gutzwiller--projected state from two given MPS representations of gaussian fermionic theories in the \ref{chap:gutz} chapter. I investigate the complexity of different approaches to achieve Gutzwiller projection for MPSs and introduce the novel algorithm which we call the Guzwiller zipper method. The performance of the algorithm is tested against two copies of a single half-filled band of spin-1/2 fermionic spinons. In a successful attempt to describe the nature of spin liquid states on quasi-1D strips of triangular and kagome-like lattices, I apply this method to two MPS of multi-band fermionic spinon theories and compare with the complexity of the traditional VMC approach. In particular, we methodically disprove the conjecture of a spinon fermi surface spin liquid for the triangular lattice. Finally, I conclude this thesis by laying out a bird's eye view of the current and future of spin liquid physics in kagome and triangular lattices and pointing out numerous possible applications of the novel Gutzwiller zipper method. In addition, I would also discuss its possible extensions to more complicated tensor networks as well as versions for higher spatial dimensions.
The search for exotic quantum spin liquid states in simple yet realistic spin models remains a central challenge in the field of frustrated quantum magnetism. Here we consider the canonical nearest-neighbor kagome Heisenberg antiferromagnet restricted to a quasi-one-dimensional strip consisting entirely of corner-sharing triangles. Using large-scale density matrix renormalization group calculations, we identify in this model an extended gapless quantum phase characterized by central charge c=2 and power-law decaying spin and bond-energy correlations which oscillate at tunably incommensurate wave vectors. We argue that this intriguing spin liquid phase can be understood as a marginal instability of a two-band spinon Fermi surface coupled to an emergent U(1) gauge field, an interpretation which we substantiate via bosonization analysis and Monte Carlo calculations on model Gutzwiller variational wave functions. Our results represent one of the first numerical demonstrations of emergent fermionic spinons in a simple SU(2) invariant nearest-neighbor Heisenberg model beyond the strictly one-dimensional (Bethe chain) limit.
The search for exotic quantum spin liquid states in simple yet realistic spin models remains a central challenge in the field of frustrated quantum magnetism. Here we consider the canonical nearest-neighbor kagome Heisenberg antiferromagnet restricted to a quasi-one-dimensional strip consisting entirely of corner-sharing triangles. Using large-scale density matrix renormalization group calculations, we identify in this model an extended gapless quantum phase characterized by central charge c = 2 and power-law decaying spin and bond-energy correlations which oscillate at tunably incommensurate wave vectors. We argue that this intriguing spin liquid phase can be understood as a marginal instability of a two-band spinon Fermi surface coupled to an emergent U(1) gauge field, an interpretation which we substantiate via bosonization analysis and Monte Carlo calculations on model Gutzwiller variational wave functions. Our results represent one of the first numerical demonstrations of emergent fermionic spinons in a simple SU(2) invariant nearest-neighbor Heisenberg model beyond the strictly one-dimensional (Bethe chain) limit.
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.
We study the ground-state properties of the quantum spin liquid (QSL) phases of the spin-1/2 antiferromagnetic Heisenberg model on the triangular lattice with nearest- (J1), next-nearest- (J2), and third-neighbor (J3) interactions by using the density-matrix renormalization group (DMRG) method. By combining parallel DMRG with SU(2) spin rotational symmetry, we are able to obtain accurate results on large cylinders with length up to Lx=48 and circumference Ly=6−−12. Our results suggest that the QSL phase of the ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model is gapped, which is characterized by the absence of a gapless mode and by short-range spin-spin and dimer-dimer correlations. In the presence of J3 interaction, we find that a critical QSL with a single gapless mode emerges. While both spin-spin and scalar chiral-chiral correlations are short ranged, dimer-dimer correlations are quasi long ranged and decay as a power law at long distances.
Motivated by recent density-matrix renormalization group (DMRG) calculations [Yan, Huse, and White, Science 332, 1173 (2011)], which claimed that the ground state of the nearest-neighbor spin-1/2 Heisenberg antiferromagnet on the kagome lattice geometry is a fully gapped spin liquid with numerical signatures of Z2 gauge structure, and a further theoretical work [Lu, Ran, and Lee, Phys. Rev. B 83, 224413 (2011)], which gave a classification of all Schwinger-fermion mean-field fully symmetric Z2 spin liquids on the kagome lattice, we have thoroughly studied Gutzwiller-projected fermionic wave functions by using quantum variational Monte Carlo techniques, hence implementing exactly the constraint of one fermion per site. In particular, we investigated the energetics of all Z2 candidates (gapped and gapless) that lie in the neighborhood of the energetically competitive U(1) gapless spin liquids. By using a state-of-the-art optimization method, we were able to conclusively show that the U(1) Dirac state is remarkably stable with respect to all Z2 spin liquids in its neighborhood, and in particular for opening a gap toward the so-called ${\mathbb{Z}}_{2}[0,\ensuremath{\pi}]\ensuremath{\beta}$ state, which was conjectured to describe the ground state obtained by the DMRG method. Finally, we also considered the addition of a small second nearest-neighbor exchange coupling of both antiferromagnetic and ferromagnetic type, and obtained similar results, namely, a U(1) Dirac spin-liquid ground state.
In this dissertation we present the results of a theoretical investigation of spin models on two-dimensional and quasi one-dimensional lattices, all unified under the concept of quantum frustrated antiferromagnetism, and all discussing various aspects of the antiferromagnetic Heisenberg model on the kagomé lattice. In the Introduction (Chapter 1), we discuss at some length such concepts as frustration and superexchange, among others, which are of common relevance in the rest of the chapters. In Chapter 2, we study the effect of Dzyaloshinskii--Moriya (DM) interactions on the zero-temperature magnetic susceptibility of systems whose low energy can be described by short-range valence bond states. Our work shows that this treatment is consistent with the experimentally observed non-vanishing susceptibility -- in the specified temperature limit -- of the spin-1/2 kagomé antiferromagnetic compound ZnCu3(OH)6Cl2, also known as herbertsmithite. Although the objective of this work is explaining the aforementioned characteristic of the experimental system, our methods are more general and we apply them to the checkerboard and Shastry-Sutherland lattices as well. In Chapter 3, we discuss our findings in the study of ghost-mediated domain wall interactions in the diamondback ladder. These domain walls are the the spin excitations -- the kinks and the antikinks -- separating the ground states along one chain of the ladder. While as individual entities an antikink is energy costly and a kink energy free, our study finds that both interact via the ghosts that they produce in the opposite side of the ladder from where they are located. Through the study of these ghosts, we find that domain walls proliferate in the system above a critical value of the system's coupling constants. It is this proliferation that makes their treatment as free, non-interacting particles impossible, so we study here their interactions both quantitatively and qualitatively, in a region where the latter are yet not very strong, namely below the critical point. Based on the calculated two-body interaction potential, domain walls interact attractively (repulsively) when separated at even (odd) distances, with a strength that decays as 1/sp, where s is their separation and p<1. We also consider higher-order interactions. In the last chapter, Chapter 4, we present our study of the spin-1 kagomé Heisenberg antiferromagnet. Our approach is to first consider an SU(2)-symmetric parent Hamiltonian with known ground states on the S=1 kagomé lattice, in which nearest-neighbor Heisenberg interactions are already present. We then enhance these interactions by an additional Heisenberg term added perturbatively in order to move the system closer to a pure Heisenberg antiferromagnet. The results of this enhancement is obtaining a description of the system in terms of an effective Hamiltonian, namely a transverse field Ising AF on the triangular lattice. Based on the particular values of this effective Hamiltonian, our system is found to be in the order-by-disorder phase.