Communication: state-to-state quantum dynamics study of the OH + CO → H + CO2 reaction in full dimensions (J = 0)

A full dimensional state-to-state quantum dynamics study is carried out for the prototypical complex-formation OH + CO → H + CO(2) reaction in the ground rovibrational initial state on the Lakin-Troya-Schatz-Harding potential energy surface by using the reactant-product decoupling method. With three heavy atoms and deep wells on the reaction path, the reaction represents a huge challenge for accurate quantum dynamics study. This state-to-state calculation is the first such a study on a four-atom reaction other than the H(2) + OH ? H(2)O + H and its isotope analogies. The product CO(2) vibrational and rotational state distributions, and product energy partitioning information are presented for ground initial rovibrational state with the total angular momentum J = 0.     

Mechanistic and kinetic study of the O + CH2OH reaction

The mechanism for the O + CH2OH reaction was investigated by various ab initio quantum chemistry methods. For the chemical activation mechanism, that is, the addition/elimination path, the couple-cluster methods including CCSD and CCSD(T) were employed with the cc-pVXZ (X = D, T, Q, 5) basis sets. For the abstraction channels, multireference methods including CASSCF, CASPT2, and MRCISD were used with the cc-pVDZ and cc-pVTZ basis sets. It has been shown that the production of H + HCOOH is the major channel in the chemical activation mechanism. The minor channels include HCO + H2O and OH + CH2O. The hydrogen abstraction by an O atom from the CH2OH radical produces either OH + CH2O or OH + HCOH. Moreover, the two abstraction reactions are essentially barrierless processes. The rate constants for the association of O with CH2OH have been calculated using the flexible transition state theory. A weak negative temperature dependence of the rate constants is found in the range 250-1000 K. Furthermore, it is estimated that the abstraction processes also play an important role in the O + CH2OH reaction. Additionally, the falloff behavior for the OCH2OH --> H + HCOOH reaction has been investigated. The present theoretical results are compared to the experimental measurements to understand the mechanism and kinetic behavior of the O + CH2OH reaction and the unimolecular reaction of the OCH2OH radical.     

Quantum mechanical rate constants for H + O2 <--> O + OH and H + O2 --> HO2 reactions

Canonical rate constants for both the forward and reverse H + O(2) <--> O + OH reactions were calculated using a quantum wave packet-based statistical model on the DMBE IV potential energy surface of Varandas and co-workers. For these bimolecular reactions, the results show reasonably good agreement with available experimental and theoretical data up to 1500 K. In addition, the capture rate for the H + O(2) --> HO(2) addition reaction at the high-pressure limit was obtained on the same potential using a time-independent quantum capture method. Excellent agreement with experimental and quasi-classical trajectory results was obtained except for at very low temperatures, where a reaction threshold was found and attributed to the centrifugal barrier of the orbital motion.     

Ultracold collisions and reactions of vibrationally excited OH radicals with oxygen atoms

We report a quantum dynamics study of O + OH (v = 1, j = 0) collisions on its ground electronic state, employing two different potential energy surfaces: the DIMKP surface by Kendrick and Pack, and the XXZLG surface by Xu et al. A time-independent quantum mechanical method based on hyperspherical coordinates has been adopted for the dynamics calculations. Energy-dependent probabilities and rate coefficients are computed for the elastic, inelastic, and reactive channels over the collision energy range E(coll) = 10(-10)-0.35 eV, for J = 0 total angular momentum. Initial state-selected reaction rate coefficients are also calculated from the J = 0 reaction probabilities by applying a J-shifting approximation, for temperatures in the range T = 10(-6)-700 K. Our results show that the dynamics of the collisional process and its outcome are strongly influenced by long-range forces, and chemical reactivity is found to be sensitive to the choice of the potential energy surface. For O + OH (v = 1, j = 0) collisions at low temperatures, vibrational relaxation of OH competes with reactive scattering. Since long-range interactions can facilitate vibrational relaxation processes, we find that the DIMKP potential (which explicitly includes van der Waals dispersion terms) favours vibrational relaxation over chemical reaction at low temperatures. On the DIMKP potential in the ultracold regime, the reaction rate coefficient for O + OH (v = 1, j = 0) is found to be a factor of thirteen lower than that for O + OH (v = 0, j = 0). This significantly high reactivity of OH (v = 0, j = 0), compared to that of OH (v = 1, j = 0), is attributed to enhancement caused by the presence of a HO(2) quasibound state (scattering resonance) with energy near the O + OH (v = 0, j = 0) dissociation threshold. In contrast, the XXZLG potential does not contain explicit van der Waals terms, being just an extrapolation by a nearly constant function at large O-OH distances. Therefore, long-range potential couplings are absent in calculations using the XXZLG surface, which does not induce vibrational relaxation as efficiently as the DIMKP potential. The XXZLG potential leads to a slightly higher reactivity (a factor of 1.4 higher) for O + OH (v = 1, j = 0) compared to that for O + OH (v = 0, j = 0) at ultracold temperatures. Overall, both potential surfaces yield comparable values of reaction rate coefficients at low temperatures for the O + OH (v = 1, j = 0) reaction.     

CH3与NO2的反应

The elementary reaction of the CH3 radical with NO2 was investigated by time-resolved FTIR spectroscopy and quantum chemical calculations. The CH3 radical was produced by laser photolysis of CH3Br or CH3I at 248 nm. Vibrationally excited products OH, HNO and CO2 were observed by the time-resolved spectroscopy for the first time. The formation of another product NO was also verified. According to these observations, the product channels leading to CH3O+NO, CH2NO+OH and HNO+H2CO were identified. The channel of CH3O+NO was the major one. The reaction mechanisms of the above channels were studied by quantum chemical calculations at CCSD(T)/6-311++G(df,p)//MP2/6-311G(d,p) level. The calculated results fit with the experimental observations well.     

Isotope branching and tunneling in O(3P)+HD-->OH+D; OD+H reactions

The O((3)P)+HD and O((3)P)+D(2) reactions are studied using quantum scattering calculations and chemically accurate potential energy surfaces developed for the O((3)P)+H(2) system by Rogers et al. [J. Phys. Chem. A 104, 2308 (2000)]. Cross sections and rate coefficients for OH and OD products are calculated using accurate quantum methods as well as the J-shifting approximation. The J-shifting approach is found to work remarkably well for both O+HD and O+D(2) collisions. The reactions are dominated by tunneling at low temperatures and for the O+HD reaction the hydrogen atom transfer leading to the OH product dominates at low temperatures. Our result for the OH/OD branching ratio is in close agreement with previous calculations over a wide range of temperatures. The computed OH/OD branching ratios are also in close agreement with experimental results of Robie et al. [Chem. Phys. Lett. 134, 579 (1987)] at temperatures above 400 K but the theoretical results do not reproduce the rapid rise in the experimental values of the branching ratio for temperatures lower than 350 K. We believe that new measurements could resolve the long-standing discrepancy between experiment and theory for this benchmark reaction.     

Theoretical and experimental study of the product branching in the reaction of acetic acid with OH radicals

The product distribution of the reaction of acetic acid, CH(3)COOH, with hydroxyl radicals, OH, was studied experimentally and theoretically. Mass-spectrometric measurements at 290 K and 2 Torr of He of the CO(2) yield versus the loss of acetic acid yielded a branching fraction of 64 +/- 14% for the abstraction of the acidic hydrogen as follows: CH(3)COOH + OH --> CH(3)COO + H(2)O --> CH(3) + CO(2) + H(2)O. A quantum chemical and theoretical kinetic analysis showed that the abstraction of the acidic hydrogen is enhanced relative to the abstraction of -CH(3) hydrogens because of the formation of a strong pre-reactive H-bonded complex, where the H-bonds are retained in the H-abstraction transition state. The potential energy surface of the reaction is explored in detail, and the reaction products of the individual channels are identified. The theoretical product branching is found to be critically dependent on the energetic and rovibrational differences between the H-abstraction transition states.     

掺Cu对MoO3-TiO2/SiO2上光促甲烷和水表面反应的影响

 在固定床环隙反应器中,借助紫外光的激发,气相甲烷和水在MoO3-TiO2/SiO2催化剂表面生成了甲醇和氢气,当在催化剂中掺杂Cu2+后,光催化剂的活性明显提高. XRD,IR,UVDRS和TPD的研究结果表明,在催化剂表面形成了具有Mo-O-Ti和Mo-O-Cu基元的高度分散物种,不但使得吸光带边明显蓝移,而且扩展了催化剂的光响应范围. 所形成的复合结构还可以优化单组分的吸光性能并促进对反应物分子的吸附活化,同时可以有效地转移光生电子和空穴. 掺杂Cu2+能够进一步延长光生电子-空穴对的寿命,进而提高反应的量子产率.     

A theoretical study of the fluoride exchange between UO2F+(aq) and UO2 2+(aq)

Experimental data on the thermodynamics and reaction mechanism of the inner-sphere fluoride exchange reaction U17O2(2+) + UO2F+ <==> U17O2F+ + UO2(2+) have been compared with different intimate reaction mechanisms using quantum chemical methods. Two models have been tested that start from the outer sphere complexes, (H2O)[U(A)O2F(OH2)4+]...[U(B)O2(OH2)5(2+)] and [U(A)O2F(OH2)4+]...[U(B)O2(OH2)5(2+)]; the geometry and energies of the intermediates and transition states along possible reaction pathways have been calculated using different ab initio methods, SCF, B3LYP and MP2. Both the experimental data and the theoretical results suggest that the fluoride exchange takes place via the formation and breaking of a U-F-U bridge that is the rate determining step. The calculated activation enthalpy DeltaH( not equal) = 30.9 kJ mol(-1) is virtually identical to the experimental value 31 kJ mol(-1); however this agreement may be a coincidence as we do not expect a larger accuracy than 10 kJ mol(-1) with the methods used. The calculations show that the fluoride bridge is formed as an insertion of U(A)O2)F(OH2)4+ into U(B)O2(OH2)5(2+) followed by a subsequent transfer of water from the first to the second coordination sphere of U(B).     

Quasi-classical Trajectory Study of the Intramolecular Isotope Effect in the Reaction O（3p）＋H2/HD

Theoretical studies of the dynamics of the reactions O（3p）＋H2/HD（ν=0, j=0）→OH＋H have been performed with quasi-classical trajectory method （QCT） on an ab initio potential surface for the lowest triplet electronic state of H2O（aA＂）. The QCT-calculated integral cross sections are in good agreement with the earlier time-dependent quantum mechanics results. The state-resolved rotational distributions reveal that the product OH rotational distributions for O＋HD have a preference for populating highly internally excited states compared with the O＋H2 reaction. Distributions of differential cross sections show that directions of scattering are strongly dependent on the choice of quantum state. The polarization dependent generalized differential cross-sections and the distributions were calculated and a pronounced isotopic effect is revealed. The calculated results indicate that the product polarization is very sensitive to the mass factor.     

13C, 18O, and D fractionation effects in the reactions of CH3OH isotopologues with Cl and OH radicals

A relative rate experiment is carried out for six isotopologues of methanol and their reactions with OH and Cl radicals. The reaction rates of CH2DOH, CHD2OH, CD3OH, (13)CH3OH, and CH3(18)OH with Cl and OH radicals are measured by long-path FTIR spectroscopy relative to CH3OH at 298 +/- 2 K and 1013 +/- 10 mbar. The OH source in the reaction chamber is photolysis of ozone to produce O((1)D) in the presence of a large excess of molecular hydrogen: O((1)D) + H2 --> OH + H. Cl is produced by the photolysis of Cl2. The FTIR spectra are fitted using a nonlinear least-squares spectral fitting method with measured high-resolution infrared spectra as references. The relative reaction rates defined as alpha = k(light)/k(heavy) are determined to be: k(OH + CH3OH)/k(OH + (13)CH3OH) = 1.031 +/- 0.020, k(OH + CH3OH)/k(OH + CH3(18)OH) = 1.017 +/- 0.012, k(OH + CH3OH)/k(OH + CH2DOH) = 1.119 +/- 0.045, k(OH + CH3OH)/k(OH + CHD2OH) = 1.326 +/- 0.021 and k(OH + CH3OH)/k(OH + CD3OH) = 2.566 +/- 0.042, k(Cl + CH3OH)/k(Cl + (13)CH3OH) = 1.055 +/- 0.016, k(Cl + CH3OH)/k(Cl + CH3(18)OH) = 1.025 +/- 0.022, k(Cl + CH3OH)/k(Cl + CH2DOH) = 1.162 +/- 0.022 and k(Cl + CH3OH)/k(Cl + CHD2OH) = 1.536 +/- 0.060, and k(Cl + CH3OH)/k(Cl + CD3OH) = 3.011 +/- 0.059. The errors represent 2sigma from the statistical analyses and do not include possible systematic errors. Ground-state potential energy hypersurfaces of the reactions were investigated in quantum chemistry calculations at the CCSD(T) level of theory with an extrapolated basis set. The (2)H, (13)C, and (18)O kinetic isotope effects of the OH and Cl reactions with CH3OH were further investigated using canonical variational transition state theory with small curvature tunneling and compared to experimental measurements as well as to those observed in CH4 and several other substituted methane species.     

Quasi-classical Trajectory Study of F+H2O→HF+OH Reaction: Influence of Barrier Height,Reactant Rotational Excitation,and Isotopic Substitution (cited: 4)

The reaction dynamics of the F+H₂O/D₂O→HF/DF+OH/OD are investigated on an ac-curate potential energy surface (PES) using a quasi-classical trajectory method. For bothisotopomers, the hydrogen/deuterium abstraction reaction is dominated by a direct rebound mechanism over a very low “reactant-like” barrier, which leads to a vibrationally hot HF/DF product with an internally cold OH/OD companion. It is shown that the lowered reaction barrier on this PES, as suggested by high-level ab initio calculations, leads to a much better agreement with the experimental reaction cross section, but has little impact on the product state distributions and mode selectivity. Our results further indicate that rotational exci-tation of the H₂O reactant leads to significant enhancement of the reactivity, suggesting a strong coupling with the reaction coordinate.     

Differential and integral cross sections for the H + O2 --> OH + O combustion reaction

We present accurate differential and integral cross sections for the H + O2 --> OH + O reaction obtained on a newly developed ab initio potential energy surface using time-independent and time-dependent quantum mechanical methods. The product angular distributions near the reaction threshold show pronounced forward and backward peaks, reflecting the complex-forming mechanism. However, the asymmetry of these peaks suggests certain nonstatistical behaviors, presumably due to some relatively short-lived resonances. The integral cross section increases monotonically with the collision energy above a reaction threshold.     

Using redundant coordinates to represent potential energy surfaces with lower-dimensional functions

We propose a method for fitting potential energy surfaces with a sum of component functions of lower dimensionality. This form facilitates quantum dynamics calculations. We show that it is possible to reduce the dimensionality of the component functions by introducing new and redundant coordinates obtained with linear transformations. The transformations are obtained from a neural network. Different coordinates are used for different component functions and the new coordinates are determined as the potential is fitted. The quality of the fits and the generality of the method are illustrated by fitting reference potential surfaces of hydrogen peroxide and of the reaction OH+H(2)-->H(2)O+H.     

Investigation of the radical product channel of the CH3COO2 + HO2 reaction in the gas phase

The reaction of CH(3)C(O)O(2) with HO(2) has been investigated at 296 K and 700 Torr using long path FTIR spectroscopy, during photolysis of Cl(2)/CH(3)CHO/CH(3)OH/air mixtures. The branching ratio for the reaction channel forming CH(3)C(O)O, OH and O(2) (reaction ) has been determined from experiments in which OH radicals were scavenged by addition of benzene to the system, with subsequent formation of phenol used as the primary diagnostic for OH radical formation. The dependence of the phenol yield on benzene concentration was found to be consistent with its formation from the OH-initiated oxidation of benzene, thereby confirming the presence of OH radicals in the system. The dependence of the phenol yield on the initial peroxy radical precursor reagent concentration ratio, [CH(3)OH](0)/[CH(3)CHO](0), is consistent with OH formation resulting mainly from the reaction of CH(3)C(O)O(2) with HO(2) in the early stages of the experiments, such that the limiting yield of phenol at high benzene concentrations is well-correlated with that of CH(3)C(O)OOH, a well-established product of the CH(3)C(O)O(2) + HO(2) reaction (via channel (3a)). However, a delayed source of phenol was also identified, which is attributed mainly to an analogous OH-forming channel of the reaction of HO(2) with HOCH(2)O(2) (reaction ), formed from the reaction of HO(2) with product HCHO. This was investigated in additional series of experiments in which Cl(2)/CH(3)OH/benzene/air and Cl(2)/HCHO/benzene/air mixtures were photolysed. The various reaction systems were fully characterised by simulations using a detailed chemical mechanism. This allowed the following branching ratios to be determined: CH(3)C(O)O(2) + HO(2)--> CH(3)C(O)OOH + O(2), k(3a)/k(3) = 0.38 +/- 0.13; --> CH(3)C(O)OH + O(3), k(3b)/k(3) = 0.12 +/- 0.04; --> CH(3)C(O)O + OH + O(2), k(3c)/k(3) = 0.43 +/- 0.10: HOCH(2)O(2) + HO(2)--> HCOOH + H(2)O + O(2), k(17b)/k(17) = 0.30 +/- 0.06; --> HOCH(2)O + OH + O(2), k(17c)/k(17) = 0.20 +/- 0.05. The results therefore provide strong evidence for significant participation of the radical-forming channels of these reactions, with the branching ratio for the title reaction being in good agreement with the value reported in one previous study. As part of this work, the kinetics of the reaction of Cl atoms with phenol (reaction (14)) have also been investigated. The rate coefficient was determined relative to the rate coefficient for the reaction of Cl with CH(3)OH, during the photolysis of mixtures of Cl(2), phenol and CH(3)OH, in either N(2) or air at 296 K and 760 Torr. A value of k(14) = (1.92 +/- 0.17) x 10(-10) cm(3) molecule(-1) s(-1) was determined from the experiments in N(2), in agreement with the literature. In air, the apparent rate coefficient was about a factor of two lower, which is interpreted in terms of regeneration of phenol from the product phenoxy radical, C(6)H(5)O, possibly via its reaction with HO(2).     

Infrared spectra of the OH+ and H2O+ cations solvated in solid argon

Infrared spectra of various OH+ and H2O+ isotopomers solvated in solid argon are presented. The OH+ and H2O+ cations were produced by co-deposition of H2O/Ar mixture with high-frequency discharged Ar at 4 K. Detailed isotopic substitution studies confirm the assignments of absorptions at 3054.9 and 3040.0 cm(-1) to the antisymmetric and symmetric H-O-H stretching vibrations of H2O+ and 2979.6 cm(-1) to the O-H stretching vibration of OH+. The frequencies of H2O+ solvated in solid argon are red-shifted, whereas the frequency of OH+ is blue-shifted with respect to the gas-phase fundamentals. On the basis of previous gas-phase studies and quantum chemical calculations, the OH+ and H2O+ cations solvated in solid argon may be regarded as the OH+-Ar5 and H2O+-Ar4 complexes isolated in the argon matrix.     

气体再燃低NOx燃烧中NO与NHi的反应机理研究

采用量子化学密度泛函理论(DFT)对NO与NH_i自由基的反应机理进行了研究,并结合经典过渡态理论对各反应速率常数进行了计算。结果表明,NO与NH₂自由基的反应体系可通过六个反应通道形成N₂+H₂O、N₂O+H₂和N₂H+OH。从能量变化和反应速率两方面考虑,产物N₂+H₂O最容易生成,其最佳反应通道为NO+NH₂→→N₂+H₂O;NO与NH自由基的反应体系可通过七个反应通道形成N₂+OH、N₂O+H和N₂H+O;其中,N₂+OH最容易生成,最佳反应通道为NO+NH→→N₂+OH。比较发现, NH比NH₂自由基更易与NO发生反应生成N₂。因此,在实际运行中改变操作条件,实现NH₂等向NH方向转化,有利于NO_x的还原。     

Measurements and automated mechanism generation modeling of OH production in photolytically initiated oxidation of the neopentyl radical

Production of OH in the reaction of the neopentyl radical with O2 has been measured by a laser photolysis/cw absorption method for various pressures and oxygen concentrations at 673, 700, and 725 K. The MIT Reaction Mechanism Generator (RMG) was used to automatically generate a model for this system, and the predicted OH concentration profiles are compared to present and literature experimental results. Several reactions significantly affect the OH profile. The experimental data provide useful constraints on the rate coefficient for the formally direct chemical activation reaction of neopentyl radical with O2 to form OH (CH3)3CCH2 + O2 --> OH + 3,3-dimethyloxetane (Rxn 1) At 673 K and 60 Torr, log k(1) (cm(3) molecule(-1) s(-1)) = -13.7 +/- 0.5. Absolute absorbance measurements on OH and I indicate that the branching ratio for R + O2 to OH is about 0.03 under these conditions. The data suggest that the ab initio neopentyl + O2 potential energy surface of Sun and Bozzelli is accurate to within 2 kcal mol(-1).