Real wave packet and quasiclassical trajectory studies of the H+ + LiH reaction

Time-dependent real wave packet (RWP) and quasiclassical trajectory (QCT) calculations have been carried out to study the H(+) + LiH reaction on the ab initio potential-energy surface of Martinazzo et al. [J. Chem. Phys., 2003, 119, 11241]. Total initial state-selected and final state-resolved reaction probabilities for the two possible reaction channels, H(2)(+) + Li and LiH + H(+), have been calculated for total angular momentum J=0 at a broad range of collision energies. Integral cross sections and thermal rate coefficients have been calculated using the QCT method and from the corresponding J=0 RWP reaction probabilities by means of a capture model. The calculated thermal rate coefficients are found to be nearly independent of temperature in the 100-500 K interval with a value of approximately 10(-9) cm(3) s(-1), which is in good agreement with estimates used in evolutionary models of early-Universe lithium chemistry. The RWP results are found to be in good agreement overall with the corresponding QCT calculations.     

Quantum chemical and conventional TST calculations of rate constants for the OH + alkane reaction

Reactions of OH with methane, ethane, propane, i-butane, and n-butane have been modeled using ab initio (MP2) and hybrid DFT (BHandHLYP) methods, and the 6-311G(d,p) basis set. Furthermore, single-point calculations at the CCSD(T) level were carried out at the optimized geometries. The rate constants have been calculated using the conventional transition-state theory (CTST). Arrhenius equations are proposed in the temperature range of 250–650 K. Hindered Internal Rotation partition functions calculations were explicitly carried out and included in the total partition functions. These corrections showed to be relevant in the determination of the pre-exponential parameters, although not so important as in the NO₃ + alkane reactions [G. Bravo-Pérez, J.R. Alvarez-Idaboy, A. Cruz-Torres, M.E. Ruíz, J. Phys. Chem. A 106 (2002) 4645]. The explicit participation of the tunnel effect has been taken into account. The calculated rate coefficients provide a very good agreement with the experimental data. The best agreement for the overall alkane + OH reactions seemed to occur when the BHandHLYP geometries and partition functions are used. For propane and i-butane, in addition to the respective secondary and tertiary H-abstraction channels, the primary one has been considered. These pathways are confirmed to be significant in spite of the large differences in activation energies between primary and secondary or primary and tertiary channels, respectively of propane and i-butane reactions and should not be disregarded.     

Quantum wave packet calculation of reaction probabilities,cross sections,and rate constants for H+ + LiH → Li + H reaction

The H⁺ + LiH → Li + H reactive scattering has been studied using a quantum real wave packet method. The state‐to‐state and state‐to‐all reaction probabilities for the entitled collision have been calculated at zero total angular momentum. The probabilities for J > 0 are estimated from the J = 0 results by using J‐shifting approximation based on the Capture model. The integral cross sections and thermal rate constants are then calculated. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Quantum wave packet calculation of reaction probabilities,cross sections,and rate constants for the C(1D) + HD reaction

The time‐dependent real wave packet method has been used to study the C(¹D) + HD reaction. The state‐to‐state and state‐to‐all reactive scattering probabilities for a broad range of energies are calculated at zero total angular momentum. The probabilities for J > 0 are estimated from accurately computed J = 0 probabilities by using the J‐shifting approximation. The integral cross sections for a large energy range, and thermal rate constants are calculated. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Cross sections and low temperature rate coefficients for the H + CH+ reaction: a quasiclassical trajectory study

The H + CH(+) reaction is studied by quasiclassical trajectory (QCT) calculations, along with phase space theory (PST) and quantum rigid rotor calculations, employing a global single-valued potential energy surface recently derived by our group. We report QCT total cross sections for each of the three channels, for low collision energies and different reactant rotational quantum numbers. At the lowest collision energies, all cross sections exhibit a capture-like behaviour, as expected from a barrierless reaction. At higher energies, there are important dynamical effects coming from the opening of new channels in the inelastic and reactive exchange collisions. The inelastic cross sections turn out to largely increase, while the reactive abstraction cross sections are declining faster than predicted by the capture theory. A large value of the reactant rotational quantum number tends to suppress these dynamical effects. The QCT rate coefficients are reported for a temperature range from 1-700 K. Below 20 K, the abstraction and exchange QCT rate coefficients are almost constant, as predicted by the capture theory. Above this temperature, the abstraction rate coefficient declines, while the exchange and inelastic rate coefficients are increasing, due to the opening of new channels. A good agreement is observed between the experimental abstraction rate coefficient and the QCT and PST ones. The QCT inelastic results are also compared with those obtained from rigid rotor close coupling (CCRR) calculations in order to check the ability of this approach to provide a reliable estimate of the inelastic rate coefficients for a reactive system without a barrier. The laws of variation as a function of temperature are found to be very similar and the curves are parallel above 20 K. However, reaction is not allowed in the rigid rotor approximation, therefore the CCRR results are about twice as large as their QCT counterparts.     

Reflected shock tube studies of high-temperature rate constants for OH + C2H2 and OH + C2H4

The reflected shock tube technique with multi-pass absorption spectrometric detection of OH-radicals at 308 nm (corresponding to a total path length of approximately 4.9 m) has been used to study the reactions, OH + C(2)H(2)--> products (1) and OH + C(2)H(4)--> C(2)H(3) + H(2)O (2). The present optical configuration gives a S/N ratio of approximately 1 at approximately 0.5-1.0 x 10(12) radicals cm(-3). Hence, kinetics experiments could be performed at [OH](0) = approximately 4-20 ppm thereby minimizing secondary reactions. OH was produced rapidly from the dissociations of either CH(3)OH or NH(2)OH (hydroxylamine). A mechanism was then used to obtain profile fits that agreed with the experiment to within <+/-5%. The derived Arrhenius expressions, in units of cm(3) molecule(-1) s(-1) are: k(1) = (1.03 +/- 0.24) x 10(-10) exp(-7212 +/- 417 K/T) for 1509-2362 K and k(2) = (10.2 +/- 5.8) x 10(-10) exp(-7411 +/- 871 K/T) for 1463-1931 K. The present study is the first ever direct measurement for reaction (1) at temperatures >1275 K while the present results extend the temperature range for (2) by approximately 700 K. These values are compared with earlier determinations and with recent theoretical calculations. The calculations agree with the present data for both reactions to within +/-10% over the entire T-range.     

Reflected shock tube studies of high-temperature rate constants for CH3 + O2, H2CO + O2, and OH + O2

The reflected shock tube technique with multipass absorption spectrometric detection of OH-radicals at 308 nm, corresponding to a total path length of approximately 2.8 m, has been used to study the reaction CH3 + O2 CH2O + OH. Experiments were performed between 1303 and 2272 K, using ppm quantities of CH3I (methyl source) and 5-10% O2, diluted with Kr as the bath gas at test pressures less than 1 atm. We have also reanalyzed our earlier ARAS measurements for the atomic channel (CH3 + O2 --> CH3O + O) and have compared both these results with other earlier studies to derive a rate expression of the Arrhenius form. The derived expressions, in units of cm3 molecule(-1) s(-1), are k = 3.11 x 10(-13) exp(-4953 K/T) over the T-range 1237-2430 K, for the OH-channel, and k = 1.253 x 10(-11) exp(-14241 K/T) over the T-range 1250-2430 K, for the O-atom channel. Since CH2O is a major product in both reactions, reliable rates for the reaction CH2O + O2 --> HCO + HO2 could be derived from [OH]t and [O]t experiments over the T-range 1587-2109 K. The combined linear least-squares fit result, k = 1.34 x 10(-8) exp(-26883 K/T) cm3 molecule(-1) s(-1), and a recent VTST calculation clearly overlap within the uncertainties in both studies. Finally, a high sensitivity for the reaction OH + O2 --> HO2 + O was noted at high temperature in the O-atom data set simulations. The values for this obtained by fitting the O-atom data sets at later times (approximately 1.2 ms) again follow the Arrhenius form, k = 2.56 x 10(-10) exp(-24145 K/T) cm3 molecule(-1) s(-1), over the T-range, 1950-2100 K.     

Time-dependent wave packet studies on the Cl + HCl hydrogen exchange reaction

The initiation of the hydrogen exchange reaction Cl((2)P)+HCl --> ClH+Cl((2)P) by excitation of the HCl molecular stretch to v=2 is studied for total angular momentum quantum number J=(1)/(2) and both even and odd parity. The calculations were performed using a time-dependent propagation from an initial quasi-bound state and employed all three relevant potential energy surfaces and the nonadiabatic couplings between them. Coriolis and spin-orbit coupling were also taken into account. The electronic and HCl rotational distributions of the products in both dissociation channels are analyzed, and the results are interpreted using features of the potential energy surfaces.     

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.     

Communication: rate coefficients from quasiclassical trajectory calculations from the reverse reaction: The Mu + H2 reaction re-visited

In a previous paper [P. G. Jambrina et al., J. Chem. Phys. 135, 034310 (2011)] various calculations of the rate coefficient for the Mu + H(2) → MuH + H reaction were presented and compared to experiment. The widely used standard quasiclassical trajectory (QCT) method was shown to overestimate the rate coefficients by several orders of magnitude over the temperature range 200-1000 K. This was attributed to a major failure of that method to describe the correct threshold for the reaction owing to the large difference in zero-point energies (ZPE) of the reactant H(2) and product MuH (～0.32 eV). In this Communication we show that by performing standard QCT calculations for the reverse reaction and then applying detailed balance, the resulting rate coefficient is in very good agreement with the other computational results that respect the ZPE, (as well as with the experiment) but which are more demanding computationally.     

Methyl vinyl ketone+OH and methacrolein+OH oxidation reactions: a master equation analysis of the pressure- and temperature-dependent rate constants

High-level electronic structure calculations and master equation analyses were carried out to obtain the pressure- and temperature-dependent rate constants of the methyl vinyl ketone+OH and methacrolein+OH reactions. The balance between the OH addition reactions at the high-pressure limit, the OH addition reactions in the fall-off region, and the pressure-independent hydrogen abstractions involved in these multiwell and multichannel systems, has been shown to be crucial to understand the pressure and temperature dependence of each global reaction. In particular, the fall-off region of the OH addition reactions contributes to the inverse temperature dependence of the rate constants in the Arrhenius plots, leading to pressure-dependent negative activation energies. The pressure dependence of the methyl vinyl ketone+OH reaction is clearly more important than in the case of the methacrolein+OH reaction owing to the weight of the hydrogen abstraction process in this second system. Comparison of the theoretical rate constants and the experimental measurements shows quite good agreement.     

Cross sections and rate constants for the C(3P)+OH(X 2Pi)-->CO(X 1Sigma+)+H(2S) reaction using a quasiclassical trajectory method

First quasiclassical trajectory calculations have been carried out for the C(3P)+OH(X 2Pi)-->CO(X 1Sigma+)+H(2S) reaction using a recent ab initio potential energy surface for the ground electronic state, X 2A', of HCO/COH. Total and state-specific integral cross sections have been determined for a wide range of collision energies (0.001-1 eV). Then, thermal and state-specific rate constants have been calculated in the 1-500 K temperature range. The thermal rate constant varies from 1.78x10(-10) cm3 s-1 at 1 K down to 5.96x10(-11) cm3 s-1 at 500 K with a maximum value of 3.39x10(-10) cm3 s-1 obtained at 7 K. Cross sections and rate constants are found to be almost independent of the rovibrational state of OH.     

Quantum and quasiclassical studies of the O(3P)+HCl-->OH+Cl(2P) reaction using benchmark potential surfaces

We have performed quantum mechanical (QM) dynamics calculations within the independent-state approximation with new benchmark triplet A" and A' surfaces [B. Ramachandran et al., J. Chem. Phys. 119, 9590 (2003)] for the rovibronic state-to-state measurements of the reaction O(3P)+HCl(v=2,j=1,6,9)-->OH(v'j')+Cl(2P) [Zhang et al., J. Chem. Phys. 94, 2704 (1991)]. The QM and experimental rotational distributions peak at similar OH(j') levels, but the QM distributions are significantly narrower than the measurements and previous quasiclassical dynamics studies. The OH(low j) populations observed in the measurements are nearly absent in the QM results. We have also performed quasiclassical trajectory with histogram binning (QCT-HB) calculations on these same benchmark surfaces. The QCT-HB rotational distributions, which are qualitatively consistent with measurements and classical dynamics studies using other surfaces, are much broader than the QM results. Application of a Gaussian binning correction (QCT-GB) dramatically narrows and shifts the QCT-HB rotational distributions to be in very good agreement with the QM results. The large QCT-GB correction stems from the special shape of the joint distribution of the classical rotational/vibrational action of OH products. We have also performed QM and QCT calculations for the transition, O+HCl(v=0,T=300 K)-->OH(v'j')+Cl from threshold to approximately 130 kcal mol(-1) collision energy as a guide for possible future hyperthermal O-atom measurements. We find in general a mixed energy release into translation and rotation consistent with a late barrier to reaction. Angular distributions at high collision energy are forward peaked, consistent with a stripping mechanism. Direct collisional excitation channel cross sections, O+HCl(v=0,T=300 K)-->O+HCl(v'=1), in the same energy range are large, comparable in magnitude to the reactive channel cross sections. Although the (3)A" state dominates most collision processes, above approximately 48 kcal mol(-1), the (3)A' state plays the major role in collisional excitation.     

Quasiclassical trajectory simulations of OH(v) + NO2 --> HONO2* --> OH(v') + NO2: capture and vibrational deactivation rate constants

Quasiclassical trajectory calculations are used to investigate the dynamics of the OH(v) + NO(2) --> HONO(2) --> OH(v') + NO(2) recombination/dissociation reaction on an analytic potential energy surface (PES) that gives good agreement with the known structure and vibrational frequencies of nitric acid. The calculated recombination rate constants depend only weakly on temperature and on the initial vibrational energy level of OH(v). The magnitude of the recombination rate constant is sensitive to the potential function describing the newly formed bond and to the switching functions in the PES that attenuate inter-mode interactions at long range. The lifetime of the nascent excited HONO(2) depends strongly not only on its internal energy but also on the identity of the initial state, in disagreement with statistical theory. This disagreement is probably due to the effects of slow intramolecular vibrational energy redistribution (IVR) from the initially excited OH stretching mode. The vibrational energy distribution of product OH(v') radicals is different from statistical distributions, a result consistent with the effects of slow IVR. Nonetheless, the trajectory results predict that vibrational deactivation of OH(v) via the HONO(2) transient complex is approximately 90% efficient, almost independent of initial OH(v) vibrational level, in qualitative agreement with recent experiments. Tests are also carried out using the HONO(2) PES, but assuming the weaker O-O bond strength found in HOONO (peroxynitrous acid). In this case, the predicted vibrational deactivation efficiencies are significantly lower and depend strongly on the initial vibrational state of OH(v), in disagreement with experiments. This disagreement suggests that the actual HOONO PES may contain more inter-mode coupling than found in the present model PES, which is based on HONO(2). For nitric acid, the measured vibrational deactivation rate constant is a useful proxy for the recombination rate, but IVR randomization of energy is not complete, suggesting that the efficacy of the proxy method must be evaluated on a case-by-case basis.     

Rate constants for the reaction OH + CO as functions of temperature and water concentration

Rate constants for the reaction OH + CO have been measured as functions of temperature (340–1220 K) and water concentration in the presence of 1 atm of argon. Results at zero water concentration yield the expression, log k_? (cm³ molecule^?1 S^?1) = ?12.96 + 4.7 × 10^?4 T, for the reaction rate constant as a function of temperature. These results are in very good agreement with previous direct measurements and in reasonable agreement with flame and shock tube measurements. Explanations are offered for the involvement of the water molecule in the present experiments and earlier measurements from this laboratory throughout the entire temperature range. Results are consistent with previous results showing little, if any, pressure effect of Ar on the reaction up to 1 atm of Ar.     

State-to-state dynamics analysis of the F + CHD3 reaction: a quasiclassical trajectory study

An exhaustive state-to-state dynamics study was performed to analyze the F + CHD3 --> FD(nu', j') + CHD2(nu) gas-phase abstraction reaction. Quasiclassical trajectory (QCT) calculations, including corrections to avoid zero-point energy leakage along the trajectories, were performed at different collision energies on an analytical potential energy surface (PES-2006) recently developed by our group. Whereas the CHD2 coproduct appears vibrationally and rotationally cold, most of the available energy appears as FD(nu') product vibrational energy, peaking at nu' = 2 and nu' = 3, with the population in the latter level growing as the energy increases. The excitation function rises from the threshold of the reaction and then levels off at higher energies, with the maximum contribution from the FD(nu' = 3) level. The state-specific FD(nu') scattering distributions correlated with the coproduct CHD2 in the nu4 = 2 and nu3 = 1 states, at different collision energies, show a steady change from backward to forward scattering as the energy increases. This similar behavior for the two coproduct vibrational states, nu4 = 2 and nu3 = 1, agrees qualitatively with the experimental measurements. Comparison with theoretical and experimental results for the isotopic analogues, F + CH4 and F + CD4, shows that the title reaction presents a direct mechanism, similar to the perdeuterated reaction, but contrasts with that of the F + CH4 reaction. These results for the dynamics of different isotopic variants, always in qualitative and sometimes in quantitative agreement with experiment, show the capacity of the PES-2006 surface to correctly describe the title reaction, even though there are differences that could be due to deficiencies of the PES but also to the known limitations of the classical treatment in the QCT method.     

A quasiclassical trajectory study of the reaction H + O2 <==> OH + O with the O2 reagent vibrationally excited

Quasiclassical trajectories have been computed on the Melius-Blint (MB) Potential Energy Surface (PES) and on the Double Many-Body Expansion (DMBE) IV PES of Pastrana et al. describing the H + O(2) <==> OH + O reaction with the nonrotating (J = 0) O(2) reagent vibrationally excited to levels v = 6, 7, 8, 9, and 10 at four temperatures: 1000, 2000, 3000, and 4000 K. The vibrational energy levels were selected by using a semiclassical Einstein-Brillouin-Keller (EBK) quantization procedure while the relative translational energy was sampled from a Boltzmann weighted distribution. The rate coefficient for the formation of the OH + O products is seen to increase monotonically with quantum number and nearly monotonically with temperature. On the MB PES, at T = 1000 K, the total rate coefficient increases by a factor of 5.2 as the initial vibrational quantum number of the O(2) diatom increases from v = 6 to v = 10. For T = 2000 K, this factor drops to 3.3, to 2.9 for T = 3000 K, and to 2.5 for T = 4000 K. On the DMBE IV PES, at T = 1000 K the total rate coefficient increases by a factor of 4.1 as the initial vibrational quantum number of the O(2) diatom increases from v = 6 to v = 10. For T = 2000 K, this factor drops to 3.5, to 2.1 for T = 3000 K, and to 2.0 for T = 4000 K. The less-direct group (defined below) of trajectories is sensitive to the initial O(2) vibrational excitation in several different temperature ranges, apparently retaining the effect of reagent vibrational excitation. The more-direct group (defined below) of trajectories does not exhibit this behavior. Reagent vibrational excitation does not increase the total rate coefficients for the title reaction more than the increase due to a simple temperature increase. The less-direct and more-direct groups of trajectories differ in their contribution to the rate coefficient for the title reaction. In particular, at T = 4000 K, the two PESs used in this work differ dramatically in the roles of the less-direct and more-direct trajectories. The behavior of the more-direct and less-direct groups of trajectories can be understood in terms of the efficiency of intramolecular vibrational energy transfer. This work utilizes the recently introduced PES Library, POTLIB 2001, which made the comparisons between the two PESs discussed in this work possible in a very straightforward way.     

Quantum mechanical rate constants for H + O2 <--> O + OH and H + O2 --> HO2 reactions

Canonical rate constants for both the forward and reverse H + O(2) <--> O + OH reactions were calculated using a quantum wave packet-based statistical model on the DMBE IV potential energy surface of Varandas and co-workers. For these bimolecular reactions, the results show reasonably good agreement with available experimental and theoretical data up to 1500 K. In addition, the capture rate for the H + O(2) --> HO(2) addition reaction at the high-pressure limit was obtained on the same potential using a time-independent quantum capture method. Excellent agreement with experimental and quasi-classical trajectory results was obtained except for at very low temperatures, where a reaction threshold was found and attributed to the centrifugal barrier of the orbital motion.     

Quantum wave packet study of the insertion reaction S+H2

S(¹D)+H₂→SH+H reaction at zero total angular momentum has been studied by using a time-dependent quantum real wave packet method. State-to-state and state-to-all reactive scattering probabilities for a broad range of energy are calculated. The probabilities show many sharp peaks that ascribed to reactive scattering resonances. The density plots of the wave function shows that the reaction presents a pure insertion pathway.     

Reflected shock tube studies of high-temperature rate constants for OH + NO2 --> HO2 + NO and OH + HO2 --> H2O + O2

The motivation for the present study comes from the preceding paper where it is suggested that accepted rate constants for OH + NO2 --> NO + HO2 are high by approximately 2. This conclusion was based on a reevaluation of heats of formation for HO2, OH, NO, and NO2 using the Active Thermochemical Table (ATcT) approach. The present experiments were performed in C2H5I/NO2 mixtures, using the reflected shock tube technique and OH-radical electronic absorption detection (at 308 nm) and using a multipass optical system. Time-dependent profile decays were fitted with a 23-step mechanism, but only OH + NO2, OH + HO2, both HO2 and NO2 dissociations, and the atom molecule reactions, O + NO2 and O + C2H4, contributed to the decay profile. Since all of the reactions except the first two are known with good accuracy, the profiles were fitted by varying only OH + NO2 and OH + HO2. The new ATcT approach was used to evaluate equilibrium constants so that back reactions were accurately taken into account. The combined rate constant from the present work and earlier work by Glaenzer and Troe (GT) is k(OH+NO2) = 2.25 x 10(-11) exp(-3831 K/T) cm3 molecule(-1) s(-1), which is a factor of 2 lower than the extrapolated direct value from Howard but agrees well with NO + HO2 --> OH + NO2 transformed with the updated equilibrium constants. Also, the rate constant for OH + HO2 suitable for combustion modeling applications over the T range (1200-1700 K) is (5 +/- 3) x 10(-11) cm3 molecule(-1) s(-1). Finally, simulating previous experimental results of GT using our updated mechanism, we suggest a constant rate for k(HO2+NO2) = (2.2 +/- 0.7) x 10(-11) cm3 molecule(-1) s(-1) over the T range 1350-1760 K.