Landau–Lifshitz–Gilbert equation

Template:Short description In physics, the Landau–Lifshitz–Gilbert equation (usually abbreviated as LLG equation), named for Lev Landau, Evgeny Lifshitz, and Thomas L. Gilbert, is a name used for a differential equation describing the dynamics (typically the precessional motion) of magnetization Template:Math in a solid. It is a modified version by Gilbert of the original equation of Landau and Lifshitz. The LLG equation is similar to the Bloch equation, but they differ in the form of the damping term. The LLG equation describes a more general scenario of magnetization dynamics beyond the simple Larmor precession. In particular, the effective field driving the precessional motion of Template:Math is not restricted to real magnetic fields; it incorporates a wide range of mechanisms including magnetic anisotropy, exchange interaction, and so on.

The various forms of the LLG equation are commonly used in micromagnetics to model the effects of a magnetic field and other magnetic interactions on ferromagnetic materials. It provides a practical way to model the time-domain behavior of magnetic elements. Recent developments generalizes the LLG equation to include the influence of spin-polarized currents in the form of spin-transfer torque.

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 H_eff is constant.

In a ferromagnet, the magnitude of the magnetization Template:Math at each spacetime point is approximated by the saturation magnetization Template:Math (although it can be smaller when averaged over a chunk of volume). The LLG equation describes the rotation of the magnetization in response to the effective field Template:Math and accounts for not only a real magnetic field but also internal magnetic interactions such as exchange and anisotropy. An earlier, but equivalent, equation (the Landau–Lifshitz equation) was introduced by Template:Harvtxt: Template:NumBlk where Template:Math is the electron gyromagnetic ratio and Template:Math is a phenomenological damping parameter, often replaced by

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

where Template:Math is a dimensionless constant called the damping factor. The effective field Template:Math is a combination of the external magnetic field, the demagnetizing field, and various internal magnetic interactions involving quantum mechanical effects, which is typically defined as the functional derivative of the magnetic free energy with respect to the local magnetization Template:Math. To solve this equation, additional conditions for the demagnetizing field must be included to accommodate the geometry of the material.

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: Template:NumBlk This is the Landau–Lifshitz–Gilbert (LLG) equation, where Template:Math is the damping parameter, which is characteristic of the material. It can be transformed into the Landau–Lifshitz equation: Template:NumBlk where

${\gamma }^{\prime }=\frac{\gamma }{1+{\gamma }^2{\eta }^2{M}_s^2}\text{and}\lambda =\frac{{\gamma }^2\eta }{1+{\gamma }^2{\eta }^2{M}_s^2}.$

In this form of the LL equation, the precessional term Template:Math depends on the damping term. This better represents the behavior of real ferromagnets when the damping is large.

Landau–Lifshitz–Gilbert–Slonczewski equation

In 1996 John 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_s$:

$\dot{𝐦}=-\gamma 𝐦\times {𝐇}_{\mathrm{e}\mathrm{f}\mathrm{f}}+\alpha 𝐦\times \dot{𝐦}+{\tau }_{\parallel }\frac{𝐦\times (𝐱\times 𝐦)}{|𝐱\times 𝐦|}+{\tau }_{\perp }\frac{𝐱\times 𝐦}{|𝐱\times 𝐦|}$

where $\alpha$ is the dimensionless damping parameter, ${\tau }_{\perp }$ and ${\tau }_{\parallel }$ are driving torques, and Template:Math is the unit vector along the polarization of the current.

Application to magnetic resonance

The response of a magnetic material to an oscillating magnetic field (e.g. the magnetic component of an electromagnetic wave) can be found using the LLG equation. Without loss of generality, the magnetization equilibrium magnetization vector ${𝐌}_s$, which is aligned with the static magnetic field ${𝐇}_s$ can be taken to be in the z-direction. One can then consider a small oscillating magnetic field $𝐡\propto e^{i\omega t}$, which causes a small change in the magnetization $\mathbf{\text{𝐦}}$. Thus, the total magnetization reads $𝐌={𝐌}_s+𝐦$, and the total magnetic field is $𝐇={𝐇}_s+𝐡$.

Filling this in to the LLG equation (neglecting the damping terms) gives

$\frac{d}{dt}({𝐌}_s+𝐦)=\gamma ({𝐌}_s+𝐦)\times ({𝐇}_s+𝐡)$.

The equilibrium magnetization does not vary over time. Additionally, the ${𝐌}_s\times {𝐇}_s$ term is zero because they are parallel. Finally, the equation is considered in leading order, so the $𝐦\times 𝐡$ term is neglected because it is the product of two small terms. Thus, the equation becomes $\frac{d𝐦}{dt}=\gamma (𝐦\times {𝐇}_s+{𝐌}_s\times 𝐡).$

The z-component of this equation is zero, meaning that the magnetization does not respond to an magnetic field oscillating in the z-direction. Assuming that $𝐦$ oscillates together with $𝐡$, the time derivative turns into $i\omega$. The x- and y-components can then be written in matrix form as

$\left(\begin{array}{cc}\gamma H_s& i\omega \\ -i\omega & \gamma H_s\end{array}\right)\left(\begin{array}{c}{m}_x\\ m_y\end{array}\right)=-\gamma M_s\left(\begin{array}{c}{h}_x\\ h_y\end{array}\right).$

Inverting this matrix gives the permeability tensor

$\left(\begin{array}{c}{m}_x\\ m_y\end{array}\right)=\frac{\gamma M_s}{{\gamma }^2{H}_s^2-{\omega }^2}\left(\begin{array}{cc}\gamma H_s& -i\omega \\ i\omega & \gamma H_s\end{array}\right)\left(\begin{array}{c}{h}_x\\ h_y\end{array}\right).$

Thus, there is a resonance at the ferromagnetic resonance frequency.

Microscopic origin of the damping

In the Landau–Lifshitz–Gilbert equation the damping is described by the phenomenological, dimensionless parameter Template:Mvar, which does not by itself specify the physical mechanism responsible for the dissipation of energy and angular momentum from the precessing magnetization. In conducting ferromagnets the dominant contribution to the intrinsic damping is generally attributed to spin–orbit coupling, which couples the magnetization to the itinerant electron system and allows energy to be transferred to the electronic degrees of freedom, where it is ultimately dissipated through electron scattering.

Kamberský mechanism (torque-correlation model)

The most widely used microscopic description of intrinsic Gilbert damping is the torque-correlation model (TCM), often referred to as the Kamberský mechanism after the Czech physicist Vladimír Kamberský, who developed the theory in a series of papers between 1970 and 1976. In this model the damping is expressed through a correlation function of the spin–orbit torque operator evaluated near the Fermi surface. When the magnetization is tilted slightly from equilibrium, spin–orbit coupling generates a non-equilibrium electron population that relaxes back towards equilibrium over a characteristic scattering time Template:Mvar, dissipating energy in the process.

The torque-correlation model can be regarded as a generalization of Kamberský's earlier breathing Fermi surface model, and it makes fewer approximations; in particular it describes damping over a wide range of temperatures and in the presence of chemical (alloy) disorder. A characteristic feature of the resulting expression is that it separates into two contributions with different dependence on the electronic scattering rate: an intraband (conductivity-like) term that grows as the scattering rate decreases, and an interband (resistivity-like) term that grows as the scattering rate increases. This partitioning has been used to interpret the non-monotonic temperature dependence of the damping observed in ferromagnetic metals such as iron, cobalt and nickel.

Kamberský later re-derived the damping formula using exact linear-response theory, expressing it in a form convenient for numerical estimates in terms of the Fermi golden rule.

Divergence problem and later developments

Calculations based on Kamberský's formula predict that the intraband contribution, which is of third order in the spin–orbit coupling parameter, diverges for a perfect crystal in the limit of zero temperature (vanishing scattering rate). Later analytical work argued that this divergence is an artifact of applying the formula beyond its regime of validity, showing that when the derivation is carried out consistently and restricted to second order in the spin–orbit parameter, no divergent intraband terms appear. In real materials the presence of defects, reflected in a finite residual resistivity, ensures that the damping parameter remains finite even as the temperature approaches zero.

The torque-correlation model has since been implemented in first-principles electronic structure calculations and extended in several directions, including tensorial and non-local formulations of the damping. It has been complemented by an equivalent formulation based on scattering theory, which relates the damping to the scattering matrix of the ferromagnet and can avoid explicit reliance on an external relaxation-time parameter.

Further reading

• Template:Cite journal This is only an abstract; the full report is "Armor Research Foundation Project No. A059, Supplementary Report, May 1, 1956", but was never published. A description of the work is given in Template:Cite journal The original abstract reads: T. L. Gilbert, Armour Research Foundation of Illinois Institute of Technology.–The gyromagnetic equation, $d𝐌/dt=\gamma 𝐌\times 𝒢$, for the motion of the magnetization field $𝐌\mathbf{(}𝐫\mathbf{)}$, in a ferromagnetic material can be derived from a variation principle, as first shown by Doering (W. Doering, Z. Naturforsch. 3a, 374 (1948)). Here $𝒢$ is the effective internal field, including the magnetic field and contributions from exchange, anisotropy, and magnetoelastic effects. Using the variational principle, the equations of motion can be recast into a Lagrangian form. This makes possible a consistent derivation of the equations of motion of the magnetization field and other fields to which it may be coupled (e.g., the displacement field of the lattice and the electromagnetic field). It also permits the introduction of viscous damping effects in a consistent manner using the Rayleigh dissipation function. It is shown that viscous damping of the magnetization fields leads to an equation of motion which reduces to the Landau-Lifshitz equation only when the damping is small. It is also shown that this Lagrangian formalism permits the introduction of damping due to disaccomodation in a consistent and very general way.
• Template:Cite journal
• Template:Cite journal
• Template:Cite book
• Template:Cite journal
• Template:Cite journal