The potential energy of a system of charged particles, characterized by particle charge e and number density n(r), is given by
U=2e2∫∣r−r′∣n(r)n(r′)drdr′+e∫n(r)ϕext(r)dr,
where ϕext(r) is the potential of an external electric field. Assume that the entropy of the system, apart from an additive constant, is given by the formula
S=−k∫n(r)lnn(r)dr;
cf. formula (3.3.13). Using these expressions, derive the equilibrium equations satisfied by the number density n(r) and the total potential ϕ(r), the latter being
ϕext(r)+e∫∣r−r′∣n(r′)dr′.