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

Abstract: 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.

• The paper introduces GC-SCF and GC-AuxH algorithms that transform fixed-electron DFT into a grand-canonical framework for realistic electrochemical modeling.
• The GC-AuxH method achieves smooth, exponential convergence and reduces computational expense by roughly 50% compared to traditional approaches.
• These innovations bridge the gap between theory and experiment, enhancing simulations of key electrochemical processes like underpotential deposition.

Grand Canonical Electronic Density-Functional Theory: Algorithms and Applications to Electrochemistry

The paper "Grand canonical electronic density-functional theory: algorithms and applications to electrochemistry" by Sundararaman, Goddard III, and Arias focuses on enhancing the applicability and efficiency of density-functional theory (DFT) for modeling electrochemical systems. It introduces novel algorithms to adapt conventional fixed-electron-number DFT methods into a grand-canonical ensemble framework, thereby allowing the number of electrons to adjust automatically in response to electrochemical potentials, a scenario frequently encountered in experimental electrochemical systems.

Algorithms and Theoretical Insights

The authors present two distinct algorithms: a self-consistent field method (GC-SCF) and a variational free energy minimization method utilizing auxiliary Hamiltonians (GC-AuxH). Both are adaptations of standard methods used in electronic structure calculations but adjusted to directly target fixed potentials.

1. GC-SCF: This method involves iterative optimization where the electron density is updated using a Pulay mixing scheme, which they modified to suit the grand-canonical ensemble's characteristics. This scheme effectively alleviates the common issue of charge sloshing in metallic systems, where electron densities oscillate but do not converge.
2. GC-AuxH: This approach employs a variational principle in which the free energy is minimized directly. This method introduces auxiliary Hamiltonians that lead to effective convergence in calculating the total energy, even with variable electron occupations. Notably, the GC-AuxH method exhibits consistently smooth and exponential convergence, outperforming GC-SCF in challenging systems, as identified through computational benchmarks.

Numerical and Computational Evaluations

The application results of these methods demonstrate a notable improvement in computational efficiency over the existing practice of manually adjusting the electron number to match potential (a process referred to as the "Loop" method). Both grand-canonical methods reduce the computational cost by approximately half, with GC-AuxH achieving more stable convergence without oscillations in the electron number or grand free energy. This stability is crucial for accurately modeling electrochemical systems close to their experimental conditions.

Implications for Electrochemical Systems

Electrochemical phenomena, particularly under-potential deposition (UPD), are highly sensitive to electron potential and electrolyte composition. By effectively modeling electron exchange between electrode surfaces and electrolytes, these algorithms provide insights that are potentially transformative for designing and predicting electrocatalytic processes. The study focuses on Cu UPD on Pt(111) and provides a computational explanation that aligns well with experimental observations regarding chloride desorption and full monolayer formation.

Future Developments

The methodologies proposed pave the way for enhanced modeling of electrochemical systems where the electron potential can vary dynamically. This advancement holds promise for bridging the gap between theoretical predictions and experimental measurements in systems of critical importance for energy storage and conversion technologies.

Conclusion

This work represents a robust contribution to the computational electrochemistry field, providing tools to address long-standing challenges associated with modeling systems at varied electrochemical potentials. The algorithms extend the reach of DFT, enabling the simulation of realistic experimental conditions, and thereby propelling theoretical investigations into practical electrochemical applications. This is expected to spur further development in both computational methods and applied electrochemistry research.