次甲基硅SiH与H2X(X=O,S)反应的热力学及动力学性质研究

在量子化学对SiH与H₂O和H₂S反应计算的基础上,运用统计热力学和Wigner校正的Eyring过渡态理论,计算了上述两反应在200～2000 K温度范围内的热力学函数、平衡常数、频率因子A和速率常数随温度的变化。计算结果表明,两反应在低温下具有热力学优势,而在高温下具有动力学优势。比较两反应的计算结果发现,在相同的温度下,SiH与H₂O反应比SiH与H₂S反应放热较多,但速率常数却较小。SiH与H₂O反应和前文报道的SiH与HF反应的比较表明,SiH与H₂O反应放热较少,而且在相同温度下,速率常数也较小。     

NH+O3→ONH+O2反应热力学和动力学研究

在量子化学对 NH自由基与臭氧 O3反应计算的基础上,应用统计热力学方法研究了 100～ 1600 K温度范围内 NH和 O3反应过程的各热力学量的变化及平衡常数,用经 Wigner校正的 Eyring过渡态理论计算了不同温度下该反应两不同反应通道的活化热力学量、反应速率常数及频率因子.计算表明,相对于反应通道 II,反应通道 I不仅有很强的反应自发性,而且在动力学上也是较容易实现的反应.     

NH＋O3→ONH＋O2反应热力学和动力学研究

在量子化学对NH自由基与臭氧O3反应计算的基础上，应用统计热力学方法研究了100～1600 K温度范围内NH和O3反应过程的各热力学量的变化及平衡常数，用经Wigner校正的Eyring过渡态理论计算了不同温度下该反应两不同反应通道的活化热力学量、反应速率常数及频率因子.计算表明，相对于反应通道II，反应通道I不仅有很强的反应自发性，而且在动力学上也是较容易实现的反应.     

锆吸氘和氚的动力学同位素效应

应用反应速率分析方法,在高真空金属系统上测定了锆在恒容体系和450～630 ℃范围内吸收氘和氚的p-t曲线。由p-t曲线可知,在同一温度下,锆吸氚的速率低于锆吸氘的速率;而随着温度的升高,由于解吸逆反应的影响,锆吸氚(氘)的反应速率降低。根据p-t曲线计算了吸气反应在不同温度的速率常数,得到锆吸氘和氚的表观活化能分别为(-25.9±0.7)和(-16.8±0.8)kJ·mol^-1。从活化能数据可以看出     

十氢化萘低温燃烧反应的动力学机理

采用量子化学方法研究了十氢化萘低温燃烧的动力学机理,获得了脱氢反应、自由基加氧反应及1,5氢迁移反应等反应的动力学参数,并在CBS-QB3水平下获得了相关物种的热力学参数,通过过渡态理论计算获得了具有紧致过渡态反应的高压极限速率常数,而无能垒反应的速率常数则由变分过渡态理论得到.基于此机理分析了十氢化萘低温反应的动力学规律和热力学机制.相比于链烷烃和单环烷烃,十氢化萘自由基加氧反应的速率常数随温度变化较快,1,5-氢迁移反应的能垒较高,揭示了物质结构对反应动力学的影响.热力学平衡常数分析结果表明,在低温下十氢化萘自由基加氧反应起主导作用.通过拟合获得了所有反应Arrhenius形式的速率常数,这些参数可用于双环烷烃低温燃烧机理的构建和优化.     

单重态CCl~2与O~3反应的热力学及动力学性质研究

在量化计算的基础上运用统计热力学和Wigner校正的Eyring过渡态理论研究了不同温度下单重态CCl~2和臭氧O~3反应的热力学及动力学性质。计算结果表明该反应在低温下具有热力学优势,而在高温下具有动力学优势。     

计算化学在正丁醚制备实验教学中的应用*

正丁醚的制备是重要的大学有机化学实验,为提高学生对该实验所涉及的反应机理和关键操作要点的深入了解,采用Gaussian计算软件对正丁醇在酸催化下和无催化剂下的反应体系进行了研究,重点考察了酸催化下反应的主、副反应方向的反应机理。结果表明无催化剂下,正丁醇在常压液相下几乎不能发生反应;在酸催化下,正丁醇发生取代反应生成醚的反应路径是优势反应通道;酸催化下正丁醇的取代和消除反应速率常数均随温度增加而迅速增大,但消除反应的反应速率随温度增加更快,温度超过420 K消除反应将变得很明显,综合考虑,制备正丁醚的反应温度应控制在130~140 ℃之间较为合适。利用计算化学以图、表和动图等形式直观、动态、量化地解释了正丁醇成醚和成烯反应的竞争,该结果有助于更好地控制该反应体系,可用作实验教材的补充内容。     

OH+C2H2与OH+C2H4加成反应的热力学和动力学

OH基与乙炔(乙烯)的反应是控制大气中OH基浓度的重要化学反应。对反应OH+C_2H_2,1975年Davis等用FP-RF技术测定了反应的速率常数。1977年,Perry又用同样方法对该反应作出了研究;结果表明,此反应的速率常数强烈地依赖于压力,与Davis等人的实验结果不符。对反应OH+C_2H_4,Atkinson等人的实验研究表明,压力在30.0kPa以下,反应的速率常数随压力而改变;而在30.0—88.4 kPa之间,与压力无关.这与前人的结果不同。两个反应的产物也因温度不同而异.此外,上述反应一般是在近101kPa下发生的,而实验结果大多是在低压下获得的。且目前尚无直接的理论计算结果。为此,我们在从头算水平上用统计热力学方法,对反应     

反应N+NH→N2+H的态-态量子动力学研究

采用在MRCI/aug-cc-pVQZ水平上构建的N₂H基态势能面, 并运用Chebyshev实波包法研究了N + NH→N₂ + H反应的量子动力学, 如反应几率、 积分截面以及产物振转态分布等. 在50~500 K温度范围内, 该反应的速率常数随着温度升高而递增, 与基于其它势能面的理论结果吻合. 然而, 在室温条件下, 所有理论计算的速率常数均显著大于实验值.     

金属卟啉的合成及其对细胞色素P-450的模拟(ⅪⅤ)——钴卟啉对环己烷的单充氧催化作用

报道了钴卟啉催化下PhIO对环己烷的单充氧反应及其反应动力学行为。探讨了反应温度和有机配体对该反应的影响。研究结果表明:在配体吡啶的作用下,钴卟啉对PhIO单充氧化环己烷具有催化活性,其催化性能与反应温度以及吡啶用量有关。反应速率随温度增加而增加,且表观反应速率常数与反应温度之间存在Arrhenius关系。反应产率随温度升高而降低。当吡啶与TPPCo(Ⅱ)用量为50:1时,TPPCo(Ⅱ)催化效果最好。     

Thermodynamic and Kinetic Studies on the SiH + XH_3 (X=N, P) Reactions

Silicon chemistry has been attracted more attention because of its applications to the production of thin silicon films and the etching of silicon wafers in micro-electronics1. Silylidyne (SiH), which plays an important role in plasma chemical vapor deposition (CVD) processes, has been investigated in experimental research2-5. The reaction mechanism of SiH insertion reaction with small molecules such as XH3 (X=N, P) was recently studied by means of M豯ler-Plesset perturbation theory6.…     

ANTA的结构、性质及其互变异构的理论研究

对3-硝基-5-氨基-1,2,4-三唑(ANTA)的三种异构体,1H-ANTA(Ⅰ),2－ANTA(Ⅱ)和4H-ANTA(Ⅲ)在,bainitio-HF/3-21G和DFT-B3LYP/3-21G势能面计算的基础上,进行6-311G^**几何参数全优化,MP2总能量和SCRF溶剂(四氢呋喃)效应计算。以振动分析和统计热力学为基础,作标题物热力学性质以及Ⅰ和Ⅱ之间的互变异构反应计算,求得分子几何,电子结构和300～1000K范围的焓、熵和热容以及Ⅰ和Ⅱ互变异构平衡常数和速率常数。发现在三种异构体中在通常温度下以Ⅱ在气相下最稳定,Ⅰ在溶液中最稳定。低温下难以发生异构化反应,温度可提高Ⅰ与Ⅱ之间的互变速率,在800K时两种异构体在气相中等量共存;大于800K时Ⅰ更为稳定。     

硅乙烯异构化反应的理论研究——从头算水平上的热力学、动力学分析

据报道,含硅碳双键的化合物已能制备,但它在通常条件下不稳定,易于发生加成、异构化反应,研究最多的两个异构化反应是:     

Theoretical study of reactions between AlH(1Σ) and HF molecule

Reaction mechanisms between AlH (¹Σ) and HF molecule are theoretically investigated. Ab initio calculations demonstrate that there are two parallel reaction channels: one is an addition reaction to give H₂AlF via the three‐membered ring transition state (TS) and the other is a dehydrogenation reaction to give AlF+H₂ via the four‐membered ring TS. The addition reaction is thermodynamically favorable and the dehydrogenation reaction is kinetically favorable. Thermodynamics and Eyring transition state theory (TST) with the Wigner correction are also used to compute the thermodynamic functions, the equilibrium constants, A factors, and the rate constants of these reaction channels at 200–1000 K. From the thermodynamics and TST calculations, it is valuable to point out that consideration on the entropy and thermal enthalpy is quite important in the study of chemical reactions on the basis of ab initio method. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 73: 417–424, 1999     

Theoretical studies of the reaction of hydroxyl radical with methyl acetate

The mechanisms and the kinetics of the OH (OD) radicals with methyl acetate CH3C(O)OCH3 are investigated theoretically. The dual-level direct dynamics method is employed in the calculation of the rate constants. The optimized geometries and frequencies and the gradients of the stationary points are calculated at the MP2/6-311G(d,p) level. The energetic information of potential energy surfaces is further refined by the multicoefficient correlation method based on QCISD (MC-QCISD) using the MP2/6-311G(d,p) geometries. Four channels are found for the title reaction. The calculated results reveal that there exists an attractive well (reactant complex) in each entrance H-abstraction channel, that is, the H-abstraction reaction makes a stepwise mechanism. The rate constants are calculated by the canonical variational transition-state theory (CVT) with the interpolated single-point energies (ISPE) approach in the temperature range of 200-1200 K. The small-curvature tunneling effect (SCT) approximation is used to evaluate the transmission coefficient. The calculated rate constants are in good agreement with the experimental ones in the measured temperature range. It is shown that the "out-of-plane hydrogen abstraction" from the methoxy end is the dominant channel at the lower temperatures, and the other two H-abstraction channels should be taken into account with the temperatures increasing. The kinetic isotope effects (KIEs) for the three H-abstraction channels and the total reaction are "inverse", and these theoretically calculated KIEs as a function of temperature are expected to be useful for the future laboratory investigation.     

Direct dynamics studies on hydrogen abstraction reactions of CH3CHFCH3 and CH3CH2CH2F with OH radicals

The hydrogen abstraction reactions of CH3CHFCH3 and CH3CH2CH2F with the OH radicals have been studied theoretically by a dual-level direct dynamics method. The geometries and frequencies of all the stationary points are optimized by means of the DFT calculation. There are complexes at the reactant side or exit route, indicating these reactions may proceed via indirect mechanisms. To improve the reaction enthalpy and potential barrier of each reaction channel, the single point energy calculation is performed by the MC-QCISD/3 method. The rate constants are evaluated by canonical variational transition state theory (CVT) with the small-curvature tunneling correction method (SCT) over a wide temperature range 200-2000 K. The canculated CVT/SCT rate constants are consistent with available experimental data. The results show that both the variation effect and the SCT contribution play an important role in the calculation of the rate constants. For reactions CH3CHFCH3 and CH3CH2CH2F with OH radicals, the channels of H-abstraction from -CHF- and -CH2- groups are the major reaction channels, respectively, at lower temperature. Furthermore, to further reveal the thermodynamics properties, the enthalpies of formation of reactants CH3CHFCH3, CH3CH2CH2F, and the product radicals CH3CFCH3, CH3CHFCH2, CH3CH2CHF, CH3CHCH2F, and CH2CH2CH2F are studied using isodesmic reactions.     

Theoretical Studies on the Abstraction Reaction of Atomic O（^3p） with Si2H6

The hydrogen abstraction reaction of O(^3P) with Si2H6 has been studied theoretially. Two transition states of ^3A″ and ^3A′ symmetries have been located for this abstraction reaction. Geometries have been optimized at the UMP2 leve with 6-311G (d) basis set. G3MP2 has been used for the final single-point energy calculation. The rate constants have been calculated over a wide temperature range of 200-3000K using canonical variational transition-state sheory (CVT) with small curvature tunneling effect(SCT). The calculated CVT/SCT rate constants match well with the experimental value.     

Unimolecular elimination of HF and HCl from chemically activated CF3CFClCH2Cl

The unimolecular reactions of CF3CFClCH2Cl molecules formed with 87 kcal mol(-1) of vibrational energy by recombination of CF3CFCl and CH2Cl radicals at room temperature have been characterized by the chemical activation technique. The 2,3-ClH and 2,3-FH elimination reactions, which have rate constants of (2.5 +/- 0.8) x 10(4) and (0.38 +/- 0.11) x 10(4) s(-1), respectively, are the major reactions. The 2,3-FCl interchange reaction was not observed. The trans (or E)-isomers of CF3CFCHCl and CF3CClCHCl are favored over the cis (or Z)-isomers. Density functional theory at the B3PW91/6-31G(d',p') level was used to evaluate thermochemistry and structures of the molecule and transition states. This information was used to calculate statistical rate constants. Matching the calculated to the experimental rate constants for the trans-isomers gave threshold energies of 62 and 63 kcal mol(-1) for HCl and HF elimination, respectively. The threshold energy for FCl interchange must be 3-4 kcal mol(-1) higher than for HF elimination. The results for CF3CFClCH2Cl are compared to those from CF3CFClCH3; the remarkable reduction in rate constants for HCl and HF elimination upon substitution of one Cl atom for one H atom is a consequence of both a lower E and higher threshold energies for CF3CFClCH2Cl.     

Quantum statistical study of O + O2 isotopic exchange reactions: cross sections and rate constants

Using a wave packet based statistical model, we compute cross sections and thermal rate constants for various isotopic variants of the O + O2 exchange reaction on a recently modified ab initio potential energy surface. The calculation predicts a highly excited rotational distribution and relatively cold vibrational distribution for the diatomic product. A small but important threshold effect was identified for the (16)O + 18O2 reaction, which is suggested to contribute to the experimentally observed negative temperature dependence of the rate ratio, k(18O + 16O2)/k(16O + 18O2). Despite reasonable agreement with quasiclassical trajectory results, however, the calculated thermal rate constants are smaller than experimental measurements by a factor from 2 to 5. The experimentally observed negative temperature dependence of the rate constants is not reproduced. Possible reasons for the theory-experiment discrepancies are discussed.