Nonadiabatic time-dependent wave packet study of the D+ + H2 reaction system

A theoretical investigation on the nonadiabatic processes of the D(+) + H(2) reaction system has been carried out by means of exact three-dimensional nonadiabatic time-dependent wave packet calculations with an extended split operator scheme (XSOS). The diabatic potential energy surface newly constructed by Kamisaka et al. (J. Chem. Phys. 2002, 116, 654) was employed in the calculations. This study provided quantum cross sections for three competing channels of the reactive charge transfer, the nonreactive charge transfer, and the reactive noncharge transfer, which contrasted markedly to many previous quantum theoretical reports on the (DH(2))(+) system restricted to the total angular momentum J = 0. These quantum theoretical cross sections derived from the ground rovibrational state of H(2) show wiggling structures and an increasing trend for both the reactive charge transfer and the nonreactive charge transfer but a decreasing trend for the reactive noncharge transfer throughout the investigated collision energy range 1.7-2.5 eV. The results also show that the channel of the reactive noncharge transfer with the largest cross section is the dominant one. A further investigation of the v-dependent behavior of the probabilities for the three channels revealed an interesting dominant trend for the reactive charge transfer and the nonreactive charge transfer at vibrational excitation v = 4 of H(2). In addition, the comparison between the centrifugal sudden (CS) and exact calculations showed the importance of the Coriolis coupling for the reactive system. The computed quantum cross sections are also compared with the experimental measurement results.     

Nonadiabatic reactant-product decoupling calculation for the F(2P(1/2)) + H2 reaction

In this paper we present a theoretical study using time-dependent nonadiabatic reactant-product decoupling method for the state-to-state reactive scattering calculation of F((2)P(1/2))+H(2) (nu=j=0) reaction on the Alexander-Stark-Werner potential energy surface. In this nonadiabatic state-to-state calculation, the full wave function is partitioned into reactant component and a sum of all product components. The reactant and product components of the wave function are solved independently. For the excited state reaction, the state-to-state reaction probabilities for J=0.5 are calculated. Comparing the state-to-state reaction probabilities, it is found that the vibrational population of the HF product is dominated by vibrational levels nu=2 and 3. The rotation specific reaction probabilities of HF product in j=1 and 2 are larger than those in other rotational levels. As the rotation quantum number j increases, the positions of the peak in the rotational reaction probability of HF product in nu=3 shift to higher collision energy.     

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.     

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.     

Full-dimensional time-dependent wave packet dynamics of H2 + D2 reaction

Collision induced dissociation (CID), four center reaction (4C), and single exchange reaction (SE) in H(2) (v(1) = high) + D(2) (v(2) = low) were studied by means of time-dependent wave packet approach within a full-dimensional model. Initial state-selected total reaction probabilities for the three competitive processes have been computed on two realistic global potential energy surfaces of Aguado-Suárez-Paniagua and Boothroyd-Martin-Keogh-Peterson (BMKP) with the total angular momentum J = 0. The role of both vibrationally excited and rotationally excited reagents was examined by varying the initial vibrational and rotational states. The vibrational excitation of the hot diatom gives an enhancement effect on the CID process, while the vibrational excitation of the cold diatom gives an inhibition effect. The rotational excitation of both reagents has a significant effect on the reaction process. The 4C and SE probabilities are at least one order of magnitude smaller than the CID probabilities over the energy range considered. Isotope substitution effects were also studied by substituting the collider D(2) by H(2) and HD on the BMKP potential energy surfaces. The CID process is most efficient for the H(2) + D(2) combination and least efficient for the H(2) + H(2) combination and is different for the 4C and SE processes.     

Vibrational mode and collision energy effects on reaction of H2CO+ with C2D4

We report the effects of collision energy (Ecol) and five different H2CO+ vibrational modes on the reaction of H2CO+ with C2D4 over the center-of-mass E(col) range from 0.1 to 2.1 eV. Properties of various complexes and transition states were also examined computationally. Seven product channels are observed. Charge transfer (CT) has the largest cross section over the entire energy range, substantially exceeding the hard sphere cross section at high energies. Competing with CT are six channels involving transfer of one or more hydrogen atoms or protons and one involving formation of propanal, followed by hydrogen elimination. Despite the existence of multiple deep wells on the potential surface, all reactions go by direct mechanisms, except at the lowest collision energies, where short-lived complexes appear to be important. Statistical complex decay appears adequate to account for the product branching at low collision energies, however, even at the lowest energies, the vibrational effects are counter to statistical expectations. The pattern of Ecol and vibrational mode effects provide insight into factors that control reaction and interchannel competition.     

Coriolis coupling effects in the dynamics of deep well reactions: application to the H(+) + D2 reaction

We present exact and estimated quantum differential and integral cross sections as well as product state distributions for the title reaction. We employ a time-dependent wavepacket method including all Coriolis couplings and also an adapted code where the helicity quantum number and with this the Coriolis couplings have been truncated. Results from helicity truncated as well as helicity conserving (HC) calculation are presented. The HC calculations fail to reproduce the exact results due to the influence of the centrifugal barrier. While the truncated calculation overestimate the exact integral cross sections they reproduce the features of the integral cross section very well. We also find that the product rotational state distributions are well reproduced if the maximum helicity state is chosen carefully. The helicity truncated calculations fail to give a good approximation of differential cross sections.     

Coriolis coupling effects in the calculation of state-to-state integral and differential cross sections for the H+D2 reaction

The quantum wavepacket parallel computational code DIFFREALWAVE is used to calculate state-to-state integral and differential cross sections for the title reaction on the BKMP2 surface in the total energy range of 0.4-1.2 eV with D2 initially in its ground vibrational-rotational state. The role of Coriolis couplings in the state-to-state quantum calculations is examined in detail. Comparison of the results from calculations including the full Coriolis coupling and those using the centrifugal sudden approximation demonstrates that both the energy dependence and the angular dependence of the calculated cross sections are extremely sensitive to the Coriolis coupling, thus emphasizing the importance of including it correctly in an accurate state-to-state calculation.     

Nonadiabatic dynamics in the CH3+HCl-->CH4+Cl(2PJ) reaction

Nonadiabatic dynamics in the title reaction have been investigated by 2+1 REMPI detection of the Cl(2P(3/2)) and Cl*(2P(1/2)) products. Reaction was initiated by photodissociation of CH(3)I at 266 nm within a single expansion of a dilute mixture of CH(3)I and HCl in argon, giving a mean collision energy of 7800 cm(-1) in the center-of-mass frame. Significant production of Cl* was observed, with careful checks made to ensure that no additional photochemical or inelastic scattering sources of Cl* perturbed the measurements. The fraction of the total yield of Cl(2P(J)) atoms formed in the J=1/2 level at this collision energy was 0.150+/-0.024, and must arise from nonadiabatic dynamics because the ground potential energy surface correlates to CH(4)+Cl(2P(3/2)) products.     

Quantum approaches for the insertion dynamics of the H+ + D2 and D+ + H2 reactive collisions

The H(+)+D(2) and D(+)+H(2) reactive collisions are studied using a recently proposed adiabatic potential energy surface of spectroscopic accuracy. The dynamics is studied using an exact wave packet method on the adiabatic surface at energies below the curve crossing occurring at approximately 1.5 eV above the threshold. It is found that the reaction is very well described by a statistical quantum method for a zero total angular momentum (J) as compared with the exact ones, while for higher J some discrepancies are found. For J >0 different centrifugal sudden approximations are proposed and compared with the exact and statistical quantum treatments. The usual centrifugal sudden approach fails by considering too high reaction barriers and too low reaction probabilities. A new statistically modified centrifugal sudden approach is considered which corrects these two failures to a rather good extent. It is also found that an adiabatic approximation for the helicities provides results in very good agreement with the statistical method, placing the reaction barrier properly. However, both statistical and adiabatic centrifugal treatments overestimate the reaction probabilities. The reaction cross sections thus obtained with the new approaches are in rather good agreement with the exact results. In spite of these deficiencies, the quantum statistical method is well adapted for describing the insertion dynamics, and it is then used to evaluate the differential cross sections.     

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.     

Quantum wave packet study of the H+ + D2 reaction on diabatic potential energy surfaces

The exact three-dimensional nonadiabatic quantum dynamics calculations were carried out for the title reaction by a time-dependent wave packet approach based on a newly constructed diabatic potential energy surface (Kamisaka et al. J. Chem. Phys. 2002, 116, 654). Three processes including those of reactive charge transfer, nonreactive charge transfer, and reactive noncharge transfer were investigated to determine the initial state-resolved probabilities and reactive cross sections. The results show that a large number of resonances can be observed in the calculated probabilities due to the deep well on adiabatic ground surface and the dominant process is the reactive noncharge-transfer process. Some interesting dynamical features such as v-dependent and j-dependent behaviors of the probabilities are also revealed. In addition, a good agreement has been achieved in the comparison between the calculated quantum cross sections from the ground rovibrational initial state and the experimental measurement data.     

Critical re-examination of the reaction Cl + H2 ⇌ HCl + H

Recent data on the system Cl + H₂

HCl + H indicate that the rate constant ratio k₁/k₂ is a factor of two smaller than the equilibrium constant K. Two earlier explanations of this discrepancy are shown to be compatible with recent experimental and theoretical work. As an alternative explanation we suggest that flow tube measurements, ostensibly of k₂, actually determine 2k₂, due to Cl + H recombination on the walls of the flow tube. After correcting for this factor of two, all kinetic studies on this system are reconciled.     

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.     

Overtone spectroscopy of H2D+ and D2H+ using laser induced reactions

The method of laser induced reaction is used to obtain high-resolution IR spectra of H2D+ and D2H+ in collision with n-H2 at a nominal temperature of 17 K. For this purpose three cw-laser systems have been coupled to a 22-pole ion trap apparatus, two commercial diode laser systems in the ranges of 6100-6600 cm(-1) and 6760-7300 cm(-1), respectively, and a high-power optical parametric oscillator tunable in the range of 2600-3200 cm(-1). In total, 27 new overtone and combination transitions have been detected for H2D+ and D2H+, as well as a weak line in the nu1 vibrational band of H2D+ (2(20)<--1(01)) at 3164.118 cm(-1). The line positions are compared to high accuracy ab initio calculations, showing small but mode-dependent differences, being largest for three vibrational quanta in the nu2 symmetric bending of H2D+. Within the experimental accuracy, the relative values of the ab initio predicted Einstein B coefficients are confirmed.     

Theoretical study of the complex-forming CH + H2 --> CH2 + H reaction

The complex-forming CH + H2 --> CH2 + H reaction is studied employing a recently developed global potential energy function. The reaction probability in the total angular momentum J = 0 limit is estimated with a four-atom quantum wave packet method and compared with classical trajectory and statistical theory results. The formation of complexes from different reactant internal states is also determined with wave packet calculations. While there is no barrier to reaction along the minimum energy path, we find that there are angular constraints to complex formation. Trajectory-based estimates of the low-pressure rate constants are made and compared with experimental results. We find that zero-point energy violation in the trajectories is a particularly severe problem for this reaction.     

Collision theory treatment of the H + H2 exchange reaction

The rate constant of the three-dimensional H + H₂ reaction is calculated accurately in the temperature range 100–600 K. By adjusting the value of an effective collision diameter one obtains good agreement of the theory with experimental data. The precise corrections to the simplest approximate kinetic theories are also computed.     

A quantum wave packet dynamics study of the N(2D) + H2 reaction

We report a dynamics study of the reaction N((2)D) + H(2) (v=0, j=0-5) --> NH + H using the time-dependent quantum wave packet method and a recently reported single-sheeted double many-body expansion potential energy surface for NH(2)(1(2)A' ') which has been modeled from accurate ab initio multireference configuration-interaction calculations. The calculated probabilities for (v=0, j=0-5) are shown to display resonance structures, a feature also visible to some extent in the calculated total cross sections for (v=0, j=0). A comparison between the calculated centrifugal-sudden and coupled-channel reaction probabilities validate the former approximation for the title system. Rate constants calculated using a uniform J-shifting scheme and averaged over a Boltzmann distribution of rotational states are shown to be in good agreement with the available experimental values. Comparisons with other theoretical results are also made.