Vikram Jadhao, Francisco J. Solis, Monica Olvera de la Cruz
For charged systems in heterogeneous dielectric media, a key obstacle for fast molecular dynamics (MD) simulations is the need to solve the Poisson equation in the media. This obstacle can be bypassed using MD methods that treat the local polarization charge density as a dynamic variable, but such approaches require access to a true free energy functional; which evaluates to the equilibrium electrostatic energy at its minimum. In this letter, we derive the needed functional. As an application, we develop a Car-Parrinello MD method for the simulation of free charges present near a spherical emulsion droplet separating two immiscible liquids with different dielectric constants. Our results show the presence of non-monotonic ionic profiles in the dielectric with lower dielectric constant.
View original:
http://arxiv.org/abs/1208.4113
No comments:
Post a Comment