Analytic Theory of Finite-size Effects in Supercell Modelings of Charged Interfaces

Abstract

The Ewald3D sum with the tinfoil boundary condition (e3dtf) evaluates the electrostatic energy of a finite unit cell inside an infinitely periodic supercell. Although it has been used as a {\it de facto} standard treatment of electrostatics for simulations of extended polar or charged interfaces, the finite-size effect on simulated properties has yet to be fully understood. There is, however, an intuitive way to quantify the average effect arising from the difference between the e3dtf and Coulomb potentials on the response of mobile charges to contact surfaces with fixed charges and/or to an applied external electric field. While any charged interface formed by mobile countercharges that compensate the fixed charges slips upon the change of the acting electric field, the distance between a pair of oppositely charged interfaces is found to be nearly stationary, which allows an analytic finite-size correction to the amount of countercharges. Application of the theory to solvated electric double layers (insulator/electrolyte interfaces) predicts that the state of complete charge compensation is invariant with respect to solvent permittivities, which is confirmed by a proper analysis of simulation data in the literature.

arXiv comment8 pages, 5 figures

References in corpus (9)

• Electric displacement as the fundamental variable in electronic-structure calculations
• Connecting local structure to interface formation: a molecular scale van der Waals theory of nonuniform liquids
• MMM2D: A fast and accurate summation method for electrostatic interactions in 2D slab geometries
• Interplay of local hydrogen-bonding and long-ranged dipolar forces in simulations of confined water
• Finite Field Methods for the Supercell Modelling of Charged Insulator-Electrolyte Interfaces
• Simulations of Coulomb systems with slab geometry using an efficient 3d Ewald summation method
• Efficient solutions of self-consistent mean field equations for dewetting and electrostatics in nonuniform liquids
• Rigorous Error Bounds for Ewald Summation of Electrostatics at Planar Interfaces
• Note: A pairwise form of the Ewald sum for non-neutral systems

Cited by in corpus (7)

• Dielectric response with short-ranged electrostatics
• Superscalability of the random batch Ewald method
• The Symmetry-Preserving Mean Field Condition for Electrostatic Correlations in Bulk
• Macroscopic surface charges from microscopic simulations
• Infinite Boundary Terms and Pairwise Interactions: A Unified Framework for Periodic Coulomb Systems
• A theory for the stabilization of polar crystal surfaces by a liquid environment
• The Finite Coulomb Lattice Sum: A Resolution of Conditional Convergence through Exact Shape and Size