Solution (source code)

= Solution

For a dilute stack, the membrane number per unit normal length is $1/d=\phi/\delta$. Dividing the fluctuation repulsion per area by the repeat distance gives its <free-energy density>
$$
f_{\rm rep}
=\frac{c(k_BT)^2}{k_cd^3}
=\frac{c(k_BT)^2}{k_c\delta^3}\phi^3.
$$

The short-range electrostatic repulsion and long-range van der Waals attraction enter at second-virial order. In general, for an effective pair energy $U_{A_H}(\Gamma)$ over relative configurations $\Gamma$, the <second virial coefficient> has the Mayer-integral form
$$
\boxed{
B_2(A_H)=\frac12\int
\left[1-e^{-U_{A_H}(\Gamma)/(k_BT)}\right]d\Gamma},
$$
with a fixed normalization by the microscopic membrane thickness making $B_2$ dimensionless here. The mean-field contribution is quadratic in membrane concentration. The required two-power free energy can therefore be written
$$
\boxed{
f(\phi)=
\frac{k_BT}{\delta^3}B_2(A_H)\phi^2
+\frac{c(k_BT)^2}{k_c\delta^3}\phi^3,
\qquad \phi\geq0}.
$$
Changes of microscopic normalization merely rescale $B_2$ by a positive constant and do not affect the transition.