已烯自由基阳离子的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-31G(d,p),6-31+G(d,p),6-311G(d,p)及6- 311+G(d,p)基组,分别对1-C_6H_(12)~+,2-C_6H_(12)~+和3-C_6H_(12)~+的各种构 象进行了几何构型优化,并在B3LYP/6-311G(d,p)水平上进行了频率分析计算,在 各优化构型上,使用B3LYP和MP2(full)方法进行了超精细结构的计算。计算的3- C_6H_(12)~+的超精细偶合常数比以往的计算结果更好;1-C_6H_(12)~+和2-C_6H_ (12)~+的超精细偶合常数目前尚无实验数据报道,本计算预言了它们的超精细偶合 常数和最稳定构型。     

已烯自由基阳离子的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-31G(d,p),6-31+G(d,p),6-311G(d,p)及6- 311+G(d,p)基组，分别对1-C_6H_(12)~+,2-C_6H_(12)~+和3-C_6H_(12)~+的各种构 象进行了几何构型优化，并在B3LYP/6-311G(d,p)水平上进行了频率分析计算，在 各优化构型上，使用B3LYP和MP2(full)方法进行了超精细结构的计算。计算的3- C_6H_(12)~+的超精细偶合常数比以往的计算结果更好；1-C_6H_(12)~+和2-C_6H_ (12)~+的超精细偶合常数目前尚无实验数据报道，本计算预言了它们的超精细偶合 常数和最稳定构型。     

丁烯和丁烷自由基阳离子的分子结构和超精细结构的 理论研究

用密度函数理论B3LYP方法和6-31G(d,p),6-311G(d,p)及6-311+G(d,p)基组，分别对1-C4H^+~8,2-C4H^+~8和C4H^+~10进行了构型优化和频率分析计算，预言1-C4H^+~8具有非平面构型，与以往报道的从头算和密度函数理论计算结果不同。在各自由基阳离子的B3LYP构型上，进行了B3LYP、MP2及MRSDCI方法的超精细偶合常数计算，得到了比以往更好的结果，特别是MP2/B3LYP计算值是至今与实验值符合得最好的理论计算结果。     

CHF_3与O(~3P)反应机理的理论研究

采用从头算MP2/6-311G(d,p)方法研究了CHF3与O(3 P)的反应机理,优化了所有反应物、产物和过渡态的几何构型,并通过振动频率分析和内禀反应坐标(IRC)方法确证了过渡态的真实性.在QCISD(T)/6-311++G(d,p)水平上精确计算了各反应物种的单点能.结果表明,标题反应共存在4类反应,分别为抽提氢反应(R1)、抽提氟反应(R2),消氟化氢反应(R3)和消氢反应(R4),在QCISD(T)/6-311++G(d,p)//MP2/6-311G(d,p)水平上,R1,R2,R3和R4反应的能垒分别为70.7,378.7,294.7和307.2kJ·mol-1,抽提氢反应为主反应通道.     

大气中OH引发的有机磷农药MTMPDA降解机理的研究

采用量子化学方法在ωB97XD/6-311++G(3df,2pd)//ωB97XD/6-31+G(d,p)和M062X/6-311++G(3df,2pd)//M062X/6-31+G(d,p)水平下对大气中MTMPDA在OH自由基引发下的反应机理及其产物的后继反应进行理论研究.结果表明MTMPDA+OH反应主要通过OH抽提与N原子相连的—CH3上的H原子、以及抽提与O原子相连的—CH3上的H原子的方式进行.此外,抽提与N原子相连的—CH3上的H原子生成的产物与大气中的O2和NO发生进一步反应并没有导致MTMPDA毒性明显地降低;而抽提与O原子相连的—CH3上的H原子的产物发生一系列后继反应可以降低MTMPDA的毒性.     

SO_2F_2与H自由基反应机理研究

用密度泛函B3PW91/6-311++g(d,p)方法研究了SO_2F_2与H自由基的反应机理,对反应势能面上的各驻点的几何构型进行了全参数优化,过渡态的真实性通过振动分析和内禀反应坐标(IRC)分析的结果得到了证实。用双水平CCSD(T)/6-311++g(d,p)//B3PW91/6-311++g(d,p)获得了SO_2F_2+H反应势能剖面图,并用传统过渡态理论计算了最佳反应通道的决速步的反应速率常数。结果表明,SO_2F_2+H可经过加成-消去或直接抽提机理生成两种裂解产物P1(SO_2F+HF)和P2(SF_2O+OH),其中通道SO_2F_2+H→TS1→IM1→TS4→P1和SO_2F_2+H→TS3→CP1→P1决速步的势垒分别为89.5和91.3 kJ·mol~(-1),势垒相差仅1.8 kJ·mol~(-1),为标题反应的两个竞争反应,且高温有利于直接抽提反应的发生。SO_2F_2分子中的O原子不可能通过抽提反应直接提取。     

戊烯自由基阳离子的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-31G(d,p)、6-31＋G(d,p)、6-311G(d,p)及6-311＋G(d,p)基组，分别对2-C5H10＋和1-C5H10＋的各种构象进行了几何构型优化，并用B3LYP/6-311G(d,p)进行了频率分析计算.计算预言1-C5H10＋具有非平面构型，与以往报导的从头算计算结论相反.在两个自由基阳离子的各种构象的B3LYP几何构型上，进行了B3LYP和UMP2(full)方法的超精细偶合常数计算，得到了比以往更好的结果.     

4,4′-联吡啶二面角和扭转势垒的理论研究

采用HF、MP2方法和密度泛函理论(包括B3LYP,BLYP和mPW1PW91)方法,分别在6 31G(d),6 311G(d),6 311+G(d),6 311++G(d,p),cc pvdz和cc pvtz水平上,计算了4,4′联吡啶分子的平衡几何构型和扭转势能面以及相应势垒,考察了各种计算方法和基组对于理论预测结果的影响.在π共轭效应和环间相邻原子的空间位阻共同作用下,2个吡啶环间有大约37°的夹角.计算得到的分子结构参数与电子衍射实验结果吻合.较小的扭转势垒说明分子可以很容易进行内转动.     

CH3F与C2H3反应机理的理论研究

采用密度泛函B3LYP/6-311G(d,p)方法对CH3F与C2H3的反应体系进行了理论研究,获得了反应的势能面信息及可能的微观机理.在QCISD(T)/6-311++G(d,p)水平上精确计算了各反应物种的单点能.结果表明,除抽提氢反应外,标题反应还存在抽提氟(R1)、消氟化氢(R2)、消氢(R3)和自由基形成(R4)四类反应.在QCISD(T)/6-311++G(d,p)//B3LYP/6-311G(d,p)水平上,R1,R2,R3和R4反应的能垒分别是163.9,152.2,209.8和224.2kJ·mol-1,相应反应能为-56.6,-164.3,-2.7和-156.0kJ·mol-1,所有反应均放热,为热力学允许的反应.     

OH-与CH2CIF反应的阴离子产物通道

理论研究了羟基负离子(OH-)与氟氯代甲烷(CH2CIF)反应的阴离子产物通道.分别在B3LYP/6-31+G(d,p)和B3LYP/6-311++G(2d,p)水平上得到反应势能面上各关键物种的优化构型,进而计算得到谐振频率和零点能.基于CCSD(T)/6-311+G(3df,3dp)水平的相对能量,描述了由质子转移和双分子亲核取代(SN2)过程生成各阴离子产物的途径.各阴离子产物途径势垒的计算结果表明质子转移过程是实验中的主要产物通道,与以往实验测量的结论相符.此外,计算还显示双分子亲核取代过程得到了非典型的阴离子产物,其中动力学效应可能会导致F-的生成.     

C6F5X＋（X=Cl,Br,I,CH3）离子的理论研究

用DFT B3LYP方法及6-311G（d,p）,6-311＋G（d,p）和LanL2dz基组,对C6F5X＋（X=Cl,Br,I,CH3）阳离子做了理论研究,优化了它们的电子基态的构型,计算了对应分子的垂直电离势（VIP）和绝热电离势（AIP）.结果表明四种离子的构型的对称点群和对应分子相同,但构型参数有明显差别.B...     

HCHO+X(X=F、Cl、Br)的反应机理

采用CCSD/6-311++G(d,p)//B3LYP/6-311++G(d,p)方法研究了HCHO与卤素原子X(X=F、Cl、Br)的反应机理. 计算结果表明, 卤素原子X(X=F、Cl、Br)主要通过直接提取HCHO中的H原子生成HCO+HX(X=F、Cl、Br). 另外还可以生成稳定的中间体, 中间体再通过卤原子夺氢和氢原子直接解离两个反应通道分别生成HCO+HX(X=F、Cl、Br)和H+XCHO(X=F、Cl、Br). 其中卤原子夺氢通道为主反应通道, HCO和HX(X=F、Cl、Br)为主要的反应产物; 且三个反应的活化能均较低, 说明此类反应很容易进行, 计算结果与实验结果符合很好. 电子密度拓扑分析显示, 在HCHO+X反应通道(b)中出现了T型结构过渡态, 结构过渡态(STS)位于能量过渡态(ETS)之后. 并且按F、Cl、Br的顺序, 结构过渡态出现得越来越晚.     

Transition state structure and energetics of the N(2)O + X (X = Cl,Br) reactions

The structural and vibrational properties of the transition state of the N(2)O + X (X = Cl,Br) reactions have been characterized by ab initio methods using density functional theory. We have employed Becke's hybrid functional (B3LYP), and transition state optimizations were performed with 6-31G(d), 6-311G(2d,2p), 6-311+G(3d,2p), and 6-311+G(3df,2p) basis sets. For the chlorine atom reaction the coupled-cluster method (CCSD(T)) with 6-31G(d) basis set was also used. All calculations resulted in transition state structures with a planar cis arrangement of atoms for both reactions. The geometrical parameters of transition states at B3LYP are very similar, and the reaction coordinates involve mainly the breaking of the N-O bond. At CCSD(T)/6-31G(d) level a contribution of the O-Cl forming bond is also observed in the reaction coordinate. In addition, several highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction energetics. All model chemistries predict exothermic reactions. The G3 and G2 methods result in the smallest deviations from experiment, 1.8 and 0 kcal mol(-1), for the enthalpies of reaction for N(2)O reaction with chlorine and bromine, respectively. The G3//B3LYP and G1 methods perform best among the composite methods in predicting energies of the transition state, with a deviation of 1.9 and 3.0 kcal mol(-1), respectively, in the activation energies for the above processes. However, the B3LYP/6-311+G(3df,2p) method gives smaller deviations of 0.4 and -1.0 kcal mol(-1), respectively. The performance of the methodologies applied in predicting transition state energies was analyzed.     

The C-H and alpha(C-X) bond dissociation enthalpies of toluene, C6H5-CH2X (X = F, Cl), and their substituted derivatives: a DFT study

The homolytic C-H bond dissociation enthalpies (BDEs) of toluene and its para- and meta-substituted derivatives have been estimated by using the (RO)B3LYP/6-311++G(2df,2p)//(U)B3LYP/6-311G(d,p) procedure. The performance of two other hybrid functionals of DFT, namely, B3PWP91 and O3LYP, has also been evaluated using the same basis sets and molecules. Our computed results are compared with the available experimental values and are found to be in good agreement. The (RO)B3LYP and (RO)O3LYP procedures are found to produce reliable BDEs for the C-H bonds in toluene and the C-X (X = F, Cl) bond in alpha-substituted toluene (C6H5-CH2X) and their substituted derivatives. The substituent effect on the BDE values has been analyzed in terms of the ground-state effect and the radical effect. The effect of polarization of the C-H bond on the substituent effect is also analyzed. The BDE(C-H) and BDE(C-X) values for alpha-substituted (X = F and Cl) toluenes with a set of para substituents are presented for the first time.     

CH3SS与XO (X＝F,Cl, Br)反应过程的电子密度拓扑研究

利用MP2/6-311＋＋G(d,p)//B3LYP/6-311＋＋G(d,p)对CH3SS与XO (X＝F, Cl, Br)的反应机理进行了研究. 着重从电子密度拓扑分析角度讨论了化学键的生成和断裂. 计算结果表明单线态反应为主要反应通道, 且由于该通道的反应能垒低、放热明显, 说明CH3SS与XO在大气中比较容易进行. 电子密度拓扑分析表明, 在单线态抽氢反应通道中存在着四元环状过渡结构, 随着反应进行, 此四元环状过渡结构通过一个T-型结构变为三元环状过渡结构, 最后环状结构消失得到产物.     

Theoretical study on the hydrogen abstraction reaction of CH3CH2F (HFC-161) with Cl atom

A direct dynamics method is employed to study the hydrogen abstraction reaction of CH₃CH₂F+Cl. Three distinct transition states are located, one for -H abstraction and two for β-H abstraction. The potential energy surface (PES) information is obtained at the QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p), CCSD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) and G2//MP2/6-311G(d,p) level. Based on the QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) results, the rate constants of the three reaction channels are evaluated by using the canonical variational transition state theory (CVT) with small-curvature tunneling (SCT) contributions over the temperature range of 220–2800 K. The calculated results indicate that -H abstraction dominates the total reaction almost over the whole temperature range.     

CH_2X(X=H,FCI)与臭氧反应机理的最子化学研究

用量化学UMP2方法，在6-311++G**基组水平上研究了CH_2X(X=H,FCI)与臭氧反 应机理，全参数优化了反应过程中反应物、中间体、过渡态和产物的内何构型，在 UQCISD（T）/6-311++G**水平上计算了它们的能量，并对它们进行了振动分析，以 确定中间体和过渡态的直实性。从CH_2X(X=H,FCI)与O_3的反应机理的研究结果看 ，它们与O_3反应的活性都比较强，相对而言，活性大小顺序为CH_2F＞CH_3＞ CH_2CI，也就是说，CH_2F自由基与臭氧间的反应活性最强，对大气臭氧的损耗将 是最大的。同时研究还发现CH_2X(X=H,FCI)系列自由基与O_3的反应都是强放热反 应。     

DFT and ab initio direct dynamics studies on the hydrogen abstraction reactions of chlorine atoms with CH(4-n)F(n) (n = 1-3)

A direct dynamics study is carried out for the hydrogen abstraction reactions Cl + CH(4-n)F(n) (n = 1-3) in the temperature range of 200-1,000 K. The minimum energy paths (MEPs) of these reactions are calculated at the BH&H-LYP/6-311G(d,p) level, and the energies along the MEPs are further refined at the QCISD(T)/6-311+G(2df,2p) and QCISD(T)/6-311+G(d,p) (single-point) level. The rate constants obtained by using the improved canonical variational transition state theory incorporating small-curvature tunneling correction (ICVT/SCT) are in good agreement with the available experimental results. It is shown that the vibrational adiabatic potential energy curves for these reactions have two barriers, a situation similar to the analogous reactions CH(3)X+Cl (X=Cl, Br). The theoretical results show that for the title reactions the variational effect should not be neglected over the whole considered temperature range, while the small-curvature tunneling effect is only important in the lower temperature range. The effects of fluorine substitution on the rate of this kind of reactions are also examined.     

Quantum theory of atoms in molecules charge-charge flux-dipole flux models for the infrared intensities of X(2)CY (X = H, F, Cl; Y = O, S) molecules

The molecular dipole moments, their derivatives, and the fundamental IR intensities of the X2CY (X = H, F, Cl; Y = O, S) molecules are determined from QTAIM atomic charges and dipoles and their fluxes at the MP2/6-311++G(3d,3p) level. Root-mean-square errors of +/-0.03 D and +/-1.4 km mol(-1) are found for the molecular dipole moments and fundamental IR intensities calculated using quantum theory of atoms in molecules (QTAIM) parameters when compared with those obtained directly from the MP2/6-311++G(3d,3p) calculations and +/-0.05 D and 51.2 km mol(-1) when compared with the experimental values. Charge (C), charge flux (CF), and dipole flux (DF) contributions are reported for all the normal vibrations of these molecules. A large negative correlation coefficient of -0.83 is calculated between the charge flux and dipole flux contributions and indicates that electronic charge transfer from one side of the molecule to the other during vibrations is accompanied by a relaxation effect with electron density polarization in the opposite direction. The characteristic substituent effect that has been observed for experimental infrared intensity parameters and core electron ionization energies has been applied to the CCFDF/QTAIM parameters of F2CO, Cl2CO, F2CS, and Cl2CS. The individual atomic charge, atomic charge flux, and atomic dipole flux contributions are seen to obey the characteristic substituent effect equation just as accurately as the total dipole moment derivative. The CH, CF, and CCl stretching normal modes of these molecules are shown to have characteristic sets of charge, charge flux, and dipole flux contributions.     

N(4S)+CH3X(X=Cl、Br)反应机理的理论研究

The reaction of N(4S)+CH3X(X=Cl、Br) was studied by the ab initio method. The geometries of the reactants, transition states and products were optimized at the MP2/6-311+G(d,p) level. The corresponding vibration frequencies were calculated at the same level. The single-point calculations for all the stationary points were carried out at the MP2/6-311++G(3df,2p) and the QCISD(T)/6-311+G(d,p) levels using the MP2/6-311+G(d,p) optimized geometries. The energies of all the stationary points were calculated by the G2MP2 method. The results of this theoretical study indicate that the reaction has three reaction channels: H abstraction reaction channel a, Cl or Br abstraction reaction channel b and substitution reaction channel c. For the N(4S)+CH3Cl reaction, reaction channel a is the main reaction channel. Reaction channels b and c may have a slight contribution in the reaction. For the N(4S)+CH3Br reaction, reaction channel a is the main reaction channel. Reaction channels b and c may have some contribution in the reaction.