Sodium chloride in supercritical water as a function of density: potentials of mean force and an equation for the dissociation constant from 723 to 1073 K and from 0 to 0.9 g/cm(3)

The potential of mean force (PMF) of sodium chloride in water has been calculated by using the ab initio classical free-energy perturbation method at five state points: at 973 K with densities of 0.2796, 0.0935, and 0.0101 g/cm (3) and at 723 K with densities of 0.0897 and 0.0098 g/cm (3). The method is based on a QM-MM model in which Na-H 2O, Cl-H 2O, and Na-Cl interactions are calculated by ab initio methods. The water-water interactions are from the polarizable TIP4P-FQ model. The logarithm of the dissociation constant (log K c) has been calculated from the PMF. These predictions, together with experimental measurements, were used to derive an equation for log K c at densities from 0 to 0.9 g/cm (3) and temperatures from 723 to 1073 K, as well as from 600 to 1073 K for densities from 0.29 g/cm (3) to 0.9 g/cm (3). Extrapolation of the present equation below 723 K for densities less than 0.29 g/cm (3) does not fit the experimental results. This is attributed to long-range changes in the local dielectric constant due to the high compressibility. Comparisons with previous predictions and simulations are presented.     

Molecular electrostatic charge models: A topographical approach

Electrostatic charge models for molecules have been developed by employing the critical topographical features of the molecular electrostatic potential (MESP ) as the “fitting” criterion. These models include one or more spherical Gaussians for incorporating the continuous electron-charge distribution in addition to the positive valued point charges representing the nuclei. The model parameters (point charges, the orbital exponents, and Gaussian centers) are optimized so as to mimic the extremal characteristics of the corresponding quantum chemical MESP . The test cases reported here include methane, ethylene, and methanol molecules. The charge models developed using the present method are seen to satisfactorily reproduce the ab initio MESP and its extremal features. © 1994 John Wiley & Sons, Inc.     

Application of Static Charge Transfer within an Ionic‐Liquid Force Field and Its Effect on Structure and Dynamics

The effects of linear scaling of the atomic charges of a reference potential on the structure, dynamics, and energetics of the ionic liquid 1,3‐dimethylimidazolium chloride are investigated. Diffusion coefficients that span over four orders of magnitude are observed between the original model and a scaled model in which the ionic charges are ±0.5 e. While the three‐dimensional structure of the liquid is less affected, the partial radial distribution functions change markedly—with the positive result that for ionic charges of ±0.7 e, an excellent agreement is observed with ab initio molecular dynamics data. Cohesive energy densities calculated from these partial‐charge models are also in better agreement with those calculated from the ab initio data. We postulate that ionic‐liquid models in which the ionic charges are assumed to be ±1 e overestimate the intermolecular attractions between ions, which results in overstructuring, slow dynamics, and increased cohesive energy densities. The use of scaled‐charge sets may be of benefit in the simulation of these systems—especially when looking at properties beyond liquid structure—thus providing an alternative to computationally expensive polarisable force fields.     

Molecular dynamics simulations of aqueous solutions of ethanolamines

We report on molecular dynamics simulations performed at constant temperature and pressure to study ethanolamines as pure components and in aqueous solutions. A new geometric integration algorithm that preserves the correct phase space volume is employed to study molecules having up to three ethanol chains. The most stable geometry, rotational barriers, and atomic charges were obtained by ab initio calculations in the gas phase. The calculated dipole moments agree well with available experimental data. The most stable conformation, due to intramolecular hydrogen bonding interactions, has a ringlike structure in one of the ethanol chains, leading to high molecular stability. All molecular dynamics simulations were performed in the liquid phase. The interaction parameters are the same for the atoms in the ethanol chains, reducing the number of variables in the potential model. Intermolecular hydrogen bonding is also analyzed, and it is shown that water associates at low water mole fractions. The force field reproduced (within 1%) the experimental liquid densities at different temperatures of pure components and aqueous solutions at 313 K. The excess and partial molar volumes are analyzed as a function of ethanolamine concentration.     

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.     

Potential energy surfaces for small alcohol dimers I: methanol and ethanol

Potential energy landscapes for homogeneous dimers of methanol and ethanol were calculated using counterpoise (CP) corrected energies at the MP26-311+G(2df,2pd) level. The landscapes were sampled at approximately 15 dimer separation distances for different relative monomer geometries, or routes, given in terms of a relative monomer yaw, pitch, and roll and the spherical angles between the monomer centers (taken as the C atom attached to the O). The 19 different routes studied for methanol and the 22 routes examined for ethanol include 607 CP corrected energies. Both landscapes can be adequately represented by site-site, pairwise-additive models, suitable for use in molecular simulations. A modified Morse potential is used for the individual pair interactions either with or without point charges to represent the monomer charge distribution. A slightly better representation of the methanol landscape is obtained using point charges, while the potential energy landscape of ethanol is slightly better without point charges. This latter representation may be computationally advantageous for molecular simulations because it avoids difficulties associated with long-range effects of point-charge-type models.     

Validation of intermolecular pair potential model of SiH4: molecular-dynamics simulation for saturated liquid density and thermal transport properties

We demonstrate a validation of the intermolecular pair potential model of SiH(4), which is constructed from ab initio molecular-orbital calculations and expressed as the sum of the exponential and the London dispersion terms. The saturated liquid densities of SiH(4) are calculated for temperatures from 100 to 225 K by molecular-dynamics (MD) simulation. The average deviation between the experiment and the MD simulation using the present potential model is 3.9%, while the deviations exceed 10% for other well-known potential models such as the five-center Lennard-Jones (LJ) model. Subsequently, the shear viscosity, the thermal conductivity, and the self-diffusion coefficient of liquid SiH(4) are calculated by an equilibrium MD simulation with the Green-Kubo formula from 100 to 225 K. The average deviations from experiment are 11.8% and 13.7% for the shear viscosity and the thermal conductivity, respectively. Comparing the present model with an empirical one-center LJ model, it turns out that the rotational energy transfer through the intermolecular potential energy, which comes from the anisotropic potential energy, plays an important role in the thermal conductivity of liquid SiH(4). These results indicate that the present intermolecular potential model has an ability to give realistic pictures for liquid SiH(4) through molecular simulations.     

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.     

The use of two-phase molecular dynamics simulations to determine the phase behavior and critical point of propane molecular models

Molecular dynamics simulations were performed to determine two-phase configurations of model propane molecules below the critical point and in the near-critical, two-phase region. A postprocessor that uses a Monte Carlo method for determination of volumes attributable to each molecule was used to obtain density histograms of the particles from which the bulk coexisting equilibrium vapor and liquid densities were determined. This method of analyzing coexisting densities in a two-phase simulation is straightforward and can be easily implemented for complex, multisite models. Various degrees of internal flexibility in the propane models have little effect on the coexisting densities at temperatures 40 K or more below the critical point, but internal flexibility (angle bending and bond vibrations) does affect the saturated liquid densities in the near-critical region, changing the critical temperature by approximately 20 K. Shorter cutoffs were also found to affect the phase dome and the location of the critical point.     

Molecular dynamics simulations of liquid methanol and methanol-water mixtures with polarizable models

A polarizable model for simulation of liquid methanol, compatible with the COS/G2 water model, has been developed using the Charge-on-Spring (COS) technique. The model consists of three point charges, with one polarizable center on the oxygen atom. The Lennard-Jones parameters on the oxygen atom together with the molecular polarizability were varied to reproduce the experimental heat of vaporization and density of liquid methanol at ambient conditions. We examined the energies of various methanol dimers in the gas phase and compared them with values obtained from ab initio calculations. The model was then used to study the thermodynamic, dynamic, structural, and dielectric properties of liquid methanol as well as of a methanol-water mixture. A microscopic picture of the structure of pure liquid methanol and of the methanol-water mixture is provided. Good agreement was found between the results from our model simulations and available experimental and ab initio calculation data. In particular, the experimental dielectric permittivity of 32 could be reproduced, which had been shown to be difficult when using nonpolarizable models.     

Low temperature second virial coefficients of H2-He obtained from quantum mechanical pair correlation functions

The previously tested ab initio potential of H₂-He has been used for calculating scattering wave functions, diagonal density matrix elements, pair correlation functions, and finally second virial coefficients, for temperatures up to 100 K. Qualitative agreement with experiment has been found at 90 K, but deviations increased with decreasing temperature up to more than a factor two at temperatures below 20 K. The most likely explanation of this strong effect is the neglect of the formation of dimer hydrogen at the very low temperatures. Some aspects concerning the reliability of the theoretical results will be discussed which are at least partly relevant for equivalent methods.     

Spectroscopic and acoustic study in binary mixtures of glycerol with several alkanols

Speeds of sound and densities of glycerol + methanol, glycerol + ethanol and glycerol + 2-propanol, were measured over the entire composition range at 298.15 K. The excess volumes, the isentropic compressibilities, molar isentropic compressibilities and excess molar isentropic compressibilities and excess speeds of sound were estimated from the densities and speeds of sound. The results indicated the presence of interactions between unlike molecules through intermolecular hydrogen bonding. The excess volumes, excess molar isentropic compressibilities and excess speeds of sound of the binary mixtures were fitted to the Redlich–Kister equation. The infrared spectra of glycerol + methanol, glycerol + ethanol and glycerol + 2-propanol have been recorded for various concentrations at room temperature. IR stretching frequencies, bandwidths and relative intensities have been estimated and analysed. Acoustic and spectroscopic measurements showed a good correlation to explain the existence of interactions between unlike molecules through intermolecular hydrogen bonding.     

Molecular dynamics study of CaCl2 in methanol

The effective potential to describe interactions between the calcium ion and methanol molecule was derived from a potential energy surface. The energy surface was constructed from more than 10,000 points obtained from ab initio calculations, performed at the MP2/6-31G(d,p) level using Gaussian98. The derived effective potential and the Palinkas-Hawlicka-Heinzinger methanol model were employed in the molecular dynamics (MD) simulation of 0.25 M CaCl(2) solution in methanol at 298 K. Average distances between ions and the methanol sites--oxygen and carbon atoms--obtained from MD simulation agreed excellently with the X-ray data. The coordination number of the cation, greater as compared with the experimental results, is in accord with its concentration dependence.     

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.     

Ab initio virial equation of state for argon using a new nonadditive three-body potential

An ab initio nonadditive three-body potential for argon has been developed using quantum-chemical calculations at the CCSD(T) and CCSDT levels of theory. Applying this potential together with a recent ab initio pair potential from the literature, the third and fourth to seventh pressure virial coefficients of argon were computed by standard numerical integration and the Mayer-sampling Monte Carlo method, respectively, for a wide temperature range. All calculated virial coefficients were fitted separately as polynomials in temperature. The results for the third virial coefficient agree with values evaluated directly from experimental data and with those computed for other nonadditive three-body potentials. We also redetermined the second and third virial coefficients from the best experimental pρT data utilizing the computed higher virial coefficients as constraints. Thus, a significantly closer agreement of the calculated third virial coefficients with the experimental data was achieved. For different orders of the virial expansion, pρT data have been calculated and compared with results from high quality measurements in the gaseous and supercritical region. The theoretically predicted pressures are within the very small experimental errors of ±0.02% for p ≤ 12 MPa in the supercritical region near room temperature, whereas for subcritical temperatures the deviations increase up to +0.3%. The computed pressure at the critical density and temperature is about 1.3% below the experimental value. At pressures between 200 MPa and 1000 MPa and at 373 K, the calculated values deviate by 1% to 9% from the experimental results.     

Vapor-liquid and vapor-solid phase equilibria for united-atom benzene models near their triple points: the importance of quadrupolar interactions

Gibbs ensemble Monte Carlo simulations were used to calculate the vapor-liquid and vapor-solid coexistence curves for benzene using two simple united-atom models. An extension of the Gibbs ensemble method that makes use of an elongated box containing a slab of the condensed phase with a vapor phase along one axis was employed for the simulations of the vapor-solid equilibria and the vapor-liquid equilibria at very low reduced temperatures. Configurational-bias and aggregation-volume-bias Monte Carlo techniques were applied to improve the sampling of particle transfers between the two simulation boxes and between the vapor and condensed-phase regions of the elongated box. An isotropic united-atom representation with six Lennard-Jones sites at the positions of the carbon atoms was used for both force fields, but one model contained three additional out-of-plane partial charge sites to explicitly represent benzene's quadrupolar interactions. Both models were fitted to reproduce the critical temperature and density of benzene and yield a fair representation of the vapor-liquid coexistence curve. In contrast, differences between the models are very large for the vapor-solid coexistence curve. In particular, the lack of explicit quadrupolar interactions for the 6-site model greatly reduces the energetic differences between liquid and solid phases, and this model yields a triple point temperature that is about a factor of 2 too low. In contrast, the 9-site model predicts a triple point of benzene at T = 253 +/- 6 K and p = 2.3 +/- 0.8 kPa in satisfactory agreement with the experimental data (T = 278.7 K and p = 4.785 kPa).     

Vapor-phase chemical equilibrium for the hydrogenation of benzene to cyclohexane from reaction-ensemble molecular simulation

The reaction-ensemble Monte-Carlo (REMC) molecular simulation method was used to study the vapor-phase chemical equilibrium for the reaction of hydrogenation of benzene to cyclohexane. A one-center modified Buckingham exponential-6 (1CMBE6) effective pair potential model (that had already been used to predict thermodynamic properties and liquid–liquid equilibria of helium+hydrogen mixtures) was used for hydrogen. Six-center modified Buckingham exponential-6 (6CMBE6) effective pair potential models (that had already been used to reproduce the saturated liquid and vapor densities, vapor pressures, second virial coefficients, and critical parameters of the six-membered ring molecules), were used for benzene and cyclohexane. No binary adjustable parameters were needed to compute the unlike-pair Buckingham exponential-6 interactions in the ternary system. Simulation results were obtained for the effect of some operating variables such as temperature (from 500 to 650 K), pressure (from 1 to 30 bar), and hydrogen to benzene feed mole ratio (from 1.5:1 to 6:1) on the reaction conversion, molar composition, and mass density of the ternary system at equilibrium. These results were found to be in excellent agreement with calculations using the predictive Soave–Redlich–Kwong (PSRK) group contribution equation of state.     

Molecular simulation of the shear viscosity and the self-diffusion coefficient of mercury along the vapor-liquid coexistence curve

In earlier work [G. Raabe and R. J. Sadus, J. Chem. Phys. 119, 6691 (2003)] we reported that the combination of an accurate two-body ab initio potential with an empirically determined multibody contribution enables the prediction of the phase coexistence properties, the heats of vaporization, and the pair distribution functions of mercury with reasonable accuracy. In this work we present molecular dynamics simulation results for the shear viscosity and self-diffusion coefficient of mercury along the vapor-liquid coexistence curve using our empirical effective potential. The comparison with experiment and calculations based on a modified Enskog theory shows that our multibody contribution yields reliable predictions of the self-diffusion coefficient at all densities. Good results are also obtained for the shear viscosity of mercury at low to moderate densities. Increasing deviations between the simulation and experimental viscosity data at high densities suggest that not only a temperature-dependent but also a density-dependent multibody contribution is necessary to account for the effect of intermolecular interactions in liquid metals. An analysis of our simulation data near the critical point yields a critical exponent of beta = 0.39, which is identical to the value obtained from the analysis of the experimental saturation densities.