Marcel D. Baer

Simulation and Theory of Ions at Atmospherically Relevant Aqueous Liquid-Air Interfaces

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.

Re-examining the properties of the aqueous vapor–liquid interface using dispersion corrected density functional theory

First-principles molecular dynamics simulations, in which the forces are computed from electronic structure calculations, have great potential to provide unique insight into structure, dynamics, electronic properties, and chemistry of interfacial systems that is not available from empirical force fields. The majority of current first-principles simulations are driven by forces derived from density functional theory with generalized gradient approximations to the exchange-correlation energy, which do not capture dispersion interactions. We have carried out first-principles molecular dynamics simulations of air-water interfaces employing a particular generalized gradient approximation to the exchange-correlation functional (BLYP), with and without empirical dispersion corrections. We assess the utility of the dispersion corrections by comparison of a variety of structural, dynamic, and thermodynamic properties of bulk and interfacial water with experimental data, as well as other first-principles and force field-based simulations.

Probing the Hydration Structure of Polarizable Halides: A Multiedge XAFS and Molecular Dynamics Study of the Iodide Anion

A comprehensive analysis of the H(2)O structure about aqueous iodide (I(-)) is reported from molecular dynamics (MD) simulation and X-ray absorption fine structure (XAFS) measurements. This study establishes the essential ingredients of an interaction potential that reproduces the experimentally determined first-solvation shell of aqueous iodide. XAFS spectra from the iodide K, L(1), and L(3) edges were corefined to establish the complete structure of the first hydration shell about aqueous iodide. Further, we have utilized molecular dynamics simulations employing both DFT (+dispersion) and empirical polarizable interaction potentials to generate an ensemble of structures that were directly compared to the XAFS data. Our results indicate that DFT-MD simulations provide a description of the molecular structure that is more consistent with the XAFS experimental data. The experimental data yield approximately 6.3 water molecules located at I-H and I-O distances of 2.65 and 3.50 Å, respectively. The differences in the two interaction potentials can be traced to the treatment of the electronic charge density in the vicinity of the iodide. The empirical polarizable interaction potential yields a significantly higher induced dipole for the aqueous iodide than the DFT study. The lower induced dipole moment from the DFT simulation produces a higher coordination number and leads to a more symmetric solvation environment than that produced by the empirical polarizable interaction potential. Furthermore, the hydrogen bonding of second-shell water with the first-shell water establishes a strong ordering of the water about the iodide surface.

Polarization- and Azimuth-Resolved Infrared Spectroscopy of Water on TiO2(110): Anisotropy and the Hydrogen-Bonding Network

We have investigated the structure and dynamics of thin water films adsorbed on TiO2(110) using infrared reflection-absorption spectroscopy (IRAS) and ab initio molecular dynamics. Infrared spectra were obtained for s- and p-polarized light with the plane of incidence parallel to the [001] and [11̅0] azimuths for water coverages ≤ 4 monolayers. The spectra indicate strong anisotropy in the water films. The vibrational densities of states predicted by the ab initio simulations for 1 and 2 monolayer coverages agree well with the observations. The results provide new insight into the structure of water on TiO2(110) and resolve a long-standing puzzle regarding the hydrogen bonding between molecules in the first and second monolayers on this surface. The results also demonstrate the capabilities of polarization- and azimuth-resolved IRAS for investigating the structure and dynamics of adsorbates on dielectric substrates.

Toward a Unified Picture of the Water Self-Ions at the Air–Water Interface: A Density Functional Theory Perspective

The propensities of the water self-ions, H3O(+) and OH(-), for the air-water interface have implications for interfacial acid-base chemistry. Despite numerous experimental and computational studies, no consensus has been reached on the question of whether or not H3O(+) and/or OH(-) prefer to be at the water surface or in the bulk. Here we report a molecular dynamics simulation study of the bulk vs interfacial behavior of H3O(+) and OH(-) that employs forces derived from density functional theory with a generalized gradient approximation exchange-correlation functional (specifically, BLYP) and empirical dispersion corrections. We computed the potential of mean force (PMF) for H3O(+) as a function of the position of the ion in the vicinity of an air-water interface. The PMF suggests that H3O(+) has equal propensity for the interface and the bulk. We compare the PMF for H3O(+) to our previously computed PMF for OH(-) adsorption, which contains a shallow minimum at the interface, and we explore how differences in solvation of each ion at the interface vs in the bulk are connected with interfacial propensity. We find that the solvation shell of H3O(+) is only slightly dependent on its position in the water slab, while OH(-) partially desolvates as it approaches the interface, and we examine how this difference in solvation behavior is manifested in the electronic structure and chemistry of the two ions.

Electrochemical Surface Potential Due to Classical Point Charge Models Drives Anion Adsorption to the Air-Water Interface.

We demonstrate that the driving forces for ion adsorption to the air-water interface for point charge models result from both cavitation and a term that is of the form of a negative electrochemical surface potential. We carefully characterize the role of the free energy due to the electrochemical surface potential computed from simple empirical models and its role in ionic adsorption within the context of dielectric continuum theory. Our research suggests that the electrochemical surface potential due to point charge models provides anions with a significant driving force for adsoprtion to the air-water interface. This is contrary to the results of ab initio simulations that indicate that the average electrostatic surface potential should favor the desorption of anions at the air-water interface. The results have profound implications for the studies of ionic distributions in the vicinity of hydrophobic surfaces and proteins.

The role of broken symmetry in solvation of a spherical cavity in classical and quantum water models

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.

Identifying Residual Structure in Intrinsically Disordered Systems: A 2D IR Spectroscopic Study of the GVGXPGVG Peptide

The peptide amide-I vibration of a proline turn encodes information on the turn structure. In this study, FTIR, two-dimensional IR spectroscopy and molecular dynamics simulations were employed to characterize the varying turn conformations that exist in the GVGX(L)PGVG family of disordered peptides. This analysis revealed that changing the size of the side chain at the X amino acid site from Gly to Ala to Val substantially alters the conformation of the peptide. To quantify this effect, proline peak shifts and intensity changes were compared to a structure-based spectroscopic model. These simulated spectra were used to assign the population of type-II β turns, bulged turns, and irregular β turns for each peptide. Of particular interest was the Val variant commonly found in the protein elastin, which contained a 25% population of irregular β turns containing two peptide hydrogen bonds to the proline C═O.

Dissociation of strong acid revisited: X-ray photoelectron spectroscopy and molecular dynamics simulations of HNO3 in water

Molecular-level insight into the dissociation of nitric acid in water is obtained from X-ray photoelectron spectroscopy and first-principles molecular dynamics (MD) simulations. Our combined studies reveal surprisingly abrupt changes in solvation configurations of undissociated nitric acid at approximately 4 M concentration. Experimentally, this is inferred from shifts of the N1s binding energy of HNO(3)(aq) as a function of concentration and is associated with variations in the local electronic structure of the nitrogen atom. It also shows up as a discontinuity in the degree of dissociation as a function of concentration, determined here from the N1s photoelectron signal intensity, which can be separately quantified for undissociated HNO(3)(aq) and dissociated NO(3)(-)(aq). Intermolecular interactions within the nitric acid solution are discussed on the basis of MD simulations, which reveal that molecular HNO(3) interacts remarkably weakly with solvating water molecules at low concentration; around 4 M there is a turnover to a more structured solvation shell, accompanied by an increase in hydrogen bonding between HNO(3) and water. We suggest that the driving force behind the more structured solvent configuration of HNO(3) is the overlap of nitric acid solvent shells that sets in around 4 M concentration.

Persistent Ion Pairing in Aqueous Hydrochloric Acid

For strong acids, like hydrochloric acid, the complete dissociation into an excess proton and conjugated base as well as the formation of independent solvated charged fragments is assumed. The existence of chloride-hydronium (Cl(-)···H3O(+)) contact ion pairs even in moderate concentration hydrochloric acid (2.5 m) demonstrates that the counterions do not behave merely as spectators. Through comparison of recent extended X-ray absorption fine structure (EXAFS) measurements to state-of-the-art density functional theory (DFT) simulations, we are able to obtain a unique view into the molecular structure of medium-to-high concentrated electrolytes. Here we report that the Cl(-)···H3O(+) contact ion pair structure persists throughout the entire concentration range studied and that these structures differ significantly from moieties studied in microsolvated hydrochloric acid gas phase clusters. Characterizing distinct populations of these ion pairs gives rise to a novel molecular level description of how to view the reaction network for acid dissociation and how it relates to our picture of acid-base equilibria.