Theoretical study of the dynamics of the H+CH4 and H+C2H6 reactions using a specific-reaction-parameter semiempirical Hamiltonian

We present a theoretical study of the reactions of hydrogen atoms with methane and ethane molecules and isotopomers. High-accuracy electronic-structure calculations have been carried out to characterize representative regions of the potential-energy surface (PES) of various reaction pathways, including H abstraction and H exchange. These ab initio calculations have been subsequently employed to derive an improved set of parameters for the modified symmetrically-orthogonalized intermediate neglect of differential overlap (MSINDO) semiempirical Hamiltonian, which are specific to the H+alkane family of reactions. The specific-reaction-parameter (SRP) Hamiltonian has then been used to perform a quasiclassical-trajectory study of both the H+CH4 and H+C2H6 reactions. The calculated values of dynamics properties of the H+CH4-->H2+CH3 reaction and isotopologues, including alkyl product speed distributions, diatomic product internal-state distributions, and cross sections, are generally in good agreement with experiment and with the results provided by the ZBB3 PES [Z. Xie et al., J. Chem. Phys. 125, 133120 (2006)]. The results of trajectories propagated with the SRP Hamiltonian for the H+C2H6-->H2+C2H5 reaction also agree with experiment. The level of agreement between the results calculated with the SRP Hamiltonian and experiment in both the H+methane and H+ethane reactions indicates that semiempirical Hamiltonians can be improved for not only a specific reaction but also a family of reactions.     

Quasiclassical trajectory study of the O(3P) + CH4 --> OH + CH3 reaction with a specific reaction parameters semiempirical Hamiltonian

We present a theoretical study of the O(3P) + CH4 --> OH + CH3 reaction using electronic structure, kinetics, and dynamics calculations. We calculate a grid of ab initio points at the PMP2/AUG-cc-pVDZ level to characterize the potential energy surface in regions of up to 1.3 eV above reagents. This grid of ab initio points is used to derive a set of specific reaction parameters (SRP) for the MSINDO semiempirical Hamiltonian. The resulting SRP-MSINDO Hamiltonian improves the quality of the standard Hamiltonian, particularly in regions of the potential energy surface beyond the minimum-energy reaction path. Quasiclassical-trajectory calculations are used to study the reaction dynamics with the original and the improved MSINDO semiempirical Hamiltonians, and a prior surface. The SRP-MSINDO semiempirical Hamiltonian yields OH rotational distributions in agreement with experimental results, improving over the results of the other surfaces. Thermal rate constants estimated with Variational Transition State Theory using the SRP-MSINDO Hamiltonian are also in agreement with experiments. Our results indicate that reparametrized semiempirical Hamiltonians are a good alternative to generating potential energy surfaces for accurate dynamics studies of polyatomic reactions.     

Ab initio and direct quasiclassical-trajectory study of the Cl + CH4-->HCl + CH3 reaction

We present an electronic structure and dynamics study of the Cl + CH(4)--> HCl + CH(3) reaction. We have characterized the stationary points of the ground-state potential-energy surface using various electronic structure methods and basis sets. Our best calculations, CCSD(T) extrapolated to the complete basis-set limit based on geometries and harmonic frequencies obtained at the CCSD(T)/aug-cc-pvtz level, are in agreement with the experimental reaction energy and indirect measurements of the barrier height. Using ab initio information, we have reparametrized a semiempirical Hamiltonian so that the predictions of the improved Hamiltonian agree with the higher-level calculations in various regions of the potential-energy surface. This improved semiempirical Hamiltonian is then used to propagate quasiclassical trajectories and characterize the reaction dynamics. The good agreement of the calculated HCl rotational and angular distributions with the experiment indicates that reparametrizing semiempirical Hamiltonians is a promising approach to derive accurate potential-energy surfaces for polyatomic reactions. However, excessive energy leakage from the initial vibrational energy of the CH(4) molecule to the reaction coordinate in the trajectory calculations calls into question the suitability of the standard quasiclassical-trajectory method to describe energy partitioning in polyatomic reactions.     

Ab initio and direct quasiclassical-trajectory study of the F+CH4-->HF+CH3 reaction

We present an electronic structure and dynamics study of the F+CH4-->HF+CH3 reaction. CCSD(T)/aug-cc-pVDZ geometry optimizations, harmonic-frequency, and energy calculations indicate that the potential-energy surface is remarkably isotropic near the transition state. In addition, while the saddle-point F-H-C angle is 180 degrees using MP2 methods, CCSD(T) geometry optimizations predict a bent transition state, with a 153 degrees F-H-C angle. We use these high-quality ab initio data to reparametrize the parameter-model 3 (PM3) semiempirical Hamiltonian so that calculations with the improved Hamiltonian and employing restricted open-shell wave functions agree with the higher accuracy data. Using this specific-reaction-parameter PM3 semiempirical Hamiltonian (SRP-PM3), we investigate the reaction dynamics by propagating quasiclassical trajectories. The results of our calculations using the SRP-PM3 Hamiltonian are compared with experiments and with the estimates of two recently reported potential-energy surfaces. The trajectory calculations using the SRP-PM3 Hamiltonian reproduce quantitatively the measured HF vibrational distributions. The calculations also agree with the experimental HF rotational distributions and capture the essential features of the excitation function. The results of the SRP semiempirical Hamiltonian developed here clearly improve over those using the two prior potential-energy surfaces and suggest that reparametrization of semiempirical Hamiltonians is a promising strategy to develop accurate potential-energy surfaces for reaction dynamics studies of polyatomic systems.     

Quantum dynamics study on the product branching for the C(3P) + C2H2 reaction: cyclic-C3H versus linear-C3H

Reduced-dimensionality quantum reactive scattering calculations for the C(3P) + C2H2 reaction have been carried out in order to understand the product branching dynamics of cyclic-C3H + H and linear-C3H + H. Our model treats only two degrees of freedom but can explicitly describe both of the C3H isomer product channels. The lowest triplet potential energy surface has been obtained by the hybrid density-functional method at the B3LYP/6-31G(d,p) level of theory. The calculated reaction probabilities were found to be dominated by resonance consistent with the complex-formation potential, and the results show that cyclic-C3H is preferentially formed via the cyclic-C3H2 intermediate produced by insertion of C(3P) into the CC bond. We have found that the isomerization from the cyclic-C3H2 to linear-C3H2 intermediate is suppressed by a barrier separating potential wells corresponding to these two intermediates. It has also been found that the energy dependence of the calculated total reaction cross section is in good agreement with the result of crossed molecular beam experiments.     

Exact quantum and TST-CEQ study of the collinear reaction O+HCl

Exact quantum calculations of reaction probabilities have been carried out using hyperspherical coordinates for the collinearr reaction O+HCl(v <1) -OH(v'<1)+Cl . A generalized LEPS potential energy surface with a barrier height of 8.12 kcal/mol has been used in the calculations. According to the calculated results we found that (1) the reaction probability oscillates with energy, (2) the reaction probability shows vibrational adiabaticity, although it is poorer than that for symmetric reaction Cl + HC1. The analysis of resonance has also been done. The reaction rate constants and average cross sections have been calculated by TST-CEQ method. The rate constants are in agreement with that by QCT and smaller than the experimental one. Finally, the threshold has been estimated and is in good agreement with that of the literature.     

A reaction surface Hamiltonian study of malonaldehyde

We report calculations using a reaction surface Hamiltonian for which the vibrations of a molecule are represented by 3N-8 normal coordinates, Q, and two large amplitude motions, s(1) and s(2). The exact form of the kinetic energy operator is derived in these coordinates. The potential surface is first represented as a quadratic in Q, the coefficients of which depend upon the values of s(1),s(2) and then extended to include up to Q(6) diagonal anharmonic terms. The vibrational energy levels are evaluated by solving the variational secular equations, using a basis of products of Hermite polynomials and appropriate functions of s(1),s(2). Our selected example is malonaldehyde (N=9) and we choose as surface parameters two OH distances of the migrating H in the internal hydrogen transfer. The reaction surface Hamiltonian is ideally suited to the study of the kind of tunneling dynamics present in malonaldehyde. Our results are in good agreement with previous calculations of the zero point tunneling splitting and in general agreement with observed data. Interpretation of our two-dimensional reaction surface states suggests that the OH stretching fundamental is incorrectly assigned in the infrared spectrum. This mode appears at a much lower frequency in our calculations due to substantial transition state character.     

Communication: theoretical exploration of Au(+)+H2, D2, and HD reactive collisions

A quasi-classical study of the endoergic Au(+)((1)S) + H(2)(X(1)Σ(g)(+)) → AuH(+) ((2)Σ(+)) + H((2)S) reaction, and isotopic variants, is performed to compare with recent experimental results [F. Li, C. S. Hinton, M. Citir, F. Liu, and P. B. Armentrout, J. Chem. Phys. 134, 024310 (2011)]. For this purpose, a new global potential energy surface has been developed based on multi-reference configuration interaction ab initio calculations. The quasi-classical trajectory results show a very good agreement with the experiments, showing the same trends for the different isotopic variants of the hydrogen molecule. It is also found that the total dissociation into three fragments, Au(+)+H+H, is the dominant reaction channel for energies above the H(2) dissociation energy. This results from a well in the entrance channel of the potential energy surface, which enhances the probability of H-Au-H insertion.     

Tunneling in the loss of hydrogen chloride from isopropyl chloride cation

The unimolecular dissociation of isopropyl chloride cation has been investigated using mass-analyzed ion kinetic energy spectrometry. The C3H6*+ ion was the only product ion in the metastable dissociation. The kinetic energy release distribution for the HCl loss was determined. Ab initio molecular orbital calculations were performed at the MP2/6-311++G(d,p) level together with single point energy calculations at the QCSID(T)/6-311++G(2d,2p) level. The calculations show that the molecular ion rearranges to an ion-dipole complex prior to loss of HCl via a transition state containing a four-membered ring. The rearrangement involves H atom transfer. On the basis of the potential energy surface obtained for the loss of HCl and Cl*, the rate constants were calculated by transition-state statistical theories with considering tunneling effect. From the calculated result, it is proposed that the observed HCl loss would occur via tunneling through the barrier for isomerization to the ion-dipole complex, CH3CHCH2*+...HCl.     

Trajectory surface hopping study of the O(3P) + ethylene reaction dynamics

Spin-orbit coupling (SOC) induced intersystem crossing (ISC) has long been believed to play a crucial role in determining the product distributions in the O(3P) + C2H4 reaction. In this paper, we present the first nonadiabatic dynamics study of the title reaction at two center-of-mass collision energies: 0.56 eV, which is barely above the H-atom abstraction barrier on the triplet surface, and 3.0 eV, which is in the hyperthermal regime. The calculations were performed using a quasiclassical trajectory surface hopping (TSH) method with the potential energy surface generated on the fly at the unrestricted B3LYP/6-31G(d,p) level of theory. To simplify our calculations, nonadiabatic transitions were only considered when the singlet surface intersects the triplet surface. At the crossing points, Landau-Zener transition probabilities were computed assuming a fixed spin-orbit coupling parameter, which was taken to be 70 cm-1 in most calculations. Comparison with a recent crossed molecular beam experiment at 0.56 eV collision energy shows qualitative agreement as to the primary product branching ratios, with the CH3 + CHO and H + CH2CHO channels accounting for over 70% of total product formation. However, our direct dynamics TSH calculations overestimate ISC so that the total triplet/singlet ratio is 25:75, compared to the observed 43:57. Smaller values of SOC reduce ISC, resulting in better agreement with the experimental product relative yields; we demonstrate that these smaller SOC values are close to being consistent with estimates based on CASSCF calculations. As the collision energy increases, ISC becomes much less important and at 3.0 eV, the triplet to singlet branching ratio is 71:29. As a result, the triplet products CH2 + CH2O, H + CH2CHO and OH + C2H3 dominate over the singlet products CH3 + CHO, H2 + CH2CO, etc.     

Crossed beams and theoretical studies of the O(3P)+CH4-->H+OCH3 reaction excitation function

The excitation function for the reaction, O(3P)+CH4-->H+OCH3, has been measured in a crossed molecular beams experiment and determined with direct dynamics calculations that use the quasiclassical trajectory method in conjunction with a recently developed semiempirical Hamiltonian. Good agreement is found between experiment and theory, enabling us to address two fundamental issues for the O(3P)+CH4 reaction that arise for all O(3P)+saturated hydrocarbon reactions: (1) the importance of triplet excited states that correlate adiabatically to ground-state reactants and products and (2) the importance of intersystem crossing processes involving the lowest singlet surface [corresponding to reaction with O(1D)]. Our results indicate that the first excited triplet surface contributes substantially to the cross section when the collision energy exceeds the reaction barrier (approximately 2 eV) by more than 0.5 eV. Although triplet-singlet crossings may occur at all energies, we have found that their effect on the excitation function is negligible for the collision energies studied-up to 1.5 eV above threshold.     

State-to-state reactive differential cross sections for the H+H2-->H2+H reaction on five different potential energy surfaces employing a new quantum wavepacket computer code: DIFFREALWAVE

State-to-state differential cross sections have been calculated for the hydrogen exchange reaction, H+H2-->H2+H, using five different high quality potential energy surfaces with the objective of examining the sensitivity of these detailed cross sections to the underlying potential energy surfaces. The calculations were performed using a new parallel computer code, DIFFREALWAVE. The code is based on the real wavepacket approach of Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)]. The calculations are parallelized over the helicity quantum number Omega' (i.e., the quantum number for the body-fixed z component of the total angular momentum) and wavepackets for each J,Omega' set are assigned to different processors, similar in spirit to the Coriolis-coupled processors approach of Goldfield and Gray [Comput. Phys. Commun. 84, 1 (1996)]. Calculations for J=0-24 have been performed to obtain converged state-to-state differential cross sections in the energy range from 0.4 to 1.2 eV. The calculations employ five different potential energy surfaces, the BKMP2 surface and a hierarchical family of four new ab initio surfaces [S. L. Mielke, et al., J. Chem. Phys. 116, 4142 (2002)]. This family of four surfaces has been calculated using three different hierarchical sets of basis functions and also an extrapolation to the complete basis set limit, the so called CCI surface. The CCI surface is the most accurate surface for the H3 system reported to date. Our calculations of differential cross sections are the first to be reported for the A2, A3, A4, and CCI surfaces. They show that there are some small differences in the cross sections obtained from the five different surfaces, particularly at higher energies. The calculations also show that the BKMP2 performs well and gives cross sections in very good agreement with the results from the CCI surface, displaying only small divergences at higher energies.     

Planar study of H2O+Cl<-->HO+HCl reactions

The first four dimensional (4D) quantum scattering calculations on the tetra-atomic H2O+Cl<-->HO+HCl reactions are reported. With respect to a full (6D) treatment, only the planar constraint and a fixed length for the HO spectator bond are imposed. This work explicitly accounts for the bending and local HO stretching vibrations in H2O, for the vibration of HCl and for the in-plane rotation of the H2O, HO and HCl molecules. The calculations are performed with the potential energy surface of Clary et al. and use a Born-Oppenheimer type separation between the motions of the light and the heavy nuclei. State-to-state cross sections are reported for a collision energy range 0-1.8 eV measured with respect to H2O+Cl. For the H2O+Cl reaction, present results agree with previous (3D) non planar calculations and confirm that excitation of the H2O stretching promotes more reactivity than excitation of the bending. New results are related to the rotation of the H2O molecule: the cross sections are maximal for planar rotational states corresponding to 10相似文献   

New ab initio potential energy surface for BrH2 and rate constants for the H + HBr → H2 + Br abstraction reaction

A global potential energy surface (PES) for the electronic ground state of the BrH(2) system was constructed based on the multireference configuration interaction (MRCI) method including the Davidson's correction using a large basis set. In addition, the spin-orbit correction were computed using the Breit-Pauli Hamiltonian and the unperturbed MRCI wavefunctions in the Br + H(2) channel and the transition state region. Adding the correction to the ground state potential, the lowest spin-orbit correlated adiabatic potential was obtained. The characters of the new potential are discussed. Accurate initial state specified rate constants for the H + HBr → H(2) + Br abstraction reaction were calculated using a time-dependent wave packet method. The predicted rate constants were found to be in excellent agreement with the available experimental values and much better than those obtained from a previous PES.     

Cl+HD反应产物的极化与态分布

利用准经典轨线理论,在BW2和G3两个势能面上,研究了Cl+HD反应的动力学.计算结果表明,产物的转动取向对势能面及反应体系的质量因子非常敏感.在BW2势能面上,计算的两个产物的转动取向强于在G3势能面上计算的结果,而无论是在BW2势能面上还是在G3势能面上,DCl产物的取向都强于HCl产物的取向.计算结果还表明,在不同的势能面上反应物的转动激发对反应的影响有着显著的不同.在BW2势能面上,反应物的初始转动激发有利于Cl+HD反应的进行;而在G3势能面上,反应物的初始转动激发消弱了反应的反应性.     

CH2Cl与OH自由基反应机理的理论研究

The reaction mechanism of CH2Cl radical with OH radical to produce HCCl+H2O,HCOCl+H2 and H2CO+HCl has been studied by using quantum chemistry ab initio calculations. The optimized geometrical parameters,and vibrational frequencies of all species were obtained at the UMP2（FC）level of theory in conjunction with 6-311++G* basis set. Besides,the zero-point energies（ZPE）,relative energies and total energies of all species were calculated using Gaussian-3（G3）model. The results of theoretical study indicate that the activated intermediate CH2ClOH is first formed through a barrierless process,followed by atoms migration,radical groups rotation and bonds fission to produce HCCl+H2O,HCOCl+H2 and H2CO+HCl,respectively. And all channels are exothermic by 72.81,338.54 and 354.08 kJ/mol. The reaction heat of reactants to H2CO+HCl is 281.27 kJ/mol more than that of reactants to HCCl+H2O. This result accords with that of experiments.     

Theoretical study of the decomposition of BCl3 induced by a H radical

We report on a theoretical study of the gas-phase decomposition of boron trichloride in the presence of hydrogen radicals using ab initio energetic calculations coupled to TST, RRKM, and VTST-VRC kinetic calculations. In particular, we present an addition-elimination mechanism (BCl(3) + H → BHCl(2) + Cl) allowing for a much more rapid consumption of BCl(3) than the direct abstraction reaction (BCl(3) + H → BCl(2) + HCl) considered up to now. At low temperatures, T ≤ 800 K, our results show that a weakly stabilized complex BHCl(3) is formed with a kinetic law compatible with the consumption rate measured in the former experiments. At higher temperatures, this complex is not stable and then easily eliminates a chlorine atom. Our work also shows that a very similar mechanism, involving the same intermediate and sharing the same transition state, allows for the elimination of HCl. A dividing coefficient between these two elimination pathways is obtained from both a potential energy surface based statistical analysis and an ab initio molecular dynamics transition path sampling simulation. It finally allows partitioning of the global consumption rate of BCl(3) in terms of the formation of (i) BHCl(3), (ii) BHCl(2) + Cl through a H addition/Cl elimination mechanism, (iii) BCl(2) + HCl through a H addition/HCl elimination mechanism, and (iv) BCl(2) + HCl through direct abstraction.     

Theoretical approach to the mechanism of reactions between halogen atoms and unsaturated hydrocarbons: the Cl + propene reaction

The potential energy surface for the Cl + propene reaction was analyzed at the MP2 level using Pople's 6-31G(d,p) and 6-311+G(d,p), and Dunning's cc-pVDZ and aug-cc-pVDZ basis sets. Two different channels for the addition reaction leading to chloroalkyl radicals and five alternative channels for the abstraction reaction leading to C(3)H(5) (.) + HCl were explored. The corresponding energy profiles were computed at the QCISD(T)/aug-cc-pVDZ//MP2/aug-cc-pVDZ level of theory. Theoretical results suggest that the previously established mechanism consisting of (1) direct abstraction and (2) addition-elimination steps is instead made up of (1) addition through an intermediate and (2) two-step abstraction processes. No direct abstraction mechanism exists on the potential energy surface. The kinetic equations derived for the new mechanism are consistent with the pressure dependence experimentally observed for this reaction.     

Quasiclassical trajectory study of the Cl+CH4 reaction dynamics on a quadratic configuration interaction with single and double excitation interpolated potential energy surface

An ab initio interpolated potential energy surface (PES) for the Cl+CH(4) reactive system has been constructed using the interpolation method of Collins and co-workers [J. Chem. Phys. 102, 5647 (1995); 108, 8302 (1998); 111, 816 (1999); Theor. Chem. Acc. 108, 313 (2002)]. The ab initio calculations have been performed using quadratic configuration interaction with single and double excitation theory to build the PES. A simple scaling all correlation technique has been used to obtain a PES which yields a barrier height and reaction energy in good agreement with high level ab initio calculations and experimental measurements. Using these interpolated PESs, a detailed quasiclassical trajectory study of integral and differential cross sections, product rovibrational populations, and internal energy distributions has been carried out for the Cl+CH(4) and Cl+CD(4) reactions, and the theoretical results have been compared with the available experimental data. It has been shown that the calculated total reaction cross sections versus collision energy for the Cl+CH(4) and Cl+CD(4) reactions is very sensitive to the barrier height. Besides, due to the zero-point energy (ZPE) leakage of the CH(4) molecule to the reaction coordinate in the quasiclassical trajectory (QCT) calculations, the reaction threshold falls below the barrier height of the PES. The ZPE leakage leads to CH(3) and HCl coproducts with internal energy below its corresponding ZPEs. We have shown that a Gaussian binning (GB) analysis of the trajectories yields excitation functions in somehow better agreement with the experimental determinations. The HCl(v'=0) and DCl(v'=0) rotational distributions are as well very sensitive to the ZPE problem. The GB correction narrows and shifts the rotational distributions to lower values of the rotational quantum numbers. However, the present QCT rotational distributions are still hotter than the experimental distributions. In both reactions the angular distributions shift from backward peaked to sideways peaked as collision energy increases, as seen in the experiments and other theoretical calculations.