Static polarizability effects on counterion distributions near charged dielectric surfaces: A coarse-grained Molecular Dynamics study employing the Drude model

Coarse-grained implicit solvent Molecular Dynamics (MD) simulations have been used to investigate the structure of the vicinal layer of polarizable counterions close to a charged interface. The classical Drude oscillator model was implemented to describe the static excess polarizability of the ions. The electrostatic layer correction with image charges (ELCIC) method was used to include the effects of the dielectric discontinuity between the aqueous solution and the bounding interfaces for the calculation of the electrostatic interactions. Cases with one or two charged bounding interfaces were investigated. The counterion density profile in the vicinity of the interfaces with different surface charge values was found to depend on the ionic polarizability. Ionic polarization effects are found to be relevant for ions with high excess polarizability near surfaces with high surface charge.