Spin-1 paramagnet
= Spin-1 paramagnet
For $N$ independent spins with $s_z\in\{-1,0,1\}$ and one-spin Hamiltonian $-\mu Bs_z$,
$$
Z=\left(1+2\cosh(\beta\mu B)\right)^N.
$$
Writing $x=\beta\mu B$, its heat capacity is
$$
C=Nk_Bx^2\frac{2(\cosh x+2)}
{(1+2\cosh x)^2}.
$$