An efficient method to establish electrostatic screening lengths of restricted primitive model electrolytes†

We present a novel, and computationally cheap, way to estimate electrostatic screening lengths from simulations of restricted primitive model (RPM) electrolytes. We demonstrate that the method is accurate by comparisons with simulated long-ranged parts of the charge density, at various Bjerrum lengths, salt concentrations and ion diameters. We find substantial underscreening in low dielectric solvent, but with an “aqueous” solvent, there is instead overscreening, the degree of which increases with ion size. Our method also offers a possible path to (future) more accurate classical density functional treatments of ionic fluids.