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.
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.
The Landau-Lifshitz-Gilbert (LLG) and Landau-Lifshitz (LL) equations play an essential role for describing the dynamics of magnetization in solids. While a quantum analog of the LL dynamics has been proposed in [Phys. Rev. Lett. 110, 147201 (2013)PRLTAO0031-900710.1103/PhysRevLett.110.147201], the corresponding quantum version of LLG remains unknown. Here, we propose such a quantum LLG equation that inherently conserves purity of the quantum state. We examine the quantum LLG dynamics of a dimer consisting of two interacting spin-1/2 particles. Our analysis reveals that, in the case of ferromagnetic coupling, the evolution of initially uncorrelated spins mirrors the classical LLG dynamics. However, in the antiferromagnetic scenario, we observe pronounced deviations from classical behavior, underscoring the unique dynamics of becoming a spinless state, which is nonlocally correlated. Moreover, when considering spins that are initially entangled, our study uncovers an unusual form of revival-type quantum correlation dynamics, which differs significantly from what is typically seen in open quantum systems.
The Landau-Lifshitz-Gilbert (LLG) equation provides the fundamental description of dissipative magnetization dynamics. A recently proposed quantum analog, the quantum Landau-Lifshitz-Gilbert (QLLG) equation, extends this framework to quantum spin systems. However, whether the QLLG equation recovers the LLG equation in the classical large-spin limit for interacting spin systems remains an open question. Here, we demonstrate this connection by deriving the LLG equation as the classical large-spin limit of the QLLG equation. We establish this correspondence analytically and verify it numerically for an interacting spin system. These results provide a quantum-to-classical description of dissipative spin dynamics, offering a microscopic foundation for the LLG equation. From numerical calculations provided here, we also conclude that quantum effects become much less noticeable for systems with larger number of spins in the Hamiltonian.
Understanding the microscopic origin of magnetic damping is essential for describing magnetization dynamics in nanoscale magnetic systems. In this work, we develop a microscopic framework to derive Gilbert damping terms from specific interactions in ferromagnetic conductors and magnetic insulators. In ferromagnetic conductors, magnetic damping is dominated by the s d exchange interaction between the localized magnetic moments and the itinerant electrons, which is effectively equivalent to an electronmagnon interaction. From the spin diffusion equation and the spin transfer torque therein, we obtain the corresponding quantum LLG equation. In magnetic insulators, where free electrons are absent, spin phonon interactions become the primary source of the damping.Starting from an effective spin phonon interaction and using linear response theory, we derive a quantum LLG equation with a Gilbert damping tensor analytically. Our results demonstrate that different microscopic interactions lead to qualitatively different structures of magnetic damping. Instead of being described by a universal scalar parameter, Gilbert damping can exhibit spatial dependence and tensorial characteristics determined by the underlying microscopic interactions.
Ground-state preparation of quantum many-body systems is a fundamental challenge across quantum science, with implications for quantum simulation, quantum computing, and the study of complex quantum matter. The exponential growth of Hilbert space and the complexity of entangled eigenstates motivate the development of new approaches that go beyond conventional strategies based on adiabatic evolution, variational optimization, imaginary-time methods, and engineered dissipative processes. Here, we demonstrate that the recently proposed quantum Landau-Lifshitz-Gilbert (QLLG) dynamics [Phys. Rev. Lett. {\bf 133}, 266704 (2024)] provides an intrinsic nonlinear dissipative mechanism for real-time ground-state preparation. Unlike engineered dissipative schemes, QLLG encodes relaxation directly into the quantum equation of motion while preserving the underlying quantum structure of the evolution. The QLLG evolution selectively suppresses excited-state contributions, driving the system toward the lowest-energy eigenstate within the accessible symmetry sector. We show that the dynamics admits a monotonic energy functional and a stable ground-state fixed point, and that for generic initial states in an N-qubit Hilbert space the characteristic convergence time scales linearly with N and inversely with the spectral gap. Numerical simulations of an interacting spin chain confirm our analytical predictions. These results establish nonlinear dissipative QLLG dynamics as a new avenue for quantum state preparation, opening a new route for quantum simulation, computation, and state engineering.
The entanglement of many localized spins (LSs) within solid magnetic materials is a topic of great basic and applied interest, particularly after becoming amenable to experimental scrutiny where recent neutron scattering experiments have witnessed macroscopic entanglement in the ground state (GS) of antiferromagnets persisting even at elevated temperatures. On the other hand, spintronics and magnonics studies assume that LSs of antiferromagnets are in unentangled N\'eel GS, as well as that they evolve, when pushed out of equilibrium by current or external fields, according to the Landau-Lifshitz-Gilbert (LLG) equation viewing LSs as classical vectors of fixed length. The prerequisite for applicability of the LLG equation is zero entanglement in the underlying many-body quantum state of LSs. In this study, we initialize quantum Heisenberg ferro- or antiferromagnetic chains hosting S=1/2,S=1, or S=5/2 LSs into an unentangled pure state and then evolve them by quantum master equations (QMEs) of Lindblad or non-Markovian type, derived by coupling LSs weakly to a bosonic bath (due to phonons in real materials) or by using additional ``reaction coordinate'' in the latter case. The time evolution is initiated by applying an external magnetic field, and entanglement of the ensuing mixed quantum states is monitored by computing its negativity. We find that non-Markovian dynamics never brings entanglement to zero, in the presence of which the vector of spin expectation value changes its length to render the LLG equation inapplicable. Conversely, Lindbladian (i.e., Markovian) dynamics ensures that entanglement goes to 0, thereby enabling quantum-to-classical transition in all cases---S=1/2,S=1, and S=5/2 ferromagnets or S=5/2 antiferromagnets---except for S=1/2 and S=1 antiferromagnets. Finally, we investigate the stability of an entangled antiferromagnetic GS upon suddenly coupling it to bosonic baths.
We examine the relation between the quantum Landau-Lifshitz equation (q-LL) [Phys. Rev. Lett. 110, 147201 (2013)] and quantum Landau-Lifshitz-Gilbert equation (q-LLG) [Phys. Rev. Lett. 133, 266704 (2024)]; two non-linear purity preserving master equations that extend classical atomistic spin dynamics into the quantum regime. While the classical LL and LLG counterparts for any number of spins are known to be equivalent, i.e., give identical spin trajectories up to a rescaling of the time parameter, the quantum formulations are equivalent only in certain cases, such as for pure states or for arbitrary single spin-21 states. Here, we demonstrate that this equivalence breaks down even at the level of a single spin, provided s≥1. Focusing on a spin-1 particle in an anisotropic crystal field, we show that the q-LL and q-LLG equations generate inequivalent time evolution. We introduce temporal rescaling misfits that quantify the inequivalence of the two types of dynamics. Although our results highlight fundamental differences in dissipation mechanisms encoded in these equations, the resulting trajectories remain qualitatively similar for this system.