A hybrid molecular dynamics/fluctuating hydrodynamics method for modelling liquids at multiple scales in space and time

A new 3D implementation of a hybrid model based on the analogy with two-phase hydrodynamics has been developed for the simulation of liquids at microscale. The idea of the method is to smoothly combine the atomistic description in the molecular dynamics zone with the Landau-Lifshitz fluctuating hydrodynamics representation in the rest of the system in the framework of macroscopic conservation laws through the use of a single "zoom-in" user-defined function s that has the meaning of a partial concentration in the two-phase analogy model. In comparison with our previous works, the implementation has been extended to full 3D simulations for a range of atomistic models in GROMACS from argon to water in equilibrium conditions with a constant or a spatially variable function s. Preliminary results of simulating the diffusion of a small peptide in water are also reported.

Two-phase flow analogy as an effective boundary condition for modelling liquids at atomistic resolution

A hybrid Molecular Dynamics/Fluctuating Hydrodynamics framework based on the analogy with two-phase hydrodynamics has been extended to dynamically tracking the feature of interest at all-atom resolution. In the model, the hydrodynamics description is used as an effective boundary condition to close the molecular dynamics solution without resorting to standard periodic boundary conditions. The approach is implemented in a popular Molecular Dynamics package GROMACS and results for two biomolecular systems are reported. A small peptide dialanine and a complete capsid of a virus porcine circovirus 2 in water are considered and shown to reproduce the structural and dynamic properties compared to those obtained in theory, purely atomistic simulations, and experiment.