Quasiclassical Trajectory Calculations of the Rate Constant of the OH

Feb 3, 2014 - quasiclassical trajectory calculations with a focus on the thermal rate constant, k(T), over the temperature range 5 to 500 K. Compariso...
1 downloads 0 Views 660KB Size
Letter pubs.acs.org/JPCL

Quasiclassical Trajectory Calculations of the Rate Constant of the OH + HBr → Br + H2O Reaction Using a Full-Dimensional Ab Initio Potential Energy Surface Over the Temperature Range 5 to 500 K Antonio G. S. de Oliveira-Filho,†,‡ Fernando R. Ornellas,† and Joel M. Bowman*,‡ †

Departamento de Química Fundamental, Instituto de Química, Universidade de São Paulo, São Paulo, São Paulo, 05508-000, Brazil Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, United States



S Supporting Information *

ABSTRACT: We report a permutationally invariant, ab initio potential energy surface (PES) for the OH + HBr → Br + H2O reaction. The PES is a fit to roughly 26 000 spin-free UCCSD(T)/cc-pVDZ-F12a energies and has no classical barrier to reaction. It is used in quasiclassical trajectory calculations with a focus on the thermal rate constant, k(T), over the temperature range 5 to 500 K. Comparisons with available experimental data over the temperature range 23 to 416 K are made using three approaches to treat the OH rotational and associated electronic partition function. All display an inverse temperature dependence of k(T) below roughly 160 K and a nearly constant temperature dependence above 160 K, in agreement with experiment. The calculated rate constant with no treatment of spin− orbit coupling is overall in the best agreement with experiment, being (probably fortuitously) within 20% of it. SECTION: Kinetics and Dynamics he OH + HBr → H2O + Br reaction is the major process by which atomic bromine, an active ozone depletion catalyst, is regenerated from the relatively stable hydrogen bromide.1,2 Thermal rate constants for this reaction have been measured over a broad temperature range (23−416 K), and revealed a strong inverse temperature dependence of the rate constant below 150 K.3−14 This surprising behavior of rapid reactivity at very low temperatures was observed, using the CRESU (Cinétique de Réaction en Encoulement Supersonique Uniforme) method.9,10 This method also has been used to study many reactions between neutral species,15,16 making this class of reactions important for the understanding of the chemistry of low-temperature molecular clouds in the interstellar medium. From a fundamental perspective, these reactions pose challenges for the understanding of the underlying dynamics responsible for this behavior. Several approaches have been used in previous theoretical work on the OH + HBr reaction,17−22 but none of these studies conclusively address the issue of the barrier height, and the role of long-range interactions and tunneling in the observed reaction rates. Clary et al.17 used a rotationally adiabatic capture theory with a long-range electrostatic interaction potential to predict a maximum for the rate constant at the temperature of 20 K, a feature that was not observed experimentally, perhaps due to the absence of measurements in the very lowtemperature limit (T < 10 K). Later, Clary et al.18 developed a London−Eyring−Polanyi−Sato (LEPS) potential that was used in a reduced dimensionality quantum scattering

T

© 2014 American Chemical Society

calculations. Based on the dependence of the reaction cross section on the OH rotational quantum number, they obtained a rate constant with a T−1/2 dependence, in accord with the earlier capture model. However, agreement with available experimental data was not good, with theory below experiment and experiment not showing the maximum found in the calculations. This potential, which has no potential barrier, was also used in trajectory calculations, conducted by Nizamov et al.,19 with a focus on experiments that determined the H2O vibrational energy distribution. In the most recent calculations, Liu et al.21 reported a direct ab initio approach to this reaction, employing reaction-path dual-level direct dynamics variational transition-state theory, and attributed the observed lowtemperature inverse T-dependence of the rate constant to tunneling through a small barrier. As noted in the Abstract, the present results find no barrier, and the disagreement with these most recent ones are focused on below. In this Letter, we report a full-dimensional potential energy surface (PES) for the title reaction. This is a permutationally invariant fit to roughly 26 000 high-level ab initio energies. We also report an accurate determination of the energetics of the stationary points of products of this reaction. Quasiclassical trajectories (QCT) were used to calculate the rate constant Received: January 6, 2014 Accepted: February 3, 2014 Published: February 3, 2014 706

dx.doi.org/10.1021/jz5000325 | J. Phys. Chem. Lett. 2014, 5, 706−712

The Journal of Physical Chemistry Letters

Letter

over the temperature range of 5−500 K. The present study also relates to recent investigations of the X + H2O → HX + OH (X = F and Cl) forward and reverse reactions,23−31 for which ab initio potential energy surfaces have been reported. Note in these cases there are reaction barriers for the diatom−diatom reaction, in contrast to the present case. This series of systems contributes greatly to the general understanding of chemical reactivity due to the larger range of reaction dynamics, compared to atom−diatom reaction, that can be studied. The PES is a fit to electronic energies calculated using spinunrestricted, explicitly correlated coupled-cluster 32−34 (UCCSD(T)-F12a) theory with restricted open-shell Hartree−Fock (ROHF) reference functions. These calculations employed the cc-pVDZ-F1235−37 orbital basis sets, together with the following choice of auxiliary basis sets: the cc-pVDZ/ OptRI38 sets were used for the resolution of the identity, the ccpVTZ/JKFIT (for H and O) and QZVPP/JKFIT39,40 (for Br) sets were used for the density fitting of the exchange and Fock operators, while the aug-cc-pVTZ/MP2FIT39 (for H and O) and cc-pVTZ-F12/MP2FIT (for Br) sets were used for the remaining two-electron integrals. In all calculations, except when otherwise indicated, the inner core of bromine was replaced by a relativistic pseudopotential.41 The full-dimensional potential energy surface was constructed using the invariant polynomial method via monomial symmetrization.42,43 The expression for the potential, in this approach, is given by V (y1 , ..., y6 ) =



Figure 1. Schematic reaction path and stationary points for the OH + HBr → H2O + Br reaction. Energies in kcal mol−1.

Cn1...n6y1n1 y6n6 [y2n2 y3n3 y4n4 y5n5 + y2n4 y3n5 y4n2 y5n3 ]

n1,..., n6

(1)

where yi = e−ri/α with α = 2a0, and ri is an internuclear distance, namely: r1 = rHH′, r2 = rHO, r3 = rHBr, r4 = rH′O, r5 = rH′Br and r6 = rOBr. The total order of polynomials involved in the expression (N = n1 + n2 + n3 + n4 + n5 + n6) does not exceed six, resulting in a total of 502 terms. The expansion coefficients in eq 1 were determined by a standard least-squares fit of 26 121 electronic energies. These configurations were first sampled using grids around the stationary points on the PES that were used to fit a tentative potential, on which batches of exploratory trajectories were calculated, at several collision energies. To ensure an accurate description of the long-range behavior, the data set extended to R values as long as 8.5 Å, where R is the distance between the centers of mass of the reactants. Configurations in regions with unphysical behavior were added to the data set, that was then fitted, providing a new tentative potential. This procedure was repeated in order to ensure that all relevant features for the description of the OH + HBr → H2O + Br were properly represented. The overall root-mean-square error (RMSE), that includes points as high as 250 kcal mol−1 above the reactants, is 0.84 kcal mol−1. Restricting the energy threshold to points up to 50 kcal mol−1 above the reactants reduces RMSE to 0.52 kcal mol−1. A schematic reaction path and the energetics of the stationary points are shown in Figure 1. As seen, there is a van der Waals well (one of two), followed by a negative energy saddle-point barrier, both in the entrance channel, followed by a van der Waals well in the product channel. A contour plot in Figure 2 displays a relaxed equipotential contour plot in which the energy, in kcal mol−1 and relative to the OH(re) + HBr(r′e ) asymptotic limit, is shown as a function of the rOH and rHBr distances, for optimized values of the other geometrical parameters. As seen, the potential is smooth and the key

Figure 2. Equipotential contour plot of the H2OBr PES in (rOH′, rH′Br) space with all other internal coordinates optimized. The zero of energy is the OH(re) + HBr(re) asymptotic limit, and the equipotential lines are separated by 5 kcal mol−1.

structures (van der Waals minima and the saddle point) on the reaction path are easily identified. From the perspective of the dynamics, this reaction profile, with a small negative-energy barrier, and large exoergicity suggests that tunneling plays a minor or no role in the dynamics. Thus, quasiclassical trajectory calculations should be reliable for both the rate constant and final-state properties of the products. Here the emphasis is on the rate constant, which as noted above, displays an interesting inverse temperature dependence at low temperatures and is nearly constant, at least up to roughly 500 K. Complementing these calculations, optimized geometries and harmonic frequencies were determined for the stationary points at the CCSD(T)-F12a/cc-pVnZ-F12 (n = D and T) levels of theory, and also using conventional coupled cluster with single, double, and perturbative estimate of connected triples, 44,45 CCSD(T), employing the aug-cc-pVnZ(PP)41,46−48 (n = T and Q) basis sets. Finally, in order to obtain accurate benchmark energetics, a series of single-point calculations at the optimized CCSD(T)F12a/aVTZ-F12 geometries was carried out for the stationary points of the PES, using a procedure similar to the one used to 707

dx.doi.org/10.1021/jz5000325 | J. Phys. Chem. Lett. 2014, 5, 706−712

The Journal of Physical Chemistry Letters

Letter

Table 1. Energies (in kcal mol−1) and Geometries (in Ångström and Degrees) of the Stationary Points on the H2OBr PES

calculate an accurate barrier height of the O + HBr reaction.49−51 The complete basis set (CBS) limit was estimated using a 1/l3max extrapolation52 of CCSD(T) energies calculated with the aug-cc-pVQZ(-PP) and the aug-cc-pV5Z(PP) basis sets.41,46−48 Additive corrections to include the effect of full iterative triple and quadruple excitations, core−valence correlation, scalar relativity and spin−orbit have been calculated. The triple excitation correction (ΔT) was defined as the difference between CCSD(T) and CCSDT53−55 calculations using the cc-pVTZ(-PP) basis set; similarly, the quadruple excitation correction (ΔQ) was calculated from CCSDT and CCSDTQ56−58 energies using the cc-pVDZ(-PP) basis set. The core−valence correlation was included at the CCSD(T) level using the aug-cc-pwCVQZ(-PP) basis set59,60 as the difference between frozen-core calculations and those with the O 1s and Br 3s3p3d electrons correlated (both in the same basis set). The second-order Douglas−Kroll−Hess (DKH) Hamiltonian61 along with the DKH-contracted basis sets,62,63 aug-cc-pVQZ-DK, was used to include scalar relativistic effects on H and O, as well as to also correct for the pseudopotential approximation. This contribution was defined as the difference between CCSD(T)/aug-cc-pVQZ(PP) and DK-CCSD(T)/aug-cc-pVQZ-DK energies within the frozen-core approximation. The inclusion of spin−orbit effects was calculated using the full Breit−Pauli operator within the state interacting approach64 using dynamically weighted state averaged complete active space self-consistent field65−67 (DWSA-CASSCF) wave functions using the aug-cc-pVTZ(-PP) basis sets in a calculation for five states, with a β−1 factor of 3 eV, employing valence active space. This additive correction was defined as the difference between nonspin−orbit and spin− orbit energies. All CCSD(T), CCSD(T)-F12a, and CASSCF calculations were carried out with the Molpro68,69 suite of programs, while CCSDT and CCSDTQ calculations were performed with the MRCC70 as interfaced to Molpro. Table 1 presents CCSD(T) energies and geometries of the relevant indicated stationary points using a variety of basis sets, along with results from the PES. As seen the results from the various basis sets all agree well and the PES reproduces the VDZ-F12 results, upon which it is based, very well. In several instances the PES is fortuitously closer to the “TZ” basis results than the VDZ-F12 ones. Focusing in particular on the saddle point in the entrance channel, we see that all results give a negative energy for this. This is one example where the PES energy of −0.52 kcal/mol is quite close to the “TZ” basis results. A comparison of corresponding harmonic frequencies is given in Supporting Information (SI), where the PES is shown to be generally in good agreement with the direct ab initio results. The frequencies shown there are substantially different from the ones reported previously for LEPS potential18 and also the lower-level direct ab initio calculations.21 Note that harmonic frequencies are not employed in a direct way in the present quasiclassical trajectory calculations. An examination of larger basis sets and extrapolated complete basis set results, spin−orbit interaction, scalar relativistic, core−valence effects, and full iterative triple and quadruple excitations is given in SI. Based on these, the best estimate for the saddle point energy is −0.61 kcal mol−1, which is 0.09 kcal mol−1 below the PES value. As expected spin−orbit coupling is largest in the region of the products, Br···H2O and the products themselves. The best estimate for the product energy is −35.43 kcal mol−1, which, as expected, is lower than

θHOH′

θOH′Br

ϕHOH′Br

−2.96 −2.91 −2.83

1.415 R-vdW (HO···HBr) 0.975 2.148 1.426 112.6 0.972 2.127 1.426 112.9 0.972 2.132 1.424 113.3

173.1 173.6 173.7

0.0 0.0 0.0

−2.85

0.972

112.8

173.5

0.0

−3.19

0.972

173.1

0.0

−2.13 −1.98 −1.67

2.316 1.421 116.5 R′-vdW (HBr···HO) 0.976 3.804 1.421 29.2 0.973 3.778 1.420 31.0 0.973 4.128 1.418 18.6

71.7 73.2 62.0

0.0 0.0 0.0

−1.67

0.973

20.3

65.9

0.0

−1.99

0.973

basis

E

rHO

aVTZ aVQZ VDZF12 VTZF12 PES

0.0 0.0 0.0

0.973 0.971 0.971

OH + HBr 1.420 1.419 1.417

0.0

0.971

1.420

0.0

0.971

aVTZ aVQZ VDZF12 VTZF12 PES aVTZ aVQZ VDZF12 VTZF12 PES aVTZ aVQZ VDZF12 VTZF12 PES aVTZ aVQZ VDZF12 VTZF12 PES aVTZ aVQZ VDZF12 VTZF12 PES

rOH′

2.129

4.036

rH′Br

1.426

1.421

20.8

66.3

0.0

−0.33 −0.45 −0.34

3.954 1.418 Saddle (HOHBr‡) 0.975 1.504 1.479 0.972 1.507 1.478 0.972 1.507 1.477

101.9 102.1 101.9

141.3 140.7 139.6

55.4 54.3 55.4

−0.55

0.972

102.0

141.7

53.3

−0.52

0.973

1.516

1.478

105.7

143.7

73.7

−35.61 −36.01 −35.28

1.551 1.487 P-vdW (Br···H2O) 0.963 0.963 3.120 0.961 0.961 3.090 0.961 0.961 3.078

104.4 104.7 104.7

61.6 61.5 61.4

−103.9 −103.9 −103.9

−36.01

0.961

0.961

3.086

104.7

61.4

−103.9

−35.39

0.960

104.8

61.9

−103.5

−32.33 −32.63 −31.86

0.962 0.959 0.959

0.960 3.070 Br + H2O 0.962 0.959 0.959

−32.61

0.959

0.959

104.4

−32.36

0.965

0.965

107.3

104.2 104.4 104.3

the values of the spin-free PES and aVQZ values, −32.36 and −32.63 kcal mol−1, respectively. This ca. 10% difference is not expected to play a significant role in the dynamics calculations of the reaction cross section. Finally, we note that we examined the widely used T1 diagnostic and find that it is below 0.02 along most of the reaction path, but it does reach a maximum of 0.07 at the vicinity of the negative-energy saddle point. This latter value is somewhat larger than the “safe” 0.02 threshold; however, the UCCSD(T) calculations in that region of the PES do not produce obviously “spurious” results and so we feel reasonably confident that the final fitted PES is smooth, realistic and reliable. 708

dx.doi.org/10.1021/jz5000325 | J. Phys. Chem. Lett. 2014, 5, 706−712

The Journal of Physical Chemistry Letters

Letter

gelOH(T ) = 2(1 + e−|A| hc / kBT )

The QCT calculations were performed using standard procedures and expressions, 71,72 employing the code VENUS.73,74 The maximum impact parameter, (bmax), is strongly dependent on the collision energy for this reaction, ranging from 15.75 Å, in the limit of low collision energies, to 3.4 Å for the most energetic collisions. The vibration−rotation energy of the reactants was determined from the initial quantum states, using semiclassical quantization of the anharmonic diatom potentials from the PES. For a given set of initial parameters (collision energy, initial states of the reactants), the reactive cross section, using Monte Carlo sampling of the impact parameter, is given by 2 σr = πbmax

This product is indicated as “RR/SO” and if spin−orbit coupling is neglected, i.e., gOH el (T) is replaced by the factor 4, the result is indicated “RR/nSO”. A more accurate procedure to obtain the partition function involves a detailed quantum mechanical treatment of the energy levels in a fully coupled calculation.75,76 In this case, the partition function is given by, Ω= 3/2states OH Q coupled (T ) = 2[

+ e−|A| hc / kBT

(2)

(3)

where P(jOH; T) and P(jHBr; T) are the Boltzmann distributions for jOH and jHBr at temperature T; and then integrating over the thermal distribution of collision energies. The final expression for the thermal rate constant is given by,



σr(Ec ; T )Ec 2

(kBT )

e−Ec / kBT dEc (4)

where kB is the Boltzmann constant, μ is the reduced mass for the collision (mOHmHBr/(mOH + mHBr)), ge is the electronic degeneracy of the reactive surface, and Qreact(T) is internal partition function of the reactants. The integral in eq 4 was calculated using a cubic spline interpolation of the averaged cross sections. The part of the reactant partition function describing the internal energies of the OH radical has received special attention in the literature, as its ground state is 2Π. Thus, there is both spin−orbit coupling and lambda-doubling to consider. Ignoring the latter, as is typically done in rate constant calculations, two versions of the rotational partition function have been used previously in rate constant calculations. The most common expression is given by a product of the spin-free rotational partition function times an electronic factor to account for the spin−orbit splitting of the Ω (the projection of the total electronic angular momentum on the molecular axis) = 1/2 and 3/2 states. Using rigid-rotor energy levels the factored partition function is given by Qrr(T) gOH el (T), where Q rr =

∑ (2J + 1)e−J(J+ 1)Bhc/k T B

J

(2J + 1)e−(EJ − EJ=1/2)/ kBT ] (7)

where EJ is the energy of a level with total angular momentum J, obtained from a numerical solution of the Schrödinger equation for the asymptotic fitted potential, calculated for states with Ω equal to 1/2 and 3/2, assuming that these states are separated by the experimental spin−orbit splitting constant,77 A. The details of these calculations and Hamiltonian are given in the SI. Results calculated using eq 7 for the partition function are indicated as “Coupled”. It could be argued that using the Coupled approach leads to an unbalanced treatment of rotational energies of OH in the partition function and as they are calculated by VENUS for the trajectories. In VENUS, as usual, they are obtained using a spinfree Hamiltonian, with the anharmonic diatom potential from the PES, and energies are obtained using standard and accurate EBK quantization. Dynamical consistency could be achieved by dynamically coupling the electronic and nuclear degrees of freedom either quantum mechanically, or with various schemes to couple nuclear and electronically degrees of freedom classically.78−80 However, that was not done in the present case, and so it might be argued even more strongly that since the present calculations are done using a single spin-free PES that any explicit consideration of spin−orbit coupling is inconsistent. In any case, the three expressions for the OH partition function will be used to calculate k(T). The thermal rate constants over the temperature range 5− 500 K obtained from QCT calculations and three different approaches for the partition function of the reactants are shown in Figure 3 along with available experimental data. The rate constant using reaction-path dual-level direct-dynamics variational transition-state theory21 is also shown. As seen, the present results with any treatment of the OH rotational partition function agree with experiment in showing a nearly constant k(T) at temperatures above roughly 180 K and a strong inverse temperature dependence below roughly 150 K. As also seen, the choice of OH-partition function does begin to play a major role below 150 K, where factors of between 2 and 4 are seen. As expected, the largest differences are at the lowest temperatures. Between 50 and 150 K, the RR/SO, RR/nSO, and Coupled results overestimate experiment on average by factors of 2.2, 1.2, and 1.5, respectively. Below 50 K, the RR/SO rate constants overestimate the measurements by a factor of 1.5, and the RR/nSO and Coupled calculations underestimate experiment by factors of 1.3 and 1.6, respectively. Overall, the best agreement with experiment is for the RR/nSO results, where the differences are roughly 20%. This level of agreement is probably fortuitously good.

σr(Ec , jOH , jHBr )P(jOH ; T )P(jHBr ; T )

ge 8kBT μπ Q react(T )

∑ J = 1/2

jOH , jHBr

k(T ) =

(2J + 1)e−(EJ − EJ=3/2)/ kBT

J = 3/2

where Nr is the number of reactive trajectories and Ntotal is the total number of trajectories. The initial state-specific cross sections were calculated for 30 values of collision energies (Ec) ranging from 0.001 to 25 kcal mol−1, with the semiclassical rotational quantum numbers ranging from 0 to 7, for both reactants in the vibrational ground state. The very low collision energy was needed to converge the rate constant to temperaturea as low as 5 K. The rate constant can be calculated by averaging the reactive cross section with respect to the population of the initial states





Ω= 1/2states

Nr Ntotal

σr(Ec ; T ) =

(6)

(5)

and 709

dx.doi.org/10.1021/jz5000325 | J. Phys. Chem. Lett. 2014, 5, 706−712

The Journal of Physical Chemistry Letters

Letter

in the QCT rate constant, however, is not evident in the experimental results, measured down to 23 K, and clearly it would be important for experiment to go below that temperature. It may be that the self-reaction OH + OH81,82 may affect to some extent the measurement of the OH + HBr rate constants at very low temperatures. To summarize, a new ab initio potential energy surface has been reported and used in quasiclassical calculations of the rate constant of the OH + HBr → Br + H2O reaction, which shows an unusual temperature dependence. The theoretical results are (finally) in good agreement with experiment. Further investigation of the product state distribution, isotope effects on the rate constants and in the energy disposal, and transitionstate theory calculation are being conducted and will be presented in a forthcoming paper. Finally, with the availability of the global PES, it would be interesting, albeit challenging, to conduct quantum calculations of the rate constant to investigate the importance of tunneling, e.g., through rotationally adiabatic barriers.



ASSOCIATED CONTENT

S Supporting Information *

Benchmark calculations of the energetics, harmonic vibrational frequencies of the stationary points, calculated and experimental thermal rate constants, and details of the determination of the energies used in the coupled partition function calculation. This material is available free of charge via the Internet at http://pubs.acs.org/.



Figure 3. Comparison of the thermal rate constants calculated for the OH + HBr reaction with the available experimental data, indicated as symbols (refs 3−14), and with the reaction-path dual-level direct dynamics variational transition-state theory calculations from Liu et al.21

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS F.R.O. acknowledges the Conselho Nacional de Desenvolví e Tecnológico (CNPq) of Brazil for academic mento Cientifico support. A.G.S.O.F. is thankful to Fundaçaõ de Amparo à Pesquisa do Estado de São Paulo (FAPESP). J.M.B. thanks the Department of Energy (DE-FG02-97ER14782) for financial support.

Note also that the three approaches give a maximum that is located at roughly 15 K for the RR/SO and RR/nSO calculations, and 50 K for the Coupled calculation. The existence of the maximum is in agreement with a prediction made by Clary et al.17 for reactions of molecules in 2Π electronic states with 1Σ polar molecules at low temperatures. Clary et al. did consider the OH + HBr reaction in the range 0 to 80 K, using adiabatic capture theory.17 The reported rate constant is in the range 30−35 × 10−11 cm3 molecule−1 s−1, which is significantly larger than the present calculated results and experiment. This is perhaps due to the assumption in the capture model of a unit reaction probability for collisions with impact parameters less than the predicted maximum one. At low collision energies, the QCT results find a smaller reaction probability owing to the formation of quasi-trapped trajectories due to the prereactive van der Waals well. A significant fraction of such trajectories emerge to reform the reactants. (Details of these trajectories will be given elsewhere.) Next consider the previous reaction-path dual-level directdynamics variational transition state theory calculations of Liu et al.21 As seen, they significantly underestimate experiment at temperatures above 60 K. They do increase very rapidly below that temperature, and no evidence of a temperature maximum is seen. Overall, the present results represent a major improvement of the theoretical modeling of the OH + HBr reaction. The level of agreement with experiment is very good. The maximum



REFERENCES

(1) Yung, Y. L.; Pinto, J. P.; Watson, R. T.; Sander, S. P. Atmospheric Bromine and Ozone Perturbations In the Lower Stratosphere. J. Atmos. Sci. 1980, 37, 339−353. (2) Wayne, R. P. Chemistry of Atmospheres; Oxford University Press: Oxford, 2000. (3) Takacs, G. A.; Glass, G. P. Reactions of Hydrogen Atoms and Hydroxyl Radicals With Hydrogen Bromide. J. Phys. Chem. 1973, 77, 1060−1064. (4) Smith, I. W. M.; Zellner, R. Rate Measurements of Reactions of OH By Resonance-Absorption 3. Reactions of OH With H2, D2 and Hydrogen and Deuterium Halides. J. Chem. Soc., Faraday Trans. 2 1974, 8, 1045−1056. (5) Husain, D.; Plane, J. M. C.; Slater, N. K. H. Kinetic Investigation of the Reactions of OH(X2Π) With the Hydrogen Halides, HCl, DC1, HBr and DBr By Time-resolved Resonance Fluorescence (A2Σ+− X2Π). J. Chem. Soc., Faraday Trans. 2 1981, 77, 1949−1962. (6) Jourdain, J. L.; Lebras, G.; Combourieu, J. EPR Determination of Absolute Rate Constants For the Reactions of H and OH Radicals with Hydrogen Bromide. Chem. Phys. Lett. 1981, 78, 483−487. 710

dx.doi.org/10.1021/jz5000325 | J. Phys. Chem. Lett. 2014, 5, 706−712

The Journal of Physical Chemistry Letters

Letter

(27) Li, J.; Jiang, B.; Guo, H. Reactant Vibrational Excitations Are More Effective than Translational Energy in Promoting an EarlyBarrier Reaction F + H2O → HF + OH. J. Am. Chem. Soc. 2013, 135, 982−985. (28) Li, J.; Dawes, R.; Guo, H. Kinetic and Dynamic Studies of the Cl(2Pu) + H2O (X1A1) → HCl(X1Σ+) + OH(X2Π) Reaction on an Ab Initio Based Full-Dimensional Global Potential Energy Surface of the Ground Electronic State of ClH2O. J. Chem. Phys. 2013, 139, 074302. (29) Guo, Y.; Zhang, M.; Xie, Y.; Schaefer, H. F., III. Communication: Some Critical Features of the Potential Energy Surface for the Cl + H2O → HCl + OH Forward and Reverse Reactions. J. Chem. Phys. 2013, 139, 041101. (30) Li, J.; Li, Y.; Guo, H. Communication: Covalent Nature of X− H2O (X = F, Cl, and Br) Interactions. J. Chem. Phys. 2013, 138, 141102. (31) Jiang, B.; Guo, H. Control of Mode/Bond Selectivity and Product Energy Disposal by the Transition State: X + H2O (X = H, F, O(3P), and Cl) Reactions. J. Am. Chem. Soc. 2013, 135, 15251−15256. (32) Adler, T. B.; Knizia, G.; Werner, H.-J. A Simple and Efficient CCSD(T)-F12 Approximation. J. Chem. Phys. 2007, 127, 221106. (33) Knizia, G.; Adler, T. B.; Werner, H.-J. Simplified CCSD(T)-F12 Methods: Theory and Benchmarks. J. Chem. Phys. 2009, 130, 054104. (34) Werner, H.-J.; Knizia, G.; Manby, F. R. Explicitly Correlated Coupled Cluster Methods with Pair-Specific Geminals. Mol. Phys. 2011, 109, 407−417. (35) Peterson, K. A.; Adler, T. B.; Werner, H.-J. Systematically Convergent Basis Sets for Explicitly Correlated Wavefunctions: The Atoms H, He, B−Ne, and Al−Ar. J. Chem. Phys. 2008, 128, 084102. (36) Peterson, K. A.; Krause, C.; Stoll, H.; Hill, J. G.; Werner, H.-J. Application of Explicitly Correlated Coupled-Cluster Methods to Molecules Containing Post-3d Main Group Elements. Mol. Phys. 2011, 109, 2607−2623. (37) Peterson, K. A.; Hill, J. G. in preparation (sets available upon request). (38) Yousaf, K. E.; Peterson, K. A. Optimized Auxiliary Basis Sets for Explicitly Correlated Methods. J. Chem. Phys. 2008, 129, 184108. (39) Weigend, F.; Köhn, A.; Hättig, C. Efficient Use of the Correlation Consistent Basis Sets in Resolution of the Identity MP2 Calculations. J. Chem. Phys. 2002, 116, 3175. (40) Weigend, F. A Fully Direct RI-HF Algorithm: Implementation, Optimised Auxiliary Basis Sets, Demonstration of Accuracy and Efficiency. Phys. Chem. Chem. Phys. 2002, 4, 4285−4291. (41) Peterson, K. A.; Figgen, D.; Goll, E.; Stoll, H.; Dolg, M. Systematically Convergent Basis Sets with Relativistic Pseudopotentials. II. Small-Core Pseudopotentials and Correlation Consistent Basis Sets for the Post-d Group 16−18 Elements. J. Chem. Phys. 2003, 119, 11113−11123. (42) Braams, B. J.; Bowman, J. M. Permutationally Invariant Potential Energy Surfaces in High Dimensionality. Int. Rev. Phys. Chem. 2009, 28, 577−606. (43) Xie, Z.; Bowman, J. M. Permutationally Invariant Polynomial Basis for Molecular Energy Surface Fitting via Monomial Symmetrization. J. Chem. Theory Comput. 2010, 6, 26−34. (44) Raghavachari, K.; Trucks, G. W.; Pople, J. A.; Head-Gordon, M. A Fifth-Order Perturbation Comparison of Electron Correlation Theories. Chem. Phys. Lett. 1989, 157, 479−483. (45) Watts, J. D.; Gauss, J.; Bartlett, R. J. Coupled-Cluster Methods with Noniterative Triple Excitations for Restricted Open-Shell Hartree-Fock and Other General Single Determinant Reference Functions. Energies and Analytical Gradients. J. Chem. Phys. 1993, 98, 8718−8733. (46) Dunning, T. H., Jr. Gaussian Basis Sets for Use in Correlated Molecular Calculations. I. The Atoms Boron Through Neon and Hydrogen. J. Chem. Phys. 1989, 90, 1007−1023. (47) Kendall, R. A.; Dunning, T. H., Jr.; Harrison, R. J. Electron Affinities of the First-Row Atoms Revisited. Systematic Basis Sets and Wave Functions. J. Chem. Phys. 1992, 96, 6796−6806.

(7) Ravishankara, A. R.; Wine, P. H.; Wells, J. R. The OH + HBr Reaction Revisited. J. Chem. Phys. 1985, 83, 447−448. (8) Ravishankara, A. R.; Wine, P. H.; Langford, A. O. Absolute Rate Constant For the Reaction OH + HBr → H2O + Br. Chem. Phys. Lett. 1979, 63, 479−484. (9) Sims, I. R.; Smith, I. W. M.; Clary, D. C.; Bocherel, P.; Rowe, B. R. Ultra-Low Temperature Kinetics of Neutral−Neutral Reactions New Experimental and Theoretical Results For OH + HBr Between 295 and 23 K. J. Chem. Phys. 1994, 101, 1748−1751. (10) Atkinson, D. B.; Jaramillo, V. I.; Smith, M. A. Low-Temperature Kinetic Behavior of the Bimolecular Reaction OH + HBr (76−242 K). J. Phys. Chem. A 1997, 101, 3356−3359. (11) Jaramillo, V. I.; Smith, M. A. Temperature-Dependent Kinetic Isotope Effects in the Gas-Phase Reaction: OH + HBr. J. Phys. Chem. A 2001, 105, 5854−5859. (12) Jaramillo, V. I.; Gougeon, S.; Le Picard, S. D.; Canosa, A.; Smith, M. A.; Rowe, B. R. A Consensus View of the Temperature Dependence of the Gas Phase Reaction: OH + HBr → H2O + Br. Int. J. Chem. Kinet. 2002, 34, 339−344. (13) Mullen, C.; Smith, M. A. Temperature Dependence and Kinetic Isotope Effects for the OH + HBr Reaction and H/D Isotopic Variants at Low Temperatures (53−135 K) Measured Using a Pulsed Supersonic Laval Nozzle Flow Reactor. J. Phys. Chem. A 2005, 109, 3893−3902. (14) Bedjanian, Y.; Riffault, V.; Le Bras, G.; Poulet, G. Kinetic Study of the Reactions of OH and OD with HBr and DBr. J. Photochem. Photobiol. A 1999, 128, 15−25. (15) Smith, I. W. M. Reactions at Very Low Temperatures: Gas Kinetics at a New Frontier. Angew. Chem., Int. Ed. 2006, 45, 2842− 2861. (16) Sabbah, H.; Biennier, L.; Sims, I. R.; Georgievskii, Y.; Klippenstein, S. J.; Smith, I. W. M. Understanding Reactivity at Very Low Temperatures: The Reactions of Oxygen Atoms with Alkenes. Science 2007, 317, 102−105. (17) Clary, D. C.; Stoecklin, T. S.; Wickham, A. G. Rate Constants for Chemical Reactions of Radicals at Low Temperatures. J. Chem. Soc., Faraday Trans. 1993, 89, 2185−2191. (18) Clary, D. C.; Nyman, G.; Hernandez, R. Mode Selective Chemistry in the Reactions of OH with HBr and HCl. J. Chem. Phys. 1994, 101, 3704−3714. (19) Nizamov, B.; Setser, D. W.; Wang, H. B.; Peslherbe, G. H.; Hase, W. L. Quasiclassical Trajectory Calculations for the OH(X2Π) and OD(X2Π) + HBr Reactions: Energy Partitioning and Rate Constants. J. Chem. Phys. 1996, 105, 9897−9911. (20) Rougeau, N.; Nyman, G.; Kubach, C. Quantum Studies of H2O + Cl → HO + HCl and H2O + Br → HO + HBr Reactions. A Comparison of Two Reduced Dimensionality Approaches. Phys. Chem. Chem. Phys. 1999, 1, 1191−1196. (21) Liu, J. Y.; Li, Z. S.; Dai, Z. W.; Huang, X. R.; Sun, C. C. Direct Ab Initio Dynamics Calculations of the Reaction Rates for the Hydrogen Abstraction OH + HBr → H2O + Br. J. Phys. Chem. A 2001, 105, 7707−7712. (22) Georgievskii, Y.; Klippenstein, S. J. Long-Range Transition State Theory. J. Chem. Phys. 2005, 122, 194103. (23) Li, G.; Zhou, L.; Li, Q.-S.; Xie, Y.; Schaefer, H. F., III. The Entrance Complex, Transition State, and Exit Complex for the F + H2O → HF + OH Reaction. Definitive Predictions. Comparison with Popular Density Functional Methods. Phys. Chem. Chem. Phys. 2012, 14, 10891−10895. (24) Li, J.; Dawes, R.; Guo, H. An Ab Initio Based Full-Dimensional Global Potential Energy Surface for FH2O(X2A′) and Dynamics for the F + H2O → HF + HO Reaction. J. Chem. Phys. 2012, 137, 094304. (25) Li, J.; Jiang, B.; Guo, H. Spin-Orbit Corrected Full-Dimensional Potential Energy Surfaces for the Two Lowest-Lying Electronic States of FH2O and Dynamics for the F + H2O → HF + OH Reaction. J. Chem. Phys. 2013, 138, 074309. (26) Li, J.; Jiang, B.; Guo, H. Enhancement of Bimolecular Reactivity by a Pre-Reaction van der Waals Complex: The Case of F + H2O → HF + HO. Chem. Sci. 2013, 4, 629−632. 711

dx.doi.org/10.1021/jz5000325 | J. Phys. Chem. Lett. 2014, 5, 706−712

The Journal of Physical Chemistry Letters

Letter

(48) Peterson, K. A.; Shepler, B. C.; Figgen, D.; Stoll, H. On the Spectroscopic and Thermochemical Properties of ClO, BrO, IO, and Their Anions. J. Phys. Chem. A 2006, 110, 13877−13883. (49) Peterson, K. A.; Feller, D.; Dixon, D. A. Chemical Accuracy in Ab Initio Thermochemistry and Spectroscopy: Current Strategies and Future Challenges. Theor. Chem. Acc. 2012, 131, 1079. (50) de Oliveira-Filho, A. G. S.; Ornellas, F. R.; Peterson, K. A. Accurate Ab Initio Potential Energy Surfaces for the 3A″ and 3A′ Electronic States of the O(3P) + HBr System. J. Chem. Phys. 2012, 136, 174316. (51) de Oliveira-Filho, A. G. S.; Ornellas, F. R.; Peterson, K. A.; Mielke, S. L. Thermal Rate Constants for the O(3P) + HBr and O(3P) + DBr Reactions: Transition-State Theory and Quantum Mechanical Calculations. J. Phys. Chem. A 2013, 117, 12703−12710. (52) Helgaker, T.; Klopper, W.; Koch, H.; Noga, J. Basis-Set Convergence of Correlated Calculations on Water. J. Chem. Phys. 1997, 106, 9639−9646. (53) Noga, J.; Bartlett, R. J. The Full CCSDT Model for Molecular Electronic Structure. J. Chem. Phys. 1987, 86, 7041−7050. (54) Scuseria, G. E.; Schaefer, H. F., III. A New Implementation of the Full CCSDT Model for Molecular Electronic Structure. Chem. Phys. Lett. 1988, 152, 382−386. (55) Kállay, M.; Surján, P. R. Higher Excitations in Coupled-Cluster Theory. J. Chem. Phys. 2001, 115, 2945−2954. (56) Kucharski, S. A.; Bartlett, R. J. The Coupled-Cluster Single, Double, Triple, and Quadruple Excitation Method. J. Chem. Phys. 1992, 97, 4282−4288. (57) Kucharski, S. A.; Bartlett, R. J. Recursive Intermediate Factorization and Complete Computational Linearization of the Coupled-Cluster Single, Double, Triple, and Quadruple Excitation Equations. Theor. Chim. Acta 1991, 80, 387−405. (58) Oliphant, N.; Adamowicz, L. Coupled-Cluster Method Truncated at Quadruples. J. Chem. Phys. 1991, 95, 6645−6651. (59) Peterson, K. A.; Dunning, T. H., Jr. Accurate Correlation Consistent Basis Sets for Molecular Core-Valence Correlation Effects: The Second Row Atoms Al−Ar, and the First Row Atoms B−Ne Revisited. J. Chem. Phys. 2002, 117, 10548−10560. (60) Peterson, K. A.; Yousaf, K. E. Molecular Core-Valence Correlation Effects Involving the Post-d Elements Ga−Rn: Benchmarks and New Pseudopotential-Based Correlation Consistent Basis Sets. J. Chem. Phys. 2010, 133, 174116. (61) Douglas, M.; Kroll, N. M. Quantum Electrodynamical Corrections to the Fine Structure of Helium. Ann. Phys. 1974, 82, 89−155. (62) Wilson, A. K.; Woon, D. E.; Peterson, K. A.; Dunning, T. H., Jr. Gaussian Basis Sets for Use in Correlated Molecular Calculations. IX. The Atoms Gallium Through Krypton. J. Chem. Phys. 1999, 110, 7667−7676. (63) de Jong, W. A.; Harrison, R. J.; Dixon, D. A. Parallel Douglas− Kroll Energy and Gradients in NWChem: Estimating Scalar Relativistic Effects Using Douglas−Kroll Contracted Basis Sets. J. Chem. Phys. 2001, 114, 48−53. (64) Berning, A.; Schweizer, M.; Werner, H.-J.; Knowles, P. J.; Palmieri, P. Spin-Orbit Matrix Elements for Internally Contracted Multireference Configuration Interaction Wavefunctions. Mol. Phys. 2000, 98, 1823−1833. (65) Deskevich, M. P.; Nesbitt, D. J.; Werner, H.-J. Dynamically Weighted Multiconfiguration Self-Consistent Field: Multistate Calculations for F + H2O → HF + OH Reaction Paths. J. Chem. Phys. 2004, 120, 7281−7289. (66) Werner, H.-J.; Knowles, P. J. A Second Order Multiconfiguration SCF Procedure with Optimum Convergence. J. Chem. Phys. 1985, 82, 5053−5063. (67) Knowles, P. J.; Werner, H.-J. An Efficient Second-Order MCSCF Method for Long Configuration Expansions. Chem. Phys. Lett. 1985, 115, 259−267. (68) Werner, H.-J.; Knowles, P. J.; Manby, F. R.; Schütz, M.; Celani, P.; Knizia, G.; Korona, T.; Lindh, R.; Mitrushenkov, A.; Rauhut, G.; et al. Molpro, version 2010.1, a package of ab initio programs, 2010.

(69) Werner, H.-J.; Knowles, P. J.; Knizia, G.; Manby, F. R.; Schütz, M. Molpro: A General-Purpose Quantum Chemistry Program Package. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2012, 2, 242−253. (70) MRCC, a string-based quantum chemical program written by M. Kállay. See www.mrcc.hu. (71) Truhlar, D. G.; Muckerman, J. T. In Atom-Molecule Collision Theory: A Guide for the Experimentalist; Bernstein, R. B., Ed.; Plenum Press: New York, 1979; pp 505−566. (72) Peslherbe, G. H.; Wang, H.; Hase, W. L. In Monte Carlo Methods in Chemical Physics; Prigogine, I., Rice, S. A., Eds.; Advances in Chemical Physics; Wiley, 2007; pp 171−201. (73) Hu, X.; Hase, W. L.; Pirraglia, T. Vectorization of the General Monte Carlo Classical Trajectory Program VENUS. J. Comput. Chem. 1991, 12, 1014−1024. (74) Hase, W. L.; Duchovic, R. J.; Hu, X.; Komornicki, A.; Lim, K. F.; Lu, D.-H.; Peslherbe, G. H.; Swamy, K. N.; Linde, S. R. V.; Varandas, A. J. C.; Wang, H.; Wolf, R. J. Quantum Chem. Program Exch. Bull. 1996, 16, 671. (75) Herzberg, G. Molecular Spectra and Molecular Structure I: Spectra of Diatomic Molecules; Van Nostrand Reinhold: New York, 1950; pp 220. (76) McQuarrie, D. A. Statistical Mechanics; University Science Books: Sausalito, CA, 2000; pp 109. (77) Maillard, J. P.; Chauville, J.; Mantz, A. W. High-Resolution Emission-Spectrum of OH in an Oxyacetylene Flame from 3.7 to 0.9 μm. J. Mol. Spectrosc. 1976, 63, 120−141. (78) Truhlar, D. G. Multiple Potential Energy Surfaces for Reactions of Species in Degenerate Electronic States. J. Chem. Phys. 1972, 56, 3189−3190. (79) Meyer, H.-D.; Miller, W. H. Classical Models for Electronic Degrees of Freedom: Derivation Via Spin Analogy and Application to F* + H2 → F + H2. J. Chem. Phys. 1979, 71, 2156−2169. (80) Tully, J. C. Molecular Dynamics with Electronic Transitions. J. Chem. Phys. 1990, 93, 1061−1071. (81) Nguyen, T. L.; Stanton, J. F. Ab Initio Thermal Rate Calculations of HO + HO = O(3P) + H2O Reaction and Isotopologues. J. Phys. Chem. A 2013, 117, 2678−2686. (82) Altinay, G.; Macdonald, R. G. Determination of the Rate Constant for the OH(X2Π) + OH(X2Π) → H2O + O(3P) Reaction Over the Temperature Range 295 to 701 K. J. Phys. Chem. A 2013, DOI: 10.1021/jp409344q.

712

dx.doi.org/10.1021/jz5000325 | J. Phys. Chem. Lett. 2014, 5, 706−712