Grand-canonical-like molecular-dynamics simulations by using an adaptive-resolution technique

Abstract

In this work, we provide a detailed theoretical analysis, supported by numerical tests, of the reliability of the adaptive resolution simulation (AdResS) technique in sampling the Grand Canonical ensemble. We demonstrate that the correct density and radial distribution functions in the hybrid region, where molecules change resolution, are two necessary conditions for considering the atomistic and coarse-grained regions in AdResS equivalent to subsystems of a full atomistic system with an accuracy up to the second order with respect to the probability distribution of the system. Moreover, we show that the work done by the thermostat and a thermodynamic force in the transition region is formally equivalent to balance the chemical potential difference between the different resolutions. From these results follows the main conclusion that the atomistic region exchanges molecules with the coarse-grained region in a Grand Canonical fashion with an accuracy up to (at least) second order. Numerical tests, for the relevant case of liquid water at ambient conditions, are carried out to strengthen the conclusions of the theoretical analysis. Finally, in order to show the computational convenience of AdResS as a Grand Canonical set up, we compare our method to the Insertion Particle Method (IMP) in its most efficient computational implementation. This fruitful combination of theoretical principles and numerical evidence candidates the adaptive resolution technique as a natural, general and efficient protocol for Grand Canonical Molecular Dynamics for the case of large systems.

References in corpus (5)

• Adaptive Resolution Molecular Dynamics Simulation: Changing the Degrees of Freedom on the Fly
• Adaptive resolution molecular dynamics simulation through coupling to an internal particle reservoir
• Coupling different levels of resolution in molecular simulations
• Coupling atomistic and continuum hydrodynamics through a mesoscopic model: application to liquid water
• Some fundamental problems for an energy conserving adaptive resolution molecular dynamics scheme

Cited by in corpus (21)

• Molecular Dynamics Simulations of Solutions at Constant Chemical Potential
• Monte Carlo adaptive resolution simulation of multicomponent molecular liquids
• Statistical Mechanics of Hamiltonian Adaptive Resolution Simulations
• Chemical potential of liquids and mixtures via Adaptive Resolution Simulation
• Molecular Dynamics in a Grand Ensemble: Bergmann-Lebowitz model and Adaptive Resolution Simulation
• Path Integral Molecular Dynamics within the Grand Canonical-like Adaptive Resolution Technique: Simulation of Liquid Water
• Open-Boundary Hamiltonian adaptive resolution. From grand canonical to non-equilibrium molecular dynamics simulations
• Formulation of Liouville's Theorem for Grand Ensemble Molecular Simulations
• Adaptive Resolution Molecular Dynamics Technique: Down to the Essential
• Open Systems out of Equilibrium: Theory and Simulation
• Probing Spatial Locality in Ionic Liquids with the Grand Canonical Adaptive Resolution Molecular Dynamics Technique
• Liouville-type equations for the n-particle distribution functions of an open system
• Adaptive coupling of a deep neural network potential to a classical force field
• Nanoscale domains in ionic liquids: A statistical mechanics definition for molecular dynamics studies
• From Classical to Quantum and Back: Hamiltonian Adaptive Resolution Path Integral, Ring Polymer, and Centroid Molecular Dynamics
• Simulation of Macromolecular Liquids with the Adaptive Resolution Molecular Dynamics Technique
• Simulation of a particle domain in a continuum/fluctuating hydrodynamics reservoir
• Adaptive Molecular Resolution Approach in Hamiltonian Form: An Asymptotic Analysis
• Molecular dynamics simulations with grand-canonical reweighting suggest cooperativity effects in RNA structure probing experiments
• 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