TTCF4LAMMPS: A toolkit for simulation of the non-equilibrium behaviour of molecular fluids at experimentally accessible shear rates

Abstract

We present TTCF4LAMMPS, a toolkit for performing non-equilibrium molecular dynamics (NEMD) simulations to study the fluid behaviour at low shear rates using the LAMMPS software. By combining direct NEMD simulations and the transient-time correlation function (TTCF) technique, we study the fluid response to shear rates spanning 15 orders of magnitude. We present two examples for simple monatomic systems: one consisting of a bulk liquid and another with a liquid layer confined between two solid walls. The small bulk system is suitable for testing on personal computers, while the larger confined system requires high-performance computing (HPC) resources. We demonstrate that the TTCF formalism can successfully detect the system response for arbitrarily weak external fields. We provide a brief mathematical explanation for this feature. Although we showcase the method for simple monatomic systems, TTCF can be readily extended to study more complex molecular fluids. Moreover, in addition to shear flows, the method can be extended to investigate elongational or mixed flows as well as thermal or electric fields. The high computational cost needed for the method is offset by the two following benefits: i) the cost is independent of the magnitude of the external field, and ii) the simulations can be made highly efficient on HPC architectures by exploiting the parallel design of the algorithm. We expect the toolkit to be useful for computational researchers striving to study the nonequilibrium behaviour of fluids under experimentally-accessible conditions. Program summaryProgram title: TTCF4LAMMPSCPC Library link to program files:https://doi.org/10.17632/hh2rkcxbrf.1Developer's repository link:https://github.com/edwardsmith999/TTCF4LAMMPSLicensing provisions: GNU General Public License 3Programming language: Python 3Nature of problem: Measuring the nonequilibrium behaviour of bulk and confined fluids under experimentally accessible strain rates in non-equilibrium molecular dynamics (NEMD) simulations.Solution method: Creating a Python-based code that utilises the transient-time correlation function method and the LAMMPS software to enable the bulk fluid properties (e.g. viscosity) and confined fluid interfacial properties (e.g. shear stress and slip velocity) to be computed at low shear rates with NEMD.