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.
Subject(s)
Algorithms , Hydrodynamics , Molecular Dynamics Simulation , Biological Transport , Diffusion , ViscosityABSTRACT
Molecular dynamics simulations were performed for the prediction of the finite-size effects of Maxwell-Stefan diffusion coefficients of molecular mixtures and a wide variety of binary Lennard-Jones systems. A strong dependency of computed diffusivities on the system size was observed. Computed diffusivities were found to increase with the number of molecules. We propose a correction for the extrapolation of Maxwell-Stefan diffusion coefficients to the thermodynamic limit, based on the study by Yeh and Hummer ( J. Phys. Chem. B , 2004 , 108 , 15873 - 15879 ). The proposed correction is a function of the viscosity of the system, the size of the simulation box, and the thermodynamic factor, which is a measure for the nonideality of the mixture. Verification is carried out for more than 200 distinct binary Lennard-Jones systems, as well as 9 binary systems of methanol, water, ethanol, acetone, methylamine, and carbon tetrachloride. Significant deviations between finite-size Maxwell-Stefan diffusivities and the corresponding diffusivities at the thermodynamic limit were found for mixtures close to demixing. In these cases, the finite-size correction can be even larger than the simulated (finite-size) Maxwell-Stefan diffusivity. Our results show that considering these finite-size effects is crucial and that the suggested correction allows for reliable computations.
ABSTRACT
Diffusion is slow. Thus, diffusion experiments are intrinsically time-consuming and laborious. Additionally, the experimental effort is multiplied for multicomponent systems as the determination of multicomponent diffusion coefficients typically requires several experiments. To reduce the experimental effort, we present the first microfluidic diffusion measurement method for multicomponent liquid systems. The measurement setup combines a microfluidic chip with Raman microspectroscopy. Excellent agreement between experimental results and literature data is achieved for the binary system cyclohexane + toluene and the ternary system 1-propanol + 1-chlorobutane + heptane. The Fick diffusion coefficients are obtained from fitting a multicomponent convection-diffusion model to the mole fractions measured in experiments. Ternary diffusion coefficients can be obtained from a single experiment; high accuracy is already obtained from two experiments. Advantages of the presented measurement method are thus short measurement times, reduced sample consumption, and less experiments for the determination of a multicomponent diffusion coefficient.
