The Symmetry-Preserving Mean Field Condition for Electrostatic Correlations in Bulk

Abstract

Accurate simulations of a condensed system of ions or polar molecules are concerned with proper handlings of the involved electrostatics. For such a Coulomb system at a charged planar interface, the Coulomb interaction averaged over the lateral directions with preserved symmetry serves as a necessary constraint in building any accurate handling that reconciles a simulated singlet charge density with the corresponding macroscopic charge/dielectric response. At present, this symmetry-preserving mean-field (SPMF) condition represented in the reciprocal space, is conjectured to be necessary for a simulated bulk system to reproduce correctly the charge structure factor of the macroscopic bulk, as well. In this work, we further examine analytically the asymptotic behavior of the charge structure factor at small wavenumbers for an arbitrary charge-charge interaction. In light of our theoretical predictions, simulations with a length of nearly 0.1 micron are carried out to demonstrate that, typical efficient methods violating the SPMF condition indeed fail to capture the exact charge correlations at small wavenumbers for both ionic and polar systems. However, for both types of systems, these existing methods can be simply amended to match the SPMF condition and subsequently to probe precisely the electrostatic correlations at all length scales.

6 pages, 3 figures

References in corpus (3)

• Simulations of ionic liquids confined by metal electrodes using periodic Green functions
• 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 (8)

• Self-Consistent Determination of Long-Range Electrostatics in Neural Network Potentials
• A random batch Ewald method for charged particles in the isothermal-isobaric ensemble
• Machine-Learning Interatomic Potentials for Long-Range Systems
• Improved random batch Ewald method in molecular dynamics simulations
• Infinite Boundary Terms and Pairwise Interactions: A Unified Framework for Periodic Coulomb Systems
• Random batch sum-of-Gaussians algorithm for molecular dynamics simulations of Yukawa systems in three dimensions
• The Finite Coulomb Lattice Sum: A Resolution of Conditional Convergence through Exact Shape and Size
• Analytic Boundary Terms for Arbitrary Crystal Geometries and Direct-Sum Evaluation of Madelung Constants in Triclinic Lattices