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.     

State-to-state reaction probabilities for the H+O2(v,j)-->O+OH(v',j') reaction on three potential energy surfaces

We report state-to-state and total reaction probabilities for J=0 and total reaction probabilities for J=2 and 4 for the title reaction, both for ground-state and initially rovibrationally excited reactants. The results for three different potential energy surfaces are compared and contrasted. The potential energy surfaces employed are the DMBE IV surface by Pastrana et al. [J. Phys. Chem. 94, 8073 (1990)], the surface by Troe and Ushakov (TU) [J. Chem. Phys. 115, 3621 (2001)], and the new XXZLG ab initio surface by Xu et al. [J. Chem. Phys. 122, 244305 (2005)]. Our results show that the total reaction probabilities from both the TU and XXZLG surfaces are much smaller in magnitude for collision energies above 1.2 eV compared to the DMBE IV surface. The three surfaces also show different behavior with regards to the effect of initial state excitation. The reactivity is increased on the XXZLG and the TU surfaces and decreased on the DMBE IV surface. Vibrational and rotational product state distributions for the XXZLG and the DMBE IV surface show different behaviors for both types of distributions. Our results show that for energies above 1.25 eV the dynamics on the DMBE IV surface are not statistical. However, there is also evidence that the dynamics on the XXZLG surface are not purely statistical for energies above the onset of the first excited product vibrational state v'=1. The magnitude of the total reaction probability is decreased for J>0 for the DMBE IV and the XXZLG surfaces for ground-state reactants. However, for initially rovibrationally excited reactants, the total reaction probability does not decrease as expected for both surfaces. As a result the total cross section averaged over all Boltzmann accessible rotational states may well be larger than the cross section reported in the literature for j=1.     

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.     

Differential and integral cross sections of the N(2D)+H2-->NH+H reaction from exact quantum and quasi-classical trajectory calculations

Exact quantum mechanical state-to-state differential and integral cross sections and their energy dependence have been calculated on an accurate NH2 potential energy surface (PES), using a newly proposed Chebyshev wave packet method. The NH product is found to have a monotonically decaying vibrational distribution and an inverted rotational distribution. The product angular distributions peak in both forward and backward directions, but with a backward bias. This backward bias is more pronounced than that observed previously on a less accurate PES. Both the differential and integral cross sections oscillate mildly with collision energy, indicating the dominance of short-lived resonances. The quantum mechanical results are compared with those obtained from quasi-classical trajectories. The agreement is generally reasonable and the discrepancies can be attributed to the neglect of quantum effects such as tunneling. Detailed analysis of the trajectories revealed that the backward bias in the differential cross section stems overwhelmingly from the fast insertion component of the reaction, augmented with some flux from the abstraction channel, particularly at high collision energies.     

Disagreement between theory and experiment in the simplest chemical reaction: collision energy dependent rotational distributions for H + D2 --> HD(nu' = 3,j') + D

We present experimental rotational distributions for the reaction H + D2 --> HD(nu' = 3,j') + D at eight different collision energies between 1.49 and 1.85 eV. We combine a previous measurement of the state-resolved excitation function for this reaction [Ayers et al., J. Chem. Phys. 119, 4662 (2003)] with the current data to produce a map of the relative reactive cross section as a function of both collision energy and rotational quantum number (an E-j' plot). To compare with the experimental data, we also present E-j' plots resulting from both time-dependent and time-independent quantum mechanical calculations carried out on the BKMP2 surface. The two calculations agree well with each other, but they produce rotational distributions significantly colder than the experiment, with the difference being more pronounced at higher collision energies. Disagreement between theory and experiment might be regarded as surprising considering the simplicity of this system; potential causes of this discrepancy are discussed.     

A time-dependent wave packet quantum scattering study of the reaction HD+ (v = 0 - 3;j0 = 1) + He --> HeH+(HeD+) + D(H)

Time-dependent wave packet quantum scattering (TWQS) calculations are presented for HD(+) (v = 0 - 3;j(0)=1) + He collisions in the center-of-mass collision energy (E(T)) range of 0.0-2.0 eV. The present TWQS approach accounts for Coriolis coupling and uses the ab initio potential energy surface of Palmieri et al. [Mol. Phys. 98, 1839 (2000)]. For a fixed total angular momentum J, the energy dependence of reaction probabilities exhibits quantum resonance structure. The resonances are more pronounced for low J values and for the HeH(+) + D channel than for the HeD(+) + H channel and are particularly prominent near threshold. The quantum effects are no longer discernable in the integral cross sections, which compare closely to quasiclassical trajectory calculations conducted on the same potential energy surface. The integral cross sections also compare well to recent state-selected experimental values over the same reactant and translational energy range. Classical impulsive dynamics and steric arguments can account for the significant isotope effect in favor of the deuteron transfer channel observed for HD(+)(v<3) and low translational energies. At higher reactant energies, angular momentum constraints favor the proton-transfer channel, and isotopic differences in the integral cross sections are no longer significant. The integral cross sections as well as the J dependence of partial cross sections exhibit a significant alignment effect in favor of collisions with the HD(+) rotational angular momentum vector perpendicular to the Jacobi R coordinate. This effect is most pronounced for the proton-transfer channel at low vibrational and translational energies.     

Quantum scattering LCAC-SW theory studies on reaction probabilities of three-dimensional H + H2 (v, j) → H2 (v′, j′) + H reaction

Upon the Liu, Siegbahn, Truhlar, Horowitz (LSTH) potential energy surface, the reaction probabilities of the three-dimensional (3-D) state-to-state H + H₂ (v, j) →H ₂(v′, j′) + H reaction are calculated with the linear combination of arrangement channels-scattering wavefunction (LCAC-SW) method. In the calculation, the vibration function of H₂ and the radial propagating wave functions are expanded by the real Gauss functions. The calculated threshold energy and the resonating structure are consistent with the results of the accurate quantum scattering calculations, which shows the accuration, simplicity and practicability of the LCAC-SW method. Project supported by the National Natural Science Fondation of China and the Doctoral Foundation of the State Education Commission of China.     

Time-dependent quantum mechanical wave packet study of the He+H(2) (+)(v,j)-->HeH(+)+H reaction

A detailed three-dimensional time-dependent quantum dynamical study of the He+H(2) (+)(v=0-3,j=0)-->HeH(+)+H reaction is reported for different vibrational v states of H(2) (+) in its ground rotational (j=0) state over a range of translational E(trans) energies on an accurate ab initio potential energy surface published by Palmieri et al. Plots of reaction probability as a function of total energy E reveal a large number of oscillations indicating the presence of a number of reactive scattering resonances. When averaged over total angular momentum J, some of the oscillations survive, indicating that they may be amenable to experimental observation. A comparison of our present results with our earlier results on the McLaughlin-Thompson-Joseph-Sathyamurthy surface and the experimental results from different research groups reveal a good deal of agreement as well as some discrepancies between theory and experiment at the level of state-selected gas phase dynamics.     

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.     

Inelastic scattering matrix elements for the nonadiabatic collision B(2P1/2)+H2(1Sigmag+,j)<-->B(2P3/2)+H2(1Sigmag+,j')

Inelastic scattering matrix elements for the nonadiabatic collision B(2P1/2)+H2(1Sigmag+,j)<-->B(2P3/2)+H2(1Sigmag+,j') are calculated using the time dependent channel packet method (CPM). The calculation employs 1 2A', 2 2A', and 1 2A" adiabatic electronic potential energy surfaces determined by numerical computation at the multireference configuration-interaction level [M. H. Alexander, J. Chem. Phys. 99, 6041 (1993)]. The 1 2A' and 2 2A', adiabatic electronic potential energy surfaces are transformed to yield diabatic electronic potential energy surfaces that, when combined with the total B+H2 rotational kinetic energy, yield a set of effective potential energy surfaces [M. H. Alexander et al., J. Chem. Phys. 103, 7956 (1995)]. Within the framework of the CPM, the number of effective potential energy surfaces used for the scattering matrix calculation is then determined by the size of the angular momentum basis used as a representation. Twenty basis vectors are employed for these calculations, and the corresponding effective potential energy surfaces are identified in the asymptotic limit by the H2 rotor quantum numbers j=0, 2, 4, 6 and B electronic states 2Pja, ja=1/2, 3/2. Scattering matrix elements are obtained from the Fourier transform of the correlation function between channel packets evolving in time on these effective potential energy surfaces. For these calculations the H2 bond length is constrained to a constant value of req=1.402 a.u. and state to state scattering matrix elements corresponding to a total angular momentum of J=1/2 are discussed for j=0<-->j'=0,2,4 and 2P1/2<-->2P1/2, 2P3/2 over a range of total energy between 0.0 and 0.01 a.u.     

Quantum calculations of the O(3P)+H2-->OH+H reaction

Quantum scattering calculations are reported for the O(3P)+H2(v=0,1) reaction using chemically accurate potential energy surfaces of 3A' and 3A" symmetry. We present state-to-state reaction cross sections and rate coefficients as well as thermal rate coefficients for the title reaction using accurate quantum calculations. Our calculations yield reaction cross sections that are in quantitative accord with results of recent crossed molecular beam experiments. Comparisons with results obtained using the J-shifting calculations show that the J-shifting approximation is quite reliable for this system. Thermal rate coefficients from the exact calculations and the J-shifting approximation agree remarkably well with experimental results. Our calculations also reproduce the markedly different OH(v'=0)/OH(v'=1) branching in O(3P)+H2(v=1) reaction, observed in experiments that use different O(3P) atom sources. In particular, we show that the branching ratio is a strong function of the kinetic energy of the O(3P) atom.     

Observation of partial wave structures in the integral cross section of the S((1)D2) + H2(j = 0) reaction

The integral cross section of the S((1)D(2)) + H(2)(j = 0) → SH + H reaction has been measured for the first time at collision energies from 0.820 down to 0.078 kJ mol(-1) in a high-resolution crossed beam experiment. The excitation function obtained exhibits a non-monotonic variation with collision energy and compares well with the results of high-level quantum calculations. In particular, the structures observed in the lower energy part, where only a few partial waves contribute, can be described in terms of the sequential opening of individual channels, consistent with the theoretical calculations.     

Coriolis coupling effects on the initial-state-resolved dynamics of the N(2D)+H2-->NH+H reaction

We present Coriolis coupling effects on the initial-state-resolved dynamics of the insertion reaction N((2)D)+H(2)(X (1)Sigma(g) (+))-->NH(X (3)Sigma(-) and a (1)Delta)+H((2)S), without and with nonadiabatic Renner-Teller (RT) interactions between the NH(2) X (2)B(1) and A (2)A(1) electronic states. We report coupled-channel (CC) Hamiltonian matrix elements, which take into account both Coriolis and RT couplings, use the real wave-packet and flux methods for calculating initial-state-resolved reaction probabilities, and contrast CC with centrifugal-sudden (CS) results. Without RT interactions, Coriolis effects are rather small up to J=40, and the CS approximation can be safely employed for calculating initial-state-resolved, integral cross sections. On the other hand, RT effects are associated with rather large Coriolis couplings, mainly near the linearity of NH(2), and the accuracy of the CS approximation thus breaks down at high collision energies, when the reaction starts on the excited A (2)A(1) surface. We also present the CC-RT distribution of the X (3)Sigma(-) and a (1)Delta electronic states of the NH products.     

Differential and integral cross sections for the H + O2 --> OH + O combustion reaction

We present accurate differential and integral cross sections for the H + O2 --> OH + O reaction obtained on a newly developed ab initio potential energy surface using time-independent and time-dependent quantum mechanical methods. The product angular distributions near the reaction threshold show pronounced forward and backward peaks, reflecting the complex-forming mechanism. However, the asymmetry of these peaks suggests certain nonstatistical behaviors, presumably due to some relatively short-lived resonances. The integral cross section increases monotonically with the collision energy above a reaction threshold.     

Renner-Teller quantum dynamics of the N(2D) + H2-->NH + H reaction

We present the Born-Oppenheimer (BO) and Renner-Teller (RT) quantum dynamics of the reaction (14)N((2)D)+(1)H(2)(X (1)Sigma(g) (+))-->NH(X (3)Sigma(-))+H((2)S), considering the NH(2) electronic states X (2)B(1) and A (2)A(1). These states correlate to the same (2)Pi(u) linear species, are coupled by RT nonadiabatic effects, and give NH(X (3)Sigma(-))+H and NH(a (1)Delta)+H, respectively. We develop the Hamiltonian matrix elements in the R embedding of the Jacobi coordinates and in the adiabatic electronic representation, using the permutation-inversion symmetry, and taking into account the nuclear-spin statistics. Collision observables are calculated via the real wave-packet (WP) and flux methods, using the potential-energy surfaces of Santoro et al. [J. Phys. Chem. A 106, 8276 (2002)]. WP snapshots show that the reaction proceeds via an insertion mechanism, and that the RT-WP avoids the A (2)A(1) potential barrier, jumping from the excited to the ground surface and giving mainly the NH(X (3)Sigma(-)) products. X (2)B(1) BO probabilities and cross sections show large tunnel effects and are approximately four to ten times larger than the A (2)A(1) ones. This implies a BO rate-constant ratio k(X (2)B(1))k(A (2)A(1)) approximately 10(5) at 300 K, i.e., a negligible BO formation of NH(a (1)Delta). When H(2) is rotationally excited, RT couplings reduce slightly the X (2)B(1) reaction observables, but enhance strongly the A (2)A(1) reactivity. These couplings are important at all collision energies, reduce the collision threshold, and increase remarkably reaction probabilities and cross sections. The RT k(A (2)A(1)) is thus approximately 3.3 order of magnitude larger than the BO value, and degeneracy-averaged, initial-state-resolved rate constants increase by approximately 13% and by approximately 47% at 300 and 500 K, respectively. Owing to an overestimation of the X (2)B(1) potential barrier, the calculated thermal rate is too low with respect to that observed, but we obtain a good agreement by shifting down the calculated cross section.     

Chemical reaction dynamics of Rydberg atoms with neutral molecules: a comparison of molecular-beam and classical trajectory results for the H(n)+D2-->HD+D(n') reaction

Recent molecular-beam experiments have probed the dynamics of the Rydberg-atom reaction, H(n)+D2-->HD+D(n) at low collision energies. It was discovered that the rotationally resolved product distribution was remarkably similar to a much more limited data set obtained at a single scattering angle for the ion-molecule reaction H++D2-->D++HD. The equivalence of these two problems would be consistent with the Fermi-independent-collider model (electron acting as a spectator) and would provide an important new avenue for the study of ion-molecule reactions. In this work, we employ a classical trajectory calculation on the ion-molecule reaction to facilitate a more extensive comparison between the two systems. The trajectory simulations tend to confirm the equivalence of the ion+molecule dynamics to that for the Rydberg-atom+molecule system. The theory reproduces the close relationship of the two experimental observations made previously. However, some differences between the Rydberg-atom experiments and the trajectory simulations are seen when comparisons are made to a broader data set. In particular, the angular distribution of the differential cross section exhibits more asymmetry in the experiment than in the theory. The potential breakdown of the classical model is discussed. The role of the "spectator" Rydberg electron is addressed and several crucial issues for future theoretical work are brought out.     

Differential cross sections and product energy distributions for the C(3P)+OH(X 2Pi)-->CO(X 1Sigma+)+H(2S) reaction using a quasiclassical trajectory method

Quasiclassical trajectory calculations have been carried out for the C((3)P)+OH(X (2)Pi)-->CO(X (1)Sigma(+))+H((2)S) reaction using a recent ab initio potential energy surface for the ground electronic state X (2)A(') of COH. Differential cross sections (DCSs), and product vibrational, rotational and translational distributions have been determined for a wide range of collision energies (0.001-1 eV). The role of excitations (rotation or vibration) of the OH reactant on these quantities has been investigated. Product vibrational, rotational, and translational distributions are found to be almost independent on the rovibrational state of OH, whereas DCSs show a weak dependence on the initial rotational state of OH. We also analyze the results using a study based on the lifetime of the intermediate complex and on the kinematic constraint associated with the mass combination.     

The effect of phonon modes on the H2(v, j)/D2(v, j)-Cu(1nn) scattering processes

We include the effect of the phonon modes originating from the three layers of Cu(1nn) surface atoms on the dynamics of incoming molecular [H(2)(v, j)/D(2)(v, j)] degrees of freedom (DOFs) through a mean-field approach, where the surface temperature is incorporated into the effective potential by considering Bose-Einstein probability (BEP) factor for the initial state distribution of the surface modes calculated within harmonic approximation. Such time and temperature dependent effective Hamiltonian is further subdivided assuming a weak coupling between the two sets of molecular DOFs, namely, (x, y, z, Z) and (X, Y), respectively, in particular, to reduce the computational cost and the corresponding coupled quantum dynamical equations of motion have been formulated in terms of Time Dependent Discrete Variable Representation (TDDVR) approach. We demonstrate the workability of TDDVR method to investigate the scattering of H(2)(v, j) on Cu(1nn) surface by calculating the reaction probabilities and scattering cross-sections. Calculated results show that the phonon modes affect (a) the state-to-state transition probabilities of the scattered H(2) molecule substantially but chemisorption and physisorption processes negligibly and (b) the reaction probability of the incoming D(2) molecule noticeably.     

Exact state-to-state quantum dynamics of the F + HD --> HF(v' = 2) + D reaction on model potential energy surfaces

In this paper, we present the results of a theoretical investigation on the dynamics of the title reaction at collision energies below 1.2 kcal/mol using rigorous quantum reactive scattering calculations. Vibrationally resolved integral and differential cross sections, as well as product rotational distributions, have been calculated using two electronically adiabatic potential energy surfaces, developed by us on the basis of semiempirical modifications of the entrance channel. In particular, we focus our attention on the role of the exothermicity and of the exit channel region of the interaction on the experimental observables. From the comparison between the theoretical results, insight about the main mechanisms governing the reaction is extracted, especially regarding the bimodal structure of the HF(v = 2) nascent rotational state distributions. A good overall agreement with molecular beam scattering experiments has been obtained.