Ivan Gladich

The Ice−Vapor Interface and the Melting Point of Ice Ih for the Polarizable POL3 Water Model

We use molecular dynamics simulations to determine the melting point of ice I(h) for the polarizable POL3 water force field (Dang, L. X. J. Chem. Phys.1992, 97, 2659). Simulations are performed on a slab of ice I(h) with two free surfaces at several different temperatures. The analysis of the time evolution of the total energy in the course of the simulations at the set of temperatures yields the melting point of the POL3 model to be T(m) = 180 ± 10 K. Moreover, the results of the simulations show that the degree of hydrogen-bond disorder occurring in the bulk of POL3 ice is larger (at the corresponding degree of undercooling) than in ice modeled by nonpolarizable water models. These results demonstrate that the POL3 water force field is rather a poor model for studying ice and ice-liquid or ice-vapor interfaces. While a number of polarizable water models have been developed over the past years, little is known about their performance in simulations of supercooled water and ice. This study thus highlights the need for testing of the existing polarizable water models over a broad range of temperatures, pressures, and phases, and developing a new polarizable water force field, reliable over larger areas of the phase diagram.

Halide affinity for the water-air interface in aqueous solutions of mixtures of sodium salts

The water-air interface plays a critical role in many physical and chemical processes of the Earths atmosphere. In particular, heavy halide ions are strongly involved in processes of fundamental importance in determining the prevalence of many atmospheric components through heterogeneous reactions at the water-air interface. In this work, molecular dynamics simulations are used to study the halide enhancements at the water-air interface in the case of mixtures of Cl(-), Br(-), and I(-) ions. The results show a pattern of enhancement directly correlated to the anion polarizability. This effect is explained in terms of the charge distribution across the slab resembling an electrical double layer. As a result, the anions with higher polarizability lower the systems potential energy by enhancing their presence at the interface.

Spectroscopic properties of benzene at the air-ice interface: A combined experimental-computational approach

A combined experimental and computational approach was used to study the spectroscopic properties of benzene at the ice-air interface at 253 and 77 K in comparison with its spectroscopic behavior in aqueous solutions. Benzene-contaminated ice samples were prepared either by shock-freezing of benzene aqueous solutions or by benzene vapor-deposition on pure ice grains and examined using UV diffuse reflectance and emission spectroscopies. Neither the absorption nor excitation nor emission spectra provided unambiguous evidence of benzene associates on the ice surface even at a higher surface coverage. Only a small increase in the fluorescence intensity in the region above 290 nm found experimentally might be associated with formation of benzene excimers perturbed by the interaction with the ice surface as shown by ADC(2) excited-state calculations. The benzene associates were found by MD simulations and ground-state DFT calculations, although not in the arrangement that corresponds to the excimer structures. Our experimental results clearly demonstrated that the energy of the S0 → S1 electronic transition of benzene is not markedly affected by the phase change or the microenvironment at the ice-air interface and its absorption is limited to the wavelengths below 268 nm. Neither benzene interactions with the water molecules of ice nor the formation of dimers and microcrystals at the air-ice interface thus causes any substantial bathochromic shift in its absorption spectrum. Such a critical evaluation of the photophysical properties of organic contaminants of snow and ice is essential for predictions and modeling of chemical processes occurring in polar regions.

Hydrogen bonding and orientation effects on the accommodation of methylamine at the air-water interface

Methylamine is an abundant amine compound detected in the atmosphere which can affect the nature of atmospheric aerosol surfaces, changing their chemical and optical properties. Molecular dynamics simulation results show that methylamine accommodation on water is close to unity with the hydrophilic head group solvated in the interfacial environment and the methyl group pointing into the air phase. A detailed analysis of the hydrogen bond network indicates stronger hydrogen bonds between water and the primary amine group at the interface, suggesting that atmospheric trace gases will likely react with the methyl group instead of the solvated amine site. These findings suggest new chemical pathways for methylamine acting on atmospheric aerosols in which the methyl group is the site of orientation specific chemistry involving its conversion into a carbonyl site providing hydrophilic groups for uptake of additional water. This conversion may explain the tendency of aged organic aerosols to form cloud condensation nuclei. At the same time, formation of NH2 radical and formaldehyde is suggested to be a new source for NH2 radicals at aerosol surfaces, other than by reaction of absorbed NH3. The results have general implications for the chemistry of other amphiphilic organics, amines in particular, at the surface of atmospherically relevant aerosols.

Adsorption, Mobility, and Self-Association of Naphthalene and 1-Methylnaphthalene at the Water–Vapor Interface

The adsorption, mobility, and self-association of naphthalene (NPH) and 1-methylnaphthalene (1MN), two of the simplest polycyclic aromatic hydrocarbons (PAHs), at the surface of liquid water at 289 K were investigated using classical molecular dynamics (MD) simulations and free energy profile calculations across the water-vapor interface. Both NPH and 1MN, which exhibit a strong preference to be adsorbed at the water-vapor interface, are found to readily self-associate at the water surface, adopting mostly configurations with distinctly nonparallel arrangement of the two monomers. The additional methyl group of 1MN represents only a minor perturbation in terms of the hydration properties, interfacial orientation, and self-association with respect to NPH. Implications of the observed self-association behavior for fluorescence spectroscopy of NPH and 1MN in aqueous interfacial environment are discussed.

Negative heat capacity of small systems in the microcanonical ensemble

We propose a new mechanism to explain the origin of negative heat capacity in small systems. We show how a particular structure of the potential energy surface combined with a microcanonical treatment results in negative heat capacity. Moreover, this effect occurs under true equilibrium conditions. An energy landscape with a sudden and spatially large opening leads to a negative heat capacity. The magnitude of this effect is related to the extent of the opening and the number of particles in the system.

A surface-stabilized ozonide triggers bromide oxidation at the aqueous solution-vapour interface

Oxidation of bromide in aqueous environments initiates the formation of molecular halogen compounds, which is important for the global tropospheric ozone budget. In the aqueous bulk, oxidation of bromide by ozone involves a [Br•OOO−] complex as intermediate. Here we report liquid jet X-ray photoelectron spectroscopy measurements that provide direct experimental evidence for the ozonide and establish its propensity for the solution-vapour interface. Theoretical calculations support these findings, showing that water stabilizes the ozonide and lowers the energy of the transition state at neutral pH. Kinetic experiments confirm the dominance of the heterogeneous oxidation route established by this precursor at low, atmospherically relevant ozone concentrations. Taken together, our results provide a strong case of different reaction kinetics and mechanisms of reactions occurring at the aqueous phase-vapour interface compared with the bulk aqueous phase.Heterogeneous oxidation of bromide in atmospheric aqueous environments has long been suspected to be accelerated at the interface between aqueous solution and air. Here, the authors provide spectroscopic, kinetic and theoretical evidence for a rate limiting, surface active ozonide formed at the interface.

Interfaces Select Specific Stereochemical Conformations: The Isomerization of Glyoxal at the Liquid Water Interface

Interfacial chemistry involving glyoxal at aerosol surfaces is postulated to catalyze aerosol growth. This chemistry remains speculative due to a lack of detailed information concerning the physicochemical behavior of glyoxal at the interface of atmospheric aerosols. Here, we report results from high-level electronic structure calculations as well as both classical and Born-Oppenheimer ab initio molecular dynamics simulations of glyoxal solvation at the air/liquid water interface. When compared to the gas phase, the trans to cis isomerization of glyoxal at the liquid water interface is found to be catalyzed; additionally, the trans conformation is selectively solvated within the bulk to a greater degree than is the cis conformation. These two processes, i.e., the catalytic effect at the water interface and the differentially selective solvation, act to enhance the concentration of the cis isomer of glyoxal at the water interface. This has important consequences for the interpretation of experiments and for the modeling of glyoxal chemistry both at the interface of water clouds and at aerosols. Broader implications of this work relate to describing the role of interfaces in selecting specific stereo molecular structures at interfacial environments.

Ab initio study of the reaction of ozone with bromide ion

Surface level ozone destruction in polar environments may be initiated by oxidation of bromide ions by ozone, ultimately leading to Br2 production. Ab initio calculations are used to support the development of atmospheric chemistry models, but errors can occur in study of the bromide-ozone reaction due to inappropriate treatment of the many-electron species and the ionic nature of the reaction. In this work, a high level ab initio study is used to take into account the electronic correlation and the polarization effects. Our results show three possible pathways for the reaction. In particular, we find that this process, though endothermic on the singlet spin state surface, can be energetically feasible on the triplet surface. The triplet surface can be reached through photoexcitation of ozone or by the spin crossing of the potential energy surface. Because this process is known to occur in the dark, it may be that it occurs after intersystem crossing to a triplet surface. This paper also provides a starting point calibration for any future ab initio calculation studies of the bromide-ozone reaction, from the gas to the condensed phase.

Tuning the Stereoselectivity and Solvation Selectivity at Interfacial and Bulk Environments by Changing Solvent Polarity: Isomerization of Glyoxal in Different Solvent Environments

Conformational isomerism plays a central role in organic synthesis and biological processes; however, effective control of isomerization processes still remains challenging and elusive. Here, we propose a novel paradigm for conformational control of isomerization in the condensed phase, in which the polarity of the solvent determines the relative concentration of conformers at the interfacial and bulk regions. By the use of state-of-the-art molecular dynamics simulations of glyoxal in different solvents, we demonstrate that the isomerization process is dipole driven: the solvent favors conformational changes toward conformers having molecular dipoles that better match its polar character. Thus, the solvent polarity modulates the conformational change, stabilizing and selectively segregating in the bulk vs the interface one conformer with respect to the others. The findings in this paper have broader implications affecting systems involving compounds with conformers of different polarity. This work suggests novel mechanisms for tuning the catalytic activity of surfaces in conformationally controlled (photo)chemical reactions and for designing a new class of molecular switches that are active in different solvent environments.