A certain magnetic system contains n independent molecules per unit volume, each of which has four energy levels given by 0,Δ−gμBB,Δ,Δ+gμBB (g is a constant). Write down the partition function, compute the Helmholtz function and hence compute the magnetization M. Hence show that the magnetic susceptibility χ is given by
χ=B→0limBμ0M=kBT(3+eΔ/kBT)2nμ0g2μB2.(20.55)