  A diffuse interface (phase field) model for an electrochemical system is
developed. We describe the minimal set of components needed to model an
electrochemical interface and present a variational derivation of the governing
equations. With a simple set of assumptions: mass and volume constraints,
Poisson's equation, ideal solution thermodynamics in the bulk, and a simple
description of the competing energies in the interface, the model captures the
charge separation associated with the equilibrium double layer at the
electrochemical interface. The decay of the electrostatic potential in the
electrolyte agrees with the classical Gouy-Chapman and Debye-H\"uckel theories.
We calculate the surface energy, surface charge, and differential capacitance
as functions of potential and find qualitative agreement between the model and
existing theories and experiments. In particular, the differential capacitance
curves exhibit complex shapes with multiple extrema, as exhibited in many
electrochemical systems.
