Molecular oxygen has a net magnetic spin $\vec{S}$ of unity, that is, $S^z$ is quantized to $-1,0$, or +1 . The Hamiltonian for an ideal gas of $N$ such molecules in a magnetic field $\vec{B} \| \hat{z}$ is
$$\mathcal{H}=\sum_{i=1}^N\left[\frac{\vec{p}_i^2}{2 m}-\mu B S_i^z\right],$$
where $\left\{\vec{p}_i\right\}$ are the center of mass momenta of the molecules. The corresponding coordinates $\left\{\vec{q}_i\right\}$ are confined to a volume $V$. (Ignore all other degrees of freedom.)
(a) Treating $\left\{\vec{p}_i, \vec{q}_i\right\}$ classically, but the spin degrees of freedom as quantized, calculate the partition function, $\tilde{Z}(T, N, V, B)$.
(b) What are the probabilities for $S_i^z$ of a specific molecule to take on values of $-1,0$, +1 at a temperature $T$ ?
(c) Find the average magnetic dipole moment, $\langle M\rangle / V$, where $M=\mu \sum_{i=1}^N S_i^z$.
(d) Calculate the zero field susceptibility $\chi=\partial\langle M\rangle /\left.\partial B\right|_{B=0}$.