No Arabic abstract
The properties of ground state of spin-$frac{1}{2}$ kagome antiferromagnetic Heisenberg (KAFH) model have attracted considerable interest in the past few decades, and recent numerical simulations reported a spin liquid phase. The nature of the spin liquid phase remains unclear. For instance, the interplay between symmetries and $Z_2$ topological order leads to different types of $Z_2$ spin liquid phases. In this paper, we develop a numerical simulation method based on symmetric projected entangled-pair states (PEPS), which is generally applicable to strongly correlated model systems in two spatial dimensions. We then apply this method to study the nature of the ground state of the KAFH model. Our results are consistent with that the ground state is a $U(1)$ Dirac spin liquid rather than a $Z_2$ spin liquid.
We study the spin-$1/2$ Heisenberg model on the triangular lattice with the antiferromagnetic first ($J_1$) and second ($J_2$) nearest-neighbor interactions using density matrix renormalization group. By studying the spin correlation function, we find a $120^{circ}$ magnetic order phase for $J_2 lesssim 0.07 J_1$ and a stripe antiferromagnetic phase for $J_2 gtrsim 0.15 J_1$. Between these two phases, we identify a spin liquid region characterized by the exponential decaying spin and dimer correlations, as well as the large spin singlet and triplet excitation gaps on finite-size systems. We find two near degenerating ground states with distinct properties in two sectors, which indicates more than one spin liquid candidates in this region. While the sector with spinon is found to respect the time reversal symmetry, the even sector without a spinon breaks such a symmetry for finite-size systems. Furthermore, we detect the signature of the fractionalization by following the evolution of different ground states with inserting spin flux into the cylinder system. Moreover, by tuning the anisotropic bond coupling, we explore the nature of the spin liquid phase and find the optimal parameter region for the gapped $Z_2$ spin liquid.
We study the zero-temperature phase diagram of the spin-$frac{1}{2}$ Heisenberg model with breathing anisotropy (i.e., with different coupling strength on the upward and downward triangles) on the kagome lattice. Our study relies on large scale tensor network simulations based on infinite projected entangled-pair state and infinite projected entangled-simplex state methods adapted to the kagome lattice. Our energy analysis suggests that the U(1) algebraic quantum spin-liquid (QSL) ground-state of the isotropic Heisenberg model is stable up to very large breathing anisotropy until it breaks down to a critical lattice-nematic phase that breaks rotational symmetry in real space through a first-order quantum phase transition. Our results also provide further insight into the recent experiment on vanadium oxyfluoride compounds which has been shown to be relevant platforms for realizing QSL in the presence of breathing anisotropy.
We report $^{35}$Cl NMR, ESR, $mu$SR and specific heat measurements on the $S=1/2$ frustrated kagome magnet kapellasite, $alpha-$Cu$_3$Zn(OH)$_6$Cl$_2$, where a gapless spin liquid phase is stabilized by a set of competing exchange interactions. Our measurements confirm the ferromagnetic character of the nearest-neighbour exchange interaction $J_1$ and give an energy scale for the competing interactions $|J| sim 10$ K. The study of the temperature-dependent ESR lineshift reveals a moderate symmetric exchange anisotropy term $D$, with $|D/J|sim 3$%. These findings validate a posteriori the use of the $J_1 - J_2 - J_d$ Heisenberg model to describe the magnetic properties of kapellasite [Bernu et al., Phys. Rev. B 87, 155107 (2013)]. We further confirm that the main deviation from this model is the severe random depletion of the magnetic kagome lattice by 27%, due to Cu/Zn site mixing, and specifically address the effect of this disorder by $^{35}$Cl NMR, performed on an oriented polycrystalline sample. Surprisingly, while being very sensitive to local structural deformations, our NMR measurements demonstrate that the system remains homogeneous with a unique spin susceptibility at high temperature, despite a variety of magnetic environments. Unconventional spin dynamics is further revealed by NMR and $mu$SR in the low-$T$, correlated, spin liquid regime, where a broad distribution of spin-lattice relaxation times is observed. We ascribe this to the presence of local low-energy modes.
We study the quantum phase diagram of the spin-$1/2$ Heisenberg model on the kagome lattice with first-, second-, and third-neighbor interactions $J_1$, $J_2$, and $J_3$ by means of density matrix renormalization group. For small $J_2$ and $J_3$, this model sustains a time-reversal invariant quantum spin liquid phase. With increasing $J_2$ and $J_3$, we find in addition a $q=(0,0)$ N{e}el phase, a chiral spin liquid phase, a valence-bond crystal phase, and a complex non-coplanar magnetically ordered state with spins forming the vertices of a cuboctahedron known as a cuboc1 phase. Both the chiral spin liquid and cuboc1 phase break time reversal symmetry in the sense of spontaneous scalar spin chirality. We show that the chiralities in the chiral spin liquid and cuboc1 are distinct, and that these two states are separated by a strong first order phase transition. The transitions from the chiral spin liquid to both the $q=(0,0)$ phase and to time-reversal symmetric spin liquid, however, are consistent with continuous quantum phase transitions.
We present a multiloop pseudofermion functional renormalization group (pffRG) approach to quantum spin systems. As a test case, we study the spin-$tfrac{1}{2}$ 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 $ell geq 5$. The multiloop flow remains stable as the infrared cutoff $Lambda$ is reduced relative to the microscopic exchange interaction $J$, allowing us to reach values of $Lambda/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 $Lambda/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.