Partial molar properties from molecular simulation using multiple linear regression

Partial molar volumes, energies, and enthalpies can be computed from <i>NpT</i>-Gibbs ensemble simulations through a post-processing procedure that leverages fluctuations in composition, total volume, and total energy of a simulation box. By recording the instantaneous box volumes <i>V</i> and instantaneous number of molecules Ni of each of <i>n</i> species for <i>M</i> frames, a large M×n matrix N is constructed, as well as the M×1 vector V. The 1×n vector of partial molar volumes V¯ may then be solved using N⋅V¯=V. A similar construction permits calculation of partial molar energies using <i>M</i> instantaneous measurements of the total energy of the simulation box, and N⋅U¯=U. Partial molar enthalpies may be derived from U¯, V¯, and pressure <i>p</i>. These properties may be used to calculate enthalpy and entropy of transfer (absorption, extraction, and adsorption) for species in complex mixtures. The method is demonstrated on three systems in the <i>NpT</i>-Gibbs ensemble: a highly compressible natural gas condensate of methane, <i>n</i>-butane, and <i>n</i>-decane, the liquid-phase adsorption of 1,5-pentanediol and ethanol onto the MFI zeolite, and a relatively incompressible mixture of ethanol, <i>n</i>-dodecane, and water at liquid-liquid equilibrium. Property predictions are compared to those from numerical differentiation of simulation data sequentially changing the composition and from equations of state. The method can also be extended to reaction equilibria in a closed system and is applied to a reactive first-principles Monte Carlo simulation of compressed nitrogen/oxygen.