Direct determination of phase behavior of square-well fluids

We have combined Gibbs ensemble Monte Carlo simulations with the aggregation volume-biased method in conjunction with the Gibbs-Duhem method to provide the first direct estimates for the vapor-solid, vapor-liquid, and liquid-solid phase coexistences of square-well fluids with three different ranges of attraction. Our results are consistent with the previous simulations and verify the notion that the vapor-liquid coexistence behavior becomes metastable for cases where the attraction well becomes smaller than 1.25 times the particle diameter. In these cases no triple point is found.     

Efficient prediction of thermodynamic properties of quadrupolar fluids from simulation of a coarse-grained model: the case of carbon dioxide

Monte Carlo simulations are presented for a coarse-grained model of real quadrupolar fluids. Molecules are represented by particles interacting with Lennard-Jones forces plus the thermally averaged quadrupole-quadrupole interaction. The properties discussed include the vapor-liquid coexistence curve, the vapor pressure along coexistence, and the surface tension. The full isotherms are also accessible over a wide range of temperatures and densities. It is shown that the critical parameters (critical temperature, density, and pressure) depend almost linearly on a quadrupolar parameter q=Q(*4)T*, where Q* is the reduced quadrupole moment of the molecule and T* the reduced temperature. The model can be applied to a variety of small quadrupolar molecules. We focus on carbon dioxide as a test case, but consider nitrogen and benzene, too. Experimental critical temperature, density, and quadrupolar moment are sufficient to fix the parameters of the model. The resulting agreement with experiments is excellent and marks a significant improvement over approaches which neglect quadrupolar effects. The same coarse-grained model was also applied in the framework of perturbation theory in the mean spherical approximation. As expected, the latter deviates from the Monte Carlo results in the critical region, but is reasonably accurate at lower temperatures.     

Transferable potentials for phase equilibria. 9. Explicit hydrogen description of benzene and five-membered and six-membered heterocyclic aromatic compounds

The explicit hydrogen version of the transferable potentials for phase equilibria (TraPPE-EH) force field is extended to benzene, pyridine, pyrimidine, pyrazine, pyridazine, thiophene, furan, pyrrole, thiazole, oxazole, isoxazole, imidazole, and pyrazole. While the Lennard-Jones parameters for carbon, hydrogen (two types), nitrogen (two types), oxygen, and sulfur are transferable for all 13 compounds, the partial charges are specific for each compound. The benzene dimer energies for sandwich, T-shape, and parallel-displaced configurations obtained for the TraPPE-EH force field compare favorably with high-level electronic structure calculations. Gibbs ensemble Monte Carlo simulations were carried out to compute the single-component vapor-liquid equilibria for benzene, pyridine, three diazenes, and eight five-membered heterocycles. The agreement with experimental data is excellent with the liquid densities and vapor pressures reproduced within 1 and 5%, respectively. The critical temperatures and normal boiling points are predicted with mean deviations of 0.8 and 1.6%, respectively.     

Phase equilibria in carbon dioxide expanded solvents: Experiments and molecular simulations

We present complementary molecular simulations and experimental results of phase equilibria for carbon dioxide expanded acetonitrile, methanol, ethanol, acetone, acetic acid, toluene, and 1-octene. The volume expansion measurements were done using a high-pressure Jerguson view cell. Molecular simulations were performed using the Gibbs ensemble Monte Carlo method. Calculations in the canonical ensemble (NVT) were performed to determine the coexistence curve of the pure solvent systems. Binary mixtures were simulated in the isobaric-isothermal distribution (NPT). Predictions of vapor-liquid equilibria of the pure components agree well with experimental data. The simulations accurately reproduced experimental data on saturated liquid and vapor densities for carbon dioxide, methanol, ethanol, acetone, acetic acid, toluene, and 1-octene. In all carbon dioxide expanded liquids (CXL's) studied, the molecular simulation results for the volume expansion of these binary mixtures were found to be as good as, and in many cases superior to, predictions based on the Peng-Robinson equation of state, demonstrating the utility of molecular simulation in the prediction of CXL phase equilibria.     

Liquid-vapor and liquid-liquid phase equilibria of the Brodholt-Sampoli-Vallauri polarizable water model

Liquid-vapor and liquid-liquid phase equilibria of the polarizable Brodholt-Sampoli-Vallauri water model have been investigated by Gibbs ensemble Monte Carlo computer simulations. The coexisting liquid and vapor densities and energy of vaporization of the model is found to be in a reasonable agreement with experimental data in the entire temperature range of liquid-vapor coexistence. The critical temperature and density of the model are found to be 615 K and 0.278 gcm(3), respectively, close to the experimental values of 647.1 K and 0.322 gcm(3). In the supercooled state two distinct liquid-liquid coexistence regions are observed. The existence of liquid-liquid phase separation of a polarizable water model is demonstrated for the first time.     

Phase behavior of elemental aluminum using monte carlo simulations

Monte Carlo simulations are presented for two models of aluminum: an embedded-atom model and an explicit many-body model. Vapor/liquid coexistence curves are determined using Gibbs ensemble Monte Carlo simulations. The normal boiling points predicted by both models are somewhat higher (by about 10%) than the experimental value. Isothermal constant-stress simulations are used to simulate solid Al from 300 K to the triple point. The solid structures are at least metastable in the face-centered cubic configuration, and the specific heat is determined to be lower than the experimental value. The melting point for the embedded-atom model determined via thermodynamic integration along a pseudo-supercritical path is approximately 20% higher than the experimental value.     

Vapor-liquid equilibrium properties for confined binary mixtures involving CO2, CH4, and N2 from Gibbs ensemble Monte Carlo simulations

The effects of solid-fluid interactions on the vapor-liquid phase diagram,coexistence density,relative volatility and vaporization enthalpy have been investigated for confined binary systems of CO 2-CH 4,CO 2-N 2 and CH 4-N 2.The Gibbs ensemble Monte Carlo(GEMC) simulation results indicate that the confinement and the solid-fluid interaction have significant influences on the vapor-liquid equilibrium properties.The confinement and the strength of the solid-fluid interaction make the p-x i phase diagram move to higher pressure regions.They also make the two-phase region become narrower for each binary mixture.The strength of the solid-fluid interactions can cause increases in the coexistence liquid and vapor densities,and cause the decrease of the relative volatility and the vaporization enthalpy for the systems studied.As the pore width is decreased,the two-phase region of the binary mixture becomes narrower.     

Transferable potentials for phase equilibria. 7. Primary, secondary, and tertiary amines, nitroalkanes and nitrobenzene, nitriles, amides, pyridine, and pyrimidine

The transferable potentials for phase equilibria (TraPPE) force fields are extended to amine, nitro, nitrile, and amide functionalities and to pyridine and pyrimidine. In many cases, the same parameters for a functional group are used for both united-atom and explicit-hydrogen representations of alkyl tails. Following the TraPPE philosophy, the nonbonded interaction parameters were fitted to the vapor-liquid coexistence curves for selected one-component systems. Coupled-decoupled configurational-bias Monte Carlo simulations in the Gibbs ensemble were applied to neat (methyl-, dimethyl-, trimethyl-, ethyl-, diethyl-, or triethyl-)amine, nitromethane, nitroethane, nitrobenzene, acetonitrile, propionitrile, acetamide, propanamide, butanamide, pyridine, and pyrimidine. Excellent agreement with experimental results was found, with the mean unsigned errors being less than 1% for both the critical temperature and the normal boiling temperature. Similarly, the liquid densities at low reduced temperatures are reproduced to within 1%, and the deviation for the critical densities is about 4%. Additional simulations were performed for the binary mixtures of methylamine + n-hexane, diethyl ether + acetonitrile, 1-propanol + acetonitrile, and nitroethane + ethanol. With the exception of the methylamine/n-hexane mixture for which the separation factor is substantially overestimated, agreement with experiment for the other three mixtures is very satisfactory.     

GEMC和GDI方法计算流体气液相平衡的比较

考察采用TraPPE联合原子和OPLS全原子力场两种分子力场, Gibbs系综蒙特卡罗(GEMC)方法和Gibbs-Duhem积分(GDI)方法计算流体气液相平衡的适用性、计算速度、计算精度等问题. 结果表明, 在采用全原子力场情况下, GDI方法比GEMC方法极大地节省了计算时间. 从计算结果来看, 两种方法各有适用范围, 在使用时可互为补充. 在给定力场的前提下, 两种方法所得到的液相密度、蒸发焓、临界温度和临界密度相差不大, 而当力场中的缺陷导致蒸发焓的计算不够准确时, 两种计算方法得到的气体的压力和密度明显不同,进而导致预测的临界压力也明显不同.     

A general perturbation approach for equation of state development: applications to simple fluids, ab initio potentials, and fullerenes

A new perturbation scheme based on the Barker-Henderson perturbation theory [J. Chem. Phys. 47, 4714 (1967)] is proposed to predict the thermodynamic properties of spherical molecules. Accurate predictions of second virial coefficients and vapor-liquid coexistence properties are obtained for a large variety of potential functions (square well, Yukawa, Sutherland, Lennard-Jones, Buckingham, Girifalco). New Gibbs ensemble Monte Carlo simulations of the generalized exp-m Buckingham potential are reported. An extension of the perturbation approach to mixtures is proposed, and excellent predictions of vapor-liquid equilibria are obtained for Lennard-Jones mixtures. The perturbation scheme can be applied to complex potential functions fitted to ab initio data to predict the properties of real molecules such as neon. The new approach can also be used as an auxiliary tool in molecular simulation studies, to efficiently optimize an intermolecular potential on macroscopic properties or match force fields based on different potential functions.     

Anisotropic united-atoms (AUA) potential for alcohols

An anisotropic united-atom (AUA4) intermolecular potential has been derived for the family of alkanols by first optimizing a set of charges to reproduce the electrostatic potential of the isolated molecules of methanol and ethanol and then by adjusting the parameters of the OH group to fit selected equilibrium properties. In particular, the proposed potential includes additional extra-atomic charges in order to improve the matching to the electrostatic field. Gibbs ensemble Monte Carlo simulations were performed to determine the phase equilibria, while the critical region was explored by means of grand canonical Monte Carlo simulations combined with histogram reweighting techniques. In order to increase the transferability of the model, only the parameters of the Lennard-Jones OH group have been fitted, the parameters of the other AUA groups are taken from previous works. Nevertheless, a good level of agreement was obtained for all compounds considered in this work. In particular, excellent results were obtained for the Henry constants calculation of different gases in alkanols.     

Improvement in molecule exchange efficiency in Gibbs ensemble Monte Carlo: development and implementation of the continuous fractional component move

The continuous fractional component Monte Carlo (CFC MC) move (J Chem Theory Comput, 2007, 3, 1451) is extended to the Gibbs ensemble. The algorithm is validated against conventional simulations for the Lennard Jones fluid and a flexible water model. The method is also used to compute the vapor-liquid coexistence densities of a model for SO(2). The CFC molecule exchange move relies on the gradual insertion and deletion of molecules in conjunction with a self-adapting bias potential. As a result, the method does not require the formation of spontaneous voids in the dense fluid phase to be successful, leading to molecule exchange acceptance probabilities that are nearly independent of temperature. For example, over 1% of the vapor-liquid molecule exchange moves are successful for water at 280 K, whereas advanced rotational and configurational bias methods have success rates of less than 0.03%. The CFC move can be combined with other Monte Carlo moves to enable efficient simulation of dense strongly associating fluids that are to this point infeasible to model with standard methods.     

Transferable potentials for phase equilibria. 8. United-atom description for thiols, sulfides, disulfides, and thiophene

An extension of the transferable potentials for phase equilibria united-atom (TraPPE-UA) force field to thiol, sulfide, and disulfide functionalities and thiophene is presented. In the TraPPE-UA force field, nonbonded interactions are governed by a Lennard-Jones plus fixed point charge functional form. Partial charges are determined through a CHELPG analysis of electrostatic potential energy surfaces derived from ab initio calculations at the HF/6-31g+(d,p) level. The Lennard-Jones well depth and size parameters for four new interaction sites, S (thiols), S(sulfides), S(disulfides), and S(thiophene), were determined by fitting simulation data to pure-component vapor-equilibrium data for methanethiol, dimethyl sulfide, dimethyl disulfide, and thiophene, respectively. Configurational-bias Monte Carlo simulations in the grand canonical ensemble combined with histogram-reweighting methods were used to calculate the vapor-liquid coexistence curves for methanethiol, ethanethiol, 2-methyl-1-propanethiol, 2-methyl-2-propanethiol, 2-butanethiol, pentanethiol, octanethiol, dimethyl sulfide, diethyl sulfide, ethylmethyl sulfide, dimethyl disulfide, diethyl disulfide, and thiophene. Excellent agreement with experiment is achieved, with unsigned errors of less than 1% for saturated liquid densities and less than 3% for critical temperatures. The normal boiling points were predicted to within 1% of experiment in most cases, although for certain molecules (pentanethiol) deviations as large as 5% were found. Additional calculations were performed to determine the pressure-composition behavior of ethanethiol+n-butane at 373.15 K and the temperature-composition behavior of 1-propanethiol+n-hexane at 1.01 bar. In each case, a good reproduction of experimental vapor-liquid equilibrium separation factors is achieved; both of the coexistence curves are somewhat shifted because of overprediction of the pure-component vapor pressures.     

Critical properties of aluminum

Gibbs ensemble Monte Carlo calculations are performed using a validated embedded-atom potential to obtain the vapor-liquid coexistence curve for elemental aluminum in good agreement with available experimental data up to the boiling point. These calculations are then extended to make a reliable prediction of the critical temperature, pressure, and density of Al, which have previously been known only with very large uncertainties. This demonstrates the ability of modern simulations to predict fundamental physical properties that are extremely difficult to measure directly.     

受限狭缝中CO_2-CH_4体系的汽液相平衡性质

受限条件下CO2-CH4体系的相平衡性质对化工工艺过程的设计具有非常重要的意义.采用Gibbs系综Monte Carlo模拟,对220K下CO2-CH4体系在主体相和受限狭缝中的相平衡性质进行了系统地研究.通过主体相模拟与实验结果比较,验证了流体分子势能参数的合理性;通过改变狭缝壁面原子的能量参数,研究了受限环境对CO2-CH4体系汽液相平衡性质的影响.与主体相相比,在硬壁狭缝中,CO2-CH4体系的露点压力增加,泡点压力降低,压力-组成相图变窄,且体系更容易达到超临界状态;在吸引狭缝中,随壁面原子能量参数的增大,CO2-CH4体系的压力-组成相图上移,临界点处CH4的摩尔分数减小,相图变窄.在体系汽液相总组成相同情况下,硬壁狭缝内体系的汽液相密度均比主体相中小;随壁面原子能量参数增大,气相密度变大、液相密度在CH4的摩尔分数较小时变大而当CH4的摩尔分数达到一定值后反而减小.在体系汽液相总组成相同时,受限环境下的汽化热比主体相的汽化热小且随壁面吸引势的增强越来越小;在主体相和硬壁狭缝中体系的汽化热随CH4含量的增加单调减小,而当壁面势能参数较大时汽化热随CH4含量增加先增大后减小.     

Influence of bond flexibility on the vapor-liquid phase equilibria of water

The authors performed Gibbs ensemble simulations on the vapor-liquid equilibrium of water to investigate the influence of incorporating intramolecular degrees of freedom in the simple point charge (SPC) water model. Results for vapor pressures, saturation densities, heats of vaporization, and the critical point for two different flexible models are compared with data for the corresponding rigid SPC and SPC/E models. They found that the introduction of internal vibrations, and also their parametrization, has an observable effect on the prediction of the vapor-liquid coexistence curve. The flexible SPC/Fw model, although optimized to describe bulk diffusion and dielectric constants at ambient conditions, gives the best prediction of saturation densities and the critical point of the examined models.     

Computational studies on thermodynamic properties, effective diameters, and free volume of argon using an ab initio potential

A quantum mechanical derived ab initio interaction potential for the argon dimer was tested in molecular simulations to reproduce the thermophysical properties of the vapor-liquid phase equilibria using the Gibbs ensemble Monte Carlo simulations as well as the liquid and supercritical equation of state using the NVT Monte Carlo simulations. The ab initio interaction potential was taken from the literature. A recently developed theory [R. Laghaei et al., J. Chem. Phys. 124, 154502 (2006)] was used to compute the effective diameters of argon in fluid phases and the results were subsequently applied in the generic van der Waals theory to compute the free volume of argon. The calculated densities of the coexisting phases, the vapor pressure, and the equation of state show excellent agreement with experimental values. The effective diameters and free volumes of argon are given over a wide range of densities and temperatures. An empirical formula was used to fit the effective diameters as a function of density and temperature. The computed free volume will be used in future investigations to calculate the transport properties of argon.     

HCFC-R22的等效Stockmayer势能模型及汽液平衡分子模拟

采用等效Stockmayer势能模型得到R22的势能参数,并应用Gibbs系综模拟汽液平衡予以考察。模拟中,采用了Metropolis抽样及周期边界条件。模拟结果表明饱和液体密度精度较高,且饱和蒸汽的密度及蒸汽压的结果也十分令人满意,说明等效Stockmayer势能模型方法的可靠性和可行性。