Sandip
Ghosh
^{a},
Rahul
Sharma
^{b},
Satrajit
Adhikari
*^{a} and
António J. C.
Varandas
*^{c}
^{a}Department of Physical Chemistry, Indian Association for the Cultivation of Science, 2A & 2B Raja S. C. Mullick Road, Jadavpur, Kolkata-700032, West Bengal, India. E-mail: pcsa@iacs.res.in
^{b}Department of Chemistry, St. Xaviers’ College, Kolkata-700016, West Bengal, India
^{c}Departamento de Química, and Centro de Química, Universidade de Coimbra, 3004-535 Coimbra, Portugal. E-mail: varandas@uc.pt

Received
13th September 2017
, Accepted 28th November 2017

First published on 28th November 2017

We perform quantum dynamics calculations for the reaction H + O_{2} → O + OH on the ground-state potential energy surfaces CHIPR [Varandas, J. Chem. Phys., 2013, 138, 134117] and DMBE IV [Pastrana et al., J. Phys. Chem. 1990, 94, 8073] using a three-dimensional time-dependent wave packet formalism based on hyperspherical coordinates. Initial rovibrational state [O_{2}(v = 0–4, j = 1–5)] dependent reaction probabilities are calculated for the case J = 0. The J-shifting scheme is employed to estimate initial state selected integral cross-sections as well as thermal rate coefficients, which is verified using a realistic extrapolation scheme. The calculated total and state-to-state rate coefficients are compared with the findings of recent experimental studies and previous theoretical calculations.

The development of an accurate quantum dynamics algorithm^{16–18} led to the evaluation of accurate integrals and differential cross-sections for atom–diatom collision processes. Various TD quantum mechanical approaches then appeared to assess reaction probabilities for a variety of reactive systems.^{19–21} On the other hand, the Coupled Channel (CC) approach in hyperspherical coordinates^{22–26} removed the awkward problem of coordinate transformation from the reactants to the products channel and made the formulation convenient for an equivalent description of all the rearrangement channels. Moreover, the extension of a reactive scattering calculation in hyperspherical coordinates on multiple electronic states for non-zero total angular momentum is straightforward to implement numerically rather than the traditional Jacobi coordinate system. The computational cost of such a CC approach scales as (JΠN_{i})^{3}, with the use of the grid method reducing the cost to , where J is the total angular momentum and N_{i} is the number of basis functions for the ith internal coordinate. Billing and Markovic^{25} formulated such a four-dimensional (4D) quantum mechanical problem of triatomic reactive scattering using hyperspherical coordinates and performed calculations for the D + H_{2} system in the simplest cases of J = 0, and J = 1. This approach was further transmuted by Adhikari and Varandas^{22} for the extension to higher J ≠ 0 situations. The state-to-state as well as total cross sections and rate coefficients have been accurately calculated by Adhikari and coworkers via the implementation of an OpenMP-MPI parallelized algorithm for the TD wave packet approach to the D^{+} + H_{2} system using both the lowest adiabatic singlet sheet (1^{1}A′)^{23,24} and lowest three singlet sheets (1^{1}A′, 2^{1}A′ and 3^{1}A′)^{27} of the double many-body expansion (DMBE) PES for H_{3}^{+} reported by Varandas and coworkers.^{28} The results were found to be in good agreement with experimental observations^{29–31} and previous theoretical calculations.^{32,33} Very recently, Adhikari and coworkers^{34} constructed accurate diabatic potential energy surfaces for the lowest three singlet states of H_{3}^{+}^{35} using first principle-based beyond Born–Oppenheimer theory^{36} and performed quantum dynamics using this hyperspherical coordinate-based scattering code on those surfaces to calculate state-to-state integral cross sections and found an overwhelming agreement with most recent available experimental data for all three competing charge transfer channels. In a similar context, Adhikari, Varandas and coworkers employed the same methodology for the O + OH reaction^{37} using both the (combined-hyperbolic-inverse-power-representation) CHIPR^{38} and DMBE IV PESs^{39} of HO_{2} and the J-shifting approximation to calculate state-to-state thermal rate coefficients. Moreover, 3D as well as coupled 3D time-dependent wave packet approaches have been employed for several other triatomic reactions in Jacobi coordinates, viz., F + H_{2},^{40,41} O(^{1}D) + N_{2},^{42} H + D_{2},^{43} O + H_{2},^{44} D + H_{2} and H + D_{2}^{45}etc. and in hyperspherical coordinates, viz., D + H_{2},^{22,25} D^{+} + H_{2},^{23,24,27,34,46} F + H_{2},^{47}etc. On the other hand, going beyond the atom–diatom to polyatomic reactive scattering is a big challenge in quantum dynamics, since the major bottleneck is to handle the exponential increase in computational cost due to the increase in mathematical dimensionalities when the number of atoms involved in the collision increases from three to four. For example, the dimensionality in terms of internal degrees of freedom increases from three for a triatomic system to six for a tetra-atomic system. On the contrary, the main computational advantage of a time-dependent approach is its smaller computational scaling with the number of basis functions (≤N^{2}vs. N^{3} in the standard CC time-independent method). Owing to its favorable computational scaling, the time-dependent approach is preferred for some large-scale computational tasks for various large molecules in reaction dynamics. In addition to that, utilization of hyperspherical coordinates is advantageous as it is able to handle all the reactant and product channels evenhandedly. To summarize, the rigorous dynamical treatment for tetra-atomic reactions is still not performed in hyperspherical coordinates, although a few calculations have been reported by using the time-dependent wave packet methodology in Jacobi coordinates for various tetra-atomic reactions, viz. H + H_{2}O,^{48–51} H + H_{2}S,^{52} H + CO_{2},^{53} OH + D_{2}^{54}etc.

The reaction between H(^{2}S) and O_{2}(^{3}Σ^{−}_{g}) is one of the most important elementary reactions in combustion chemistry. It is considered to be the rate-determining step in the combustion of hydrogen and hydrocarbons, and consumes ∼80% of O_{2} producing hydrocarbon-air flames in the atmosphere.^{55} Moreover, the title reaction is considered to be a prototype for complex-forming chemical reactions.^{56} Owing to such important attributes, many electronic structure calculations have been performed to construct a potential energy surface (PES) for the title system, and thereafter reactive scattering calculations are carried out to obtain temperature-dependent rate coefficients. Consequently, a kinetics study of accurate rate coefficients for the title reaction is highly demanding. Recall that the reverse reaction is also very important in interstellar chemistry, where it has a primary role in the destruction of ozone.

Numerous experimental measurements of the rate coefficient for the title reaction have been accomplished. Cohen and Westberg^{57} in 1983 reported a rate constant for the temperature range 250–2000 K. Later on, Baulch et al.^{58} obtained rate coefficients for various combustion processes including the title reaction, and Warnatz^{59} accumulated data for the same reaction. In turn, Shin and Michael^{60} carried out experiments by employing a laser photolysis-shock tube technique to check the rate constant in the temperature range 1103 ≤ T/K ≤ 2055, which were the first direct rate constant measurements in such a high temperature range. Owing to the endothermic nature of the reaction, the temperature-dependent rate constant is found to be Arrhenius in nature. To have a complete understanding of the reaction mechanism, several theoretical calculations have been performed to construct the potential energy surface with different levels of accuracy. Several versions of the electronically adiabatic PES have been reported, starting with Melius and Blint^{61} (MB). Subsequently, Pastrana et al.^{39} reported the DMBE IV surface by using dynamical simulations. In turn, Kendrick and Pack^{62} reported a diatomics-in-molecules type (DIMKP) PES. There are other available forms of ab initio surfaces by Troe and Ushakov (TU)^{63} and XXZLG PES of Xu et al.^{64,65} For the present calculations, we have employed the CHIPR^{38} PES, the most recent and likely accurate PES for the title system, which utilizes ab initio data from the XXZLG PES^{64} and the CHIPR method of Varandas.^{66} Dynamical calculations on the popular DMBE IV PES^{39} are also reported.

On the computational side, a large number of quantum dynamics studies have been carried out for the title reaction. Recently, Teixidor and Varandas employed the CHIPR^{38} PES of HO_{2} for dynamics calculations with the ABC time-independent quantum scattering program based on hyperspherical coordinates. They were able to evaluate thermal rate constants by the J-shifting technique,^{67} which were compared with the theoretical calculations^{68} on XXZLG and DIMKP PESs as well as experimental findings.^{57–59} The same authors reported similar calculations^{69} on the CHIPR PES for the X + O_{2} reactive system, with X being atomic hydrogen isotopes (Mu, H, D and T). Previously, Varandas^{70,71} utilized the trajectory calculations on the DMBE IV PES and calculated the thermal rate coefficients for the title reaction via an extrapolation method,^{70} whereas Pack et al.^{72,73} performed a TI quantum dynamics study on the same PES using the adiabatically adjusting principal-axis hyperspherical coordinates (APH) method. In turn, Duchovic et al. studied the H + O_{2} reaction on the MB and DMBE IV PESs by using quasiclassical trajectory (QCT) calculations.^{74} A while later, Lin et al.^{75} calculated the state-to-state reaction probabilities of the same reaction on the XXZLG PES with time-independent as well as time-dependent quantum mechanical methods. Moreover, the title reaction was extensively investigated by Bergueno et al.^{76} using time-dependent wave packet (TDWP) and statistical methods to calculate reaction probabilities and integral cross sections. Additionally, Guo and coworkers^{77} employed a fully Coriolis coupled time-dependent quantum mechanical methodology to calculate reaction probabilities for J ≠ 0 situations, whereas Sun et al.^{78} predicted the nonstatistical nature of the reaction. On the other hand, the effect of initial rotational state excitation on the reactivity of the title reaction has been studied by Guo, Lendvay and coworkers^{79} by utilizing the statistical phase space theory, the QCT method and the quantum wave packet approach. In turn, Viel et al.^{80} utilized a Lanczos formalism based on Green′s function and the DMBE IV PES to obtain J = 0 reaction probabilities and then, the J-shifting approach to calculate thermal rate coefficients. Meijer and Goldfield,^{81–83} in a series of studies using TD quantum mechanics with body fixed (BF) Jacobi coordinates on the DMBE IV PES, reported total reaction probabilities and cross sections for J = 0 and J > 0 situations, which highlighted the importance of the close-coupling approach. In turn, Sultanov and Balakrishnan^{84} performed TI quantum calculations on both DMBE IV and TU PESs for computing J = 0 cumulative reaction probabilities and thereafter J-shifted rate coefficients. A while later, Lin and coworkers,^{85} performed both the TD and TI quantum dynamics calculations on the XXZLG and DMBE IV PESs, whereas Hankel et al.^{86} calculated reaction probabilities for both the J = 0 and J > 0 situations by using a TD quantum method and the DMBE IV, TU and XXZLG PESs. In turn, Honvault et al.^{87} presented both the differential and integral cross sections for non-zero total angular momenta up to 50. Moreover, Quéméner et al.^{68} performed quantum dynamics calculations on the XXZLG and DIMKP PESs to study the initial rovibrational excitation of the reactant O_{2} molecule on the state-to-state reaction probabilities as well as J-shifted rate coefficients, and compared the calculated integral cross sections with available experimental measurements.^{88–90}

To summarize, we have investigated the following reaction for the formation of OH(^{2}Π):

H(^{2}S) + O_{2}(^{3}Σ^{−}_{g}) → OH(^{2}Π) + O(^{3}P), |

(1) |

(2) |

cosΘ′ = −sinβcosγ, | (3) |

sinξ = cosβ/sinΘ′, | (4) |

Moreover, A_{Kμ} obeys the orthonormality condition,

(5) |

(6) |

ϕ
_{
vj
}(r) is a diatomic Morse wave function, χ(R) is the translational wave function corresponding to the A + BC relative motion, and J(Rrη|ρθϕ) is the Jacobian for the transformation from Jacobi to hyperspherical coordinates.^{22,23}

The Hamiltonian operator describing the three-particle system can be written in terms of Johnson's hypersherical coordinates:^{91}

(7) |

Upon substituting the wave packet and the Hamiltonian into the TD Schrödinger equation (TDSE), the following set of coupled equations can be obtained:

(8) |

(9) |

(10) |

(11) |

(12) |

(13) |

(14) |

(15) |

When the wave packet passes through the projection region and gets almost absorbed at the boundary, the scattering amplitudes are Fourier transformed from time to energy space:

(16) |

(17) |

h^{±}(k_{v′j′}R*) = −n_{l}(k_{v′j′}R*) ± ij_{l}(k_{v′j′}R*), | (18) |

(19) |

(20) |

(21) |

(22) |

(23) |

The initial state selected cross sections by J-shifting can be obtained as:^{68}

(24) |

E^{‡}_{JK} = B^{‡}J(J + 1) + (A^{‡} − B^{‡})K^{2}, | (25) |

The thermal rate coefficient can be calculated by simply applying the J-shifting approximation as:

(26) |

(27) |

(28) |

(29) |

(30) |

On the other hand, considering the general form of the integral cross section,

(31) |

(32) |

(33) |

(34) |

V(ρ,θ,ϕ) = V_{0}(ρ,θ,ϕ) − iV_{Im}(ρ). | (35) |

where ρ

Panel-(a) of Fig. 1 depicts the total reaction probabilities calculated as a function of collision energy for the reaction H + O_{2}(v = 0; j = 1) → OH(v′,j′) + O at J = 0 performed on the adiabatic CHIPR PES, whereas panel-(b) displays the same, but on the DMBE IV PES. In panel-(b), the previous time-dependent quantum mechanical calculations on the DMBE IV PES for the reaction probability at J = 0 are shown over the collision energy range 0.0–1.8 eV by Hankel and coworkers^{86} and Abu Bajeh et al.^{89} While comparing the results from the CHIPR PES with DMBE IV, all the reaction profiles are seen to be in close agreement except some apparent differences due to resonances. Specifically, the CHIPR reaction probability profile shows a threshold value close to ∼0.67 eV, and gradually acquires higher values from there onwards, a pattern also consistent with the previous calculations performed on the DMBE IV PES. The existence of such a threshold value is due to the appearance of an effective barrier in the exit channel which corresponds to the opening of the OH(v′,j′) + O products channel. Note also the fact that the magnitude of the profile for the DMBE IV case is in good accordance with the previous theoretical calculations^{86,89} except for some peaks and troughs of the profiles up to ∼1.8 eV; after a collision energy of ∼2.0 eV, the present profile tends to slightly decrease. On the other hand, Fig. 2 depicts the calculated cumulative reaction probability (CRP) on the CHIPR PES for the H + O_{2} reaction as a function of total energy due to all possible initial rovibrational states of the reactant, O_{2}. Also included for comparison in this plot are the results extracted from the time-independent quantum scattering calculations of Quéméner et al.^{68} on the DIMKP PES (magenta curve) and the results of Teixidor and Varandas^{93} on the CHIPR PES (blue curve). It may be noted that we have calculated the CRP by involving reaction probabilities for all the initial rovibrational states (v, j) of an O_{2} diatom as well as all the final rovibrational states (v′, j′) of the OH diatom considered in our calculation [see Table 1], where the state-to-state reaction probability profiles at various collision energies are well converged with respect to final rotational (j′) as well as final vibrational (v′) quantum numbers. Additionally, the convergence profiles for state-to-state reaction probabilities at four different collision energies, viz., E_{c} = 0.50 eV, 1.0 eV, 1.5 eV and 2.0 eV are depicted in the ESI,† which clearly indicate that the reaction probability profiles are gradually converging as j′ increases from 0 to 15 and v′ increases from 0 to 10. It is quite noteworthy that the trends of the other theoretical profiles obtained from the DIMKP and CHIPR PESs are quite similar to the presently calculated one on the CHIPR PES. Indeed, despite some mismatch in the magnitudes of different profiles, the overall features of all the profiles remain the same.

Fig. 1 Total reaction probability for the reaction H + O_{2} (v = 0, j = 1) → O + OH (v′,j′) compared with previous theoretical investigations.^{86,89} Panel-(a): on the CHIPR PES. Panel-(b): on the DMBE IV PES. |

Fig. 2 Cumulative reaction probability (N^{J=0}(E)) as a function of the total energy for the reaction H + O_{2} → O + OH in comparison with previous theoretical investigations on the CHIPR^{93} and DIMKP^{68} PESs. |

Fig. 3 shows the initial rovibrational state selected total reaction probabilities calculated on the CHIPR PES as a function of collision energy. It is quite evident from the profiles of the reaction probability for various initial rotational states of the O_{2} molecule that the vibrational excitation enhances the reaction probability to a large extent, whereas rotational excitation does not have that much influence on the total reaction probability. With the increase of initial vibrational state from v = 0 to v = 1, 2 and 3, the profiles for total reaction probability are seen to shift towards the left, which indicates a decrease in the threshold collision energy due to smaller effective barriers for higher vibrational states. Notably, for the initial state v = 3, such a reaction threshold diminishes to ∼0.1 eV. On the other hand, with the increase in initial rotational excitation from j = 1 to j = 3 and 5, the profiles shift slightly towards lower collision energies which causes only small alterations on the overall magnitude of the reaction probability. As it is visible, the trend of the profiles so obtained is more-or-less consistent with the calculations of Quéméner et al.^{68}

Fig. 3 Initial rotational state resolved total reaction probability for the H + O_{2} (v = 0, 1, 2, 3; j = 1, 3, 5) reaction as a function of collision energy. |

Once the reaction probability is calculated for zero total angular momentum, J = 0, it is easy to calculate the integral cross section by invoking the J-shifting scheme. The integral cross sections for the initial v = 0, j = 1 rovibrational state of the O_{2} molecule in CHIPR and DMBE IV PES over the collision energy range 0.0–2.5 eV are displayed in Fig. 4. Other theoretical calculations by Varandas^{70} on the DMBE IV PES,^{39} and Quéméner et al. on the DIMKP PES as well as XXZLG PES^{68} are also shown together along with the available experimental findings by Keβler et al.,^{88} Abu Bajeh et al.^{89} and Seeger et al.^{90} The current profiles of integral cross section (both CHIPR and DMBE IV) show a threshold value of collision energy of ∼0.6 eV. Our calculated profile on the CHIPR PES shows a broad maximum at around ∼1.7 eV, which is in accordance with the general trend of experimental measurements by Keβler et al.^{88} and Seeger et al.,^{90} whereas the maximum seems to appear around ∼1.6 eV for the DMBE IV case. Notably, the present profile calculated on the CHIPR PES generally attains higher magnitude compared to the previous calculation by Quéméner et al., who performed TD quantum mechanical scattering calculations at J = 0 based on hyperspherical coordinates and invoked the J-shifting approximation to evaluate the integral cross-section. On the other hand, the calculated integral cross section on the DMBE IV PES is in quite good agreement with the previous QCT calculation by Varandas.^{70} In summary, the present result on the CHIPR PES shows good agreement with the experimental measurements of Keβler et al. and Seeger et al. especially in the lower collision energy range, whereas the profile for DMBE IV matches well with the same experimental profiles^{88,90} in the higher energy regime. Interestingly, our profiles also agree in shape with the curve reported by Varandas,^{70} who proposed to use quasiclassical trajectories to “extrapolate” accurate quantum mechanical results for J = 0, thence an alternative to the J-shifting scheme. Some quantitative differences between his results and the ones reported here may naturally be attributed to both the methodology and distinct potential energy surfaces that have been utilized.

Fig. 4 Integral cross section for the H + O_{2} (v = 0, j = 1) reaction as a function of collision energy in comparison with previous theoretical results and experimental measurements. |

While preforming the dynamics on the CHIPR PES for higher rovibrational states, viz. v = 1, 2, 3 and j = 1, 3, 5 of the O_{2} molecule, the state-specific integral cross sections shown in Fig. 5 seem to be largely influenced by the vibrational excitation from v = 0 to v = 1, 2 and 3, whereas the rotational excitation has little contribution to the magnitude of the cross section. Generally, the pattern of all the profiles is retained for all those rovibrational states, but the threshold collision energy gradually decreases with the higher vibrational state. Thus, we can say that with the increase in vibrational excitation, the title reaction tends gradually to be barrierless in nature. On the other hand, it is clearly visible that the initial rotational excitation of O_{2} does not alter the magnitude of the integral cross section, except for v = 0 where the profiles for j = 3 and 5 show a larger magnitude compared to the j = 1 case in the lower collision energy regime (0.0–1.5 eV). In addition, the threshold value decreases minutely with the increase in initial rotational state for all those initial vibrational states. Once again, the trends of our calculated profiles are in accordance with the previous calculations of Quéméner et al.^{68} Expectedly, by applying J-shifting approximation, the narrower and sharper resonance structures in the reaction probability at J = 0 get transmuted to ICS profiles. Those reaction probabilities for various initial rovibrational states at J = 0 show pronounced resonance structures (see Fig. 3) originating from the long-lived resonance states of HO_{2} due to its deep potential well.

Fig. 5 Initial rovibrational state resolved integral cross section for the H + O_{2} (v = 0, 1, 2, 3; j = 1, 3, 5) reaction as a function of the collision energy. |

While exploring the reliability of Bowman's J-shifting scheme,^{67} we have used Varandas's^{70} extrapolation scheme to calculate integral cross sections. Fig. 6 and 7 depict the results obtained from both methods: (a) J-shifting^{67}versus (b) extrapolation of quantum mechanical J = 0 reaction probabilities with the help of QCT calculations for specific J values.^{70}Fig. 6 shows that our reported integral cross sections on the CHIPR PES by J-shifting are in excellent agreement with the integral cross sections obtained via the extrapolation method. On the other hand, Fig. 7 also depicts that the cross sections calculated by these two approaches for the DMBE IV PES are in close agreement with each other, except for the fact that the J-shifted cross sections slightly overestimate the extrapolated results for E_{coll} beyond 1.5 eV. Of course, one is led to think that the extrapolation scheme is probably more reliable for high collision energies, where the correspondence principle suggests that quantum results should approach their classical counterparts. It is perhaps not surprising that unlike the J-shifting case, the extrapolated cross section curve exhibits a profound oscillatory pattern, which basically originates from the rich resonance structure of the P^{J=0} (v = 0,j = 1) profile. Thence, it should be viewed on an averaged sense. Since extrapolation of quantum results via classical mechanics is a quite natural thing to do for molecules of any size and it agrees quite well with the J-shifting results, the present theoretical prediction of experimental cross sections using both such methods appears to corroborate our better outcome when compared with other theoretical calculations.

Fig. 7 Same as Fig. 6, but on the DMBE IV PES. |

Employing the J = 0 reaction probability profile and using the J-shifting scheme, one may estimate the temperature-dependent rate coefficient. Fig. 8 depicts the Arrhenius plot of the Boltzmann averaged rate coefficient over the initial rovibrational states (v = 0–4; j = 1, 3, 5) as a function of inverse temperature for the temperature range 200–2000 K. Also included for comparison in Fig. 8 are previous results calculated by Teixidor and Varandas^{93} on the CHIPR PES and experimental measurements of Baulch et al.,^{58} Cohen and Westberg et al.^{57} and Warnatz.^{59} As shown, our results agree well with the experimental findings of Baulch et al.,^{58} particularly at the high temperature regime. Regarding the calculated profile of Teixidor and Varandas,^{93} no notable differences are observed except for the fact that ours lies slightly lower in magnitude compared to theirs over an intermediate temperature range of ∼500–1000 K.

To illustrate the effect of rovibrational excitation on the reactivity, Fig. 9 displays the thermal rate coefficients for various initial rovibrational states of O_{2} (v = 0–4; j = 1, 3, 5). From the Arrhenius plot of such state-specific thermal rate coefficients, one observes that rotational excitation does not show a large influence on the rate constant. Conversely, as we go to higher vibrational states, the profile of the rate coefficient changes drastically in the low temperature regime, whereas the effect is less pronounced at higher temperatures. In turn, the less negative slope of the rate coefficient profile for higher initial vibrational states indicates that the title reaction is evolving towards being barrierless in nature. In this case too, the trends of our calculated profiles are more-or-less consistent with the previous TI calculations by Teixidor and Varandas^{93} and Quéméner et al.^{68} Moreover, it is quite evident that as we increase the vibrational state of the O_{2} molecule from 0 to 4 at a particular rotational state, the difference among the profiles diminishes with increasing vibrational excitation. Additionally, the difference between the profiles seems to decrease further for higher rotational states. Since the rate coefficient does not change much for vibrational states above v = 4, the contribution from those rovibrational states on reactivity may be considered as an accurate estimation of the thermal rate coefficient.

Fig. 9 Comparison of total rate constants for various initial rovibrational states of O_{2}, viz., v = 0–4; j = 1, 3, 5. |

The integral cross section and rate coefficient profiles for various initial rovibrational states have shown that the J-shifting approximation^{67} and the extrapolation technique^{70} are cost-effective approaches to obtain reliable estimation of such attributes. Quantum mechanical calculations using the fully close-coupling approach^{23} are computationally very demanding for all higher values of total angular momentum that are required to obtain converged results. Although the thermal rate coefficient reported here shows a slight deviation from the available experimental and theoretical results in the low temperature regime, the overall agreement is quite encouraging. On the other hand, the cross-section profile is also in good agreement with the experimental findings. We believe that the deviations can be overcome by performing fully close-coupled calculations for larger J values, an endeavor that we plan to tackle in future work by performing converged calculations for J > 0 cases with our parallelized quantum scattering code.

- A. Kuppermann, G. C. Schatz and M. Baer, J. Chem. Phys., 1976, 65, 4596–4623 CrossRef CAS .
- F. Webster and J. C. Light, J. Chem. Phys., 1989, 90, 265–299 CrossRef CAS .
- G. C. Schatz, Chem. Phys. Lett., 1988, 150, 92–98 CrossRef CAS .
- J. Z. H. Zhang and W. H. Miller, J. Chem. Phys., 1989, 91, 1528–1547 CrossRef CAS .
- M. Mladenovic, M. Zhao, D. G. Truhlar, D. W. Schwenke, Y. Sun and D. J. Kouri, J. Phys. Chem., 1988, 92, 7035–7038 CrossRef CAS .
- Y. T. Lee, Science, 1987, 236, 793 CAS .
- D. A. V. Kliner, D. E. Adelman and R. N. Zare, J. Chem. Phys., 1991, 95, 1648 CrossRef CAS .
- P. Casvecchia, Rep. Prog. Phys., 2000, 63, 355–414 CrossRef .
- X. Liu, J. J. Lin, S. A. Harich, G. C. Schatz and X. Yang, Science, 2000, 289, 1536–1538 CrossRef CAS PubMed .
- B. K. Kendrick, J. Chem. Phys., 2000, 112, 5679–5704 CrossRef CAS .
- D. Sokolovski and J. F. Castillo, Phys. Chem. Chem. Phys., 2000, 2, 507–512 RSC .
- A. J. H. M. Meijer and E. M. Goldfield, J. Chem. Phys., 1999, 110, 870–880 CrossRef CAS .
- D. H. Zhang, S. Y. Lee and M. Baer, J. Chem. Phys., 2000, 112, 9802–9809 CrossRef CAS .
- W. Zhu, T. Peng and J. Z. H. Zhang, J. Chem. Phys., 1997, 106, 1742–1748 CrossRef CAS .
- S. C. Althorpe, D. J. Kouri and D. K. Hoffmann, J. Chem. Phys., 1997, 106, 7629–7636 CrossRef CAS .
- J. C. Juanes-Marcos and S. C. Althorpe, Chem. Phys. Lett., 2003, 381, 743–750 CrossRef CAS .
- H. Gao, M. Sneha, F. Bouakline, S. C. Althorpe and R. N. Zare, J. Phys. Chem. A, 2015, 119, 12036–12042 CrossRef CAS PubMed .
- S. K. Gray and G. G. Balint-Kurti, J. Chem. Phys., 1998, 108, 950–962 CrossRef CAS .
- S. Y. Lin and H. Guo, Phys. Rev. A: At., Mol., Opt. Phys., 2006, 74, 022703/1–8 CAS .
- S. G. Carrasco and O. Roncero, J. Chem. Phys., 2006, 125, 054102/1–14 Search PubMed .
- Z. Sun, H. Guo and D. H. Zhang, J. Chem. Phys., 2010, 132, 084112/1–11 CAS .
- S. Adhikari and A. J. C. Varandas, Comput. Phys. Commun., 2013, 184, 270–283 CrossRef CAS .
- T. Sahoo, S. Ghosh, S. Adhikari, R. Sharma and A. J. C. Varandas, J. Phys. Chem. A, 2014, 118, 4837–4850 CrossRef CAS PubMed .
- T. Sahoo, S. Ghosh, S. Adhikari, R. Sharma and A. J. C. Varandas, J. Chem. Phys., 2015, 142, 024304/1–9 CrossRef CAS PubMed .
- G. D. Billing and N. Markovic, J. Chem. Phys., 1994, 100, 1085–1093 CrossRef .
- J. D. Kress, R. T. Pack and G. A. Parker, Chem. Phys. Lett., 1990, 170, 306–310 CrossRef CAS .
- S. Ghosh, T. Sahoo, S. Adhikari, R. Sharma and A. J. C. Varandas, J. Phys. Chem. A, 2015, 119, 12392–12403 CrossRef CAS PubMed .
- L. P. Viegas, A. Alijah and A. J. C. Varandas, J. Chem. Phys., 2007, 126, 074309/1–9 CrossRef CAS PubMed .
- F. C. Fehsenfeld, D. B. Dunkin, E. E. Ferguson and D. L. Albritton, Astrophys. J., 1973, 183, L25–L26 CrossRef CAS .
- D. Gerlich, Adv. Chem. Phys., 1992, 82, 1–176 CrossRef CAS .
- C. Schlier, U. Nowotny and E. Teloy, Chem. Phys., 1987, 111, 401–408 CrossRef CAS .
- T. S. Chu and K. L. Han, J. Phys. Chem. A, 2005, 109, 2050–2056 CrossRef CAS PubMed .
- P. Honvault and Y. Scribano, J. Phys. Chem. A, 2013, 117, 9778–9784 CrossRef CAS PubMed .
- S. Ghosh, S. Mukherjee, B. Mukherjee, S. Mandal, R. Sharma, P. Chaudhury and S. Adhikari, J. Chem. Phys., 2017, 147, 074105/1–16 CAS .
- S. Mukherjee, D. Mukhopadhyay and S. Adhikari, J. Chem. Phys., 2014, 141, 204306/1–13 CrossRef CAS PubMed .
- B. Sarkar and S. Adhikari, J. Chem. Phys., 2006, 124, 074101/1–18 CrossRef CAS PubMed .
- S. Ghosh, R. Sharma, S. Adhikari and A. J. C. Varandas, Chem. Phys. Lett., 2017, 675, 85–91 CrossRef CAS .
- A. J. C. Varandas, J. Chem. Phys., 2013, 138, 134117/1–7 CAS .
- M. R. Pastrana, L. A. M. Quintales, J. Brandao and A. J. C. Varandas, J. Phys. Chem., 1990, 94, 8073–8080 CrossRef CAS .
- D. Neuhauser, M. Baer, R. S. Judson and D. J. Kouri, Chem. Phys. Lett., 1990, 169, 372–379 CrossRef CAS .
- Y. Zhang, T.-X. Xie, K. L. Han and J. Z. H. Zhang, J. Chem. Phys., 2003, 119, 12921–12925 CrossRef CAS .
- T. S. Chu, T. X. Xie and K. L. Han, J. Chem. Phys., 2004, 121, 9352–9360 CrossRef CAS PubMed .
- R. F. Lu, T. S. Chu, Y. Zhang, K. L. Han, A. J. C. Varandas and J. Z. H. Zhang, J. Chem. Phys., 2006, 125, 133108/1–6 CAS .
- T. S. Chu, X. Zhang and K. L. Han, J. Chem. Phys., 2005, 122, 214301/1–6 CAS .
- T. S. Chu, A. J. C. Varandas and K. L. Han, Chem. Phys. Lett., 2009, 471, 222–228 CrossRef CAS .
- N. Markovic and G. D. Billing, Chem. Phys., 1995, 191, 247–260 CrossRef CAS .
- N. Markovic and G. D. Billing, Mol. Phys., 2000, 78, 1771–1781 CrossRef .
- D. H. Zhang and J. C. Light, J. Chem. Phys., 1996, 104, 4544–4553 CrossRef CAS .
- D. H. Zhang, M. Yang and S. Y. Lee, J. Chem. Phys., 2002, 117, 10067–10072 CrossRef CAS .
- D. H. Zhang and J. Z. H. Zhang, J. Chem. Phys., 1994, 101, 1146–1156 CrossRef CAS .
- N. Balakrishnan and G. D. Billing, J. Chem. Phys., 1994, 101, 2785–2792 CrossRef CAS .
- J. Qi, D. Lu, H. Song, J. Li and M. Yang, J. Chem. Phys., 2017, 146, 124303/1–8 CrossRef CAS PubMed .
- P. Sun, J. Chen, S. Liu and D. H. Zhang, Chem. Phys. Lett., 2017, 683, 352–356 CrossRef CAS .
- P. Y. Zhang, R. F. Lu, T. S. Chu and K. L. Han, J. Chem. Phys., 2010, 133, 174316/1–7 CAS .
- W. C. Gardiner, Combustion Chemistry, Springer, Berlin, 1984 Search PubMed .
- H. Guo, Int. Rev. Phys. Chem., 2012, 31, 1–68 CrossRef CAS .
- N. Cohen and K. R. Westberg, J. Phys. Chem. Ref. Data, 1983, 12, 531–590 CrossRef CAS .
- D. L. Baulch, C. J. Cobos, R. A. Cox, P. Frank, G. Hayman, T. Just, J. A. Kerr, T. Murrells, M. J. Pilling, J. Troe, R. W. Walker and J. Warnatz, J. Phys. Chem. Ref. Data, 1994, 23, 847–1033 CrossRef CAS .
- J. Warnatz, in Combustion Chemistry, ed. W. C. Gardiner, Springer-Verlag, New York, 1984 Search PubMed .
- K. S. Shin and J. V. Michael, J. Chem. Phys., 1991, 95, 262–273 CrossRef CAS .
- C. F. Melius and R. J. Blint, Chem. Phys. Lett., 1979, 64, 183–189 CrossRef CAS .
- B. Kendrick and R. T. Pack, J. Chem. Phys., 1995, 102, 1994–2012 CrossRef CAS .
- J. Troe and V. G. Ushakov, J. Chem. Phys., 2001, 115, 3621–3628 CrossRef CAS .
- C. Xu, D. Xie, D. H. Zhang, S. Y. Lin and H. Guo, J. Chem. Phys., 2005, 122, 244305 CrossRef PubMed .
- D. Xie, C. Xu, T. S. Ho, H. Rabitz, G. Lendvay, S. Y. Lin and H. Guo, J. Chem. Phys., 2007, 126, 074315 CrossRef PubMed .
- A. J. C. Varandas, J. Chem. Phys., 2013, 138, 054120 CrossRef CAS PubMed .
- J. Bowman, J. Phys. Chem., 1991, 95, 4960–4968 CrossRef CAS .
- G. Quéméner, B. K. Kendrick and N. Balakrishnan, J. Chem. Phys., 2010, 132, 014302 CrossRef PubMed .
- M. M. Teixidor and A. J. C. Varandas, Chem. Phys. Lett., 2015, 638, 61–65 CrossRef CAS .
- A. J. C. Varandas, Mol. Phys., 1995, 85, 1159–1164 CrossRef CAS .
- A. J. C. Varandas, Chem. Phys. Lett., 1995, 235, 111–118 CrossRef CAS .
- R. T. Pack, E. A. Butcher and G. A. Parker, J. Chem. Phys., 1993, 99, 9310–9313 CrossRef CAS .
- R. T. Pack, E. A. Butcher and G. A. Parker, J. Chem. Phys., 1995, 102, 5998–6012 CrossRef CAS .
- R. J. Duchovic and M. A. Parker, J. Phys. Chem. A, 2005, 109, 5883–5896 CrossRef CAS PubMed .
- S. Y. Lin, H. Guo, P. Honvault and D. Xie, J. Phys. Chem. B, 2006, 110, 23641–23643 CrossRef CAS PubMed .
- P. Bargueno, T. Gonzalez-Lezana, P. Larregaray, L. Bonnet and J. C. Rayez, Phys. Chem. Chem. Phys., 2007, 9, 1127–1137 RSC .
- S. Y. Lin, Z. Sun, H. Guo, D. H. Zhang, P. Honvault, D. Xie and S. Y. Lee, J. Phys. Chem. A, 2008, 112, 602–611 CrossRef CAS PubMed .
- Z. Sun, D. H. Zhang, C. Xu, S. Zhou, D. Xie, G. Lendvay, S. Y. Lee, S. Y. Lin and H. Guo, J. Am. Chem. Soc., 2008, 130, 14962–14963 CrossRef CAS PubMed .
- S. Y. Lin, H. Guo, G. Lendvay and D. Xie, Phys. Chem. Chem. Phys., 2009, 11, 4715–4721 RSC .
- A. Viel, C. Leforestier and W. H. Miller, J. Chem. Phys., 1998, 108, 3489–3497 CrossRef CAS .
- A. J. H. M. Meijer and E. M. Goldfield, J. Chem. Phys., 1998, 108, 5404–5413 CrossRef CAS .
- A. J. H. M. Meijer and E. M. Goldfield, J. Chem. Phys., 1999, 110, 870–880 CrossRef CAS .
- E. M. Goldfield and A. J. H. M. Meijer, J. Chem. Phys., 2000, 113, 11055–11062 CrossRef CAS .
- R. A. Sultanov and N. Balakrishnan, J. Phys. Chem. A, 2004, 108, 8759–8764 CrossRef CAS .
- S. Y. Lin, E. J. Rackham and H. Guo, J. Phys. Chem. A, 2006, 110, 1534–1540 CrossRef CAS PubMed .
- M. Hankel, S. C. Smith and A. J. H. M. Meijer, J. Chem. Phys., 2007, 127, 064316/1–9 CrossRef CAS PubMed .
- P. Honvault, S. Y. Lin, d. Xie and H. Guo, J. Phys. Chem. A, 2007, 111, 5349–5352 CrossRef CAS PubMed .
- K. Keβler and K. Kleinermanns, J. Chem. Phys., 1992, 97, 374–377 CrossRef .
- M. A. Bajeh, E. M. Goldfield, A. Hanf, C. Kappel, A. J. H. M. Meijer, H. R. Volpp and J. Wolfrum, J. Phys. Chem. A, 2001, 105, 3359–3364 CrossRef CAS .
- S. Seeger, V. Sick, H. R. Volpp and J. Wolfrum, Isr. J. Chem., 1994, 34, 5–18 CrossRef CAS .
- B. R. Johnson, J. Chem. Phys., 1983, 79, 1916–1925 CrossRef CAS .
- G. D. Billing and N. Markovic, J. Chem. Phys., 1993, 99, 2674–2681 CrossRef CAS .
- M. M. Teixidor and A. J. C. Varandas, J. Chem. Phys., 2015, 142, 014309/1–10 CrossRef CAS PubMed .
- D. Neuhasuer and M. Baer, J. Chem. Phys., 1989, 90, 4351–4355 CrossRef CAS .
- D. Kosloff and R. Kosloff, J. Comput. Phys., 1986, 52, 35–53 CrossRef .

## Footnote |

† Electronic supplementary information (ESI) available. See DOI: 10.1039/c7cp06254k |

This journal is © the Owner Societies 2018 |