The O + O2 reaction: quantum detailed probabilities and thermal rate coefficients

The detailed quantum probabilities of the O + O₂ reactive system have been computed at zero total angular momentum using the time-independent quantum program ABC thanks to the restructuring of the code and its implementation on the EGEE production Grid. Their main features are discussed and out of them J-shifting thermal rate coefficients have been computed to compare with the experiment and quasiclassical trajectory results over a wide range of temperatures.     

Constraints at the transition state of the D + H2 reaction: quantum bottlenecks vs. stereodynamics

This article presents a quasiclassical trajectory method for the calculation of cumulative reaction probabilities by sampling of the helicity quantum number of the reagents (k). The method is applied to the D + H(2) reaction at various total angular momentum (J) values, and the helicity-resolved quasiclassical cumulative reaction probabilities are compared to their quantum mechanical counterparts. The agreement between the two sets of results is fairly good. In particular, k-dependent, J-independent reaction thresholds found with quantum methods are reproduced by the quasiclassical calculations. The shift of these thresholds with increasing k, which has been previously attributed to the quantum bottleneck states taking part in the reaction, is revisited and discussed also in terms of the reaction stereodynamics.     

A state-to-state dynamical study of the Br + H(2) reaction: comparison of quantum and classical trajectory results

We present a detailed theoretical investigation of the dynamics corresponding to the strongly endothermic Br + H(2) (v = 0-1, j = 0) → H + HBr reaction in the 0.85 to 1.9 eV total energy range. State-averaged and state-to-state results obtained through time-independent wave packet (TIWP) and time-independent quantum mechanical (TIQM) calculations and quasiclassical trajectories (QCT) are compared and analyzed. The agreement in the results obtained with both quantum mechanical results is very good overall. However, although QCT calculations reproduce the general features, their agreement with the QM results is sometimes only qualitative. The analysis of the mechanism based on state-averaged results turns out to be deceptive and conveys an oversimplified picture of the reaction consistent with a direct-rebound mechanism. Consideration of state-to-state processes, in contrast, unveils the existence of multiple mechanisms that give rise to a succession of maxima in the differential cross section (DCS). Such mechanisms correlate with different sets of partial waves and display similar collision times when analyzed through the time-dependent DCS.     

Quantum dynamics of the H + O2 --> O + OH reaction on an accurate ab initio potential energy surface

We report exact time-dependent and time-independent quantum mechanical studies of the title reaction on an accurate ab initio potential energy surface of Xu et al. (J. Chem. Phys. 2005, 122, 24305). The J = 0 reaction probabilities for several reactant states show sharp resonance structures superimposed on relatively low backgrounds, and they are remarkably different from existing quantum results on an earlier potential energy surface (DMBE-IV). The new findings reported here suggest that our current understanding of this important reaction might require significant revision.     

Competition between dissociation and exchange processes in a collinear A + BC collision: Comparison of quantum and classical results

Previous calculations on a reactive/dissociative H + HD model system have been extended to higher collision energies. Exact quantum dissociation probabilities are now available in the 1–10 eV total energy range for initial vibrational quantum numbers v = 0,1 and 2. Comparison with quasiclassical results shows the absence of quantum tails at dissociation threshold, but large quantum effects at higher energies.     

Isotopic replacement in ionic systems: the 4He2+ + 3He-->3He 4He+ + 4He reaction

Full quantum dynamics calculations have been carried out for the ionic reaction (4)He(2) (+)+(3)He and state-to-state reactive probabilities have been obtained using both time-dependent and time-independent approaches. An accurate ab initio potential-energy surface has been employed for the present quantum dynamics and the two sets of results are shown to be in agreement with each other. The results for zero total angular momentum suggest a marked presence of atom exchange (isotopic replacement) reaction with probabilities as high as 60%. The reaction probabilities are only weakly dependent on the initial vibrational state of the reactants, while they are slightly more sensitive to the degree of rotational excitation. A brief discussion of the results for selected higher total angular momentum values is also presented, while the l-shifting approximation [S. K. Gray et al., Phys. Chem. Chem. Phys. 1, 1141 (1999)] has been used to provide estimates of the total reaction rates for the title process. Such rates are found to be large enough to possibly become experimentally accessible.     

Quantum mechanical coupled-channel collision-induced dissociation calculations with hyperspherical coordinates

A time-independent coupled-channel method, using hyperspherical coordinates, has been developed for calculating quantum mechanical collision-induced dissociation probabilities for collinear atom-diatom systems in which the exchange reaction can also occur. The results for a model potential energy surface are compared with quasi-classical trajectory calculations and discussed.     

Wave packet and quasiclassical trajectory calculations for the N(D) + H2 reaction and its isotopic variants

The N(²D) + H₂(v = 0, j = 0) reaction and its HD and D₂ isotopic variants have been studied by means of quantum mechanical real wave packet and wave packet with split operator and quasiclassical trajectory methodologies on the potential energy surface of Ho et al. [J. Chem. Phys. 119 (2003) 6]. Total initial state-selected and final state-resolved reaction probabilities and product rotational distributions have been calculated for total angular momentum J = 0 in a broad range of collision energies. The real wave packet results are in very good agreement with the corresponding split operator wave packet calculations. A reasonable overall good agreement has been found between the wave packet and quasiclassical trajectory results. Integral cross-sections and thermal rate constants have been calculated from the wave packet reaction probabilities by means of standard J-shifting, refined J-shifting and uniform J-shifting methods in combination with the centrifugal sudden approximation for J > 0. Comparisons with available exact wave packet, quasiclassical trajectory and experimental results are made and discussed.     

Time-dependent reactive scattering for the system H- + D2 <--> HD + D- and comparison with H- + H2 <--> H2 + H-

This work presents results of quantum mechanical calculations of reaction probabilities for the ion-neutral molecule collisions H- + D2 <--> HD + D-. Time-dependent wave packet propagations for total angular momentum J not equal to 0, including the full Coriolis coupling, are performed. The calculated state-to-state reaction probabilities using product Jacobi coordinates are compared with energy-resolved reaction probabilities calculated with the flux-operator using reactant Jacobi coordinates and with time-independent calculations. Differences between nearly converged integral cross sections and those using the J-shifting method and centrifugal sudden approximation and comparison with experimental results will be presented.     

Quasiclassical determination of reaction probabilities as a function of the total angular momentum

This article presents a quasiclassical trajectory (QCT) method to determine the reaction probability as a function of the total angular momentum J for any given value of the initial rotational angular momentum j. The proposed method is based on a discrete sampling of the total and orbital angular momenta for each trajectory and on the development of equations that have a clear counterpart in the quantum-mechanical (QM) case. The reliability of the method is illustrated by comparing QCT and time-dependent wave-packet QM results for the H+D(2)(upsilon=0,j=4,10) reaction. The small discrepancies between both sets of calculations, when they exist, indicate some genuine quantum effects. In addition, a procedure to extract the reaction probabilities as a function of J when trajectories are calculated in the usual way using a continuous distribution of impact parameters is also described.     

Quantum mechanical scattering calculations on a spline-fitted ab initio surface: the He + H+2 (ν = 0, 1, 2) → HeH+ + H reaction

All-channel time-dependent quantum mechanical reaction probabilities are reported for the collinear He + H⁺₂(ν = 0, 1, 2) → HeH⁺ + H reaction at a total energy of 1.2 eV on previously reported diatomics-in-molecule (DIM) and spline fitted ab initio (SAI) surfaces. These results are in agreement with the previous quasiclassical trajectory results in that there is vibrational enhancement of the reaction probability on the SAI surface but not on the DIM surface. This agreement lends support to our previously drawn conclusion that small differences in the potential-energy surface can lead to substantially different dynamic results.     

S(1D)+D2反应的立体动力学理论研究

在ab initio 势能面上利用非含时量子反应散射理论研究了S(1D)+D2反应在散射能为22.18 kJ·mol-1时的立体动力学. 计算的微分反应截面呈现弱反对称的前向后向散射, 这一结果与以前的准经典结果和实验结果相一致. 态分辨的极化依赖的微分反应截面和极化参数显示产物极化并非各向异性分布, 这意味着在反应过程中由于势能面上存在一个深势阱从而形成了长寿命的复合物.     

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.     

A comparison of attack angle dependence of exchange channel of reaction H + HS (v = 0, 1; j = 0) on 3A″ and 3A′ surfaces

In this theoretical work, we report quasiclassical dynamics predictions for the attack angle‐dependence exchange processes for the H + HS (v = 0, 1; j = 0) reaction by using the new triplet 3A″ and 3A′ potential energy surfaces, respectively. The calculated quasiclassical reaction probabilities of exchange reaction channel of reaction H(D)′ + H(D)S for J = 0, 10, 20, 30, 40 are in good agreement with quantum wave packet results over the collision energy range from 0.1 to 2.0 eV on 3A″ surfaces. The attack angle dependence reaction probability of the title reactions at J = 0 are calculated, respectively, on the two surfaces. The reaction probability was found to be strongly dependent on the attack angle. It may be ascribe to the significant difference of the effective potential barrier height in the two reactions. Besides, the reaction probabilities of exchange reaction channel of reaction H(D)′ + H(D)S for J = 0, 10, 20, 30, 40 are also predicted on 3A′ surfaces. © 2013 Wiley Periodicals, Inc.     

Exact quantum dynamics of N(2D) + H2 --> NH + H reaction: cross-sections, rate constants, and dependence on reactant rotation

Using an exact Chebyshev wave packet method, initial state-specified (upsilon(i)=0, j(i)=0,2) integral cross-sections and rate constants are obtained for the title reaction on the latest ab initio potential energy surface. Reaction probabilities up to J=29 are dependent on the reactant rotation and show mild oscillations superimposed on a broad background. Due to a barrier in the entrance channel, the cross sections increase with energy with clear thresholds and the rate constants vary with temperature in the Arrhenius form. The calculated canonical rate constant is in good agreement with the experimental measurements. Our results also indicate that the quasiclassical trajectory method underestimates the rate due to the neglect of tunneling, while the quantum statistical approach overestimates because of the short lifetime of the reaction intermediate.     

An innovative computational comparison of exact and centrifugal sudden quantum properties of the N + N2 reaction

The development of an innovative computational strategy suited to provide an accurate quantum evaluation of the detailed properties of the N + N(2) exchange reaction has been undertaken by carrying out an extended theoretical study of such reaction. To this end exact and approximate quantum calculations (based on both time-independent and time-dependent techniques) of state-specific and state-to-state probabilities of the title reaction have been performed by considering values of the total angular momentum quantum number up to 20, values of total energy up to 2.3 eV and by making a combined use of both high throughput and high performance computing platforms. The comparison of the results obtained from calculations performed by taking into account the full Coriolis coupling of the allowed helicity states with those obtained when neglecting the Coriolis coupling or even a model energy shift treatment has allowed us to find out when a workflow managing the distribution of the jobs can replace exact treatments with approximate ones and for what type of properties this is possible.     

The H + Li2 bimolecular exchange reaction: dynamical and kinetical properties at J = 0

For the first time in the literature, rigorous time-independent quantum scattering formalism was applied, by means of the ABC program, to the H + Li(2) → LiH + Li reaction. The state-to-state probabilities as a function of the total energy have been computed at zero total angular momentum (J = 0) allowing us to evaluate the effect of vibrational/rotational excitation on the reaction promotion/inhibition, the energetic distribution of products, and the temperature dependence of the J-shifting thermal rate coefficients.     

Ar+ ArH+ Reactive Collisions of Astrophysical Interest: The Case of 36Ar

The reactive collision between ³⁶Ar and the ³⁶ArH⁺ species has been investigated by means of quantum mechanical (QM), quasiclassical trajectories (QCT) and statistical quantum mechanical (SQM) approaches. Reaction probabilities, cross sections as a function of the energy and rate constants in terms of the temperature have been obtained. Cumulative distributions as a function of the collision time and the inspection of selected QCT corresponding to specific dynamical mechanisms have been analysed. Predictions by means of the SQM method are in good agreement with the QM results, thus supporting the complex-forming nature of the process.     

Rate coefficient calculations on the three-dimensional D + HCl → DCl + H reaction using a new quantum-dynamical approximation

Rate coefficient calculations for the three-dimensional D + HCl→DC1 + H reaction are reported. A new quantum method is used, which involves applying the energy sudden approximation to the entrance channel and the centrifugal sudden approximation to the exit channel. Good agreement with quasiclassical results is obtained.