“…Furthermore, many materials interfaces are characterized by a discontinuity in the permittivity, resulting in induced surface polarization charges. This induced charge may significantly affect the ionic distributions [21,22,23,24,25,26] and dynamics [27,28,29], so that developing efficient methods for modeling spatially varying dielectrics has received considerable research attention [30,31,32,33,34,35,36]. Whereas this is typically costly, the special situation of planar interfaces (as well as systems that can be approximated as such) can be treated highly efficiently and accurately via the wellknown image-charge method (ICM) [37].…”