Nonperturbative semiclassical spin dynamics for ordered quantum magnets

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.