TY - JOUR
T1 - Incorporating surface polarization effects into large-scale coarse-grained Molecular Dynamics simulation
AU - Nguyen, Trung Dac
AU - Li, Honghao
AU - Bagchi, Debarshee
AU - Solis, Francisco J.
AU - Olvera de la Cruz, Monica
N1 - Funding Information:
TDN thanks Meng Shen for helpful discussion on the energy functional approach and Yamil Colón for help with SSAGES setup. TDN was supported by the Midwest Integrated Center for Computational Materials (MICCoM) . HL was supported by the National Science Foundation through Grant No. DMR-1611076 . MOdlC thanks the computational support of the Sherman Fairchild Foundation .
Publisher Copyright:
© 2019 Elsevier B.V.
PY - 2019/8
Y1 - 2019/8
N2 - We implement and extend three methods for incorporating surface polarization effects into coarse-grained Molecular Dynamics (MD) simulations for LAMMPS: the boundary element method using the generalized minimum residual (GMRES) solver for solving the Poisson equation, the induced charge computation method using the successive over relaxations scheme, and the direct optimization of an energy functional of induced charge density. Our implementations are validated against the analytical results for several cases and against published results on electrolyte and polyelectrolyte adsorption on charged interfaces. We have also employed the mentioned methods to examine the effects of surface polarization due to dielectric mismatch on the conformational behavior and relaxation times of the Rouse modes of a highly charged polyelectrolyte with explicit multivalent counterions in a spherical droplet surrounded by a low dielectric medium. We have found that dielectric mismatch coupled with spatial confinement results in noticeable changes in the radius of gyration of the collapsed conformation as well as the work required to compress and stretch the chain. Dielectric confinement is therefore believed to play an important role in the packaging of polyelectrolytes, which was generally overlooked in previous studies. Finally, the parallel performance of the implemented methods for practical simulations is characterized for a model electrolyte system. For polarizable interfaces composed of more than 1000 mesh points, the induced charge computation method is found to be the most efficient.
AB - We implement and extend three methods for incorporating surface polarization effects into coarse-grained Molecular Dynamics (MD) simulations for LAMMPS: the boundary element method using the generalized minimum residual (GMRES) solver for solving the Poisson equation, the induced charge computation method using the successive over relaxations scheme, and the direct optimization of an energy functional of induced charge density. Our implementations are validated against the analytical results for several cases and against published results on electrolyte and polyelectrolyte adsorption on charged interfaces. We have also employed the mentioned methods to examine the effects of surface polarization due to dielectric mismatch on the conformational behavior and relaxation times of the Rouse modes of a highly charged polyelectrolyte with explicit multivalent counterions in a spherical droplet surrounded by a low dielectric medium. We have found that dielectric mismatch coupled with spatial confinement results in noticeable changes in the radius of gyration of the collapsed conformation as well as the work required to compress and stretch the chain. Dielectric confinement is therefore believed to play an important role in the packaging of polyelectrolytes, which was generally overlooked in previous studies. Finally, the parallel performance of the implemented methods for practical simulations is characterized for a model electrolyte system. For polarizable interfaces composed of more than 1000 mesh points, the induced charge computation method is found to be the most efficient.
KW - Coarse-grained Molecular Dynamics
KW - Electrostatics
KW - LAMMPS
KW - Surface polarization
UR - http://www.scopus.com/inward/record.url?scp=85064595728&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85064595728&partnerID=8YFLogxK
U2 - 10.1016/j.cpc.2019.03.006
DO - 10.1016/j.cpc.2019.03.006
M3 - Article
AN - SCOPUS:85064595728
SN - 0010-4655
VL - 241
SP - 80
EP - 91
JO - Computer Physics Communications
JF - Computer Physics Communications
ER -