Micromagnetics

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

Micromagnetics originated from a 1935 paper

by Lev Landau and Evgeny Lifshitz on antidomain walls.{{rp|133}}{{cite journal | last=Aharoni | first=Amikam | title=Micromagnetics: past, present and future | journal=Physica B: Condensed Matter | volume=306 | issue=1–4 | date=2001 | doi=10.1016/S0921-4526(01)00954-1 | pages=1–9| bibcode=2001PhyB..306....1A }}{{cite journal | last=Brown | first=William Fuller | title=Domains, micromagnetics, and beyond: Reminiscences and assessments | journal=Journal of Applied Physics | volume=49 | issue=3 | date=1978-03-01 | issn=0021-8979 | doi=10.1063/1.324811 | pages=1937–1942| bibcode=1978JAP....49.1937B }}{{cite journal | last1=Bar'yakhtar | first1=V. G. | last2=Ivanov | first2=B. A. | title=The Landau-Lifshitz equation: 80 years of history, advances, and prospects | journal=Low Temperature Physics | volume=41 | issue=9 | date=2015-09-01 | issn=1063-777X | doi=10.1063/1.4931649 | pages=663–669| bibcode=2015LTP....41..663B }}{{cite journal| last1=Landau | first1=L | last2=Lifshitz | first2=E| journal=Sov. Phys. | volume=8 | date=1935 | title=On the theory of magnetic permeability dispersion in ferromagnetic solids | pages=153–166}}{{rp|440}}

Micromagnetics was then expanded upon by William Fuller Brown Jr. in several works in 1940-1941{{rp|133}}{{primary source inline| reason= Brown is essentially a primary source for this material, which makes it difficult to take anything without performing original research. |date=June 2024}}{{cite journal | last=Brown | first=William Fuller | title=Theory of the Approach to Magnetic Saturation | journal=Physical Review | volume=58 | issue=8 | date=1940-10-15 | issn=0031-899X | doi=10.1103/PhysRev.58.736 | pages=736–743| bibcode=1940PhRv...58..736B }}{{cite journal | last=Brown | first=William Fuller | title=The Effect of Dislocations on Magnetization Near Saturation | journal=Physical Review | volume=60 | issue=2 | date=1941-07-15 | issn=0031-899X | doi=10.1103/PhysRev.60.139 | pages=139–147| bibcode=1941PhRv...60..139B }}{{rp|440}} using energy expressions taken from a 1938 paper by William Cronk Elmore.{{cite journal | last=Elmore | first=W. C. | title=The Magnetic Structure of Cobalt | journal=Physical Review | volume=53 | issue=9 | date=1938-05-01 | issn=0031-899X | doi=10.1103/PhysRev.53.757 | pages=757–764}}

According to D. Wei, Brown introduced the name "micromagnetics" in 1958.{{cite book | last=Wei | first=Dan | title=Micromagnetics and Recording Materials | publisher=Springer Science & Business Media | publication-place=Heidelberg; New York | date=2012-04-28 | isbn=978-3-642-28577-6 | oclc=779865001 | page=}}{{rp|41}}{{cite journal | last=Brown | first=William Fuller | title=Rigorous Approach to the Theory of Ferromagnetic Microstructure | journal=Journal of Applied Physics | volume=29 | issue=3 | date=1958-03-01 | issn=0021-8979 | doi=10.1063/1.1723183 | pages=470–471| bibcode=1958JAP....29..470B }}

The field prior to 1960 was summarised in Brown's book Micromagnetics.{{rp|41}}

In the 1970s computational methods were developed for the analysis of recording media due to the introduction of personal computers.{{rp|44}}

Static micromagnetics

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

The static equilibria are found by minimizing the magnetic energy,{{rp|138}}

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

subject to the constraint $|\mathbf{M}|=M_{s}$ or $|\mathbf{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:{{cite book | last=Aharoni | first=Amikam | title=Introduction to the Theory of Ferromagnetism | publisher=Oxford University Press | publication-place=Oxford; New York | date=2007 |orig-date=2000 | isbn=978-0-19-850809-0 | page=135 |edition=second}}{{cite book | last1=Hubert | first1=A. | last2=Schäfer | first2=R. | title=Magnetic Domains: The Analysis of Magnetic Microstructures | publisher=Springer | year=1998 | isbn=978-3-540-64108-7 | url=https://books.google.com/books?id=1BLBWbPXrUgC&pg=PR19 | access-date=2024-06-23}}{{rp|101–104}}

: $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; $m_{x}$ , $m_{y}$ and $m_{z}$ are the components of $\mathbf{m}$ ;

and the integral is performed over the volume of the sample.

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

The exchange term is isotropic,

so any direction is equally acceptable.{{rp|83}}

= Anisotropy energy =

{{Main|Magnetic anisotropy|Anisotropy energy|Magnetocrystalline anisotropy}}

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

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

where $F_{\text{anis}}$ , the anisotropy energy density, is a function of the orientation of the magnetization. Minimum-energy directions for $F_{\text{anis}}$ are called easy axes.

Time-reversal symmetry ensures that $F_{\text{anis}}$ is an even function of $\mathbf{m}$ .{{rp|108}} The simplest such function is

: $F_\text{anis}(\mathbf{m}) = -K_1 m_z^2,$

where K₁ is called the anisotropy constant. In this approximation, called uniaxial anisotropy, the easy axis is the $z$ axis.{{rp|85}}

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 is written as:{{rp|174}}{{rp|109}}

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

where $\mathbf{H}_{\text{a}}$ is the applied field and $\mu_0$ is the vacuum permeability.

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

= Energy of the demagnetizing field =

File:Micromagnetics by Zureks.svg

The demagnetizing field is the magnetic field created by the magnetic sample upon itself. The associated energy is:{{rp|110}}

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

where $\mathbf{H}_{\text{d}}$ is the demagnetizing field. The field satisfies

: $\nabla\times\mathbf{H}_\text{d} = 0$

and hence can be written as the gradient of a potential $\mathbf{H}_\text{d} = -\nabla U$ . This field depends on the magnetic configuration itself, and it can be found by solving

: $\nabla^2 U_{\text{in}} = \nabla\cdot\mathbf{M}$

inside of the body and

: $\nabla^2 U_{\text{out}} =0$

outside of the body.

These are supplemented with the boundary conditions on the surface of the body

: $U_{\text{out}} =U_{\text{in}},\quad \frac{\partial U_{\text{in}}}{\partial \mathbf{n}} - \frac{\partial U_{\text{out}}}{\partial \mathbf{n}}=\mathbf{M}\cdot\mathbf{n}$

where $\mathbf{n}$ is the unit normal to the surface. Furthermore, the potential satisfies the condition that $|rU|$ and $|r^2\nabla U|$ remain bounded as $r\to\infty$ .{{rp|109–111}} The solution of these equations (cf. magnetostatics) is:

: $U(\mathbf{r}) = \frac{1}{4\pi} \left(-\int_V \frac{\nabla'\cdot\mathbf{M}(\mathbf{r}')}$

\mathbf{r}-\mathbf{r}'

\mathrm{d}V + \int_{\partial V}\frac{\mathbf{n}\cdot\mathbf{M}(\mathbf{r}')}

\mathbf{r}-\mathbf{r}'

\mathrm{d}S\right).

The quantity $-\nabla\cdot\mathbf{M}$ is often called the volume charge density, and $\mathbf{M}\cdot\mathbf{n}$ is called the surface charge density.{{rp|125–126}}{{rp|110}}{{cite journal | last=Kruzík | first=Martin | last2=Prohl | first2=Andreas | title=Recent Developments in the Modeling, Analysis, and Numerics of Ferromagnetism | journal=SIAM Review | volume=48 | issue=3 | date=2006 | issn=0036-1445 | doi=10.1137/S0036144504446187 | pages=439–483}}{{rp|441}}

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.{{citation needed|date=June 2024}} The static equilibrium then is a compromise that minimizes the total magnetic energy, although it may not minimize individually any particular term.

= Dzyaloshinskii–Moriya Interaction Energy =

This interaction arises when a crystal lacks inversion symmetry, encouraging the magnetization to be perpendicular to its neighbours. It directly competes with the exchange energy. It is modelled with the energy contribution{{cite journal | last1=Hoffmann | first1=Markus | last2=Zimmermann | first2=Bernd | last3=Müller | first3=Gideon P. | last4=Schürhoff | first4=Daniel | last5=Kiselev | first5=Nikolai S. | last6=Melcher | first6=Christof | last7=Blügel | first7=Stefan | title=Antiskyrmions stabilized at interfaces by anisotropic Dzyaloshinskii-Moriya interactions | journal=Nature Communications | volume=8 | issue=1 | date=2017-08-21 | issn=2041-1723 | pmid=28827700 | pmc=5566362 | doi=10.1038/s41467-017-00313-0 | page=308| arxiv=1702.07573 | bibcode=2017NatCo...8..308H }}

$E_\text{DMI} = \int_{V}\mathbf{D}:(\nabla \mathbf{m}\times \mathbf{m})$

where $\mathbf{D}$ is the spiralization tensor,

that depends upon the crystal class.{{cite journal | last1=Cortés-Ortuño | first1=David | last2=Beg | first2=Marijan | last3=Nehruji | first3=Vanessa | last4=Breth | first4=Leoni | last5=Pepper | first5=Ryan | last6=Kluyver | first6=Thomas | last7=Downing | first7=Gary | last8=Hesjedal | first8=Thorsten | last9=Hatton | first9=Peter | last10=Lancaster | first10=Tom | last11=Hertel | first11=Riccardo | last12=Hovorka | first12=Ondrej | last13=Fangohr | first13=Hans | title=Proposal for a micromagnetic standard problem for materials with Dzyaloshinskii–Moriya interaction | journal=New Journal of Physics | volume=20 | issue=11 | date=2018-11-12 | issn=1367-2630 | doi=10.1088/1367-2630/aaea1c | page=113015| arxiv=1803.11174 | bibcode=2018NJPh...20k3015C }} For bulk DMI,

: $E_\text{DMI} = \int_{V}D \mathbf{m}\cdot(\nabla \times \mathbf{m}),$

and for a thin film in the $x-y$ plane

interfacial DMI takes the form

: $E_\text{DMI} = \int_{V}D(\mathbf{m}\cdot\nabla m_{z} - m_{z}\nabla\cdot\mathbf{m}),$

and for materials with symmetry class $D_{2d}$ the energy contribution is

: $E_\text{DMI} = \int_{V}D \mathbf{m}\cdot\left(\frac{\partial \mathbf{m}}{\partial x}\times \hat{x} - \frac{\partial \mathbf{m}}{\partial y}\times \hat{y}\right).$

This term is important for the formation of magnetic skyrmions.

= 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 $\mathbf{m}$ .

For a simple small-strain model, one can assume this strain to be isochoric and fully

isotropic in the lateral direction, yielding the deviatoric ansatz{{rp|128}} {{cite book | last1=Miyazaki | first1=Terunobu | last2=Jin | first2=Hanmin | title=The Physics of Ferromagnetism | publisher=Springer Science & Business Media | publication-place=Berlin New York | date=2012-08-22 | isbn=978-3-642-25583-0}}{{rp|250–251}}

$\mathbf{\varepsilon}_0(\mathbf{m}) = \frac{3}{2} \lambda_{\text{s}}\, \left[\mathbf{m}\otimes \mathbf{m} - \frac{1}{3}\mathbf{1}\right]$

where the material parameter $\lambda_{\text{s}}$ is the isotropic magnetostrictive

constant. The elastic

energy density is assumed 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{{rp|138}}

$E_\text{m-e} = \frac{1}{2} \int_{V}[\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 $\lambda$ and $\mu$ ).

Taking into account the constant length of $\mathbf{m}$ , we obtain the invariant-based representation

E_\text{m-e} = \int_{V}

\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.{{rp|181–182}} This is especially important if the sample is 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 terms of the local effective field acting on it.

= Effective field =

The effective field is the local field felt by the magnetization. The only real fields however are the magnetostatic field and the applied field. 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$

Since m is a unit vector, dm is always perpendicular to m. Then the above definition leaves unspecified the component of H_eff that is parallel to m.{{cite book | last=Brown Jr. | first=W.F. | title=Micromagnetics | publisher=Interscience Publishers | series=Interscience tracts on physics and astronomy | year=1963 | volume=18 | page=7 | url=https://books.google.com/books?id=KvuXnAEACAAJ | access-date=2024-06-14 }} 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 (excluding the DMI and magnetoelastic contributions):{{rp|178}}

: $\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 =

File:Damped Magnetization Precession.jpg

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:{{rp|181}}{{rp|462}}

: $\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 $\gamma$ is the electron gyromagnetic ratio and $\alpha$ the Gilbert damping constant.

It can be shown that this is mathematically equivalent to the following Landau-Lifshitz (or explicit) form:{{Cite journal |last=Abert |first=Claas |date=June 2019 |title=Micromagnetics and spintronics: models and numerical methods |journal=The European Physical Journal B |language=en |volume=92 |issue=6 |pages=120 |doi=10.1140/epjb/e2019-90599-6 |bibcode=2019EPJB...92..120A |issn=1434-6028|doi-access=free |arxiv=1810.12365 }}{{rp|181–182}}

: $\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}),

where $\alpha$ is the Gilbert Damping constant, characterizing how quickly the damping term takes away energy from the system ( $\alpha$ = 0, no damping, permanent precession).

These equations preserve the constraint $|\mathbf{m}| = 1$ , as{{rp|181}}

: $\frac{\mathrm{d}}{\mathrm{d}t}|\mathbf{m}|^2 = 2\mathbf{m}\cdot\frac{\partial \mathbf{m}}{\partial t}=0.$

Applications

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.

{{cite journal

| last1 = Miehe

| first1 = Christian

| last2 = Ethiraj

| first2 = Gautam

| date = 2011-10-15

| title = A geometrically consistent incremental variational formulation for phase field models in micromagnetics

| journal = Computer Methods in Applied Mechanics and Engineering

| publisher = Elsevier

| volume = 245–246

| issue =

| pages = 331–347

| doi = 10.1016/j.cma.2012.03.021

|bibcode = 2012CMAME.245..331M }}

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;{{cite arXiv |last=Komineas |first=Stavros|author2=Papanicolaou, Nikos |authorlink= |eprint=0712.3684v1 |title=Dynamics of vortex-antivortex pairs in ferromagnets |class=cond-mat.mtrl-sci |year=2007}} or even 3d-Bloch points,{{cite journal|last=Thiaville|first=André|author2=García, José |author3=Dittrich, Rok |author4=Miltat, Jacques |author5= Schrefl, Thomas |title=Micromagnetic study of Bloch-point-mediated vortex core reversal|journal=Physical Review B|date=March 2003|volume=67|issue=9|pages=094410|doi=10.1103/PhysRevB.67.094410|bibcode = 2003PhRvB..67i4410T |url=https://digital.csic.es/bitstream/10261/25225/1/Thiaville%2c%20A.%20et%20al%20Phys.Rev.B_67_2003.pdf|hdl=10261/25225|hdl-access=free}}{{cite journal|last=Döring|first=W.|title=Point Singularities in Micromagnetism|journal=Journal of Applied Physics|year=1968|volume=39|issue=2|pages=1006–1007|doi=10.1063/1.1656144|bibcode = 1968JAP....39.1006D }} 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 are thought{{by whom|date=June 2024}} to be used as information carriers, to apply the most recent, and already studied, propositions in information technology.

Another application that has emerged in the last decade is the application of micromagnetics towards neuronal stimulation. In this discipline, numerical methods such as finite-element analysis are used to analyze the electric/magnetic fields generated by the stimulation apparatus; then the results are validated or explored further using in-vivo or in-vitro neuronal stimulation. Several distinct set of neurons have been studied using this methodology including retinal neurons, cochlear neurons,{{cite journal|last=Mukesh|first=S.|title=Modeling intracochlear magnetic stimulation: a Finite-Element Analysis|journal=IEEE Transactions on Neural Systems and Rehabilitation Engineering|year=2017|volume=25|issue=8|pages=1353–1362|doi=10.1109/TNSRE.2016.2624275 |pmid=27831887 |pmc=5750049 |doi-access=free}} vestibular neurons, and cortical neurons of embryonic rats.{{cite conference|last=Mukesh|first=S.|title=Magnetic Stimulation of Dissociated Cortical Neurons on a Planar Mulitelectrode Array|conference=2019 9th International IEEE/EMBS Conference on Neural Engineering (NER)|year=2019|pages=1758–761|doi=10.1109/NER.2019.8717125 }}

Footnotes and references

External links

[http://www.ctcms.nist.gov/mumag/mumag.org.html μMAG -- Micromagnetic Modeling Activity Group].
[http://math.nist.gov/oommf OOMMF -- Micromagnetic Modeling Tool].
[https://mumax.github.io MuMax -- GPU-accelerated Micromagnetic Modeling Tool].

Category:Dynamical systems

Category:Magnetic ordering

Category:Magnetostatics