Computing the Helmholtz Capacitance of Charged Insulator-Electrolyte Interfaces from the Supercell Polarization

Abstract

Supercell modelling of an electrical double layer (EDL) at electrified solid-electrolyte interfaces is a challenge. The net polarization of EDLs arising from the fixed chemical composition setup leads to uncompensated EDLs under periodic boundary condition and convolutes the calculation of the Helmholtz capacitance [Zhang and Sprik, Phys. Rev. B, 94, 245309 (2016)]. Here we provide a new formula based on the supercell polarization at zero electric field $\bar{E}=0$ (i.e. standard Ewald boundary condition) to calculate the Helmholtz capacitance of charged insulator-electrolyte interfaces and validate it using atomistic simulations. Results are shown to be independent of the supercell size. This formula gives a shortcut to compute the Helmholtz capacitance without locating the zero net charge state of EDL and applies directly to any standard molecular dynamics code where the electrostatic interactions are treated by the Ewald summation or its variants.

References in corpus (5)

• Water orientation and hydrogen-bond structure at the fluorite/water interface
• Computing the dielectric constant of liquid water at constant dielectric displacement
• Finite Field Methods for the Supercell Modelling of Charged Insulator-Electrolyte Interfaces
• Note: On the dielectric constant of nanoconfined water
• Charge compensation at the interface between the polar NaCl(111) surface and a NaCl aqueous solution

Cited by in corpus (3)

• Coupling of surface chemistry and electric double layer at TiO${}_2$ electrochemical interfaces
• Stabilization of AgI's polar surfaces by the aqueous environment, and its implications for ice formation
• Macroscopic surface charges from microscopic simulations