Dynamics of the Rydberg electron in H*+D2-->D*+HD reactive collisions

Experimental crossed-beam studies carried out previously have indicated that the dynamics of the Rydberg-atom-molecule reaction H*+D2-->D*+HD are very similar to those of the corresponding ion-molecule reaction H++D2-->D++HD. The equivalence of the cross sections for these related systems would open up a new approach to the experimental study of ion-molecule reactions. However, a recent experimental and theoretical study has brought to light some important qualitative differences between the Rydberg-atom reaction and the ion-molecule reaction; in particular, the experimental cross section for the Rydberg-atom reaction exhibits a higher degree of forward-backward scattering asymmetry than predicted by a quasiclassical trajectory study of the ion-molecule reaction. In this paper, the authors consider the dynamics of the Rydberg-electron over the course of a reactive collision and the implications of these dynamics for the Rydberg-atom-molecule crossed-beam experiment. Using an approach based on perturbation theory, they estimate the attenuation of the experimental signal due to the Rydberg-electron dynamics as a function of the scattering angle. They show that at least part of the experimental asymmetry can be ascribed to this angle dependent attenuation. Their results offer general insight into the practical aspects of the experimental study of ion-molecule reactions by means of their Rydberg-atom counterparts.     

Angular distributions of H-induced HD and D2 desorptions from the Si(100) surfaces

We measured angular distributions of HD and D2 molecules desorbed via the reactions H+DSi(100)-->HD [abstraction (ABS)] and H+DSi(100)-->D2 [adsorption-induced-desorption (AID)], respectively. It was found that the angular distribution of HD molecules desorbed along ABS is broader than that of D2 molecules desorbed along AID, i.e., the former could be fit with cos(2.0+/-0.2) theta, while the latter with cos(5.0+/-0.5) theta. This difference of the angular distributions between the two reaction paths suggests that their dynamic mechanisms are different. The observed cos2 theta distribution for the ABS reaction was reproduced by the classical trajectory calculations over the London-Eyring-Polanyi-Sato potential-energy surfaces. The simulation suggests that the HD desorption along the ABS path takes place along the direction of Si-D bonds, but the apparent angular distribution is comprised of multiple components reflecting the different orientations of D-occupied Si dimers in the (2 x 1) and (1 x 2) double domain structures.     

On the dynamics of the H+ +D2(v=0,j=0)-->HD+D + reaction: a comparison between theory and experiment

The H+ +D2(v=0,j=0)-->HD+D + reaction has been theoretically investigated by means of a time independent exact quantum mechanical approach, a quantum wave packet calculation within an adiabatic centrifugal sudden approximation, a statistical quantum model, and a quasiclassical trajectory calculation. Besides reaction probabilities as a function of collision energy at different values of the total angular momentum, J, special emphasis has been made at two specific collision energies, 0.1 and 0.524 eV. The occurrence of distinctive dynamical behavior at these two energies is analyzed in some detail. An extensive comparison with previous experimental measurements on the Rydberg H atom with D2 molecules has been carried out at the higher collision energy. In particular, the present theoretical results have been employed to perform simulations of the experimental kinetic energy spectra.     

New, unexpected, and dominant mechanisms in the hydrogen exchange reaction

A quasiclassical trajectory study of the state specific H+D(2)(upsilon = 0,j = 0) --> HD(upsilon' = 0,j' = 0) + D reaction at a collision energy of 1.85 eV (total energy of 2.04 eV) found that the scattering is governed by two unexpected and dominant new mechanisms, and not by direct recoil as is generally assumed. The new mechanisms involve strong interaction with the sloping potential around the conical intersection, an area of the potential energy surface not previously considered to have much effect upon reactive scattering. Initial investigations indicate that more than 50% of reactive scattering could be the result of these new mechanisms at this collision energy. Features in the corresponding quantum mechanical results can be attributed to these new (classical) reaction mechanisms.     

Nonadiabatic effects in the H+D2 reaction

The state-to-state dynamics of the H+D2 reaction is studied by the reactant-product decoupling method using the double many-body expansion potential energy surface. Two approaches are compared: one uses only the lowest adiabatic sheet while the other employs both coupled diabatic sheets. Rotational distributions for the reaction H+D2 (upsilon = 0, j = 0)-->HD(upsilon' = 3, j')+D are obtained at eight different collision energies between 1.49 and 1.85 eV; no significant difference are found between the two approaches. Initial state-selected total reaction probabilities and integral cross sections are also given for energies ranging from 0.25 up to 2.0 eV with extremely small differences being observed between the two sets of results, thus showing that the nonadiabatic effects in the title reaction are negligible at least for small energies below 2.0 eV.     

Time dependent quantum dynamics study of the O++H2(v=0,j=0)-->OH++H ion-molecule reaction and isotopic variants (D2,HD)

The time dependent real wave packet method using the helicity decoupling approximation was used to calculate the cross section evolution with collision energy (excitation function) of the O++H2(v=0,j=0)-->OH++H reaction and its isotopic variants with D2 and HD, using the best available ab initio analytical potential energy surface. The comparison of the calculated excitation functions with exact quantum results and experimental data showed that the present quantum dynamics approach is a very useful tool for the study of the selected and related systems, in a quite wide collision energy interval (approximately 0.0-1.1 eV), involving a much lower computational cost than the quantum exact methods and without a significant loss of accuracy in the cross sections.     

Quasiclassical trajectory study of the C(1D) + H2 reaction and isotopomeric variants: kinetic isotope effect and CD/CH branching ratio

The recently proposed ab initio single-sheeted double many-body expansion potential energy for the methylene molecule has been used to perform quasiclassical trajectory (QCT) calculations for the title reaction. Thermal and initial state-specific (v = 0, j = 0) rate constants for the C((1)D) + H(2)/HD/D(2) reactions have been obtained over a wide range of temperatures. Cross sections for the reaction C((1)D) + H(2) and its deuterated isotopes have also been calculated, as well as the CD/CH branching ratios for the C((1)D) + HD reaction. It is found that the CD + H product channel in the C((1)D) + HD reaction is preferred relative to the CH + D channel. The estimated rate constants are predicted to be in the order k(H2) > k(HD) > k(D2) and the calculated cross sections and rate constants compared with available theoretical and experimental data.     

Macroscopic evidences for non-Rice-Ramsperger-Kassel effects in the reaction between H3O+ and D2O: the occurrence of nonstatistical isotopic branching ratio

The dynamics of the isotopic scrambling in the energized and metastable complex D2O-H3O+ has been studied using classical molecular dynamics (MD) trajectories starting from regions of phase space corresponding to an already formed collisional complex. The simulations cover the range of internal energies spanned by gas phase collision experiments. Rate constants for the isotopic exchange and the complex dissociation have been computed; the isotopic branching ratio R=[HD2O+]/[H2DO+] has also been obtained from MD simulations and found to deviate substantially from an equivalent prediction based on a previously proposed kinetic scheme. This finding suggests the possibility that details of the reaction dynamics play a role in defining the isotopic branching ratio. The analysis of trajectory results indicated a relatively long lifetime for the collisional complex and the presence of multiple time scales for the exchange process, with a large fraction of the exchange events being separated only by a single oxygen-oxygen vibration or half of it. The occurrence of these fast consecutive jumps and their different probabilities as a function of the relative direction between first and second jumps suggest the presence of ballistic motion in the complex following each reactive event. This can be explained on the basis of overlapping regions in phase space and it is used to provide an explanation of the difference between kinetic and MD branching ratios.     

Thermochemistry and accurate quantum reaction rate calculations for H2/HD/D2 + CH3

Accurate quantum-mechanical results for thermodynamic data, cumulative reaction probabilities (for J = 0), thermal rate constants, and kinetic isotope effects for the three isotopic reactions H2 + CH3 --> CH4 + H, HD + CH3 --> CH4 + D, and D2 + CH3 --> CH(3)D + D are presented. The calculations are performed using flux correlation functions and the multiconfigurational time-dependent Hartree (MCTDH) method to propagate wave packets employing a Shephard interpolated potential energy surface based on high-level ab initio calculations. The calculated exothermicity for the H2 + CH3 --> CH4 + H reaction agrees to within 0.2 kcal/mol with experimentally deduced values. For the H2 + CH3 --> CH4 + H and D2 + CH3 --> CH(3)D + D reactions, experimental rate constants from several groups are available. In comparing to these, we typically find agreement to within a factor of 2 or better. The kinetic isotope effect for the rate of the H2 + CH3 --> CH4 + H reaction compared to those for the HD + CH3 --> CH4 + D and D2 + CH3 --> CH(3)D + D reactions agree with experimental results to within 25% for all data points. Transition state theory is found to predict the kinetic isotope effect accurately when the mass of the transferred atom is unchanged. On the other hand, if the mass of the transferred atom differs between the isotopic reactions, transition state theory fails in the low-temperature regime (T < 400 K), due to the neglect of the tunneling effect.     

Cross sections of the O+ + H2 --> OH+ + H ion-molecule reaction and isotopic variants (D2, HD): quasiclassical trajectory study and comparison with experiments

A dynamics study [cross section and microscopic mechanism versus collision energy (E(T))] of the reaction O+ + H2 --> OH+ + H, which plays an important role in Earth's ionosphere and interstellar chemistry, was conducted using the quasiclassical trajectory method, employing an analytical potential energy surface (PES) recently derived by our group [R. Martinez et al., J. Chem. Phys. 120, 4705 (2004)]. Experimental excitation functions for the title reaction, as well as its isotopic variants with D2 and HD, were near-quantitatively reproduced in the calculations in the very broad collision energy range explored (E(T) = 0.01-6.0 eV). Intramolecular and intermolecular isotopic effects were also examined, yielding data in good agreement with experimental results. The reaction occurs via two microscopic mechanisms (direct and nondirect abstraction). The results were satisfactorily interpreted based on the reaction probability and the maximum impact parameter dependences with E(T), and considering the influence of the collinear [OHH]+ absolute minimum of the PES on the evolution from reactants to products. The agreement between theory and experiment suggests that the reaction mainly occurs through the lowest energy PES and nonadiabatic processes are not very important in the wide collision energy range analyzed. Hence, the PES used to describe this reaction is suitable for both kinetics and dynamics studies.     

Tunneling chemical reactions D + H2 --> DH + H and D + DH --> D2 + H in solid D2-H2 and HD-H2 mixtures: an electron-spin-resonance study

Tunneling chemical reactions D + H2 --> DH + H and D + DH --> D2 + H in solid HD-H2 and D2-H2 mixtures were studied in the temperature range between 4 and 8 K. These reactions were initiated by UV photolysis of DI molecules doped in these solids for 30 s and followed by measuring the time course of electron-spin-resonance (ESR) intensities of D and H atoms. ESR intensity of D atoms produced by the photolysis decreases but that of H atoms increases with time. Time course of the D and H intensities has the fast and slow processes. The fast process, which finishes within approximately 300 s after the photolysis, is assigned to the reaction of D atom with one of its nearest-neighboring H2 molecules, D(H2)n(HD)(12-n) --> H(H2)(n-1)(HD)(13-n) or D(H2)n(D2)(12-n) --> H(HD)(H2)(n-1)(D2)(12-n) for 12 > or = n > or = 1. Rate constant for the D + H2 reaction between neighboring D atom-H2 molecule pair is determined to be (7.5 +/- 0.7) x 10(-3) s(-1) in solid HD-H2 and (1.3+/-0.3) x 10(-2) s(-1) in D2-H2 at 4.1 K, which is very close to that calculated based on the theory of chemical reaction in gas phase by Hancock et al. [J. Chem. Phys. 91, 3492 (1989)] and Takayanagi and Sato [J. Chem. Phys. 92, 2862 (1990)]. This rate constant was found to be independent of temperature up to 7 K within experimental error of +/-30%. The slow process is assigned to the reaction of D atom produced in a cage fully surrounded by HD or D2 molecules, D(HD)12 or D(D2)12. This D atom undergoes the D + DH reaction with one of its nearest-neighboring HD molecules in solid HD-H2 or diffuses to the neighbor of H2 molecules to allow the D + H2 reaction in solid HD-H2 and D2-H2. The former is the main channel in solid HD-H2 below 6 K where D atoms diffuse very slowly, whereas the latter dominates over the former above 6 K. Rate for the reactions in the slow process is independent of temperature below 6 K but increases with the increase in temperature above 6 K. We found that the increase is due to the increase in hopping rate of D atoms to the neighbor of H2 molecules. Rate constant for the D + DH reaction was found to be independent of temperature up to 7 K as well.     

Effect of Coriolis coupling in chemical reaction dynamics

It is essential to evaluate the role of Coriolis coupling effect in molecular reaction dynamics. Here we consider Coriolis coupling effect in quantum reactive scattering calculations in the context of both adiabaticity and nonadiabaticity, with particular emphasis on examining the role of Coriolis coupling effect in reaction dynamics of triatomic molecular systems. We present the results of our own calculations by the time-dependent quantum wave packet approach for H + D2 and F(2P3/2,2P1/2) + H2 as well as for the ion-molecule collisions of He + H2 +, D(-) + H2, H(-) + D2, and D+ + H2, after reviewing in detail other related research efforts on this issue.     

Theoretical study of the reaction mechanism of HCN+ and CH4 of relevance to Titan's ion chemistry

Titan is the largest satellite of Saturn. In its atmosphere, CH4 is the most abundant neutral after nitrogen. In this paper, the complex doublet potential-energy surface related to the reaction between HCN+ and CH4 is investigated at the B3LYP/6-311G(d,p), CCSD(T)/6-311G++(3df,2pd)(single-point), and QCISD/6-311G(d,p) computational levels. A total of seven products are located on the PES. The initial association of HCN+ with CH4 is found to be a prereaction complex 1 (HCNHCH3(+)) without barrier. Starting from 1, the most feasible pathway is the direct H-abstraction process (the internal C-H bond dissociation) leading to the product P1 (HCNH++CH3). By C-C addition, prereaction complex 1 can form intermediate 2 (HNCHCH3(+)) and then lead to the product P2 (CH3CNH++H). The rate-controlling step of this process is only 25.6 kcal/mol. It makes the Path P2 (1) R --> 1 --> TS1/2 --> 2 --> TS2/P2 --> P2 another possible way for the reaction. P3 (HCNCH3(+) + H), P5 (cNCHCH2(+) + H2), and P6 (NCCH3(+) + H2) are exothermic products, but they have higher barriers (more than 40.0 kcal/mol); P4 (H + HCN + CH3(+)) and P7 (H + H2 + HCCNH+) are endothermic products. They should be discovered under different experimental or interstellar conditions. The present study may be helpful for investigating the analogous ion-molecule reaction in Titan's atmosphere.     

Influence of rotation and isotope effects on the dynamics of the N((2)D)+H(2) reactive system and of its deuterated variants

Integral cross sections and thermal rate constants have been calculated for the N((2)D)+H(2) reaction and its isotopic variants N((2)D)+D(2) and the two-channel N((2)D)+HD by means of quasiclassical trajectory and statistical quantum-mechanical model methods on the latest ab initio potential-energy surface [T.-S. Ho et al., J. Chem. Phys. 119, 3063 (2003)]. The effect of rotational excitation of the diatom on the dynamics of these reactions has been investigated and interesting discrepancies between the classical and statistical model calculations have been found. Whereas a net effect of reagent rotation on reactivity is always observed in the classical calculations, only a very slight effect is observed in the case of the asymmetric N((2)D)+HD reaction for the statistical quantum-mechanical method. The thermal rate constants calculated on this Potential-Energy Surface using quasiclassical trajectory and statistical model methods are in good agreement with the experimental determinations, although the latter are somewhat larger. A reevaluation of the collinear barrier of the potential surface used in the present study seems timely. Further theoretical and experimental studies are needed for a full understanding of the dynamics of the title reaction.     

An experimental and quasiclassical trajectory study of the rovibrationally state-selected reactions: HD+(v=0-15,j=1)+He-->HeH+(HeD+)+D(H)

The absolute integral cross sections for the formation of HeH+ and HeD+ from the collisions of HD+(v,j=1)+He have been examined over a broad range of vibrational energy levels v=0-13 at the center-of-mass collision energies (ET) of 0.6 and 1.4 eV using the vacuum ultraviolet (VUV) pulsed field ionization photoelectron secondary ion coincidence method. The ET dependencies of the integral cross sections for products HeH+ and HeD+ from HD+(v=0-4)+He collisions in the ET range of 0-3 eV have also been measured using the VUV photoionization guided ion beam mass spectrometric technique, in which vibrationally selected HD+(v) reactant ions were prepared via excitation of selected autoionization resonances of HD. At low total energies, a pronounced isotope effect is observed in absolute integral cross sections for the HeH++D and HeD++H channels with significant favoring of the deuteron transfer channel. As v is increased in the range of v=0-9, the integral cross sections of the HeH++D channel are found to approach those of HeD++H. The observed velocity distributions of products HeD+ and HeH+ are consistent with an impulsive or spectator-stripping mechanism. Detailed quasiclassical trajectory (QCT) calculations are also presented for HD+(v,j=1)+He collisions at the same energies of the experiment. The QCT calculations were performed on the most accurate ab initio potential energy surface available. If the zero-point energy of the reaction products is taken into account, the QCT cross sections for products HeH+ and HeD+ from HD+(v)+He are found to be significantly lower than the experimental results at ET values near the reaction thresholds. The agreement between the experimental and QCT cross sections improves with translational energy. Except for prethreshold reactivity, QCT calculations ignoring the zero-point energy in the products are generally in good agreement with experimental absolute cross sections. The experimental HeH+/HeD+ branching ratios for the HD+(v=0-9)+He collisions are generally consistent with QCT predictions. The observed isotope effects can be rationalized on the basis of differences in thermochemical thresholds and angular momentum conservation constraints.     

Direct dynamics trajectory study of vibrational effects: can polanyi rules be generalized to a polyatomic system?

We present an ab initio direct dynamics trajectory study of the hydrogen abstraction reaction: H2CO+ + CD4 --> H2COD+ + CD3, with methane excited in two different distortion modes (nu4 and nu2). The trajectory simulations were able to reproduce experimental results and for the first time show how vibrational enhancement originates in reaction of small polyatomic species. Roughly equal contributions from two vibrational enhancement mechanisms were found. The "distortion" mechanism correlates the vibrational effects with vibration-induced reactant distortions, and the "velocity" mechanism correlates vibrational effects with vibrational velocities of the constituent atoms. This reaction has a reactant-like transition state and, thus, would correspond to an "early" barrier system in the context of the well-known Polanyi rules for predicting effects of vibration and collision energy. Straightforward application of these rules would predict that vibration should be ineffective in driving reaction, in disagreement with both experiment and trajectory results. Using the trajectories for guidance, we were able to construct a two-dimensional cut through the reaction potential energy surface that does suggest a predictive, Polanyi-type rule.     

Quantum-classical wave packet calculations for the O(1D) + N2O → NO + NO/N2 + O2 reaction

Mixed quantum-classical calculations have been carried out for the O(¹D) + N₂O reaction with an emphasis on the effect of the relative translational energy as well as initial vibrational state of N₂O on the NO + NO/N₂ + O₂ product branching. The calculations were done within a planar constraint using a five-dimensional analytical potential energy surface previously developed by our group. Three vibrational coordinates in the N₂O molecule were treated with a quantum wave packet technique while other two degrees of freedom, translational and angular motions of O(¹D) with respect to N₂O, were described with classical mechanics. We have found that the initial orientation angle significantly affects the NO + NO/N₂ + O₂ product branching similar to our previous classical trajectory result using the same potential surface. It has been also found that the branching ratio decreases as the translational energy increases except for a low energy region. Excitation of the initial vibrational state of the N₂O reactant does not largely affect the reaction dynamics.     

Ab initio analytical potential energy surface and quasiclassical trajectory study of the O(+)((4)S)+H(2)(X (1)Sigma(g) (+))-->OH(+)(X (3)Sigma(-))+H((2)S) reaction and isotopic variants

An analytical potential energy surface (PES) representation of the O(+)((4)S)+H(2)(X (1)Sigma(g) (+)) system was developed by fitting around 600 CCSD(T)/cc-pVQZ ab initio points. Rate constant calculations for this reaction and its isotopic variants (D(2) and HD) were performed using the quasiclassical trajectory (QCT) method, obtaining a good agreement with experimental data. Calculations conducted to determine the cross section of the title reaction, considering collision energies (E(T)) below 0.3 eV, also led to good accord with experiments. This PES appears to be suitable for kinetics and dynamics studies. Moreover, the QCT results show that, although the hypotheses of a widely used capture model are not satisfied, the resulting expression for the cross section can be applied within a suitable E(T) interval, due to errors cancellation. This could be a general situation regarding the application of this simple model to ion-molecule processes.     

Ab initio molecular dynamics study of the hydrolysis reaction of diborane

The hydrolysis reaction of the diborane molecule in aqueous solution has been studied by a series of Car-Parrinello Molecular Dynamics simulations in the Blue Moon Ensemble. The total reaction has been divided into two parts: one dealing with the breaking of B(2)H(6) molecule and the formation of a BH(4)(-) ion, a H(2)BOH molecule and a H(+) ion; the second leads to the formation of two hydrogen molecules and another H(2)BOH molecule, starting from BH(4)(-), two water molecules and a H(+) ion. The total reaction studied in this work has been B(2)H(6) + 2H(2)O --> 2H(2)BOH + 2H(2). We have described both structurally and electronically the reagents and the products through the radial distribution functions and the Wannier Function Center positions calculations, with attention to the solvent effects on the compounds. The free energy barrier value for the first part of the reaction and a detailed mechanisms for both parts have been reported. An interesting behavior of BH(3) and H(2) molecules in solution has been observed. They form a quite stable three center bond between the electron pair of the hydrogen molecule and the empty orbital of the boron atom in BH(3), which has been described from both a structural and electronic point of view.