Landau–Lifshitz–Gilbert equation

Landau–Lifshitz–Gilbert equation
Name	Landau–Lifshitz–Gilbert equation
Type	Partial differential equation
Field	Condensed matter physics
Discovered by	Lev Landau, Evgeny Lifshitz, T. L. Gilbert
Statement	Describes the precessional dynamics of magnetization in a ferromagnet.

Landau–Lifshitz–Gilbert equation. This fundamental equation of micromagnetics describes the time evolution of magnetization in ferromagnetic materials. It was first introduced by Lev Landau and Evgeny Lifshitz in 1935 to explain magnetic damping without involving thermodynamic concepts. The form commonly used today, incorporating a phenomenological damping term, was later developed by T. L. Gilbert in the 1950s, leading to the modern unified name.

Mathematical formulation

The standard form of the equation is a partial differential equation for the magnetization vector field \(\mathbf{M}(\mathbf{r}, t)\). It is given by: \[ \frac{\partial \mathbf{M}}{\partial t} = -\gamma \mathbf{M} \times \mathbf{H}_\text{eff} + \frac{\alpha}{M_s} \left( \mathbf{M} \times \frac{\partial \mathbf{M}}{\partial t} \right), \] where \(\gamma\) is the gyromagnetic ratio, \(\mathbf{H}_\text{eff}\) is the effective magnetic field, \(M_s\) is the saturation magnetization, and \(\alpha\) is the dimensionless Gilbert damping parameter. The effective field is derived from the functional derivative of the total magnetic energy density, which includes contributions from exchange interaction, magnetocrystalline anisotropy, Zeeman energy, and demagnetizing field. An algebraically equivalent form, often used in numerical simulations, eliminates the implicit time derivative on the right-hand side. This equation is central to the theoretical framework established by John H. Van Vleck and Louis Néel.

Physical interpretation

The equation describes the precessional motion of magnetization around the effective field, analogous to the Larmor precession of a magnetic moment in an external field. The first term represents this conservative torque, which conserves the magnitude of \(\mathbf{M}\). The second term, introduced by T. L. Gilbert, provides a phenomenological damping torque that drives the magnetization toward alignment with \(\mathbf{H}_\text{eff}\), dissipating energy. This damping is crucial for modeling realistic magnetic relaxation processes observed in experiments, such as those conducted at Bell Labs. The competition between precession and damping governs the dynamics leading to equilibrium, a concept also explored in the context of the Bloch equations for nuclear magnetization.

Derivation and relation to other equations

The equation can be derived from the Lagrangian mechanics of a continuous magnetic system or from the Landau-Lifshitz energy functional using variational principles. It is closely related to the Landau–Lifshitz equation, which uses a different damping form. In the limit of small damping (\(\alpha \ll 1\)), both forms are equivalent. The equation is also a specialization of the more general torque equation for angular momentum dynamics. Its connection to the Brown's equation in micromagnetics is foundational, as it provides the dynamic counterpart to William Fuller Brown Jr.'s static equilibrium theory. Comparisons are often made to the Schrödinger equation in quantum mechanics for its fundamental role in its domain.

Applications in magnetism

The primary application is in modeling the dynamics of magnetic domain walls and the switching of magnetic nanoparticles. It is essential for simulating the behavior of magnetic random-access memory (MRAM) devices and spin-transfer torque phenomena studied at institutions like IBM and Intel. The equation is used to analyze ferromagnetic resonance experiments, which probe material properties like damping. It also underpins the design of magnetic sensors and recording heads for hard disk drives, technologies advanced by companies such as Seagate Technology. Research into skyrmions and spintronics at facilities like the Max Planck Institute relies heavily on simulations using this equation.

Numerical solution and computational aspects

Solving the equation requires sophisticated numerical methods due to its nonlinearity and the stiffness of the associated partial differential equation. Common approaches include finite difference methods, such as those implemented in the OOMMF software package from the National Institute of Standards and Technology, and finite element methods. Explicit time-integration schemes like the Runge–Kutta methods are often used, but stability concerns necessitate implicit or semi-implicit methods, especially for stiff problems. High-performance computing resources at laboratories like Lawrence Berkeley National Laboratory are employed for large-scale micromagnetic simulations. The computational challenge involves accurately resolving fast precessional motion and long relaxation timescales simultaneously.

Limitations and extensions

The standard equation assumes a purely phenomenological, local damping form, which may not capture all microscopic dissipation mechanisms in materials like yttrium iron garnet. It does not inherently include effects from spin-polarized currents, which led to the development of the Landau–Lifshitz–Gilbert–Slonczewski equation incorporating spin-transfer torque. Further extensions include adding field-like torque terms and stochastic fields for thermal effects, formalized in the Stratonovich calculus for Brownian motion. For ultrafast demagnetization processes, observed with femtosecond laser pulses, coupling to other degrees of freedom via models like the Bloch–Bloembergen equations or atomistic spin dynamics simulations is necessary. Current research at CERN and Argonne National Laboratory explores these quantum-thermodynamic limits.

Category:Equations Category:Condensed matter physics Category:Magnetism