Using thermal tensor-network approach, we investigate the spin Seebeck effect (SSE) of the triangular lattice quantum antiferromagnet hosting spin supersolid phase. We focus on the low-temperature scaling behaviors of the normalized spin current across the interface. For the 1D Heisenberg chain, we find a negative spinon spin current in the bulk with algebraic temperature scaling; at low fields, boundary effects induce a second sign reversal at lower temperatures. These benchmark results are consistent with field-theoretical analysis. On the triangular lattice, spin frustration dramatically enhances the low-temperature SSE, with distinct spin-current signatures-particularly the sign reversal and characteristic temperature dependence-distinguishing different spin states. Remarkably, we discover a persistent, negative spin current in the spin supersolid phase, which saturates to a nonzero value in the low-temperature limit and can be ascribed to the Goldstone-mode-mediated spin supercurrent. Moreover, a universal scaling T^{d/z} is found at the U(1)-symmetric polarization quantum critical points. These distinct quantum spin transport traits provide sensitive spin current probes for spin supersolid states in quantum magnets such as Na_{2}BaCo(PO_{4})_{2}. Furthermore, our results also establish spin supersolids as a promising quantum platform for spin caloritronics in the ultralow-temperature regime.
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 Na
2BaNi(PO
4)
2 and K
2Ni(SeO
3)
2.
Low-dimensional quantum magnets, particularly those with strong spin frustration, are characterized by their notable spin fluctuations. Nuclear magnetic resonance (NMR) serves as a sensitive probe of low-energy fluctuations that offers valuable insight into rich magnetic phases and emergent phenomena in quantum magnets. Although experimentally accessible, the numerical simulation of NMR relaxation rates, specifically the spin-lattice relaxation rate
1/T1, remains a significant challenge. Analytical continuation based on Monte Carlo calculations are hampered by the notorious negative sign for frustrated systems, and the real-time simulations incur significant costs to capture low-energy fluctuations. Here we propose computing the relaxation rate using thermal tensor networks (TTNs), which provides a streamlined approach by calculating its imaginary-time proxy. We showcase the accuracy and versatility of our methodology by applying it to one-dimensional spin chains and two-dimensional lattices, where we find that the critical exponents
η and
zν can be extracted from the low-temperature scalings of the simulated
1/T1 near quantum critical points. Our results also provide insights into the low-dimensional and frustrated magnetic materials, elucidating universal scaling behaviors in the Ising chain compound CoNb
2O
6 and revealing the renormalized classical behaviors in the triangular-lattice antiferromagnet Ba
8CoNb
6O
24. We apply the approach to effective model of the family of frustrated magnets AYbCh
2 (A = Na, K, Cs, and Ch = O, S, Se), and find dramatic changes from spin ordered to the proposed quantum spin liquid phase. Overall, with high reliability and accuracy, the TTN methodology offers a systematic strategy for studying the intricate dynamics observed across a broad spectrum of quantum magnets and related fields.
Abstract Prototypical models and their material incarnations are cornerstones to the understanding of quantum magnetism. Here we show theoretically that the recently synthesized magnetic compound Na 2 BaCo(PO 4 ) 2 (NBCP) is a rare, nearly ideal material realization of the S = 1/2 triangular-lattice antiferromagnet with significant easy-axis spin exchange anisotropy. By combining the automatic parameter searching and tensor-network simulations, we establish a microscopic model description of this material with realistic model parameters, which can not only fit well the experimental thermodynamic data but also reproduce the measured magnetization curves without further adjustment of parameters. According to the established model, the NBCP hosts a spin supersolid state that breaks both the lattice translation symmetry and the spin rotational symmetry. Such a state is a spin analog of the long-sought supersolid state, thought to exist in solid Helium and optical lattice systems, and share similar traits. The NBCP therefore represents an ideal material-based platform to explore the physics of supersolidity as well as its quantum and thermal melting.
Based on the recent experiments on the triangular lattice antiferromagnet Na 2 BaCo ( PO 4 ) 2 , we propose the easy-axis XXZ spin-1/2 model on the triangular lattice, that exhibits a quantum spin supersolid, to be a precursory Dirac spin liquid. Despite the presence of a three-sublattice magnetic order as a spin supersolid, we suggest that this system is close to a Dirac spin liquid by exploring its spectroscopic response. The physical consequence is examined from the spectroscopic response, and we establish the continuous spectra near the M point in addition to the K point excitation from the spinon continuum on top of the three-sublattice order. Moreover, the satellite peaks were predicted at the mid-points connecting the Γ and K points. This proposal offers a plausible understanding of the recent inelastic neutron scattering measurement in Na 2 BaCo ( PO 4 ) 2 and could inspire further research in relevant models and materials, such as K 2 Co ( SeO 3 ) 2 and Rb 2 Co ( SeO 3 ) 2 , and even more anisotropic magnets like PrMgAl 11 O 19 . Published by the American Physical Society 2024
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 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.