Skip to content

Magnetism

A bar magnet has two poles: like poles repel and opposite poles attract. Cutting it in half gives two smaller magnets, each with both poles, rather than an isolated north or south pole. No isolated magnetic monopole has been observed.

Different materials respond differently. Magnetite can be naturally magnetized; some materials become magnetized near another magnet; others respond so weakly that the effect is hard to notice.

For calculations, use fields and currents instead of treating the poles as separate magnetic charges. Moving charges produce magnetic fields and experience magnetic forces. Intrinsic magnetic moments, such as electron spin, also contribute to magnetism in matter.

The magnetic field B⃗(x,y,z)\vec B(x,y,z) is a vector field with three components. It is also called magnetic induction or magnetic flux density, and its SI unit is the tesla:

1 T=1 NA⋅m.1\ \text{T}=1\ \frac{\text{N}}{\text{A}\cdot\text{m}}.

Field lines are tangent to B⃗\vec B, with closer spacing representing a stronger field. They do not begin or end on magnetic charges. Around a straight wire they are circles; around a current loop they resemble the field of a bar magnet. They cannot cross where the field has a well-defined nonzero direction. The field exists between the drawn lines too, and a two-dimensional sketch only shows part of a three-dimensional field.


A small current element I dℓ⃗′I\,d\vec\ell' produces a field at an observation point. Define R⃗=r⃗−r⃗′\vec R=\vec r-\vec r' from the source element to that point, with R=∣R⃗∣R=\lvert\vec R\rvert. Then

dB⃗=μ04πI dℓ⃗′×R^R2=μ0I4πdℓ⃗′×R⃗R3.d\vec B=\frac{\mu_0}{4\pi}\frac{I\,d\vec\ell'\times\hat R}{R^2} =\frac{\mu_0 I}{4\pi}\frac{d\vec\ell'\times\vec R}{R^3}.

Here dℓ⃗′d\vec\ell' points along conventional current, and the permeability of free space is approximately

μ0≈4π×10−7 T⋅m/A.\mu_0\approx4\pi\times10^{-7}\ \text{T}\cdot\text{m/A}. A current element and displacement vector to an observation point, showing the perpendicular magnetic field contribution

Magnetic fields obey superposition. Add the contributions from every segment of the source circuit:

B⃗(r⃗)=μ0I4π∫Cdℓ⃗′×R^R2.\vec B(\vec r)=\frac{\mu_0 I}{4\pi}\int_C\frac{d\vec\ell'\times\hat R}{R^2}.

This is the magnetostatic form: currents and charge distributions are steady. Use it directly for currents in vacuum, or when magnetic effects of the surrounding material can be neglected. A single moving point charge is not a steady current distribution; its general field requires the time-dependent theory, not a direct substitution into this wire formula.

Example. An infinitely long wire lies on the xx-axis and carries current II in the positive xx direction. Find the field at P=(0,a,0)P=(0,a,0), where a>0a>0.

For a source element at (x,0,0)(x,0,0),

dℓ⃗′=dx x^,R⃗=−xx^+ay^,dℓ⃗′×R⃗=a dx z^.d\vec\ell'=dx\,\hat x,\qquad \vec R=-x\hat x+a\hat y,\qquad d\vec\ell'\times\vec R=a\,dx\,\hat z.

Every contribution points out of the page. Integrating along the wire,

B⃗=μ0I4π∫−∞∞a dx(x2+a2)3/2 z^.\vec B=\frac{\mu_0 I}{4\pi}\int_{-\infty}^{\infty} \frac{a\,dx}{(x^2+a^2)^{3/2}}\,\hat z.

Use x=atan⁡θx=a\tan\theta. The integrand becomes cos⁡θ dθ/a\cos\theta\,d\theta/a, so

B⃗=μ0I4πa∫−π/2π/2cos⁡θ dθ z^=μ0I2πa z^.\vec B=\frac{\mu_0 I}{4\pi a}\int_{-\pi/2}^{\pi/2}\cos\theta\,d\theta\,\hat z =\frac{\mu_0 I}{2\pi a}\,\hat z.

At perpendicular distance rr from the wire, the magnitude is B=μ0I/(2πr)B=\mu_0 I/(2\pi r), with direction tangent to a circle around the wire.

Example. A circular loop of radius aa lies in the xyxy-plane and carries current II counterclockwise when viewed from positive zz. Find its field on the zz-axis.

Opposite elements give equal and opposite transverse components, leaving only the axial component. Every source element is distance R=a2+z2R=\sqrt{a^2+z^2} from the observation point, and

dBz=μ0I4πa dℓ′(a2+z2)3/2.dB_z=\frac{\mu_0 I}{4\pi}\frac{a\,d\ell'}{(a^2+z^2)^{3/2}}.

Since ∫dℓ′=2πa\int d\ell'=2\pi a,

B⃗(z)=μ0Ia22(a2+z2)3/2 z^.\vec B(z)=\frac{\mu_0 I a^2}{2(a^2+z^2)^{3/2}}\,\hat z.

At the center,

B⃗(0)=μ0I2a z^.\vec B(0)=\frac{\mu_0 I}{2a}\,\hat z.

Reverse the current and the field reverses. For NN identical closely stacked turns, multiply by NN.

For a surface current density K⃗\vec K, measured in A/m\text{A/m}, or a volume current density J⃗\vec J, measured in A/m2\text{A/m}^2, replace the wire element by the appropriate distributed-current element:

B⃗(r⃗)=μ04π∫SK⃗(r⃗′)×R^R2 dA′,\vec B(\vec r)=\frac{\mu_0}{4\pi}\int_S \frac{\vec K(\vec r')\times\hat R}{R^2}\,dA', B⃗(r⃗)=μ04π∫VJ⃗(r⃗′)×R^R2 dV′.\vec B(\vec r)=\frac{\mu_0}{4\pi}\int_V \frac{\vec J(\vec r')\times\hat R}{R^2}\,dV'.

These integrals can be difficult. Check symmetry before committing to a direct integration: Ampère’s law may give the same field much faster.


Magnetic flux measures the field passing through an oriented surface:

ΦB=∫SB⃗⋅dA⃗.\Phi_B=\int_S\vec B\cdot d\vec A.

For a uniform field through a flat surface,

ΦB=BAcos⁡θ,\Phi_B=BA\cos\theta,

where θ\theta is the angle between the field and the surface normal, not the surface itself. The unit is the weber: 1 Wb=1 T⋅m21\ \text{Wb}=1\ \text{T}\cdot\text{m}^2.

Gauss’s law for magnetism states that the net flux through any closed surface is zero:

∮SB⃗⋅dA⃗=0,∇⋅B⃗=0.\oint_S\vec B\cdot d\vec A=0, \qquad \nabla\cdot\vec B=0.

Whatever flux enters a closed surface also leaves it. An open surface can have nonzero flux; the closed-surface condition is essential. This is the field-law statement that there are no magnetic monopoles.

A divergence-free magnetic field can be represented using a vector potential A⃗\vec A:

B⃗=∇×A⃗.\vec B=\nabla\times\vec A.

The identity ∇⋅(∇×A⃗)=0\nabla\cdot(\nabla\times\vec A)=0 makes this consistent with Gauss’s law. Unlike electric potential, A⃗\vec A is a vector, not a scalar. It is not unique: adding ∇f\nabla f leaves its curl unchanged.


For steady currents in vacuum,

∮CB⃗⋅dℓ⃗=μ0Ienc,Ienc=∫SJ⃗⋅dA⃗.\oint_C\vec B\cdot d\vec\ell=\mu_0 I_{\text{enc}}, \qquad I_{\text{enc}}=\int_S\vec J\cdot d\vec A.

Choose a direction around the Amperian loop first. Curl your right-hand fingers in that direction; your thumb gives the positive surface normal. Currents crossing in that direction count positively, and currents crossing the other way count negatively. Use the algebraic sum of enclosed currents.

The law is true for any closed loop in magnetostatics, but it only makes the calculation simple when symmetry tells you enough about B⃗\vec B. Zero enclosed current means zero circulation, not necessarily zero field everywhere on the loop. Time-dependent fields require the additional displacement-current term; the formulas here use the steady-current limit.

Example. An infinite sheet lies in the xyxy-plane and carries uniform surface current K⃗=Kx^\vec K=K\hat x. Find the field on either side.

Symmetry and the right-hand rule give equal field magnitudes: −y^-\hat y above the sheet and +y^+\hat y below it. Choose a rectangular loop in the yzyz-plane with length LL parallel to the field on each side.

The two parallel sides contribute BLBL each. The other two sides contribute zero, and the enclosed current is KLKL. Thus

2BL=μ0KL,B=μ0K2.2BL=\mu_0 KL,\qquad B=\frac{\mu_0 K}{2}.

Therefore,

B⃗(z>0)=−μ0K2y^,B⃗(z<0)=+μ0K2y^.\vec B(z>0)=-\frac{\mu_0 K}{2}\hat y, \qquad \vec B(z<0)=+\frac{\mu_0 K}{2}\hat y.

The field magnitude does not decrease with distance for this ideal infinite sheet.

Example. An infinitely long cylindrical wire of radius RR carries total current II uniformly through its cross-section. Find B(r)B(r) inside and outside the wire.

The current density is J=I/(πR2)J=I/(\pi R^2). Cylindrical symmetry makes the field tangent to circles centered on the axis and constant around a circle of radius rr.

Inside, only the current within radius rr is enclosed:

Ienc=Jπr2=Ir2R2,B(2πr)=μ0Ir2R2.I_{\text{enc}}=J\pi r^2=I\frac{r^2}{R^2}, \qquad B(2\pi r)=\mu_0 I\frac{r^2}{R^2}.

Outside, all the current is enclosed. Combining the two regions,

B(r)={μ0Ir2πR2,0≤r≤R,μ0I2πr,r≥R.B(r)= \begin{cases} \dfrac{\mu_0 I r}{2\pi R^2},&0\le r\le R,\\[6pt] \dfrac{\mu_0 I}{2\pi r},&r\ge R. \end{cases}

The field starts at zero, grows linearly inside, and then falls as 1/r1/r outside. Both expressions give μ0I/(2πR)\mu_0 I/(2\pi R) at the surface. Its direction follows the right-hand rule around the current.

A tightly wound solenoid acts like many circular current loops stacked together. Define the turn density n=N/Ln=N/L, the number of turns per unit length.

Solenoid winding with the right-hand grip rule relating current direction to the axial magnetic field

Proof (Field of an ideal solenoid). Model the winding as an infinitely long cylindrical current sheet. Symmetry gives an axial field. Rectangular Amperian loops with both long sides outside show that the external axial field is constant; requiring the solenoid’s field to vanish far away makes that constant zero. Loops entirely inside similarly give a uniform interior field.

Now choose a rectangular loop with one long side of length ℓ\ell inside and the other outside. The short sides are perpendicular to the field. The surface cuts through nℓn\ell turns, so

Bℓ=μ0(nℓ)I.B\ell=\mu_0(n\ell)I.

Therefore,

Binside=μ0nI,Boutside=0.B_{\text{inside}}=\mu_0 nI,\qquad B_{\text{outside}}=0.

The direction follows your right thumb when your fingers curl along the winding current. A finite solenoid has a nonzero external field and end effects; the ideal result is a good approximation deep inside a long solenoid.


Magnetic force and charged-particle motion

Section titled “Magnetic force and charged-particle motion”

The magnetic part of the Lorentz force is

F⃗B=qv⃗×B⃗,FB=∣q∣vBsin⁡θ.\vec F_B=q\vec v\times\vec B, \qquad F_B=\lvert q\rvert vB\sin\theta.

Point your right index finger along v⃗\vec v and your middle finger along B⃗\vec B. Your thumb gives the force on a positive charge. Reverse it for a negative charge.

The force is perpendicular to both v⃗\vec v and B⃗\vec B, so

PB=F⃗B⋅v⃗=0.P_B=\vec F_B\cdot\vec v=0.

A magnetic field alone changes a particle’s direction, not its speed or kinetic energy. The full Lorentz force is q(E⃗+v⃗×B⃗)q(\vec E+\vec v\times\vec B); an electric field can do work, so the no-work statement applies only to the magnetic part.

If v⃗⊥B⃗\vec v\perp\vec B in a uniform field, the magnetic force supplies the centripetal force:

∣q∣vB=mv2r,r=mv∣q∣B=p∣q∣B.\lvert q\rvert vB=\frac{mv^2}{r},\qquad r=\frac{mv}{\lvert q\rvert B}=\frac{p}{\lvert q\rvert B}.

For nonrelativistic motion, the angular-frequency magnitude and period are

ωc=∣q∣Bm,T=2πm∣q∣B.\omega_c=\frac{\lvert q\rvert B}{m}, \qquad T=\frac{2\pi m}{\lvert q\rvert B}.

The period does not depend on speed: faster particles trace proportionally larger circles. This speed independence is nonrelativistic; at relativistic speeds the period includes a factor of γ\gamma. The sign of qq determines the direction of rotation, not the sign of the radius or period.

Resolve the initial velocity into components parallel and perpendicular to the uniform field:

v⃗=v⃗∥+v⃗⊥.\vec v=\vec v_{\parallel}+\vec v_{\perp}.

The parallel component is unchanged because its cross product with B⃗\vec B vanishes. The perpendicular component undergoes uniform circular motion. Combining them gives a helix with

r=mv⊥∣q∣B,T=2πm∣q∣B,pitch=∣v∥∣T.r=\frac{mv_{\perp}}{\lvert q\rvert B}, \qquad T=\frac{2\pi m}{\lvert q\rvert B}, \qquad \text{pitch}=\lvert v_{\parallel}\rvert T.

If v∥=0v_{\parallel}=0, the helix becomes a circle. If v⊥=0v_{\perp}=0, the path is straight along the field.


The force on a wire comes from the forces on its moving charge carriers. For a straight segment of length LL and cross-sectional area AA, the number of carriers is ncALn_cAL. Using their drift velocity and J⃗=ncqv⃗d\vec J=n_cq\vec v_d gives

F⃗=(ncAL)qv⃗d×B⃗=IL⃗×B⃗.\vec F=(n_cAL)q\vec v_d\times\vec B =I\vec L\times\vec B.

Here L⃗\vec L points along conventional current, so the same formula works for positive or negative charge carriers. In magnitude,

F=ILBsin⁡θ.F=ILB\sin\theta.

For a curved wire or a nonuniform field, add the forces on individual elements:

dF⃗=I dℓ⃗×B⃗,F⃗=I∫Cdℓ⃗×B⃗.d\vec F=I\,d\vec\ell\times\vec B, \qquad \vec F=I\int_Cd\vec\ell\times\vec B.

Use the applied field at the wire, not the wire’s own singular idealized field. The force direction follows the same cross-product rule as for a positive charge: index finger along current, middle finger along field, thumb along force. Do not interchange current and field.

A closed loop in a uniform field has zero net force because

F⃗=I(∮Cdℓ⃗)×B⃗=0.\vec F=I\left(\oint_Cd\vec\ell\right)\times\vec B=0.

The forces can still produce a torque.

Opposite magnetic forces on the sides of a rectangular current loop forming a torque about its center

Proof (Torque on a rectangular loop). Let a rectangular loop have sides aa and bb, so its area is A=abA=ab. First put the field in the plane of the loop, parallel to the sides of length aa. Those sides feel no force. The two sides of length bb feel opposite forces of magnitude IBbIBb.

Each force has lever arm a/2a/2 about the central pivot, giving

τ=(IBb)a2+(IBb)a2=IBab=IAB.\tau=(IBb)\frac a2+(IBb)\frac a2=IBab=IAB.

When the loop’s normal makes angle θ\theta with the field, the effective lever arm is reduced by sin⁡θ\sin\theta:

τ=IABsin⁡θ.\tau=IAB\sin\theta.

The angle is measured from the normal, not from the plane of the loop. The torque is largest when the field lies in the loop’s plane and zero when the normal is aligned with the field.

The opposite forces on a current loop can turn a rotor. With a fixed current direction, the magnetic torque tries to align the loop’s magnetic moment with the field; it does not keep driving the rotation in the same sense through a full turn.

DC motor with a current loop, brushes, and split-ring commutator between magnetic poles

A split-ring commutator reverses the current every half-turn. This reverses the loop’s magnetic moment at the appropriate time, keeping the driving torque in the same rotational direction. Inertia carries the rotor through the orientations where the torque is momentarily zero.


A small current loop behaves as a magnetic dipole. For a planar loop, define

m⃗=IAn^.\vec m=IA\hat n.

Curl your right-hand fingers along the current; your thumb gives n^\hat n. This works for any planar loop shape, not just a rectangle. For NN identical aligned turns,

m⃗=NIAn^.\vec m=NIA\hat n.

The unit is A⋅m2\text{A}\cdot\text{m}^2. Magnetic moment is often written μ⃗\vec\mu, but we use m⃗\vec m here to distinguish it from permeability μ\mu. Plain mm in the particle-motion section denotes mass.

In a uniform external field,

τ⃗=m⃗×B⃗,U=−m⃗⋅B⃗=−mBcos⁡θ.\vec\tau=\vec m\times\vec B, \qquad U=-\vec m\cdot\vec B=-mB\cos\theta.

The net force is zero, but the torque tends to align the moment with the field. Parallel alignment minimizes the energy; antiparallel alignment has zero torque but is unstable.

In a nonuniform field, a small dipole with fixed moment can experience a net force:

F⃗=∇(m⃗⋅B⃗).\vec F=\nabla(\vec m\cdot\vec B).

For m⃗=mx^\vec m=m\hat x and a field along x^\hat x on the dipole’s path,

Fx=mdBxdx.F_x=m\frac{dB_x}{dx}.

An aligned fixed dipole is pulled toward stronger field. These energy and force expressions use an externally imposed field and a fixed moment; do not apply the fixed-moment derivative blindly to an induced moment that itself changes with position.

For the circular loop of radius aa, m=Iπa2m=I\pi a^2. Far along its positive axis, where z≫az\gg a,

B⃗(z)≈μ0Ia22z3z^=μ0m2πz3z^.\vec B(z)\approx\frac{\mu_0 I a^2}{2z^3}\hat z =\frac{\mu_0 m}{2\pi z^3}\hat z.

At any point far from a localized dipole,

B⃗(r⃗)≈μ04πr3[3(m⃗⋅r^)r^−m⃗].\vec B(\vec r)\approx\frac{\mu_0}{4\pi r^3} \left[3(\vec m\cdot\hat r)\hat r-\vec m\right].

If the moment points along positive zz and θ\theta is the polar angle from that axis, this is

B⃗≈μ0m4πr3(2cos⁡θ r^+sin⁡θ θ^).\vec B\approx\frac{\mu_0 m}{4\pi r^3} \left(2\cos\theta\,\hat r+\sin\theta\,\hat\theta\right).

On the equatorial plane the field points opposite m⃗\vec m, with half the axial magnitude at the same distance. The dipole approximation requires distance much larger than the source’s size.

For a localized steady current distribution, the general magnetic moment is

m⃗=12∫Vr⃗×J⃗(r⃗) dV.\vec m=\frac12\int_V\vec r\times\vec J(\vec r)\,dV.

Atoms and molecules can have magnetic moments from electron orbital motion and intrinsic electron spin. Contributions often cancel, especially in filled shells. Spin is intrinsic angular momentum, not a little charged sphere literally rotating about an axis.

A magnet can be modeled as many microscopic dipoles. To find its total force or torque in an external field, add the contributions from those dipoles. Their average magnetic moment per unit volume is the magnetization:

M⃗=∑im⃗iΔV.\vec M=\frac{\sum_i\vec m_i}{\Delta V}.

The averaging volume is small on the scale of the object but contains many atoms. Magnetization has units A/m\text{A/m}.

In a paramagnetic material, an applied field weakly favors alignment of microscopic moments along the field. Thermal motion prevents full alignment, so the average magnetization is usually small. Removing the field removes the preferred direction, and the bulk magnetization normally disappears.

In materials such as iron, nickel, and cobalt, neighboring moments can strongly favor parallel alignment. Regions with aligned moments are called magnetic domains. An unmagnetized sample can contain strongly magnetized domains pointing in different directions, with little net magnetization.

An applied field favors domains oriented along it. Domain walls move and moments rotate, producing a much larger response than ordinary paramagnetism.

Magnetic domains before, during, and after application of a magnetic field, showing increased alignment and remanence
  • Hard magnetic materials resist demagnetization and can retain substantial alignment after the field is removed. They are useful for permanent magnets.
  • Soft magnetic materials are readily magnetized and demagnetized: their domain configuration changes under relatively small applied fields. This is not simply thermal randomization of all the moments.

An applied field also changes the orbital motion of electrons, inducing a magnetic response opposite to the applied field. This diamagnetic contribution occurs in all materials, but it is often hidden by stronger paramagnetic or ferromagnetic effects.

In materials such as copper and water, the diamagnetic response dominates. It is usually weak, although strong nonuniform fields can make its mechanical effects noticeable. Do not interpret the response as every electron simply beginning the same classical circular orbit; the orbital picture is a model for the induced opposing moment.

A superconductor supports persistent current without electrical resistance. In the Meissner state, screening currents near its surface expel magnetic flux from the bulk, so B⃗≈0\vec B\approx0 well inside. The field penetrates a thin surface layer rather than stopping at a mathematically sharp boundary.

Magnetic field lines diverted around a superconducting body by screening currents

Type I superconductors lose superconductivity above a critical field. Type II superconductors also have a Meissner state below a lower critical field; between their lower and upper critical fields they enter a mixed state in which flux penetrates in vortices. Thus, “type II always freezes the field inside” is not a general rule. Zero resistance alone does not imply flux expulsion; the Meissner effect is an additional property. See the superconductivity discussion in OpenStax for this distinction.


The microscopic current-loop model of magnetization can be replaced by equivalent bound currents. These describe the magnetic effect of the material’s dipoles, rather than a transport current supplied through a wire. The latter is a free current.

For magnetization M⃗\vec M, the bound surface and volume current densities are

K⃗b=M⃗×n^,J⃗b=∇×M⃗,\vec K_b=\vec M\times\hat n, \qquad \vec J_b=\nabla\times\vec M,

where n^\hat n points outward from the material. If M⃗\vec M is uniform inside the material, J⃗b=0\vec J_b=0 there: neighboring microscopic loops cancel internally, leaving an equivalent surface current.

Example. A cylinder has uniform magnetization M⃗=Mz^\vec M=M\hat z. Find its equivalent bound currents.

The volume current is zero because ∇×M⃗=0\nabla\times\vec M=0. On the curved surface, n^=r^\hat n=\hat r, so

K⃗b=Mz^×r^=Mϕ^.\vec K_b=M\hat z\times\hat r=M\hat\phi.

On the end faces, n^=±z^\hat n=\pm\hat z and K⃗b=0\vec K_b=0. The cylinder is therefore magnetically equivalent to an azimuthal surface current, like a solenoid winding.

Ampère’s law counts both types of current:

∮CB⃗⋅dℓ⃗=μ0(If,enc+Ib,enc).\oint_C\vec B\cdot d\vec\ell =\mu_0(I_{f,\text{enc}}+I_{b,\text{enc}}).

The corresponding magnetization circulation gives the enclosed bound current, with surface contributions included at material boundaries:

Ib,enc=∮CM⃗⋅dℓ⃗.I_{b,\text{enc}}=\oint_C\vec M\cdot d\vec\ell.

Moving this contribution to the left motivates the magnetic field strength H⃗\vec H:

H⃗=B⃗μ0−M⃗,∮CH⃗⋅dℓ⃗=If,enc.\vec H=\frac{\vec B}{\mu_0}-\vec M, \qquad \oint_C\vec H\cdot d\vec\ell=I_{f,\text{enc}}.

Both H⃗\vec H and M⃗\vec M have units A/m\text{A/m}. The field acting in the magnetic force law is still B⃗\vec B, not H⃗\vec H. In vacuum, M⃗=0\vec M=0 and B⃗=μ0H⃗\vec B=\mu_0\vec H.

For a linear, isotropic magnetic material,

M⃗=χmH⃗.\vec M=\chi_m\vec H.

The dimensionless constant χm\chi_m is the magnetic susceptibility. Substituting into the definition of H⃗\vec H gives

B⃗=μ0(H⃗+M⃗)=μ0(1+χm)H⃗=μH⃗,\vec B=\mu_0(\vec H+\vec M) =\mu_0(1+\chi_m)\vec H=\mu\vec H,

where

μ=μ0(1+χm),μr=μμ0=1+χm.\mu=\mu_0(1+\chi_m),\qquad \mu_r=\frac{\mu}{\mu_0}=1+\chi_m.

Paramagnets have χm>0\chi_m>0; diamagnets have χm<0\chi_m<0. For many ordinary weakly magnetic materials, ∣χm∣≪1\lvert\chi_m\rvert\ll1, so μ\mu is close to μ0\mu_0. The magnitude depends on material and temperature; a single small numerical range does not describe every paramagnet.

Ferromagnets can have a very large response, but a single constant χm\chi_m or μ\mu generally does not describe them. Their response is nonlinear and depends on the magnetization history.

Increasing and then decreasing the applied field does not take a ferromagnet through the same sequence of domain configurations. A plot of MM against HH, or BB against HH, traces a hysteresis loop.

Magnetization versus field-strength hysteresis loop showing saturation, remanent magnetization, and coercive field

At large fields, the magnetization approaches saturation. After the field returns to zero, some magnetization can remain: remanence. A reverse field is needed to reduce the magnetization to zero; its magnitude is the coercive field for the MM–HH loop. Hard magnets have high coercivity, while soft magnets have low coercivity. Distinguish an MM–HH graph from a BB–HH graph when reading its intercepts.