O(3P)+O2H→OH+O2反应机理的密度泛函理论研究

用密度泛函理论方法研究了O(3P)与O2H反应生成羟基和氧分子的反应机理. 在PW91/6-31+G水平上用梯度解析技术全自由度优化上述反应物、产物和反应路径上的中间体及过渡态几何构型, 并通过频率振动分析加以确认, 计算IRC反应路径及中间体异构化过程, 确定了此反应的可能反应通道. 结果表明: 该反应是多通道多步骤的强放热反应. 首先形成顺式或反式O3H富能中间体, 此过程无能垒; 然后跨过一个能垒分解成产物OH和O2. 通道IM1→TS1比IM2→TS2克服的能垒要大, 反应放热372.822 kJ*mol-1. IM1TS3IM2 可相互转化.     

Nonadiabatic quantum dynamics of C(1D)+H2→CH+H: coupled-channel calculations including Renner-Teller and Coriolis terms

The Renner-Teller (RT) coupled-channel dynamics for the C((1)D)+H(2)(X(1)Σ(g) (+))→CH(X(2)Π)+H((2)S) reaction has been investigated for the first time, considering the first two singlet states ??(1)A' and b(1)A' of CH(2) dissociating into the products and RT couplings, evaluated through the ab initio matrix elements of the electronic angular momentum. We have obtained initial-state-resolved probabilities, cross sections and thermal rate constants via the real wavepacket method for both coupled electronic states. In contrast to the N((2)D)+H(2)(X(1)Σ(g)(+)) system, RT effects tend to reduce probabilities, cross sections, and rate constants in the low energy range compared to Born-Oppenheimer (BO) ones, due to the presence of a repulsive RT barrier in the effective potentials and to long-lived resonances. Furthermore, contrary to BO results, the rate constants have a positive temperature dependence in the 100-400 K range. The two-state RT rate constant at 300 K, lower than the BO one, remains inside the error bars of the experimental value.     

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.     

O(3P)+O2H——OH+O2反应机理的密度泛函理论研究

用密度泛函理论方法研究了O(³P)与O₂H反应生成羟基和氧分子的反应机理.在PW91/6-31+G^*水平上用梯度解析技术全自由度优化上述反应物、产物和反应路径上的中间体及过渡态几何构型,并通过频率振动分析加以确认,计算IRC反应路径及中间体异构化过程,确定了此反应的可能反应通道.结果表明:该反应是多通道多步骤的强放热反应.首先形成顺式或反式O₃H富能中间体,此过程无能垒;然后跨过一个能垒分解成产物OH和O₂.通道IM1→TS1比IM2→TS2克服的能垒要大,反应放热372.822kJ.mol^-1.IM1TS3IM2可相互转化.     

Abstraction versus insertion in the reaction O(1D2)+H2→OH(υ″ = 2) + H

The product hydroxyl radical arising from the reaction O(¹D₂)+H₂→OH+H was detected by LIF following excitation of the off-diagonal transition OH(A²Σ⁺, υ′=1←X²Π, υ″=2) in the region 348–357 nm. The rotational population distribution in υ″=2 appears to be inverted and quite similar to that previously reported for υ″=0 and 1. Because rotationally cool OH was not observed, there is no evidence for the existence of an abstractive pathway in which the subject reaction occurs without the initial formation, via insertion, of a chemically activated HOH collision complex.     

Quasiclassical trajectory study of the reaction O(3P) + HI → OH + I

Three-dimensional quasiclassical trajectory calculations were carried out for the reaction of oxygen atoms O(³P) with hydrogen iodide molecules (HI and DI) for the temperature range 200–550 K, using a LEPS potential-energy surface. The calculated results include reaction cross sections, rate constants, kinetic isotope effects, the influence of vibrational and rotational excitation of the reactants on the dynamics, and the product energy partitioning and angular distribution. The calculated results are in good agreement with the available experimental results. The dynamics of the O + HI reaction is discussed in view of the associated mass combination H + LH′ (H and H′ are heavy atoms and L is a light atom), and in relation to earlier trajectory results for the reactions O + HCl and O + HBr.     

HNCO+OH→H2O+NCO的反应机理

采用从头算分子轨道法(UHF/6-31G**水平,并用MP4加以相关能校正)研究了HNCO+OH→H2O+NCO反应机理.同时用Morokuma数值法获得了反应途径即内禀反应坐标(IRC).沿着IRC,运用反应途径哈密顿理论,获得反应途径动态学信息.在此基础上,根据过渡态理论和相应隧道效应校正,计算了在不同温度下的反应速率常数,得到了和实验相一致的结果.计算结果表明,此反应是一步直接型的抽提H反应.     

Nonadiabatic quantum reactive scattering of the OH(A 2Σ+) + D2

The seams of conical intersection exist between the ground (1 (2)A(')) and the first-excited (2 (2)A(')) electronic potential energy surfaces (PESs) of OH(A (2)Σ(+),X (2)Π) + H(2) system. This intersection induces the nonadiabatic quenching of OH(A (2)Σ(+)) by D(2). We present nonadiabatic quantum dynamics study for OH(A (2)Σ(+)) + D(2) on new five-dimensional coplanar PESs. The ab initio calculations of PESs are based on multireference configuration interaction (MRCI)/aug-cc-pVQZ level. A back-propagation neural network is utilized to fit the PESs and nonadiabatic coupling. High degrees of rotational excitation of quenched OH(X (2)Π) products are found in nonreactive quenching channel, and the quenched D(2) products are vibrationally excited up to quantum number v(2) (')=8. The theoretical results of nonadiabatic time-dependent wave-packet calculation are in good agreement with the existing experimental data.     

State-to-state reactive scattering in six dimensions using reactant-product decoupling: OH + H2 → H2O + H (J = 0)

We extend to full dimensionality a recently developed wave packet method [M. T. Cvitas? and S. C. Althorpe, J. Phys. Chem. A 113, 4557 (2009)] for computing the state-to-state quantum dynamics of AB + CD → ABC + D reactions and also increase the computational efficiency of the method. This is done by introducing a new set of product coordinates, by applying the Crank-Nicholson approximation to the angular kinetic energy part of the split-operator propagator and by using a symmetry-adapted basis-to-grid transformation to evaluate integrals over the potential energy surface. The newly extended method is tested on the benchmark OH + H(2) → H(2)O + H reaction, where it allows us to obtain accurately converged state-to-state reaction probabilities (on the Wu-Schatz-Fang-Lendvay-Harding potential energy surface) with modest computational effort. These methodological advances will make possible efficient calculations of state-to-state differential cross sections on this system in the near future.     

CH(^4∑^-)+H2O→CH2(^3B1)+OH的直接动态学研究

用从头算方法计算了反应CH(^4∑^-)+H2O→CH2(^3B1)+OH的反应途径。在此基础上, 计算沿反应途径的动态学性质和正则变分过渡态理论的速率常数, 并进行隧道效应校正。结果表明, 电子相关能对反应活化位垒影响较大; 反应存在返回效应和隧道效应, 用正则变分过渡态方法和小曲率近似的隧道校正是有效的。     

Quantum chemical and master equation study of OH + CH2O → H2O + CHO reaction rates in supercritical CO2 environment

We investigated the reaction rates of OH + CH₂O → H₂O + CHO at CO₂ pressures of up to 1000 atm with and without CO₂ molecule included in a reactive complex. Both mechanisms begin with formation of the hydrogen-bonded prereactive complexes. Our ab initio calculations indicate a possibility of catalytic effect, predicting an activation barrier that one order of magnitude lower when the CO₂ molecule is involved. To verify this effect, we use the Rice–Ramsperger–Kassel–Marcus theory and solve unimolecular master equations in the steady-state approximation. We assume the equilibrium between prereactive complexes and reactants and compare the bimolecular reaction rates for the two mechanisms. The catalyzed reaction mechanism is found to be faster at higher CO₂ pressures and lower temperatures, when prereactive complexes have nonnegligible concentration. Therefore, this catalytic effect may be important for this reactive process in room temperature supercritical CO₂ solvent, but is unlikely to play a role during oxy-combustion.     

Quantum dynamics of the S+OH→SO+H reaction

First accurate quantum mechanical scattering calculations have been carried out for the S((3)P)+OH(X?(2)Π)→SO(X?(3)Σ(-))+H((2)S) reaction using a recent ab initio potential energy surface for the ground electronic state, X?(2)A("), of HSO. Total and state-to-state reaction probabilities for a total angular momentum J=0 have been determined for collision energies up to 0.5 eV. A rate constant has been calculated by means of the J-shifting approach in the 10-400 K temperature range. Vibrational and rotational product distributions show no specific behavior and are consistent with a mixture of direct and indirect reaction mechanisms.     

CHF3…H2O complex revisited: a matrix isolation and ab initio study

Structural Chemistry - The structural and spectroscopic features of the CHF3…H2O complex have been investigated using high-level ab initio calculations and IR matrix isolation spectroscopy....     

State-to-state quantum dynamics of the H((2)S) + O2(ã(1)Δ(g)) → O((3)P)+OH(X(2)Π) reaction on the first excited state of HO2(Ã(2)A')

State-to-state differential and integral cross sections for the title reaction were calculated using an exact wave packet method on a recently developed ab initio potential energy surface of the first excited state HO(2)(?(2)A'). The calculation results indicate that the reaction is dominated by highly rotationally excited OH products scattered in both the forward and backward directions, consistent with the formation of a long-lived HO(2) reaction intermediate. However, a statistical model was found to overestimate the integral cross sections, due apparently to dynamical bottlenecks. In addition, a unique feature in the OH + O exit channel potential promotes rotational excitation of the departing OH product by exerting a torque force. The role of the title reaction in high temperature combustion is also discussed.     

Ab initio potential energy surface and quantum dynamics for the H + CH4 → H2 + CH3 reaction

A new full-dimensional potential energy surface for the title reaction has been constructed using the modified Shepard interpolation scheme. Energies and derivatives were calculated using the UCCSD(T) method with aug-cc-pVTZ and 6-311++G(3df,2pd) basis sets, respectively. A total number of 30,000 data points were selected from a huge number of molecular configurations sampled by trajectory method. Quantum dynamical calculations showed that the potential energy surface is well converged for the number of data points for collision energy up to 2.5 eV. Total reaction probabilities and integral cross sections were calculated on the present surface, as well as on the ZBB3 and EG-2008 surfaces for the title reaction. Satisfactory agreements were achieved between the present and the ZBB3 potential energy surfaces, indicating we are approaching the final stage to obtain a global potential energy surface of quantitative accuracy for this benchmark polyatomic system. Our calculations also showed that the EG-2008 surface is less accurate than the present and ZBB3 surfaces, particularly in high energy region.     

Calculated energy barriers for the identity SN2 reaction H2O + CH3OH2+ → +H2OCH3 + OH2 in the gas phase,in water clusters,and in aqueous solution

The bimolecular nucleophilic substitution reaction H₂O + CH₃OH₂⁺ → ⁺H₂OCH₃ + OH₂ has been studied using various quantum chemical methods. Accurate barriers for the reaction in the gas phase are presented and discussed. The effect of microsolvation by water molecules in small clusters has been investigated. Extrapolation of the barrier obtained in the small clusters, using a linear relationship between the activation energy and the proton affinity of water clusters, gives a barrier for the reaction in aqueous solution which is in good agreement with that obtained in separate model calculations (polarized continuum model of a super molecule with the first solvation shell included).     

State-to-state quantum dynamics of the N(4S) + OH(X2Π) → H(2S) + NO(X 2Π) reaction

Quantum state-to-state dynamics for the N((4)S) + OH(X(?2)Π) → H((2)S) + NO(X(?2)Π) reaction is reported on an accurate ab initio potential energy surface of the lowest triplet electronic state (a(3)A(")) of HNO∕HON. It was found that the reaction is dominated by long-lived resonances supported by the HNO and HON wells. Significant non-reactive scattering was observed, indicating substantial deviations from the statistical limit. Due to the large exothermicity of the reaction, the NO product has hot internal state distributions: its rotational state distribution is inverted and peaks near the highest accessible rotational level; and its vibrational state distribution extends to υ = 10 and decays monotonically with the vibrational quantum number. In particular, the predicted product vibrational distribution is in reasonably good agreement with experiment. The calculated differential cross section is dominated by scattering in both the forward and backward directions, consistent with the formation of reaction intermediates.