# Landau–Lifshitz–Gilbert equation

In physics, the Landau–Lifshitz–Gilbert equation, named for Lev Landau, Evgeny Lifshitz, and T. L. Gilbert, is a name used for a differential equation describing the precessional motion of magnetization M in a solid. It is a modification by Gilbert of the original equation of Landau and Lifshitz.

The various forms of the equation are commonly used in micromagnetics to model the effects of a magnetic field on ferromagnetic materials. In particular it can be used to model the time domain behavior of magnetic elements due to a magnetic field.[1] An additional term was added to the equation to describe the effect of spin polarized current on magnets.[2]

## Landau–Lifshitz equation

The terms of the Landau–Lifshitz–Gilbert equation: precession (red) and damping (blue). The trajectory of the magnetization (dotted spiral) is drawn under the simplifying assumption that the effective field Heff is constant.

In a ferromagnet, the magnetization M can vary internally but at each point its magnitude is equal to the saturation magnetization Ms. The Landau–Lifshitz–Gilbert equation predicts the rotation of the magnetization in response to torques. An earlier, but equivalent, equation (the Landau–Lifshitz equation) was introduced by Landau & Lifshitz (1935):[3][4][5]

${\displaystyle {\frac {d\mathbf {M} }{dt}}=-\gamma \mathbf {M} \times \mathbf {H_{\mathrm {eff} }} -\lambda \mathbf {M} \times \left(\mathbf {M} \times \mathbf {H_{\mathrm {eff} }} \right)}$

(1)

where γ is the electron gyromagnetic ratio. and λ is a phenomenological damping parameter, often replaced by

${\displaystyle \lambda =\alpha {\frac {\gamma }{M_{\mathrm {s} }}},}$

where α is a dimensionless constant called the damping factor. The effective field Heff is a combination of the external magnetic field, the demagnetizing field (magnetic field due to the magnetization), and some quantum mechanical effects. To solve this equation, additional equations for the demagnetizing field must be included.

Using the methods of irreversible statistical mechanics, numerous authors have independently obtained the Landau–Lifshitz equation.[6]

## Landau–Lifshitz–Gilbert equation

In 1955 Gilbert replaced the damping term in the Landau–Lifshitz (LL) equation by one that depends on the time derivative of the magnetization:

${\displaystyle {\frac {d\mathbf {M} }{dt}}=-\gamma \left(\mathbf {M} \times \mathbf {H} _{\mathrm {eff} }-\eta \mathbf {M} \times {\frac {d\mathbf {M} }{dt}}\right)}$

(2b)

This is the Landau–Lifshitz–Gilbert (LLG) equation, where η is the damping parameter, which is characteristic of the material. It can be transformed into the Landau–Lifshitz equation:[3]

${\displaystyle {\frac {d\mathbf {M} }{dt}}=-\gamma '\mathbf {M} \times \mathbf {H} _{\mathrm {eff} }-\lambda \mathbf {M} \times (\mathbf {M} \times \mathbf {H} _{\mathrm {eff} })}$

(2a)

where

${\displaystyle \gamma '={\frac {\gamma }{1+\gamma ^{2}\eta ^{2}M_{s}^{2}}}\qquad {\text{and}}\qquad \lambda ={\frac {\gamma ^{2}\eta }{1+\gamma ^{2}\eta ^{2}M_{s}^{2}}}.}$

In this form of the LL equation, the precessional term γ' depends on the damping term. This better represents the behavior of real ferromagnets when the damping is large.[7]

## Landau–Lifshitz–Gilbert–Slonczewski equation

In 1996 Slonczewski expanded the model to account for the spin-transfer torque, i.e. the torque induced upon the magnetization by spin-polarized current flowing through the ferromagnet. This is commonly written in terms of the unit moment defined by m = M / MS:

${\displaystyle {\dot {\mathbf {m} }}=-\gamma \mathbf {m} \times \mathbf {H} _{\mathrm {eff} }+\alpha \mathbf {m} \times {\dot {\mathbf {m} }}+\tau _{\parallel }{\frac {\mathbf {m} \times (\mathbf {x} \times \mathbf {m} )}{\left|\mathbf {x} \times \mathbf {m} \right|}}+\tau _{\perp }{\frac {\mathbf {x} \times \mathbf {m} }{\left|\mathbf {x} \times \mathbf {m} \right|}}}$

where ${\displaystyle \alpha }$ is the unitless damping parameter, ${\displaystyle \tau _{\perp }}$ and ${\displaystyle \tau _{\parallel }}$ are driving torques, and x is the unit vector along the polarization of the current.[8][9]

## References and footnotes

1. ^ Yang, Bo. "Numerical Studies of Dynamical Micromagnetics". Retrieved 8 August 2011.
2. ^
3. ^ a b Aharoni 1996
4. ^ Brown, Jr. 1978
5. ^ Chikazumi 1997
6. ^ T. Iwata, J. Magn. Magn. Mater. 31–34, 1013 (1983); T. Iwata, J. Magn. Magn. Mater. 59, 215 (1986); V.G. Baryakhtar, Zh. Eksp. Teor. Fiz. 87, 1501 (1984); S. Barta (unpublished, 1999); W. M. Saslow, J. Appl. Phys. 105, 07D315 (2009).
7. ^ For details of Kelly's non-resonant experiment, and of Gilbert's analysis (which led to Gilbert's modifying the damping term), see T. L. Gilbert and J. M. Kelly, "Anomalous rotational damping in ferromagnetic sheets", Conf. Magnetism and Magnetic Materials, Pittsburgh, PA, June 14–16, 1955 (New York: American Institute of Electrical Engineers, Oct. 1955, pp. 253–263). http://people.physics.tamu.edu/saslow/MMMConf55_253GilbertKelly.pdf Text references to Figures 5 and 6 should have been to Tables 1 and 2. Gilbert could not fit Kelly's experiments with fixed usual gyromagnetic ratio γ and a frequency-dependent λ=αγ, but could fit that data for a fixed Gilbert gyromagnetic ratio γG=γ/(1+α2) and a frequency-dependent α. Values of α as large as 9 were required, indicating very broad absorption, and thus a relatively low-quality sample. Modern samples, when analyzed from resonance absorption, give α's on the order of 0.05 or less.
8. ^ Slonczewski, John C. (1996). "Current-driven excitation of magnetic multilayers". Journal of Magnetism and Magnetic Materials. 159 (1): –1–L7. Bibcode:1996JMMM..159L...1S. doi:10.1016/0304-8853(96)00062-5. Retrieved 2014-09-25.
9. ^ Wolf, S. A. (16 November 2001). "Spintronics: A Spin-Based Electronics Vision for the Future". Science. 294 (5546): 1488–1495. Bibcode:2001Sci...294.1488W. doi:10.1126/science.1065389. PMID 11711666.