Numerical solution of quantum Landau-Lifshitz-Gilbert equation

The classical Landau–Lifshitz–Gilbert (LLG) equation has long served as a cornerstone for modeling magnetization dynamics in magnetic systems, yet its classical nature limits its applicability to inherently quantum phenomena such as entanglement and nonlocal correlations. Inspired by the need to incorporate quantum effects into spin dynamics, recently a quantum generalization of the LLG equation is proposed [Phys. Rev. Lett. 133, 266704 (2024)] which captures essential quantum behavior in many-body systems. In this work, we develop a robust numerical methodology tailored to this quantum LLG framework that not only handles the complexity of quantum many-body systems but also preserves the intrinsic mathematical structures and physical properties dictated by the equation. We apply the proposed method to a class of quantum systems with a moderate number of spins that host host topological states of matter, and demonstrate rich quantum behavior, including the emergence of long-time entangled states. This approach opens a pathway toward reliable simulations of quantum magnetism beyond classical approximations, potentially leading to new discoveries.