Determination of the effective dielectric constant from the accurate solution of the Poisson equation

Abstract

Constant dielectric (CD) and distance-dependent dielectric (DDD) functions are the most popular and widespread in the Molecular Mechanism simulations of large molecular systems. In this article, we present a simple procedure to derive an effective dielectric constant, εout,eff, for these two methods based on numerical solutions of the Poisson equation. It was found that because of the very approximate of the CD and DDD models there is no universal εout,eff which will work equally well for all molecular systems. For example, different MD trajectories of the same molecule can produce different optimal εout,effs. The DDD function was found to yield better agreement with the numerical solutions of the Poisson equation than a CD model does. The reason is that a DDD function gives a better description of the electrostatic interactions short distances the atoms. Another interesting finding of this study is that under certain conditions εout,eff can take negative values for a system of two atoms at a limited distance range, However, in principle, there is nothing to prevent the εout,eff from taking negative values for specific conformations of some molecules.