Investigation of dynamical excitations is difficult but crucial to the understanding of many exotic quantum phenomena discovered in quantum materials. This is particularly true for highly frustrated quantum antiferromagnets whose dynamical properties deviate strongly from theoretical predictions made based on the spin-wave or other approximations. Here, we present a large-scale numerical calculation on the dynamical correlation functions of spin-1/2 triangular Heisenberg model using a state-of-the-art tensor network renormalization group method. The calculated results allow us to gain for the first time a comprehensive picture on the nature of spin excitation spectra in this highly frustrated quantum system. It provides a quantitative account for all the key features of the dynamical spectra disclosed by inelastic neutron scattering measurements for Ba_{3}CoSb_{2}O_{9}, revealing the importance of the interplay between low- and high-energy excitations and its renormalization effect to the low-energy magnon bands and high-energy continuums. We identify the longitudinal Higgs modes in the intermediate-energy scale and predict the energy and momentum dependence of spectral functions along the three principal axes that can be verified by polarized neutron scattering experiments. Furthermore, we find that the spin excitation spectra weakly depend on the anisotropic ratio of the antiferromagnetic interaction.
Investigation of dynamical excitations is difficult but crucial to the understanding of many exotic quantum phenomena discovered in quantum materials. This is particularly true for highly frustrated quantum antiferromagnets whose dynamical properties deviate strongly from theoretical predictions made based on the spin-wave or other approximations. Here, we present a large-scale numerical calculation on the dynamical correlation functions of spin-1/2 triangular Heisenberg model using a state-of-the-art tensor network renormalization group method. The calculated results allow us to gain for the first time a comprehensive picture on the nature of spin excitation spectra in this highly frustrated quantum system. It provides a quantitative account for all the key features of the dynamical spectra disclosed by inelastic neutron scattering measurements for Ba_{3}CoSb_{2}O_{9}, revealing the importance of the interplay between low- and high-energy excitations and its renormalization effect to the low-energy magnon bands and high-energy continuums. We identify the longitudinal Higgs modes in the intermediate-energy scale and predict the energy and momentum dependence of spectral functions along the three principal axes that can be verified by polarized neutron scattering experiments. Furthermore, we find that the spin excitation spectra weakly depend on the anisotropic ratio of the antiferromagnetic interaction.
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 investigate the quantum dynamics of the spins on different Heisenberg antiferromagnetic spin lattice systems. Firstly, we applied the coupled-cluster method to the spin-1/2 antiferromagnetic \emph{XXZ} model on a square lattice by employing an approximation which contains two-body long-range correlations and high-order four-body local correlations. Improvement is found for the ground-state energy, sublattice magnetization, and the critical anisotropy when comparing with the approximation including the two-body correlations alone. We also obtain the full excitation spectrum which is in good agreement with the quantum Monte Carlo results and the high-order spin-wave theory.Secondly, we study the longitudinal excitations of quantum antiferromagnets on a triangular lattice by a recently proposed microscopic many-body approach based on magnon-density waves. We calculate the full longitudinal excitation spectra of the antiferromagnetic Heisenberg model for a general spin quantum number in the isotropic limit. Similar to the square lattice model, we find that, at the center of the first hexagonal Brillouin zone Γ(q=0) and at the magnetic ordering wavevectors ±[Q=(4π/3,0)], the excitation spectra become gapless in the thermodynamic limit, due to the slow, logarithmic divergence of the structure factor. However, these longitudinal modes on two-dimensional models may be considered as quasi-gapped, as any finite-size effect or small anisotropy will induce a large energy gap, when compared with the counterpart of the transverse spin-wave excitations. We have also investigated the excited states of the quasi-one-dimensional quantum antiferromagnets on hexagonal lattices, including the longitudinal modes based on the magnon-density waves. A model Hamiltonian with a uniaxial single-ion anisotropy is first studied by a spin-wave theory based on the one-boson method; the ground state thus obtained is employed for the study of the longitudinal modes. The full energy spectra of both the transverse modes (i.e., magnons) and the longitudinal modes are obtained as functions of the nearest-neighbor coupling and the anisotropy constants. We have found two longitudinal modes due to the non-collinear nature of the triangular antiferromagnetic order, similar to that of the phenomenological field theory approach by Affleck. The excitation energy gaps due to the anisotropy and the energy gaps of the longitudinal modes without anisotropy are then investigated. We then compares our results for the longitudinal energy gaps at the magnetic wavevectors with the experimental results for several antiferromagnetic compounds with both integer and non-integer spin quantum numbers, and we find good agreements after the higher-order contributions are included in our calculations.
Neutron scattering is central to identifying quantum states of magnetic materials. In the search for quantum spin liquids, broad spectral features of inelastic spectra have been cited as evidence for spinon excitations, but can also arise from magnon excitations excitations in the presence of quenched disorder and strong magnon interactions. We develop a new approach to this problem, based on the adiabatic continuity in the XXZ Heisenberg model on geometrically frustrating (GF) lattices as a function of the model's anisotropy. Using this approach, we identify universal features and energies of finite-temperature spin correlators. Focusing on the kagome lattice, we show that the low-energy spin spectral function contains robust, momentum-independent peaks with frequencies: ω1≈3.4T∗ and ω2≈6.3T∗, where the ``hidden energy scale'' T∗ is the characteristic scale of a low-temperature peak in the heat capacity, at which many GF magnets also display spin-glass freezing. We show that the spectral features at low energies ω≲T∗ arise from single-magnon scattering and identify the magnetizations of the respective excitations. We explore the evolution of the spectral features with temperature and discuss extensions to other GF lattices. Our results provide a sharp spectroscopic criterion for interpreting neutron scattering in kagome and other GF quantum magnets.
Abstract Fifty years after Anderson’s resonating valence-bond proposal, the spin-1/2 triangular-lattice Heisenberg antiferromagnet (TLHAF) remains the ultimate platform to explore highly entangled quantum spin states in proximity to magnetic order. Yb-based delafossites are ideal candidate TLHAF materials, which allow experimental access to the full range of applied in-plane magnetic fields. We perform a systematic neutron scattering study of CsYbSe 2 , first proving the Heisenberg character of the interactions and quantifying the second-neighbor coupling. We then measure the complex evolution of the excitation spectrum, finding extensive continuum features near the 120°-ordered state, throughout the 1/3-magnetization plateau and beyond this up to saturation. We perform cylinder matrix-product-state (MPS) calculations to obtain an unbiased numerical benchmark for the TLHAF and spectacular agreement with the experimental spectra. The measured and calculated longitudinal spectral functions reflect the role of multi-magnon bound and scattering states. These results provide valuable insight into unconventional field-induced spin excitations in frustrated quantum materials.
In ordered quantum magnets where interactions between elementary excitations dominate over their kinetic energy, perturbative approaches often fail, making non-perturbative methods essential to capture spectral features such as bound states and the redistribution of weight within excitation continua. Although an increasing number of experiments report anomalous spin excitation continua in such systems, their microscopic interpretation remains an open challenge. Here, we investigate the spin dynamics of the triangular-lattice antiferromagnet in its 1/3-plateau phase using two complementary non-perturbative approaches: exact diagonalization in a truncated Hilbert space for a gas of elementary excitations (THED) and matrix product state (MPS) simulations. Alongside cross-validation between these methods, we benchmark our results against inelastic neutron scattering (INS) data. The THED analysis confirms the presence of two-magnon bound states and identifies the anomalous scattering continuum observed in both MPS and INS as a two-magnon resonance, arising from hybridization between the bound state and the two-magnon continuum. Furthermore, THED reveals bound states overlapping with the continuum, enriching the interpretation of continuum anomalies. More broadly, THED provides a robust framework for investigating anomalous spin excitation continua and bound-state effects in other materials with gapped spectra. Its combination of accuracy and computational efficiency makes it a powerful tool for extracting reliable microscopic models in semiclassical regimes.
Abstract In ordered quantum magnets where interactions between elementary excitations dominate over their kinetic energy, perturbative approaches often fail, making non-perturbative methods essential to capture spectral features such as bound states and the redistribution of weight within excitation continua. Although an increasing number of experiments report anomalous spin excitation continua in such systems, their microscopic interpretation remains an open challenge. Here, we investigate the spin dynamics of the triangular-lattice antiferromagnet in its 1/3-plateau phase using two complementary non-perturbative approaches: exact diagonalization in a truncated Hilbert space for a gas of elementary excitations (THED) and matrix product state (MPS) simulations. Alongside cross-validation between these methods, we benchmark our results against inelastic neutron scattering (INS) data. The THED analysis confirms the presence of two-magnon bound states and identifies the anomalous scattering continuum observed in both MPS and INS as a two-magnon resonance, arising from hybridization between the bound state and the two-magnon continuum. Furthermore, THED reveals bound states overlapping with the continuum, enriching the interpretation of continuum anomalies. More broadly, THED provides a robust framework for investigating anomalous spin excitation continua and bound-state effects in other materials with gapped spectra. Its combination of accuracy and computational efficiency makes it a powerful tool for extracting reliable microscopic models in semiclassical regimes.