Reaction pathway for the nonadiabatic reaction of Ca(4s3d 1D)+H2-->CaH(X 2Sigma+)+H

The reaction pathway and the nascent CaH product distribution in the reaction Ca(4s3d (1)D)+H(2)-->CaH(X (2)Sigma(+))+H are obtained using a pump-probe technique. The Ca atom is first prepared in the 3 (1)D state by a two-photon absorption, and then in brief time delay the laser-induced fluorescence of the reaction product CaH is monitored. The CaH(v=0,1) distributions appear to be single peaked, as characterized by Boltzmann rotational temperature of 807+/-38 K (v=0) and 684+/-77 K (v=1). The vibrational population ratio of CaH(v=0)/CaH(v=1) is determined to be 3.3+/-0.1, while the v=2 population is not detectable. The fractions of the available energy partitioning into rotation, vibration, and translation are estimated to be 0.36+/-0.05, 0.28+/-0.04, and 0.36+/-0.05, respectively. With the aid of the potential energy surfaces calculations, the current reaction should favor a near C(2v) collision configuration. The temperature dependence measurement yields a positive slope, indicative of the reaction occurrence without any potential barrier. The colliding species are anticipated to follow an attractive 1B(2) (or 2A') surface and then transit nonadiabatically to the reactive ground state surface.     

Reaction pathway and potential barrier for the CaH product in the reaction of Ca(4s4p1P1) + H2 --> CaH(X2Sigma+) + H

The nascent CaH product in the reaction Ca(4s4p1P1) + H2 --> CaH(X2Sigma+) + H is obtained using a pump-probe technique. The CaH(v = 0,1) distributions, with a population ratio of CaH(v = 0)/CaH(v = 1) = 2.7+/-0.2, may be characterized by low Boltzmann rotational temperature. According to Arrhenius theory, the temperature dependence measurement yields a potential barrier of 3820+/-480 cm(-1) for the current reaction. As a result of the potential energy surfaces (PES) calculations, the reaction pathway favors a Ca insertion into the H2 bond along a (near) C2v geometric approach. As the H2 bond is elongated, the configurational mixing between the orbital components of the 4p and nearby low-lying 3d state with the same symmetry makes significant the nonadiabatic transition between the 5A' and 2A' surface in the repulsive limbs. Therefore, the collision species are anticipated to track along the 5A' surface, then undergo nonadiabatic transition to the inner limb of the 2A' surface, and finally cross to the reactive 1A' surface. The observed energy barrier probably accounts for the energy requirement to surmount the repulsive hill in the entrance. The findings of the nascent CaH product distributions may be reasonably interpreted from the nature of the intermediate structure and lifetime after the 2A'-1A' surface transition. The distinct product distributions between the Ca(4 1P1) and Mg(3 1P1) reactions with H2 may also be realized with the aid of the PES calculations.     

Reaction dynamics of Si(3PJ)+O2-->SiO(X 1 Sigma+)+O studied by a crossed-beam laser-induced fluorescence technique

Oxidation reaction of the ground state Si atom was studied by using a crossed molecular beam technique at 13.0 kJ/mol of collision energy. The Si atomic beam was generated by laser vaporization and crossed with the oxygen molecular beam at right angle. Products at the crossing region were detected by the laser-induced fluorescence (LIF). The LIF of SiO(A 1 Pi-X 1 Sigma+) was used to determine the vibrational state distribution of the electronic ground state, SiO(X 1 Sigma+). The determined distribution was inverted with the maximum population at v"=4, and in good agreement with the recent quasiclassical trajectory calculation on the singlet potential energy surface. The agreement suggested that an abstraction mechanism is dominant at the collision energy studied here. One of the counterproducts, O(3PJ), was also observed by the vacuum ultraviolet LIF and the distribution of the spin-orbit levels were determined. The formation of O(3PJ) was consistent with the significant population of v"=7 and 8 states of SiO, which could be explained by the presence of the triplet product channel with higher exothermicity.     

Differential cross sections and product energy distributions for the C(3P)+OH(X 2Pi)-->CO(X 1Sigma+)+H(2S) reaction using a quasiclassical trajectory method

Quasiclassical trajectory calculations have been carried out for the C((3)P)+OH(X (2)Pi)-->CO(X (1)Sigma(+))+H((2)S) reaction using a recent ab initio potential energy surface for the ground electronic state X (2)A(') of COH. Differential cross sections (DCSs), and product vibrational, rotational and translational distributions have been determined for a wide range of collision energies (0.001-1 eV). The role of excitations (rotation or vibration) of the OH reactant on these quantities has been investigated. Product vibrational, rotational, and translational distributions are found to be almost independent on the rovibrational state of OH, whereas DCSs show a weak dependence on the initial rotational state of OH. We also analyze the results using a study based on the lifetime of the intermediate complex and on the kinematic constraint associated with the mass combination.     

Spin-orbit effect in the energy pooling reaction O2(a 1Delta)+O2(a 1Delta)-->O2(b 1Sigma)+O2(X 3Sigma)

Five-dimensional nonadiabatic quantum dynamics studies have been carried out on two new potential energy surfaces of S(2)((1)A(')) and T(7)((3)A(")) states for the title oxygen molecules collision with coplanar configurations, along with the spin-orbit coupling between them. The ab initio calculations are based on complete active state second-order perturbation theory with the 6-31+G(d) basis set. The calculated spin-orbit induced transition probability as a function of collision energy is found to be very small for this energy pooling reaction. The rate constant obtained from a uniform J-shifting approach is compared with the existing theoretical and experimental data, and the spin-orbit effect is also discussed in this electronic energy-transfer process.     

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

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

Potential surfaces and dynamics of the O(3P)+H2O(X1A1)-->OH(X2pi)+OH(X2pi) reaction

We present global potential energy surfaces for the three lowest triplet states in O(3P)+H2O(X1A1) collisions and present results of classical dynamics calculations on the O(3P)+H2O(X1A1)-->OH(X2pi)+OH(X2pi) reaction using these surfaces. The surfaces are spline-based fits of approximately 20,000 fixed geometry ab initio calculations at the complete-active-space self-consistent field+second-order perturbation theory (CASSCF+MP2) level with a O(4s3p2d1f)/H(3s2p) one electron basis set. Computed rate constants compare well to measurements in the 1000-2500 K range using these surfaces. We also compute the total, rovibrationally resolved, and differential angular cross sections at fixed collision velocities from near threshold at approximately 4 km s(-1) (16.9 kcal mol(-1) collision energy) to 11 km s(-1) (122.5 kcal mol(-1) collision energy), and we compare these computed cross sections to available space-based and laboratory data. A major finding of the present work is that above approximately 40 kcal mol(-1) collision energy rovibrationally excited OH(X2pi) products are a significant and perhaps dominant contributor to the observed 1-5 micro spectral emission from O(3P)+H2O(X1A1) collisions. Another important result is that OH(X2pi) products are formed in two distinct rovibrational distributions. The "active" OH products are formed with the reagent O atom, and their rovibrational distributions are extremely hot. The remaining "spectator" OH is relatively rovibrationally cold. For the active OH, rotational energy is dominant at all collision velocities, but the opposite holds for the spectator OH. Summed over both OH products, below approximately 50 kcal mol(-1) collision energy, vibration dominates the OH internal energy, and above approximately 50 kcal mol(-1) rotation is greater than vibrational energy. As the collision energy increases, energy is diverted from vibration to mostly translational energy. We note that the present fitted surfaces can also be used to investigate direct collisional excitation of H2O(X1A1) by O(3P) and also OH(X2pi)+OH(X2pi) collisions.     

Full six-dimensional nonadiabatic quantum dynamics calculation for the energy pooling reaction O(2)(a (1)Delta)+O(2)(a (1)Delta)-->O(2)(b (1)Sigma)+O(2)(X (3)Sigma)

Six new potential energy surfaces of four singlet states and two triplet states for the title oxygen molecule reaction along with the spin-orbit coupling among them have been constructed from the complete active space second-order perturbation theory with a 6-311+G(d) basis. Accurate integral cross sections are calculated with a full six-dimensional nonadiabatic time-dependent quantum wave packet method. The thermal rate constant based on the integral cross sections agrees well with the result of the experimental measurements, and the intersystem crossing effects are also discussed in this electronic energy-transfer process.     

Quantum mechanical investigations of the N(4S)+O2(X 3Sigmag -)-->NO(X 2Pi)+O(3P) reaction

The reaction between energetic nitrogen atoms and oxygen molecules has received important attention in connection with nitric oxide chemistry in the lower thermosphere. We report time-independent quantum mechanical calculations of the N(4S)+O2-->NO+O reaction employing the X 2A' and a 4A' electronic potential energy surfaces of Sayos et al. [J. Chem. Phys. 117, 670 (2002)]. We confirm the production of highly vibrationally excited NO molecules, consistent with previous semiclassical and more recent time-dependent quantum wave packet studies. Calculations are carried out for total angular momentum quantum number J=0 and cross sections and rate coefficients are extracted using the J-shifting approximation. The results are in good agreement with available experimental and theoretical data.     

Quantum real wave-packet dynamics of the N(4S)+NO(X2Pi)-->N2(X1Sigma(g)+)+O(3P) reaction on the ground and first excited triplet potential energy surfaces: rate constants, cross sections, and product distributions

The reaction N+NO-->N(2)+O was studied by means of the time-dependent real wave-packet (WP) method and the J-shifting approximation. We consider the ground 1 (3)A(") and first excited 1 (3)A(') triplet states, which correlate with both reactants and products, using analytical potential energy surfaces (PESs) recently developed in our group. This work extends our previous quantum dynamics study, and probabilities, cross sections, and rate constants were calculated and interpreted on the basis of the different shapes of the PESs (barrierless 1 (3)A(") and with barrier 1 (3)A(') surfaces, respectively). The WP rate constant (k(1)) shows a weak dependence on T(200-2500 K), as the dominant contribution to reactivity is provided by the barrierless ground PES. There is a good agreement of WP k(1) with the measurements and variational transition state theory (VTST) data, and also between the WP and VTST k(1)(1 (3)A(")) results. Nevertheless, there is a large discrepancy between the WP and VTST k(1)(1 (3)A(')) results. Product state distributions were also calculated for the much more reactive 1 (3)A(") PES. There is an excellent agreement with the experimental average fraction of vibrational energy in N(2)(25+/-3%), the only measured dynamics property of this reaction.     

Evidence for the production of N4 via the N2 A3Sigma(u)+ + N2 A3Sigma(u)+ energy pooling reaction

We report mass spectrometric evidence supporting our proposed mechanistic pathway for the production of N4 through the energy pooling reaction N2 A3Sigma(u)+ + N2 A3Sigma(u)+. N2 A3Sigma(u)+ is generated from the quenching of resonantly excited xenon in a mixture of xenon, 15N2, and 14N2 that is illuminated with xenon resonant lamps (147 nm). Mass spectra are periodically taken of the mixture. Over time, we observe significant isotopic scrambling of the 15N2 and 14N2, generating 15N14N in concentrations approaching 10% (approximately 2 Torr) of the initial 15N2 concentration. Though we do not observe the direct formation of N4, the isotopic ratios indicate that an excited complex (15N2(14)N2) exists long enough so that scrambling of the nitrogen atoms can occur, offering a possible route to the formation of tetrahedral nitrogen (1Td N4).     

Crossed-beam radical-radical reaction dynamics of O(3P)+C3H3-->H(2S)+C3H2O

The radical-radical oxidation reaction, O(3P)+C3H3 (propargyl)-->H(2S)+C3H2O (propynal), was investigated using vacuum-ultraviolet laser-induced fluorescence spectroscopy in a crossed-beam configuration, together with ab initio and statistical calculations. The barrierless addition of O(3P) to C3H3 is calculated to form energy-rich addition complexes on the lowest doublet potential energy surface, which subsequently undergo direct decomposition steps leading to the major reaction products, H+C3H(2)O (propynal). According to the nascent H-atom Doppler-profile analysis, the average translational energy of the products and the fraction of the average transitional energy to the total available energy were determined to be 5.09+/-0.36 kcal/mol and 0.077, respectively. On the basis of a comparison with statistical prior calculations, the reaction mechanism and the significant internal excitation of the polyatomic propynal product can be rationalized in terms of the formation of highly activated, short-lived addition-complex intermediates and the adiabaticity of the excess available energy along the reaction coordinate.     

Imaging the pair-correlated excitation function: The F+CH4-->HF(v')+CH3(nu=0) reaction

The velocity map ion imaging technique was applied to measure the reaction excitation function for the first time. It was found that the "raw" excitation function was significantly distorted by the density-to-flux transformation of the title reaction. Through a systematic investigation, possible reasons for such a dramatic effect are outlined. In addition, the state-resolved, pair-correlated excitation functions and branching ratios are presented. Effects of imperfect time slicing in the time-sliced velocity imaging technique in general are also discussed.     

Seven-dimensional quantum dynamics study of the O(3P)+CH4 reaction

The initial state selected time-dependent wave packet calculations have been carried out to study the title reaction with seven degrees of freedom included by restricting the nonreacting CH(3) group under C(3V) symmetry and the CH bond length in the group. Total reaction probabilities as well as integral cross sections were calculated for the ground and four vibrationally excited reagent states. Our calculation shows that the reactivity is very small for the reaction for collision energy up to 1.0 eV for all the initial states. Initial vibration excitation of CH(4), in particular, the CH stretch excitation, enhances the reactivity, but only part of the excitation energy deposited can be used to reduce the reaction threshold. The rate constant for the ground initial state agrees rather well with that from a recent quasiclassical trajectory study and is larger than that from the semirigid vibrating rotor target calculations, in particular, in the low temperature region. On the other hand, the thermal rate constant calculated from the integral cross sections for these five vibrational states is about a factor of 20 smaller than that from the multiconfiguration time-dependent Hartree study.     

CH2((X)3B1)自由基与N2O反应的研究

The reaction dynamics of methylene radicalCH2((x)3B1)with N2O was investigated by Time Resolved Fourier Transform Infrared Spectroscopy(TR-FTIRS). Pure CH2((x)3B1) radical was produced via laser photolysis of ketene at 351 nm.. Nascent vibrationally excited products CO , NO and HCN were observed.. Some reaction pathways which may lead to these products were proposed and a possible reaction mechanism was outlined..     

The effect of spin-orbit coupling on fast neutral chemical reaction O(3P)+CH3-->CH3O

The effect of nonadiabatic transitions through the spin-orbit couplings has been investigated on the fast neutral reaction, O((3)P)+CH(3)-->CH(3)O. Adiabatic potential energies and the spin-orbit coupling terms have been evaluated for the four electronic states of CH(3)O ((2)E, (2)A(2), (4)E, and (4)A(2)) that correlate with the O((3)P)+CH(3) asymptote, as a function of CO distance and OCH angle under the C(3v) symmetry, by ab initio electronic structure calculations using multireference internally contracted single and double excitation configuration interaction method with the 6-311G(2df,2pd) basis sets. Multistate quantum reactive scattering calculations have been carried out with the use of thus obtained potential energies and spin-orbit coupling matrices, based on the generalized R-matrix propagation method. The calculated thermal rate constants show a slight positive dependence on temperature in a range between 50 and 2000 K, supporting the previous experimental results. It is shown that the spin-orbit coupled excited states give rise to reflections over the centrifugal barrier due to the quantum interference. Classical capture calculations yield larger rate constants due to the neglect of quantum reflections. It is concluded that the effect of nonadiabatic transitions is of minor importance on the overall reactivity in this reaction.     

Potential energy surface, kinetics, and dynamics study of the Cl+CH4-->HCl+CH3 reaction

A modified and recalibrated potential energy surface for the gas-phase Cl+CH4-->HCl+CH3 reaction is reported and tested. It is completely symmetric with respect to the permutation of the four methane hydrogen atoms and is calibrated with respect to updated experimental and theoretical stationary point properties and experimental forward thermal rate constants. From the kinetics point of view, the forward and reverse thermal rate constants and the activation energies were calculated using the variational transition-state theory with semiclassical transmission coefficients over a wide temperature range of 150-2500 K. The theoretical results reproduce the available experimental data, with a small curvature of the Arrhenius plot which indicates the role of tunneling in this hydrogen abstraction reaction. A dynamics study was also performed on this PES using quasiclassical trajectory (QCT) calculations, including corrections to avoid zero-point energy leakage along the trajectories. First, we found a noticeable internal energy in the coproduct methyl radical, both in the ground-state [CH4 (v=0)] and vibrationally excited [CH4 (v=1)] reactions. This CH3 internal energy was directly precluded in some experiments or oversimplified in previous theoretical studies using pseudotriatomic models. Second, our QCT calculations give HCl rotational distributions slightly hotter than those in experiment, but correctly describing the experimental trend of decreasing the HCl product rotation excitation in going from HCl (v'=0) to HCl (v'=1) for the CH4 (v=1) reaction. Third, the state specific scattering distributions present a reasonable agreement with experiment, although they tend to make the reaction more forward and backward scattered than found experimentally probably because of the hotter rotational distribution and the deficiencies of the QCT methods.