O(3P) + CO2 collisions at hyperthermal energies: dynamics of nonreactive scattering, oxygen isotope exchange, and oxygen-atom abstraction

The dynamics of O((3)P) + CO(2) collisions at hyperthermal energies were investigated experimentally and theoretically. Crossed-molecular-beams experiments at = 98.8 kcal mol(-1) were performed with isotopically labeled (12)C(18)O(2) to distinguish products of nonreactive scattering from those of reactive scattering. The following product channels were observed: elastic and inelastic scattering ((16)O((3)P) + (12)C(18)O(2)), isotope exchange ((18)O + (16)O(12)C(18)O), and oxygen-atom abstraction ((18)O(16)O + (12)C(18)O). Stationary points on the two lowest triplet potential energy surfaces of the O((3)P) + CO(2) system were characterized at the CCSD(T)/aug-cc-pVTZ level of theory and by means of W4 theory, which represents an approximation to the relativistic basis set limit, full-configuration-interaction (FCI) energy. The calculations predict a planar CO(3)(C(2v), (3)A') intermediate that lies 16.3 kcal mol(-1) (W4 FCI excluding zero point energy) above reactants and is approached by a C(2v) transition state with energy 24.08 kcal mol(-1). Quasi-classical trajectory (QCT) calculations with collision energies in the range 23-150 kcal mol(-1) were performed at the B3LYP/6-311G(d) and BMK/6-311G(d) levels. Both reactive channels observed in the experiment were predicted by these calculations. In the isotope exchange reaction, the experimental center-of-mass (c.m.) angular distribution, T(θ(c.m.)), of the (16)O(12)C(18)O products peaked along the initial CO(2) direction (backward relative to the direction of the reagent O atoms), with a smaller isotropic component. The product translational energy distribution, P(E(T)), had a relatively low average of = 35 kcal mol(-1), indicating that the (16)O(12)C(18)O products were formed with substantial internal energy. The QCT calculations give c.m. P(E(T)) and T(θ(c.m.)) distributions and a relative product yield that agree qualitatively with the experimental results, and the trajectories indicate that exchange occurs through a short-lived CO(3)* intermediate. A low yield for the abstraction reaction was seen in both the experiment and the theory. Experimentally, a fast and weak (16)O(18)O product signal from an abstraction reaction was observed, which could only be detected in the forward direction. A small number of QCT trajectories leading to abstraction were observed to occur primarily via a transient CO(3) intermediate, albeit only at high collision energies (149 kcal mol(-1)). The oxygen isotope exchange mechanism for CO(2) in collisions with ground state O atoms is a newly discovered pathway through which oxygen isotopes may be cycled in the upper atmosphere, where O((3)P) atoms with hyperthermal translational energies can be generated by photodissociation of O(3) and O(2).     

Experimental and theoretical investigations of the inelastic and reactive scattering dynamics of O(3p) + D2

This paper presents a combined experimental and theoretical study of the dynamics of O((3)P) + D(2) collisions, with emphasis on a center-of-mass (c.m.) collision energy of 25 kcal mol(-1). The experiments were conducted with a crossed-molecular-beams apparatus, employing a laser detonation source to produce hyperthermal atomic oxygen and mass spectrometric detection to measure the product angular and time-of-flight distributions. The novel beam source, which enabled these experiments to be conducted, contributed unique challenges to the experiments and to the analysis, so the experimental methods and approach to the analysis are discussed in detail. Three different levels of theory were used: (1) quasiclassical trajectories (QCT), (2) time-independent quantum scattering calculations based on high-quality potential surfaces for the two lower-energy triplet states, and (3) trajectory-surface-hopping (TSH) studies that couple the triplet surfaces with the lowest singlet surface using a spin-orbit Hamiltonian derived from ab-initio calculations. The latter calculations explore the importance of intersystem crossing in the dynamics. Both experiment and theory show that inelastically scattered O atoms scatter almost exclusively in the forward direction, with little or no loss of translational energy. For the reaction, O((3)P) + D(2) --> OD + D, the experiment shows that, on average, approximately 50% of the available energy goes into product translation and that the OD product angular distributions are largely backward-peaked. These results may be interpreted in light of the QCT and TSH calculations, leading to the conclusion that the reaction occurs mainly on triplet potential energy surfaces with, at most, minor intersystem crossing to a singlet surface. Reaction on either of the two low-lying reactive triplet surfaces proceeds through a rebound mechanism in which the angular distributions are backward-peaked and the product OD is both vibrationally and rotationally excited. The quantum scattering results are in good agreement with QCT calculations, indicating that quantum effects are relatively small for this reaction at a collision energy of 25 kcal mol(-1).     

Dynamics study of the reaction OH- + C2H2-->C2H- + H2O with crossed beams and density-functional theory calculations

The proton transfer reaction between OH- and C2H2, the sole reactive process observed over the collision energy range from 0.37 to 1.40 eV, has been studied using the crossed beam technique and density-functional theory (DFT) calculations. The center of mass flux distributions of the product C2H- ions at three different energies are highly asymmetric, characteristic of a direct process occurring on a time scale much less than a rotational period of any transient intermediate. The maxima in the flux distributions correspond to product velocities and directions close to those of the precursor acetylene reactants. The reaction quantitatively transforms the entire exothermicity into internal excitation of the products, consistent with an energy release motif in which the proton is transferred early, in a configuration in which the forming bond is extended. This picture is supported by DFT calculations showing that the first electrostatically bound intermediate on the reaction pathway is the productlike C2H- H2O species. Most of the incremental translational energy in the two higher collision energy experiments appears in product translational energy, and provides an example of induced repulsive energy release characteristic of the heavy+light-heavy mass combination.     

State-resolved reactive scattering by slice imaging: a new view of the Cl+C2H6 reaction

We present state-resolved crossed beam scattering results for the reaction Cl+C2H6-->HCl+C2H5, obtained using direct current slice imaging. The HCl (v=0,J=2) image, recorded at a collision energy of 6.7+/-0.6 kcalmol, shows strongly coupled angular and translational energy distributions revealing features of the reaction not seen in previous studies. The overall distribution is mainly forward scattered with respect to the Cl beam, with a translational energy distribution peaking near the collision energy. However, there is a substantial backscattered contribution that is very different. It shows a sharp peak at 8.0 kcalmol, but extends to much lower energy, implying substantial internal excitation in the ethyl radical coproduct. These results provide new insight into the reaction, and they are considered in terms of alternative models of the dynamics. This work represents the first genuine crossed-beam study in which a product other than the methyl radical was detected with quantum state specificity, showing the promise of the approach generally for high resolution state-resolved reactive scattering.     

Vibrational excitation of SF6 scattering from graphite

We have observed vibrationally excited sulfur hexafluoride molecules in direct inelastic scattering from hot graphite surfaces. The vibrational temperature for the scattered flux has been determined by probing the effect of internal temperature on electron-induced fragmentation observed in mass spectra. The vibrational excitation depends on incident translational energy, E_tr, and a maximum temperature increase of 50 K is reached in direct scattering at E_tr = 2.5 eV. No effect of surface temperature has been observed at 950–1400 K. Inelastic angular distributions are reproduced by a collision complex model, and the experimental results are related to existing models for vibrational excitation.     

Reaction dynamics of OH+(3Sigma-)+C2H2 studied with crossed beams and density functional theory calculations

The reactions between OH+(3Sigma-) and C2H2 have been studied using crossed ion and molecular beams and density functional theory calculations. Both charge transfer and proton transfer channels are observed. Products formed by carbon-carbon bond cleavage analogous to those formed in the isoelectronic O(3P)+C2H2 reaction, e.g., 3CH2 + HCO+, are not observed. The center of mass flux distributions of both product ions at three different energies are highly asymmetric, with maxima close to the velocity and direction of the precursor acetylene beam, characteristic of direct reactions. The internal energy distributions of the charge transfer products are independent of collision energy and are peaked at the reaction exothermicity, inconsistent with either the existence of favorable Franck-Condon factors or energy resonance. In proton transfer, almost the entire reaction exothermicity is transformed into product internal excitation, consistent with mixed energy release in which the proton is transferred with both the breaking and forming bonds extended. Most of the incremental translational energy in the two higher-energy experiments appears in product translational energy, providing an example of induced repulsive energy release.     

Reactive quenching of OD A 2Σ+ by H2: translational energy distributions for H- and D-atom product channels

The H- and D-atom products from collisional quenching of OD A (2)Σ(+) by H(2) are characterized through Doppler spectroscopy using two-photon (2 (2)S ←← 1 (2)S) laser-induced fluorescence. Partial deuteration enables separation of the channel forming H + HOD products, which accounts for 75% of reactive quenching events, from the D + H(2)O product channel. The Doppler profiles, along with those reported previously for other isotopic variants, are transformed into product translational energy distributions using a robust fitting procedure based on discrete velocity basis functions. The product translational energy distribution for the H-atom channel is strongly peaked at low energy (below 0.5 eV) with a long tail extending to the energetic limit. By contrast, the D-atom channel exhibits a small peak at low translational energy with a distinctive secondary peak at higher translational energy (approximately 1.8 eV) before falling off to higher energy. In both cases, most of the available energy flows into internal excitation of the water products. Similar distributions are obtained upon reanalysis of D- and H-atom Doppler profiles, respectively, from reactive quenching of OH A (2)Σ(+) by D(2). The sum of the translational energy distributions for H- and D-atom channels is remarkably similar to that obtained for OH A (2)Σ(+) + H(2), where the two channels cannot be distinguished from one another. The product translational energy distributions from reactive quenching are compared with those obtained from a previous experiment performed at higher collision energy, quasiclassical trajectory calculations of the post-quenching dynamics, and a statistical model.     

Crossed beams and theoretical studies of the dynamics of hyperthermal collisions between Ar and ethane

Crossed molecular beams experiments and classical trajectory calculations have been used to study the dynamics of Ar+ethane collisions at hyperthermal collision energies. Experimental time-of-flight and angular distributions of ethane molecules that scatter into the backward hemisphere (with respect to their original direction in the center-of-mass frame) have been collected. Translational energy distributions, derived from the time-of-flight distributions, reveal that a substantial fraction of the collisions transfer abnormally large amounts of energy to internal excitation of ethane. The flux of the scattered ethane molecules increased only slightly from directly backward scattering to sideways scattering. Theoretical calculations show angular and translational energy distributions which are in reasonable agreement with the experimental results. These calculations have been used to examine the microscopic mechanism for large energy transfer collisions ("supercollisions"). Collinear ("head-on") or perpendicular ("side-on") approaches of Ar to the C-C axis of ethane do not promote energy transfer as much as bent approaches, and collisions in which the H atom is "sandwiched" in a bent Ar...H-C configuration lead to the largest energy transfer. The sensitivity of collisional energy transfer to the intramolecular potential energy of ethane has also been examined.     

Reaction dynamics of CN + O2 --> NCO + O(3P2)

We have used oxygen Rydberg time-of-flight spectroscopy to carry out a crossed molecular beam study of the CN + O2 reaction at collision energies of 3.1 and 4.1 kcal/mol. The O(3P2) products were tagged by excitation to high-n Rydberg levels and subsequently field ionized at a detector. The translational energy distributions were broad, indicating that the NCO is formed with a wide range of internal excitation, and the angular distribution was forward-backward symmetric, indicating the participation of NCOO intermediates with lifetimes comparable to or longer than their rotational periods. Rice-Ramsperger-Kassel-Marcus modeling of the dissociation of NCOO to NCO + O suggests that Do(NC-OO) > or = 38 kcal/mol, which is consistent with several theoretical calculations. Implications for the competing CO + NO channel are discussed.     

Rovibrational product state distribution for inelastic H+D2 collisions

Experimental measurements of rovibrational product state distributions for the inelastic scattering process H + D2(nu=0,j)-->H + D2(nu' = 1,2,j') are presented and compared with the results of quasiclassical and quantum mechanical calculations. Agreement between theory and experiment is almost quantitative. Two subtle trends are found: the relative amount of energy in product rotational excitation decreases slightly with increasing collision energy and increases slightly with increasing product vibrational excitation. These trends are the reverse of what has been found for reactive scattering in which the opposite trends are much more pronounced.     

Dynamics of the reaction H2+(He,H)HeH+. Influence of various forms of reactant energy on the total and differential cross section

Results of quasiclassical trajectory calculations of reactive processes between He atoms and H₂⁺ (υ, J) molecular ions in the collision energy interval 0.5–5.0 eV (c.m.) for a large number of selected υ, J combinations are analyzed with respect to the influence of the initial translational, vibrational, and rotational energy on the total and differential reaction cross sections. Vibrational energy is more effective in promoting the reaction than translational energy. Small rotational excitation has a negligible effect, whereas high rotational excitation has a similar influence on the reaction cross sections as the vibrational excitation of the same magnitude.     

Monte Carlo simulation of RRKM unimolecular decomposition in molecular beam experiments. III. Application to angular distributions and branching fractions for the system K + CsF

The experiments on the reactive systems K + CsF and K + RbF by Stolte et al. are unique in that product angular distributions and branching fractions with variation of both translational and rotational energies have been obtained. The results for the reaction K + CsF with translational energy variation have now been re-evaluated using an algorithm recently published by us. A satisfactory fit to the experimental data is found by a long-lived collision complex model including strict angular momentum conservation in both the reactive and the non-reactive channels. No “osculating” behaviour has to be added to the model. Angular momentum limits have to be included in the model to reach agreement with the measured branching fractions and cross sections. This indicates that the internal dynamics of the complexes is important. The calculation shows that the reactants giving reactive complexes have large rotational and vibrational energies, while their translational energy is nearly average. The internal energy of reactive pairs of molecules decreases with increasing translational energy. The question how rapidly the branching fraction varies with translational energy is thus more complicated to answer than believed previously.     

Experimental and theoretical studies of charge transfer and deuterium ion transfer between D2O+ and C2H4

The charge transfer and deuterium ion transfer reactions between D(2)O(+) and C(2)H(4) have been studied using the crossed beam technique at relative collision energies below one electron volt and by density functional theory (DFT) calculations. Both direct and rearrangement charge transfer processes are observed, forming C(2)H(4) (+) and C(2)H(3)D(+), respectively. Independent of collision energy, deuterium ion transfer accounts for approximately 20% of the reactive collisions. Between 22 and 36 % of charge transfer collisions occur with rearrangement. In both charge transfer processes, comparison of the internal energy distributions of products with the photoelectron spectrum of C(2)H(4) shows that Franck-Condon factors determine energy disposal in these channels. DFT calculations provide evidence for transient intermediates that undergo H/D migration with rearrangement, but with minimal modification of the product energy distributions determined by long range electron transfer. The cross section for charge transfer with rearrangement is approximately 10(3) larger than predicted from the Rice-Ramsperger-Kassel-Marcus isomerization rate in transient complexes, suggesting a nonstatistical mechanism for H/D exchange. DFT calculations suggest that reactive trajectories for deuterium ion transfer follow a pathway in which a deuterium atom from D(2)O(+) approaches the pi-cloud of ethylene along the perpendicular bisector of the C-C bond. The product kinetic energy distributions exhibit structure consistent with vibrational motion of the D-atom in the bridged C(2)H(4)D(+) product perpendicular to the C-C bond. The reaction quantitatively transforms the reaction exothermicity into internal excitation of the products, consistent with mixed energy release in which the deuterium ion is transferred in a configuration in which both the breaking and the forming bonds are extended.     

Crossed molecular beam studies on the reaction dynamics of O(1D)+N2O

The reaction of oxygen atom in its first singlet excited state with nitrous oxide was investigated under the crossed molecular beam condition. This reaction has two major product channels, NO+NO and N2+O2. The product translational energy distributions and angular distributions of both channels were determined. Using oxygen-18 isotope labeled O(1D) reactant, the newly formed NO can be distinguished from the remaining NO that was contained in the reactant N2O. Both channels have asymmetric and forward-biased angular distributions, suggesting that there is no long-lived collision complex with lifetime longer than its rotational period. The translational energy release of the N2+O2 channel (fT = 0.57) is much higher than that of the NO+NO channel (fT = 0.31). The product energy partitioning into translational, rotational, and vibrational degrees of freedom is discussed to learn more about the reaction mechanism. The branching ratio between the two product channels was estimated. The 46N2O product of the isotope exchange channel, 18O+44N2O-->16O+46N2O, was below the detection limit and therefore, the upper limit of its yield was estimated to be 0.8%.     

Imaging H2O photofragments in the predissociation of the HCl-H2O hydrogen-bonded dimer

The state-to-state vibrational predissociation (VP) dynamics of the hydrogen-bonded HCl-H(2)O dimer was studied following excitation of the dimer's HCl stretch by detecting the H(2)O fragment. Velocity map imaging (VMI) and resonance-enhanced multiphoton ionization (REMPI) were used to determine pair-correlated product energy distributions. Following vibrational excitation of the HCl stretch of the dimer, H(2)O fragments were detected by 2 + 1 REMPI via the C (1)B(1) (000) ← X (1)A(1) (000) transition. REMPI spectra clearly show H(2)O from dissociation produced in the ground vibrational state. The fragments' center-of-mass (c.m.) translational energy distributions were determined from images of selected rotational states of H(2)O and were converted to rotational state distributions of the HCl cofragment. The distributions were consistent with the previously measured dissociation energy of D(0) = 1334 ± 10 cm(-1) and show a clear preference for rotational levels in the HCl fragment that minimize translational energy release. The usefulness of 2 + 1 REMPI detection of water fragments is discussed.     

Molecular beam scattering studies of the reaction D + H2 (v = 0) and D + H2 (v = 1) → HD + H

The reaction D + H₂ → HD + H has been investigated in two molecular beam scattering experiments. Angular and time-of-flight distributions have been measured for the initial vibrational ground state (v = 0) at a most probable collision energy of E_cm = 1.5 eV and for the first vibrational excited state (v = 1) at E_cm = 0.28 eV with the same apparatus. Results for the ground-state experiment are compared with quasiclassical trajectory calculations(QCT) on the LSTH-hypersurface transformed into the laboratory system and averaged over the apparatus distributions. The agreement isquite satisfactory. At this high collision energy the HD products are no longer scattered in a backward direction but in a wide angular region concentrated about θ = 90° in the center-of-mass system. The absolute reactive cross section has been determined and the agreement with the theoretical value from QCT calculations is within the experimental error. The high sensitivity of the experiment to different properties of the doubly differential cross section has also been demonstrated. A preliminary evaluation of the experiment with initial vibrational excitation (v = 1) shows that the HD-product molecules are preferably backward scattered and the change of internal energy is small supporting the concept of a reaction which is adiabatic with respect to the internal degrees of freedom.     

Dynamics of collisions of hydroxyl radicals with fluorinated self-assembled monolayers

We present a classical trajectory study of the dynamics of collisions between OH radicals and fluorinated self-assembled monolayers (F-SAMs). The gas/surface interaction potential required in the simulations has been derived from high-level ab initio calculations (focal-point-CCSD(T)/aug-cc-pVQZ) of various approaches of OH to a model fluorinated alkane. The two lowest-energy doublet potential energy surfaces considered in the electronic structure calculations have been averaged to produce a pairwise analytic potential. This analytic potential has been subsequently employed to propagate classical trajectories of collisions between OH and F-SAMs at initial conditions relevant to recent experiments on related systems. The calculated rotational distributions of the inelastically scattered OH agree well with the experiment, which serves to validate the accuracy of the simulations. Investigation of the dynamics of energy transfer for different initial rotational states of OH indicates that an increase in the initial rotation of OH results in increases in both the final average OH rotational and translational energy and in a slight decrease in the amount of energy transferred to the surface. Analysis of the dynamics as a function of the desorption angle of OH from the surface shows that while there is a correlation between the final scattering angle and OH’s amount of final translational energy, the amount of rotational energy in OH is largely independent of the desorption angle. The mechanism of the collisions is found to be mostly direct; in about 90% of most trajectories, OH only collides with the surface once before desorbing, which exemplifies the rigidity of fluorinated monolayer surfaces and their inability to efficiently accommodate gas species.     

Experimental study of electron impact excitation of multiply charged ions

We report on measurements of angular differential cross sections for the excitation of multiply charged ions by electron impact. An ion beam is crossed by an electron beam; electrons which are inelastically scattered at different angles are identified by their energy loss due to the excitation process. Absolute excitation cross sections are obtained by comparing the signals of the elastic and the inelastic electron-ion scattering. Results obtained for the 3s→3p excitation of Ar⁷⁺ are discussed.     

Quantum state-resolved energy transfer dynamics at gas-liquid interfaces: IR laser studies of CO2 scattering from perfluorinated liquids

An apparatus for detailed study of quantum state-resolved inelastic energy transfer dynamics at the gas-liquid interface is described. The approach relies on supersonic jet-cooled molecular beams impinging on a continuously renewable liquid surface in a vacuum and exploits sub-Doppler high-resolution laser absorption methods to probe rotational, vibrational, and translational distributions in the scattered flux. First results are presented for skimmed beams of jet-cooled CO(2) (T(beam) approximately 15 K) colliding at normal incidence with a liquid perfluoropolyether (PFPE) surface at E(inc) = 10.6(8) kcal/mol. The experiment uses a tunable Pb-salt diode laser for direct absorption on the CO(2) nu(3) asymmetric stretch. Measured rotational distributions in both 00(0)0 and 01(1)0 vibrational manifolds indicate CO(2) inelastically scatters from the liquid surface into a clearly non-Boltzmann distribution, revealing nonequilibrium dynamics with average rotational energies in excess of the liquid (T(s) = 300 K). Furthermore, high-resolution analysis of the absorption profiles reveals that Doppler widths correspond to temperatures significantly warmer than T(s) and increase systematically with the J rotational state. These rotational and translational distributions are consistent with two distinct gas-liquid collision pathways: (i) a T approximately 300 K component due to trapping-desorption (TD) and (ii) a much hotter distribution (T approximately 750 K) due to "prompt" impulsive scattering (IS) from the gas-liquid interface. By way of contrast, vibrational populations in the CO(2) bending mode are inefficiently excited by scattering from the liquid, presumably reflecting much slower T-V collisional energy transfer rates.     

Imaging the state-specific vibrational predissociation of the C2H2-NH3 hydrogen-bonded dimer

The state-to-state vibrational predissociation (VP) dynamics of the hydrogen-bonded ammonia-acetylene dimer were studied following excitation in the asymmetric CH stretch. Velocity map imaging (VMI) and resonance-enhanced multiphoton ionization (REMPI) were used to determine pair-correlated product energy distributions. Following vibrational excitation of the asymmetric CH stretch fundamental, ammonia fragments were detected by 2 + 1 REMPI via the B1E' <-- X1A1' and C'1A1' <-- X1A1' transitions. The fragments' center-of-mass (c.m.) translational energy distributions were determined from images of selected rotational levels of ammonia with one or two quanta in the symmetric bend (nu2 umbrella mode) and were converted to rotational-state distributions of the acetylene co-fragment. The latter is always generated with one or two quanta of bending excitation. All the distributions could be fit well when using a dimer dissociation energy of D0 = 900 +/- 10 cm(-1). Only channels with maximum translational energy <150 cm(-1) are observed. The rotational excitation in the ammonia fragments is modest and can be fit by temperatures of 150 +/- 50 and 50 +/- 20 K for 1nu2 and 2nu2, respectively. The rotational distributions in the acetylene co-fragment pair-correlated with specific rovibrational states of ammonia appear statistical as well. The vibrational-state distributions, however, show distinct state specificity among channels with low translational energy release. The predominant channel is NH3(1nu2) + C2H2(2nu4 or 1nu4 + 1nu5), where nu4 and nu5 are the trans- and cis-bend vibrations of acetylene, respectively. A second observed channel, with much lower population, is NH3(2nu2) + C2H2(1nu4). No products are generated in which the ammonia is in the vibrational ground state or the asymmetric bend (1nu4) state, nor is acetylene ever generated in the ground vibrational state or with CC stretch excitation. The angular momentum (AM) model of McCaffery and Marsh is used to estimate impact parameters in the internal collisions that give rise to the observed rotational distributions. These calculations show that dissociation takes place from bent geometries, which can also explain the propensity to excite fragment bending levels. The low recoil velocities associated with the observed channels facilitate energy exchange in the exit channel, which results in statistical-like fragment rotational distributions.