{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2025,10,13]],"date-time":"2025-10-13T15:25:56Z","timestamp":1760369156924,"version":"build-2065373602"},"reference-count":49,"publisher":"MDPI AG","issue":"11","license":[{"start":{"date-parts":[[2013,10,24]],"date-time":"2013-10-24T00:00:00Z","timestamp":1382572800000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/3.0\/"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["Entropy"],"abstract":"<jats:p>Coarse-grained models of soft matter are usually combined with implicit solvent models that take the electrostatic polarizability into account via a dielectric background. In biophysical or nanoscale simulations that include water, this constant can vary greatly within the system. Performing molecular dynamics or other simulations that need to compute exact electrostatic interactions between charges in those systems is computationally demanding. We review here several algorithms developed by us that perform exactly this task. For planar dielectric surfaces in partial periodic boundary conditions, the arising image charges can be either treated with the MMM2D algorithm in a very efficient and accurate way or with the electrostatic layer correction term, which enables the user to use his favorite 3D periodic Coulomb solver. Arbitrarily-shaped interfaces can be dealt with using induced surface charges with the induced charge calculation (ICC*) algorithm. Finally, the local electrostatics algorithm, MEMD(Maxwell Equations Molecular Dynamics), even allows one to employ a smoothly varying dielectric constant in the systems. We introduce the concepts of these three algorithms and an extension for the inclusion of boundaries that are to be held fixed at a constant potential (metal conditions). For each method, we present a showcase application to highlight the importance of dielectric interfaces.<\/jats:p>","DOI":"10.3390\/e15114569","type":"journal-article","created":{"date-parts":[[2013,10,24]],"date-time":"2013-10-24T12:03:57Z","timestamp":1382616237000},"page":"4569-4588","update-policy":"https:\/\/doi.org\/10.3390\/mdpi_crossmark_policy","source":"Crossref","is-referenced-by-count":42,"title":["Efficient Algorithms for Electrostatic Interactions Including Dielectric Contrasts"],"prefix":"10.3390","volume":"15","author":[{"given":"Axel","family":"Arnold","sequence":"first","affiliation":[{"name":"Institute for Computational Physics, University of Stuttgart, Allmandring 3, Stuttgart 70569, Germany"}]},{"given":"Konrad","family":"Breitsprecher","sequence":"additional","affiliation":[{"name":"Institute for Computational Physics, University of Stuttgart, Allmandring 3, Stuttgart 70569, Germany"}]},{"given":"Florian","family":"Fahrenberger","sequence":"additional","affiliation":[{"name":"Institute for Computational Physics, University of Stuttgart, Allmandring 3, Stuttgart 70569, Germany"}]},{"given":"Stefan","family":"Kesselheim","sequence":"additional","affiliation":[{"name":"Institute for Computational Physics, University of Stuttgart, Allmandring 3, Stuttgart 70569, Germany"}]},{"given":"Olaf","family":"Lenz","sequence":"additional","affiliation":[{"name":"Institute for Computational Physics, University of Stuttgart, Allmandring 3, Stuttgart 70569, Germany"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-2739-310X","authenticated-orcid":false,"given":"Christian","family":"Holm","sequence":"additional","affiliation":[{"name":"Institute for Computational Physics, University of Stuttgart, Allmandring 3, Stuttgart 70569, Germany"}]}],"member":"1968","published-online":{"date-parts":[[2013,10,24]]},"reference":[{"key":"ref_1","doi-asserted-by":"crossref","first-page":"253","DOI":"10.1002\/andp.19213690304","article-title":"Die Berechnung optischer und elektrostatischer Gitterpotentiale","volume":"369","author":"Ewald","year":"1921","journal-title":"Ann. Phys."},{"key":"ref_2","doi-asserted-by":"crossref","first-page":"1924","DOI":"10.1063\/1.441285","article-title":"Electrostatic potentials and fields in infinite point charge lattices","volume":"74","author":"Heyes","year":"1981","journal-title":"J. Chem. Phys."},{"key":"ref_3","doi-asserted-by":"crossref","first-page":"27","DOI":"10.1098\/rspa.1980.0135","article-title":"Simulation of electrostatic systems in periodic boundary conditions. I. Lattice sums and dielectric constants","volume":"373","author":"Perram","year":"1980","journal-title":"Proc. R. Soc. Lond. A"},{"key":"ref_4","doi-asserted-by":"crossref","first-page":"57","DOI":"10.1098\/rspa.1980.0136","article-title":"Simulation of electrostatic systems in periodic boundary conditions. II. Equivalence of boundary conditions","volume":"373","author":"Perram","year":"1980","journal-title":"Proc. R. Soc. Lond. A"},{"key":"ref_5","doi-asserted-by":"crossref","unstructured":"Hockney, R.W., and Eastwood, J.W. (1988). Computer Simulation Using Particles, IOP.","DOI":"10.1201\/9781439822050"},{"key":"ref_6","doi-asserted-by":"crossref","first-page":"10089","DOI":"10.1063\/1.464397","article-title":"Particle Mesh Ewald: An N log(N) method for Ewald sums in large systems","volume":"98","author":"Darden","year":"1993","journal-title":"J. Chem. Phys."},{"key":"ref_7","doi-asserted-by":"crossref","first-page":"8577","DOI":"10.1063\/1.470117","article-title":"A smooth Particle Mesh Ewald method","volume":"103","author":"Essmann","year":"1995","journal-title":"J. Chem. Phys."},{"key":"ref_8","doi-asserted-by":"crossref","first-page":"7678","DOI":"10.1063\/1.477414","article-title":"How to mesh up Ewald sums. I. A theoretical and numerical comparison of various particle mesh routines","volume":"109","author":"Deserno","year":"1998","journal-title":"J. Chem. Phys."},{"key":"ref_9","doi-asserted-by":"crossref","first-page":"167","DOI":"10.1016\/S0010-4655(97)00090-8","article-title":"A comparison of Ewald summation techniques for planar surfaces","volume":"107","author":"Widmann","year":"1997","journal-title":"Comput. Phys. Commun."},{"key":"ref_10","doi-asserted-by":"crossref","first-page":"3155","DOI":"10.1063\/1.479595","article-title":"Ewald summation for systems with slab geometry","volume":"111","author":"Yeh","year":"1999","journal-title":"J. Chem. Phys."},{"key":"ref_11","doi-asserted-by":"crossref","first-page":"2496","DOI":"10.1063\/1.1491955","article-title":"Electrostatics in periodic slab geometries I","volume":"117","author":"Arnold","year":"2002","journal-title":"J. Chem. Phys."},{"key":"ref_12","doi-asserted-by":"crossref","first-page":"2503","DOI":"10.1063\/1.1491955","article-title":"Electrostatics in periodic slab geometries II","volume":"117","author":"Arnold","year":"2002","journal-title":"J. Chem. Phys."},{"key":"ref_13","doi-asserted-by":"crossref","first-page":"214103","DOI":"10.1063\/1.2918365","article-title":"Boundary conditions in local electrostatics algorithms","volume":"128","author":"Levrel","year":"2008","journal-title":"J. Chem. Phys."},{"key":"ref_14","doi-asserted-by":"crossref","first-page":"214102","DOI":"10.1063\/1.2918361","article-title":"Local monte carlo for electrostatics in anisotropic and nonperiodic geometries","volume":"128","author":"Thompson","year":"2008","journal-title":"J. Chem. Phys."},{"key":"ref_15","doi-asserted-by":"crossref","first-page":"1089","DOI":"10.1080\/00268978800101621","article-title":"Electrostatic potentials for thin layers","volume":"65","author":"Smith","year":"1988","journal-title":"Mol. Phys."},{"key":"ref_16","doi-asserted-by":"crossref","first-page":"154723","DOI":"10.1063\/1.2790428","article-title":"ICMMM2D: An accurate method to include planar dielectric interfaces via image charge summation","volume":"127","author":"Tyagi","year":"2007","journal-title":"J. Chem. Phys."},{"key":"ref_17","doi-asserted-by":"crossref","first-page":"204102","DOI":"10.1063\/1.3021064","article-title":"Electrostatic layer correction with image charges: A linear scaling method to treat slab 2D + h systems with dielectric interfaces","volume":"129","author":"Tyagi","year":"2008","journal-title":"J. Chem. Phys."},{"key":"ref_18","doi-asserted-by":"crossref","first-page":"154112","DOI":"10.1063\/1.3376011","article-title":"An iterative, fast, linear-scaling method for computing induced charges on arbitrary dielectric boundaries","volume":"132","author":"Tyagi","year":"2010","journal-title":"J. Chem. Phys."},{"key":"ref_19","doi-asserted-by":"crossref","first-page":"33","DOI":"10.1016\/j.cpc.2010.08.014","article-title":"Applying ICC* to DNA translocation. Effect of dielectric boundaries","volume":"182","author":"Kesselheim","year":"2011","journal-title":"Comput. Phys. Commun."},{"key":"ref_20","doi-asserted-by":"crossref","first-page":"196402","DOI":"10.1103\/PhysRevLett.88.196402","article-title":"Local simulation algorithms for coulombic interactions","volume":"88","author":"Maggs","year":"2002","journal-title":"Phys. Rev. Lett."},{"key":"ref_21","doi-asserted-by":"crossref","first-page":"3999","DOI":"10.1088\/0953-8984\/16\/38\/017","article-title":"Coulomb interactions via local dynamics: A molecular-dynamics algorithm","volume":"16","author":"Pasichnyk","year":"2004","journal-title":"J. Phys. Condens. Matter"},{"key":"ref_22","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1007\/978-3-642-32979-1_1","article-title":"ESPResSo 3.1\u2014Molecular Dynamics Software for Coarse-Grained Models","volume":"Volume 89","author":"Griebel","year":"2013","journal-title":"Meshfree Methods for Partial Differential Equations VI"},{"key":"ref_23","doi-asserted-by":"crossref","first-page":"704","DOI":"10.1016\/j.cpc.2005.10.005","article-title":"ESPResSo\u2014An extensible simulation package for research on soft matter systems","volume":"174","author":"Limbach","year":"2006","journal-title":"Comput. Phys. Commun."},{"key":"ref_24","first-page":"59","article-title":"Efficient Methods to Compute Long Range Interactions for Soft Matter Systems","volume":"Volume II","author":"Holm","year":"2005","journal-title":"Advanced Computer Simulation Approaches for Soft Matter Sciences II"},{"key":"ref_25","doi-asserted-by":"crossref","unstructured":"Jackson, J.D. (1999). Classical Electrodynamics, Wiley. [3rd ed.].","DOI":"10.1119\/1.19136"},{"key":"ref_26","doi-asserted-by":"crossref","first-page":"327","DOI":"10.1016\/S0010-4655(02)00586-6","article-title":"MMM2D: A fast and accurate summation method for electrostatic interactions in 2D slab geometries","volume":"148","author":"Arnold","year":"2002","journal-title":"Comput. Phys. Commun."},{"key":"ref_27","doi-asserted-by":"crossref","first-page":"324","DOI":"10.1016\/S0009-2614(02)00131-8","article-title":"A novel method for calculating electrostatic interactions in 2D periodic slab geometries","volume":"354","author":"Arnold","year":"2002","journal-title":"Chem. Phys. Lett."},{"key":"ref_28","doi-asserted-by":"crossref","first-page":"094107","DOI":"10.1063\/1.3216473","article-title":"Simulations of non-neutral slab systems with long-range electrostatic interactions in two-dimensional periodic boundary conditions","volume":"131","author":"Ballenegger","year":"2009","journal-title":"J. Chem. Phys."},{"key":"ref_29","doi-asserted-by":"crossref","unstructured":"Katsikadelis, J.T. (2002). Boundary Elements: Theory and Applications: Theory and Applications, Elsevier.","DOI":"10.1016\/B978-008044107-8\/50006-7"},{"key":"ref_30","doi-asserted-by":"crossref","first-page":"9480","DOI":"10.1039\/c2sm25615k","article-title":"Effects of dielectric mismatch and chain flexibility on the translocation barriers of charged macromolecules through solid state nanopores","volume":"8","author":"Kesselheim","year":"2012","journal-title":"Soft Matter"},{"key":"ref_31","doi-asserted-by":"crossref","first-page":"89","DOI":"10.1021\/nl052107w","article-title":"Salt dependence of ion transport and DNA translocation through solid-state nanopores","volume":"6","author":"Dekker","year":"2006","journal-title":"Nano Lett."},{"key":"ref_32","doi-asserted-by":"crossref","first-page":"209","DOI":"10.1038\/nnano.2007.27","article-title":"Solid-state nanopores","volume":"2","author":"Dekker","year":"2007","journal-title":"Nat. Nanotechnol."},{"key":"ref_33","doi-asserted-by":"crossref","first-page":"048102","DOI":"10.1103\/PhysRevLett.94.048102","article-title":"Asymmetric diffusion through synthetic nanopores","volume":"94","author":"Siwy","year":"2005","journal-title":"Phys. Rev. Lett."},{"key":"ref_34","doi-asserted-by":"crossref","first-page":"105105","DOI":"10.1063\/1.2358705","article-title":"Optical tweezers for force measurements on DNA in nanopores","volume":"77","author":"Keyser","year":"2006","journal-title":"Rev. Sci. Instrum."},{"key":"ref_35","doi-asserted-by":"crossref","first-page":"160","DOI":"10.1038\/nnano.2009.379","article-title":"Electrostatic focusing of unlabelled DNA into nanoscale pores using a salt gradient","volume":"5","author":"Wanunu","year":"2009","journal-title":"Nat. Nanotechnol."},{"key":"ref_36","unstructured":"Bastian, P., Blatt, M., Dedner, A., Engwer, C., Kl\u00f6fkorn, R., Kuttanikkad, S., Ohlberger, M., and Sander, O. (2006, January 12\u201314). The Distributed and Unified Numerics Environment (DUNE). Proceedings of the 19th Symposium on Simulation Technique, Hannover, Germany."},{"key":"ref_37","first-page":"294","article-title":"Generic implementation of finite element methods in the Distributed and Unified Numerics Environment (DUNE)","volume":"2","author":"Bastian","year":"2010","journal-title":"Kybernetika"},{"key":"ref_38","doi-asserted-by":"crossref","first-page":"7694","DOI":"10.1063\/1.477415","article-title":"How to mesh up Ewald sums. II. An accurate error estimate for the Particle-Particle-Particle-Mesh algorithm","volume":"109","author":"Deserno","year":"1998","journal-title":"J. Chem. Phys."},{"key":"ref_39","doi-asserted-by":"crossref","first-page":"7687","DOI":"10.1021\/jp3008877","article-title":"New coarse-grained models of imidazolium ionic liquids for bulk and interfacial molecular simulations","volume":"116","author":"Merlet","year":"2012","journal-title":"J. Phys. Chem. C"},{"key":"ref_40","doi-asserted-by":"crossref","first-page":"022201","DOI":"10.1088\/0953-8984\/23\/2\/022201","article-title":"Superionic state in double-layer capacitors with nanoporous electrodes","volume":"23","author":"Kondrat","year":"2011","journal-title":"J. Phys. Condens. Matter"},{"key":"ref_41","doi-asserted-by":"crossref","first-page":"4549","DOI":"10.1021\/jp809900w","article-title":"Microstructure and capacitance of the electrical double layers at the interface of ionic liquids and planar electrodes","volume":"113","author":"Feng","year":"2009","journal-title":"J. Phys. Chem. C"},{"key":"ref_42","doi-asserted-by":"crossref","first-page":"3108","DOI":"10.1063\/1.1642587","article-title":"Auxilary field Monte Carlo for charged particles","volume":"120","author":"Maggs","year":"2004","journal-title":"J. Chem. Phys."},{"key":"ref_43","doi-asserted-by":"crossref","first-page":"170201","DOI":"10.1103\/PhysRevLett.93.170201","article-title":"Local molecular dynamics with Coulombic interactions","volume":"93","author":"Rottler","year":"2004","journal-title":"Phys. Rev. Lett."},{"key":"ref_44","unstructured":"Fahrenberger, F., and Holm, C. Computing Coulomb Interaction in Inhomogeneous Dielectric Media via a Local Electrostatics Lattice Algorithm. Available online: http:\/\/arxiv.org\/abs\/1309.7859."},{"key":"ref_45","doi-asserted-by":"crossref","first-page":"1761","DOI":"10.1021\/jp077401m","article-title":"Simulating van der Waals interactions in water\/hydrocarbon-based complex fluids","volume":"112","author":"Pasichnyk","year":"2008","journal-title":"J. Phys. Chem. B"},{"key":"ref_46","doi-asserted-by":"crossref","first-page":"166102","DOI":"10.1103\/PhysRevLett.107.166102","article-title":"Dielectric profile of interfacial water and its effect on double-layer capacitance","volume":"107","author":"Bonthuis","year":"2011","journal-title":"Phys. Rev. Lett."},{"key":"ref_47","doi-asserted-by":"crossref","first-page":"67","DOI":"10.1016\/0021-9797(79)90222-4","article-title":"Electrophoretic mobility of a cylinder with high surface charge density","volume":"71","author":"Drift","year":"1979","journal-title":"J. Colloid Interface Sci."},{"key":"ref_48","doi-asserted-by":"crossref","first-page":"7679","DOI":"10.1021\/la2051564","article-title":"Profile of the static permittivity tensor of water at interfaces: Consequences for capacitance, hydration interaction and ion adsorption","volume":"28","author":"Bonthuis","year":"2012","journal-title":"Langmuir"},{"key":"ref_49","doi-asserted-by":"crossref","unstructured":"Arnold, A., Bolten, M., Dachsel, H., Fahrenberger, F., G\u00e4hler, F., Halver, R., Heber, F., Hofmann, M., Holm, C., and Iseringhausen, J. (2013). Comparison of scalable fast methods for long-range interactions. Phys. Rev. E, submitted.","DOI":"10.1103\/PhysRevE.88.063308"}],"container-title":["Entropy"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/www.mdpi.com\/1099-4300\/15\/11\/4569\/pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,10,11]],"date-time":"2025-10-11T21:50:07Z","timestamp":1760219407000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.mdpi.com\/1099-4300\/15\/11\/4569"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2013,10,24]]},"references-count":49,"journal-issue":{"issue":"11","published-online":{"date-parts":[[2013,11]]}},"alternative-id":["e15114569"],"URL":"https:\/\/doi.org\/10.3390\/e15114569","relation":{},"ISSN":["1099-4300"],"issn-type":[{"type":"electronic","value":"1099-4300"}],"subject":[],"published":{"date-parts":[[2013,10,24]]}}}