Multimagnon dynamics and thermalization in the $S=1$ easy-axis ferromagnetic chain

TL;DR

This paper investigates the dynamics, interactions, and thermalization of multimagnon states in an S=1 ferromagnetic chain with easy-axis anisotropy, combining numerical simulations and analytical insights relevant for experiments and ultracold atomic systems.

Contribution

It provides a detailed analysis of multimagnon behavior, including binding energies and nonequilibrium dynamics, using advanced numerical methods and analytical approximations.

Findings

Magnon binding energy decreases with the number of magnons.

Predictions for ultracold atomic experiments on magnon thermalization.

Observation of behaviors similar to many-body quantum scars.

Abstract

Quasiparticles are physically motivated mathematical constructs for simplifying the seemingly complicated many-body description of solids. A complete understanding of their dynamics and the nature of the effective interactions between them provides rich information on real material properties at the microscopic level. In this work, we explore the dynamics and interactions of magnon quasiparticles in a ferromagnetic spin-1 Heisenberg chain with easy-axis onsite anisotropy, a model relevant for the explanation of recent terahertz optics experiments on NiNb$_2$O$_6$ [P. Chauhan et al., Phys. Rev. Lett. 124, 037203 (2020)],and nonequilibrium dynamics in ultracold atomic settings [W.C. Chung et al., Phys. Rev. Lett. 126, 163203 (2021)]. We build a picture for the properties of clouds of a few magnons with the help of exact diagonalization and density matrix renormalization group calculations supported by physically motivated Jastrow wavefunctions. We show how the binding energy of magnons effectively reduces with their number and explain how this energy scale is of direct relevance for dynamical magnetic susceptibility measurements. This understanding is used to make predictions for ultracold-atomic platforms which are ideally suited to study the thermalization of multimagnon states. We simulate the nonequilibrium dynamics of these chains using the matrix product state based time-evolution block decimation algorithm and explore the dependence of revivals and thermalization on magnon density and easy-axis onsite anisotropy (which controls the strength of effective magnon interactions). We observe behaviors akin to those reported for many-body quantum scars which we explain with an analytic approximation that is accurate in the limit of small anisotropy.