Vapour liquid equilibria of mixtures from theNpTplus test particle method

Solubility of Carbon Dioxide in Pentaerythritol Hexanoate: Molecular Dynamics Simulation of a Refrigerant-Lubricant Oil System

We have investigated the solubility and the solvation structure between a refrigerant (carbon dioxide, CO2) and a lubricant oil (pentaerythritol hexanoate, PEC6) by molecular dynamics simulations. First, to investigate the solubility, we calculated the vapor-liquid equilibrium pressure. The chemical potential of the liquid phase and the gas phase were calculated, and the equilibrium state was obtained from the crossing point of these chemical potentials. The equilibrium pressures agreed well with experimental data over a wide range of temperatures and mole fractions of CO2. Second, the solvation structure was also investigated on a molecular scale. We found the following characteristics. First, the tails of the lubricant oil are relatively rigid inside the ester groups but flexible beyond. Second, CO2 molecules barely enter the lubricant core as delimited by the ester groups. Third, the double-bonded oxygen atoms of the ester groups are good sorption sites for CO2. Fourth, only a few CO2 molecules are attached to more than one carbonyl oxygen simultaneously. Finally, there is also significant unspecific sorption of CO2 in the alkane tail region. These results indicate that increasing the size of the rigid lubricant core would probably decrease the solubility, whereas increasing the number of polar groups would increase it.

Microscopic structure of liquid 1-1-1-2-tetrafluoroethane (R134a) from Monte Carlo simulation

1-1-1-2-tetrafluoroethane (R134a) is one of the most commonly used refrigerants. Its thermophysical properties are important for evaluating the performance of refrigeration cycles. These can be obtained via computer simulation, with an insight into the microscopic structure of the liquid, which is not accessible to experiment. In this paper, vapour-liquid equilibrium properties of R134a and its liquid microscopic structure are investigated using coupled-decoupled configurational-bias Monte Carlo simulation in the Gibbs ensemble, with a recent potential [J. Phys. Chem. B 2009, 113, 178]. We find that the simulations agree well with the experimental data, except at the vicinity of the critical region. Liquid R134a packs like liquid argon, with a coordination number in the first solvation shell of 12 at 260 K. The nearest neighbours prefer to be localized in three different spaces around the central molecule, in such a manner that the dipole moments are in a parallel alignment. Analysis of the pair interaction energy shows clear association of R134a molecules, but no evidence for C-HF type hydrogen bonding is found. The above findings should be of relevance to a broad range of fluoroalkanes.

Gibbs ensemble Monte Carlo simulations of binary mixtures of methane, difluoromethane, and carbon dioxide

Gibbs ensemble Monte Carlo simulations were used to study the vapor-liquid equilibrium of binary mixtures of carbon dioxide + methane and carbon dioxide + difluoromethane. The potential forcefields we employ are all atomistic models, and have not previously been mixed together to study the vapor-liquid equilibrium of the binary mixtures. In addition, we characterize the microscopic structure of these liquid mixtures. In carbon dioxide + methane at 230 K and 56 bar, the microscopic structure of carbon dioxide in the mixture is the same as that in the pure liquid. In carbon dioxide + difluoromethane at 283 K and 56 bar, the presence of carbon dioxide does not noticeably perturb the liquid structure of difluoromethane, but the structure of carbon dioxide is subtly changed, due to a strong interaction between it and difluoromethane. The simulations in the isobaric-isothermal (NPT) ensemble agree well with the experimental data, except at the two extreme regions of the pressure range. The good agreement of most simulated state points with experimental data encourages one to develop more accurate potentials for predicting the thermodynamic properties of these systems as well as other complicated systems, which are less amenable to measurement by experiment.

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.

Determination of surface tension in binary mixtures using transition-matrix Monte Carlo

We present a methodology based on grand-canonical transition-matrix Monte Carlo and finite-size scaling analysis to calculate surface tensions in binary mixtures. In particular, mixture transition-matrix Monte Carlo is first used to calculate apparent, system-size-dependent free-energy barriers separating coexisting fluid phases. Finite-size scaling is then used to extrapolate these values to the infinitely large system limit to determine the true thermodynamic surface tension. A key distinction of the methodology is that it yields the entire isothermal surface-tension curve for a binary mixture in a relatively small number of simulations. We demonstrate the utility of the method by calculating surface-tension curves for three binary Lennard-Jones mixtures. While we have only examined the surface tension of simple fluids in this work, the method is general and can be extended to molecular fluids as well as to determine interfacial tensions of liquid-liquid interfaces.

Vapour-liquid equilibrium of the OPLS, optimized potentials for liquid simulations, model for binary systems of alkanes and alkanes + alcohols

The NpT - Gibbs ensemble Monte Carlo computer simulation method was applied to predict the vapour-liquid equilibrium (VLE) behavior of the binary systems ethane + pentane at 277.55 K and 310.95 K, ethane + hexane at 298.15 K, propane + methanol at 313.15 K and propane + ethanol at 325.15 K and 425.15 K. The optimized potentials for the liquid simulating (OPLS) model were used to describe the interactions of alkanes and alcohols. The simulated VLE predictions are compared with experimental data available for the pressure and phase composition of the analyzed binary systems. The agreement between the experimental data and the simulation results is found to be generally good, although slightly better for system in which both components were nonpolar.