A energia de interação de um momento magnético $\vec{\mu}$ com o campo $\vec{B} = B_z \hat{k}$ é dada por $E = -\vec{\mu} \cdot \vec{B}$. Com a quantização, temos $2s+1$ níveis igualmente espaçados:
\[ E_m = g \mu_B B_z m \quad \text{onde} \quad m \in \{-s, \dots, +s\} \]
A função de partição para uma única partícula ($Z_1$) é a soma sobre todos os estados de spin acessíveis, resultando em uma série geométrica finita:
\[ Z_1(x) = \sum_{m=-s}^{s} e^{-\beta g \mu_B B_z m} = \frac{\sinh(ax)}{\sinh(bx)} \]
Onde definimos as variáveis adimensionais: $x = \frac{g \mu_B s B_z}{k_B T}$, $a = 1 + \frac{1}{2s}$ e $b = \frac{1}{2s}$. Note que $x$ representa a razão entre energia magnética e térmica.