Green's function for a spherical dielectric discontinuity and its application to simulation

Per Linse, Leo Lue

Research output: Contribution to journalArticle

6 Citations (Scopus)

Abstract

We present rapidly convergent expressions for the Green's function of the Poisson equation for spherically symmetric systems where the dielectric constant varies discontinuously in the radial direction. These expressions are used in Monte Carlo simulations of various electrolyte systems, and their efficiency is assessed. With only the leading term of the expansion included, a precision of the polarization energy of 0.01 kJ mol−1 or better was achieved, which is smaller than the statistical uncertainty of a typical simulation. The inclusion of the dielectric inhomogeneity leads to a 2.5-fold increase of the computational effort, which is modest for this type of model. The simulations are performed on six types of systems having either (i) a uniform surface charge distribution, (ii) a uniform volume charge distribution, or (iii) mobile ions, which were neutralized by mobile counterions. The ion density distributions are investigated for different dielectric conditions. These spatial distributions are discussed in terms of the importance of (i) the direct mean-field Coulomb interaction, (ii) the surface charge polarization at the dielectric discontinuity, and/or (iii) the change in the attractive Coulomb correlations.
LanguageEnglish
Article number044903
Number of pages12
JournalJournal of Chemical Physics
Volume140
Issue number4
DOIs
Publication statusPublished - Jan 2014

Fingerprint

Green's function
discontinuity
Green's functions
Charge distribution
Surface charge
charge distribution
Ions
Polarization
simulation
Poisson equation
polarization
Coulomb interactions
Spatial distribution
Electrolytes
density distribution
spatial distribution
inhomogeneity
Permittivity
interactions
electrolytes

Keywords

  • dielectrics
  • dielectric constant
  • surface charge
  • Green's function methods

Cite this

@article{f5671ac4ce9e488ab3565185bb5432f5,
title = "Green's function for a spherical dielectric discontinuity and its application to simulation",
abstract = "We present rapidly convergent expressions for the Green's function of the Poisson equation for spherically symmetric systems where the dielectric constant varies discontinuously in the radial direction. These expressions are used in Monte Carlo simulations of various electrolyte systems, and their efficiency is assessed. With only the leading term of the expansion included, a precision of the polarization energy of 0.01 kJ mol−1 or better was achieved, which is smaller than the statistical uncertainty of a typical simulation. The inclusion of the dielectric inhomogeneity leads to a 2.5-fold increase of the computational effort, which is modest for this type of model. The simulations are performed on six types of systems having either (i) a uniform surface charge distribution, (ii) a uniform volume charge distribution, or (iii) mobile ions, which were neutralized by mobile counterions. The ion density distributions are investigated for different dielectric conditions. These spatial distributions are discussed in terms of the importance of (i) the direct mean-field Coulomb interaction, (ii) the surface charge polarization at the dielectric discontinuity, and/or (iii) the change in the attractive Coulomb correlations.",
keywords = "dielectrics, dielectric constant, surface charge, Green's function methods",
author = "Per Linse and Leo Lue",
year = "2014",
month = "1",
doi = "10.1063/1.4862148",
language = "English",
volume = "140",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
number = "4",

}

Green's function for a spherical dielectric discontinuity and its application to simulation. / Linse, Per; Lue, Leo.

In: Journal of Chemical Physics, Vol. 140, No. 4, 044903, 01.2014.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Green's function for a spherical dielectric discontinuity and its application to simulation

AU - Linse, Per

AU - Lue, Leo

PY - 2014/1

Y1 - 2014/1

N2 - We present rapidly convergent expressions for the Green's function of the Poisson equation for spherically symmetric systems where the dielectric constant varies discontinuously in the radial direction. These expressions are used in Monte Carlo simulations of various electrolyte systems, and their efficiency is assessed. With only the leading term of the expansion included, a precision of the polarization energy of 0.01 kJ mol−1 or better was achieved, which is smaller than the statistical uncertainty of a typical simulation. The inclusion of the dielectric inhomogeneity leads to a 2.5-fold increase of the computational effort, which is modest for this type of model. The simulations are performed on six types of systems having either (i) a uniform surface charge distribution, (ii) a uniform volume charge distribution, or (iii) mobile ions, which were neutralized by mobile counterions. The ion density distributions are investigated for different dielectric conditions. These spatial distributions are discussed in terms of the importance of (i) the direct mean-field Coulomb interaction, (ii) the surface charge polarization at the dielectric discontinuity, and/or (iii) the change in the attractive Coulomb correlations.

AB - We present rapidly convergent expressions for the Green's function of the Poisson equation for spherically symmetric systems where the dielectric constant varies discontinuously in the radial direction. These expressions are used in Monte Carlo simulations of various electrolyte systems, and their efficiency is assessed. With only the leading term of the expansion included, a precision of the polarization energy of 0.01 kJ mol−1 or better was achieved, which is smaller than the statistical uncertainty of a typical simulation. The inclusion of the dielectric inhomogeneity leads to a 2.5-fold increase of the computational effort, which is modest for this type of model. The simulations are performed on six types of systems having either (i) a uniform surface charge distribution, (ii) a uniform volume charge distribution, or (iii) mobile ions, which were neutralized by mobile counterions. The ion density distributions are investigated for different dielectric conditions. These spatial distributions are discussed in terms of the importance of (i) the direct mean-field Coulomb interaction, (ii) the surface charge polarization at the dielectric discontinuity, and/or (iii) the change in the attractive Coulomb correlations.

KW - dielectrics

KW - dielectric constant

KW - surface charge

KW - Green's function methods

UR - http://scitation.aip.org/content/aip/journal/jcp

U2 - 10.1063/1.4862148

DO - 10.1063/1.4862148

M3 - Article

VL - 140

JO - Journal of Chemical Physics

T2 - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 4

M1 - 044903

ER -