乙醇醛光解离机理的理论研究

采用多参考态方法, 在CASPT2//CASSCF/6-311+G(2df, 2p) 水平上计算了乙醇醛(HOCH2CHO)分子在三个最低电子态(S0、S1和T1)上驻点的电子结构和解离势能面。结合势能面交叉点,探讨了HOCH2CHO与波长有关的光解离机理,分析了可能的光解离产物。结果表明, 在实验光解波长240 – 400 nm的激发下,HOCH2CHO分子主要发生S1态上的解离反应或通过S0和S1态之间的振动相互作用驰豫到基态,随之发生基态解离反应。C-C键断裂生成基态光解产物HOCH2 (2A′)+ HCO (2A′)是最主要的反应途径;而在一定波长下,生成CH3OH + CO的基态协同反应、脱醛基氢及脱羟基通道都是能量上可行的反应途径。本文的计算结果和实验观察一致。     

叠氮化氰的光解离机理

采用多参考态方法,在MRCI+Q//CAS(10,9)/6-311+G(2df)水平上对叠氮化氰(N3CN)的光解离机理进行理论研究.优化得到基态(S0)和低激发态(S1、S2、T1)势能面上的极小点、过渡态、内转换交叉点(IC-S1/S0)和隙间窜跃交叉点(ISC-S1/T1)的结构和能量,构建反应势能面.在MRCI+Q//CAS(10,9)水平上计算N3CN的垂直激发能,并和实验值进行对比.结果表明,在S0、S1、S2和T1态势能面上,N—N键断裂生成N2+NCN是主要解离途径,而C—N键断裂通道是次要通道.实验观测到220 nm处的吸收峰对应分子由S0态到S1态的激发,对应主要光解离产物为NCN[a1△g];而在275 nm处的吸收峰则对应分子被激发到T1态,然后直接生成基态产物NCN[X3Σg-].我们的理论结果与实验测量符合得很好.     

HS+HO2气相反应机理及主通道速率常数的理论研究

采用CCSD(T)/6-311++G(3df, 2pd)//B3LYP/6-311+G(2df, 2p)双水平计算方法构建了HO₂+HS反应体系的单、三重态反应势能面,并对该反应主通道的速率常数进行了研究。研究结果表明,标题反应经历了八条反应通道,其中三重态反应通道R1是标题反应主通道。此通道包含路径Path 1 (R → ³IM1 → ³TS1 → P1(³O₂+H₂S))和Path 1a (R → ³IM1a → ³TS1a → P1(³O₂+H₂S))两条路径。利用经典过渡态理论(TST)与变分过渡态理论(CVT)并结合小曲率隧道效应模型(SCT),分别计算了主路径Path 1和Path 1a在200-800 K温度范围内的速率常数k^TST、k^CVT和k^CVT/SCT,在此温度区间内路径Path 1和Path 1a具有负温度系数效应。速率常数计算结果显示,对主路径Path 1和Path 1a而言,变分效应在计算温度段内有一定影响,与此同时量子力学隧道效应在低温段有显著影响。路径Path 1和Path 1a的CVT/SCT速率常数的三参数表达式分别为k₁^CVT/SCT(200-800 K) = 1.54×10^-5T^-2.70exp(1154/T) cm³ ·molecule^-1·s^-1和k_1a^CVT/SCT(200-800 K) = 5.82×10^-8T^-1.84exp(1388/T) cm³·molecule^-1·s^-1。     

N-甲硝胺异构化和分解反应机理和动力学的理论研究

The complex potential energy surface and reaction mechanisms for the unimolecular isomerization and decomposition of methyl-nitramine (CH3NHNO2) were theoretically probed at the QCISD(T)/6-311+G*//B3LYP/6-311+G* level of theory. The results demonstrated that there are four low-lying energy channels: (i) the NN bond fission pathway; (ii) a sequence of isomerization reactions via CH3NN(OH)O; (IS2a); (iii) the HONO elimination pathway; (iv) the isomerization and the dissociation reactions via CH3NHONO (IS3). The rate constants of each initial step (rate-determining step) for these channels were calculated using the canonical transition state theory. The Arrhenius expressions of the channels over the temperature range 298-2000 K are k6(T)=1014:8e-46:0=RT , k7(T)=1013:7e-42:1=RT , k8(T)=1013:6e-51:8=RT and k9(T)=1015:6e-54:3=RT s-1, respectively. The calculated overall rate constants is 6.9￡10-4 at 543 K, which is in good agreement with the experimental data. Based on the analysis of the rate constants, the dominant pathway is the isomerization reaction to form CH3NN(OH)O at low temperatures, while the NN bond fission and the isomerization reaction to produce CH3NHONO are expected to be competitive with the isomerization reaction to form CH3NN(OH)O at high temperatures.     

Ab initio Study on Mechanism of Forming a Silicic Bis-Heterocyclic Compound Between Dimethylmethylenesilylene and Ethene

The mechanism of the cycloaddition reaction of forming a silicic bis-heterocyclic compound between singlet dimethylmethylenesilylene (Me₂C=Si:) and ethene has been investigated with the CCSD(T)//MP2/6-31G^* method. From the potential energy profile, it can be predicted that, this reaction has one dominant channel. The presented rule of this dominant channel: the 3p unoccupied orbital of Si in dimethylmethylenesilylene and the π orbital of ethene forming the π→p donor-acceptor bond, resulting in the formation of three-membered ring intermediate (INT1); INT1 then isomerizes to a four-membered ring silylene (P2), which is driven by ring-enlargement effect; due to sp³ hybridization of Si atom in P2, P2 further combines with ethene to form a silicic bis-heterocyclic compound.     

Theoretical study of the reaction of Ni with OCS

The reactivity of Ni⁺ with OCS on both doublet and quartet potential energy surfaces (PES) has been investigated at the B3LYP/6-311+G(d) level. The object of this investigation was the elucidation of the reaction mechanism. The calculated results indicated that both the CS and CO bond activations proceed via an insertion–elimination mechanism. Intersystem crossing between the doublet and quartet surfaces may occur along both the CS and CO bond activation branches. The ground states of NiS⁺ and NiO⁺ were found to be quartets, whereas NiCO⁺ and NiCS⁺ have doublet ground states. The CS bond activation is energetically much more favorable than the CO bond activation. All theoretical results are in line with early experiments.     

Theoretical study of the reaction of Cu with OCS

The reactivity of Cu⁺ with OCS on both singlet and triplet potential energy surfaces (PES) has been investigated at the UB3LYP/6-311+G(d) level. The object of this investigation was the elucidation of the reaction mechanism. The calculated results indicated that both the C–S and C–O bond activations proceed via an insertion–elimination mechanism. Intersystem crossing between the singlet and triplet surfaces may occur along both the C–S and C–O bond activation branches. The ground states of CuS⁺ and CuO⁺ were found to be triplets, whereas CuCO⁺ and CuCS⁺ have singlet ground states. The C–S bond activation is energetically much more favorable than the C–O bond activation. All theoretical results are in line with early experiments.     

正一溴代烷烃的紫外光解动力学研究

利用共振增强多光子电离飞行时间质谱(REMPI-TOFMS),研究了长链正一溴代烷烃R_Br(R为正烷烃基)(C₂H5Br,n-C₃H₇Br,n-C₄H₉Br)在234及267nm附近的光解动力学.溴碎片来源于R_Br的直接解离:R_Br→R+Br(²P_3/2)/Br^*(2P1/2).根据测定的离子信号强度,得到了Br^*与Br的分支比N(Br^*)/N(Br)及相应的相对量子产额(Br^*)和(Br).(Br^*)与激光波长及分子结构显示了一定的依赖关系.将实验结果用CH3Br的解离模型进行拟合,得到了长链R_Br的光解动力学行为的定性解释.     

Quasi-classical Trajectory Study of F-I-H20→HF-I-OH Reaction： Influence of Barrier Height,Reactant Rotational Excitation,and Isotopic Substitution

The reaction dynamics of the F＋H20/D20→HF/DF＋OH/OD are investigated on an accurate potential energy surface （PES） using a quasi-classical trajectory method. For both isotopomers, 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 excitation of the H20 reactant leads to significant enhancement of the reactivity, suggesting a strong coupling with the reaction coordinate.     

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.     

HOD分子在(~C)1B1态的光解动力学:产物OH/OD的分支比以及OD键解离能

Photodissociation of jet-cooled HOD via the ? state around 124 nm has been studied using the H(D)-atom Rydberg tagging time-of-flight technique. Rotational state resolved action spectrum and the product translational energy distribution spectra have been recorded for both D+OH and H+OD dissociation channels. Product channel OH/OD branching ratios for the individual ?- X rotational transition have been determined. A comparison is also given with the B- X and ?- X transitions. In addition, the dissociation energy of the OD bond in HOD has been determined accurately to be 41751.3±5 cm^-1.     

Density Functional Theory Study of Mechanism of Cycloaddition Reaction Between Dimethyl-Silylene Carbene and Acetone

The mechanism of the cycloaddition reaction between singlet dimethyl-silylene carbene and acetone has been investigated with density functional theory, From the potential energy profile, it can be predicted that the reaction has two competitive dominant reaction pathways. The presented rule of this reaction: the [2+2] cycloaddition effect between the πorbital of dimethyl-silylene carbene and the π orbital of π-bonded compounds leads to the formation of a twisty four-membered ring intermediate and a planar four-membered ring product; The unsaturated property of C atom from carbene in the planar four-membered ring product,resulting in the generation of CH₃-transfer product and silicic bis-heterocyclic compound.     

Ab initio Study of Mechanism of Forming Germanic Bis-Heterocyclic Compound Between Dichloro-Germylene Carbene (Cl2Ge=C:) and Formaldehyde

The mechanism of the cycloaddition reaction of forming germanic bis-heterocyclic compound between singlet dichloro-germylene carbene and formaldehyde has been investigated with CCSD(T)//MP2/6-31G^* method, from the potential energy profile, we predict that the re-action has two competitive dominant reaction pathways. The presented rule of this reaction: the 2p unoccupied orbital of the C atom in dichloro-germylene carbene insert the π orbital of formaldehyde from oxygen side, resulting in the formation of intermediate. In the interme-diate and between two reactants, because of the two bonding π orbital in dichloro-germylene carbene and formaldehyde have occurred [2+2] cycloaddition reaction, forming two four-membered ring compounds in which Ge and O are in the opposite orientation and in the syn-position, respectively. Because of the unsaturated property of C atom from carbene in the two four-membered ring compounds, they further reacts with formaldehyde, resulting in the generation of two germanic bis-heterocyclic compounds.     

BH+2与C2H2反应势能面的量子化学研究

用B3LYP/6-311G(d,p) 密度泛函方法和高级电子相关的CCSD(T)/6-311G(2df,p)偶合簇法研究BH+2与C2H2反应势能面. 结果表明 该势能面上存在(a) H2B+*C2H2, (b) HBCHCH2+, (c) H2BCCH+2和(d) H2*BHCCH+四种异构体, 其中(b)能量最低且在动力学上最稳定, (a), (c)和(d)在动力学上均不稳定; BH+2通过对C2H2的分步亲电加成以及随后的氢迁移和H2消除等反应形成离解产物HBCCH++H2, 该反应不需要活化能且大量放热. 计算结果有助于深入了解BH+2等缺电子硼氢正离子的反应行为.     

气相中Y^＋活化CS2中的C—S键

以Y+与CS2反应作为第二行前过渡金属离子与CS2反应的范例体系.采用密度泛函UB3LYP方法,对于Y+采用Stuttgart赝势基组,对于CS2采用6-311+G(2d)基组,计算研究了Y+离子在基态和激发态时与CS2气相反应的机理.并用UCCSD(T)方法在相同的基组水平上对各驻点作了单点能量校正.结果表明Y+离子与CS2的反应是插入-消去反应,在反应过程中会发生系间窜越,并且找到了两个势能面的能量最低交义点.     

Interpolated potential energy surface for abstraction and exchange reactions of NH 3 + H and deuterated analogues

An ab initio interpolated potential energy surface for the hydrogen abstraction and exchange reactions between ammonia and a hydrogen atom is reported. The interpolation is constructed over a set of data points calculated at the unrestricted coupled cluster approximation, using single and double excitations, and including the triple excitations non-iteratively. New data point selection methods were used to improve the convergence and accuracy of the interpolated surface.     

基本不变量神经网络解析梯度方法的研究

提升势能面的运行速度对于动力学模拟至关重要. 相对于计算简单、 但耗时更长的数值梯度计算, 直接求解势能面梯度的解析公式能够大幅提高势能面的运行效率. 本文发展了基本不变量神经网络解析梯度的生成方法. 计算解析梯度的代码可以通过程序自动生成. 对大量数据点进行测试后, 证明了该方法可以得到正确的势能面梯度输出结果. 通过测试不同势能面的调用时间, 发现采用解析梯度方法能够带来10倍以上的性能提升. 随着体系的增大, 这种性能提升也会越明显.