https://doi.org/10.1140/epjst/e2015-02415-x
Review
Open boundary molecular dynamics
1 Departamento Física Teórica de la Materia Condensada, Universidad Autónoma de Madrid, Campus de Cantoblanco, 28049 Madrid, Spain
2 Laboratory for Molecular Modeling, National Institute of Chemistry, Hajdrihova 19, 1001 Ljubljana, Slovenia
a e-mail: rafael.delgado@uam.es
b e-mail: praprot@cmm.ki.si
Received: 17 March 2015
Revised: 5 May 2015
Published online: 22 June 2015
This contribution analyzes several strategies and combination of methodologies to perform molecular dynamic simulations in open systems. Here, the term open indicates that the total system has boundaries where transfer of mass, momentum and energy can take place. This formalism, which we call Open Boundary Molecular Dynamics (OBMD), can act as interface of different schemes, such as Adaptive Resolution Scheme (AdResS) and Hybrid continuum-particle dynamics to link atomistic, coarse-grained (CG) and continuum (Eulerian) fluid dynamics in the general framework of fluctuating Navier-Stokes equations. The core domain of the simulation box is solved using all-atom descriptions. The CG layer introduced using AdResS is located at the outer part of the open box to make feasible the insertion of large molecules into the system. Communications between the molecular system and the outer world are carried out in the outer layers, called buffers. These coupling preserve momentum and mass conservation laws and can thus be linked with Eulerian hydro- dynamic solvers. In its simpler form, OBMD allows, however, to impose a local pressure tensor and a heat flux across the system's boundaries. For a one component molecular system, the external normal pressure and temperature determine the external chemical potential and thus the independent parameters of a grand-canonical ensemble simulation. Extended ensembles under non-equilibrium stationary states can also be simulated as well as time dependent forcings (e.g. oscillatory rheology). To illustrate the robustness of the combined OBMD-AdResS method, we present simulations of star-polymer melts at equilibrium and in sheared flow.
© EDP Sciences, Springer-Verlag, 2015