Ab initio molecular dynamics simulations of the liquid/vapor interface of sulfuric acid solutions.

Ab initio molecular dynamics simulations of the liquid-vapor interface are presented for thin slabs of 72 water molecules containing a single molecule of sulfuric acid. Trajectories in the 306-330 K range are calculated for two functionals with double- and triple-ζ quality basis sets. Comparisons are made between BLYP and HCTH/120 results for the slab simulations and for bulk simulations of one H(2)SO(4) in a periodic box with 63 waters. Good agreement is found with the available experimental data and the results of other relevant AIMD studies with respect to ionization of the acid, size of the coordination shells, partitioning of the ions with the hydronium exhibiting a surface preference and the anions in the interior, and the orientational distributions for the hydronium ions and for the surface/subsurface water molecules. The major differences in the performance of the two functionals are attributable to the greater basicity of the anion oxygen atoms with the HCTH functional and the more structured aqueous solution with BLYP. The enhanced basicity results in larger aqueous coordination shells for the anion oxygens. The structuring of the BLYP aqueous solution is observed in the corrugation of the water density profile, the higher first peak in g(OO)(r), and a smaller water self-diffusion constant. This structuring with the BLYP functional yields anion hydrogen bonds that endure longer and where the dissociated ions more rapidly and directly segregate in the slab. The simulations indicate that aqueous surfaces containing ionizable diprotic acids can be modeled with rather modest sized systems and be informative.

Vibrational spectra of α-glucose, ß-glucose, and sucrose: anharmonic calculations and experiment.

The anharmonic vibrational spectra of α-D-glucose, ß-D-glucose, and sucrose are computed by the vibrational self-consistent field (VSCF) method, using potential energy surfaces from electronic structure theory, for the lowest energy conformers that correspond to the gas phase and to the crystalline phase, respectively. The results are compared with ultraviolet-infrared (UV-IR) spectra of phenyl ß-D-glucopyranoside in a molecular beam, with literature results for sugars in matrices and with new experimental data for the crystalline state. Car-Parrinello dynamics simulations are also used to study temperature effects on the spectra of α-D-glucose and ß-D-glucose and to predict their vibrational spectra at 50, 150, and 300 K. The effects of temperature on the spectral features are analyzed and compared with results of the VSCF calculations conducted at 0 K. The main results include: (i) new potential surfaces, constructed from Hartree-Fock, adjusted to fit harmonic frequencies from Møller-Plesset (MP2) calculations, that give very good agreement with gas phase, matrix, and solid state spectra; (ii) computed infrared spectra of the crystalline solid of α-glucose, which are substantially improved by including mimic groups that represent the effect of the solid environment on the sugar; and (iii) identification of a small number of combination-mode transitions, which are predicted to be strong enough for experimental observation. The results are used to assess the role of anharmonic effects in the spectra of the sugars in isolation and in the solid state and to discuss the spectroscopic accuracy of potentials from different electronic structure methods.

Sugar-salt and sugar-salt-water complexes: structure and dynamics of glucose-KNO3-(H2O)n.

Molecular dynamics (MD) simulations are carried out for the complex of glucose with KNO(3) and for complexes of the type glucose-KNO(3)-(H(2)O)(n), for n < or = 11. Structure and dynamic properties of the systems are explored. The MD simulations are carried out using primarily the DL_POLY/OPLS force field, and global and local minimum energy structures of some of the systems are compared with ab initio calculations. The main findings include: (1) complexation with KNO(3) leads to an "inverse anomeric effect", with the beta-glucose complex more stable than the alpha-glucose by approximately 1.74 kcal mol(-1); (2) as temperature is increased to 600 K, the KNO(3) remains undissociated in the 1 : 1 complex, with the K(+) hooked to the equilibrium site, and the NO(3)(-) bound to it, undergoing large-amplitude bending/torsional motions; (3) for n > or = 3 water molecules added to the system, charge separation into K(+) and NO(3)(-) ions takes place; (4) for the sugar-water system with n = 11 water molecules all hydroxyl groups are hydrated with the glucose adopting a surface position, indicative of a surfactant property of the sugar; and (5) comparison of DL_POLY with MP2/TZP structure predictions indicates that the empirical force field predicts global and local minimum structures reasonably well, but errs in giving the energy rankings of the different minima. The implications of the results on the effects of salts on saccharides are discussed.

Clathrate hydrates with hydrogen-bonding guests.

Clathrate hydrates (CHs) are inclusion compounds in which "tetrahedrally" bonded H(2)O forms a crystalline host lattice composed of a periodic array of cages. The structure is stabilized by guest particles which occupy the cages and interact with cage walls via van der Waals interactions. A host of atoms or small molecules can act as guests; here the focus is on guests that are capable of strong to intermediate H-bonding to water (small ethers, H(2)S, etc.) but nevertheless "choose" this hydrate crystal form in which H-bonding is absent from the equilibrium crystal structure. These CHs can form by exposure of ice to guest molecules at temperatures as low as 100-150 K, at the (low) guest saturation pressure. This is in contrast to the "normal" CHs whose formation typically requires temperatures well above 200 K and at least moderate pressures. The experimental part of this study addresses formation kinetics of CHs with H-bonding guests, as well as transformation kinetics between different CH forms, studied by CH infrared spectroscopy. The accompanying computational study suggests that the unique properties of this family of CHs are due to exceptional richness of the host lattice in point defects, caused by defect stabilization by H-bonding of water to the guests.

Hydroxide impurity in ice.

The role of hydroxide ions in proton mobility in water ice I is studied using both computational methods and isotopic exchange experiments. A strong influence of base adsorbate at the ice nanoparticle surface on proton activity is observed experimentally. Trace doping of the ice surface with base adsorbates stops proton activity, whereas the activity is restarted if larger amounts of adsorbate are present. On the basis of density functional theory calculations for ice slabs, these observations are rationalized by the existence of traps strongly binding hydroxide ions both at the surface and in the ice interior. These strongly binding sites stop proton mobility by immobilization of hydroxide ions; however, in the presence of an abundance of OH(-), the deepest traps are saturated and thus the proton activity is restored. A rather broad energy distribution is calculated for the hydroxide located at different lattice sites.

Hydroxide trapped in the interior of ice: a computational study.

A new configuration is proposed for trapping of the hydroxide impurity in ice and in this "off-the-lattice" configuration, OH(-) accepts four hydrogen bonds, and, since its H atom is pointing towards a cavity in the structure, it does not donate any hydrogen bond and this configuration is proposed to account for relatively low proton activity in hydroxide-rich ice systems, as observed in isotopic exchange experiments.

The single-crystal, basal face of ice I(h) investigated with sum frequency generation.

Sum frequency generation spectroscopy has been used to investigate the hydrogen-bonded region of single-crystal, hexagonal ice in the temperature range of 113-178 K. The temperature and polarization dependences of the signal are used in conjunction with a recent theoretical model to suggest an interpretation of the bluest and reddest of the hydrogen-bonded peaks. The reddest feature is associated with strong hydrogen bonding; the dynamic polarizability of this feature is primarily parallel to the surface. It is assigned to a cooperative motion among the companion to the free-OH and four-coordinate oscillators hydrogen bonded to dangling lone-pair molecules on the surface. The bluest hydrogen-bonded feature is similarly assigned to a cooperative motion of the OH stretch of dangling lone-pair molecules and of four-coordinate molecules in the lower half bilayer that are hydrogen bonded to free-OH molecules. Reconstruction induced strain is present at as low as 113 K. These results provide a richer picture of the ice surface than has heretofore been possible.

Evidence for the surface origin of point defects in ice: control of interior proton activity by adsorbates.

Spectroscopic studies are presented of H-D isotopic exchange in the interior of ice nanocrystals. The exchange process is dominated by ionic and orientational defects long viewed as governing the electrical properties of ice. A new finding that interior exchange rates can be controlled by acidic and basic adsorbates is evidence that the defects originate at the ice surface. In particular, it is argued that interior isotopic exchange is a reflection of proton concentrations equilibrated at the ice surface.

Autoionization at the surface of neat water: is the top layer pH neutral, basic, or acidic?

Autoionization of water which gives rise to its pH is one of the key properties of aqueous systems. Surfaces of water and aqueous electrolyte solutions are traditionally viewed as devoid of inorganic ions; however, recent molecular simulations and spectroscopic experiments show the presence of certain ions including hydronium in the topmost layer. This raises the question of what is the pH (defined using proton concentration in the topmost layer) of the surface of neat water. Microscopic simulations and measurements with atomistic resolution show that the water surface is acidic due to a strong propensity of hydronium (but not of hydroxide) for the surface. In contrast, macroscopic experiments, such as zeta potential and titration measurements, indicate a negatively charged water surface interpreted in terms of preferential adsorption of OH(-). Here we review recent simulations and experiments characterizing autoionization at the surface of liquid water and ice crystals in an attempt to present and discuss in detail, if not fully resolve, this controversy.

Water surface is acidic.

Water autoionization reaction 2H2O --> H3O- + OH- is a textbook process of basic importance, resulting in pH = 7 for pure water. However, pH of pure water surface is shown to be significantly lower, the reduction being caused by proton stabilization at the surface. The evidence presented here includes ab initio and classical molecular dynamics simulations of water slabs with solvated H3O+ and OH- ions, density functional studies of (H2O)(48)H+ clusters, and spectroscopic isotopic-exchange data for D2O substitutional impurities at the surface and in the interior of ice nanocrystals. Because H3O+ does, but OH- does not, display preference for surface sites, the H2O surface is predicted to be acidic with pH < 4.8. For similar reasons, the strength of some weak acids, such as carbonic acid, is expected to increase at the surface. Enhanced surface acidity can have a significant impact on aqueous surface chemistry, e.g., in the atmosphere.

HCl solvation at the surface and within methanol clusters/nanoparticles II: Evidence for molecular wires.

Condensed-phase solvation of HCl on and within methanol nanoparticles was investigated by Fourier transform infrared (FTIR) spectroscopy, on-the-fly molecular dynamics as implemented in the density functional code Quickstep (which is part of the CP2K package), and ab initio calculations. Adsorption and solvation stages are identified and assigned with the help of calculated infrared spectra obtained from the simulations. The results have been further checked with MP2-level ab initio calculations. The range of acid solvation states extends from the single-coordinated slightly stretched HCl to proton-sharing with Zundel-like methanol O...H+...X- states, and finally to MeOH2+...Cl- units with full proton transfer. Furthermore, once the proton moves to methanol, it is mobilized along methanol molecular chains. Since the proton dynamics reflects the evolving local structures, the "proton" spectra display broad bands usually with underlying continua.

Rates and mechanisms of conversion of ice nanocrystals to hydrates of HCl and HBr: acid diffusion in the ionic hydrates.

This FTIR study focuses on solid-state chemistry associated with formation and interconversion of the ionic HX (X = Cl, Br) hydrates. Kinetic data are reported for conversions of ice nanocrystal arrays exposed to the saturation pressure of the acids in the 110 approximately 125 K range. The product is amorphous acid dihydrate in the case of HBr, and amorphous monohydrate for HCl. The rate-determining step is identified as HX diffusion through the hydrate product crust toward the interfacial reaction zone, rather than diffusion through ice, as commonly believed. Slowing of the conversion process is thus observed with increasing thickness of the crust. The diffusion coefficient (D(e)) and activation energy values for HX diffusion through the hydrates were evaluated with the help of the shrinking-core model. Hydrate crystallization occurs as a separate step, upon heating above 130 K. Subsequently, rates of reversible transitions between crystal di- and monohydrates were observed upon exposure to acid vapor and acid evacuation. In conversion from di- to monohydrate, the rate slows after fast formation of several layers; subsequently, diffusion through the product crust appears to be the rate-controlling step. The activation energy for HBr diffusion through crystal dihydrate is found to be significantly higher than that for the amorphous analogue. Conjecture is offered for a molecular mechanism of HX transport through the crystal hydrate, based on (i) spectroscopic/computational evidence for the presence of molecular HX bonded to X(-) in each of the ionic hydrate phases and (ii) the relative E(a) values found for HBr and HCl diffusion. Monte Carlo modeling suggests acid transport to the reaction zone along boundaries between "nanocrystallites" generated by multiple hydrate nucleation events at the particle surfaces. The reverse conversion, of crystalline monohydrate particles to the dihydrate phase, as well as dihydrate to trihydrate, displays nearly constant rate throughout the particle conversion; suggesting desorption of HX from the particle surface as the rate-limiting factor. Like for D(e), the activation energies for desorption were found to be approximately 20% greater for HCl than HBr for related hydrate phases.

Calculation of infrared absorption spectra using Gaussian variational wave packets.

Formulas are developed for calculation of infrared absorption spectra with the help of Gaussian variational wave packets. The excitation is viewed as occurring to a linear combination of basis states corresponding to a Gaussian multiplied by Hermite polynomials of coordinates. Application of the formulas is examplified on HDO isolated in a D2O matrix. A simple approximation is proposed for the fundamental transitions.