Solution (source code)

= Solution

The surface <pressure> follows by differentiating the logarithm of the <canonical partition function> with respect to container area:
$$
\Pi=k_BT(\partial_a\ln Q)_{N,T}
=\frac{Nk_BT}{a-2Na_d}
=\frac{k_BT}{a_d}\frac{\phi}{1-2\phi}.
$$
Its formal divergence is therefore \b[$\phi=1/2$]. This lies outside the low-density regime used to derive the expression. It is a model pole, not an actual geometrical close-packing fraction.