首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Global analytic potential energy surfaces for O((3)P) + H(2)O((1)A(1)) collisions, including the OH + OH hydrogen abstraction and H + OOH hydrogen elimination channels, are presented. Ab initio electronic structure calculations were performed at the CASSCF + MP2 level with an O(4s3p2d1f)/H(3s2p) one electron basis set. Approximately 10(5) geometries were used to fit the three lowest triplet adiabatic states corresponding to the triply degenerate O((3)P) + H(2)O((1)A(1)) reactants. Transition state theory rate constant and total cross section calculations using classical trajectories to collision energies up to 120?kcal mol(-1) (~11?km s(-1) collision velocity) were performed and show good agreement with experimental data. Flux-velocity contour maps are presented at selected energies for H(2)O collisional excitation, OH + OH, and H + OOH channels to further investigate the dynamics, especially the competition and distinct dynamics of the two reactive channels. There are large differences in the contributions of each of the triplet surfaces to the reactive channels, especially at higher energies. The present surfaces should support quantitative modeling of O((3)P) + H(2)O((1)A(1)) collision processes up to ~150?kcal mol(-1).  相似文献   

2.
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.  相似文献   

3.
The mechanism of the reactions of W and W(+) with the water molecule have been studied for several lower-lying electronic states of tungsten centers at the CCSD(T)/6-311G(d,p)+SDD and B3LYP/6-31G(d,p)+SDD levels of theory. It is shown that these reactions are essentially multistate processes, during which lower-lying electronic states of the systems cross several times. They start with the formation of initial prereaction M(H(2)O) complexes with M-H(2)O bonding energies of 9.6 and 48.2 kcal/mol for M = W and W(+), followed by insertion of the metal center into an O-H bond with 20.0 and 53.3 kcal/mol barriers for neutral and cationic systems, respectively. The overall process of M + H(2)O --> t-HM(OH) is calculated to be highly exothermic, 48.4 and 48.8 kcal/mol for M = W and W(+). From the HM(OH) intermediate the reaction may proceed via several different channels, among which the stepwise HM(OH) --> HMO + H --> (H)(2)MO and concerted HM(OH) --> (H)(2)MO pathways are more favorable and can compete (energetically) with each other. For the neutral system (M = W), the concerted process is the most favorable, whereas for the charged system (M = W(+)), the stepwise pathway is slightly more favorable. From the energetically most favorable intermediate (H)(2)MO the reactions proceed via H(2)-molecule formation with a 53.1 kcal/mol activation barrier for the neutral system. For the cationic system, H-H formation and dissociation is an almost barrierless process. The overall reaction of W and W(+) with the water molecule leading to H(2) + MO formation is found to be exothermic by 48.2 and 39.8 kcal/mol, respectively. In the gas phase with the collision-less conditions the reactions W((7)S) + H(2)O --> H(2) + WO((3)Sigma(+)), and W(+)((6)D) + H(2)O --> H(2) + WO(+)((4)Sigma(+)) are expected to proceed via a 10.4 and 5.1 kcal/mol overall energy barrier corresponding to the first O-H dissociation at the TS1. On the basis of these PESs, we predict kinetic rate constants for the reactions of W and W(+) with H(2)O.  相似文献   

4.
The dynamics of the O(1D) + HCl --> OH + Cl(2P) reaction are investigated by a crossed molecular beam ion-imaging method and quasiclassical trajectory calculations on the three ab initio potential energy surfaces, the ground 1(1)A' and two excited (1(1)A' and 2(1)A') states. The scattering experiment was carried out at collision energies of 4.2, 4.5, and 6.4 kcal/mol. The observed doubly differential cross sections (DCSs) for the Cl(2P) product exhibit almost no collision energy dependence over this inspected energy range. The nearly forward-backward symmetric DCS indicates that the reaction proceeds predominantly on the ground-state potential energy surface at these energies. Variation of the forward-backward asymmetry with collision energy is interpreted using an osculating complex model. Although the potential energy surfaces obtained by CASSCF-MRCI ab initio calculations exhibit relatively low potential barriers of 1.6 and 6.5 kcal/mol for 1(1)A' and 2(1)A', respectively, the dynamics calculations indicate that contributions of these excited states are small at the collision energies lower than 15.0 kcal/mol. Theoretical DCSs calculated for the ground-state reaction pathway agree well with the observed ones. These experimental and theoretical results suggest that the titled reaction at collision energies less than 6.5 kcal/mol is predominantly via the ground electronic state.  相似文献   

5.
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相似文献   

6.
The H-atom abstraction reaction, O((3)P) + CH(4) → OH + CH(3), has been studied at a hyperthermal collision energy of 64 kcal mol(-1) by two crossed-molecular-beams techniques. The OH products were detected with a rotatable mass spectrometer employing electron-impact ionization, and the CH(3) products were detected with the combination of resonance-enhanced multiphoton ionization (REMPI) and time-sliced ion velocity-map imaging. The OH products are mainly formed through a stripping mechanism, in which the reagent O atom approaches the CH(4) molecule at large impact parameters and the OH product is scattered in the forward direction: roughly the same direction as the reagent O atoms. Most of the available energy is partitioned into product translation. The dominance of the stripping mechanism is a unique feature of such H-atom abstraction reactions at hyperthermal collision energies. In the hyperthermal reaction of O((3)P) with CH(4), the H-atom abstraction reaction pathway accounts for 70% of the reactive collisions, while the H-atom elimination pathway to produce OCH(3) + H accounts for the other 30%.  相似文献   

7.
The product state-resolved dynamics of the reactions H+H(2)O/D(2)O-->OH/OD((2)Pi(Omega);v',N',f )+H(2)/HD have been explored at center-of-mass collision energies around 1.2, 1.4, and 2.5 eV. The experiments employ pulsed laser photolysis coupled with polarized Doppler-resolved laser induced fluorescence detection of the OH/OD radical products. The populations in the OH spin-orbit states at a collision energy of 1.2 eV have been determined for the H+H(2)O reaction, and for low rotational levels they are shown to deviate from the statistical limit. For the H+D(2)O reaction at the highest collision energy studied the OD((2)Pi(3/2),v'=0,N'=1,A') angular distributions show scattering over a wide range of angles with a preference towards the forward direction. The kinetic energy release distributions obtained at 2.5 eV also indicate that the HD coproducts are born with significantly more internal excitation than at 1.4 eV. The OD((2)Pi(3/2),v'=0,N'=1,A') angular and kinetic energy release distributions are almost identical to those of their spin-orbit excited OD((2)Pi(1/2),v'=0,N'=1,A') counterpart. The data are compared with previous experimental measurements at similar collision energies, and with the results of previously published quasiclassical trajectory and quantum mechanical calculations employing the most recently developed potential energy surface. Product OH/OD spin-orbit effects in the reaction are discussed with reference to simple models.  相似文献   

8.
The dynamics of O((3)P) + CO(2) collisions at hyperthermal energies were investigated experimentally and theoretically. Crossed-molecular-beams experiments at = 98.8 kcal mol(-1) were performed with isotopically labeled (12)C(18)O(2) to distinguish products of nonreactive scattering from those of reactive scattering. The following product channels were observed: elastic and inelastic scattering ((16)O((3)P) + (12)C(18)O(2)), isotope exchange ((18)O + (16)O(12)C(18)O), and oxygen-atom abstraction ((18)O(16)O + (12)C(18)O). Stationary points on the two lowest triplet potential energy surfaces of the O((3)P) + CO(2) system were characterized at the CCSD(T)/aug-cc-pVTZ level of theory and by means of W4 theory, which represents an approximation to the relativistic basis set limit, full-configuration-interaction (FCI) energy. The calculations predict a planar CO(3)(C(2v), (3)A') intermediate that lies 16.3 kcal mol(-1) (W4 FCI excluding zero point energy) above reactants and is approached by a C(2v) transition state with energy 24.08 kcal mol(-1). Quasi-classical trajectory (QCT) calculations with collision energies in the range 23-150 kcal mol(-1) were performed at the B3LYP/6-311G(d) and BMK/6-311G(d) levels. Both reactive channels observed in the experiment were predicted by these calculations. In the isotope exchange reaction, the experimental center-of-mass (c.m.) angular distribution, T(θ(c.m.)), of the (16)O(12)C(18)O products peaked along the initial CO(2) direction (backward relative to the direction of the reagent O atoms), with a smaller isotropic component. The product translational energy distribution, P(E(T)), had a relatively low average of = 35 kcal mol(-1), indicating that the (16)O(12)C(18)O products were formed with substantial internal energy. The QCT calculations give c.m. P(E(T)) and T(θ(c.m.)) distributions and a relative product yield that agree qualitatively with the experimental results, and the trajectories indicate that exchange occurs through a short-lived CO(3)* intermediate. A low yield for the abstraction reaction was seen in both the experiment and the theory. Experimentally, a fast and weak (16)O(18)O product signal from an abstraction reaction was observed, which could only be detected in the forward direction. A small number of QCT trajectories leading to abstraction were observed to occur primarily via a transient CO(3) intermediate, albeit only at high collision energies (149 kcal mol(-1)). The oxygen isotope exchange mechanism for CO(2) in collisions with ground state O atoms is a newly discovered pathway through which oxygen isotopes may be cycled in the upper atmosphere, where O((3)P) atoms with hyperthermal translational energies can be generated by photodissociation of O(3) and O(2).  相似文献   

9.
We have determined limits on the cross section for both electronically nonadiabatic excitation and quenching in the Cl((2)P(j)) + D(2) system. Our experiment incorporates crossed-molecular-beam scattering with state-selective Cl((2)P(12,32)) detection and velocity-mapped ion imaging. By colliding atomic chlorine with D(2), we address the propensity for collisions that result in a change of the spin-orbit level of atomic chlorine either through electronically nonadiabatic spin-orbit excitation Cl((2)P(32)) + D(2)-->Cl(*)((2)P(12)) + D(2) or through electronically nonadiabatic spin-orbit quenching Cl(*)((2)P(12)) + D(2)-->Cl((2)P(32)) + D(2). In the first part of this report, we estimate an upper limit for the electronically nonadiabatic spin-orbit excitation cross section at a collision energy of 5.3 kcal/mol, which lies above the energy of the reaction barrier (4.9 kcal/mol). Our analysis and simulation of the experimental data determine an upper limit for the excitation cross section as sigma(NA)< or =0.012 A(2). In the second part of this paper we investigate the propensity for electronically nonadiabatic spin-orbit quenching of Cl(*) following a collision with D(2) or He. We perform these experiments at collision energies above and below the energy of the reaction barrier. By comparing the amount of scattered Cl(*) in our images to the amount of Cl(*) lost from the atomic beam we obtain the maximum cross section for electronically nonadiabatic quenching as sigma(NA)< or =15(-15) (+44) A(2) for a collision energy of 7.6 kcal/mol. Our experiments show the probability for electronically nonadiabatic quenching in Cl(*) + D(2) to be indistinguishable to that for the kinematically identical system of Cl(*) + He.  相似文献   

10.
The integral cross section of the S((1)D(2)) + H(2)(j = 0) → SH + H reaction has been measured for the first time at collision energies from 0.820 down to 0.078 kJ mol(-1) in a high-resolution crossed beam experiment. The excitation function obtained exhibits a non-monotonic variation with collision energy and compares well with the results of high-level quantum calculations. In particular, the structures observed in the lower energy part, where only a few partial waves contribute, can be described in terms of the sequential opening of individual channels, consistent with the theoretical calculations.  相似文献   

11.
We present results of time-dependent quantum mechanics (TDQM) and quasiclassical trajectory (QCT) studies of the excitation function for O(3P) + H2(v = 0-3,j = 0) --> OH + H from threshold to 30 kcal/mol collision energy using benchmark potential energy surfaces [Rogers et al., J. Phys. Chem. A 104, 2308 (2000)]. For H2(v = 0) there is excellent agreement between quantum and classical results. The TDQM results show that the reactive threshold drops from 10 kcal/mol for v = 0 to 6 for v = 1, 5 for v = 2 and 4 for v = 3, suggesting a much slower increase in rate constant with vibrational excitation above v = 1 than below. For H2(v > 0), the classical results are larger than the quantum results by a factor approximately 2 near threshold, but the agreement monotonically improves until they are within approximately 10% near 30 kcal/mol collision energy. We believe these differences arise from stronger vibrational adiabaticity in the quantum dynamics, an effect examined before for this system at lower energies. We have also computed QCT OH(v',j') state-resolved cross sections and angular distributions. The QCT state-resolved OH(v') cross sections peak at the same vibrational quantum number as the H2 reagent. The OH rotational distributions are also quite hot and tend to cluster around high rotational quantum numbers. However, the dynamics seem to dictate a cutoff in the energy going into OH rotation indicating an angular momentum constraint. The state-resolved OH distributions were fit to probability functions based on conventional information theory extended to include an energy gap law for product vibrations.  相似文献   

12.
A crossed molecular beam study is presented for the O((1)D(2))+HCl-->OH+Cl((2)P(J)) reaction at the collision energy of 6 kcal mol(-1). State-resolved doubly differential cross sections are obtained for the Cl((2)P(J=3/2) ) and Cl*((2)P(J=1/2) ) products by velocity-map ion imaging. Both products are slightly more forward scattered, which suggests a reaction mechanism without a long-lived intermediate in the ground electronic state. A small fraction (23 %) of the energy release into the translational degree of freedom indicates strong internal excitation of the counterpart OH radical. The contribution of the electronic excited states of O--HCl to the overall reaction is also examined from the doubly differential cross sections.  相似文献   

13.
Reliable thermochemical data for the reaction SO3 + H2O<-->SO3 x H2O (1a) are of crucial importance for an adequate modeling of the homogeneous H2SO4 formation in the atmosphere. We report on high-level quantum chemical calculations to predict the binding energy of the SO3 x H2O complex. The electronic binding energy is accurately computed to De = 40.9+/-1.0 kJ/mol = 9.8+/-0.2 kcal/mol. By using harmonic frequencies from density functional theory calculations (B3LYP/cc-pVTZ and TPSS/def2-TZVP), zero-point and thermal energies were calculated. From these data, we estimate D0 = -Delta H(1a)0(0 K) = 7.7+/-0.5 kcal/mol and Delta H(1a)0(298 K) = -8.3+/-1.0 kcal/mol.  相似文献   

14.
15.
Mechanisms for the deamination reaction of cytosine with H 2O/OH (-) and 2H 2O/OH (-) to produce uracil were investigated using ab initio calculations. Optimized geometries of reactants, transition states, intermediates, and products were determined at MP2 and B3LYP using the 6-31G(d) basis set and at B3LYP/6-31+G(d) levels of theory. Single point energies were also determined at MP2/G3MP2Large and G3MP2 levels of theory. Thermodynamic properties (Delta E, Delta H, and Delta G), activation energies, enthalpies, and free energies of activation were calculated for each reaction pathway investigated. Intrinsic reaction coordinate (IRC) analysis was performed to characterize the transition states on the potential energy surface. Seven pathways for the deamination reaction were found. All pathways produce an initial tetrahedral intermediate followed by several conformational changes. The final intermediate for all pathways dissociates to product via a 1-3 proton shift. The activation energy for the rate-determining step, the formation of the tetrahedral intermediate for pathway D, the only pathway that can lead to uracil, is 115.3 kJ mol (-1) at the G3MP2 level of theory, in excellent agreement with the experimental value (117 +/- 4 kJ mol (-1)).  相似文献   

16.
Reaction dynamics for a microsolvated SN2 reaction OH-(H2O)+CH3Cl have been investigated by means of the direct ab initio molecular dynamics method. The relative center-of-mass collision energies were chosen as 10, 15, and 25 kcal/mol. Three reaction channels were found as products. These are (1) a channel leading to complete dissociation (the products are CH3OH+Cl- +H2O: denoted by channel I), (2) a solvation channel (the products are Cl-(H2O)+CH3OH: channel II), and (3) a complex formation channel (the products are CH3OH...H2O+Cl-: channel III). The branching ratios for the three channels were drastically changed as a function of center-of-mass collision energy. The ratio of complete dissociation channel (channel I) increased with increasing collision energy, whereas that of channel III decreased. The solvation channel (channel II) was minor at all collision energies. The selectivity of the reaction channels and the mechanism are discussed on the basis of the theoretical results.  相似文献   

17.
Quantum mechanical wave packet calculations are carried out for the H((2)S) + FO((2)II) --> OH((2)II) + F((2)P) reaction on the adiabatic potential energy surface of the ground (3)A' triplet state. The state-to-state and state-to-all reaction probabilities for total angular momentum J = 0 have been calculated. The probabilities for J > 0 have been estimated from the J = 0 results by using J-shifting approximation based on a capture model. Then, the integral cross sections and initial state-selected rate constants have been calculated. The calculations show that the initial state-selected reaction probabilities are dominated by many sharp peaks. The reaction cross section does not manifest any sharp oscillations and the initial state-selected rate constants are sensitive to the temperature.  相似文献   

18.
Radical recombination reactions are important in the combustion of fuel oils. Shale oil contains alkylated heteroaromatic species, the simplest example of which is the 2-thienylmethyl radical. The ab initio potential energy surface for the reaction of the 2-thienylmethyl radical with the HO(2) radical has been examined. Seventeen product channels corresponding to either addition/elimination or direct hydrogen abstraction have been characterized for the first time. Direct hydrogen abstract from HO(2) proceeds via a weakly bound van der Waals complex, which leads to 2-methylthiophene, 2-methylene-2,3-dihydrothiophene, or 2-methylene-2,5-dihydrothiophene depending upon the 2-thienylmethyl radical reaction site. The addition pathway for the two radical reactants is barrierless with the formation of three adducts, as distinguished by HO(2) reaction at three different sites on the 2-thienylmethyl radical. The addition is exothermic by 37-55 kcal mol(-1) relative to the entrance channel, and these excess energies are available to promote further decomposition or rearrangement of the adducts, leading to nascent products such as H, OH, H(2)O, and CH(2)O. The reaction surfaces are characterized by relatively low barriers (most lower than 10 kcal mol(-1)). Upon the basis of a careful analysis of the overall barrier heights and reaction exothermicities, the formations of O(2), OH, and H(2)O are likely to be important pathways in the radical recombination reactions of 2-thienylmethyl + HO(2).  相似文献   

19.
We have performed quantum mechanical (QM) dynamics calculations within the independent-state approximation with new benchmark triplet A" and A' surfaces [B. Ramachandran et al., J. Chem. Phys. 119, 9590 (2003)] for the rovibronic state-to-state measurements of the reaction O(3P)+HCl(v=2,j=1,6,9)-->OH(v'j')+Cl(2P) [Zhang et al., J. Chem. Phys. 94, 2704 (1991)]. The QM and experimental rotational distributions peak at similar OH(j') levels, but the QM distributions are significantly narrower than the measurements and previous quasiclassical dynamics studies. The OH(low j) populations observed in the measurements are nearly absent in the QM results. We have also performed quasiclassical trajectory with histogram binning (QCT-HB) calculations on these same benchmark surfaces. The QCT-HB rotational distributions, which are qualitatively consistent with measurements and classical dynamics studies using other surfaces, are much broader than the QM results. Application of a Gaussian binning correction (QCT-GB) dramatically narrows and shifts the QCT-HB rotational distributions to be in very good agreement with the QM results. The large QCT-GB correction stems from the special shape of the joint distribution of the classical rotational/vibrational action of OH products. We have also performed QM and QCT calculations for the transition, O+HCl(v=0,T=300 K)-->OH(v'j')+Cl from threshold to approximately 130 kcal mol(-1) collision energy as a guide for possible future hyperthermal O-atom measurements. We find in general a mixed energy release into translation and rotation consistent with a late barrier to reaction. Angular distributions at high collision energy are forward peaked, consistent with a stripping mechanism. Direct collisional excitation channel cross sections, O+HCl(v=0,T=300 K)-->O+HCl(v'=1), in the same energy range are large, comparable in magnitude to the reactive channel cross sections. Although the (3)A" state dominates most collision processes, above approximately 48 kcal mol(-1), the (3)A' state plays the major role in collisional excitation.  相似文献   

20.
This article describes an ab initio investigation on the potential surfaces for one of the simplest hydrogen atom abstraction reactions, that is, HO + H2O <--> H2O + OH. In accord with the findings in the previously reported theoretical investigations, two types of the hydrogen-bonding complexes [HOH--OH] and [H2O--HO] were located on the potential energy surface. The water molecule acts as a hydrogen donor in the [HOH--OH] complex, while the OH radical acts as a hydrogen donor in the [H2O--HO] complex. The energy evaluations at the MP2(FC) basis set limit, as well as those through the CBS-APNO procedure, have provided estimates for enthalpies of association for these complexes at 298 K as -2.1 approximately -2.3 and -4.1 approximately -4.3 kcal/mol, respectively. The IRC calculations have suggested that the [H2O--HO] complex should be located along the reaction coordinate for the hydrogen abstraction. Our best estimate for the classical barrier height for the hydrogen abstraction is 7.8 kcal/mol, which was obtained from the CBS-APNO energy evaluations. After fitting the CBS-APNO potential energy curve to a symmetrical Eckart function, the rate constants were calculated by using the transition state theory including the tunneling correction. Our estimates for the Arrhenius parameters in the temperature region from 300 to 420 K show quite reasonable agreement with the experimentally derived values.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号