Revealing Local Field Effects at the Electrical Double Layer with Efficient Open Boundary Simulations under Potential Control

Margherita Buraschi,Andrew P. Horsfield,Clotilde S. Cucinotta
2023-11-06
Abstract:A major challenge in modelling interfacial processes in electrochemical (EC) devices is performing simulations at constant potential. This requires an open-boundary description of the electrons, so that they can enter and leave the computational cell. To enable realistic modelling of EC processes under potential control we have interfaced Density Functional Theory with the Hairy Probe method in the weak coupling limit (DOI: <a class="link-https link-external" data-doi="10.1103/PhysRevB.97.045116" href="https://doi.org/10.1103/PhysRevB.97.045116" rel="external noopener nofollow">https://doi.org/10.1103/PhysRevB.97.045116</a>). Our implementation was systematically tested using simple parallel-plate capacitor models with pristine surfaces and a single layer of adsorbed water molecules. Remarkably, our code's efficiency is comparable with a standard DFT calculation. We reveal that local field effects at the electrical double layer induced by the change of applied potential can significantly affect the energies of chemical steps in heterogeneous electrocatalysis. Our results demonstrate the importance of an explicit modelling of the applied potential in a simulation and provide an efficient tool to control this critical parameter.
Chemical Physics
What problem does this paper attempt to address?
The paper primarily aims to address the simulation issues of the Electrified Interface (EI) in Electrochemical (EC) devices, particularly how to effectively simulate open boundaries under constant potential. Specifically, the research focuses on the following aspects: 1. **Background and Challenges**: - The functionality of electrochemical devices largely depends on the EI, i.e., the formation of the Electrical Double Layer (EDL). - To reasonably design electrochemical devices, an atomic-level understanding of the EDL structure and its response to applied potential is required. - Although first-principles calculation methods can achieve electronic-level resolution, incorporating potential into computer simulations remains a challenge. 2. **Problems with Existing Methods**: - Most methods use Density Functional Theory (DFT) to simulate under the canonical ensemble (constant temperature, constant volume, constant number of particles). - The Computational Hydrogen Electrode (CHE) method, although successful in replicating the volcano plot of the hydrogen evolution reaction, cannot capture the potential-induced charging of the interface under operating conditions. - There is a conceptual difference between fixed charge simulation and fixed potential simulation, and a series of limitations exist, such as the need for a large number of samples or repeated simulations, and the lack of direct control over the electrode potential. 3. **Proposed Solution**: - Researchers have developed a novel and efficient method—Hairy-Probes Density Functional Theory (HP-DFT) to introduce electronic open boundaries. - This method combines the Hairy Probes (HP) formalism with the traditional Kohn-Sham DFT formalism, enabling the calculation of electronic structure and forces under applied potential control. - HP-DFT can achieve potential-controlled simulations while maintaining computational efficiency close to standard DFT. 4. **Research Findings**: - HP-DFT successfully simulated the parallel plate capacitor model and correctly described its physical properties, demonstrating the method's suitability for describing the behavior of open boundary systems under fixed potential. - HP-DFT was applied to study a simple water molecule bilayer/platinum interface system, revealing that surface polarization significantly affects coverage, structure, and adsorption energy even in non-electrochemical steps such as water adsorption. - The results indicate that a water bilayer does not accurately represent the behavior of a full water layer at the electrochemical interface, and systems with low coverage are more stable at zero bias. In summary, this paper addresses the challenge of simulating electrochemical interfaces under constant potential by proposing the efficient HP-DFT method and demonstrates its potential application in simulating electrochemical processes.