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

Program summary: Program title: TTCF4LAMMPS CPC Library link to program files: https://doi.org/10.17632/hh2rkcxbrf.1 Developer's repository link: https://github.com/edwardsmith999/TTCF4LAMMPS Licensing provisions: GNU General Public License 3 Programming language: Python 3 Nature 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.Data availability: All data used in this work are freely reproducible using the program provided. The source files are also available at the link https://github.com/edwardsmith999/TTCF4LAMMPS.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.We thank the Australian Research Council for a grant obtained through the Discovery Projects Scheme (Grant No. DP200100422) and The Royal Society for support via International Exchanges, Grant No. IES/R3/170/233. J.P.E. was supported by the Royal Academy of Engineering (RAEng) through their Research Fellowships scheme. D.D. was supported through a Shell/RAEng Research Chair in Complex Engineering Interfaces. The authors acknowledge the Swinburne OzSTAR Supercomputing facility and the Imperial College London Research Computing Service (DOI:https://doi.org/10.14469/hpc/2232) for providing computational resources for this work. We thank Debra Bernhardt and Stephen Sanderson (University of Queensland) for useful discussions regarding the implementation of SLLOD in LAMMPS

Measuring heat flux beyond Fourier's law

We use nonequilibrium molecular dynamics (NEMD) to explore the effect of shear flow on heat flux. By simulating a simple fluid in a channel bounded by tethered atoms, the heat flux is computed for two systems: a temperature driven one with no flow and a wall driven, Couette flow system. The results for the temperature driven system give the Fourier's law thermal conductivity, which is shown to agree well with experiments. Through comparison of the two systems, we quantify the additional components of the heat flux parallel and normal to the walls due to shear flow. To compute the heat flux in the flow direction, the Irving-Kirkwood equations are integrated over a volume, giving the so-called volume average form, and they are also manipulated to get expressions for the surface averaged and method of planes forms. The method of planes and volume average forms are shown to give equivalent results for the heat flux when using small volumes. The heat flux in the flow direction is obtained consistently over a range of simulations, and it is shown to vary linearly with strain rate, as predicted by theory. The additional strain rate dependent component of the heat flux normal to the wall is obtained by fitting the strain rate dependence of the heat flux to the expected form. As a result, the additional terms in the thermal conductivity tensor quantified in this work should be experimentally testable.Comment: 16 pages double column, 11 figures, one tabl