---
title: Atomistic Spin-Lattice Dynamics (ASLD)
url: https://www.emergentmind.com/topics/atomistic-spin-lattice-dynamics-asld
type: topic
---

# Atomistic Spin-Lattice Dynamics (ASLD)

Atomistic Spin-Lattice Dynamics (ASLD) is an atomistic framework in which magnetic moments or spins and lattice degrees of freedom are evolved together in time, rather than treated as separate or adiabatically decoupled subsystems. In the literature, ASLD appears both as a microscopic electronically mediated theory—where spin-spin, lattice-lattice, and spin-lattice couplings are obtained after integrating out the electrons—and as a practical simulation methodology that combines atomistic spin dynamics with molecular dynamics or Langevin lattice dynamics using first-principles, phenomenological, or data-driven parameterizations [1505.08005, 1804.03119]. Its central purpose is to resolve feedback between magnetic fluctuations, phonons, strain, and defects, including regimes where non-adiabatic scattering, magneto-elastic renormalization, or angular-momentum exchange are essential [1612.08503, 1802.02934].

## 1. Conceptual foundations

A microscopic formulation of ASLD treats localized magnetization $\mathbf M_i$ and ionic displacement $\mathbf Q_i$ on the same footing. In the electronically mediated theory, the starting point is a local exchange coupling between electron spin and magnetic moment together with a local coupling between electronic charge and lattice displacement. After integrating out the electrons, one obtains an effective action containing three interaction classes: spin-spin, lattice-lattice, and spin-lattice couplings [1505.08005]. In discrete form, the resulting minimal bilinear Hamiltonian is

$$
\mathcal H_{\text{MQ}}= -\frac{1}{2}\sum_{ij}\Bigl(\mathbf Q_i\cdot[T^{cc}_{ij}\cdot \mathbf Q_j + T^{cs}_{ij}\cdot \mathbf M_j] + \mathbf M_i\cdot[T^{sc}_{ij}\cdot \mathbf Q_j + T^{ss}_{ij}\cdot \mathbf M_j]\Bigr).
$$

This construction is notable because the spin-lattice term is not introduced as a post hoc correction. It is an electronically mediated susceptibility, alongside the more familiar interatomic force constants and spin-spin interactions, and it is explicitly constrained by time-reversal and inversion symmetry [1505.08005].

A computationally efficient first-principles version of ASLD reformulates the same idea as a coupled classical dynamics problem. In that setting, atomistic spin dynamics provides the spin evolution, molecular dynamics provides the atomic motion, and a spin-lattice Hamiltonian couples the two through displacement-dependent magnetic interactions [1804.03119]. This formulation generalizes separate ASD and MD treatments into a single framework that can be parameterized from density functional theory rather than fitted phenomenologically [1804.03119].

An important conceptual revision came from work showing that the conventional combined molecular and spin dynamics framework, with exchange interactions depending on atomic positions, still lacks a channel for transferring angular momentum between spins and lattice. The proposed remedy is to encode spin-orbit physics through local magnetic anisotropies induced by symmetry breaking from phonons or defects, thereby extending ASLD from an energy-coupled model to one that also supports angular-momentum exchange [1612.08503].

## 2. Hamiltonian structure and coupling mechanisms

A widely used decomposition writes the total Hamiltonian as

$$
\mathcal H_{\rm SLD}=\mathcal H_{\rm S}+\mathcal H_{\rm L}+\mathcal H_{\rm LS},
$$

where $\mathcal H_{\rm S}$ is the spin Hamiltonian, $\mathcal H_{\rm L}$ the lattice Hamiltonian, and $\mathcal H_{\rm LS}$ the spin-lattice coupling [1804.03119]. In first-principles ASLD for solids, the bilinear magnetic term is expanded in atomic displacements, producing a hierarchy of couplings. The leading term is the first-order spin-lattice or exchange-striction term,

$$
\mathcal H_{\rm SSL} = -\frac{1}{2}\sum_{ijk} A_{ijk}\cdot \mathbf{u}_k\,(\mathbf{m}_i\cdot \mathbf{m}_j),
$$

with $A_{ijk}$ obtained from derivatives of the exchange tensor with respect to displacement [1804.03119]. This makes explicit the reciprocal relation between magnetism and structure: atomic motion modulates magnetic exchange, and spin correlations generate forces on atoms.

A more direct realization of the same idea appears in paramagnetic CrN, where the magnetic subsystem is represented by a classical Heisenberg model with distance-dependent exchange,

$$
\mathcal{H} = -\sum_{i\neq j} J_{ij}(R_{ij})\, \hat{\mathbf{S}}_i \cdot \hat{\mathbf{S}}_j,
$$

and the dependence of $J_{ij}$ on the Cr–Cr distance $R_{ij}$ is the essential magneto-lattice coupling channel [1802.02934]. In that material, the exchange is short-ranged and can be captured by first- and second-neighbor Cr–Cr terms [1802.02934].

Not all coupling channels are equivalent. A spin-orbit-inspired route introduces local anisotropy terms that emerge from symmetry breaking in the atomic environment. In the anisotropy-enhanced framework, the additional contribution is

$$
\mathcal{H}_{\text{anis}} = - C_1 \sum_{i=1}^N \mathbf{K}_i \cdot \mathbf{S}_i - C_2 \sum_{i=1}^N \mathbf{S}_i^\intercal \cdot \mathbf{\Lambda}_i \cdot \mathbf{S}_i,
$$

where $\mathbf K_i$ and $\mathbf \Lambda_i$ are determined by the instantaneous local atomic environment through a symmetry descriptor $\rho_i$ [1612.08503]. A related microscopic coupling used in coupled spin and lattice dynamics is a two-site pseudo-dipolar term motivated by spin-orbit interaction,

$$
\mathcal{H}_c = -\sum_{i,j} f(r_{ij}) \left[ (\mathbf{S}_i\cdot \hat{\mathbf{r}}_{ij}) (\mathbf{S}_j\cdot \hat{\mathbf{r}}_{ij}) -\frac13 \mathbf{S}_i\cdot\mathbf{S}_j \right],
$$

chosen so that the Curie temperature is preserved under changes in coupling strength, net anisotropy is removed on a cubic lattice, and a spurious uniform translation is avoided when the magnetization is saturated [2010.00642].

## 3. Equations of motion and computational realization

The coupled dynamics are typically written as simultaneous equations for positions, velocities, and spins. In one representative formulation,

$$
\frac{\partial \mathbf{r}_i}{\partial t}=\mathbf{v}_i,\qquad
\frac{\partial \mathbf{v}_i}{\partial t}=-\eta \mathbf{v}_i+\frac{\mathbf{F}_i}{m_i},\qquad
\frac{\partial \mathbf{S}_i}{\partial t}=-\gamma \mathbf{S}_i\times \mathbf{H}_i,
$$

with forces and effective fields obtained from derivatives of the total Hamiltonian with respect to positions and spins [2010.00642]. Other implementations replace the purely precessional spin equation by stochastic Landau–Lifshitz–Gilbert dynamics and use Langevin terms for both lattice and spin subsystems [2101.07332].

A first-principles realization of fully coupled ASLD in paramagnetic CrN alternates atomistic spin dynamics with ab initio molecular dynamics. The workflow begins with pre-equilibration by adiabatic DLM-AIMD, then computes distance-dependent $J_{ij}$, generates an initial spin configuration by Monte Carlo, and iterates a coupled loop in which the current spins are used in a Born–Oppenheimer AIMD step for $1$ fs, the new atomic positions update $J_{ij}(R_{ij})$, and the spins are then propagated for $1$ fs using ASD [1802.02934]. The force calculations use spin-polarized, noncollinear DFT in VASP with constrained spin directions through the Ma–Dudarev constrained-moment method with Lagrange multipliers [1802.02934].

The numerical integrator depends on the coupling structure. A Suzuki–Trotter decomposition is used in self-consistent spin-lattice dynamics with pseudo-dipolar coupling, including a Cayley transform for norm-conserving spin updates [2010.00642]. When nonlinear anisotropy terms are added, the conventional Suzuki–Trotter scheme becomes inapplicable in its standard form, and a hybrid integrator combining Suzuki–Trotter decomposition with the iterative method of Krech et al. is used; in that case the timestep is reduced from $1$ fs to $0.1$ fs to maintain microcanonical energy conservation [1612.08503]. Data-driven large-scale ASLD for $\alpha$-Fe is implemented in LAMMPS using the SPIN package and a massively parallel symplectic ASLD algorithm, with the total potential energy surface decomposed into a collective atomic spin model and a SNAP-based machine-learning interatomic potential [2101.07332].

## 4. Non-adiabaticity, damping, and angular-momentum exchange

A recurrent simplification in early coupled spin-lattice models is the assumption that coordinate-dependent exchange alone is sufficient to mediate full relaxation. The literature summarized here shows that this is not generally the case. In the conventional MD-SD Hamiltonian, exchange depending on positions allows energy transfer between lattice and spins, but because of rotational symmetry it does not allow the lattice and spin subsystems to exchange angular momentum. A heat bath connected only to the lattice can therefore thermalize the atoms while leaving the spins trapped [1612.08503]. The anisotropy-based extension remedies this by introducing symmetry-breaking torques that preserve energy, linear momentum, and total angular momentum while allowing angular momentum to move between spin and lattice sectors [1612.08503].

The Einstein–de Haas effect provides a direct benchmark. In simulations of a freely rotating prolate spheroidal Fe nanocluster containing $505$ Fe atoms, a modified version of SPILADY with the dynamic anisotropy correction of Perera et al. shows that the model conserves total angular momentum while spin and lattice angular momenta are not conserved separately. The rate of angular momentum transfer is proportional to the anisotropy strength $C_2$, and the added anisotropy allows full spin-lattice relaxation on a timescale of approximately $100$ ps [2211.03706].

ASLD also supplies a microscopic route to effective damping without imposing a phenomenological Gilbert term on the spin subsystem. In the pseudo-dipolar coupling model, the effective damping parameter increases with both coupling strength and temperature; the reported dependence on coupling is approximately quadratic, and at low temperature damping values on the order of $\alpha \sim 10^{-5}$ are obtained [2010.00642]. On the spin side, a generalized atomistic equation with a damping tensor and a moment-of-inertia tensor,

$$
\dot\mathbf{M} = \mathbf{M}\times\left( -\gamma\mathbf{B} +\hat{\mathbf G}\dot\mathbf{M} +\hat{\mathbf I}\ddot\mathbf{M} \right),
$$

extends standard atomistic spin dynamics into the femtosecond regime and provides a first-principles foundation for damping and nutational effects, although it does not itself constitute a full spin-lattice model [1110.3387].

## 5. Materials results and experimental observables

One of the clearest demonstrations of ASLD’s non-adiabatic content is paramagnetic CrN above the magnetic ordering temperature. In fully coupled ASD-AIMD, the dynamic coupling between spin fluctuations and lattice vibrations strongly broadens phonon spectra at $300\,\mathrm{K}$ relative to the adiabatic DLM-AIMD control, whereas at $1000\,\mathrm{K}$ the difference is much smaller. The effect is strongest for acoustic phonons, and correlation analysis shows that antiferromagnetically aligned nearest-neighbor Cr pairs tend to have shorter distances while ferromagnetically aligned pairs tend to have longer distances; this correlation is strong at $300\,\mathrm{K}$ and weaker at $1000\,\mathrm{K}$ [1802.02934]. The result is a microscopic explanation of the anomalous thermal conductivity of paramagnetic CrN: near the transition, dynamic magneto-lattice coupling shortens acoustic phonon lifetimes, while farther above the transition the coupling weakens and the extra spin-induced phonon scattering diminishes [1802.02934].

In bcc Fe, first-principles ASLD based on a displacement expansion of the exchange tensor produces coupling-induced modifications of both magnon and phonon spectra, and the same framework exhibits dissipation-free coupled motion in small magnetic clusters such as dimers, trimers, and quadmers [1804.03119]. A data-driven extension for $\alpha$-iron couples a collective spin Hamiltonian with a SNAP machine-learning interatomic potential and yields quantitative predictions of bulk modulus, magnetization, and specific heat across the ferromagnetic-paramagnetic phase transition. Under fixed-volume conditions the Curie temperature is predicted near $716$ K, while pressure-controlled and pressure- and magnetization-controlled protocols recover thermal softening and the specific-heat anomaly more faithfully [2101.07332].

ASLD has also been propagated into electron microscopy observables. In body-centered cubic iron at $300$ K, an extension of TACAW to ASLD produces momentum-resolved EELS from unified finite-temperature trajectories containing both atomic displacements and spin configurations. The resulting “Full” spectrum is not the sum of separate “Phonons (ASLD)” and “Magnons (ASLD)” signals; instead it shows interference and energy redistribution in $(\mathbf q,E)$-space. With an annular dark-field detector of inner and outer collection angles $2$ and $7$ mrad, phonons dominate below $100$ meV, magnetic intensity becomes visible above approximately $100$ meV, and energy-resolved EELS improves the signal-to-noise ratio by roughly three orders of magnitude relative to a no-energy-resolution configuration [2508.07073]. This establishes that phonon and magnon excitations in the coupled system are not generally additive observables.

## 6. Software, scaling, and methodological boundaries

ASLD practice spans full spin-lattice codes, spin-only frameworks that supply the magnetic subsystem, and accelerator-oriented infrastructure. UppASD is used as the ASLD engine in coupled phonon-magnon EELS simulations, where it supplies unified spin-lattice trajectories through coupled Langevin equations [2508.07073]. SPILADY has been modified to include dynamic anisotropy for Einstein–de Haas simulations in Fe nanoclusters [2211.03706]. LAMMPS, through the SPIN package and a massively parallel symplectic ASLD algorithm, supports data-driven magneto-elastic simulations with machine-learning interatomic potentials [2101.07332]. Spirit provides the atomistic spin side of the ecosystem, with an extended Heisenberg-type spin-lattice Hamiltonian, Monte Carlo, LLG solvers, geodesic nudged elastic band, minimum mode following, harmonic transition-state theory, and CPU/GPU parallelization, but it does not propagate atomic coordinates in the cited work [1901.11350]. SpinX supplies a JAX-based atomistic spin framework with dense, sparse, FFT-based, pair-list, and reciprocal-space dipolar backends, plus deterministic and stochastic LLG, Monte Carlo, spectroscopy, and string/GNEB calculations; performance benchmarks report peak throughput exceeding $10$ billion spin-site operations per second on a single accelerator and aggregate single-node workloads of over $1$ billion atomic spins [2605.26016].

Post-processing and spectroscopy are themselves major computational bottlenecks. For real-space dynamical spin-spin correlations, GPU kernel fusion reformulates the pair-correlation construction as GEMM, yielding up to $39\times$ speedup over a baseline GPU code and up to $25\times$ over a $64$-core CPU GEMM implementation. Fusing inner-product work into the GEMM kernel with CUTLASS improves performance by $26\%\text{–}33\%$ over a cuBLAS-plus-Thrust implementation, while on-the-fly epilogue computation of the phase matrix enables larger systems; on four A100 GPUs, reported speedups range from $2.9\times$ to $3.6\times$ over one GPU [2308.07487].

Several methodological boundaries recur across the literature. Most implementations remain semiclassical, with classical spins and nuclei and electrons integrated out [1804.03119]. Phenomenological ingredients persist where first-principles extraction is difficult, particularly for local anisotropy parameters and some spin-orbit-mediated couplings [1612.08503]. Fixed-magnitude spin approximations can limit accuracy at strong spin disorder or elevated pressure [2101.07332]. Explicit quantum spin effects are commonly neglected, although path-integral approaches have been proposed to encode quantum corrections as effective anisotropy-like fields within ASD-style sampling [2303.00602]. Three simplifications are shown to fail in specific contexts: isotropic exchange alone does not enable spin-lattice angular-momentum transfer [1612.08503], adiabatic magnetic disorder misses non-adiabatic phonon lifetime broadening in CrN [1802.02934], and separate phonon and magnon spectra need not add to the full signal in EELS [2508.07073]. Together, these results define ASLD less as a single algorithm than as a class of coupled atomistic theories and simulation strategies designed to resolve the mutual dynamics of magnetic and structural degrees of freedom.

Source: https://www.emergentmind.com/topics/atomistic-spin-lattice-dynamics-asld