Molecular dynamics simulations of D2O ice photodesorption

Molecular dynamics (MD) calculations have been performed to study the ultraviolet (UV) photodissociation of D(2)O in an amorphous D(2)O ice surface at 10, 20, 60, and 90 K, in order to investigate the influence of isotope effects on the photodesorption processes. As for H(2)O, the main processes after UV photodissociation are trapping and desorption of either fragments or D(2)O molecules. Trapping mainly takes place in the deeper monolayers of the ice, whereas desorption occurs in the uppermost layers. There are three desorption processes: D atom, OD radical, and D(2)O molecule photodesorption. D(2)O desorption takes places either by direct desorption of a recombined D(2)O molecule, or when an energetic D atom produced by photodissociation kicks a surrounding D(2)O molecule out of the surface by transferring part of its momentum. Desorption probabilities are calculated for photoexcitation of D(2)O in the top four monolayers and are compared quantitatively with those for H(2)O obtained from previous MD simulations of UV photodissociation of amorphous water ice at different ice temperatures [Arasa et al., J. Chem. Phys. 132, 184510 (2010)]. The main conclusions are the same, but the average D atom photodesorption probability is smaller than that of the H atom (by about a factor of 0.9) because D has lower kinetic energy than H, whereas the average OD radical photodesorption probability is larger than that of OH (by about a factor of 2.5-2.9 depending on ice temperature) because OD has higher translational energy than OH for every ice temperature studied. The average D(2)O photodesorption probability is larger than that of H(2)O (by about a factor of 1.4-2.3 depending on ice temperature), and this is entirely due to a larger contribution of the D(2)O kick-out mechanism. This is an isotope effect: the kick-out mechanism is more efficient for D(2)O ice, because the D atom formed after D(2)O photodissociation has a larger momentum than photogenerated H atoms from H(2)O, and D transfers momentum more easily to D(2)O than H to H(2)O. The total (OD + D(2)O) yield has been compared with experiments and the total (OH + H(2)O) yield from previous simulations. We find better agreement when we compare experimental yields with calculated yields for D(2)O ice than when we compare with calculated yields for H(2)O ice.     

Overtone spectroscopy of H2O clusters in the V(OH) = 2 manifold: infrared-ultraviolet vibrationally mediated dissociation studies

Spectroscopy and predissociation dynamics of (H2O)2 and Ar-H2O are investigated with vibrationally mediated dissociation (VMD) techniques, wherein upsilon(OH) = 2 overtones of the complexes are selectively prepared with direct infrared pumping, followed by 193 nm photolysis of the excited H2O molecules. As a function of relative laser timing, the photolysis breaks H2O into OH and H fragments either (i) directly inside the complex or (ii) after the complex undergoes vibrational predissociation, with the nascent quantum state distribution of the OH photofragment probed via laser-induced fluorescence. This capability provides the first rotationally resolved spectroscopic analysis of (H2O)2 in the first overtone region and vibrational predissociation dynamics of water dimer and Ar-water clusters. The sensitivity of the VMD approach permits several upsilon(OH) = 2 overtone bands to be observed, the spectroscopic assignment of which is discussed in the context of recent anharmonic theoretical calculations.     

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.     

高振动激发态吡嗪碰撞传能的QCT计算研究

用准经典轨线(QCT)方法计算了高振动激发态吡嗪(C4N2H4)与N2、O2、NH3、基态吡嗪之间的碰撞传能. C4N2H4通过计算发现, 高振动激发态C4N2H4与N2、O2碰撞发生的主要是V-V传能, 与NH3碰撞发生的主要是V-R传能, 与基态C4N2H4碰撞发生的主要是V-V(R)传能. 通过比较高振动激发态C4N2H4、C6F6、C6H6与其基态分子的碰撞传能， 发现此类碰撞传能中, 若分子的对称性高， 则V-V传能更容易实现.     

Quasiclassical Trajectory Study of Collisional Energy Transfer between Highly Excited C6F6 and N2 ,O2 and Ground State C6F6

Quasiclassical trajectory calculation (QCT) is used frequently for studying collisional energy transfer between highly vibrationally excited molecules and bath gases. In this paper, the QCT of the energy transfer between highly vibrationally excited C6F6 and N2 ,O2 and ground state C6F6 were performed. The results indicate that highly vibrationally excited C6F6 transferred vibrational energy to vibrational distribution of N2, O2 and ground state C6F6, so they are V-V energy transfer. Especially it is mainly V-V resonance energy transfer between excited C6F6 andground state C6F6, excited C6F6 transfers more vibrational energy to ground state C6F6 than to N2 and O2. The values of QCT, - (ΔEvib) of excited C6F6 are smaller than those of experiments.     

Acceleration of the reaction OH + CO → H + CO2 by vibrational excitation of OH

The collision complex formed from a vibrationally excited reactant undergoes redissociation to the reactant, intramolecular vibrational relaxation (randomization of vibrational energy), or chemical reaction to the products. If attractive interaction between the reactants is large, efficient vibrational relaxation in the complex prevents redissociation to the reactants with the initial vibrational energy, and the complex decomposes to the reactants with low vibrational energy or converts to the products. In this paper, we have studied the branching ratios between the intramolecular vibrational relaxation and chemical reaction of an adduct HO(v)-CO formed from OH(X(2)Π(i)) in different vibrational levels v = 0-4 and CO. OH(v = 0-4) generated in a gaseous mixture of O(3)/H(2)/CO/He irradiated at 266 nm was detected with laser-induced fluorescence (LIF) via the A(2)Σ(+)-X(2)Π(i) transition, and H atoms were probed by the two-photon excited LIF technique. From the kinetic analysis of the time-resolved LIF intensities of OH(v) and H, we have found that the intramolecular vibrational relaxation is mainly governed by a single quantum change, HO(v)-CO → HO(v-1)-CO, followed by redissociation to OH(v-1) and CO. With the vibrational quantum number v, chemical process from the adduct to H + CO(2) is accelerated, and vibrational relaxation is decelerated. The countertrend is elucidated by the competition between chemical reaction and vibrational relaxation in the adduct HOCO.     

Hydrogen-bond disruption by vibrational excitations in water

An excitation of the OH-stretch nu(OH) of water has unique disruptive effects on the local hydrogen bonding. The disruption is not an immediate vibrational predissociation, which is frequently the case with hydrogen-bonded clusters, but instead is a delayed disruption caused by a burst of energy from a vibrationally excited water molecule. The disruptive effects are the result of a fragile hydrogen-bonding network subjected to a large amount of vibrational energy released in a short time by the relaxation of nu(OH) stretching and delta(H2O) bending excitations. The energy of a single nu(OH) vibration distributed over one, two, or three (classical) water molecules would be enough to raise the local temperature to 1100, 700, or 570 K, respectively. Our understanding of the properties of the metastable water state having this excess energy in nearby hydrogen bonds, termed H2O*, has emerged as a result of experiments where a femtosecond IR pulse is used to pump nu(OH), which is probed by either Raman or IR spectroscopy. These experiments show that the H2O* spectrum is blue-shifted and narrowed, and the spectrum looks very much like supercritical water at approximately 600 K, which is consistent with the temperature estimates above. The H2O* is created within approximately 400 fs after nu(OH) excitation, and it relaxes with an 0.8 ps lifetime by re-formation of the disrupted hydrogen-bond network. Vibrationally excited H2O* with one quantum of excitation in the stretching mode has the same 0.8 ps lifetime, suggesting it also relaxes by hydrogen-bond re-formation.     

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.     

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.     

Quantum dynamics of proton migration in H2O dications: H2+ formation on ultrafast timescales

Irradiation of isolated water molecules by few-cycle pulses of intense infrared laser light can give rise to ultrafast rearrangement resulting in formation of the H(2) (+) ion. Such unimolecular reactions occur on the potential energy surface of the H(2)O(2+) dication that is accessed when peak laser intensities in the 10(15) W cm(-2) range and pulse durations as short as 9-10 fs are used; ion yields of ~1.5% relative to the H(2)O(+) ion are measured. We also study such reactions by means of time-dependent wavepacket dynamics on an ab initio potential energy surface of the dication and show that a proton, generated from O-H bond rupture, migrates towards the H-atom, and forms vibrationally excited H(2)(+) in a well-defined spatial zone.     

Ionization dynamics of a water dimer: specific reaction selectivity

Ionization dynamics of a water dimer have been investigated by means of a direct ab initio molecular dynamics (MD) method. Two electronic state potential energy surfaces of (H(2)O)(2)(+) (ground and first excited states, (2)A' and (2)A') were examined as cationic states of (H(2)O)(2)(+). Three intermediate complexes were found as product channels. One is a proton transfer channel where a proton of H(2)O(+) is transferred into the H(2)O and then a complex composed of H(3)O(+)(OH) was formed. The second is a face-to-face complex channel denoted by (H(2)O-OH(2))(+) where the oxygen-oxygen atoms directly bind each other. Both water molecules are equivalent to each other. The third one is a dynamical complex where H(2)O(+) and H(2)O interact weakly and vibrate largely with a large intermolecular amplitude motion. The dynamics calculations showed that in the ionization to the (2)A' state, a proton transfer complex H(3)O(+)(OH) is only formed as a long-lived complex. On the other hand, in the ionization to the (2)A' state, two complexes, the face-to-face and dynamical complexes, were found as product channels. The proton of H(2)O(+) was transferred to H(2)O within 25-50 fs at the (2)A' state, meaning that the proton transfer on the ground state is a very fast process. On the other hand, the decay process on the first excited state is a slow process due to the molecular rotation. The mechanism of the ionization dynamics of (H(2)O)(2) was discussed on the basis of theoretical results.     

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.     

The OH radical-H2O molecular interaction potential

The OH radical is one of the most important oxidants in the atmosphere due to its high reactivity. The study of hydrogen-bonded complexes of OH with the water molecules is a topic of significant current interest. In this work, we present the development of a new analytical functional form for the interaction potential between the rigid OH radical and H(2)O molecules. To do this we fit a selected functional form to a set of high level ab initio data. Since there is a low-lying excited state for the H(2)O.OH complex, the impact of the excited state on the chemical behavior of the OH radical can be very important. We perform a potential energy surface scan using the CCSD(T)/aug-cc-pVTZ level of electronic structure theory for both excited and ground states. To model the physics of the unpaired electron in the OH radical, we develop a tensor polarizability generalization of the Thole-type all-atom polarizable rigid potential for the OH radical, which effectively describes the interaction of OH with H(2)O for both ground and excited states. The stationary points of (H(2)O)(n)OH clusters were identified as a benchmark of the potential.     

Theoretical determination of rate constants for vibrational relaxation and reaction of OH(X 2Pi, v = 1) with O(3P) atoms

Collisions of the vibrationally excited OH(v = 1) molecule with atomic oxygen are investigated theoretically using a coupled-states, statistical capture (CS-ST) model. Vibrational relaxation can occur by inelastic scattering, and the vibrationally excited molecule can also be removed by reaction to form O(2) in both the ground (X (3)Sigma(g)(-)) and first excited (a (1)Delta(g)) state. In the former case, reaction occurs on the lowest potential energy surface of (2)A(") symmetry, and, in the latter case, by reaction on the lowest potential energy surface of (2)A(') symmetry. We report new ab initio potential energy surfaces for both these states in the product and reactant regions necessary for application of the coupled-states, statistical method. Comparison with exact, reactive scattering calculations within the J-shifting approximation indicate that the CS-ST rate constants for removal of OH(v = 1) can be expected to be reasonably accurate. Our calculated rate constants at 300 K agree well with the experimental results of Khachatrian and Dagdigian [Chem. Phys. Lett. 415, 1 (2005)]. Reaction to yield O(2) (X (3)Sigma(g)(-)) is the dominant removal pathway. At subthermal temperatures, the rate constants for the various vibrational quenching processes all increase down to T approximately = 60 K and then decrease at lower temperature.     

Dissociative recombination of H(3)(+) in the ground and excited vibrational states

The article presents calculated dissociative recombination (DR) rate coefficients for H(3) (+). The previous theoretical work on H(3) (+) was performed using the adiabatic hyperspherical approximation to calculate the target ion vibrational states and it considered just a limited number of ionic rotational states. In this study, we use accurate vibrational wave functions and a larger number of possible rotational states of the H(3) (+) ground vibrational level. The DR rate coefficient obtained is found to agree better with the experimental data from storage ring experiments than the previous theoretical calculation. We present evidence that excited rotational states could be playing an important role in those experiments for collision energies above 10 meV. The DR rate coefficients calculated separately for ortho- and para-H(3) (+) are predicted to differ significantly at low energy, a result consistent with a recent experiment. We also present DR rate coefficients for vibrationally excited initial states of H(3) (+), which are found to be somewhat larger than the rate coefficient for the ground vibrational level.     

Electronic quenching of OH A 2Sigma+ radicals in single collision events with H2 and D2: a comprehensive quantum state distribution of the OH X 2Pi products

A pump-probe laser-induced fluorescence technique has been used to examine the nascent OH X (2)Pi product state distribution arising from non-reactive quenching of electronically excited OH A (2)Sigma(+) by molecular hydrogen and deuterium under single-collision conditions. The OH X (2)Pi products were detected in v'=0, 1 and 2; the distribution peaks in v'=0 and decreases monotonically with increasing vibrational excitation. In all vibrational levels probed, the OH X (2)Pi products are found to be highly rotationally excited, the distribution peaking at N'=15 when H(2) was used as the collision partner and N'=17 for D(2). A marked propensity for production of Pi(A') Lambda-doublet levels was observed, while both OH X (2)Pi spin-orbit manifolds were equally populated. These observations are interpreted as dynamical signatures of the nonadiabatic passage of the OH + H(2)/D(2) system through the seams of conical intersection that couple the excited state (2 (2)A') and ground state (1 (2)A') surfaces.     

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.     

Discrete variable approaches to tetratomic molecules: part II: application to H2O2 and H2CO

We have carried out large-scale calculations for accurate vibrational energy levels of formaldehyde and hydrogen peroxide. The discrete variable representations of the radial and angular coordinates are employed together with the contraction scheme resulting from several diagonalization/truncation steps. The global potential energy surface due to Carter et al. [J. Mol. Spectrosc. 90 (1997) 729] is used for H2CO and due to Koput et al. [J. Phys. Chem. A 102 (1998) 6325] for H2O2. For both molecules, the calculated vibrational energy levels are characterized by combining vibrationally averaged geometries and expectation values of rotational constants with several adiabatic projection schemes for automatic quantum number assignments. The energy levels of H2CO involving the excited v2 and v3 vibrations appear as resonances beyond the zero-order picture consisting of uncoupled 3D stretching and 2D bending modes. The torsional energy levels of H2O2 are studied in great detail and different energy patterns occurring below and above the cis barrier are discussed. Our full dimensional calculations for H2O2 have shown that the OH triad levels, 2vOH, are symmetry adapted local mode states.     

Quasiclassical trajectory study of the vibrational quenching of hydroxyl radicals through collision with O atoms

The collisional removal of vibrationally excited OH radicals by O atoms is studied by the quasiclassical trajectory method. To evaluate the effect of different topological features on the scattering processes two different global potential energy surfaces, DMBE IV and TU, are used. Results for reactive, exchange, and inelastic scattering probabilities are reported for central collisions (with zero total angular momentum) with a fixed relative translational energy for vibrational levels of OH ranging from nu=1 to v=8. Vibrational state distributions of product molecules are also compared on the two potential energy surfaces. Both surfaces predict higher probabilities for reaction than for exchange or inelastic scattering. The vibrational state distributions of the product diatomic molecules are different on the two surfaces. In particular, the two surfaces give substantially different probabilities for multiquantum OH vibrational relaxation transitions OH(v)+O-->OH(v')+O.     

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.