  The Poisson-Boltzmann equation is often presented via a variational
formulation based on the electrostatic potential. However, the functional has
the defect of being non-convex. It can not be used as a local minimization
principle while coupled to other dynamic degrees of freedom. We formulate a
convex dual functional which is numerically equivalent at its minimum and which
is more suited to local optimization.
