USHER: an algorithm for particle insertion in dense fluids

Abstract

The insertion of solvent particles in molecular dynamics simulations of complex fluids is required in many situations involving open systems, but this challenging task has been scarcely explored in the literature. We propose a simple and fast algorithm (USHER) that inserts the new solvent particles at locations where the potential energy has the desired prespecified value. For instance, this value may be set equal to the system's excess energy per particle, in such way that the inserted particles are energetically indistinguishable from the other particles present. During the search for the insertion site, the USHER algorithm uses a steepest descent iterator with a displacement whose magnitude is adapted to the local features of the energy landscape. The only adjustable parameter in the algorithm is the maximum displacement and we show that its optimal value can be extracted from an analysis of the structure of the potential energy landscape. We present insertion tests in periodic and non-periodic systems filled with a Lennard-Jones fluid whose density ranges from moderate values to high values.

10 pages (Latex), 8 figures (postscript); J. Chem. Phys. (in press) 2003

References in corpus (1)

• Continuum-particle hybrid coupling for mass, momentum and energy transfers in unsteady fluid flow

Cited by in corpus (9)

• Molecular Dynamics Simulations of Solutions at Constant Chemical Potential
• Multiscale modelling of liquids with molecular specificity
• Coupling atomistic and continuum hydrodynamics through a mesoscopic model: application to liquid water
• Coupling Lattice Boltzmann and Molecular Dynamics models for dense fluids
• Comparison of Molecular Dynamics with Hybrid Continuum-Molecular Dynamics for a Single Tethered Polymer in a Solvent
• Energy controlled insertion of polar molecules in dense fluids
• Moving Contact Lines: Linking Molecular Dynamics and Continuum-Scale Modeling
• Determination of the chemical potential using energy-biased sampling
• Multiscale Simulation of Surface Nanostructure Effect on Bubble Nucleation