GCMe: Efficient implementation of the Gaussian core model with smeared electrostatic interactions for molecular dynamics simulations of soft matter systems

Abstract

In recent years, molecular dynamics (MD) simulations have emerged as a pivotal tool for understanding the structure, dynamics, and phase behavior in charged soft matter systems. To explore phenomena across greater length and time scales in MD simulations, molecules are often coarse-grained for better computational performance. However, commonly-used force fields represent particles as hard-core interaction centers with point charges, which often overemphasizes the packing effect and short-range electrostatics, especially in systems with bulky deformable organic molecules and systems with strong coarse-graining. This underscores the need for an efficient soft-core model to physically capture the effective interactions between coarse-grained particles. To this end, we implement a soft-core model uniting the Gaussian core model with smeared electrostatic interactions that is phenomenologically equivalent to recent theoretical models. We first parameterize it generically using water as the model solvent. Then, we benchmark its performance in the OpenMM toolkit for different boundary conditions to highlight a computational speedup of up to $34\times$ compared to commonly used force fields and existing implementations. Finally, we demonstrate its utility by investigating how boundary polarizability affects the adsorption behavior of a polyelectrolyte solution on perfectly conducting and nonmetal electrodes.

References in corpus (14)

• The Electrostatic Screening Length in Concentrated Electrolytes Increases with Concentration
• Complex coacervation: A field theoretic simulation study of polyelectrolyte complexation
• Phase diagram of the Gaussian-core model
• Cooperative Dynamics in Unentangled Polymer Fluids
• A First Principle Approach to Rescale the Dynamics of Simulated Coarse-Grained Macromolecular Liquids
• Screening properties of Gaussian electrolyte models, with application to dissipative particle dynamics
• Intermolecular effects in the center-of-mass dynamics of unentangled polymer fluids
• Ultrasoft primitive model of polyionic solutions: structure, aggregation, and dynamics
• Theoretical Reconstruction of Realistic Dynamics of Highly Coarse-Grained cis-1,4-Polybutadiene Melts
• Asymmetric electrolytes near structured dielectric interfaces
• Modulation of ionic conduction using polarizable surfaces
• MATILDA.FT, a Mesoscale Simulation Package for Inhomogeneous Soft Matter
• Mediation of Long-Range Attraction Selectively between Negatively-Charged Colloids on Surfaces by Solvation
• Electric Fields Near Undulating Dielectric Membranes