# Demagnetizing field

The demagnetizing field, also called the stray field, is the magnetic field (H-field)[1] generated by the magnetization in a magnet. The total magnetic field in a region containing magnets is the sum of the demagnetizing fields of the magnets and the magnetic field due to any free currents or displacement currents. The term demagnetizing field reflects its tendency to act on the magnetization so as to reduce the total magnetic moment. It gives rise to shape anisotropy in ferromagnets with a single magnetic domain and to magnetic domains in larger ferromagnets.

The demagnetizing field of an arbitrary shaped object is very difficult to calculate even for the simple case of uniform magnetization. For the special case of ellipsoids (which includes spheres) the demagnetization field is linearly related to the magnetization by a geometry dependent constant called the demagnetizing factor. Since the magnetization of a sample at a given location depends on the total magnetic field at that point, the demagnetization factor must be used in order to accurately determine how a magnetic material responds to a magnetic field. (See magnetic hysteresis.)

## Magnetostatic principles

### Maxwell's equations

Main article: Maxwell's equations

In general the demagnetizing field is a function of position H(r). It is derived from the magnetostatic equations for a body with no electric currents.[2] These are Ampère's law

$\nabla\times\mathbf{H} = 0,$[3]

(1)

and Gauss's law

$\nabla\cdot\mathbf{B} = 0.$[4]

(2)

The magnetic field and flux density are related by[5][6]

$\mathbf{B} = \mu_0 \left(\mathbf{M}+\mathbf{H}\right),$[7]

(3)

where $\mu_0$ is the permeability of vacuum and M is the magnetisation.

### The magnetic potential

Main article: magnetic potential

The general solution of the first equation can be expressed as the gradient of a scalar potential U(r):

$\mathbf{H} = -\nabla U.$[5][6]

(4)

Inside the magnetic body, the potential Uin is determined by substituting (3) and (4) in (2):

$\nabla^2 U_\text{in} = \nabla\cdot\mathbf{M}.$[8]

(5)

Outside the body, where the magnetization is zero,

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

(6)

At the surface of the magnet, there are two continuity requirements:[5]

• The component of H parallel to the surface must be continuous (no jump in value at the surface).
• The component of B perpendicular to the surface must be continuous.

This leads to the following boundary conditions at the surface of the magnet:

\begin{align} U_\text{in} &= U_\text{out}\\ \frac{\partial U_\text{in}}{\partial n} &= \frac{\partial U_\text{out}}{\partial n} + \mathbf{M}\cdot\mathbf{n}. \end{align}

(7)

Here n is the surface normal and $\scriptstyle \left(\partial/ \partial n\right)$ is the derivative with respect to distance from the surface.[9]

The outer potential Uout must also be regular at infinity: both |r U| and |r2 U| must be bounded as r goes to infinity. This ensures that the magnetic energy is finite.[10] Sufficiently far away, the magnetic field looks like the field of a magnetic dipole with the same moment as the finite body.

### Uniqueness of the demagnetizing field

Any two potentials that satisfy equations (5), (6) and (7), along with regularity at infinity, are identical. The demagnetizing field Hd is the gradient of this potential (equation 4).

### Energy

The energy of the demagnetizing field is completely determined by an integral over the volume V of the magnet:

$E = -\frac{1}{2}\int_\text{magnet} \mathbf{M}\cdot\mathbf{H}_\text{d} dV$

(7)

Suppose there are two magnets with magnetizations M1 and M2. The energy of the first magnet in the demagnetizing field Hd(2) of the second is

$E = \int_\text{magnet 1} \mathbf{M}_1\cdot\mathbf{H}_\text{d}^{(2)} dV.$

(8)

The reciprocity theorem states that[9]

$\int_\text{magnet 1} \mathbf{M}_1\cdot\mathbf{H}_\text{d}^{(2)} dV = \int_\text{magnet 2} \mathbf{M}_2\cdot\mathbf{H}_\text{d}^{(1)} dV.$

(9)

### Magnetic charge and the pole-avoidance principle

Formally, the solution of the equations for the potential is

$U(\mathbf{r}) = -\frac{1}{4\pi}\int_\text{volume} \frac{\nabla'\cdot\mathbf{M\left(r'\right)}}{|\mathbf{r}-\mathbf{r}'|}dV' + \frac{1}{4\pi}\int_\text{surface} \frac{\mathbf{n}\cdot\mathbf{M\left(r'\right)}}{|\mathbf{r}-\mathbf{r}'|}dS',$

(10)

where r is the variable to be integrated over the volume of the body in the first integral and the surface in the second, and is the gradient with respect to this variable.[9]

Qualitatively, the negative of the divergence of the magnetization − ∇ · M (called a volume pole) is analogous to a bulk bound electric charge in the body while n · M (called a surface pole) is analogous to a bound surface electric charge. Although the magnetic charges do not exist, it can be useful to think of them in this way. In particular, the arrangement of magnetization that reduces the magnetic energy can often be understood in terms of the pole-avoidance principle, which states that the magnetization tries to reduce the poles as much as possible.[9]

## Effect on magnetization

### Single domain

Illustration of the magnetic charges at the surface of a single-domain ferromagnet. The arrows indicate the direction of magnetization. The thickness of the colored region indicates the surface charge density.

One way to remove the magnetic poles inside a ferromagnet is to make the magnetization uniform. This occurs in single-domain ferromagnets. This still leaves the surface poles, so division into domains reduces the poles further. However, very small ferromagnets are kept uniformly magnetized by the exchange interaction.

The concentration of poles depends on the direction of magnetization (see the figure). If the magnetization is along the longest axis, the poles are spread across a smaller surface, so the energy is lower. This is a form of magnetic anisotropy called shape anisotropy.

### Multiple domains

Illustration of a magnet with four magnetic closure domains. The magnetic charges contributed by each domain are pictured at one domain wall. The charges balance, so the total charge is zero.

If the ferromagnet is large enough, its magnetization can divide into domains. It is then possible to have the magnetization parallel to the surface. Within each domain the magnetization is uniform, so there are no volume poles, but there are surface poles at the interfaces (domain walls) between domains. However, these poles vanish if the magnetic moments on each side of the domain wall meet the wall at the same angle (so that the components n · M are the same but opposite in sign). Domains configured this way are called closure domains.

## Demagnetizing factor

An arbitrary shaped magnetic object has a total magnetic field that varies with location inside the object and can be quite difficult to calculate. This makes it very difficult to determine the magnetic properties of a material such as, for instance, how the magnetization of a material varies with the magnetic field. For a uniformly magnetized sphere in a uniform magnetic field H0 the internal magnetic field H is uniform:

$\mathbf{H}=\mathbf{H}_0 - \frac{\gamma}{4\pi} \mathbf{M}_0,$

(11)

where M0 is the magnetization of the sphere and γ is called the demagnetizing factor and equals 4π/3 for a sphere.[5][6][11]

This equation can be generalized to include ellipsoids having principal axes in x, y, and z directions such that each component has a relationship of the form:[6]

$H_k=(H_0)_k - \frac{\gamma_k}{4\pi} (M_0)_k, \qquad k = x,y,z.$

(12)

Other important examples are an infinite plate (an ellipsoid with two of its axes going to infinity) which has γ = 4π in a direction normal to the plate and zero otherwise and an infinite cylinder (an ellipsoid with one of its axes tending toward infinity with the other two being the same) which has γ = 2π along its axis.[12] The demagnetizing factors are the principal values of the depolarization tensor, which gives both the internal and external values of the fields induced in ellipsoidal bodies by applied electric or magnetic fields.[13][14]

## Notes and references

1. ^ In this article the term 'magnetic field' is used for the magnetic 'H field' while 'magnetic flux density' is used for the magnetic 'B-field'.
2. ^ If there are electric currents in the system, they can be calculated separately and added to the solutions of these equations.
3. ^ In words, the curl of the magnetic field is zero.
4. ^ In words, the divergence of the magnetic flux density is zero.
5. ^ a b c d Jackson 1975, chapter 5
6. ^ a b c d Nayfeh & Brussel 1985, chapter 9
8. ^ The symbol 2 · is the Laplace operator.
9. ^ a b c d Aharoni 1996, chapter 6
10. ^ Brown, Jr. 1962
11. ^ Griffiths 1999, chapter 6
12. ^ For tables or equations for the magnetizing factors of the general ellipsoid see Osborn, J. A. (1945). "Demagnetizing Factors of the General Ellipsoid". Physical Review 67 (11–12): 351–7. Bibcode:1945PhRv...67..351O. doi:10.1103/PhysRev.67.351.
13. ^ Solivérez, C. E. (1981). "Magnetostatics of Anisotropic Ellipsoidal Bodies". IEEE Transactions on Magnetics 17 (3): 1363–4. Bibcode:1981ITM....17.1363S. doi:10.1109/TMAG.1981.1061200.
14. ^ Solivérez, C. E. (2012). Electrostática y Magnetostática de Cuerpos Elipsoidales: Formalismo del Tensor Depolarización (Electrostatics and Magnetostatics of Ellipsoidal Bodies: The Depolarization Tensor Method). Información Libre. ISBN 978-987-28304-0-3.