We use the matrix product state formalism to construct stationary scattering states of elementary excitations in generic one-dimensional quantum lattice systems. Our method is applied to the spin-1 Heisenberg antiferromagnet, for which we calculate the full magnon-magnon S matrix for arbitrary momenta and spin, the two-particle contribution to the spectral function, and higher order corrections to the magnetization curve. As our method provides an accurate microscopic representation of the interaction between elementary excitations, we envisage the description of low-energy dynamics of one-dimensional spin chains in terms of these particlelike excitations.
We use the matrix product state formalism to construct stationary scattering states of elementary excitations in generic one-dimensional quantum lattice systems. Our method is applied to the spin-1 Heisenberg antiferromagnet, for which we calculate the full magnon-magnon S matrix for arbitrary momenta and spin, the two-particle contribution to the spectral function, and higher order corrections to the magnetization curve. As our method provides an accurate microscopic representation of the interaction between elementary excitations, we envisage the description of low-energy dynamics of one-dimensional spin chains in terms of these particlelike excitations.
Interactions between elementary excitations in quasi-one-dimensional antiferromagnets are of experimental relevance and their quantitative theoretical treatment has been a theoretical challenge for many years. Using matrix product states, one can explicitly determine the wave functions of the one- and two-particle excitations, and, consequently, the contributions to dynamical correlations. We apply this framework to the (nonintegrable) frustrated dimerized spin-1/2 chain, a model for generic spin-Peierls systems, where low-energy quasiparticle excitations are bound states of topological solitons. The spin structure factor involving two quasiparticle scattering states is obtained in the thermodynamic limit with full momentum and frequency resolution. This allows very subtle features in the two-particle spectral function to be revealed which, we argue, could be seen, e.g., in inelastic neutron scattering of spin-Peierls compounds under a change of the external pressure.
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.
Using a variational ansatz based on matrix product states, the authors present a methodology for determine the elementary excitations and their low energy scattering matrices of general quasi-one-dimensional spin systems. From these scattering matrices, the authors map the system unto an effective integrable model and then employ the thermodynamic Bethe ansatz so as to determine thermodynamic properties of the spin system. As a demonstration of the method, the authors apply it to a Heisenberg spin ladder.
We introduce a simple and efficient variation of the tangent-space excitation Ansatz used to compute elementary excitation spectra of one-dimensional quantum lattice systems using matrix product states (MPS). A small basis for the excitation tensors is formed based on a single diagonalization analogous to a single site DMRG step but for multiple states. Once overlap and Hamiltonian matrix elements are found, obtaining the excitation for any momentum only requires diagonalization of a tiny matrix, akin to a non-orthogonal band-theory diagonalization. The approach is based on an infinite MPS description of the ground state, and we introduce an extremely simple alternative to variational uniform matrix product states (VUMPS) based on finite system DMRG. For the S=1 S = 1 Heisenberg chain, our method—site basis excitation Ansatz (SBEA)—efficiently produces the one-magnon dispersion with high accuracy. We also examine the role of MPS gauge choices, finding that not imposing a gauge condition—leaving the basis nonorthogonal—is crucial for the approach, whereas imposing a left-orthonormal gauge (as in prior work) severely hampers convergence. We also show how one can construct Wannier excitations, analogous to the Wannier functions of band theory, where one Wannier excitation, translated to all sites, can reconstruct the single magnon modes exactly for all momenta.
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.