Search NASA⌕ Search

SEARCH · Search NASA

Results for “Landau–Lifshitz–Gilbert dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

11 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Proximity Spin–Orbit Torque on a Two-Dimensional Magnet within van der Waals Heterostructure: Current-Driven Antiferromagnet-to-Ferromagnet Reversible Nonequilibrium Phase Transition in Bilayer CrI 3

The recently discovered two-dimensional magnetic insulator CrI 3 is an intriguing case for basic research and spintronic applications since it is a ferromagnet in the bulk but an antiferromagnet in bilayer form, with its magnetic ordering amenable to external manipulations. Using the first-principles quantum transport approach, we predict that injecting unpolarized charge current parallel to the interface of the bilayer-CrI 3 /monolayer-TaSe 2 van der Waals (vdW) heterostructure will induce spin–orbit torque and thereby drive the dynamics of magnetization on the first monolayer of CrI 3 in direct contact with TaSe 2 . By combining the calculated complex angular dependence of spin–orbit torque with the Landau-Lifshitz-Gilbert equation for classical dynamics of magnetization, we demonstrate that current pulses can switch the direction of magnetization on the first monolayer to become parallel to that of the second monolayer, thereby converting CrI 3 from antiferromagnet to ferromagnet while not requiring any external magnetic field. We explain the mechanism of this reversible current-driven nonequilibrium phase transition by showing that first monolayer of CrI 3 carries current due to evanescent wave functions injected by metallic transition metal dichalcogenide TaSe 2 , while concurrently acquiring strong spin–orbit coupling via such a proximity effect, whereas the second monolayer of CrI 3 remains insulating. Furthermore, the transition can be detected by passing vertical read current through the vdW heterostructure, encapsulated by a bilayer of hexagonal boron nitride and sandwiched between graphite electrodes, where we find a tunneling magnetoresistance of ≃240%.

2D magnetic materials↗
DOE OSTI2025

OpenFerro v0.1.0

OpenFerro is a Python package for on-lattice atomistic dynamics simulation of ferroic materials. OpenFerro is based on JAX, a high-performance linear algebra package supporting auto-differentiation and GPU acceleration. OpenFerro is designed to minimize the effort required to build on-lattice Hamiltonian models, and to perform molecular dynamics (MD) and Landau-Lifshitz-Gilbert simulations. Unlike existing codes, OpenFerro provides a unified interface to model different types of local order parameters.

Xie, Pinchen [Lawrence Berkeley National Laborator↗
DOE OSTI2026

Machine-learning modeling of magnetization dynamics in quasi-equilibrium and driven metallic spin systems

Here, we present a perspective on recent progress in machine-learning (ML) force-field approaches for large-scale Landau–Lifshitz–Gilbert (LLG) simulations of metallic spin systems. Building on a generalization of the Behler–Parrinello (BP) architecture originally developed for quantum molecular dynamics, we develop scalable and transferable ML models that faithfully capture the complex, environment-dependent electron-mediated exchange fields characteristic of itinerant magnets. A central ingredient of this framework is the implementation of symmetry-aware magnetic descriptors based on group-theoretical bispectrum formalisms. Leveraging these ML force fields, LLG simulations faithfully reproduce hallmark non-collinear magnetic orders—such as the 120° and tetrahedral states—on the triangular lattice, and successfully capture the complex spin textures emerging in the mixed-phase states of a square-lattice double-exchange model under thermal quench. We further discuss a generalized potential theory that extends the BP formalism to incorporate both conservative and nonconservative electronic torques, thereby enabling ML models to learn nonequilibrium exchange fields from computationally demanding microscopic approaches such as nonequilibrium Green’s-function techniques. This extension yields quantitatively accurate predictions of voltage-driven domain-wall motion and establishes a foundation for quantum-accurate, multiscale modeling of nonequilibrium spin dynamics and spintronic functionalities.

Descriptors↗
DOE OSTI2021

Quantum Spin Torque Driven Transmutation of an Antiferromagnetic Mott Insulator

The basic model of spin-transfer torque (STT) in antiferromagnetic spintronics considers the exchange of angular momentum between quantum spins of flowing electrons and noncollinear-to-them localized spins treated as classical vectors. These vectors are assumed to realize Néel order in equilibrium, ↑ ↓ ∙∙∙ ↑ ↓, and their STT-driven dynamics is described by the Landau-Lifshitz-Gilbert (LLG) equation. However, many experimentally employed materials (such as archetypal NiO) are strongly electron-correlated antiferromagnetic Mott insulators (AFMIs) whose localized spins form a ground state quite different from the unentangled Néel state | ↑ ↓ ∙∙∙ ↑ ↓$\rangle$. The true ground state is entangled by quantum spin fluctuations, leading to the expectation value of all localized spins being zero, so that LLG dynamics of classical vectors of fixed length rotating due to STT cannot even be initiated. Instead, a fully quantum treatment of both conduction electrons and localized spins is necessary to capture the exchange of spin angular momentum between them, denoted as quantum STT. We use a recently developed time-dependent density matrix renormalization group approach to quantum STT to predict how injection of a spin-polarized current pulse into a normal metal layer coupled to an AFMI overlayer via exchange interaction and possibly small interlayer hopping—mimicking, e.g., topological-insulator/NiO bilayer employed experimentally—will induce a nonzero expectation value of AFMI localized spins. This new nonequilibrium phase is a spatially inhomogeneous ferromagnet with a zigzag profile of localized spins. The total spin absorbed by AFMI increases with electron-electron repulsion in AFMIs, as well as when the two layers do not exchange any charge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Linearized frequency domain Landau-Lifshitz-Gilbert equation formulation

We present a general finite element linearized Landau-Lifshitz-Gilbert equation (LLGE) solver for magnetic systems under weak time-harmonic excitation field. The linearized LLGE is obtained by assuming a small deviation around the equilibrium state of the magnetic system. Inserting such expansion into LLGE and keeping only first order terms gives the linearized LLGE, which gives a frequency domain solution for the complex magnetization amplitudes under an external time-harmonic applied field of a given frequency. We solve the linear system with an iterative solver using generalized minimal residual method. We construct a preconditioner matrix to effectively solve the linear system. The validity, effectiveness, speed, and scalability of the linear solver are demonstrated via numerical examples.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Machine-learning force-field models for dynamical simulations of metallic magnets

We review recent advances in machine-learning (ML) force-field methods for Landau–Lifshitz–Gilbert simulations of itinerant electron magnets, focusing on their scalability and transferability. Built on the principle of locality, a deep neural-network model is developed to efficiently and accurately predict electron-mediated forces governing spin dynamics. Symmetry-aware descriptors constructed through a group-theoretical approach ensure rigorous incorporation of both lattice and spin-rotation symmetries. The framework is demonstrated using the prototypical s-d exchange model widely employed in spintronics. ML-enabled large-scale simulations reveal novel nonequilibrium phenomena, including anomalous coarsening of tetrahedral spin order on the triangular lattice and the freezing of phase-separation dynamics in lightly hole-doped, strong-coupling square-lattice systems. These results establish ML force-field frameworks as scalable, accurate, and versatile tools for modeling nonequilibrium spin dynamics in itinerant magnets.

Artificial neural networks↗
DOE OSTI2023

Control of magnon–photon coupling by spin torque

We demonstrate the influence of damping and field-like torques in the magnon–photon coupling process by classically integrating the generalized Landau–Lifshitz–Gilbert equation with RLC equation in which a phase correlation between dynamic magnetization and microwave current through combined Ampère and Faraday effects are considered. We show that the gap between two hybridized modes can be controlled in samples with damping parameter in the order of 10 –3 by changing the direction of the dc current density J if a certain threshold is reached. Furthermore, our results suggest that an experimental realization of the proposed magnon–photon coupling control mechanism is feasible in yttrium iron garnet/Pt hybrid structures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Skyrmions Near Defects

We study the impact of an exchange-reducing defect on a skyrmion for a thin film of finite thickness. Attraction of the skyrmion to a defect is demonstrated in a lattice model by computing the micromagnetic energy accounting for the exchange, Dzyaloshinskii-Moriya interaction, magnetic anisotropy, and dipole-dipole coupling. The spiraling dynamics of the skyrmion towards the defect is illustrated by solving numerically the full Landau-Lifshitz-Gilbert equations on a lattice and, independently, the Thiele equation, with the two methods in agreement with each other. We find that the skyrmion deforms when approaching the defect. Here, deformation is small in the parameter space far from the phase boundary that determines stability of skyrmions. It increases dramatically near the phase boundary, leading to the transformation of the skyrmion by the defect into a snake-like magnetic domain.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Ultrafast Light-Induced Magnetoelectric Effect in van der Waals Magnetic Semiconductor Heterostructures

Atomic-scale heterostructures of van der Waals (vdW) magnets and semiconductors provide a unique environment for exploring magnetic dynamics. In contrast to typical photothermal excitation of precessional magnetization dynamics by a pump laser pulse, we find that ultrafast optical excitation of a WS2/CrGeTe3 (CGT) bilayer produces an opposite sign of magnetic torque compared to an isolated CGT film. Experimental observations by time-resolved magneto-optic Kerr effect (TR-MOKE) and theoretical analysis by density functional theory (DFT) and Landau-Lifshitz-Gilbert (LLG) simulations support a mechanism in which charge transfer of photoexcited carriers across the interface alters the perpendicular magnetic anisotropy, which in turn generates a torque on the magnetic layer to trigger precessional magnetization dynamics. These results provide new avenues for ultrafast manipulation of magnetization in vdW heterostructures with type-II band alignments. Lastly, we show that optically-generated spin currents from WS2 into CGT can also trigger precessional dynamics via angular momentum transfer.

Zhou, Wenyi [The Ohio State University]↗
DOE OSTI2020

Domain wall diode based on functionally graded Dzyaloshinskii–Moriya interaction

In this work, we present a general approach for studying the dynamics of domain walls in biaxial ferromagnetic stripes with functionally graded Dzyaloshinskii–Moriya interaction (DMI). By engineering the spatial profile of the DMI parameter we propose the concept of a diode, which implements filtering of domain walls of certain topological charge and helicity. We base our study on phenomenological Landau–Lifshitz–Gilbert equations with additional Zhang–Li spin-transfer terms using a collective variable approach. In the effective equations of motion the gradients of DMI play the role of a driving force which competes with current driving. All analytical predictions are confirmed by numerical simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Coarsening dynamics of Ising-nematic order in a frustrated Heisenberg antiferromagnet

We study the phase ordering dynamics of the classical antiferromagnetic 𝐽 1 −𝐽 2 (nearest-neighbor and next-nearest-neighbor couplings) Heisenberg model on the square lattice in the strong frustration regime (𝐽 2 /𝐽 1 > 1/2). While thermal fluctuations preclude any long-range magnetic order at finite temperatures, the system exhibits a long-range spin-driven nematic phase at low temperatures. The transition into the nematic phase is further shown to belong to the two-dimensional Ising universality class based on the critical exponents near the phase transition. Our large-scale stochastic Landau-Lifshitz-Gilbert simulations find a two-stage phase ordering when the system is quenched from a high-temperature paramagnetic state into the nematic phase. In the early stage, collinear alignments of spins lead to a locally saturated Ising-nematic order. Once domains of well-defined Ising order are developed, the late-stage relaxation is dominated by curvature-driven domain coarsening, as described by the Allen-Cahn equation. The characteristic size of Ising-nematic domains scales as the square root of time, similar to the kinetic Ising model described by the time-dependent Ginzburg-Landau theory. Our results confirm that the late-stage ordering kinetics of the spin-driven nematic, which is a vestigial order of the frustrated Heisenberg model, belongs to the dynamical universality class of a nonconserved Ising order. Interestingly, the system shows no violation of the superuniversality hypothesis under weak bond disorder. The dynamic scaling invariance is preserved in the presence of weak bond disorder. Here, we also discuss possible applications of our results to materials for which vestigial Ising-nematic order is realized.

Antiferromagnets↗