COS与O2反应的密度泛函理论研究

用量子化学密度泛函理论（DFT）方法，对COS与O2的反应进行了理论研究．在UB3LYP／6—31G^＊，UB3LYP／6—311＋＋G^＊＊水平上，优化了反应势能面上各驻点（反应物、产物、中间体和过渡态）的几何构型，在UB3LYP／6—31G^＊水平上通过内禀反应坐标（IRC）计算和振动分析，对过渡态进行了确认．在CCSD（T）／6—311＋＋G（2d，2p）水平上进行了单点能量计算，并确定了反应机理．研究结果表明，反应主要产物为CO2和SO.     

C2F51分子C—I键解离的自旋-轨道从头算研究

采用考虑相对论效应的6—311G^＊＊全电子基组与多参考微扰理论,计算了该分子的包含自旋-轨道耦合效应的垂直激发能和基态、激发态C—I键解离势能曲线．理论计算发现,势能曲线3^3A＂与1^1A＂,2^1A＇出现交叉,交叉区域在C—I键长为0．241nm附近;基态1^A＇到激发态3^3A＂（^3Q0）的垂直激发能为4．658eV,与实验值4．662eV非常吻合．讨论了C2F5I分子作为碘激光介质的可行性．     

平面型四核铜簇合物Cu4(CH2XMen)4桥键性质的密度泛函研究

采用密度泛函B3LYP方法，用LanL2DZ和6—31G^＊基组分别优化了平面型四核铜簇合物Cu4（CH2SiMe3）4和Cu4（CH2XMe2）4（X=P，As）的几何构型，并对B3LYP／LanL2DZ方法优化的结构进行了红外振动频率计算和自然键轨道分析．结果表明，簇合物均呈笼状结构，Cu—C—Cu三中心桥键之间电子的离域增强了Cu簇合物的稳定性，配位C原子的C—H平伏键与C—Cu配位键之间存在σ-超共轭效应．     

X-H (X = C, N, O, Si, P, S) 键离解能计算的密度泛函方法研究

使用了不同密度泛函方法计算X-H (X = C, N, O, Si, P, S) 键离解能，并分析不同密度泛函方法的计算精度。研究发现大多数密度泛函方法包括B3LYP, B3P86, B3PW91, G96LYP, PBE1PBE，和BH&HLYP都明显低估键离解能13-25 kJ/mol。该现象与是否使用无限基组无关，因为即使使用无限基组键离解能仍然被低估。因此密度泛函方法不适合用于键离解能的估算。其中B3P86方法的偏差最小。进一步分析表明，使用限制性开壳层计算并无任何优势，在大多数情况下非限制性开壳层计算实际上比限制性开壳层计算要好。最后，我们发现了密度泛函方法对键离解能的低估是系统的，因此建议利用校准后的UDFT/6-311++G(d, p)方法计算化学键离解能。     

（CH3）2S与HOCl分子间的卤键和氢键相互作用

在DFT-B3LYP/6-311＋＋G＊＊水平上分别求得（CH3）2S…ClOH卤键复合物和（CH3）2S…HOCl氢键复合物势能面上的稳定构型.频率分析表明,与单体HOCl相比,在两种复合物中,10Cl—11O和12H—11O键伸缩振动频率发生显著的红移.经MP2/6-311＋＋G＊＊水平计算的含基组重叠误差（BSSE）校正的气相中相互作用能分别为-11.69和-24.16kJ·mol^-1.自然键轨道理论（NBO）分析表明,在（CH3）2S…ClOH卤键复合物中,引起10Cl—11O键变长的因素包括两种电荷转移：（i）孤对电子LP（1S）1→σ＊（10Cl—11O）;（ii）孤对电子LP（1S）2→σ＊（10Cl—11O）,其中孤对电子LP（1S）2→σ＊（10Cl—11O）转移占主要作用,总的结果是使σ＊（10Cl—11O）的自然布居数增加0.14035e,同时11O原子的再杂化使其与10Cl成键时s成分增加,即具有与电荷转移作用同样的“拉长效应”;在（CH3）2S…HOCl氢键复合物中也存在类似的电荷转移,但是11O原子的再杂化不同于前者.自然键共振理论（NRT）进行键序分析表明,在卤键复合物和氢键复合物中,10Cl—11O和12H—11O键的键序都减小.通过分子中原子理论（AIM）分析了复合物中卤键和氢键的电子密度拓扑性质.     

硝酸甲与硝酸乙酯的结构、振动频率和热力学性质的密度函数 理论(DFT)研究

用密度函数理论(DFT)的BLYP和B3LYP方法,取6-31G,6-31G^*,6-31G^*^*,6-311G,6-311G^*和6-311G^*^*六种基组,对硝酸甲酯和硝酸乙酯的几何构型和红外振动频率进行了计算研究.结果表明,B3LYP方法在采用极化基组(6-31G^*,6-31G^*^*,6-311G^*和6-311G^*^*)时计算得到的结果均较好,适用于硝酸酯类化合物的研究.而BLYP方法无论采用何种基组均不适用;运用校正后的B3LYP/6-31G^*频率(校正因子0.975)计算得到的热力学性质(C^o~p,H^o和S^o)与实验结果较吻合。     

硝酸甲与硝酸乙酯的结构、振动频率和热力学性质的密度函数 理论(DFT)研究

用密度函数理论(DFT)的BLYP和B3LYP方法,取6-31G,6-31G^*,6-31G^*^*,6-311G,6-311G^*和6-311G^*^*六种基组,对硝酸甲酯和硝酸乙酯的几何构型和红外振动频率进行了计算研究.结果表明,B3LYP方法在采用极化基组(6-31G^*,6-31G^*^*,6-311G^*和6-311G^*^*)时计算得到的结果均较好,适用于硝酸酯类化合物的研究.而BLYP方法无论采用何种基组均不适用;运用校正后的B3LYP/6-31G^*频率(校正因子0.975)计算得到的热力学性质(C^o~p,H^o和S^o)与实验结果较吻合。     

惰性气体体系HeAuF和NeAuF的理论预测

近年来,含惰性气体原子的体系受到越来越多的关注,特别是NgMX类体系,其中Ng表示惰性气体原子,M是金属原子Cu、Ag、Au,X为卤素原子F、Cl、Br．已有的实验和理论研究结果表明,这类体系为线性构型,并且表现出较强的Ng-M共价键的特征,含Au的体系共价键最强,其次是含Cu的,含Ag的作用最弱．对惰性气体原子而言,共价键强弱的顺序依次为：Ar〈Kr〈Xe．因此Xe—Au共价键是目前发现最强的,而Ar-Ag是最弱的．本文旨在探讨是否存在He—Au和Ne-Au共价键．我们将采用高精度量子化学方法对HeAuF和NeAuF体系的结构和稳定性进行理论研究．在MP2和CCSD（X）水平上优化了HeAuF和NeAuF的平衡构型．对Au原子,采用Stuttgart和Koeln的有效势（ECP60MWB）和价电子基组（9s9p6d4f）,对He、Ne和F采用aug-cc-pVTZ基组（定为基组A）．在CCSD（T）优化结构的基础上采用基组A和B（其中对He、Ne和F采用aug-cc-pVQZ基组）计算了NgAuF→Ng＋AuF的离解能．CCSD（T）计算结果表明He—Au的键长为0．184nm,和vdW半径极限（0．217nm）相比,更接近共价半径极限（0．167-0．187nm）;而Ne—Au的键长为0．248nin,明显长于共价半径极限（0．197mm）,更接近vdW半径极限（O．231nm）．另外从同类体系键长规律来看,He—Au的键长依次小于Ar-Au、Kr-Au和Xe—Au,而Ne-Au却长于Ar-Au和Kr-Au的键长．对比Ng—Au谐性振动频率发现,频率大小依次为Ne—Au〈He—Au〈Ar-Au〈Kr—Au〈Xe-Au．在CCSD（T）／基组A水平下得到的HeAuF体系He—Au键的离解能为24kJ·mol^-1,明显高于典型vdW体系ArNaCl的Ar-Na离解能（8kJ·mol^-1）;而NeAuF体系离解为Ne和AuF的离解能只有7kJ·mol^-1．结果表明HeAuF体系可能存在较弱的He—Au共价键,而NeAuF不是一个稳定体系．为了更精确确定HeAuF的稳定性,我们采用了更大的基组B’（在基组B基础上对Au采用cc—pVTZ-PP基组）对HeAuF体系进行进一步计算．在CCSD（T）／基组B’水平下得到的He-Au键的离解能为26kJ·mol^-1．自然键轨道（NBO）分析表明,HeAuF体系中有明显的从He到AuF电荷转移0．06e,而NeAuF体系的电荷转移仅为0．02e,ArAuF,KrAuF和XeAuF相应的值分别为0．11、0．13和0．26e．该趋势与另外一类惰性气体体系HNgF非常类似．另外从体系价层轨道的电子密度图发现,HeAuF的价层轨道4σ有明显的He的1s和Au的5Pz轨道的重叠,而NeAuF未出现Ne和Au轨道重叠的情况．     

木质素模化物键离解能的理论研究

采用密度泛函理论B3P86方法,在6-31G(d,p)基组水平上,对木质素结构中的6种连接方式(β-O-4、α-O-4、4-O-5、β-1、α-1、5-5)的63个木质素模化物的醚键(C-O)和C-C键的键离解能E_B进行了理论计算研究。分析了不同取代基对键离解能的影响以及键长与键离解能的相关性。计算结果表明,C-O键的键离解能通常比C-C键的小,在各种醚键中C_α-O键的平均键离解能最小,为182.7 kJ/mol;其次是β-O-4连接中的C_β-O键,苯环和烷烃基上的取代基对醚键的键离解能有较强的弱化作用,C-O键的键长和键离解能的相关性较差。与C-O键相比,C-C键的键离解能受苯环上取代基的影响很小,而烷烃基上的取代基对C-C键的键离解能有较大的影响,C-C键的键离解能和键长之间存在较强的线性关系,C-C键的键长越长,其键离解能越小。     

1, 2-二硒方酸气相酸性和芳香性的理论研究

在RHF/6-311G^*^*水平优化得到1,2-二硒方酸(3,4-二羟基-3-环丁烯-1,2-二硒酮)三种平面构象异构体的平衡几何构型。进一步用MP2(full)/6-311G^*//RHF/6-311G^*^*方法计算三种异构体的单点能量,发现ZZ型异构体是能量最低构象,且ZZ和ZE型能量非常接近。用优化的最稳定构象ZZ型异构体在RHF/6-311G^*^*//RHF/6-311G^*^*,RHF/6-311+G^*^*//RHF/6-311+G^*^*,MP2(full)/6-311+G^*^*//RHF/6-311+G^*^*和B3LYP/6-311+G^*^*//B3LYP/6-311+G^*^*水平计算其气相酸性[ΔGⅲ~(~2~9~8~K~)]和同键反应芳香性稳定化能(HASE)。用基团加和法(groupincrementapproach)在RHF/6-311+G^*^*//RHF/6-311+G^*^*和B3LYP/6-311+G^*^*//B3LYP/6-311+G^*^*水平计算其磁化率增量(Λ)。计算结果指出标题化合物的键长发生了平均化,同键反应芳香性稳定化能和磁化率增量均为负值,表明它具有芳香性,实现了标题化合物芳香性的几何、能量和磁性的判定。     

DFT study of the C-Cl bond dissociation enthalpies and electronic structure of substituted chlorobenzene compounds

Quantum chemical calculations were used to estimate the bond dissociation energies (BDEs) for 13 substituted chlorobenzene compounds. These compounds were studied by employing the hybrid density functional theory methods (B3LYP, B3PW1, B3P86) with 6-31G** and 6-311G** basis sets. It was demonstrated that B3P86/6-311G** method is the best method for computing the reliable BDEs for substituted chlorobenzene compounds which contain the C-Cl bond. It was found that the C-Cl BDE depends strongly on a computational method and basis set used. Substitution effect on the C-Cl BDE of substituted chlorobenzene compounds is further discussed. It is shown that the effects of substitution on the C-Cl BDE of substituted chlorobenzene compounds are very insignificant. Frontier orbital energy gap of studied compounds was also investigated. From the data on frontier orbital energies gap, we estimated the relative thermal stability of substituted chlorobenzene compounds.     

Density Functional Calculations of C-NO2 Bond Dissociation Energies for Nitroalkanes Molecules

Bond dissociation energies for the removal of nitrogen dioxide group in some nitroalkane energetic materials have been calculated by using the three hybrid density functional theory （DFT） methods B3LYP, B3PW91 and B3P86 with 6-31g^＊＊ and 6-311g^＊＊ basis sets. The computed BDEs have been compared with the available experimental results. It is found that the B3P86 method with 6-31g^＊＊ and 6-311g^＊＊ basis sets can obtain satisfactory bond dissociation energies （BDEs）, which are in extraordinary agreement with the experimental data. Considering the smaller mean absolute deviation and maximum difference, the reliable B3P86/6-311g^＊＊ method was recommended to compute the BDEs for the removal of nitrogen dioxide group in the nitroalkane energetic materials. Using the method, the BDEs of 8 other nitroalkane energetic materials have been calculated and the maximum difference from experimental value is 1.76 kcal.mo1^-1 （for the BDE of tC4Hg-NOz）, which further proves the reliability of B3P86/6-311g^＊＊ method. In addition, it is noted that the BDEs of C-NO2 bond change slightly for main chain nitroalkane compounds with the maximum difference of only 3.43 kcal mo1^-1.     

Density Functional Calculations of C–NO_2 Bond Dissociation Energies for Nitroalkanes Molecules

Bond dissociation energies for the removal of nitrogen dioxide group in some nit- roalkane energetic materials have been calculated by using the three hybrid density functional theory (DFT) methods B3LYP, B3PW91 and B3P86 with 6-31g** and 6-311g** basis sets. The computed BDEs have been compared with the available experimental results. It is found that the B3P86 method with 6-31g** and 6-311g** basis sets can obtain satisfactory bond dissociation energies (BDEs), which are in extraordinary agreement with the experimental data. Considering the smaller mean absolute deviation and maximum difference, the reliable B3P86/6-311g** method was recommended to compute the BDEs for the removal of nitrogen dioxide group in the nitroalkane energetic materials. Using the method, the BDEs of 8 other nitroalkane energetic materials have been calculated and the maximum difference from experimental value is 1.76 kcal·mol-1 (for the BDE of tC4H9–NO2), which further proves the reliability of B3P86/6-311g** method. In addition, it is noted that the BDEs of C–NO2 bond change slightly for main chain nitroalkane compounds with the maximum difference of only 3.43 kcal mol-1.     

DFT Calculations and NBO Analysis of 2-Chloroethylethyldichlorosilane Unimolecular Elimination Kinetics in the Gas Phase

Ab initio molecular orbital calculations have been used to investigate the thermal decomposition kinetics of 2-chloroethylethyldichlorosilane at the B3LYP/6-311+G**,B3PW91/6-311+G**,and MPW1PW91/6-311+G** levels of theory.Among these methods,the results(activation parameters) obtained using the B3LYP/6-311+G** level are in good agreement with the available experimental data.The calculated data imply that in the unimolecular β-elimination reactions of the studied compound in the gas phase,the polarization of C(1)-Cl(3) and C(1)-H(4) bonds in the sense of C(1)δ+-Cl(3)δ-and C(1)δ+-H(4)δ-,respectively,is a determining factor in the gas phase elimination reactions 1,2 and 3.Analysis of bond order,natural bond orbital charges,bond indexes,synchro-nicity parameters,and IRC calculations suggest the elimination of 2-chloroethylethyldichlorosilane via reactions 1～3 can be described as concerted and slightly asynchronous.The transition state structures of these reactions are a four-membered cyclic structure.     

Theoretical Study of the N-NO2 Bond Dissociation Energies for Energetic Materials with Density Functional Theory

The N-NO2 bond dissociation energies （BDEs） for 7 energetic materials were computed by means of accurate density functional theory （B3LYP, B3PW91 and B3P86） with 6-31G＊＊ and 6-311G＊＊ basis sets. By comparing the computed energies and experimental results, we find that the B3P86/6-311G＊＊ method can give good results of BDE, which has the mean absolute deviation of 1.30kcal/mol. In addition, substituent effects were also taken into account. It is noted that the Hammett constants of substituent groups are related to the BDEs of the N-NO2 bond and the bond dissociation energies of the energetic materials studied decrease when increasing the number of NO2 group.     

Bond Dissociation Energies and Electronic Structures in a Series of Peroxy Radicals: A Theoretical Study

Quantum chemical calculations are performed to estimate the bond dissociation energies (BDEs) for 18 peroxy radicals. Since DFT methods are researched to have low basis sets sensitivity, these radicals are studied by utilizing the hybrid density functional theory (DFT) (B3LYP, B3P86, B3PW91 and PBE1PBE) in conjunction with the 6‐311G** basis set and the complete basis set (CBS‐Q) method. On the basis of comparisons of the computational results and the experimental values, we evaluate the effectiveness of above methods. It is demonstrated that CBS‐Q method is the best method for computing the reliable BDEs of C—OO bond, with the average absolute errors of 2.1 kcal/mol. So CBS‐Q method is suitable to predict accurate BDEs of C‐OO bond for peroxy compounds. The computational energy gaps between the HOMO and LUMO of studied compounds are almost identical from the point of view of stability and substantial HOMO‐LUMO gaps for all molecules suggest their electronic stability. In addition, substituent effect on the C—OO BDE of peroxy radicals is analyzed. It is noted that the effects of substitution on the C—OO BDE of peroxy radicals are significant. Our results will shed lights on future theoretical and experimental work.     

PCM study of the solvent and substituent effects on bond dissociation energies of the C?NO bond

Quantum chemical calculations are used to estimate the equilibrium C? NO bond dissociation energies (BDEs) for eight X? NO molecule (X = CCl₃, C₆F₅, CH₃, CH₃CH₂, iC₃H₇, tC₄H₉, CH₂CHCH₂, and C₆H₅CH₂). These compounds are studied by employing the hybrid density functional theory (B3LYP, B3PW91, B3P86) methods together with 6‐31G** and 6‐311G** basis sets and the complete basis set (CBS‐QB3) method. The obtained results are compared with the available experimental results. It is demonstrated that B3P86/6‐31G** and CBS‐QB3 methods are accurate for computing the reliable BDEs for the X? NO molecule. Considering the inevitably computational cost of CBS‐QB3 method and the reliability of the B3P86 calculations, B3P86 method with 6‐31G** basis set may be more suitable to calculate the BDEs of the C? NO bond. The solvent effects on the BDEs of the C? NO bond are analyzed and it is shown that the C? NO BDEs in a vacuum computed by using B3PW91/6‐311G** method are the closest to the computed values in acetontrile and the average solvent effect is 1.48 kcal/mol. Subsequently, the substituent effects of the BDEs of the C? NO bond are further analyzed and it is found that electron denoting group stabilizes the radical and as a result BDE decreases; whereas electron withdrawing group stabilizes the group state of the molecule and thus increases the BDE from the parent molecule. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

密度泛函法计算C1-C14正构烃生成焓及C-C键裂解能

采用密度泛函理论(DFT)方法研究了费托石脑油裂解反应中涉及到C₁-C₁₄正构烃和自由基中间体的生成焓及其C-C键解离能(BDE)。 结果表明,在所有评价的密度泛函理论方法（B97-1、BB1K、B1B95、MPWB1K和MPW1B95）中,MPW1B95/6-311G(d,p)方法计算最精确。 以此方法为基准,进一步对高碳烃及其裂解产物的标准生成焓和C-C键解离能进行了预测。 与可得到的实验数据相比,MPW1B95/6-311G(d,p)方法预测的烃和自由基的平均生成焓分别为0.8和2.7 kJ/mol,C-C键解离能的平均绝对误差只有3.1 kJ/mol,表明此方法不仅可准确计算正构烃标准生成焓和C-C键解离能,而且还能正确预测C-C键解离能变化趋势。     

PCM Study of Bond Dissociation Energies of the S–NO Bond: A DFT Study

Quantum chemical calculations are used to estimate the equilibrium S–NO bond dissociation energies (BDEs) for seven S-nitro-N-acetyl-d,l-penicillamine dipeptides (SNAP-two peptides) in acetonitrile solution. These compounds were studied by employing the hybrid density functional theory (B3LYP, B3P86 and B3PW91) methods together with the 6-31G** basis set. The obtained results are compared with available experimental results. It is demonstrated that B3PW91 method is the best method to compute the bond dissociation energies of SNAP-two peptides. The substituent and solvent effects of the S–NO BDEs are further analyzed. The results show that S–NO BDE increases with the increment of isoelectric points of substituted groups. In addition, the S–NO BDE decreases due to the inclusion of solvent effects. Furthermore, SNAP-two peptides and the other NO-donors are compared.     

Investigation of correlation between impact sensitivities and bond dissociation energies in some triazole energetic compounds

In this article, density functional theory has been utilized to study on the correlation between impact sensitivities h _50% and the bond dissociation energies (BDEs) of nine triazole energetic explosives. By employing B3LYP and B3P86 method with the 6-311G** basis set, all the molecules have been fully optimized. The BDEs for removal of the NO₂ group in these compounds have also been calculated at the same level. Computed results show that BDEs calculated by B3LYP method are all less than those by B3P86 method. The relationship between the impact sensitivities and the weakest C–NO₂ bond dissociation energy (BDE) values have been investigated. The results indicate a good linear correlation between the impact sensitivity h _50% and the ratio (BDE/E) of the weakest BDE to the total energy E.