Motivated by recent experimental evidence for spin supersolid states in triangular-lattice compounds, we numerically investigate the dynamical properties of magnetic field-induced phases in the spin-1/2 easy-axis triangular antiferromagnetic Heisenberg model in the presence of magnetic impurities. In both weak- and strong-field spin supersolid states, the gapless Goldstone mode at the $K$ points remains robust against impurities, which is a direct manifestation of spin superfluidity. By contrast, at the same impurity density, impurities induce a splitting of the magnon bands in the conventional magnetic state, the so-called up-up-down state. In addition, the finite superfluid stiffness probed by the twisted phase in the spin supersolid states is consistent with the excitation spectrum. We argue that the excitation spectrum with impurities provides direct spectroscopic evidence for dissipationless spin dynamics in the spin supersolid states, which is experimentally accessible via inelastic neutron scattering.
Motivated by recent experimental evidence for spin supersolid states in triangular-lattice compounds, we numerically investigate the dynamical properties of magnetic field-induced phases in the spin-1/2 easy-axis triangular antiferromagnetic Heisenberg model in the presence of magnetic impurities. In both weak- and strong-field spin supersolid states, the gapless Goldstone mode at the K points remains robust against impurities, which is a direct manifestation of spin superfluidity. By contrast, at the same impurity density, impurities induce a splitting of the magnon bands in the conventional magnetic state, the so-called up-up-down state. In addition, the finite superfluid stiffness probed by the twisted phase in the spin supersolid states is consistent with the excitation spectrum. We argue that the excitation spectrum with impurities provides direct spectroscopic evidence for dissipationless spin dynamics in the spin supersolid states, which is experimentally accessible via inelastic neutron scattering.
We present a systematic numerical study of the spin-1 antiferromagnetic Heisenberg model on the triangular lattice in an out-of-plane magnetic field, using Density Matrix Renormalization Group (DMRG) methods. By mapping out the quantum phase diagram as a function of the single-ion anisotropy Dz and magnetic field, we identify distinct dipolar and quadrupolar spin supersolid states, characterized by spontaneous U(1) symmetry breaking with finite spin superfluid stiffness coexisting with longitudinal translational symmetry breaking. At zero field, the dipolar spin supersolid with a 'Y'-type spin configuration persists down to Dz=0, whereas the quadrupolar spin supersolid prevails at large Dz. At intermediate fields, the phase diagram is dominated by an up-up-down phase. At high fields below saturation, a quadrupolar spin superfluid emerges in the large-Dz regime, whereas a dipolar spin supersolid with a 'V'-type spin configuration dominates at small Dz. These phases are characterized through their order parameters and spin superfluid stiffness using calculations on various system sizes. Furthermore, the dynamical spin structure factor is obtained across the phase diagram, where characteristic spectral signatures of different phases are observed, including the gapless Goldstone mode and the roton-like minima. These features are directly accessible to inelastic neutron scattering experiments. Our results provide a theoretical understanding of the interplay between frustrations, anisotropy, and Zeeman interactions in driving distinct spin supersolid phases in the spin-1 system, which are relevant to various triangular-lattice antiferromagnets such as Na2BaNi(PO4)2 and K2Ni(SeO3)2.
We employ tensor network renormalization to explore the dynamical spectra of the easy-axis triangular-lattice antiferromagnet (TLAF) in a magnetic field. Our analysis identifies two distinct low-energy magnon excitations: a gapless Goldstone mode and a gapped mode. At zero field, the spectra display two nearly degenerate roton modes near the M point. With the increase of the magnetic field within the Y-shape superfluid phase, these modes diverge, with the roton excitation vanishing from the Goldstone mode branch, suggesting that the roton dip in this mode may just result from the energy-level repulsion imposed by the roton excitation in the gapped mode. Moreover, the in-plane spectral function shows substantial weight in high energies in the same spin excitation channel where the low-energy roton excitation appears. However, these roton excitations are absent in the V-shape supersolid phase.
We employ tensor-network renormalization to explore the dynamical spectra of the easy-axis triangular-lattice antiferromagnet (TLAF) in a magnetic field. Our analysis identifies two distinct low-energy magnon excitations: a gapless Goldstone mode and a gapped mode. At zero field, the spectra display two nearly degenerate roton modes near the M point. With the increase of the magnetic field within the Y-shaped superfluid phase, these modes diverge, with the roton excitation vanishing from the Goldstone mode branch, suggesting that the roton dip in this mode may just result from the energy-level repulsion imposed by the roton excitation in the gapped mode. Moreover, the in-plane spectral function shows substantial weight in high energies in the same spin excitation channel where the low-energy roton excitation appears. However, these roton excitations are absent in the V-shaped supersolid phase.
Motivated by recent experiments, we investigate the spin- 1 2 XXZ model on the triangular lattice with strong Ising anisotropy, combining large-scale numerical simulations and analytical methods to uncover unconventional spin dynamics at T = 0 . First, we compute the dynamical spin structure factor using density matrix renormalization group (DMRG) simulations and find excellent agreement with inelastic neutron scattering data on the layered compound K 2 Co ( SeO 3 ) 2 . The low-energy spectrum reveals a rotonlike minimum at the M point, absent in linear spin-wave theory, accompanied by peak intensity and a broad continuum above it. Near the Γ point, we observe an approximately linear dispersion with vanishing spectral weight. Second, we compare two analytical frameworks that reproduce the observed features. The first is a hard-core boson approach, which includes (i) an effective staggered boson model at zero magnetic field, (ii) perturbation theory applied to the one-third magnetization plateau, and (iii) a self-consistent mean-field Schwinger boson theory (SBT). The second framework is based on a variational supersolid quantum dimer model (QDM) ansatz, combined with a single-mode approximation. The SBT captures the broad continuum, the M -point minimum, and linear dispersion at Γ , whereas the QDM reproduces the roton minimum and linear dispersion at finite momentum near Γ . Remarkably, both the QDM wave function and the DMRG ground state exhibit nearly identical structure factors with pronounced transverse photonlike excitations. Together, our comprehensive theoretical and numerical analysis elucidates the microscopic origin of supersolid excitations in the XXZ triangular-lattice model and their proximity to a spin-liquid phase observed experimentally.
We investigate the dynamical signatures of the supersolid phase in the triangular-lattice hardcore Bose-Hubbard model via large-scale quantum Monte Carlo simulations and linear spin-wave theory. We reveal a distinct momentum-space separation in the spectral weights. The transverse spectrum shows gapless Goldstone modes at the Γ and K points, while the longitudinal spectrum features a gapless mode at K and a roton-like minimum at M. Directly in the imaginary-time domain, the supersolid is unambiguously characterized by the coexistence of transverse spectral lines at Γ and K. Furthermore, we show that strong quantum fluctuations induce a breakdown of the sharp quasiparticle picture, leading to asymmetric line shapes and broad continua. Our work establishes the definitive dynamical signatures of the hardcore-boson supersolid, providing key insights into analogous spin supersolids in triangular quantum magnets.
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.