Boson magnetism: consider a gas of non-interacting spin 1 bosons, each subject to a Hamiltonian
$$
\mathcal{H}_1\left(\vec{p}, s_z\right)=\frac{\vec{p}^2}{2 m}-\mu_0 s_z B,
$$
where $\mu_0=e \hbar / m c$, and $s_z$ takes three possible values of $(-1,0,+1)$. (The orbital effect, $\vec{p} \rightarrow \vec{p}-e \vec{A}$, has been ignored.)
(a) In a grand canonical ensemble of chemical potential $\mu$, what are the average occupation numbers $\left\{\left(n_{+}(\vec{k})\right\rangle,\left\langle n_0(\vec{k})\right\rangle,\left\langle n_{-}(\vec{k})\right\rangle\right\}$ of one-particle states of wavenumber $\vec{k}=\vec{p} / \hbar$ ?
(b) Calculate the average total numbers $\left\{N_{+}, N_0, N_{-}\right\}$of bosons with the three possible values of $s_z$ in terms of the functions $f_m^{+}(z)$.
(c) Write down the expression for the magnetization $M(T, \mu)=\mu_0\left(N_{+}-N_{-}\right)$, and by expanding the result for small $B$ find the zero-field susceptibility $\chi(T, \mu)=$ $\partial M /\left.\partial B\right|_{B=0}$.
To find the behavior of $\chi(T, n)$, where $n=N / V$ is the total density, proceed as follows:
(d) For $B=0$, find the high-temperature expansion for $z(\beta, n)=\mathrm{e}^{\beta \mu}$, correct to second order in $n$. Hence obtain the first correction from quantum statistics to $\chi(T, n)$ at high temperatures.
(c) Find the temperature $T_c(n, B=0)$ of Bose-Einstein condensation. What happens to $\chi(T, n)$ on approaching $T_c(n)$ from the high-temperature side?
(f) What is the chemical potential $\mu$ for $T<T_c(n)$, at a small but finite value of $B$ ? Which one-particle state has a macroscopic occupation number?
(g) Using the result in (f), find the spontaneous magnetization
$$
\bar{M}(T, n)=\lim _{B \rightarrow 0} M(T, n, B)
$$
$* * * * * * * *$