We present an ab initio molecular dynamics simulation of the aqueous liquid-vapor interface. Having successfully stabilized a region of bulk water in the center of a water slab, we were able to reproduce and further quantify the experimentally observed abundance of surface "acceptor-only"(19%) and "single-donor"(66%) moieties as well as substantial surface relaxation approaching the liquid-vapor interface. Examination of the orientational dynamics points to a faster relaxation in the interfacial region. Furthermore, the average value of the dipole decreases and the average value of the highest occupied molecular orbital for each water molecule increases approaching the liquid-vapor interface. Our results support the idea that the surface contains, on average, far more reactive states than the bulk.

We present herein a comprehensive density functional theory study toward assessing the accuracy of two popular gradient-corrected exchange correlation functionals on the structure and density of liquid water at near ambient conditions in the isobaric-isothermal ensemble. Our results indicate that both PBE and BLYP functionals under predict the density and over structure the liquid. Adding the dispersion correction due to Grimme (1, 2) improves the predicted densities for both BLYP and PBE in a significant manner. Moreover, the addition of the dispersion correction for BLYP yields an oxygen-oxygen radial distribution function in excellent agreement with experiment. Thus, we conclude that one can obtain a very satisfactory model for water using BLYP and a correction for dispersion.

A series of first principles molecular dynamics and Monte Carlo simulations were carried out for liquid water to investigate the reproducibility of different sampling approaches. These simulations include Car−Parrinello molecular dynamics simulations using the program cpmd with different values of the fictitious electron mass in the microcanonical and canonical ensembles, Born−Oppenheimer molecular dynamics using the programs cpmd and cp2k in the microcanonical ensemble, and Metropolis Monte Carlo using cp2k in the canonical ensemble. With the exception of one simulation for 128 water molecules, all other simulations were carried out for systems consisting of 64 molecules. Although the simulations yield somewhat fortuitous agreement in structural properties, analysis of other properties demonstrate that one should exercise caution when assuming the reproducibility of Car−Parrinello and Born−Oppenheimer molecular dynamics simulations for small system sizes in the microcanonical ensemble. In contrast, the molecular dynamics and Monte Carlo simulations in the canonical ensemble appear to be more reliable. Furthermore, in the case of canonical Car−Parrinello molecular dynamics simulations the application of Nosé−Hoover chain thermostats allows the use of larger fictitious electron masses. For the Becke−Lee−Yang−Parr exchange and correlation energy functionals and norm-conserving Troullier−Martins or Goedecker−Teter−Hutter pseudopotentials, these simulations at a fixed density of 1.0 g/cm3 and a temperature close to 315 K point to an overstructured liquid with a height of the first peak in the oxygen−oxygen radial distribution function of about 3.0, an underestimated value of the classical constant-volume heat capacity of about 70 J/(mol K), and an underestimated self-diffusion constant of about 0.04 Å2/ps.

The structure, dynamical and electronic properties of liquid water utilizing different hybrid density functionals were tested within the plane wave framework of first principles molecular dynamics simulations. The computational approach, which employs modified functionals with shortranged Hartree-Fock exchange, was first tested in calculations of the structural and bonding properties of the water dimer and cyclic water trimer.Liquid water simulations were performed at the state point of 350 K at the experimental density. Simulations included three different hybrid functionals, a meta functional, four gradient corrected functionals, the local density and Hartree-Fock approximation. It is found that hybrid functionals are superior in reproducing the experimental structure and dynamical * corresponding author 1 properties as measured by the radial distribution function and self diffusion constant when compared to the pure density functionals. The local density and Hartree-Fock approximations show strongly over-and understructured liquids, respectively. Hydrogen bond analysis shows that the hybrid functionals give slightly smaller averaged numbers of hydrogen bonds and similar hydrogen bond populations as pure density functionals. The average molecular dipole moments in the liquid from the three hybrid functionals are lower than from the corresponding pure density functionals.

We have resolved the inconsistency in quantifying the surface potential at the liquid-vapor interface when using explicit ab initio electronic charge density and effective atomic partial charge models of liquid water. This is related, in part, to the fact that the resulting electric potentials from partial-charge models and ab initio charge distributions are quite different except for those regions of space between the molecules. We show that the electrostatic surface potential from a quantum mechanical charge distribution compares well to high-energy electron diffraction and electron holography measurements, as opposed to the comparison with electrochemical measurements. We suggest that certain regions of space be excluded when comparing computed surface potentials with electrochemical measurements. This work describes a novel interpretation of ab initio computed surface potentials through high-energy electron holography measurements as useful benchmarks toward a better understanding of electrochemistry.

Chemistry occurring at or near the surface of aqueous droplets and thin films in the atmosphere influences air quality and climate. Molecular dynamics simulations are becoming increasingly useful for gaining atomic-scale insight into the structure and reactivity of aqueous interfaces in the atmosphere. Here we review simulation studies of atmospherically relevant aqueous liquid-air interfaces, with an emphasis on ions that play important roles in the chemistry of atmospheric aerosols. In addition to surveying results from simulation studies, we discuss challenges to the refinement and experimental validation of the methodology for simulating ion adsorption to the air-water interface and recent advances in elucidating the driving forces for adsorption. We also review the recent development of a dielectric continuum theory capable of reproducing simulation and experimental data on ion behavior at aqueous interfaces.

Insertion of a hard sphere cavity in liquid water breaks translational symmetry and generates an electrostatic potential difference between the region near the cavity and the bulk. Here, we clarify the physical interpretation of this potential and its calculation. We also show that the electrostatic potential in the center of small, medium, and large cavities depends very sensitively on the form of the assumed molecular interactions for different classical simple point-charge models and quantum mechanical DFT-based interaction potentials, as reflected in their description of donor and acceptor hydrogen bonds near the cavity. These differences can significantly affect the magnitude of the scalar electrostatic potential. We argue that the result of these studies will have direct consequences toward our understanding of the thermodynamics of ion solvation through the cavity charging process.

This research addresses a comprehensive particle-based simulation study of the structural, dynamic, and electronic properties of the liquid-vapor interface of water utilizing both ab initio (based on density functional theory) and empirical (fixed charge and polarizable) models. Numerous properties such as interfacial width, hydrogen bond populations, dipole moments, and correlation times will be characterized with identical schemes to draw useful conclusions on the strengths and weakness of the proposed models for interfacial water. Our findings indicate that all models considered in this study yield similar results for the radial distribution functions, hydrogen bond populations, and orientational relaxation times. Significant differences in the models appear when examining both the dipole moments and surface relaxation near the aqueous liquid-vapor interface. Here, the ab initio interaction potential predicts a significant decrease in the molecular dipole moment and expansion in the oxygen-oxygen distance as one approaches the interface in accordance with recent experiments. All classical polarizable interaction potentials show a less dramatic drop in the molecular dipole moment, and all empirical interaction potentials studied yield an oxygen-oxygen contraction as the interface is approached.