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.
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.
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.
We study the ground-state phase diagram of the quantum
spin-1/2 Heisenberg model on the kagome lattice with first-
(J1<0), second-
(J2<0), and third-neighbor interactions
(Jd>0) by means of analytical low-energy field theory and numerical density-matrix renormalization group (DMRG) studies. The results offer a consistent picture of the
Jd-dominant regime in terms of three sets of spin chains weakly coupled by the ferromagnetic interchain interactions
J1,2. When either
J1 or
J2 is much stronger than the other one, the model is found to support one of two cuboctohedral phases, cuboc1, and cuboc2. These cuboc states host noncoplanar long-ranged magnetic order and possess finite scalar spin chirality. However, in the compensated regime ${J}_{1}\ensuremath{\simeq}{J}_{2}$, a valence bond crystal phase emerges between the two cuboc phases. We find excellent agreement between an analytical theory based on coupled spin chains and unbiased DMRG calculations, including at a very detailed level of comparison of the structure of the valence bond crystal state. To our knowledge, this is the first such comprehensive understanding of a highly frustrated two-dimensional quantum antiferromagnet. We find no evidence of either the one-dimensional gapless spin liquid or the chiral spin liquids, which were previously suggested by parton mean-field theories.
Motivated by recent experiments on the Heisenberg
S=1/2 quantum spin liquid candidate material kapellasite, we classify all possible chiral (time-reversal symmetry breaking) spin liquids with fermionic spinons on the kagome lattice. We obtain the phase diagram for the physically relevant extended Heisenberg model, comparing the energies of a wide range of microscopic variational wave functions. We propose that, at low temperature, kapellasite exhibits a gapless chiral spin liquid phase with spinon Fermi surfaces. This two-dimensional state inherits many properties of the nearby one-dimensional phase of decoupled antiferromagnetic spin chains, but also shows some remarkable differences. We discuss the spin structure factors and other physical properties.
We analyze the antiferromagnetic spin-1/2 XXZ model on the kagome lattice at finite external magnetic field with the help of a non-perturbative zero-temperature renormalization group (RG) technique. The exact nature of the ground and excited state properties (e.g. gapped or gapless spectrum etc) of this system are still debated. Approximate methods have typically been adopted towards understanding the low-energy spectrum. Following the work of Kumar et al (2014 Phys. Rev. B 90 174409), we use a Jordan–Wigner transformation to map the spin problem into one of spinless fermions (spinons) in the presence of a statistical gauge field, and with nearest-neighbor interactions. While the work of Kumar et al was confined mostly to the plateau at 1/3-filling (magnetization per site) in the XY regime, we analyze the role of inter-spinon interactions in shaping the phases around this plateau in the entire XXZ model. The RG phase diagram obtained contains three spin liquid phases whose position is determined as a function of the exchange anisotropy and the energy scale for fluctuations arising from spinon scattering. Two of these spins liquids are topologically ordered states of matter with gapped, degenerate states on the torus. The gap for one of these phases corresponds to the one-spinon band gap of the Azbel–Hofstadter spectrum for the XY part of the Hamiltonian, while the other arises from two-spinon interactions. The Heisenberg point of this problem is found to lie within the interaction gapped spin liquid phase, in broad agreement with a recent experimental finding. The third phase is an algebraic spin liquid with a gapless Dirac spectrum for spinon excitations, and possess properties that show departures from the Fermi liquid paradigm. The three phase boundaries correspond to critical theories, and meet at a SU (2)-symmetric multicritical point. This special critical point agrees well with the gap-closing transition point predicted by Kumar et al . We discuss the relevance of our findings to various recent experiments, as well as results obtained from other theoretical analyses.