Curie's law

${\displaystyle M=\chi H,\chi =\frac{C}{T},}$

where

${\displaystyle \chi >0}$ is the (volume) magnetic susceptibility,

${\displaystyle M}$ is the magnitude of the resulting magnetization (A/m),

${\displaystyle H}$ is the magnitude of the applied magnetic field (A/m),

${\displaystyle T}$ is absolute temperature (K),

${\displaystyle C}$ is a material-specific Curie constant (K).

Pierre Curie discovered this relation, now known as Curie's law, by fitting data from experiment. It only holds for high temperatures and weak magnetic fields. As the derivations below show, the magnetization saturates in the opposite limit of low temperatures and strong fields. If the Curie constant is null, other magnetic effects dominate, like Langevin diamagnetism or Van Vleck paramagnetism.

Derivation with quantum mechanics

A simple model of a paramagnet concentrates on the particles which compose it which do not interact with each other. Each particle has a magnetic moment given by ${\displaystyle \overrightarrow{\mu }}$. The energy of a magnetic moment in a magnetic field is given by

${\displaystyle E=-\mathit{\mu }\cdot \mathbf{B},}$

where ${\displaystyle \mathbf{B}={\mu }_0(\mathbf{H}+\mathbf{M})}$ is the magnetic field density, measured in teslas (T).

Two-state (spin-1/2) particles

To simplify the calculation, we are going to work with a 2-state particle: it may either align its magnetic moment with the magnetic field or against it. So the only possible values of magnetic moment are then ${\displaystyle \mu }$ and ${\displaystyle -\mu }$. If so, then such a particle has only two possible energies, ${\displaystyle -\mu B}$ when it is aligned with the field and ${\displaystyle +\mu B}$ when it is oriented opposite to the field.

The extent to which the magnetic moments are aligned with the field can be calculated from the partition function. For a single particle, this is

${\displaystyle Z_1=\sum _{n=0,1}{e}^{-E_n\beta }=e^{\mu B\beta }+e^{-\mu B\beta }=2\cosh (\mu B\beta ).}$

The partition function for a set of N such particles, if they do not interact with each other, is

${\displaystyle Z=Z_1^N,}$

and the free energy is therefore

${\displaystyle G=-\frac{1}{\beta }\log Z=-N{k}_{\mathrm{B}}T\log Z_1.}$

The magnetization is the negative derivative of the free energy with respect to the applied field, and so the magnetization per unit volume is

${\displaystyle M=n\mu \tanh \frac{\mu B}{k_{\mathrm{B}}T},}$

where n is the number density of magnetic moments. The formula above is known as the Langevin paramagnetic equation. Pierre Curie found an approximation to this law that applies to the relatively high temperatures and low magnetic fields used in his experiments. As temperature increases and magnetic field decreases, the argument of the hyperbolic tangent decreases. In the Curie regime,

${\displaystyle \frac{\mu B}{k_{\mathrm{B}}T}\ll 1.}$

Moreover, if ${\displaystyle |x|\ll 1}$, then

${\displaystyle \tanh x\approx x,}$

so the magnetization is small, and we can write ${\displaystyle B\approx {\mu }_{0}H}$, and thus

${\displaystyle M\approx \frac{{\mu }_0{\mu }^{2}n}{k_{\mathrm{B}}}\frac{H}{T}.}$

In this regime, the magnetic susceptibility given by

${\displaystyle \chi =\frac{\mathrm{\partial }M}{\mathrm{\partial }H}\approx \frac{M}{H}}$

yields

${\displaystyle \chi (T\to \mathrm{\infty })=\frac{C}{T},}$

with a Curie constant given by ${\displaystyle C={\mu }_{0}n{\mu }^2/k_{\mathrm{B}}}$, in kelvins (K).

In the regime of low temperatures or high fields, ${\displaystyle M}$ tends to a maximum value of ${\displaystyle n\mu }$, corresponding to all the particles being completely aligned with the field. Since this calculation doesn't describe the electrons embedded deep within the Fermi surface, forbidden by the Pauli exclusion principle to flip their spins, it does not exemplify the quantum statistics of the problem at low temperatures. Using the Fermi–Dirac distribution, one will find that at low temperatures ${\displaystyle M}$ is linearly dependent on the magnetic field, so that the magnetic susceptibility saturates to a constant.

General case

When the particles have an arbitrary spin (any number of spin states), the formula is a bit more complicated. At low magnetic fields or high temperature, the spin follows Curie's law, with

${\displaystyle C=\frac{{\mu }_0{\mu }_{\text{B}}^2}{3k_{\mathrm{B}}}n{g}^{2}J(J+1),}$

where ${\displaystyle J}$ is the total angular momentum quantum number, and ${\displaystyle g}$ is the g-factor (such that ${\displaystyle \mu =gJ{\mu }_{\text{B}}}$ is the magnetic moment). For a two-level system with magnetic moment ${\displaystyle \mu }$, the formula reduces to $${\displaystyle C=\frac{1}{k_{\mathrm{B}}}n{\mu }_0{\mu }^2,}$$ as above, while the corresponding expressions in Gaussian units are $${\displaystyle C=\frac{{\mu }_{\mathrm{B}}^2}{3k_{\mathrm{B}}}n{g}^{2}J(J+1),}$$ $${\displaystyle C=\frac{1}{k_{\mathrm{B}}}n{\mu }^2.}$$

For this more general formula and its derivation (including high field, low temperature) see the article Brillouin function. As the spin approaches infinity, the formula for the magnetization approaches the classical value derived in the following section.