Grand Canonical Ensemble Approaches in CP2K for Modeling Electrochemistry at Constant Electrode Potentials.

作者：Ziwei Chai, S. Luber · 发表于：Journal of Chemical Theory and Computation · 年份：2024 · DOI：10.1021/acs.jctc.4c00671 · 被引用次数：11 · 研究领域：Medicine

In electrochemical experiments, the number of electrons of the electrode immersed in the electrolyte is usually variable. Additionally, the numbers of adsorbed substances on the surface of the electrode, the solvent molecules, and counter charge ions in the near-surface region can also vary. Treating electrochemical solid-liquid interfaces with the typical fixed electron number density functional theory (DFT) approach tends to be a challenge. This can be addressed by using grand canonical ensemble approaches. We present the implementation of two grand canonical ensemble approaches in the open-source computational chemistry software CP2K that go beyond the existing canonical ensemble paradigm. The first approach is based on implicit solvent models and explicit atomistic solute (electrode with/without adsorbed species) models, and includes two recent developments: (a) grand canonical self-consistent field (GC-SCF) method (J. Chem. Phys. 2017, 146, 114104) allowing the electron number of the system to fluctuate naturally and accordingly with the experimental electrode potential, (b) planar counter charge (J. Chem. Phys. 2019, 150, 041722, Phys. Rev. B 2003, 68, 245416) salt model completely screening the net charge of the electrode model. In contrast with previous studies, in our implementation, the work function (WF) (absolute electrode potential if the potential drop at the electrolyte-vacuum interface is omitted) is the constrained quantity during an SCF optimization instead ...