OCTP: A Tool for On-the-Fly Calculation of Transport Properties of Fluids with the Order- n Algorithm in LAMMPS

Seyed Hossein Jamali, Ludger Wolff, Tim M. Becker, Mariëtte De Groen, Mahinder Ramdin, Remco Hartkamp, André Bardow, Thijs J.H. Vlugt, Othonas A. Moultos

Research output: Contribution to journalArticleScientificpeer-review

17 Citations (Scopus)
84 Downloads (Pure)

Abstract

We present a new plugin for LAMMPS for on-the-fly computation of transport properties (OCTP) in equilibrium molecular dynamics. OCTP computes the self- and Maxwell-Stefan diffusivities, bulk and shear viscosities, and thermal conductivities of pure fluids and mixtures in a single simulation. OCTP is the first implementation in LAMMPS that uses the Einstein relations combined with the order-n algorithm for the efficient sampling of dynamic variables. OCTP has low computational requirements and is easy to use because it follows the native input file format of LAMMPS. A tool for calculating the radial distribution function (RDF) of the fluid beyond the cutoff radius, while taking into account the system size effects, is also part of the new plugin. The RDFs computed from OCTP are needed to obtain the thermodynamic factor, which relates Maxwell-Stefan and Fick diffusivities. To demonstrate the efficiency of the new plugin, the transport properties of an equimolar mixture of water-methanol were computed at 298 K and 1 bar.

Original languageEnglish
Pages (from-to)1290-1294
JournalJournal of Chemical Information and Modeling
Volume59
Issue number4
DOIs
Publication statusPublished - 2019

Fingerprint Dive into the research topics of 'OCTP: A Tool for On-the-Fly Calculation of Transport Properties of Fluids with the Order- n Algorithm in LAMMPS'. Together they form a unique fingerprint.

Cite this