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.     

Adsorption of carbon dioxide of 1-site and 3-site models in pillared clays: a Gibbs ensemble Monte Carlo simulation

Adsorption behavior of carbon dioxide confined in pillared clays is analyzed by using constant pressure Gibbs ensemble Monte Carlo (GEMC) method. In our simulation, 1-site and 3-site models are used to represent carbon dioxide. At the 1-site model, carbon dioxide is described as a Lennard-Jones (LJ) sphere, while at the 3-site model, carbon dioxide is modeled as a three-sites linear chain represented by EPM2 potential considering the quadrapolar effect. The potential model from Yi et al. for pillared clays is used to emphasize its quasi two-dimensional structure. Comparing the calculated results from the 1-site and the 3-site models at T=228.15 and 258.15 K, we observe that the adsorption amount from the two models is the same basically. However, the local density presents a significant difference, because the shoulder in the main peak near the wall from 3-site model can reflect the orientation of carbon dioxide. Accordingly, in the systematical investigation to explore the effect of porosity and pore width on the adsorption of carbon dioxide in pillared clays, the 3-site model was only used. We observe that for a narrow pore of H=1.02 nm, each isotherm shape displays type I curve, suggesting that it is not inflected by the porosity. However, for the larger pores of H=1.70 and H=2.38 nm, the increase of the porosity alters the shape of adsorption isotherms from a simple linear relation to the first order jump, indicating that the porosity is of very important factor to affect adsorption and phase behavior of fluids confined in pillared clays. The excess adsorptions of carbon dioxide at supercritical temperatures of T=323.15 and 348.15 K are also investigated. We find that the maximum exists for each excess isotherm, and the optimal pressure corresponding to the maximum increases with the pore width. However, the porosity has no significant effect on the optimal pressure.     

Gibbs系综Monte Carlo模拟甲烷的吸附平衡

在263、298和313 K下,对甲烷在1.91 nm的活性炭孔中的吸附平衡进行了Gibbs系综Monte Carlo(GEMC)模拟的研究.改进了GEMC方法,使之可用于模拟指定压力下的吸附平衡.通过改进的GEMC模拟,得到了在1.91 nm的活性炭中甲烷在263、298和313 K时的吸附等温线;发现263 K时的超额吸附量要大于298 K、313 K时的超额吸附量; 且不同温度下的超额吸附等温线均存在一最大超额吸附.263 K时,超额吸附量在5.0 MPa时出现最大值;而298 K、313 K时超额吸附量则在7.0 MPa时出现最大值.此工作为不同温度下天然气吸附存贮过程的开发及设计提供了依据.     

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.     

Prediction of the bubble point pressure for the binary mixture of ethanol and 1,1,1,2,3,3,3-heptafluoropropane from Gibbs ensemble Monte Carlo simulations using the TraPPE force field

Configurational-bias Monte Carlo simulations in the Gibbs ensemble using the TraPPE force field were carried out to predict the pressure–composition diagrams for the binary mixture of ethanol and 1,1,1,2,3,3,3-heptafluoropropane at 283.17 and 343.13 K. A new approach is introduced that allows one to scale predictions at one temperature based on the differences in Gibbs free energies of transfer between experiment and simulation obtained at another temperature. A detailed analysis of the molecular structure and hydrogen bonding for this fluid mixture is provided.     

Ethanoled gasoline bubble pressure determination: Experimental and Monte Carlo modeling

In this work, we present some experimental and modeling studies of ethanoled gasoline bubble pressures (ethanol + gasoline blends) at various temperatures and ethanol contents. Modelings are carried out using Monte Carlo simulations in a specific bubble-point pseudo ensemble and using the AUA4 force field. This method is first validated on the prediction of binary mixture bubble pressures (ethanol + n-hexane, ethanol + propylene, ethanol + toluene, ethanol + isooctane). It is shown that a good accuracy is reached without introducing empirical binary interaction parameter, demonstrating the predictivity of the approach. Then, simulations of ethanoled gasolines have been performed. The molecular representation of the gasoline is obtained using a lumping scheme from the detailed composition of a commercial gasoline. Simulation results are compared to experimental bubble pressures measured in this work on this commercial gasoline in which various proportions of ethanol have been added. From a qualitative point of view, the azeotropic behavior of such fuels is observed both experimentally and by simulations. From a quantitative point of view, an average deviation of 15% between experimental and simulation data is found. Such results show that Monte Carlo simulation using an accurate force field is an efficient method to predict phase equilibrium of complex mixtures such as oxygenated gasolines. This methodology can thus be seen as an efficient tool that can be used by engineers for fuel formulation or for equation of state or process model calibration.     

Calculating Gibbs free energies of transfer from Gibbs ensemble Monte Carlo simulations

The partitioning of the ternary systems n-pentane/n-heptane/(helium or argon) at ambient conditions is investigated using configurational-bias Monte Carlo simulations in the Gibbs ensemble. The results demonstrate that this approach yields very precise partition constants and free energies of transfer. Simulations are carried out to study the dependence of the n-pentane partitioning with respect to the carrier gas, the system size, and the overall solute concentrations. None of the changes of variables, within the ranges used here, has a significant effect on the alkane partitioning. However, chemical potentials calculated via Widom's ghost particle insertions show a strong number dependence for phases containing relatively few molecules of a given type. This problem originates from the fact that the chemical potential is calculated for a concentration of real particles plus one ghost particle that is systematically larger than the equilibrium concentration. A simple correction term is suggested to account for this problem. Received: 13 May 1998 / Accepted: 18 June 1998 / Published online: 4 September 1998     

Multibaric-multithermal ensemble molecular dynamics simulations

We present new generalized-ensemble molecular dynamics simulation algorithms, which we refer to as the multibaric-multithermal molecular dynamics. We describe three algorithms based on (1) the Nosé thermostat and the Andersen barostat, (2) the Nosé-Poincaré thermostat and the Andersen barostat, and (3) the Gaussian thermostat and the Andersen barostat. The multibaric-multithermal simulations perform random walks widely both in the potential-energy space and in the volume space. Therefore, one can calculate isobaric-isothermal ensemble averages in wide ranges of temperature and pressure from only one simulation run. We test the effectiveness of the multibaric-multithermal algorithm by applying it to a Lennard-Jones 12-6 potential system.     

Hybrid Monte Carlo with adaptive temperature in mixed-canonical ensemble: Efficient conformational analysis of RNA

A hybrid Monte Carlo method with adaptive temperature choice is presented that exactly generates the distribution of a mixed-canonical ensemble composed of two canonical ensembles at low and high temperature. The analysis of resulting Markov chains with the reweighting technique shows an efficient sampling of the canonical distribution at low temperature whereas the high temperature component facilitates conformational transitions, which allows shorter simulation times. The algorithm is tested by comparing analytical and numerical results for the small n-butane molecule before simulations are performed for a triribonucleotide. Sampling the complex multiminima energy landscape of this small RNA segment, we observe enforced crossing of energy barriers. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1689–1697, 1998     

Extension of the exp-6 model to the simulation of vapor-liquid equilibria of primary alcohols and their mixtures

Configurational-biased Gibbs ensemble Monte Carlo simulations were performed to obtain the phase behavior of the homologous series of primary alcohols from ethanol to 1-heptanol. Molecular interactions in these systems are modeled by a newly developed exp-6 potential in combination with a Coulombic intermolecular potential. Some of exp-6 potential parameters required to describe these alcohols were taken from the previous literature data reported for methanol and n-alkanes. The oxygen's potential parameters were optimized to fit the coexistence curve of these alcohols to the experimental data. Simulated values of saturated liquid and vapor densities, vapor pressures and critical constants of the alcohols are in good agreement with experimental data. The efficiency of the new model in the prediction of binary phase diagram of water/ethanol and n-hexane/1-propanol mixtures is also evaluated. The calculated mole fractions in the vapor and liquid phases of these binary mixtures also show satisfactory agreement with the experimental data.     

Water in nanopores. I. Coexistence curves from Gibbs ensemble Monte Carlo simulations

Coexistence curves of water in cylindrical and slitlike nanopores of different size and water-substrate interaction strength were simulated in the Gibbs ensemble. The two-phase coexistence regions cover a wide range of pore filling level and temperature, including ambient temperature. Five different kinds of two-phase coexistence are observed. A single liquid-vapor coexistence is observed in hydrophobic and moderately hydrophilic pores. Surface transitions split from the main liquid-vapor coexistence region, when the water-substrate interaction becomes comparable or stronger than the water-water pair interaction. In this case prewetting, one and two layering transitions were observed. The critical temperature of the first layering transition decreases with strengthening water-substrate interaction towards the critical temperature expected for two-dimensional systems and is not sensitive to the variation of pore size and shape. Liquid-vapor phase transition in a pore with a wall which is already covered with two water layers is most typical for hydrophilic pores. The critical temperature of this transition is very sensitive to the pore size, in contrast to the liquid-vapor critical temperature in hydrophobic pores. The observed rich phase behavior of water in pores evidences that the knowledge of coexistence curves is of crucial importance for the analysis of experimental results and a prerequiste of meaningful simulations.     

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

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

Phase behavior of C60 by computer simulation using ab‐initio interaction potential

A first‐principles intermolecular potential recently proposed by Pacheco and Ramalho [Phys Rev Lett 1997, 79, 3873–3876] has been used with the Gibbs ensemble and Gibbs–Duhem integration Monte Carlo methods to simulate the vapor–liquid and fluid–solid coexistence properties of C₆₀. The critical properties were calculated by fitting the results to the laws of rectilinear diameters and order parameter scaling. The triple‐point properties were determined from the limiting behavior of the Gibbs ensemble vapor–liquid simulations at the lowest temperature range. A stable liquid phase is predicted for temperatures between 1570±20 and 2006±27 K and densities between 0.444±0.003 and 1.05±0.01 nm^?3. The estimated critical and triple‐point pressures are, respectively, 35±6 and 5±16 bars. We show for the first time, to our knowledge, that it is possible, strictly by computer simulation, to estimate a triple point for C₆₀ in accordance with the predictions of theoretical methods and the basic concepts of thermodynamics. The liquid and fluid radial distribution functions indicate the presence of solid or glasslike features. This may support the suggestion of a more cooperative interaction of clusters in C₆₀. A comparison of our results with the data obtained by other authors is presented and discussed. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem 84: 375–387, 2001     

Monte Carlo simulations and analysis of scattering from neutral and polyelectrolyte polymer and polymer-like systems

Methods for analyzing small-angle scattering data from complex liquids with particles with many internal degrees of freedom have progressed substantially during the past years. This is mainly due to the use of Monte Carlo simulation techniques for obtaining scattering functions and the use of an approach in which the obtained scattering functions are subsequently parameterized, so that they can be used for data fitting. The present paper reviews recent applications of the Monte Carlo technique for obtaining parameterized scattering functions, and the application of these with emphasis on semiflexible polymers and wormlike micelle systems with and without electrostatic interactions.     

Conformational behavior of antineoplastic peptides Dolastatin 10 and Dolastatin 15 from Monte Carlo and molecular dynamics simulations

In this work, we report the results of a combined Monte Carlo (MC) and molecular dynamics (MD) approach applied to the conformational study of natural potent mitosis inhibitor Dolastatin 15. A GRID and CLogP analysis has been performed with the aim to provide some hints to the interpretation of the different behavior of Dolastatin 15 and Dolastatin 10, which has been examined in a former study. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Monte Carlo simulations of biomolecules: The MC module in CHARMM

We describe the implementation of a general and flexible Monte Carlo (MC) module for the program CHARMM, which is used widely for modeling biomolecular systems with empirical energy functions. Construction and use of an almost arbitrary move set with only a few commands is made possible by providing several predefined types of moves that can be combined. Sampling can be enhanced by noncanonical acceptance criteria, automatic optimization of step sizes, and energy minimization. A systematic procedure for improving MC move sets is introduced and applied to simulations of two peptides. The resulting move sets allow MC to sample the configuration spaces of these systems much more rapidly than Langevin dynamics. The rate of convergence of the difference in free energy between ethane and methanol in explicit solvent is also examined, and comparable performances are observed for MC and the Nosé-Hoover algorithm. Its ease of use combined with its sampling efficiency make the MC module in CHARMM an attractive alternative for exploring the behavior of biomolecular systems.     

Dynamics of polyethylene melts studied by Monte Carlo simulations on a high coordination lattice

A detailed comparison is made between the experiment, prior simulations by other groups, and our simulation based on a newly designed dynamic Monte Carlo algorithm, on the dynamics of polyethylene (PE) melts. The new algorithm, namely, noncross random two-bead move has been developed on a high coordination lattice (the 2nnd lattice) for studying the dynamics of realistic polymers. The chain length (molecular weight) in our simulation ranges from C40 (562 Da) to C324 (4538 Da). The effects of finite chain length have been confirmed and significant non-Gaussian statistics evidently results in nonstandard static and dynamic properties of short PE chains. The diffusion coefficients scale with molecular weight (M) to the −1.7 power for short chains and −2.2 for longer chains, which coincides very well with experimental results. No pure Rouse scaling in diffusion has been observed. The transitional molecular weight to the entanglement regime is around 1500 Da. The detailed mean square displacements of middle bead (g1) are presented for several chain lengths. The reptation-like slowdown can be clearly observed only above M ∼ 2400 Da. The slope 0.25 predicted by the theory for the intermediate regime is missing; instead a slope close to 0.4 appears, indicating that additional relaxation mechanism exists in this transitional region. The relaxation times extracted by fitting the autocorrelation function of end-to-end vectors with reptation model scale with M to 2.5 for long chains, which seemingly conflicts with the scaling of diffusion. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 2556–2571, 2006     

Dissipative particle dynamics simulations in the grand canonical ensemble: applications to polymer brushes.

We have used the dissipative particle dynamics (DPD) method in the grand canonical ensemble to study the compression of grafted polymer brushes in good solvent conditions. The force-distance profiles calculated from DPD simulations in the grand canonical ensemble are in very good agreement with the self-consistent field (SCF) theoretical models and with experimental results for two polystyrene brush layers grafted onto mica surfaces in toluene.