Molecular-dynamics study of photodissociation of water in crystalline and amorphous ices

We present the results of classical dynamics calculations performed to study the photodissociation of water in crystalline and amorphous ice surfaces at a surface temperature of 10 K. A modified form of a recently developed potential model for the photodissociation of a water molecule in ice [S. Andersson et al., Chem. Phys. Lett. 408, 415 (2005)] is used. Dissociation in the top six monolayers is considered. Desorption of H(2)O has a low probability (less than 0.5% yield per absorbed photon) for both types of ice. The final outcome strongly depends on the original position of the photodissociated molecule. For molecules in the first bilayer of crystalline ice and the corresponding layers in amorphous ice, desorption of H atoms dominates. In the second bilayer H atom desorption, trapping of the H and OH fragments in the ice, and recombination of H and OH are of roughly equal importance. Deeper into the ice H atom desorption becomes less important and trapping and recombination dominate. Motion of the photofragments is somewhat more restricted in amorphous ice. The distribution of distances traveled by H atoms in the ice peaks at 6-7 Angstroms with a tail going to about 60 Angstroms for both types of ice. The mobility of OH radicals is low within the ice with most probable distances traveled of 2 and 1 Angstrom for crystalline and amorphous ices, respectively. OH is, however, quite mobile on top of the surface, where it has been found to travel more than 80 Angstroms. Simulated absorption spectra of crystalline ice, amorphous ice, and liquid water are found to be in very good agreement with the experiments. The outcomes of photodissociation in crystalline and amorphous ices are overall similar, but with some intriguing differences in detail. The probability of H atoms desorbing is 40% higher from amorphous than from crystalline ice and the kinetic-energy distribution of the H atoms is on average 30% hotter for amorphous ice. In contrast, the probability of desorption of OH radicals from crystalline ice is much higher than that from amorphous ice.     

A theoretical and experimental study on translational and internal energies of H2O and OH from the 157 nm irradiation of amorphous solid water at 90 K

The photodesorption of H(2)O in its vibrational ground state, and of OH radicals in their ground and first excited vibrational states, following 157 nm photoexcitation of amorphous solid water has been studied using molecular dynamics simulations and detected experimentally by resonance-enhanced multiphoton ionization techniques. There is good agreement between the simulated and measured energy distributions. In addition, signals of H(+) and OH(+) were detected in the experiments. These are inferred to originate from vibrationally excited H(2)O molecules that are ejected from the surface by two distinct mechanisms: a direct desorption mechanism and desorption induced by secondary recombination of photoproducts at the ice surface. This is the first reported experimental evidence of photodesorption of vibrationally excited H(2)O molecules from water ice.     

Infrared spectroscopic observation of the group 13 metal hydroxides, M(OH)1,2,3 (M =Al, Ga, In, and Tl) and HAl(OH)2

Reactions of laser-ablated Al, Ga, In, and Tl atoms with H2O2 and with H2 + O2 mixtures diluted in argon give new absorptions in the O-H and M-O stretching and O-H bending regions, which are assigned to the metal mono-, di-, and trihydroxide molecules. Isotopic substitutions (D2O2, 18O2, 16,18O2, HD, and D2) confirm the assignments, and DFT calculations reproduce the experimental results. Infrared spectra for the Al(OH)(OD) molecule verify the calculated C2v structure. The trihydroxide molecules increase on annealing from the spontaneous reaction with a second H2O2 molecule. Aluminum atom reactions with the H2 + O2 mixtures favor the HAl(OH)2 product, suggesting that AlH3 generated by UV irradiation combines with O2 to form HAl(OH)2.     

Selective OD bond dissociation of HOD: photodissociation of vibrationally excited HOD in the 5nu(OD) state

Exclusively selective OD bond dissociation of HOD has been demonstrated by the ultraviolet photodissociation at 243.1 nm through the fourth overtone state of the OD stretching mode (5nu(OD)). Branching ratio between the OH and OD bond dissociation channels has been determined by detecting H and D atoms, utilizing a (2+1) resonance-enhanced multiphoton ionization (REMPI) process. The OD bond dissociation has been solely observed with the branching ratio phi(D+OH)/phi(H+OD)>12, which has been determined by the detection limit for the H atom. Time-dependent wave-packet calculations reveal two important features for the highly selective OD bond dissociation: (1) strong local-mode character of the 5nu(OD) state and (2) limitation of the total excitation energy lower than the saddle point between the OH and OD dissociation channels in the A state. Additionally, the recoil velocity and angular distribution of the nascent D atom are roughly evaluated by analyzing the Doppler-resolved REMPI spectrum. Based on these results, the dynamics of the selective OD dissociation has been discussed in detail.     

Isotopic effect on the cage-induced quenching of OH(A)/OD(A) inside small argon clusters

In this paper we report on the isotopic effect on the cage-induced excited-state quenching inside small Ar(m) clusters (m<10(2)) solvated in large Ne(N) clusters (N approximately 7.5x10(3)). Excited OH(A)/OD(A) fragments are produced by photodissociation of H2O and D2O molecules and the quenching agents are correspondingly H or D atoms. The decrease of the fluorescence yield with the size of the cluster m>m0 is observed in both cases and it is attributed to the formation of the cage of argon atoms around the doped molecule. Interestingly, more atoms are needed to induce the fluorescence quenching of OD*(A) fragments, m0=21+/-3, compared to the electronically excited state quenching of OH*(A) molecules, 11+/-2. A diffusion model containing two free parameters, the quenching cross section sigmaq and the number of argon atoms forming the cage m0, explains the effect in terms of the residence time of the hydrogen atom inside the cage. We suggest that the melting of the doped rare gas clusters is responsible for the different predissociation dynamics. The quenching cross section obtained from the experimental data is in good agreement with former experiments.     

Hot precursor reactions during the collisions of gas-phase oxygen atoms with deuterium chemisorbed on Pt(100)

We utilized direct rate measurements and temperature programmed desorption to investigate reactions that occur during the collisions of gaseous oxygen atoms with deuterium-covered Pt(100). We find that both D2O and D2 desorb promptly when an oxygen atom beam impinges upon D-covered Pt(100) held at surface temperatures ranging from 90 to 150 K, and estimate effective cross sections of 12 and 1.8 A2, respectively, for the production of gaseous D2O and D2 at 90 K. The yields of D2O and D2 that desorb at 90 K are about 13% and 2%, respectively, of the initial D atom coverage, though most of the D2O product molecules (approximately 80%) thermalize to the surface rather than desorb at the surface temperatures studied. Increasing the surface temperature from 90 to 150 K causes the D2O desorption rate to decay more quickly during O atom exposures to the surface and results in lower yields of gaseous D2O. We attribute the production of D2O and D2 in these experiments to reactions involving intermediates that are not thermally accommodated to the surface, so-called hot precursors. The results are consistent with the production of hot D2O involving first the generation of hot OD groups from the reaction O*+D(a)-->OD*, where the asterisk denotes a hot precursor, followed by the parallel pathways OD*+D(a)-->D2O* and OD*+OD(a)-->D2O*+O(a). The final reaction contributes significantly to hot D2O production only later in the reaction period when thermalized OD groups have accumulated on the surface, and it becomes less important at higher temperature due to depletion of the OD(a) concentration by thermally activated D2O production.     

Surface abundance change in vacuum ultraviolet photodissociation of CO2 and H2O mixture ices

Photodissociation of amorphous ice films of carbon dioxide and water co-adsorbed at 90 K was carried out at 157 nm using oxygen-16 and -18 isotopomers with a time-of-flight photofragment mass spectrometer. O((3)P(J)) atoms, OH (v = 0) radicals, and CO (v = 0,1) molecules were detected as photofragments. CO is produced directly from the photodissociation of CO(2). Two different adsorption states of CO(2), i.e., physisorbed CO(2) on the surface of amorphous solid water and trapped CO(2) in the pores of the film, are clearly distinguished by the translational and internal energy distributions of the CO molecules. The O atom and OH radical are produced from the photodissociation of H(2)O. Since the absorption cross section of CO(2) is smaller than that of H(2)O at 157 nm, the CO(2) surface abundance is relatively increased after prolonged photoirradiation of the mixed ice film, resulting in the formation of a heterogeneously layered structure in the mixed ice at low temperatures. Astrophysical implications are discussed.     

Intracluster stereochemistry in van der Waals complexes: steric effects in ultraviolet photodissociation of state-selected Ar-HOD/H2O

High-resolution IR-UV multiple resonance methods are employed to elucidate the photodissociation dynamics of quantum state-selected Ar-HOD and Ar-H(2)O van der Waals clusters. A single mode pulsed OPO operating in the region of the OH second overtone is used to prepare individual rovibrational states that are selectively photodissociated at specific excimer wavelengths. Subsequent fluorescence excitation of the resulting OH (OD) fragments yields dynamical information on the photofragmentation event and any resulting intracluster collisions. This technique is used to characterize spectroscopically the Pi(1(01)), nu(OH)=3<--Sigma(0(00)), v(OH)=0 overtone band of the Ar-HOD complex with an origin at 10648.27 cm(-1). The effects of Ar complexation on the dissociation dynamics are inferred by comparison of the OD photofragment quantum state distributions resulting from dissociation of single rovibrational states of the complex with those from isolated HOD photodissociation. The important role played by the initial internal state of the complex is demonstrated by comparison of the current Ar-HOD data with previously published results for the Ar-H(2)O Sigma(0(00))[03(-)> state. We interpret the dramatic differences in the dynamics of the two systems as manifestations of the nodal structure of the vibrational state in the parent complex and the way in which it governs the collision probability between the Ar atom and the escaping photofragments.     

Photoproduct channels from BrCD2CD2OH at 193 nm and the HDO + vinyl products from the CD2CD2OH radical intermediate

We present the results of our product branching studies of the OH + C(2)D(4) reaction, beginning at the CD(2)CD(2)OH radical intermediate of the reaction, which is generated by the photodissociation of the precursor molecule BrCD(2)CD(2)OH at 193 nm. Using a crossed laser-molecular beam scattering apparatus with tunable photoionization detection, and a velocity map imaging apparatus with VUV photoionization, we detect the products of the major primary photodissociation channel (Br and CD(2)CD(2)OH), and of the secondary dissociation of vibrationally excited CD(2)CD(2)OH radicals (OH, C(2)D(4)/CD(2)O, C(2)D(3), CD(2)H, and CD(2)CDOH). We also characterize two additional photodissociation channels, which generate HBr + CD(2)CD(2)O and DBr + CD(2)CDOH, and measure the branching ratio between the C-Br bond fission, HBr elimination, and DBr elimination primary photodissociation channels as 0.99:0.0064:0.0046. The velocity distribution of the signal at m/e = 30 upon 10.5 eV photoionization allows us to identify the signal from the vinyl (C(2)D(3)) product, assigned to a frustrated dissociation toward OH + ethene followed by D-atom abstraction. The relative amount of vinyl and Br atom signal shows the quantum yield of this HDO + C(2)D(3) product channel is reduced by a factor of 0.77 ± 0.33 from that measured for the undeuterated system. However, because the vibrational energy distribution of the deuterated radicals is lower than that of the undeuterated radicals, the observed reduction in the water + vinyl product quantum yield likely reflects the smaller fraction of radicals that dissociate in the deuterated system, not the effect of quantum tunneling. We compare these results to predictions from statistical transition state theory and prior classical trajectory calculations on the OH + ethene potential energy surface that evidenced a roaming channel to produce water + vinyl products and consider how the branching to the water + vinyl channel might be sensitive to the angular momentum of the β-hydroxyethyl radicals.     

Infrared and Raman line shapes of dilute HOD in liquid H2O and D2O from 10 to 90 degrees C

A combined electronic structure/molecular dynamics approach was used to calculate infrared and isotropic Raman spectra for the OH or OD stretches of dilute HOD in D2O or H2O, respectively. The quantities needed to compute the infrared and Raman spectra were obtained from density functional theory calculations performed on clusters, generated from liquid-state configurations, containing an HOD molecule along with 4-9 solvent water molecules. The frequency, transition dipole, and isotropic transition polarizability were each empirically related to the electric field due to the solvent along the OH (or OD) bond, calculated on the H (or D) atom of interest. The frequency and transition dipole moment of the OH (or OD) stretch of the HOD molecule were found to be very sensitive to its instantaneous solvent environment, as opposed to the isotropic transition polarizability, which was found to be relatively insensitive to environment. Infrared and isotropic Raman spectra were computed within a molecular dynamics simulation by using the empirical relationships and semiclassical expressions for the line shapes. The line shapes agree well with experiment over a temperature range from 10 to 90 degrees C.     

Segregation of hydroxide ions to an ice surface

Hydroxide ions that are initially buried within an ice film segregate to the ice film surface at elevated temperatures. This process was observed by conducting experiments with an ice film constructed with a bottom H(2)O layer and an upper D(2)O layer, with an excess of hydroxide ions trapped at the H(2)O/D(2)O interface as they were generated by Na hydrolysis. The transport of hydroxide ions from the interfacial layer to the surface was examined as a function of time using a low energy sputtering method. The progress of the H/D exchange reaction in surface water molecules was also monitored with the Cs(+) reactive ion scattering technique. At 90 K, only a small portion of buried hydroxide ions moved to the surface in the form of OD(-) species. This was due to hydroxide transport via proton hopping through a D(2)O layer, 3 BL thick, in the surface region. At 135 K, at which point water self-diffusion is active in the ice film, the majority of the buried hydroxide ions segregated to the surface after ～1 h. Both OH(-) and OD(-) species were produced at the surface, at an OH(-)/OD(-) population ratio ≥1. Based on kinetic measurements for the transport of OH(-) and OD(-) species and the H/D exchange of surface water molecules, we concluded that the major transport channel for hydroxide ions in this regime is the migration of molecular hydroxide species. H/D exchange reactions also occur between surface hydroxide ions and water molecules. No evidence was observed for the occurrence of the hop-and-turn process at 135 K, although it is known as an important mechanism of proton transport in ice.     

UV photodesorption of interstellar CO ice analogues: from subsurface excitation to surface desorption

Carbon monoxide is after H(2) the most abundant molecule identified in the interstellar medium (ISM), and is used as a major tracer for the gas phase physical conditions. Accreted at the surface of water-rich icy grains, CO is considered to be the starting point of a complex organic--presumably prebiotic--chemistry. Non-thermal desorption processes, and especially photodesorption by UV photons, are seen as the main cause that drives the gas-to-ice CO balance in the colder parts of the ISM. The process is known to be efficient and wavelength-dependent, but, the underlying mechanism and the physical-chemical parameters governing the photodesorption are still largely unknown. Using monochromatized photons from a synchrotron beamline, we reveal that the molecular mechanism responsible for CO photoejection is an indirect, (sub)surface-located process. The local environment of the molecules plays a key role in the photodesorption efficiency, and is quenched by at least an order of magnitude for CO interacting with a water ice surface.     

Wetting of mixed OHH(2)O layers on Pt(111)

We describe the effect of growth temperature and OHH(2)O composition on the wetting behavior of Pt(111). Changes to the desorption rate of ice films were measured and correlated to the film morphology using low energy electron diffraction and thermal desorption of chloroform to measure the area of multilayer ice and monolayer OHH(2)O exposed. Thin ice films roughen, forming bare (radical39 x radical39)R16 degrees water monolayer and ice clusters. The size of the clusters depends on growth temperature and determines their kinetic stability, with the desorption rate decreasing when larger clusters are formed by growth at high temperature. Continuous films of more than approximately 50 layers thick stabilize an ordered incommensurate ice film that does not dewet. OH coadsorption pins the first layer into registry with Pt, forming an ordered hexagonal (OH+H(2)O) structure with all the H atoms involved in hydrogen bonding. Although this layer has a similar honeycomb OH(x) skeleton to ice Ih, it is unable to reconstruct to match the bulk ice lattice parameter and does not form a stable wetting layer. Water aggregates to expose bare monolayer (OH+H(2)O), forming bulk ice crystallites whose size depend on preparation temperature. Increasing the proportion of water in the first layer provides free OH groups which stabilize the multilayer. The factors influencing multilayer wetting are discussed using density functional theory calculations to compare water adsorption on top of (OH+H(2)O) and on simple models for commensurate water structures. We show that both the (OH+H(2)O) structure and "H-down" water layers are poor proton acceptors, bonding to the first layer being enhanced by the presence of free OH groups. Formation of an ordered ice multilayer requires a water-metal interaction sufficient to wet the surface, but not so strong as to prevent the first layer relaxing to stabilize the interface between the metal and bulk ice.     

Hydrogen atom formation from the photodissociation of water ice at 193 nm

The TOF spectra of photofragment hydrogen atoms from the 193 nm photodissociation of amorphous ice at 90-140 K have been measured. The spectra consist of both a fast and a slow components that are characterized by average translational energies of 2k(B)T(trans)=0.39+/-0.04 eV (2300+/-200 K) and 0.02 eV (120+/-20 K), respectively. The incident laser power dependency of the hydrogen atom production suggests one-photon process. The electronic excitation energy of a branched cluster, (H(2)O)(6+1), has been theoretically calculated, where (H(2)O)(6+1) is a (H(2)O)(6) cyclic cluster attached by a water molecule with the hydrogen bond. The photoabsorption of this branched cluster is expected to appear at around 200 nm. The source of the hydrogen atoms is attributed to the photodissociation of the ice surface that is attached by water molecules with the hydrogen bond. Atmospheric implications are estimated for the photodissociation of the ice particles (Noctilucent clouds) at 190-230 nm in the region between 80 and 85 km altitude.     

Release of hydrogen molecules from the photodissociation of amorphous solid water and polycrystalline ice at 157 and 193 nm

The production of H(2) in highly excited vibrational and rotational states (v=0-5, J=0-17) from the 157 nm photodissociation of amorphous solid water ice films at 100 K was observed directly using resonance-enhanced multiphoton ionization. Weaker signals from H(2)(v=2,3 and 4) were obtained from 157 nm photolysis of polycrystalline ice, but H(2)(v=0 and 1) populations in this case were below the detection limit. The H(2) products show two distinct formation mechanisms. Endothermic abstraction of a hydrogen atom from H(2)O by a photolytically produced H atom yields vibrationally cold H(2) products, whereas exothermic recombination of two H-atom photoproducts yields H(2) molecules with a highly excited vibrational distribution and non-Boltzmann rotational population distributions as has been predicted previously by both quantum-mechanical and molecular dynamics calculations.     

Site-dependent electron-stimulated reactions in water films on TiO2(110)

Electron-stimulated reactions in thin [<3 ML (monolayer)] water films adsorbed on TiO(2)(110) are investigated. Irradiation with 100 eV electrons results in electron-stimulated dissociation and electron-stimulated desorption (ESD) of adsorbed water molecules. The molecular water ESD yield increases linearly with water coverage theta for 0< or =theta< or =1 ML and 11 ML, the water ESD yield per additional water molecule adsorbed (i.e., the slope of the ESD yield versus coverage) is 3.5 times larger than for theta<1 ML. In contrast, the number of water molecules dissociated per incident electron increases linearly for theta< or =2 ML without changing slope at theta=1 ML. The total electron-stimulated sputtering rate, as measured by postirradiation temperature programmed desorption of the remaining water, is larger for theta>1 ML due to the increased water ESD for those coverages. The water ESD yields versus electron energy (for 5-50 eV) are qualitatively similar for 1, 2, and 40 ML water films. In each case, the observed ESD threshold is at approximately 10 eV and the yield increases monotonically with increasing electron energy. The results indicate that excitations in the adsorbed water layer are primarily responsible for the ESD in thin water films on TiO(2)(110). Experiments on "isotopically layered" films with D(2)O adsorbed on the Ti(4+) sites (D(2)O(Ti)) and H(2)O adsorbed on the bridging oxygen atoms (H(2)O(BBO)) demonstrate that increasing the water coverage above 1 ML rapidly suppresses the electron-stimulated desorption of D(2)O(Ti) and D atoms, despite the fact that the total water ESD and atomic hydrogen ESD yields increase with increasing coverage. The coverage dependence of the electron-stimulated reactions is probably related to the different bonding geometries for H(2)O(Ti) and H(2)O(BBO) and its influence on the desorption probability of the reaction products.     

Photochemical charge transfer and trapping at the interface between an organic adlayer and an oxide semiconductor

Identification of charge transfer and trapping sites on semiconducting oxide surfaces is of fundamental importance in furthering the field of heterogeneous photocatalysts. Using scanning tunneling microscopy, electron energy loss spectroscopy, and photodesorption, we observed both electron trapping and hole transfer events on the (110) surface of TiO2 rutile. UV irradiation of a saturated monolayer of trimethyl acetate (TMA) on TiO2(110) at room temperature resulted in hole transfer to the carboxylate group, followed by (CH3)3C-COO bond cleavage and desorption of CO2 and isobutene/isobutane. Hole transfer to TMA proceeded in the absence of a gas-phase electron scavenger (which is typically O2) because the accompanying photogenerated electrons could be trapped at the surface as Ti3+ cations bound to bridging OH groups. The extent of electron trapping, gauged by electron spectroscopy, correlated directly with the yields of photodesorption fragments resulting from the hole transfer channel. Charge at the Ti3+ sites was titrated in the dark via a reaction between O2 and the Ti3+-OH groups.     

Methane activation by titanium monoxide molecules: a matrix isolation infrared spectroscopic and theoretical study

Reactions of titanium monoxides with methane have been investigated using matrix isolation infrared spectroscopy and theoretical calculations. Titanium derivatives of several simple oxyhydrocarbons have been prepared and identified. The titanium monoxide molecules prepared by laser evaporation of bulk TiO2 target reacted with methane to form the TiO(CH4) complex in solid argon, which was predicted to have C3v symmetry with the oxygen atom coordinated to one hydrogen atom of the methane molecule. The complex rearranged to the CH3Ti(O)H titano-acetaldehyde molecule upon visible (lambda > 500 nm) irradiation. The titano-acetaldehyde molecule sustained further photochemical rearrangement to the CH2Ti(H)OH titano-vinyl alcohol molecule, which was characterized to be a simple carbene complex involving agostic bonding. The CH2Ti(H)OH molecule reacted with a second methane to form the (CH3)2Ti(H)OH titano-isopropyl alcohol molecule spontaneously on annealing. The (CH3)2Ti(H)OH molecule also can be produced via UV photon-induced rearrangement of the CH3Ti(O)H(CH4) complex.     

Why does argon bind to deuterium? Isotope effects and structures of Ar x H5O2(+) complexes

Recently, we reported the spectrum of Ar x D4HO2(+) [McCunn; et, al. J. Phys. Chem. B 2008, 112, 321], and here, we extend that work to include the Ar x H4DO2(+) isotopologue in order to explore why the Ar atom has a much greater propensity for attachment to a dangling OD group than it does for OH, even when many more of the latter binding sites are available. Calculated (MP2/6-311+G(d,p) level of theory/basis) harmonic frequencies reproduce the observed multiplet patterns of OH and OD stretches and confirm the presence of various isomers arising from the different Ar binding sites. The preferential bonding of Ar to OD is traced to changes in the frequencies of the wag and rock modes of the H5O2(+) moiety rather than to shifts in the oscillator that directly binds the Ar atom.     

Full-dimensional quantum dynamics study of exchange processes for the D + H2O and D + HOD reactions

The exchange processes of D + H(2)O and D + HOD reactions are studied using initial state-selected time-dependent wave packet approach in full dimension. The total reaction probabilities for different partial waves, together with the integral cross sections, are obtained both by the centrifugal sudden (CS) approximation and exact coupled-channel (CC) calculations, for the H(2)O(HOD) reactant initially in the ground rovibrational state. In the CC calculations, small resonance peaks in the reaction probabilities and quick diminishing of the resonance peaks with the increase of total angular momenta J do not lead to clear step-like features just above the threshold in the cross sections for the title reactions, which are different in other isotopically substituted reactions where the hydrogen atom was included as the reactant instead of the deuterium atom [B. Fu, Y. Zhou, and D. H. Zhang, Chem. Sci. 3, 270 (2012); B. Fu and D. H. Zhang, J. Phys. Chem. A 116, 820 (2012)]. It is interesting that the shape resonance-induced features resulting from the reaction tunneling are significantly diminished accordingly in the reactions of the deuterium atom and H(2)O or HOD, owing to the weaker tunneling capability of the reagent deuterium atom in the title reactions than the reagent hydrogen atom in other reactions. In the CS calculations, the resonance peaks persist in many partial waves but cannot survive the partial-wave summations. The cross sections for the D(') + H(2)O → D(')OH + H and D(') + HOD → D(')OD + H reactions are substantially larger than those for the D(') + HOD → HOD(') + D reaction, indicating that the D(')/H exchange reactions are much more favored than the D(')/D exchange.