Simulations of Coulomb systems with slab geometry using an efficient 3d Ewald summation method
arXiv:1801.05735 · doi:10.1063/1.4945560
Abstract
We present a new approach to efficiently simulate electrolytes confined between infinite charged walls using a 3d Ewald summation method. The optimal performance is achieved by separating the electrostatic potential produced by the charged walls from the electrostatic potential of electrolyte. The electric field produced by 3d periodic images of the plates is constant, with the field produced by the transverse images of the charged plates canceling out. We show that under suitable renormalization, the non-neutral electrolyte confined between charged plates can be simulated using 3d Ewald summation with a correction that accounts for the conditional convergence of the resulting lattice sum. The new algorithm is at least an order of magnitude more rapid than the usual simulation methods for the slab geometry and can be further sped up by adopting Particle-Particle Particle-Mesh (P 3 M ) approach.
References in corpus (1)
Cited by in corpus (9)
- Charge neutrality breakdown in confined aqueous electrolytes: theory and simulation
- Simulations of Coulomb systems confined by polarizable surfaces using periodic Green functions
- Simulations of Polyelectrolyte Adsorption to a Dielectric Like-Charged Surface
- Simulations of ionic liquids confined by metal electrodes using periodic Green functions
- Lattice Model of an Ionic Liquid at an Electrified Interface
- Note: A pairwise form of the Ewald sum for non-neutral systems
- Widom insertion method in simulations with Ewald summation
- Adsorption isotherms of charged nanoparticles
- Electrolytes in regimes of strong confinements: surface charge modulations, osmotic equilibrium and electroneutrality