Search

First principles Monte Carlo simulations in the Gibbs and isobaric–isothermal ensembles were performed to map the vapor–liquid coexistence curves of methanol and methane described by Kohn–Sham density functional theory using the Becke–Lee–Yang–Parr (BLYP) exchange and correlation functionals with the Grimme correction term for dispersive (D2) interactions. The simulations indicate that the BLYP-D2 description with the TZV2P basis set underpredicts the saturated vapor densities and overpredicts the saturated liquid densities and critical and boiling temperatures for both compounds. Although the deviations are quite large, these results present a significant improvement over the BLYP functional without the correction term, which misses the experimental results by a larger extent in the opposite direction. Simulations at one temperature indicate that use of the larger QZV3P basis set may lead to improved saturated vapor densities, but not to significant changes in the liquid density. [ABSTRACT FROM AUTHOR]

Gibbs ensemble Monte Carlo simulations were run with an efficient mixed-basis electronic structure method to explore the phase equilibria of water from first principles using Kohn-Sham density functional theory. The Perdew-Burke-Ernzerhof exchange/correlation density functional gives a higher critical temperature (700 K) and boiling point (480 K) than experiment, although good agreement is found for the saturated liquid densities. A systematic increase in the size of the basis set for the Becke-Lee-Yang-Parr exchange/correlation density functional from a double-ζ to quadruple-ζ split valence leads to further deviations from experiment on the saturated liquid and vapor densities, while the intermediate basis set gives the best results for the heat of vaporization at T = 423 K. Analysis of the liquid structure for all simulations shows changes that can partially be explained by the different densities at a given temperature, and both density functionals show a similar temperature dependence of the liquid structure. [ABSTRACT FROM AUTHOR]

Abstract: Gibbs ensemble Monte Carlo simulations were carried out to investigate the properties of a frozen-electron-density (or Kim–Gordon, KG) model of water along the vapor–liquid coexistence curve. Because of its theoretical basis, such a KG model provides for seamless coupling to Kohn–Sham density functional theory in mixed quantum mechanics/molecular mechanics implementations. The Gibbs ensemble simulations indicate limited transferability of a KG model to other state points. Specifically, a KG model that was parameterized by Barker and Sprik to the properties of liquid water at 300K, yields saturated vapor pressures and a critical temperature that are significantly under- and overestimated, respectively. [Copyright &y& Elsevier]