Last updated

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.


Micromagnetics as a field (i.e., that deals specifically with the behaviour of ferromagnetic 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 vectorm = M/Ms, also called reduced magnetization.

The static equilibria are found by minimizing the magnetic energy,


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:

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:

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


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:

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. Micromagnetics by Zureks.svg
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:

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

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

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, one can assume this strain to be isochoric and fully isotropic in the lateral direction, yielding the deviatoric ansatz

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 . A quadratic form for the magnetoelastic energy is

where 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

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 terms 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:

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:

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:

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. Damped Magnetization Precession.jpg
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:

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:


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.

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, [5] vestibular neurons, and cortical neurons of embryonic rats. [6]

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. Bibcode:2012CMAME.245..331M. doi:10.1016/j.cma.2012.03.021.
  2. Komineas, Stavros; Papanicolaou, Nikos (2007). "Dynamics of vortex-antivortex pairs in ferromagnets". arXiv: 0712.3684v1 [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" (PDF). Physical Review B. 67 (9): 094410. Bibcode:2003PhRvB..67i4410T. doi:10.1103/PhysRevB.67.094410. hdl: 10261/25225 .
  4. 1 2 Döring, W. (1968). "Point Singularities in Micromagnetism". Journal of Applied Physics. 39 (2): 1006–1007. Bibcode:1968JAP....39.1006D. doi:10.1063/1.1656144.
  5. Mukesh, S. (2017). "Modeling intracochlear magnetic stimulation: a Finite-Element Analysis". IEEE Transactions on Neural Systems and Rehabilitation Engineering. 25 (8): 1353–1362. doi: 10.1109/TNSRE.2016.2624275 .
  6. Mukesh, S. (2019). Magnetic Stimulation of Dissociated Cortical Neurons on a Planar Mulitelectrode Array. 2019 9th International IEEE/EMBS Conference on Neural Engineering (NER). pp. 1758–761. doi:10.1109/NER.2019.8717125.

Further reading

Related Research Articles

In physics, the cross section is a measure of the probability that a specific process will take place when some kind of radiant excitation intersects a localized phenomenon. For example, the Rutherford cross-section is a measure of probability that an alpha particle will be deflected by a given angle during a collision with an atomic nucleus. Cross section is typically denoted σ (sigma) and is expressed in units of area, more specifically in barns. In a way, it can be thought of as the size of the object that the excitation must hit in order for the process to occur, but more exactly, it is a parameter of a stochastic process.

Maxwells equations Equations describing classical electromagnetism

Maxwell's equations are a set of coupled partial differential equations that, together with the Lorentz force law, form the foundation of classical electromagnetism, classical optics, and electric circuits. The equations provide a mathematical model for electric, optical, and radio technologies, such as power generation, electric motors, wireless communication, lenses, radar etc. They describe how electric and magnetic fields are generated by charges, currents, and changes of the fields. The equations are named after the physicist and mathematician James Clerk Maxwell, who, in 1861 and 1862, published an early form of the equations that included the Lorentz force law. Maxwell first used the equations to propose that light is an electromagnetic phenomenon. The modern form of the equations in their most common formulation is credited to Oliver Heaviside.

Paramagnetism Weak, attractive magnetism possessed by most elements and some compounds

Paramagnetism is a form of magnetism whereby some materials are weakly attracted by an externally applied magnetic field, and form internal, induced magnetic fields in the direction of the applied magnetic field. In contrast with this behavior, diamagnetic materials are repelled by magnetic fields and form induced magnetic fields in the direction opposite to that of the applied magnetic field. Paramagnetic materials include most chemical elements and some compounds; they have a relative magnetic permeability slightly greater than 1 and hence are attracted to magnetic fields. The magnetic moment induced by the applied field is linear in the field strength and rather weak. It typically requires a sensitive analytical balance to detect the effect and modern measurements on paramagnetic materials are often conducted with a SQUID magnetometer.

Ampères circuital law

In classical electromagnetism, Ampère's circuital law relates the integrated magnetic field around a closed loop to the electric current passing through the loop. James Clerk Maxwell derived it using hydrodynamics in his 1861 published paper "On Physical Lines of Force" In 1865 he generalized the equation to apply to time-varying currents by adding the displacement current term, resulting in the modern form of the law, sometimes called the Ampère–Maxwell law, which is one of Maxwell's equations which form the basis of classical electromagnetism.

Magnetic moment Physical quantity; measured in ampere square metre

The magnetic moment is the magnetic strength and orientation of a magnet or other object that produces a magnetic field. Examples of objects that have magnetic moments include: loops of electric current, permanent magnets, elementary particles, various molecules, and many astronomical objects.

Displacement current Physical quantity in electromagnetism

In electromagnetism, displacement current density is the quantity D/∂t appearing in Maxwell's equations that is defined in terms of the rate of change of D, the electric displacement field. Displacement current density has the same units as electric current density, and it is a source of the magnetic field just as actual current is. However it is not an electric current of moving charges, but a time-varying electric field. In physical materials, there is also a contribution from the slight motion of charges bound in atoms, called dielectric polarization.

In physics and engineering, a constitutive equation or constitutive relation is a relation between two physical quantities that is specific to a material or substance, and approximates the response of that material to external stimuli, usually as applied fields or forces. They are combined with other equations governing physical laws to solve physical problems; for example in fluid mechanics the flow of a fluid in a pipe, in solid state physics the response of a crystal to an electric field, or in structural analysis, the connection between applied stresses or forces to strains or deformations.

Magnetostatics Branch of physics concerned with magnetic behavior in systems with steady electric currents

Magnetostatics is the study of magnetic fields in systems where the currents are steady. It is the magnetic analogue of electrostatics, where the charges are stationary. The magnetization need not be static; the equations of magnetostatics can be used to predict fast magnetic switching events that occur on time scales of nanoseconds or less. Magnetostatics is even a good approximation when the currents are not static — as long as the currents do not alternate rapidly. Magnetostatics is widely used in applications of micromagnetics such as models of magnetic storage devices as in computer memory.

Magnetization Physical quantity, density of magnetic moment per volume

In classical electromagnetism, magnetization is the vector field that expresses the density of permanent or induced magnetic dipole moments in a magnetic material. Movement within this field is described by direction and is either Axial or Diametric. The origin of the magnetic moments responsible for magnetization can be either microscopic electric currents resulting from the motion of electrons in atoms, or the spin of the electrons or the nuclei. Net magnetization results from the response of a material to an external magnetic field. Paramagnetic materials have a weak induced magnetization in a magnetic field, which disappears when the magnetic field is removed. Ferromagnetic and ferrimagnetic materials have strong magnetization in a magnetic field, and can be magnetized to have magnetization in the absence of an external field, becoming a permanent magnet. Magnetization is not necessarily uniform within a material, but may vary between different points. Magnetization also describes how a material responds to an applied magnetic field as well as the way the material changes the magnetic field, and can be used to calculate the forces that result from those interactions. It can be compared to electric polarization, which is the measure of the corresponding response of a material to an electric field in electrostatics. Physicists and engineers usually define magnetization as the quantity of magnetic moment per unit volume. It is represented by a pseudovector M.

The electromagnetic wave equation is a second-order partial differential equation that describes the propagation of electromagnetic waves through a medium or in a vacuum. It is a three-dimensional form of the wave equation. The homogeneous form of the equation, written in terms of either the electric field E or the magnetic field B, takes the form:

Maxwell stress tensor

The Maxwell stress tensor is a symmetric second-order tensor used in classical electromagnetism to represent the interaction between electromagnetic forces and mechanical momentum. In simple situations, such as a point charge moving freely in a homogeneous magnetic field, it is easy to calculate the forces on the charge from the Lorentz force law. When the situation becomes more complicated, this ordinary procedure can become impractically difficult, with equations spanning multiple lines. It is therefore convenient to collect many of these terms in the Maxwell stress tensor, and to use tensor arithmetic to find the answer to the problem at hand.

Covariant formulation of classical electromagnetism

The covariant formulation of classical electromagnetism refers to ways of writing the laws of classical electromagnetism in a form that is manifestly invariant under Lorentz transformations, in the formalism of special relativity using rectilinear inertial coordinate systems. These expressions both make it simple to prove that the laws of classical electromagnetism take the same form in any inertial coordinate system, and also provide a way to translate the fields and forces from one frame to another. However, this is not as general as Maxwell's equations in curved spacetime or non-rectilinear coordinate systems.

Inhomogeneous electromagnetic wave equation Wave equations describing the propagation of electromagnetic waves generated by nonzero source charges and currents

In electromagnetism and applications, an inhomogeneous electromagnetic wave equation, or nonhomogeneous electromagnetic wave equation, is one of a set of wave equations describing the propagation of electromagnetic waves generated by nonzero source charges and currents. The source terms in the wave equations make the partial differential equations inhomogeneous, if the source terms are zero the equations reduce to the homogeneous electromagnetic wave equations. The equations follow from Maxwell's equations.

Dielectric loss quantifies a dielectric material's inherent dissipation of electromagnetic energy. It can be parameterized in terms of either the loss angleδ or the corresponding loss tangent tan δ. Both refer to the phasor in the complex plane whose real and imaginary parts are the resistive (lossy) component of an electromagnetic field and its reactive (lossless) counterpart.

Mathematical descriptions of the electromagnetic field Formulations of electromagnetism

There are various mathematical descriptions of the electromagnetic field that are used in the study of electromagnetism, one of the four fundamental interactions of nature. In this article, several approaches are discussed, although the equations are in terms of electric and magnetic fields, potentials, and charges with currents, generally speaking.

Effective medium approximations (EMA) or effective medium theory (EMT) pertain to analytical or theoretical modeling that describes the macroscopic properties of composite materials. EMAs or EMTs are developed from averaging the multiple values of the constituents that directly make up the composite material. At the constituent level, the values of the materials vary and are inhomogeneous. Precise calculation of the many constituent values is nearly impossible. However, theories have been developed that can produce acceptable approximations which in turn describe useful parameters including the effective permittivity and permeability of the materials as a whole. In this sense, effective medium approximations are descriptions of a medium based on the properties and the relative fractions of its components and are derived from calculations, and effective medium theory. There are two widely used formulae.

Retarded potential Type of potential in electrodynamics

In electrodynamics, the retarded potentials are the electromagnetic potentials for the electromagnetic field generated by time-varying electric current or charge distributions in the past. The fields propagate at the speed of light c, so the delay of the fields connecting cause and effect at earlier and later times is an important factor: the signal takes a finite time to propagate from a point in the charge or current distribution to another point in space, see figure below.

When an electromagnetic wave travels through a medium in which it gets attenuated, it undergoes exponential decay as described by the Beer–Lambert law. However, there are many possible ways to characterize the wave and how quickly it is attenuated. This article describes the mathematical relationships among:

In physics, magnetization dynamics is the branch of solid-state physics that describes the evolution of the magnetization of a material.

The Wannier equation describes a quantum mechanical eigenvalue problem in solids where an electron in a conduction band and an electronic vacancy within a valence band attract each other via the Coulomb interaction. For one electron and one hole, this problem is analogous to the Schrödinger equation of the hydrogen atom; and the bound-state solutions are called excitons. When an exciton's radius extends over several unit cells, it is referred to as a Wannier exciton in contrast to Frenkel excitons whose size is comparable with the unit cell. An excited solid typically contains many electrons and holes; this modifies the Wannier equation considerably. The resulting generalized Wannier equation can be determined from the homogeneous part of the semiconductor Bloch equations or the semiconductor luminescence equations.