No Arabic abstract
Temperature-dependent dynamical spin correlations, which can be readily accessed via a variety of experimental techniques, hold the potential of offering a unique fingerprint of quantum spin liquids and other intriguing dynamical states. In this work we present an in-depth study of the temperature-dependent dynamical spin structure factor $S({bf q}, omega)$ of the antiferromagnetic (AFM) Heisenberg spin-1/2 model on the kagome lattice with additional Dzyaloshinskii--Moriya (DM) interactions. Using the finite-temperature Lanczos method on lattices with up to $N = 30$ sites we find that even without DM interactions, chiral low-energy spin fluctuations of the $120^circ$ AFM order parameter dominate the dynamical response. This leads to a nontrivial frequency dependence of $S({bf q}, omega)$ and the appearance of a pronounced low-frequency mode at the M point of the extended Brillouin zone. Adding an out-of-plane DM interactions $D^z$ gives rise to an anisotropic dynamical response, a softening of in-plane spin fluctuations, and, ultimately, the onset of a coplanar AFM ground-state order at $D^z > 0.1 J$. Our results are in very good agreement with existing inelastic neutron scattering and temperature-dependent NMR spin-lattice relaxation rate ($1/T_1$) data on the paradigmatic kagome AFM herbertsmithite, where the effect of its small $D^z$ on the dynamical spin correlations is shown to be rather small, as well as with $1/T_1$ data on the novel kagome AFM YCu$_3$(OH)$_6$Cl$_3$, where its substantial $D^z approx 0.25 J$ interaction is found to strongly affect the spin dynamics.
We determine dynamical response functions of the S=1/2 Heisenberg quantum antiferromagnet on the kagome lattice based on large-scale exact diagonalizations combined with a continued fraction technique. The dynamical spin structure factor has important spectral weight predominantly along the boundary of the extended Brillouin zone and energy scans reveal broad response extending over a range of 2 sim 3J concomitant with pronounced intensity at lowest available energies. Dispersive features are largely absent. Dynamical singlet correlations -- which are relevant for inelastic light probes -- reveal a similar broad response, with a high intensity at low frequencies omega/J lesssim 0.2J. These low energy singlet excitations do however not seem to favor a specific valence bond crystal, but instead spread over many symmetry allowed eigenstates.
We believe that a necessary first step in understanding the ground state properties of the spin-${scriptstylefrac{1}{2}}$ kagome Heisenberg antiferromagnet is a better understanding of this models very large number of low energy singlet states. A description of the low energy states that is both accurate and amenable for numerical work may ultimately prove to have greater value than knowing only what these properties are, in particular when these turn on the delicate balance of many small energies. We demonstrate how this program would be implemented using the basis of spin-singlet dimerized states, though other bases that have been proposed may serve the same purpose. The quality of a basis is evaluated by its participation in all the low energy singlets, not just the ground state. From an experimental perspective, and again in light of the small energy scales involved, methods that can deliver all the low energy states promise more robust predictions than methods that only refine a fraction of these states.
A clear thermal Hall signal ($kappa_{xy}$) was observed in the spin liquid phase of the $S=1/2$ kagome antiferromagnet Ca kapellasite (CaCu$_3$(OH)$_6$Cl$_2cdot 0.6$H$_2$O). We found that $kappa_{xy}$ is well reproduced, both qualitatively and quantitatively, using the Schwinger-boson mean-field theory with the Dzyaloshinskii--Moriya interaction of $D/J sim 0.1$. In particular, $kappa_{xy}$ values of Ca kapellasite and those of another kagome antiferromagnet, volborthite, converge to one single curve in simulations modeled using Schwinger bosons, indicating a common temperature dependence of $kappa_{xy}$ for the spins of a kagome antiferromagnet.
We investigate the spin-1/2 Heisenberg antiferromagnet on the kagome lattice with breathing anisotropy (i.e. with weak and strong triangular units), constructing an improved simplex Resonating Valence Bond (RVB) ansatz by successive applications (up to three times) of local quantum gates which implement a filtering operation on the bare nearest-neighbor RVB state. The resulting Projected Entangled Pair State involves a small number of variational parameters (only one at each level of application) and preserves full lattice and spin-rotation symmetries. Despite its simple analytic form, the simplex RVB provides very good variational energies at strong and even intermediate breathing anisotropy. We show that it carries $Z_2$ topological order which does not fade away under the first few applications of the quantum gates, suggesting that the RVB topological spin liquid becomes a competing ground state candidate for the kagome antiferromagnet at large breathing anisotropy.
The nature of the ground state of the spin $S=1/2$ Heisenberg antiferromagnet on the kagome lattice with breathing anisotropy (i.e., with different superexchange couplings $J_{vartriangle}$ and $J_{triangledown}$ within elementary up- and down-pointing triangles) is investigated within the framework of Gutzwiller projected fermionic wave functions and Monte Carlo methods. We analyze the stability of the U(1) Dirac spin liquid with respect to the presence of fermionic pairing that leads to a gapped $mathbb{Z}_{2}$ spin liquid. For several values of the ratio $J_{triangledown}/J_{vartriangle}$, the size scaling of the energy gain due to the pairing fields and the variational parameters are reported. Our results show that the energy gain of the gapped spin liquid with respect to the gapless state either vanishes for large enough system size or scales to zero in the thermodynamic limit. Similarly, the optimized pairing amplitudes (responsible for opening the spin gap) are shown to vanish in the thermodynamic limit. Our outcome is corroborated by the application of one and two Lanczos steps to the gapless and gapped wave functions, for which no energy gain of the gapped state is detected when improving the quality of the variational states. Finally, we discuss the competition with the simplex $mathbb{Z}_{2}$ resonating-valence-bond spin liquid, valence-bond crystal, and nematic states in the strongly anisotropic regime, i.e., $J_{triangledown} ll J_{vartriangle}$.