Solution
= Solution
Impose the fixed <conserved order parameter>
$$
\int\phi\,d\mathbf r=\Phi
$$
with a <Lagrange multiplier> $\mu_0$. Stationarity of
$$
\mathcal F[\phi]-\mu_0\int\phi\,d\mathbf r
$$
under every local variation gives
$$
\frac{\delta\mathcal F}{\delta\phi(\mathbf r)}
=\mu(\mathbf r)=\mu_0.
$$
\b[Thus the <chemical potential> is spatially constant at constrained equilibrium.]