Simulation of charged systems in heterogeneous dielectric media via a true energy functional

Vikram Jadhao, Francisco Solis, Monica Olvera De La Cruz

Research output: Contribution to journalArticle

50 Citations (Scopus)

Abstract

For charged systems in heterogeneous dielectric media, a key obstacle for 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, one that 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 nonmonotonic ionic profiles in the dielectric with a lower dielectric constant.

Original languageEnglish (US)
Article number223905
JournalPhysical Review Letters
Volume109
Issue number22
DOIs
StatePublished - Nov 27 2012

Fingerprint

molecular dynamics
permittivity
simulation
Poisson equation
emulsions
energy
free energy
electrostatics
polarization
liquids
profiles

ASJC Scopus subject areas

  • Physics and Astronomy(all)

Cite this

Simulation of charged systems in heterogeneous dielectric media via a true energy functional. / Jadhao, Vikram; Solis, Francisco; De La Cruz, Monica Olvera.

In: Physical Review Letters, Vol. 109, No. 22, 223905, 27.11.2012.

Research output: Contribution to journalArticle

@article{70c0dce5da6f49508c41e156581d5d80,
title = "Simulation of charged systems in heterogeneous dielectric media via a true energy functional",
abstract = "For charged systems in heterogeneous dielectric media, a key obstacle for 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, one that 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 nonmonotonic ionic profiles in the dielectric with a lower dielectric constant.",
author = "Vikram Jadhao and Francisco Solis and {De La Cruz}, {Monica Olvera}",
year = "2012",
month = "11",
day = "27",
doi = "10.1103/PhysRevLett.109.223905",
language = "English (US)",
volume = "109",
journal = "Physical Review Letters",
issn = "0031-9007",
publisher = "American Physical Society",
number = "22",

}

TY - JOUR

T1 - Simulation of charged systems in heterogeneous dielectric media via a true energy functional

AU - Jadhao, Vikram

AU - Solis, Francisco

AU - De La Cruz, Monica Olvera

PY - 2012/11/27

Y1 - 2012/11/27

N2 - For charged systems in heterogeneous dielectric media, a key obstacle for 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, one that 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 nonmonotonic ionic profiles in the dielectric with a lower dielectric constant.

AB - For charged systems in heterogeneous dielectric media, a key obstacle for 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, one that 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 nonmonotonic ionic profiles in the dielectric with a lower dielectric constant.

UR - http://www.scopus.com/inward/record.url?scp=84870209083&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84870209083&partnerID=8YFLogxK

U2 - 10.1103/PhysRevLett.109.223905

DO - 10.1103/PhysRevLett.109.223905

M3 - Article

VL - 109

JO - Physical Review Letters

JF - Physical Review Letters

SN - 0031-9007

IS - 22

M1 - 223905

ER -