First-principles calculations combining density-functional theory and continuum solvation models enable realistic theoretical modeling and design of electrochemical systems. When a reaction proceeds in such systems, the number of electrons in the portion of the system treated quantum mechanically changes continuously, with a balancing charge appearing in the continuum electrolyte. A grand-canonical ensemble of electrons at a chemical potential set by the electrode potential is therefore the ideal description of such systems that directly mimics the experimental condition. We present two distinct algorithms, a self-consistent field method (GC-SCF) and a direct variational free energy minimization method using auxiliary Hamiltonians (GC-AuxH), to solve the Kohn-Sham equations of electronic density-functional theory directly in the grand canonical ensemble at fixed potential. Both methods substantially improve performance compared to a sequence of conventional fixed-number calculations targeting the desired potential, with the GC-AuxH method additionally exhibiting reliable and smooth exponential convergence of the grand free energy. Finally, we apply grand-canonical DFT to the under-potential deposition of copper on platinum from chloride-containing electrolytes and show that chloride desorption, not partial copper monolayer formation, is responsible for the second voltammetric peak.

/pdf/grand-canonical-electronic-density-functional-theory-51vtskpzj4.pdf

Grand canonical electronic density-functional theory: algorithms and applications to electrochemistry

Corrosion Mechanisms for Zinc Exposed to the Atmosphere

A detailed atomic-scale description of the electrochemical interface is essential to the understanding of electrochemical energy transformations. In this work, we investigate the charge of solvated protons at the Pt(111) | H2O and Al(111) | H2O interfaces. Using semi-local density-functional theory as well as hybrid functionals and embedded correlated wavefunction methods as higher-level benchmarks, we show that the effective charge of a solvated proton in the electrochemical double layer or outer Helmholtz plane at all levels of theory is fractional, when the solvated proton and solvent band edges are aligned correctly with the Fermi level of the metal (EF). The observed fractional charge in the absence of frontier band misalignment arises from a significant overlap between the proton and the electron density from the metal surface, and results in an energetic difference between protons in bulk solution and those in the outer Helmholtz plane.

Sheng Sun

Papers

Potential and distance dependent charge transfer: Hybrid first-principles/continuum calculations