 Solution

Write ${ϵ}_k=E_0(k+\frac{1}{2})$ and use Lagrange multipliers $\alpha ,\beta$ for the particle-number and energy constraints. The stationary equations areThus $n_k=e^{-\alpha }{e}^{-\beta {ϵ}_k}$. Finite normalization over infinitely many levels requires $\beta >0$. Put $q=e^{-\beta E_0}$, so $0<q<1$. The geometric series givesThereforeFor the energy constraint, use $n_k=N(1-q)q^k$ and $\sum k{q}^k=q/(1-q{)}^2$, obtained by differentiating the geometric series. This yieldsRecognizing the hyperbolic cotangent givesThe strict positivity of all occupations requires $E>N{E}_0/2$; the paper's positivity assumptions alone do not guarantee this feasibility condition. For admissible energy the multiplier is uniquely $\beta =E_0^{-1}\ln [(2E+N{E}_0)/(2E-N{E}_0)]>0$. Equality $E=N{E}_0/2$ is only the limiting state with all occupation in level zero, outside the stipulated strictly positive family; lower energy is impossible.

For rigor beyond formal multiplier equations, $p_k=n_k/N$ has fixed mean $E/(N{E}_0)-1/2$. The objective differs from $N$ times its information entropy only by the fixed constant $N-N\ln N$. The maximum entropy distribution on the nonnegative integers is uniquely geometric at that mean, confirming the stationary family and its maximum character.