Partial molar properties from molecular simulation using multiple linear regression

Author/Creator ORCID





Citation of Original Publication

Josephson, Tyler R.; Singh, Ramanish; Minkara, Mona S.; Fetisov, Evgenii O.; Siepmann, J. Ilja; Partial molar properties from molecular simulation using multiple linear regression; Molecular Physics : An International Journal at the Interface Between Chemistry and Physics, Volume 117, 2019 - Issue 23-24;


This item is likely protected under Title 17 of the U.S. Copyright Law. Unless on a Creative Commons license, for uses protected by Copyright Law, contact the copyright holder or the author.
Public Domain Mark 1.0
This work was written as part of one of the author's official duties as an Employee of the United States Government and is therefore a work of the United States Government. In accordance with 17 U.S.C. 105, no copyright protection is available for such works under U.S. Law.



Partial molar volumes, energies, and enthalpies can be computed from NpT-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 V and instantaneous number of molecules Ni of each of n species for M 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 M 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 p. 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 NpT-Gibbs ensemble: a highly compressible natural gas condensate of methane, n-butane, and n-decane, the liquid-phase adsorption of 1,5-pentanediol and ethanol onto the MFI zeolite, and a relatively incompressible mixture of ethanol, n-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.