Adsorption and Solvation of HCl into Ice Surfaces

Some fundamental properties and reactions of ice surfaces at low temperatures

Ice surfaces offer a unique chemical environment in which reactions occur quite differently from those in liquid water or gas phases. In this article, we examine the basic properties of ice surfaces below the surface premelting temperature and discuss some of the recent investigations carried out on reactions at the ice surfaces. The static and dynamic properties of an ice surface as a reaction medium, such as its structure, molecule diffusion and proton transfer dynamics, and the surface preference of hydronium and hydroxide ions, are discussed in relation to the reactivity of the surface.

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.

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.

H/D isotopic exchange between water molecules at ice surfaces

H/D isotopic exchange between H(2)O and D(2)O molecules was studied at the surface of ice films at 90-140 K by the technique of Cs(+) reactive ion scattering. Ice films were deposited on a Ru(0001) substrate in different compositions of H(2)O and D(2)O and in various structures to study the kinetics of isotopic exchange. H/D exchange was very slow on an ice film at 95-100 K, even when H(2)O and D(2)O were uniformly mixed in the film. At 140 K, H/D exchange occurred in a time scale of several minutes on the uniform mixture film. Kinetic measurement gave the rate coefficient for the exchange reaction, k(140 K)=1.6(+/-0.3) x 10(-19) cm(2) molecule(-1) s(-1) and k(100 K)< or =5.7(+/-0.5) x 10(-21) cm(2) molecule(-1) s(-1) and the Arrhenius activation energy, E(a)> or =9.8 kJ mol(-1). Addition of HCl on the film to provide excess protons greatly accelerated the isotopic exchange reaction such that it went to completion very quickly at the surface. The rapid reaction, however, was confined within the first bilayer (BL) of the surface and did not readily propagate to the underlying sublayer. The isotopic exchange in the vertical direction was almost completely blocked at 95 K, and it slowly occurred only to a depth of 3 BLs from the surface at 140 K. Thus, the proton transfer was highly directional. The lateral proton transfer at the surface was attributed to the increased mobility of protonic defects at the molecularly disordered and activated surface. The slow, vertical proton transfer was probably assisted by self-diffusion of water molecules.

Diffusion of HDO in Pure and Acid-Doped Ice Films

In these experiments, a few bilayers of D(2)O were vapor-deposited on a pure crystalline H(2)O ice film or an ice film doped with a small amount of HCl. Upon deposition, H/D isotopic exchange quickly converted the D(2)O layer into an HDO-rich mixture layer. Infrared absorption spectroscopy followed the changes of the HDO from the initial HDO mixture layer to HDO isolated in the H(2)O ice film. This was possible because isolated HDO in H(2)O ice has a unique, sharp peak in the O-D stretch region that can be distinguished from the broad peak due to the initial HDO mixture layer. The absorbance of isolated HDO displayed first-order kinetics and was attributed to diffusion of HDO from the HDO-rich mixture layer into the underlying H(2)O ice film. While negligible diffusion was observed for pure ice films and for ice films with HCl concentrations up to 1 x 10(-4) mole fraction, diffusion of HDO occurred for higher concentrations of (2-20) x 10(-4) mole fraction HCl with a concentration-independent rate constant. The diffusion under these conditions followed Arrhenius behavior for T = 135-145 K yielding E(a) = 25 +/- 5 kJ/mol. The mechanism for the HDO diffusion involves either (i) molecular self-diffusion or (ii) long-range H/D diffusion by a series of multiple proton hop and orientational turn steps. While these spectroscopic results compare favorably with recent studies of molecular self-diffusion in low-temperature ice films, the diffusion results from all the ice film studies at low temperatures (ca. T < 170 K) differ from earlier bulk ice studies at higher temperatures (ca. T > 220 K). A comparison and discussion of the various diffusion studies are included in this report.

Fast thermal desorption spectroscopy study of morphology and vaporization kinetics of polycrystalline ice films

Fast thermal desorption spectroscopy was used to investigate the vaporization kinetics of thin (50-100 nm) H(2)O(18) and HDO tracer layers from 2-5 microm thick polycrystalline H(2)O(16) ice films at temperatures ranging from -15 to -2 degrees C. The isothermal desorption spectra of tracer species demonstrate two distinct peaks, alpha and beta, which we attribute to the vaporization of H(2)O(18) initially trapped at or near the grain boundaries and in the crystallites of the polycrystalline ice, respectively. We show that the diffusive transport of the H(2)O(18) and HDO tracer molecules in the bulk of the H(2)O(16) film is slow as compared to the film vaporization. Thus, the two peaks in the isothermal spectra are due to unequal vaporization rates of H(2)O(18) from grain boundary grooves and from the crystallites and, therefore, can be used to determine independently the vaporization rate of the single crystal part of the film and rate of thermal etching of the film. Our analysis of the tracer vaporization kinetics demonstrates that the vaporization coefficient of single crystal ice is significantly greater than those predicted by the classical vaporization mechanism at temperatures near ice melting point. We discuss surface morphological dynamics and the bulk transport phenomena in single crystal and polycrystalline ice near 0 degrees C.

Adsorption, Ionization, and Migration of Hydrogen Chloride on Ice Films at Temperatures between 100 and 140 K

Adsorption of hydrogen chloride (HCl) on water ice films is studied in the temperature range of 100-140 K by using Cs+ reactive ion scattering (Cs+ RIS), low energy sputtering (LES), and temperature-programmed-desorption mass spectrometry (TPDMS). At 100 K, HCl on ice partially dissociates to hydronium and chloride ions and the undissociated HCl exists in two distinct molecular states (alpha- and beta-states). Upon heating of the ice films, HCl molecules in the alpha-state desorb at 135-150 K, whereas those in the beta-state first become ionized and then desorb via recombinative reaction of ions at 170 K. An adsorption kinetics study reveals that HCl adsorption into the ionized state is slightly favored over adsorption into the molecular states at 100 K, leading to earlier saturation of the ionized state. Between the two molecular states, the beta-state is formed first, and the alpha-state appears only at high HCl coverage. At 140 K, ionic dissociation of HCl is completed. The resulting hydronium ion can migrate into the underlying sublayer to a depth <4 bilayers, suggesting that the migration is assisted by self-diffusion of water molecules near the surface. When HCl is covered by a water overlayer at 100 K, its ionization efficiency is enhanced, but a substantial portion of HCl remains undissociated as molecules or contact ion pairs. The observation suggests that three-dimensional surrounding by water molecules does not guarantee ionic dissociation of HCl. Complete ionization of HCl requires additional thermal energy to separate the hydronium and chloride ions.

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.

Adsorption of HCl on the Water Ice Surface Studied by X-ray Absorption Spectroscopy

The adsorption state of HCl at 20 and 90 K on crystalline water ice films deposited under ultrahigh vacuum at 150 K has been studied by X-ray absorption spectroscopy at the O1s K-edge and Cl2p L-edge. We show that HCl dissociates at temperatures as low as 20 K, in agreement with the prediction of a spontaneous ionization of HCl on ice. Comparison between the rate of saturation of the "dangling" hydrogen bonds and the chlorine uptake indicates that hydrogen bonding of HCl with the surface native water "dangling" groups only accounts for a small part of the ionization events (20% at 90 K). A further mechanism drives the rest of the dissociation/solvation process. We suggest that the weakening of the ice surface hydrogen-bond network after the initial HCl adsorption phase facilitates the generation of new dissociation/solvation sites, which increases the uptake capacity of ice. These results also emphasize the necessity to take into account not only a single dissociation event but its catalyzing effect on the subsequent events when modeling the uptake of hydrogen-bonding molecules on the ice surface.

The vaporization rate of ice at temperatures near its melting point

The first study of free vaporization kinetics of ice at temperatures near its melting point is reported. The experimental approach employed is based on a unique combination of thermal desorption spectroscopy, microcalorimetry, and time-of-flight mass spectrometry, making it possible to overcome challenges associated with the introduction of volatile solids into a high vacuum environment. Measurements of the vaporization rate of polycrystalline ice demonstrate that the vaporization kinetics deviate dramatically from those predicted by a simple mobile precursor mechanism. The vaporization rate follows Arrhenius behavior from -40 to 0 degrees C with an effective activation energy of 50+/-4 kJ/mol, which is significantly higher than the value predicted by the simple mobile precursor mechanism. Extrapolation of earlier measurements conducted below -40 degrees C yields a value of approximately 0.02 at 0 degrees C for the vaporization coefficient alphav. In contrast, experimentally determined vaporization coefficient is found to be 0.7+/-0.3 and shows a weak dependence on temperature up to the bulk melting point. The role of possible surface phase transitions in the mechanisms of release and uptake of H2O and other chemical species by ice surfaces is discussed.

Monte Carlo simulation for the formation of a mixed crystal from two solids in contact

The study focuses on nucleation and growth of a binary mixed crystal phase from two pure crystals in contact. Monte Carlo simulations of this process are conducted, with the dynamics proceeding via activated atom-vacancy exchanges. Intermolecular interactions, ranging up to next-nearest neighbors, are of size typical of hydrogen bonded systems. The process is driven by the formation of strong AB bonds at the expense of weaker AA and BB bonds. In the resulting model, the material is channeled and transported through the mixed phase crust along antiphase boundaries. The flow of molecules through the channels is directed, due to molecular energy lowering via gradual acquisition of an increasing number of nearest neighbors of the second species. On the other hand, defect motion is quasirandom. The model accounts partially for the t(1/alpha) (alpha>3) time dependence observed for conversion of nanoparticles of HBr dihydrate to monohydrate, by exposure to acid adsorbate.

Ionization and solvation of HCl adsorbed on the D[sub 2]O-ice surface

The interaction of HCl with the D(2)O-ice surface has been investigated in the temperature range 15-200 K by utilizing time-of-flight secondary ion mass spectroscopy, temperature-programmed desorption, and x-ray photoelectron spectroscopy. The intensities of sputtered H(+)(D(2)O) and Cl(-) ions (the H(+) ions) are increased (decreased) markedly above 40 K due to the hydrogen bond formation between the HCl and D(2)O molecules. The HCl molecules which form ionic hydrates undergo H/D exchange at 110-140 K and a considerable fraction of them dissolves into the bulk above 140 K. The neutral hydrates of HCl should coexist as evidenced by the desorption of HCl above 170 K. They are incorporated completely in the D(2)O layer up to 140 K. The HCl molecules embedded in the thick D(2)O layer dissolve into the bulk, and the ionic hydrate tends to segregate to the surface above 150 K.