We compute the momentum-resolved dynamical spin structure factor $S(k,\ensuremath{\omega})$ of the SU(4) Heisenberg model on the honeycomb lattice assuming the $\ensuremath{\pi}$-flux Dirac spin liquid ground state by using two methods: (i) variationally using Gutzwiller projected particle-hole excitations of the $\ensuremath{\pi}$-flux Fermi sea, and (ii) in the noninteracting parton mean-field picture. The two approaches produce qualitatively similar results. Based on this analogy, we argue that the energy spectrum of the projected excitations is a gapless continuum of fractional excitations. Quantitatively, the Gutzwiller projection shifts the weight from higher to lower energies, thus emphasizing the lower edge of the continuum. In the mean-field approach, we obtained the
1/distance4 decay of the spin correlation function, and the local correlations show ${S}_{\text{MF}}^{33}(\ensuremath{\omega})\ensuremath{\propto}{\ensuremath{\omega}}^{3}$ behavior.
By combining the density matrix renormalization group (DMRG) method with Gutzwiller projected wave functions, we study the SU(4) symmetric spin-orbital model on the honeycomb lattice. We find that the ground states can be well described by a Gutzwiller projected $\ensuremath{\pi}$-flux state with Dirac-type gapless excitations at one quarter filling. Although these Dirac points are gapped by emergent gauge fluxes on finite cylinders, they govern the critical behavior in the thermodynamic limit. By inserting a $\ensuremath{\theta}=\ensuremath{\pi}$ spin flux to twist the boundary condition, we can shift the gapless sector to the ground state, which provides compelling evidence for the presence of a gapless Dirac spin-orbital liquid.
We perform large-scale density matrix renormalization group simulations of the SU(4) Heisenberg model on the honeycomb lattice to address the long-standing question of its ground state in an unbiased and quantitatively controlled manner. We find reliable numerical evidence that the ground state is a gapped spin-orbital liquid, presumably with a Z_{4} topological order, characterized by a finite topological entanglement entropy close to ln(4), the absence of both SU(4) and lattice symmetry breaking, and a variationally optimized ground-state energy well below the previously proposed π-flux variational state. By exploiting full SU(4) symmetry and keeping up to 12 800 SU(4) multiplets, corresponding to more than 1×10^{6} U(1) states, we achieve unprecedented accuracy for two-dimensional SU(4) quantum magnets. Finite-size scaling of energies and entanglement entropies supports a robust gapped phase in the two-dimensional limit, while a gapless critical state on narrow cylinders is identified as a proximate remnant of a Dirac spin-orbital liquid. Our results find the SU(4) honeycomb Heisenberg model a realization of a gapped topological spin-orbital liquid and provide convincing numerical evidence for topological order in a highly symmetric two-dimensional quantum magnet.
Understanding the nature of the excitation spectrum in quantum spin liquids is of fundamental importance, in particular for the experimental detection of candidate materials.However, current theoretical and numerical techniques have limited capabilities, especially in obtaining the dynamical structure factor, which gives a crucial characterization of the ultimate nature of the quantum state and may be directly assessed by inelastic neutron scattering.In this work, we investigate the low-energy properties of the S ¼ 1=2 Heisenberg model on the triangular lattice, including both nearest-neighbor J 1 and next-nearestneighbor J 2 superexchanges, by a dynamical variational Monte Carlo approach that allows accurate results on spin models.For J 2 ¼ 0, our calculations are compatible with the existence of a well-defined magnon in the whole Brillouin zone, with gapless excitations at K points (i.e., at the corners of the Brillouin zone).The strong renormalization of the magnon branch (also including rotonlike minima around the M points, i.e., midpoints of the border zone) is described by our Gutzwiller-projected state, where Abrikosov fermions are subject to a nontrivial magnetic π flux threading half of the triangular plaquettes.When increasing the frustrating ratio J 2 =J 1 , we detect a progressive softening of the magnon branch at M, which eventually becomes gapless within the spin-liquid phase.This feature is captured by the band structure of the unprojected wave function (with two Dirac points for each spin component).In addition, we observe an intense signal at low energies around the K points, which cannot be understood within the unprojected picture and emerges only when the Gutzwiller projection is considered, suggesting the relevance of gauge fields for the low-energy physics of spin liquids.
We compute the dynamical spin structure factor $S(k,\ensuremath{\omega})$ of the SU(3) Heisenberg chain variationally using a truncated Hilbert space spanned by the Gutzwiller projected particle-hole excitations of the Fermi sea, introduced by Dalla Piazza et al. [Nat. Phys. 11, 62 (2015)], with a modified importance sampling. We check the reliability of the method by comparing the $S(k,\ensuremath{\omega})$ to exact diagonalization results for 18 sites and to the two-soliton continuum of the Bethe ansatz for 72 sites. We get an excellent agreement in both cases. Detailed analysis of the finite-size effects shows that the method captures the critical Wess-Zumino-Witten SU(3)
1 behavior and reproduces the correct exponent, with the exception of the size dependence of the weight of the bottom of the conformal tower. We also calculate the single-mode approximation for the
SU(N) Heisenberg model and determine the velocity of excitations. Finally, we apply the method to the SU(3) Haldane-Shastry model and find that the variational method gives the exact wave function for the lowest excitation at $k=\ifmmode\pm\else\textpm\fi{}2\ensuremath{\pi}/3$.
A prominent feature of quantum spin liquids is fractionalization of the spin degree of freedom. Fractionalized excitations have their own dynamics in different energy scales, and hence, affect finite-temperature (
T) properties in a peculiar manner even in the paramagnetic state harboring the quantum spin liquid state. We here present a comprehensive theoretical study of the spin dynamics in a wide
T range for the Kitaev model on a honeycomb lattice, whose ground state is a quantum spin liquid. In this model, the fractionalization occurs to break up quantum spins into itinerant matter fermions and localized
Z2 fluxes, which results in two crossovers at very different
T scales. Extending the previous study for the isotropic coupling case [J. Yoshitake, J. Nasu, and Y. Motome, Phys. Rev. Lett. 117, 157203 (2016)], we calculate the dynamical spin structure factor $S(\mathbf{q},\ensuremath{\omega})$, the NMR relaxation rate
1/T1, and the magnetic susceptibility $\ensuremath{\chi}$ while changing the anisotropy in the exchange coupling constants, by using the dynamical mean-field theory based on a Majorana fermion representation. We describe the details of the methodology including the continuous-time quantum Monte Carlo method for computing dynamical spin correlations and the maximum entropy method for analytic continuation. We confirm that the combined method provides accurate results in a wide
T range including the region where the spins are fractionalized. We find that also in the anisotropic cases the system exhibits peculiar behaviors below the high-
T crossover whose temperature is comparable to the average of the exchange constants: $S(\mathbf{q},\ensuremath{\omega})$ shows an inelastic response at the energy scale of the averaged exchange constant,
1/T1 continues to grow even though the equal-time spin correlations are saturated and almost
T independent, and $\ensuremath{\chi}$ deviates from the Curie-Weiss behavior. In particular, when the exchange interaction in one direction is stronger than the other two, the dynamical quantities exhibit qualitatively different
T dependences from the isotropic case at low
T, reflecting the opposite parity between the flux-free ground state and the flux-excited state, and a larger energy cost for flipping a spin in the strong interaction direction. On the other hand, when the exchange anisotropy is in the opposite way, the results are qualitatively similar to those in the isotropic case. All these behaviors manifest the spin fractionalization in the paramagnetic region. Among them, the dichotomy between the static and dynamical spin correlations is unusual behavior hardly seen in conventional magnets. We discuss the relation between the dichotomy and the spatial configuration of the
Z2 gauge fields. Our results could stimulate further experimental and theoretical analyses of candidate materials for the Kitaev quantum spin liquids.