Abstract: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 of the Fermi energy. The chemical potential of electrons (negative WF) is a natural variable of the grand potential in the GC ensemble of electronic states, and this method can easily achieve stable SCF convergence and obtain an electronic structure that precisely corresponds to a user-specified WF. The second approach referred to as the GC DFT molecular dynamics (DFT-MD) simulation scheme (Phys. Rev. Lett. 2002, 88, 213002, J. Chem. Phys. 2005, 122, 234505, J. Am. Chem. Soc. 2004, 126 (12), 3928-3938) is based on fully explicit modeling the solvent molecules and the ions and is used to calculate the electron chemical potential corresponding to an equilibrium electrochemical half-reaction (M(n+m)+ + ne- ⇌ Mm+) which involves DFT-MD, by allowing the number of electrons to vary during the DFT-MD simulation process. This opens the way for forefront electrochemical calculations in CP2K for a broad range of systems.

Widom insertion method in simulations with Ewald summation

Simulations of Coulomb systems with slab geometry using an efficient 3d Ewald summation method

Non-Ewald methods for evaluating the electrostatic interactions of charge systems: similarity and difference

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

Analytic gradients for the electrostatic embedding QM/MM in periodic boundary conditions using particle-mesh Ewald sums and electrostatic potential fitted charge operators

Kawamata-Viehweg vanishing and quint-canonical map of a Complex threefold

Lekner summations and Ewald summations for quasi-two dimensional systems

Ewald summation for ferroelectric perovskites with charges and dipoles

Latent Ewald summation for machine learning of long-range interactions

A formula and numerical study on Ewald 1D summation

A Simple and Efficient Lattice Summation Method for Metallic Electrodes in Constant Potential Molecular Dynamics Simulation

Grand-canonical molecular dynamics simulations powered by a hybrid 4D nonequilibrium MD/MC method: Implementation in LAMMPS and applications to electrolyte solutions

The ENUF Method -- Ewald Summation based on Non-Uniform Fast Fourier Transform: Implementation, Parallelization, and Application

Polarizable Continuum Models and Green's Function GW Formalism: On the Dynamics of the Solvent Electrons

Chemical potential of quadrupolar two-centre Lennard-Jones fluids by gradual insertion

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

Two Improved Electronegativity Equalization Methods for Charge Distribution in Large Scale Non-Uniform System

Pressure Calculation in Polar and Charged Systems using Ewald Summation: Results for the Extended Simple Point Charge Model of Water

On Validity of Constant pH Simulations

Ewald summation for ion-dipole mixture under dielectric confinement

Molecular Theories and Simulation of Ions and Polar Molecules in Water