A multiscale algorithm to drastically reduce computational times when simulating liquid/solid interaction at atomic resolution with realistic hydrodynamics effects

Abstract

A hybrid multiscale model that incorporates both continuum fluid dynamics effects including thermal fluctuations and molecular dynamics is presented for the challenging 3D liquid-solid interface problem of gap water flows between the moving tip of the Atomic Force Microscope (AFM) cantilever and a material substrate surface. Highlights of the method include the all-atom resolution in the vicinity of material surfaces, thereby avoiding the empiricism of interface boundary treatments of continuum mechanics approaches, with considerable computational savings in comparison with direct single-scale approaches such as Non-Equilibrium Molecular Dynamics (NEMD) methods. All components of the multiscale method are systematically described and verified. To optimise the coupling of the continuum fluid dynamics solution with the molecular dynamics solution a simple iterative method is used. Predictions of the multiscale method are compared with the reference equilibrium molecular dynamics solutions and discussed in the context of available AFM measurements in the literature. Computational performance of the suggested multiscale method is analysed in comparison with that of the NEMD method for a relevant range of the AFM problem parameters. In particular, in comparison with the single-resolution MD method, the implemented multiscale model leads to an eight order of magnitude acceleration of the solution for a slowly moving AFM tip typical of the existing experiments. Perspectives of further development and application of the suggested multiscale method are discussed. The current model is implemented in the open-source GROMACS software, and the source code of which is also provided.

More information

Identifiers

ORCID for Ivan Korotkin: orcid.org/0000-0002-5023-3684

Catalogue record

Export record