World Library  
Flag as Inappropriate
Email this Article


Article Id: WHEBN0004153112
Reproduction Date:

Title: Micromagnetics  
Author: World Heritage Encyclopedia
Language: English
Subject: Magnetic domain, Landau–Lifshitz–Gilbert equation, Magnetostatics, Magnetism, Magnonics
Publisher: World Heritage Encyclopedia


Micromagnetics is a field of physics dealing with the prediction of magnetic behaviors at sub-micrometer length scales. The length scales considered are large enough for the atomic structure of the material to be ignored (the continuum approximation), yet small enough to resolve magnetic structures such as domain walls or vortices.

Micromagnetics can deal with static equilibria, by minimizing the magnetic energy, and with dynamic behavior, by solving the time-dependent dynamical equation.


  • History 1
  • Static micromagnetics 2
    • Exchange energy 2.1
    • Anisotropy energy 2.2
    • Zeeman energy 2.3
    • Energy of the demagnetizing field 2.4
    • Magnetoelastic Energy 2.5
  • Dynamic micromagnetics 3
    • Effective field 3.1
    • Landau-Lifshitz-Gilbert equation 3.2
  • Applications 4
  • See also 5
  • Footnotes and references 6
  • Further reading 7
  • External links 8


Micromagnetics as a field (i.e., that deals specifically with the behaviour of (ferro)magnetic materials at sub-micrometer length scales) was introduced in 1963 when William Fuller Brown, Jr. published a paper on antiparallel domain wall structures. Until comparatively recently computational micromagnetics has been prohibitively expensive in terms of computational power, but smaller problems are now solvable on a modern desktop PC.

Static micromagnetics

The purpose of static micromagnetics is to solve for the spatial distribution of the magnetization M at equilibrium. In most cases, as the temperature is much lower than the Curie temperature of the material considered, the modulus |M| of the magnetization is assumed to be everywhere equal to the saturation magnetization Ms. The problem then consists in finding the spatial orientation of the magnetization, which is given by the magnetization direction vector m = M/Ms, also called reduced magnetization.

The static equilibria are found by minimizing the magnetic energy,

E = E_\text{exch} + E_\text{anis} + E_\text{Z} + E_\text{demag}+E_\text{m-e},

subject to the constraint |M|=Ms or |m|=1.

The contributions to this energy are the following:

Exchange energy

The exchange energy is a phenomenological continuum description of the quantum-mechanical exchange interaction. It is written as:

E_\text{exch} = A \int_V \left((\nabla m_x)^2 + (\nabla m_y)^2 + (\nabla m_z)^2\right) \mathrm{d}V

where A is the exchange constant; mx, my and mz are the components of m; and the integral is performed over the volume of the sample.

The exchange energy tends to favor configurations where the magnetization varies only slowly across the sample. This energy is minimized when the magnetization is perfectly uniform.

Anisotropy energy

Magnetic anisotropy arises due to a combination of crystal structure and spin-orbit interaction. It can be generally written as:

E_\text{anis} = \int_V F_\text{anis}(\mathbf{m}) \mathrm{d}V

where Fanis, the anisotropy energy density, is a function of the orientation of the magnetization. Minimum-energy directions for Fanis are called easy axes.

Time-reversal symmetry ensures that Fanis is an even function of m. The simplest such function is

F_\text{anis}(\mathbf{m}) = -K m_z^2.

where K is called the anisotropy constant. In this approximation, called uniaxial anisotropy, the easy axis is the z direction.

The anisotropy energy favors magnetic configurations where the magnetization is everywhere aligned along an easy axis.

Zeeman energy

The Zeeman energy is the interaction energy between the magnetization and any externally applied field. It's written as:

E_\text{Z} = -\mu_0 \int_V \mathbf{M}\cdot\mathbf{H}_\text{a} \mathrm{d}V

where Ha is the applied field and µ0 is the vacuum permeability.

The Zeeman energy favors alignment of the magnetization parallel to the applied field.

Energy of the demagnetizing field

Example of micromagnetic configuration. Compared to a uniform state, the flux closure structure lowers the energy of the demagnetizing field, at the expense of some exchange energy.

The demagnetizing field is the magnetic field created by the magnetic sample upon itself. The associated energy is:

E_\text{demag} = -\frac{\mu_0}{2} \int_V \mathbf{M}\cdot\mathbf{H}_\text{d} \mathrm{d}V

where Hd is the demagnetizing field. This field depends on the magnetic configuration itself, and it can be found by solving:

\nabla\cdot\mathbf{H}_\text{d} = -\nabla\cdot\mathbf{M}
\nabla\times\mathbf{H}_\text{d} = 0

where −∇·M is sometimes called magnetic charge density. The solution of these equations (c.f. magnetostatics) is:

\mathbf{H}_\text{d} = -\frac{1}{4\pi} \int_V \nabla\cdot\mathbf{M} \frac{\mathbf{r}}{r^3} \mathrm{d}V

where r is the vector going from the current integration point to the point where Hd is being calculated.

It is worth noting that the magnetic charge density can be infinite at the edges of the sample, due to M changing discontinuously from a finite value inside to zero outside of the sample. This is usually dealt with by using suitable boundary conditions on the edge of the sample.

The energy of the demagnetizing field favors magnetic configurations that minimize magnetic charges. In particular, on the edges of the sample, the magnetization tends to run parallel to the surface. In most cases it is not possible to minimize this energy term at the same time as the others. The static equilibrium then is a compromise that minimizes the total magnetic energy, although it may not minimize individually any particular term.

Magnetoelastic Energy

The magnetoelastic energy describes the energy storage due to elastic lattice distortions. It may be neglected if magnetoelastic coupled effects are neglected. There exists a preferred local distortion of the crystalline solid associated with the magnetization director m, . For a simple model, on can assume this strain to be isochoric and fully isotropic in the lateral direction, yielding the deviatoric ansatz

\mathbf{\varepsilon}_0(\mathbf{m}) = \frac{3}{2} E\, [\mathbf{m}\otimes \mathbf{m} - \frac{1}{3}\mathbf{1}]

where the material parameter E > 0 is the magnetostrictive constant. Clearly, E is the strain induced by the magnetization in the direction m. With this ansatz at hand, we consider the elastic energy density to be a function of the elastic, stress-producing strains \mathbf{\varepsilon}_e := \mathbf{\varepsilon} -\mathbf{\varepsilon}_0. A quadratic form for the magnetoelastic energy is

E_\text{m-e} = \frac{1}{2} [\mathbf{\varepsilon} -\mathbf{\varepsilon}_0(\mathbf{m})] : \mathbb{C} : [\mathbf{\varepsilon} -\mathbf{\varepsilon}_0(\mathbf{m})]

where \mathbb{C} :=\lambda \mathbf{1}\otimes \mathbf{1} + 2\mu \mathbb{I} is the fourth-order elasticity tensor. Here the elastic response is assumed to be isotropic (based on the two Lamé constants λ and μ). Taking into account the constant length of m, we obtain the invariant-based representation

E_\text{m-e}= \frac{\lambda}{2} \mbox{tr}^2[\mathbf{\varepsilon}] + \mu \, \mbox{tr}[\mathbf{\varepsilon}^2] - 3\mu E \big\{ \mbox{tr}[\mathbf{\varepsilon}(\mathbf{m}\otimes\mathbf{m})] - \frac{1}{3}\mbox{tr}[\mathbf{\varepsilon}] \big\} .

This energy term contributes to magnetostriction.

Dynamic micromagnetics

The purpose of dynamic micromagnetics is to predict the time evolution of the magnetic configuration of a sample subject to some non-steady conditions such as the application of a field pulse or an AC field. This is done by solving the Landau-Lifshitz-Gilbert equation, which is a partial differential equation describing the evolution of the magnetization in term of the local effective field acting on it.

Effective field

The effective field is the local field felt by the magnetization. It can be described informally as the derivative of the magnetic energy density with respect to the orientation of the magnetization, as in:

\mathbf{H}_\mathrm{eff} = - \frac{1}{\mu_0 M_s} \frac{\mathrm{d}^2 E}{\mathrm{d}\mathbf{m}\mathrm{d}V}

where dE/dV is the energy density. In variational terms, a change dm of the magnetization and the associated change dE of the magnetic energy are related by:

\mathrm{d}E = -\mu_0 M_s \int_V (\mathrm{d}\mathbf{m})\cdot\mathbf{H}_\text{eff}\,\mathrm{d}V

It should be noted that, since m is a unit vector, dm is always perpendicular to m. Then the above definition leaves unspecified the component of Heff that is parallel to m. This is usually not a problem, as this component has no effect on the magnetization dynamics.

From the expression of the different contributions to the magnetic energy, the effective field can be found to be:

\mathbf{H}_\mathrm{eff} = \frac{2A}{\mu_0 M_s} \nabla^2 \mathbf{m} - \frac{1}{\mu_0 M_s} \frac{\partial F_\text{anis}}{\partial \mathbf{m}} + \mathbf{H}_\text{a} + \mathbf{H}_\text{d}

Landau-Lifshitz-Gilbert 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.

This is the equation of motion of the magnetization. It describes a Larmor precession of the magnetization around the effective field, with an additional damping term arising from the coupling of the magnetic system to the environment. The equation can be written in the so-called Gilbert form (or implicit form) as:

\frac{\partial \mathbf m}{\partial t} = - |\gamma| \mathbf{m} \times \mathbf{H}_\mathrm{eff} + \alpha \mathbf{m}\times\frac{\partial \mathbf{m}} {\partial t}

where γ is the electron gyromagnetic ratio and α the Gilbert damping constant.

It can be shown that this is mathematically equivalent to the following Landau-Lifshitz (or explicit) form:

\frac{\partial\mathbf m}{\partial t} = - \frac{|\gamma|}{1+\alpha^2} \mathbf{m} \times \mathbf{H}_\mathrm{eff} - \frac{\alpha|\gamma|}{1+\alpha^2} \mathbf{m}\times(\mathbf{m}\times\mathbf{H}_\text{eff})


The interaction of micromagnetics with mechanics is also of interest in the context of industrial applications that deal with magneto-acoustic resonance such as in hypersound speakers, high frequency magnetostrictive transducers etc. FEM simulations taking into account the effect of magnetostriction into micromagnetics are of importance. Such simulations use models described above within a finite element framework.[1]

Apart from conventional magnetic domains and domain-walls, the theory also treats the statics and dynamics of topological line and point configurations, e.g. magnetic vortex and antivortex states;[2] or even 3d-Bloch points,[3][4] where, for example, the magnetization leads radially into all directions from the origin, or into topologically equivalent configurations. Thus in space, and also in time, nano- (and even pico-)scales are used.

The corresponding topological quantum numbers[4] are thought to be used as information carriers, to apply the most recent, and already studied, propositions in information technology.

See also

Footnotes and references

  1. ^ Miehe, Christian; Ethiraj, Gautam (2011-10-15). "A geometrically consistent incremental variational formulation for phase field models in micromagnetics". Computer Methods in Applied Mechanics and Engineering (Elsevier). 245–246: 331–347.  
  2. ^ Komineas, S. N. (2007). "Dynamics of vortex-antivortex pairs in ferromagnets". arXiv:0712.36841 [cond-mat.mtrl-sci].
  3. ^ Thiaville, André; García, José; Dittrich, Rok; Miltat, Jacques; Schrefl, Thomas (March 2003). "Micromagnetic study of Bloch-point-mediated vortex core reversal". Physical Review B 67 (9).  
  4. ^ a b Döring, W. (1968). "Point Singularities in Micromagnetism". Journal of Applied Physics 39 (2): 1006.  

Further reading

  • Brown, William Fuller, Jr. (1963). Micromagnetics. New York: Wiley.  
  • Gilbert, Thomas L. (2004). "A Phenomenological Theory of Damping in Ferromagnetic Materials". IEEE Transactions on Magnetics 40 (6): 3443–3449.  
  • Kruzik Martin, Prohl Andreas (2006). "Recent Developments in the Modeling, Analysis, and Numerics of Ferromagnetism". SIAM Review 48 (3): 439–483.  
  • Maugin, Gérard A. (1988). Continuum mechanics of electromagnetic solids. Amsterdam: North-Holland.  
  • Pohl, Andreas (2001). Computational micromagnetism (1. Aufl. ed.). Stuttgart: Teubner.  
  • Tiersten, H. F. (1964). "Coupled Magnetomechanical Equations for Magnetically Saturated Insulators". Journal of Mathematical Physics 5 (9): 1298.  

External links

  • µMAG -- Micromagnetic Modeling Activity Group.
  • OOMMF -- Micromagnetic Modeling Tool.
  • MuMax -- GPU-accelerated Micromagnetic Modeling Tool.
This article was sourced from Creative Commons Attribution-ShareAlike License; additional terms may apply. World Heritage Encyclopedia content is assembled from numerous content providers, Open Access Publishing, and in compliance with The Fair Access to Science and Technology Research Act (FASTR), Wikimedia Foundation, Inc., Public Library of Science, The Encyclopedia of Life, Open Book Publishers (OBP), PubMed, U.S. National Library of Medicine, National Center for Biotechnology Information, U.S. National Library of Medicine, National Institutes of Health (NIH), U.S. Department of Health & Human Services, and, which sources content from all federal, state, local, tribal, and territorial government publication portals (.gov, .mil, .edu). Funding for and content contributors is made possible from the U.S. Congress, E-Government Act of 2002.
Crowd sourced content that is contributed to World Heritage Encyclopedia is peer reviewed and edited by our editorial staff to ensure quality scholarly research articles.
By using this site, you agree to the Terms of Use and Privacy Policy. World Heritage Encyclopedia™ is a registered trademark of the World Public Library Association, a non-profit organization.

Copyright © World Library Foundation. All rights reserved. eBooks from Project Gutenberg are sponsored by the World Library Foundation,
a 501c(4) Member's Support Non-Profit Organization, and is NOT affiliated with any governmental agency or department.