O(1D)和CD4反应的准经典轨线

在新的全域势能面上, 用准经典轨线方法细致地研究了O(¹D)+CD₄多通道化学反应的动力学．这个势能面是用交换不变多项式方法基于MRC+Q/aug-cc-pVTZ从头算点拟合得到的．通过计算得到了产物OD+CD₃、D+CD₂OD/CD₃O和D₂+DCOD/D₂CO的分支比、平动能分布以及角度分布,结果显示理论与实验吻合得较好, 从而说明了这个反应的同位素取代效应很小． 研究表明,O(¹D)+CD₄反应是经过陷入的抽取机理发生的： 最初主要通过D原子的抽取,并不是之前人们认为的直接C-D键的插入形成CD3OD中间物后再进而解离成各个产物通道．     

Termolecular chemistry facilitated by radical-radical recombinations and its impact on flame speed predictions

Recent theoretical studies have shown that termolecular chemistry can be facilitated through reactions of flame radicals (H, O, and OH) or O₂ with highly-energized collision complexes (either radical or stable species) formed in exothermic reactions. In this work, radical-radical recombination reaction induced termolecular chemistry and its impact on combustion modeling was studied. Two recombination reactions, H + CH₃ + M → CH₄ + M and H + OH + M → H₂O + M, were analyzed using ab-initio master equation analyses guided by quasiclassical trajectory results. The dynamics results and the master equation calculations indicate that CH₄^? and H₂O^? (formed in the two radical-radical reactions outlined above) react rapidly with flame radicals and O₂ at rates that are competitive with collisional cooling. The addition of these processes into conventional combustion modeling requires two modifications: the inclusion of the new nonthermal termolecular reaction rates and the simultaneous reduction of the competing recombination reaction rates. The former is described with newly derived Arrhenius expressions based on quasiclassical trajectories, and the latter is achieved by perturbing the recombination reaction rate during the simulation. Kinetic modeling was used to gauge the impact of including this nonthermal chemistry for H₂/CH₄-air laminar flames speeds. Inclusion of this nonthermal chemistry has a noticeable impact on simulated flame speeds. The procedure developed here can be utilized to properly quantify the effects of such nonthermal reactions in macroscopic kinetic models.     

Quasiclassical dynamics of light + heavy-heavy atom reactions: the reaction X + F2 → XF + F (X = Mu,H, D,T)

Quasiclassical trajectory calculations have been carried out for the collinear reaction X + F₂(v = 0, 1) → XF + F (X = Mu, H, D, T) using the best extended LEPS surface No II of Jonathan et al.. The behaviour of individual trajectories has been used to find the regions of the potential surface that are important for determining the total reaction probability. The bottleneck close to the saddle point and repulsive wall of the exit valley are found to be the important dynamical regions of the surface. Comparison is made with other quasiclassical calculations for light + heavy-heavy atom reactions and with quantum calculations.     

Stereodynamics of chemical reactions: quasi-classical,quantum and mixed quantum-classical theories

In this review, some benchmark works by Han and coworkers on the stereodynamics of typical chemical reactions, triatomic reactions H + D₂, Cl + H₂ and O + H₂ and polyatomic reaction Cl+CH₄/CD₄, are presented by using the quasi-classical, quantum and mixed quantum-classical methods. The product alignment and orientation in these A+BC model reactions are discussed in detail. We have also compared our theoretical results with experimental measurements and demonstrated that our theoretical results are in good agreement with the experimental results. Quasi-classical trajectory (QCT) method ignores some quantum effects like the tunneling effect and zero-point energy. The quantum method will be very time-consuming. Moreover, the mixed quantum-classical method can take into account some quantum effects and hence is expected to be applicable to large systems and widely used in chemical stereodynamics studies.     

D+CD4→CD3+D2反应的四维量子散射计算

运用约化维数量子动力学理论，利用含时波包法，对反应D+CD₄→CD_{3+D2进行了四维量子散射计算.将反应多原子CD4看作双原子D—CD 3，反应D+CD4→CD3+D2看作单原子-双原 子反应，把体系的反应简化为四维散射问题. 波函数的传播采用分裂算符法，为避免格点边界处含时波函数的边界反射，采用了光学吸收 势法，在格点边界处引入光学势，消除边界反射.根据CD4分子的C3v对称性， 选取了Jordan和Gilbert提出的半经验势能面.计算结果表明，反应概率随平动能的变化图像 ，呈现出显著的量子共振特性，这是很多提取反应的共同特征.而不同振动态下的反应概率 随平动能的变化表明，随振动量子数的增大，反应概率有明显提高，且反应阈能明显降低， 这说明反应分子的振动能对分子的碰撞反应有重要贡献.而对基态和第一振动激发态时散射 截面的计算，也证明了这一结论.同时，还分别通过计算量子数j，k，m对反应概率的影响， 对该反应的空间取向效应进行了研究，并与H+CH4→CH3+H2反应进行了比较. 关键词： 含时波包 量子散射 反应概率}     

Roaming

In this review the roaming pathway in chemical reaction dynamics is described and illustrated. After a brief review and illustration of this pathway, an update on the chemical reactions where this pathway has been reported is given. Then the focus is on two recent systems where new aspects of roaming have been reported, namely roaming-isomerisation in CH₃NO₂ unimolecular dissociation and multi-electronic-state roaming in NO₃ photodissociation. Full-dimensional potentials for these systems have been developed and used in quasiclassical trajectory calculations. Selected results from these studies are reviewed. So-called roaming saddle points and associated intrinsic reaction coordinate paths have been determined for these systems. These will be reviewed and investigated by making relevant comparisons with quasiclassical trajectories. Manifestations of roaming in quantum scattering and resonance calculations of the MgH + H → Mg + H₂ reaction are reviewed. That is followed by a review of several statistical theories of roaming. The review concludes with some ‘editorial’ comments about the challenges that roaming presents to approximate theories and comments about possible future research directions.     

Analyses of the semi-classical wavepacket approach to chemical reactions: the F + H2 → HF + H reaction

Mixed quantum-classical calculations using propagation of 3D wavepackets in hyperspherical coordinates have been performed for the F + H₂ → HF + H reaction using two different potential energy surfaces. Vibrationally resolved cross-sections show good agreement with those obtained from accurate quantum mechanical calculations.     

氘代最简单Criegee中间体CD2OO的振转光谱研究

本文采用振动自洽场/虚组态相关(VSCF/VCI)方法计算了氘代最简单Criegee中间体CD₂OO的振动和转动光谱. 计算得到的基频振动频率和转动常数与已有实验吻合. 这些数据可以用于未来该体系的光谱研究,特别是振动激发的转动常数对于实验光谱指认非常重要. 另外,不同来源的光谱强度,包括本文理论计算,文献中在NEVPT2和B3LYP水平上的计算结果以及实验结果,互相之间均不符合.     

Reaction mechanism and rate constants of the CH+CH4 reaction: a theoretical study

Ab initio and density functional calculations have been performed to elucidate the mechanism of CH radical insertion into methane. The results show that the reaction can be viewed to occur via two stages. On the first stage, the CH radical approaches methane without large structural changes to acquire proper positioning for the subsequent stage, where H-migration occurs from CH₄ to CH, along with a C–C bond formation. Where the first stage ends and the second begins, a tight transition state was located using the B3LYP/6-311G(d,p) and MP4(SDQ)/6-311++G(d,p) methods. Using a rigid rotor – harmonic oscillator approach within transition state theory, we show that at the MP5/6-311++G(d,p)//MP4(SDQ)/6-311++G(d,p) level the calculated rate constants are in a reasonably good agreement with experiment in a broad temperature range of 145–581 K. Even at low temperatures, the insertion reaction bottleneck is found about the location of the tight transition state, rather than at long separations between the CH and CH₄ reactants. In addition, high level CCSD(T)-F12/CBS calculations of the remainder of the C₂H₅ potential energy surface predict the CH+CH₄ reaction to proceed via the initial insertion step to the ethyl radical which then can emit a hydrogen atom to form highly exothermic C₂H₄+H products.     

The vibrational levels of methane obtained from analyses of high-resolution spectra

Methane and its tetrahedral isotopologues are spherical-top molecules whose high-resolution rovibrational spectra can only be analyzed in detail, thanks to sophisticated symmetry-adapted tensorial models. However, the effective Hamiltonian parameters of such models do not give direct access to the positions of the vibrational sublevels. In this paper, we present a calculation of the vibrational level positions for ¹²CH₄, ¹³CH₄, ¹²CD₄ and ¹³CD₄ performed using the effective Hamiltonian parameters obtained through recent analyses. We also include the results of a re-analysis of the octad system of ¹²CH₄ performed with a higher order of the development which slightly improves the previous work on this polyad.     

Quasiclassical calculation of the chemical reaction Ba+CH3I

This paper reports the quasiclassical trajectory (QCT) calculations on extended London–Eyring–Polanyi–Sato (LEPS) potential energy surface (PES) for the reaction of Ba atom with CH₃I. In the calculations, the Sato parameters needed in construction of the reaction PES are obtained from simulation of experimental results of BaI product vibrational distributions. The parameters obtained are further tested by comparison of the calculated BaI vibrational distributions with another research group's experimental results at different collision energy. In addition, the reaction cross-section versus collision energy, product BaI angular distribution and alignment are obtained by QCT calculations on the PES. The possible reaction configurations are also discussed upon the calculations.     

Calculation of accurate imaginary frequencies and tunnelling coefficients for hydrogen abstraction reactions using IRCmax

The effect of level of theory on the imaginary frequency and corresponding tunnelling coefficients has been studied for a test set of hydrogen abstraction reactions: ?CH₂X + CH₃Y → CH₃X + ?CH₂Y for (X,Y) = (H,H), (H,CN), (H,F), (H,Li) and (F,Li). It is found that the imaginary frequency is very sensitive to the level of theory used, with Hartree-Fock (HF) methods severely overestimating the imaginary frequency compared with high-level CCSD(T)/6-311G(d,p) calculations. The errors for the other methods are smaller but nonetheless significant, with MP2 methods overestimating the imaginary frequency and density functional theory (DFT) methods underestimating it. In the case of the HF methods, this leads to errors in the tunnelling coefficient of several orders of magnitude, while for the better DFT and MP2 methods errors of a factor of 2–3 are observed. To address this problem, an IRCmax procedure for estimating the imaginary frequency has been developed and it is found that IRCmax imaginary frequencies calculated with CCSD(T)/6-311G(d,p) single points along a low-level HF/6-31G(d) minimum energy path provide excellent approximations to the high-level values, at a fraction of the computational cost.     

O~+ + H_2及其同位素取代反应的立体动力学研究

运用准经典轨线方法,基于RODRIGO势能面,对碰撞能为20kcal/mol时,O++H2及其同位素取代反应的立体动力学性质进行了理论研究,对k-j′两矢量相关和k-k′-j′三矢量相关的分布函数、极化微分反应截面,以及产物转动取向参数进行了详细的讨论.结果表明,O++H2→OH++H,O++DH→OD++H和O++TH→OT++H反应的立体动力学性质对体系的质量因数非常敏感.     

F和CH3OH多通道反应的速率常数与分支比的研究

F和CH₃OH有两个夺氢反应通道,分别生成HF+CH₃O 和HF+CH₂OH. 尽管这两个通道都没有能垒,但前一个通道即生成HF+CH₃O的反应分支比远远高于期望的统计平均值(四分之一). 不同实验测得的分支比不仅相去甚远,而且定量上与早期由过渡态理论(稳定点信息在MP2以及G2理论水平下计算得到)得到的计算结果也不符合. 此前在CCSD(T)-F12a/AVDZ水平上计算得到了121000个几何构型的能量,采用对易不变多项式结合神经网络的方法拟合得到了该体系的全维高精度势能面. 本文采用该势能面,结合准经典轨线动力学方法,对该反应的反应速率常数和反应分支比进行了理论研究,得到的结果与实验吻合. 由于反应没有能垒,理论计算结果表明反应速率常数随温度升高而有微弱的下降趋势.     

碰撞能和同位素取代对H+BrF→HBr+F反应立体动力学影响的理论研究

基于拟合得到的London-Eyring-Polanyi-Sato势能面,运用准经典轨线方法对H+BrF→HBr+F反应的立体动力学性质进行了理论研究.计算了反映矢量相关的角分布和对光诱导的双分子反应实验敏感的四个极化微分散射截面.结果表明:随着碰撞能的增加,产物的转动极化变强,并且产物的后向散射占主导地位.通过比较D+BrF→DBr+F和H+BrF→HBr+F反应的产物极化,揭示了明显的同位素效应.     

Anharmonic force field for methanol

An ab initio quartic anharmonic force field for methanol has been calculated at the equilibrium position using the CCSD(T) method for the structure and the harmonic potential energy surface, and the MP4(SDQ) method for the anharmonic part of the surface. A triple zeta basis set was employed with symmetrized curvilinear internal valence coordinates in all calculations. The internal coordinate force field constants have been transformed into force constants in the dimensionless normal coordinate representation for various isotopomers. Vibrational term values for CH₃OH, CH₃OD, CD₃OH, and CD₃OD have been obtained using second order perturbation theory. Particular care has been devoted to the inclusion of Fermi resonance interactions between different vibrational states. A good accuracy has been achieved in the calculation of the fundamentals for all the isotopomers, the mean absolute error being 5.8 cm^?1.     

From quantum chemistry to dissociation kinetics: what we need to know †

The relationship between rate constants for dissociation and the reverse association reactions and their potential energy surfaces is illustrated. The reaction systems e^? + SF₆ ? SF₆ ^? →SF₅ ^? + F, H + CH₃ ?CH₄, 2 CF₂ ? C₂F₄, H + O₂ →HO₂, HO + O ?HO₂ ? H + O₂, and C + HO →CHO are chosen as representative examples. The necessity to know precise thermochemical data is emphasised. The interplay between attractive and anisotropic components of the potentials influences the rate constants. Spin–orbit and electronic–rotational coupling in reactions between electronic open-shell radicals so far generally has been neglected, but is shown to have a marked influence on low temperature rate constants.     

Quantum mechanical computational studies of chemical reactions : III. Collinear A + BC reaction with some model potential energy surfaces

The effects of the location of the energy barrier and of changing vibrational frequency along the reaction path on the reaction dynamics of a collinear A + BC → AB + C reaction were studied using a series of ‘diagnostic’ model potential energy surfaces. In accord with the classical trajectory study, vibration was found markedly more effective than translation in promoting reaction on a surface with a barrier located at the exit valley, and vice versa on a surface with a barrier placed in the entry valley. Considerable vibrational population inversion in the reaction product was found for an ‘early-down-hill’-type surface with a barrier placed in the entry valley. The energy release was exceedingly speecific, particularly so for the F + H₂/F + D₂ collinear collision on a realistic surface. Resonances in reactive molecular collisions were evident whenever the static effect of changing vibrational frequency along the reaction path gave rise to the potential wells necessary to support these quasibound states.     

A global potential energy surface for H2S+(X 4A′′) and quasi-classical trajectory study of the S+(4S) + H2(X1Σ+g) reaction

A novel global potential energy surface for H₂S⁺(X?⁴A″) based on accurate ab initio calculations is presented. Energies are calculated at the multi-reference configuration interaction level with Davidson correction using aug-cc-pVQZ basis set plus core-polarisation high-exponent d functions. A grid of 4552 points is used for the least-square fitting procedure in the frame of a many-body expansion. The topographical features of the new potential energy surface are here discussed in detail. Such a surface is then employed for dynamic studies of the S(⁴S) + H₂(X?¹Σ⁺_g) →SH⁺(X?³Σ^?) + H(²S) reaction using the quasi-classical trajectory method. State specific trajectories are calculated, for both ground and ro-vibrationally excited initial states of H₂(X?¹Σ⁺_g). Corrections to the zero point energy leakage of the classical calculations are also presented. Calculated reaction cross sections and rate constants are here reported and compared with available literature.