A theoretical study of the H2SO4+H2O → HSO4 −+H3O+ reaction at the surface of aqueous aerosols

The surface region of sulfate aerosols (supercooled aqueous concentrated sulfuric acid solutions) is the likely site of a number of important heterogeneous reactions in various locations in the atmosphere, but the surface region ionic composition is not known. As a first step in exploring this issue, the first acid ionization reaction for sulfuric acid, H₂SO₄ + H₂O HSO₄^– + H₃O⁺, is studied via electronic structure calculations at the Hartree–Fock level on an H₂SO₄ molecule embedded in the surface region of a cluster containing 33 water molecules. An initial H₂SO₄ configuration is selected which could produce H₃O⁺ readily available for heterogeneous reactions, but which involves reduced solvation and is consistent with no dangling OH bonds for H₂SO₄. It is found that at 0 K and with zero-point energy included, the proton transfer is endothermic by 3.4 kcal/mol. This result is discussed in the context of reactions on sulfate aerosol surfaces and, further, more complex calculations.Contribution to the Jacopo Tomasi Honorary Issue     

Dynamical studies of the reaction H~++He+He→HHe~++He

The classical trajectory method is employed to calculate the rate coefficient kr for the reaction H++He+He-HHe++He at the temperatures ranging 200-350 K, based on an ab initio potential energy surface. The results show that kr is strongly dependent on the temperature, which can be well fitted by the function kr=ATDr-3 with A=4.192x10-31 cm6/s and the reaction dimension Dr=2.706. The product molecules HHe+ are found in high vibrational states.     

Comment on the possible role of the reaction H+H2O→H2+OH in the radiolysis of water at high temperatures

In a recent paper (Radiation Physics and Chemistry, 2005, vol. 74, pp. 210) it was suggested that the anomalous increase of molecular hydrogen radiolysis yields observed in high-temperature water is explained by a high activation energy for the reaction H+H₂O→H₂+OH. In this comment we present thermodynamic arguments to demonstrate that this reaction cannot be as fast as suggested. A best estimate for the rate constant is 2.2×10³ M⁻¹ s⁻¹ at 300 °C. Central to this argument is an estimate of the OH radical hydration free energy vs. temperature, ΔG_hyd(OH)=0.0278t−18.4 kJ/mole (t in °C, equidensity standard states), which is based on analogy with the hydration free energy of water and of hydrogen peroxide.     

Mean potential statistical theory of the H+ + D2 → HD + D+ reaction

The reactive collision process H(+) + D(2)(ν = 0, j = 0) → HD + D(+) is theoretically analyzed for collision energies ranging from threshold up to 1.3 eV. It is assumed that the reaction takes place via formation of a collision complex. In calculations, a statistical theory is used, based on a mean isotropic potential deduced from a full potential energy surface. Calculated integral cross sections, opacity functions, and rotational distributions of the HD products are compared with recent statistical and quantum mechanical calculations performed using a full potential energy surface. Satisfactory agreement between the results obtained using the two statistical methods is found, both of which however overestimate the existing quantum mechanical predictions. The effects due to the presence of identical particles are also discussed.     

Reply to comment on the possible role of the reaction H+H2O→H2+OH in the radiolysis of water at high temperatures

In reply to “Comment on the possible role of reaction H+H₂O→H₂+OH in the radiolysis of water at high temperatures” (Bartels, 2009 Comment on the possible role of the reaction H+H₂O→H₂+OH in the radiolysis of water at high temperatures. Radiat. Phys. Chem. 78, 191–194) we present an alternative thermodynamic estimation of the reaction rate constant k. Based on the non-symmetric standard state convention we have calculated that the Gibbs energy of reaction Δ_rG=57.26 kJ mol^?1 and the reaction rate constant k=7.23×10^?5 M^?1 s^?1 at ambient temperature. Re-analysis of the thermodynamic estimation (Bartels, 2009 Comment on the possible role of the reaction H+H₂O→H₂+OH in the radiolysis of water at high temperatures. Radiat. Phys. Chem. 78, 191–194) showed that the upper limit for the rate constant at 573 K is k=1.75×10⁴ M^?1 s^?1 compared to the value predicted by the diffusion-kinetic modelling (3.18±1.25)×10⁴ M^?1 s^?1 (Swiatla-Wojcik, D., Buxton, G.V., 2005. On the possible role of the reaction H+H₂O→H₂+OH in the radiolysis of water at high temperatures. Radiat. Phys. Chem. 74(3–4), 210–219). The presented thermodynamic evaluation of k(573) is based on the assumption that k can be calculated from Δ_rG and the rate constant of the reverse reaction which, as discussed, are both uncertain at high temperatures.     

Acceleration of the reaction OH + CO → H + CO2 by vibrational excitation of OH

The collision complex formed from a vibrationally excited reactant undergoes redissociation to the reactant, intramolecular vibrational relaxation (randomization of vibrational energy), or chemical reaction to the products. If attractive interaction between the reactants is large, efficient vibrational relaxation in the complex prevents redissociation to the reactants with the initial vibrational energy, and the complex decomposes to the reactants with low vibrational energy or converts to the products. In this paper, we have studied the branching ratios between the intramolecular vibrational relaxation and chemical reaction of an adduct HO(v)-CO formed from OH(X(2)Π(i)) in different vibrational levels v = 0-4 and CO. OH(v = 0-4) generated in a gaseous mixture of O(3)/H(2)/CO/He irradiated at 266 nm was detected with laser-induced fluorescence (LIF) via the A(2)Σ(+)-X(2)Π(i) transition, and H atoms were probed by the two-photon excited LIF technique. From the kinetic analysis of the time-resolved LIF intensities of OH(v) and H, we have found that the intramolecular vibrational relaxation is mainly governed by a single quantum change, HO(v)-CO → HO(v-1)-CO, followed by redissociation to OH(v-1) and CO. With the vibrational quantum number v, chemical process from the adduct to H + CO(2) is accelerated, and vibrational relaxation is decelerated. The countertrend is elucidated by the competition between chemical reaction and vibrational relaxation in the adduct HOCO.     

Quantum dynamics of the S+OH→SO+H reaction

First accurate quantum mechanical scattering calculations have been carried out for the S((3)P)+OH(X?(2)Π)→SO(X?(3)Σ(-))+H((2)S) reaction using a recent ab initio potential energy surface for the ground electronic state, X?(2)A("), of HSO. Total and state-to-state reaction probabilities for a total angular momentum J=0 have been determined for collision energies up to 0.5 eV. A rate constant has been calculated by means of the J-shifting approach in the 10-400 K temperature range. Vibrational and rotational product distributions show no specific behavior and are consistent with a mixture of direct and indirect reaction mechanisms.     

Quasi-classical trajectory studies of the stereodynamics of the reaction O + HCl → ClO + H

The stereodynamics of the O + HCl → ClO + H reaction are investigated by quasi-classical trajectory (QCT) method. The calculations are carried out on the ground 1 ¹ A′ potential energy surface (PES). The orientation and alignments of the product rotational angular momentum for the title reaction are reported. The influence of collision energy on the product vector properties is also studied in the present work. Four (2π/σ)(dσ₀₀/dω _t ), (2π/σ)(dσ₂₀/dω _t ), (2π/σ)(dσ₂₂₊/dω _t ), and (2π / σ)(dσ₂₁₋/dω _t ), and have been calculated in the center of mass frame.     

Comparative study of the photoassociation reaction of He + X+→ HeX+ (X = H,D): Including multiphoton transitions and dissociations

Using the time-dependent quantum wave packet method, the photoassociation (PA) processes of He + H ⁺→ HeH⁺ and He + D ⁺→ HeD⁺, driven by the sin²-shaped femtosecond laser pulse in the electronic ground state, including multiphoton transitions and dissociations, are investigated for a wide range of initial collision momenta spanning from 1 to 4 a.u. (or for the collision energy roughly in the ranges of 0.009∼0.148 eV and 0.006∼0.089 eV for HeH⁺ and HeD⁺ systems, respectively). It is found that, at some collision momenta, multiphoton transitions to deeply bound states are inevitable to occur and can greatly decrease the PA probability of the target state that selected is the vibrational state v = 6. For the dissociation process, the higher-order (two- and three-photon) dissociations, measured from the target state, tend to be significant at relative high collision energies, which implies that above-threshold dissociations may also be an important loss mechanism in the PA process. In addition, it is also shown that the higher-order dissociation is much stronger for HeH⁺ systems than that for HeD⁺ systems at a given collision momentum, and could be enhanced by the strong transitions among deeply bound states.     

Kinetics and dynamics of the NH3 + H → NH2 + H2 reaction using transition state methods, quasi-classical trajectories, and quantum-mechanical scattering

On a recent analytical potential energy surface developed by two of the authors, an exhaustive kinetics study, using variational transition state theory with multidimensional tunneling effect, and dynamics study, using both quasi-classical trajectory and full-dimensional quantum scattering methods, was carried out to understand the reactivity of the NH(3) + H → NH(2) + H(2) gas-phase reaction. Initial state-selected time-dependent wave packet calculations using a full-dimensional model were performed, where the total reaction probabilities were calculated for the initial ground vibrational state and for four excited vibrational states of ammonia. Thermal rate constants were calculated for the temperature range 200-2000 K using the three methods and compared with available experimental data. We found that (a) the total reaction probabilities are very small, (b) the symmetric and asymmetric N-H stretch excitations enhance the reactivity, (c) the quantum-mechanical calculated thermal rate constants are about one order of magnitude smaller than the transition state theory results, which reproduce the experimental evidence, and (d) quasi-classical trajectory calculations, which were performed with the main goal of analyzing the influence of the zero-point energy problem on the final dynamics results, reproduce the quantum scattering calculations on the same surface.     

Theoretical studies on the reaction path and dynamics of the reaction CH_3+H_2O→CH_4+OH

The reaction path, the dynamical properties along the reaction path and CVT rate constants are computed by the ab initio MO method, the reaction path Hamiltonian theory and the variational transition state theory. The results show that the effect of the electron correlation energy on activation barrier is large, the recrossing and tunneling effects exist in the reaction.     

Ab initio study of the reaction path of the reaction HNCO+CN

HNCO is a convenient photolytic source of NCO and NH radicals for laboratory kinetics studies of elementary reaction[1] and plays an important role in the combustion and atmosphere chemistry. It can re- move deleterious compounds rapidly from exhausted ga…     

Rate constant for the reaction C2H5 + HBr → C2H6 + Br

RRKM theory has been employed to analyze the kinetics of the title reaction, in particular, the once-controversial negative activation energy. Stationary points along the reaction coordinate were characterized with coupled cluster theory combined with basis set extrapolation to the complete basis set limit. A shallow minimum, bound by 9.7 kJ?mol(-1) relative to C(2)H(5) + HBr, was located, with a very small energy barrier to dissociation to Br + C(2)H(6). The transition state is tight compared to the adduct. The influence of vibrational anharmonicity on the kinetics and thermochemistry of the title reaction were explored quantitatively. With adjustment of the adduct binding energy by ～4 kJ?mol(-1), the computed rate constants may be brought into agreement with most experimental data in the literature, including new room-temperature results described here. There are indications that at temperatures above those studied experimentally, the activation energy may switch from negative to positive.     

Dynamics and mode-selected reaction of NCO+H_2→HNCO+H

The reaction path of the reaction NCO+H2→HNCO + H has been traced by Fukui's theory and the ab initio method. On this basis, the dynamical properties along the reaction path, canonical variational theory (CVT) rate constants and vibrational-mode-selected rate constants have been computed. The results show that the effect of the electron correlation energy on the activation barrier is large, and tiros the correction by MP4 method is effective; the results also show that the recrossing and tunneling effects exist, and thus the corrections by the variational transition state theory (VTST) and the small curvature (SC) approximation method are also effective. In the reaction, the coupling and energy transfer between mode 8(7) and reaction path are strong, so the rate is effectively enhanced while these modes, especially H2 stretching, are vibrationally excited.     

Ab initio potential energy surface and quantum dynamics for the H + CH4 → H2 + CH3 reaction

A new full-dimensional potential energy surface for the title reaction has been constructed using the modified Shepard interpolation scheme. Energies and derivatives were calculated using the UCCSD(T) method with aug-cc-pVTZ and 6-311++G(3df,2pd) basis sets, respectively. A total number of 30,000 data points were selected from a huge number of molecular configurations sampled by trajectory method. Quantum dynamical calculations showed that the potential energy surface is well converged for the number of data points for collision energy up to 2.5 eV. Total reaction probabilities and integral cross sections were calculated on the present surface, as well as on the ZBB3 and EG-2008 surfaces for the title reaction. Satisfactory agreements were achieved between the present and the ZBB3 potential energy surfaces, indicating we are approaching the final stage to obtain a global potential energy surface of quantitative accuracy for this benchmark polyatomic system. Our calculations also showed that the EG-2008 surface is less accurate than the present and ZBB3 surfaces, particularly in high energy region.     

Ab initio study on the dynamical properties of the hydrogen abstraction reaction NH2 + OH → NH + H2O

 The geometry of the transition state of the title reaction was optimized at the unrestricted Hartree–Fock, the spin-unrestricted second-order M?ller–Plesset, and the spin-unrestricted quadratic configuration interaction with all single and double substitutions levels of theory. The changes in the geometry, the bound vibrational modes, and the potential energy along the minimum energy path are discussed. Variational transition-state theory rate constants calculated with the tunneling and curvature effect correction agree very well with the experimental values. Received: 23 April 1999 / Accepted: 9 June 1999 / Published online: 15 December 1999     

Quasi-classical trajectory study of the dynamics of the Cl + CH4→ HCl + CH3 reaction

We present an on-the-fly classical trajectory study of the Cl + CH(4)→ HCl + CH(3) reaction using a specific reaction parameter (SRP) AM1 Hamiltonian that was previously optimized for the Cl + ethane reaction [S. J. Greaves et al., J. Phys Chem A, 2008, 112, 9387]. The SRP-AM1 Hamiltonian is shown to be a good model for the potential energy surface of the title reaction. Calculated differential cross sections, obtained from trajectories propagated with the SRP-AM1 Hamiltonian compare favourably with experimental results for this system. Analysis of the vibrational modes of the methyl radical shows different scattering distributions for ground and vibrationally excited products.     

Quantum instanton calculation of rate constants for the C2H6 + H → C2H5 + H2 reaction: anharmonicity and kinetic isotope effects

Thermal rate constants and kinetic isotope effects for the title reaction are calculated by using the quantum instanton approximation within the full dimensional Cartesian coordinates. The obtained results are in good agreement with experimental measurements at high temperatures. The detailed investigation reveals that the anharmonicity of the hindered internal rotation motion does not influence the rate too much compared to its harmonic oscillator approximation. However, the motion of the nonreactive methyl group in C(2)H(6) significantly enhances the rates compared to its rigid case, which makes conventional reduced-dimensionality calculations a challenge. In addition, the temperature dependence of kinetic isotope effects is also revealed.