反应N+NH→N2+H的态-态量子动力学研究

采用在MRCI/aug-cc-pVQZ水平上构建的N₂H基态势能面, 并运用Chebyshev实波包法研究了N + NH→N₂ + H反应的量子动力学, 如反应几率、 积分截面以及产物振转态分布等. 在50~500 K温度范围内, 该反应的速率常数随着温度升高而递增, 与基于其它势能面的理论结果吻合. 然而, 在室温条件下, 所有理论计算的速率常数均显著大于实验值.     

Reactions of N+, N2+, and N3+ with NO from 300 to 1400 K

Rate constants have been measured from 300 to 1400 K in a selected ion flow tube (SIFT) and a high temperature flowing afterglow for the reactions of N+, N2+ and N3+ with NO. In all of the systems, the rate constants are substantially less than the collision rate constant. Comparing the high temperature results to kinetics studies as a function of translational energy show that all types of energy (translational, rotational, and vibrational) affect the reactivity approximately equally for all three ions. Branching ratios have also been measured at 300 and 500 K in a SIFT for the N+ and N3+ reactions. An increase in the N2+ product at the expense of NO+ nondissociative charge transfer product occurs at 500 K with N+. The branching ratios for the reaction of N3+ with NO have also been measured in the SIFT, showing that only nondissociative charge transfer giving NO+ occurs up to 500 K. The current results are discussed in the context of the many previous studies of these ions in the literature.     

Experimental and theoretical study of the reaction of the ethynyl radical with nitrous oxide, C2H + N2O

We investigated the rate constants and reaction mechanism of the gas phase reaction between the ethynyl radical and nitrous oxide (C(2)H + N(2)O) using both experimental methods and electronic structure calculations. A pulsed-laser photolysis/chemiluminescence technique was used to determine the absolute rate coefficient over the temperature range 570 K to 836 K. In this experimental temperature range, the measured temperature dependence of the overall rate constants can be expressed as: k(T) (C(2)H + N(2)O) = 2.93 × 10(-11) exp((-4000 ± 1100) K/T) cm(3) s(-1) (95% statistical confidence). Portions of the C(2)H + N(2)O potential energy surface (PES), containing low-energy pathways, were constructed using the composite G3B3 method. A multi-step reaction route leading to the products HCCO + N(2) is clearly preferred. The high selectivity between product channels favouring N(2) formation occurs very early. The pathway corresponds to the addition of the terminal C atom of C(2)H to the terminal N atom of N(2)O. Refined calculations using the coupled-cluster theory whose electronic energies were extrapolated to the complete basis set limit CCSD(T)/CBS led to an energy barrier of 6.0 kcal mol(-1) for the entrance channel. The overall rate constant was also determined by application of transition-state theory and Rice-Ramsperger-Kassel-Marcus (RRKM) statistical analyses to the PES. The computed rate constants have similar temperature dependence to the experimental values, though were somewhat lower.     

Theoretical analysis of the rate constants for the interstellar reaction N+OH→NO+H

The title reaction, a key elementary process involved in the chemistry of molecular clouds, has been theoretically studied over the 5–600 K temperature range. Rate constants calculations have been carried out using the full version of the statistical adiabatic channel model in conjunction with a potential energy surface that has been derived from recent ab initio quantum chemical data. By using various switching functions, the influence of the attenuation of the bound-complex bending frequency upon N? OH bond elongation on the temperature dependence of the reaction was investigated. The rate constants exhibit a slightly positive temperature dependence with a calculated rate constant value at 300 K in very good agreement with the measured value. A comparison with the available experimental data between 250 and 515 K suggests that recrossing trajectories might occur with increasing importance as the temperature increases. However, the nonstatistical recrossing effects are expected to be of minor importance at interstellar temperatures such that the rate constants over the 5–200 K temperature range are given by k = 8.41 × 10^?12 T^+0.30 cm³ molecule^?1 s^?1. The rate constant calculated at 10 K is consistent with that derived in the astrochemical modeling of the L134N dark cloud. Rate constants for individual quantum states are also presented. © 1995 John Wiley & Sons, Inc.     

Direct dynamics study on the hydrogen abstraction reactions N2H4 + R-->N2H3 + RH (R=NH2, CH3)

We present a direct ab initio dynamics study on the hydrogen abstraction reactions N(2)H(4)+R-->N(2)H(3)+RH (R=NH(2),CH(3)), which are predicted to have six possible reaction channels for NH(2) abstraction and four for CH(3) abstraction caused by the different N(2)H(4) isomers and various attacking orientations of foreign radical to N(2)H(4). The structures and frequencies at the stationary points and the points along the minimum energy paths (MEPs) of all reaction channels are obtained at the UMP2(full)6-31+G(d,p) level of theory. Energetic information of stationary points and the points along the MEPs is further refined by means of MC-QCISD method. The rate constants of these channels are calculated using the improved canonical variational transition-state theory with the small-curvature tunneling correction (ICVT/SCT) method. The calculated results show that the favorable reaction channels are channels (n1) and (n4) as well as (c1) and (c3) (refer to Scheme 1) in the whole temperature range. The total ICVT/SCT rate constants of all channels for the two reactions at the MC-QCISDUMP2(full)6-31+G(d,p) level are both in good agreement with the available experimental data, and corresponding three-parameter expressions of k(ICVTSCT) in 220-3000 K are fitted as 6.46 x 10(-15)(T298)(3.60) exp(-386T) cm(3) mol(-1) s(-1) for NH(2) abstraction and 1.04 x 10(-14)(T298)(4.00) exp(-2037T) cm(3) mol(-1) s(-1) for CH(3) abstraction. Additionally, the long range interaction between the H atom of X-H bond in foreign radicals and the lone pair on the nonreactive N atom of the transition states is further discussed to explain the various transition-state numbers of the two similar hydrogen abstraction reactions.     

Quantum Mechanics Rate Constant for the N+ND Reaction

We present nonadiabatic quantum dynamical calculations on the two coupled potential en-ergy surfaces (1²A′ and 2²A′) [J. Theor. Comput. Chem. 8, 849 (2009)] for the reaction. Initial state-resolved reaction probabilities and cross sections for the N+ND→N₂+D reaction and N′+ND→N′D+N reaction for collision energies of 5 meV to 1.0 eV are determined, re-spectively. It is found that the N+ND→N₂+D reaction is dominated in the N+ND reaction.In addition, we obtained the rate constants for the N+ND→N₂+D reaction which demand further experimental investigations.     

On the chaperon mechanism: application to ClO + ClO (+N2) --> ClOOCl (+N2)

The dynamics of the ClO + ClO (+N(2)) radical complex (or chaperon) mechanism is studied by electronic structure methods and quasi-classical trajectory calculations. The geometries and frequencies of the stationary points on the potential energy surface (PES) are optimized at the B3LYP/6-311+G(3df) level of theory, and the energies are refined at the CCSD(T)/6-311+G(3df) (single-point) level of theory. Basis set superposition error (BSSE) corrections are applied to obtain 1.5 kcal mol(-1) for the binding energy of the ClO.N(2) van der Waals (VDW) complex. A model PES is developed and used in quasi-classical trajectory calculations to obtain the capture rate constant and nascent energy distributions of ClOOCl* produced via the chaperon mechanism. A range of VDW binding energies from 1.5 to 9.0 kcal mol(-1) are investigated. The anisotropic PES for the ClO.N(2) complex and a separable anharmonic oscillator approximation are used to estimate the equilibrium constant for formation of the VDW complex. Rate constants, branching ratios to produce ClOOCl, and nascent energy distributions of excited ClOOCl* are discussed with respect to the VDW binding energy and temperature. Interestingly, even for weak VDW binding energies, the N(2) usually carries away enough energy to stabilize the nascent ClOOCl*, although the VDW equilibrium constant is small. For stronger binding energies, the stabilization efficiency is reduced, but the capture rate constant is increased commensurately. The resulting rate constants for forming ClOOCl* from the title reaction are only weakly dependent on the VDW binding energy and temperature. As a result, the relative importance of the chaperon mechanism is mostly dependent on the VDW equilibrium constant. For the calculated ClO.N(2) binding energy of 1.5 kcal mol(-1), the VDW equilibrium constant is small, and the chaperon mechanism is only important at very high pressures.     

Mechanisms of electron-ion recombination of N2H+/N2D+ and HCO+/DCO+ ions: temperature dependence and isotopic effect

The temperature dependencies of the rate coefficients, alpha(e), for electron-ion dissociative recombination (DR) of N2H+/N2D+ and HCO+/DCO+ ions with electrons have been measured over the range 100-500 K. Also, optical emissions have been detected at approximately 100 K from the N2(B3(pi)g) electronically excited products of N2H+/N2D+ recombination. The measurements were carried out using the classic FALP technique combined with an optical monochromator. For N2H+, there was no variation of alpha(e) with temperature above 200 K, with an average value of alpha(e)(N2H+) = 2.8 x 10(-7) cm3 s(-1). The temperature variation for T approximately 100-300 K observed for alpha(e)(HCO+) is similar to that of N2H+ ions for T approximately 300-500 K. The smaller rate coefficient measured for DCO+ and N2D+ ions shows the influence of an isotope effect. The substantial enhancement of the vibrational level, upsilon' = 6, from the N2B state for N2H+ recombination over N2D+ recombination is consistent with previous result at 300 K and implies the influence of a tunneling mechanism of DR.     

Thermal rate constant calculation of the NF + F reactive system multiple arrangements

In this work we analyzed the multiple channels of the reaction NF+F through the evaluation of thermal rate constants with both Wigner and Eckart tunneling corrections. Minimum energy paths and intrinsic reaction coordinates of the systems were obtained and accurately studied in order to ensure the consistency of our results. Specifically, we investigated the NF + F = N + F(2), NF + F = NF + F, and NF(2) = NF + F, reactive systems. As experimental data are available for the latter reaction, we were able to conclude that our thermal rate constants are in agreement for a wide range of temperatures. The here performed study is relevant to the understanding of the decomposition process of nitrogen trifluoride (NF(3)).     

Ab initio chemical kinetics for the OH+HNCN reaction

The kinetics and mechanism of the reaction of the cyanomidyl radical (HNCN) with the hydroxyl radical (OH) have been investigated by ab initio calculations with rate constants prediction. The single and triplet potential energy surfaces of this reaction have been calculated by single-point calculations at the CCSD(T)/6-311+G(3df,2p) level based on geometries optimized at the B3LYP/6-311+G(3df,2p) and CCSD/6-311++G(d,p) levels. The rate constants for various product channels in the temperature range of 300-3000 K are predicted by variational transition-state and Rice-Ramsperger-Kassel-Marcus (RRKM) theories. The predicted total rate constants can be represented by the expressions ktotal=2.66 x 10(+2)xT-4.50 exp(-239/T) in which T=300-1000 K and 1.38x10(-20)xT2.78 exp(1578/T) cm3 molecule(-1) s(-1) where T=1000-3000 K. The branching ratios of primary channels are predicted: k1 for forming singlet HON(H)CN accounts for 0.32-0.28, and k4 for forming singlet HONCNH accounts for 0.68-0.17 in the temperature range of 300-800 K. k2+k7 for producing H2O+NCN accounts for 0.55-0.99 in the high-temperature range of 800-3000 K. The branching ratios of k3 for producing HCN+HNO, k6 for producing H2N+NCO, k8 for forming 3HN(OH)CN, k9 for producing CNOH+3NH, and k5+k10 for producing NH2+NCO are negligible. The rate constants for key individual product channels are provided in a table for different temperature and pressure conditions.     

Vibrational energy transfer in N2-N2 collisions: a new semiclassical study

The vibrational energy relaxation in collisions between N2 molecules in the low- and medium-lying vibrationally excited levels was revisited using the semiclassical coupled-state method and the use of two different potential-energy surfaces having the same short-range potential recently determined from ab initio calculations but with different long-range interactions. Compared to the data reported in the classical work by Billing and Fisher [Chem. Phys. 43, 395 (1979)], the newly calculated vibration-to-translation rate constant K(1,0 / 0,0) is in much better agreement with the available experimental data over a large temperature interval, from T = 200 K up to T = 6000 K. Nevertheless, as far as the vibration-to-translation exchanges are concerned, the lower-temperature regime remains quite critical in that the new rate constants do not completely account for the rate constant curvature suggested by the experiments for temperatures lower than T = 500 K. The dependence of the state-selected vibration-to-vibration rate constants, K(v,v-delta v / 0,1), both upon the vibrational quantum number v and the gas temperature are calculated. The substantial deviations from previously found behaviors could have major consequences for the vibrational kinetic modeling of N2-containing gas mixtures.     

A dual-level ab initio and hybrid density functional theory dynamics study on the unimolecular decomposition reaction C2H5O-->CH2O + CH3

We present a direct ab initio and hybrid density functional theory dynamics study of the thermal rate constants of the unimolecular decomposition reaction of C2H5O-->CH2O + CH3 at a high-pressure limit. MPW1K/6-31+G(d,p), MP2/6-31+G(d,p), and MP2(full)/6-31G(d) methods were employed to optimize the geometries of all stationary points and to calculate the minimum energy path (MEP). The energies of all the stationary points were refined at a series of multicoefficient and multilevel methods. Among all methods, the QCISD(T)/aug-cc-pVTZ energies are in good agreement with the available experimental data. The rate constants were evaluated based on the energetics from the QCISD(T)/aug-cc-pVTZ//MPW1K/6-31+G(d,p) level of theory using both microcanonical variational transition state theory (microVT) and RRKM theory with the Eckart tunneling correction in the temperature range of 300-2500 K. The calculated rate constants at the QCISD(T)/aug-cc-pVTZ/MPW1K/6-31+G(d,p) level of theory are in good consistent with experimental data. The fitted three-parameter Arrhenius expression from the microVT/Eckart rate constants in the temperature range 200-2500 K is k = 2.52 x 10(12)T(0.41)e(-8894.0/T) s(-1). The falloff curves of pressure-dependent rate constants are performed using master-equation method within the temperature range of 391-471 K. The calculated results are in good agreement with the available experimental data.     

Ab initio studies of ClO(x) reactions: prediction of the rate constants of ClO+NO2 for the forward and reverse processes.

The potential-energy surface for the reaction of ClO with NO2 has been constructed at the CCSD(T)/6-311+G(3df)//B3LYP/6-311+G(3df) level of theory. Six ClNO3 isomers are located; these are ClONO2, pc-ClOONO, pt-ClOONO, OClNO2, pt-OClONO, pc-OClONO, with predicted energies relative to the reactants of -25.6, -0.5, 1.0, 1.9, 12.2 and 13.6 kcal mol-1, and heats of formation at 0 K of 7.8, 32.9, 34.4, 35.5, 45.6 and 47.0 kcal mol-1, respectively. Isomerizations among them are also discussed. The rate constants for the low-energy pathways have been computed by statistical theory calculations. For the association reaction producing exclusively ClONO2, the predicted low- and high-pressure-limit rate constants in N2 for the temperature range of 200-600 K can be represented by: (N2)=3.19 x 10-17 T-5.54 exp(-384 K/T) cm6 molecule-2 s-1 and =3.33 x 10-7 T-1.48 exp(-18 K/T) cm3 molecule-1 s-1. The predicted low- and high-pressure-limit rate constants for the decomposition of ClONO2 in N2 at 200-600 K can be expressed, respectively, by =6.08 x 1013 T-6.54 exp(-13813 K/T) cm3 molecule-1 s-1 and =4.59 x 1023 T-2.43 exp(-13437 K/T) s-1. The predicted values compare satisfactorily with available experimental data. The reverse Cl+NO3 reaction was found to be independent of the pressure, giving exclusively ClO+NO2; the predicted rate constant can be expressed as k(Cl+NO3)=1.19 x 10-9 T-0.60 exp(58 K/T) cm3 molecule-1 s-1..     

Seven dimensional quantum dynamics study of the H2+NH2-->H+NH3 reaction

Initial state-selected time-dependent wave packet dynamics calculations have been performed for the H2+NH2-->H+NH3 reaction using a seven dimensional model on an analytical potential energy surface based on the one developed by Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The model assumes that the two spectator NH bonds are fixed at their equilibrium values and nonreactive NH2 group keeps C2v symmetry and the rotation-vibration coupling in NH2 is neglected. The total reaction probabilities are calculated when the two reactants are initially at their ground states, when the NH2 bending mode is excited, and when H2 is on its first vibrational excited state, with total angular momentum J=0. The converged cross sections for the reaction are also reported for these initial states. Thermal rate constants and equilibrium constants are calculated for the temperature range of 200-2000 K and compared with transition state theory results and the available experimental data. The study shows that (a) the reaction is dominated by ground-state reactivity and the main contribution to the thermal rate constants is thought to come from this state, (b) the excitation energy of H2 was used to enhance reactivity while the excitation of the NH2 bending mode hampers the reaction, (c) the calculated thermal rate constants are very close to the experimental data and transition state theory results at high and middle temperature, while they are ten times higher than that of transition state theory at low temperature (T=200 K), and (d) the equilibrium constants results indicate that the approximations applied may have different roles in the forward and reverse reactions.     

Potential energy surface for asymmetrically substituted reactions of type CWXYZ + A. Kinetics study

The CHClF(2) + Cl --> ClH + CClF(2) gas-phase abstraction reaction was chosen as a model of asymmetrically substituted polyatomic reactions of type CWXYZ + A --> products. The analytical potential energy surface for this reaction was constructed with suitable functional forms to represent vibrational modes, and calibrated with respect to experimental thermal rate constants which are available over the temperature range 296-410 K. On this surface, the thermal rate constants were calculated using variational transition-state theory with semiclassical transmission coefficients over a wider temperature range, 200-2500 K, therefore obtaining kinetics information at lower and higher temperatures than are experimentally available. This surface was also used to analyze dynamical features, such as tunneling and reaction-path curvature. In the former, the influence of the tunneling factor was important since a light hydrogen atom passes through the barrier. In the latter, it was found that vibrational excitation of the C-F and C-Cl stretching modes can be expected in the exit channel.     

Quantum statistical study of O + O2 isotopic exchange reactions: cross sections and rate constants

Using a wave packet based statistical model, we compute cross sections and thermal rate constants for various isotopic variants of the O + O2 exchange reaction on a recently modified ab initio potential energy surface. The calculation predicts a highly excited rotational distribution and relatively cold vibrational distribution for the diatomic product. A small but important threshold effect was identified for the (16)O + 18O2 reaction, which is suggested to contribute to the experimentally observed negative temperature dependence of the rate ratio, k(18O + 16O2)/k(16O + 18O2). Despite reasonable agreement with quasiclassical trajectory results, however, the calculated thermal rate constants are smaller than experimental measurements by a factor from 2 to 5. The experimentally observed negative temperature dependence of the rate constants is not reproduced. Possible reasons for the theory-experiment discrepancies are discussed.     

Mechanistic and kinetic investigations of N2H4 + OH reaction

The reaction of N₂H₄ with OH has been investigated by quantum chemical methods. The results show that hydrogen abstraction mechanism is more feasible than substitution mechanism thermodynamically. The calculated rate constants agree with the available experimental data. The calculated results show that the variational effect is small at lower temperature region, while it becomes significant at higher temperature region. On the other hand, the small‐curvature tunneling effect may play an important role in the temperature range 220?3000 K. Moreover, the calculated rate constants show negative temperature dependence at the temperatures below 500 K, which is in accordance with Vaghjiani's report that slightly negative temperature dependence is found over the temperature range of 258?637 K. The mechanism of the major product (N₂H₃) with OH has also been investigated theoretically to understand the title reaction thoroughly. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Theoretical study on the OH + CH3NHCOOCH3 reaction

The multiple-channel reactions OH + CH3NHC(O)OCH3 --> products are investigated by direct dynamics method. The optimized geometries, frequencies, and minimum energy path are all obtained at the MP2/6-311+G(d,p) level, and energetic information is further refined by the BMC-CCSD (single-point) method. The rate constants for every reaction channels, R1, R2, R3, and R4, are calculated by canonical variational transition state theory with small-curvature tunneling correction over the temperature range 200-1000 K. The total rate constants are in good agreement with the available experimental data and the two-parameter expression k(T) = 3.95 x 10(-12) exp(15.41/T) cm3 molecule(-1) s(-1) over the temperature range 200-1000 K is given. Our calculations indicate that hydrogen abstraction channels R1 and R2 are the major channels due to the smaller barrier height among four channels considered, and the other two channels to yield CH3NC(O)OCH3 + H2O and CH3NHC(O)(OH)OCH3 + H2O are minor channels over the whole temperature range.     

Direct dynamics study of the hydrogen abstraction reaction of CF3CH2Cl + Cl → CF3CHCl + HCl

The hydrogen abstraction reaction of Cl atoms with CF₃CH₂Cl (HCFC‐133a) is investigated by using density function theory and ab initio approach, and the rate constants are calculated by using the dual‐level direct dynamics method. Optimized geometries and frequencies of reactants, transition state, and products are computed at the B3LYP/6‐311+G(2d,2p) level. To refine the energetic information along the minimum energy path, single‐point energy calculations are carried out at the G3(MP2) level of theory. The interpolated single‐point energy method is employed to correct the energy profiles for the title reaction. The rate constants are evaluated by using the canonical variational transition state theory with a small‐curvature tunneling correction over a wide range of temperature, 200–2000 K. The variational effect for the reaction is moderate at low temperatures and very small at high temperatures. However, the tunneling correction has an important contribution in the lower temperature range. The agreement between calculated rate constants and available experimental values is good at lower temperatures but diverges significantly at higher temperatures. © 2012 Wiley Periodicals, Inc. Int J Chem Kinet 44: 661–667, 2012     

Dual-level direct dynamics studies on the reaction Cl + CHBr(2)Cl

Theoretical investigations are carried out on the multichannel reaction CHBr(2)Cl + Cl by means of direct dynamics methods. The minimum energy path (MEP) is obtained at the BH&H-LYP/6-311G(d,p) level, and energetic information is further refined at the CCSD(T)/6-311+G(2df,2p) (single-point) level. The rate constants for three reaction channels, H-abstraction, Br-abstraction, and Cl-abstraction, are calculated by using the improved canonical variational transition state theory (ICVT) incorporating with the small-curvature tunneling (SCT) correction. The theoretical overall rate constants are in good agreement with the available experimental data and are found to be k=2.58 x 10(-15) T(1.18) exp(-861.17/T) cm(3)molecule(-1)s(-1) over the temperature range 200--2400 K. For the title reaction, H-abstraction reaction channel is the major channel at the lower temperatures, while as the temperature increases, the contribution of Br-abstraction reaction channel should be taken into account. At 2180 K, the rate constants of these two pathways are equal. Cl-abstraction reaction channel is minor channel over the whole temperature region.