Extension of the NpT + test particle method for the calculation of phase equilibria of nitrogen + ethane

Computing the starting state for Gibbs-Duhem integration

Gibbs-Duhem integration implies the numerical integration of a Clapeyron equation. To start the numerical integration, an initial coexistence point and a corresponding initial slope of the Clapeyron equation are needed. In order to apply Gibbs-Duhem integration to all kinds of systems at diverse physical conditions, one has to investigate and assess the available methods that can be used to compute these initial values. This publication focuses on vapor-liquid equilibria in binary mixtures comprising chain molecules. The initial coexistence point is either computed with the NVbeta Gibbs ensemble or with the Npbeta+test molecule method with overlapping distributions, which is introduced in this publication. Although computationally demanding, the Npbeta+test molecule method with overlapping distributions is applicable at conditions where the NVbeta Gibbs ensemble fails. We investigated three methods that can be employed to compute the initial slope of the Clapeyron equation. The Widom method and the overlapping-distributions difference method provide correct values for the initial slope. The difference method does only provide the correct answer in special cases. The possibility to judge the reliability of the results makes the overlapping-distributions difference method the safest route to the initial slope. Gibbs-Duhem integration requires the frequent computation of the slope of the Clapeyron equation. This slope depends on ensemble averages of the composition. A new bias method for efficient sampling of the composition in a semigrand-canonical simulation of chain molecules is presented. This bias method considerably enhances the composition sampling in systems comprising chain molecules of different sizes.

An advanced Gibbs-Duhem integration method: Theory and applications

The conventional Gibbs-Duhem integration method is very convenient for the prediction of phase equilibria of both pure components and mixtures. However, it turns out to be inefficient. The method requires a number of lengthy simulations to predict the state conditions at which phase coexistence occurs. This number is not known from the outset of the numerical integration process. Furthermore, the molecular configurations generated during the simulations are merely used to predict the coexistence condition and not the liquid- and vapor-phase densities and mole fractions at coexistence. In this publication, an advanced Gibbs-Duhem integration method is presented that overcomes above-mentioned disadvantage and inefficiency. The advanced method is a combination of Gibbs-Duhem integration and multiple-histogram reweighting. Application of multiple-histogram reweighting enables the substitution of the unknown number of simulations by a fixed and predetermined number. The advanced method has a retroactive nature; a current simulation improves the predictions of previously computed coexistence points as well. The advanced Gibbs-Duhem integration method has been applied for the prediction of vapor-liquid equilibria of a number of binary mixtures. The method turned out to be very convenient, much faster than the conventional method, and provided smooth simulation results. As the employed force fields perfectly predict pure-component vapor-liquid equilibria, the binary simulations were very well suitable for testing the performance of different sets of combining rules. Employing Lorentz-Hudson-McCoubrey combining rules for interactions between unlike molecules, as opposed to Lorentz-Berthelot combining rules for all interactions, considerably improved the agreement between experimental and simulated data.