Molecular dynamics of open systems: construction of a mean-field particle reservoir

Abstract

The simulation of open molecular systems requires explicit or implicit reservoirs of energy and particles. Whereas full atomistic resolution is desired in the region of interest, there is some freedom in the implementation of the reservoirs. Here, we construct a combined, explicit reservoir by interfacing the atomistic region with regions of point-like, non-interacting particles (tracers) embedded in a thermodynamic mean field. The tracer molecules acquire atomistic resolution upon entering the atomistic region and equilibrate with this environment, while atomistic molecules become tracers governed by an effective mean-field potential after crossing the atomistic boundary. The approach is extensively tested on thermodynamic, structural, and dynamic properties of liquid water. Conceptual and numerical advantages of the procedure as well as new perspectives are highlighted and discussed.

accepted as communication in Advanced Theory and Simulations

References in corpus (2)

• Molecular Dynamics Simulations of Solutions at Constant Chemical Potential
• Nanoscale domains in ionic liquids: A statistical mechanics definition for molecular dynamics studies

Cited by in corpus (12)

• The Seven Deadly Sins: when computing crystal nucleation rates, the devil is in the details
• Neural functional theory for inhomogeneous fluids: Fundamentals and applications
• Perspective: How to overcome dynamical density functional theory
• Open Systems out of Equilibrium: Theory and Simulation
• Thermodynamic relations at the coupling boundary in adaptive resolution simulations for open systems
• Why hyperdensity functionals describe any equilibrium observable
• Finite-size corrections for the static structure factor of a liquid slab with open boundaries
• Simulation of a particle domain in a continuum/fluctuating hydrodynamics reservoir
• Nonequilibrium induced by reservoirs: Physico-mathematical models and numerical tests
• Continuous demixing transition of binary liquids: finite-size scaling from the analysis of sub-systems
• Determining the chemical potential via universal density functional learning
• Dynamics of systems with varying number of particles: from Liouville equations to general master equations for open systems