偶氮四唑非金属盐的分子结构和热分解参数与其撞击感度间关联关系的研究

采用DFT 方法对偶氮四唑胍盐(GZT)等6 种偶氮四唑非金属盐进行了B3LYP/6-31G 水平的几何结构全优化, 计算了分子的最低空轨道能量(E_LUMO)、最高占有轨道能量(E_HOMO)及原子上的电荷分布等分子结构参数. 研究了偶氮四唑非金属盐的分子结构参数(包括氧平衡)和热分解参数(热分解温度和热分解活化能等)与其撞击感度之间的相关性. 结果表明, 偶氮四唑非金属盐的氧平衡越低, 撞击感度越低; 热分解温度和热分解活化能越低, 撞击感度越高; 取代基团上所带正电荷越大, 撞击感度越低.     

椅式(8,8)单壁碳纳米管内偶氮苯的顺反异构化

采用ONIOM(B3LYP/6-31+G*:UFF)方法对偶氮苯受限于椅式单壁碳纳米管CNT(8,8)内的结构、电子光谱和热致顺反异构化势能面进行了计算.结果表明,反式偶氮苯进入CNT(8,8)碳纳米管内是一个放热过程.由于碳纳米管的限制作用,偶氮苯分子的苯环绕CN键有一定的旋转,反式偶氮苯的平面结构发生扭曲,但其它结构参数变化不明显.受限于碳纳米管内的偶氮苯的顺反异构体能量差比非受限状态下增加了8.1 kJ.mol-1,表明碳纳米管的限制作用对偶氮苯顺反异构体的相对热稳定性有一定影响.光谱计算表明,受限于碳纳米管内的偶氮苯最低的三个单线态吸收波长仅蓝移1-5 nm.异构化势能面计算发现,偶氮苯在CNT(8,8)碳纳米管内发生反式到顺式异构化的能垒增加,由顺式回复到反式异构体的能垒无明显变化,表明CNT(8,8)碳纳米管的限制作用可以抑制偶氮苯从反式到顺式异构体的热致异构化过程.此外,受限状态下的偶氮苯主要通过CNN键角反转发生热致异构化.     

单重态CCl2与O3反应机理的理论研究

用量子化学从头算方法,研究了单重态CCl_2与O_3反应的机理.在HF/6-31G(d)水平上用梯度解析技术全参数优化上述反应的反应物、中间体、过渡态和产物构型,MP2/6-31G(d)//HF/6-31G(d)方法计算能量.给出了有关化合物的结构数据.结果表明:CCl_2与O_3首先生成富能中间体CCl_2O_3,然后中间体裂解生成CCl_2O和O_2.该反应为强放热反应,放出的热量为516.88kJ·mol~(-1)[MP2/6-31G(d)//HF/6-31G(d)].通过内禀反应坐标(IRC)计算,获得了沿反应途径的势能剖面.     

苯和苯胺类硝基衍生物热解机理和撞击感度的理论研究

在DFT-B3LYP/6-31G*水平获得苯和苯胺类硝基衍生物的全优化分子几何和电子结构.通过非限制性(U)B3LYP/6-31G*计算求得标题物各化学键离解能(BDE).用UHF-PM3 MO方法求得引发键C-NO2键均裂反应的活化能(Ea).以静态指标(键集居数、前线轨道能级差和硝基上净电荷)和动态理论指标(BDE和Ea)阐明了热解引发机理,关联了实验撞击感度.运用SPSS程序关联静态和动态理论指标,表明它们均可以用作预示标题物的热解引发机理和撞击感度.     

高能化合物热解机理和撞击感度的理论研究-苯和苯胺类硝基衍生物

在DFT-B3LYP/6-31G*水平求得苯和苯胺类硝基衍生物的全优化分子几何和电子结构。通过非限制性 (U) B3LYP/6-31G*计算求得标题物各化学键离解能(BDE)。用UHF-PM3 MO方法求得引发键C-NO2键均裂反应的活化能(Ea)。以静态指标(键集居数、前线轨道能级差和硝基上净电荷)和动态理论指标(BDE和Ea) 阐明了热解引发机理,关联了实验撞击感度。运用SPSS程序关联静态和动态理论指标,表明它们可平行或等价地用作预示标题物的热解引发机理和撞击感度。     

不同取代基巯基偶氮苯甲酸TiO2光催化降解的密度泛函理论研究

以TiO2为光催化剂,对含六种不同取代基的巯基偶氮苯甲酸进行了光催化降解实验,并采用量子化学密度泛函理论B3LYP/6-31G**计算了巯基偶氮苯甲酸的电子结构,研究了巯基偶氮苯甲酸中推、拉电子取代基对其光催化降解活性的影响.结果表明,拉电子基(-COOH或-SO3H)的引入使巯基偶氮苯甲酸的偶极矩增大,绝对电负性减小,最高占有轨道能量升高,HOMO-LUMO能隙和C-N键的键级降低,从而提高了巯基偶氮苯甲酸的光催化降解活性,而分子中推电子基团的影响则相反.     

液晶作为气相色谱固定液的研究

1888年F.Reinitzer发现了液晶,直到二十世纪五十年代人们才对液晶的应用作了尝试。1963年H.Kelker首先提出用对.对’-二甲氧基氧化偶氮苯液晶作为气相色谱固定液分离二甲苯异构体.近年来液晶固定相在苯二取代化合物、多环芳烃、甾族差向异构体、取代萘类、氮杂环化合物、顺反异构体、对映异构体等的分离研究取得了较大进展.     

1,2,3-三唑类含能化合物的合成研究进展

1,2,3-三唑类含能化合物具有密度高、生成焓高及热稳定性好等特点,近年来受到科研人员的广泛研究和报道.较全面地综述了单环1,2,3-三唑、多环1,2,3-三唑和稠环1,2,3-三唑类含能化合物的研究成果,对其合成方法、感度和爆轰性能进行了总结.结果表明,1,2,3-三唑含能化合物具有良好的爆轰性能和热稳定性,在含能材料领域具有重要的研究价值和应用潜力.最后,对1,2,3-三唑含能化合物的发展现状进行总结与展望,提出了各类1,2,3-三唑含能化合物未来可能的应用方向,以期为从事含能化合物研究的科研人员提供一定的参考.     

特殊氢方法预测丙氨酸-α-四肽构象稳定性

使用B3LYP/6-31G*方法优化得到了丙氨酸-α-四肽分子的48个稳定构象.定义了丙氨酸-α-多肽中的特殊氢原子,对构象中与特殊氢原子有关的主要非键相互作用进行分析,提出了使用与特殊氢原子有关的非键相互作用来预测多肽构象的相对稳定性,并将之称为特殊氢方法.基于丙氨酸-α-二肽分子的5个构象和三肽分子的7个构象,确定了特殊氢方法中的与特殊氢原子有关的非键相互作用参数.利用特殊氢方法预测丙氨酸-α-四肽分子的48个构象的相对稳定性,与B3LYP/6-31G*方法的相对能量比较,得到了满意的结果.特殊氢方法得到的相对能量(Y)和B3LYP/6-31G*方法得到的相对能量(X)的线性相关方程为Y=0.9296X 2.2041,相关系数R=0.9532,标准偏差为3.0kJ/mol,偏差绝对值≤4.18kJ/mol的构象占87.5%.     

含能材料的密度、爆速、爆压和静电感度的理论研究

用密度泛函理论(DFT) B3LYP方法, 在6-31G*基组水平下, 全优化计算了系列硝胺类和硝基芳烃类爆炸物的几何构型, 用Monte-Carlo方法和自编程序, 基于0.001 e•bohr^－3等电子密度面所包围的体积空间求得分子平均摩尔体积(V)和理论密度(ρ). 用Kamlet-Jacobs方程基于理论密度(ρ)和PM3计算生成焓(ΔH_f)估算标题物的爆速(D)和爆压(p), 发现多环硝胺类化合物的爆轰性能优于芳烃硝基类化合物, 故此, 在寻求高能量密度材料(HEDM)时, 我们应特别关注多环硝胺化合物. 与ρ和D文献值比较, 表明本理论计算方法和结果是适用可靠的. 将爆速(D)和爆压(p)计算值与静电感度实验值(E_ES)进行比较和关联, 发现: 若对化合物进行细致分类讨论, 则它们之间存在较好的线性关系. 据此建议, 在含能材料分子设计中, 可通过理论计算爆轰性质(D或p)去预估难以定量求得或尚未合成的含能材料的静电火花感度值(E_ES). 此外, 我们还讨论了取代基对ρ, D和p的影响, 也有助于分子设计.     

Conformational study of the structure of 12-crown-4-alkali metal cation complexes

A conformational search was performed for the 12-crown-4 (12c4)-alkali metal cation complexes using two different methods, one of them is the CONFLEX method, whereby eight conformations were predicted. Computations were performed for the eight predicted conformations at the HF/6-31+G*, MP2/6-31+G*//HF/6-31+G*, B3LYP/6-31+G*, MP2/6-31+G*//B3LYP/6-31+G*, and MP2/6-31+G* levels. The calculated energies predict a C4 conformation for the 12c4-Na+, -K+, -Rb+, and -Cs+ complexes and a C(s) conformation for the 12c4-Li+ complex to be the lowest energy conformations. For most of the conformations considered, the relative energies, with respect to the C4 conformation, at the MP2/6-31+G*//B3LYP/6-31+G* are overestimated, compared to those at the MP2/6-31+G* level, the highest level of theory considerd in this report, by 0.2 kcal/mol. Larger relative energy differences are attributed to larger differences between the B3LYP and MP2 optimized geomtries. Binding enthalpies (BEs) were calculated at the above-mentioned levels for the eight conformations. The agreement between the calculated and experimental BEs is discussed.     

吲哚分子振动光谱的密度泛函理论研究

用多种密度泛函理论(DFT)方法(BLYP/6-31G^*^*,B3LYP/6-31G^*^*,B3PW91/6-31G^*^*和SVWN/6-31G^*^*)对吲哚分子的平衡几何构型进行了优化。在优化构型的基础上计算了吲哚分子的谐力场、振动基频和红外光谱强度。计算得到的振动频率与实验值比较平均偏差对四种计算方法(BLYP/6-31G^*^*,P3LYP/6-31G^*^*,B3PW91/6-31G^*^*和SVWN/6-31G^*^*)分别为16.3,40.5,45.1和26.4cm^-^1。BLYP/6-31G^*^*理论力场被用于吲哚分子的简正坐标分析计算中。根据振动率的势能分布(PEDs)对此分子的振动基频进行了理论归属。     

Quantum chemical studies on structure, conformation and isomerization of nitroso, nitro substituted benzene and 1,3-cyclopentadiene

The molecular structure, conformational stability and isomerization of nitroso, nitro substituted benzene and 1,3-cyclopentadiene in gas phase have been investigated using ab initio and density functional theory methods. The molecular geometries and energetics of possible conformers were obtained by employing MP2, B3LYP and B3PW91 levels of theory implementing 6-31G* basis set. The relative stabilities of the conformations were evaluated from the energy differences of the structure. Chemical hardness (η) and chemical potential (μ) were calculated at HF/6-31G* level of theory for all the positional and geometrical isomers to study the maximum hardness principle. Each optimized structure has been tested against the imaginary frequencies at MP2/6-31G* level of theory in order to be sure they are located at energy minimum.     

An Ab Initio Molecular Orbital Study of Valence Bond Isomers of Silabenzene and Phosphabenzene

Ab initio molecular orbital calculation at HF/6-31G*, HF/6-31G**, HF/6-311G**, HF/6-311++G**, RMP2-FC/6-31G*, and B3LYP/6-31G* levels of theory for geometry optimization and MP4(SDQ)/6-31G* for a single point total energy calculation are reported for silabenzene ( 7 ), phosphabenzene ( 8 ) and 16 valence bond isomers of silabenzene and phosphabenzene ( 9-24 ). The calculated energy difference (19.78 kcal mol m 1 ) between silabenzene and the most stable valence bond isomer of silabenzene (1-silabenzvalene, 9 ) is much smaller than the difference (73.60 kcal mol m 1 ) between benzene and benzvalene ( 2 ). The energy difference between phosphabenzene and the most stable valence bond isomer of phosphabenzene (1-phosphabenzvalene, 17 ) is calculated to be 43.29 kcal mol m 1 .     

Accuracy of the three-body fragment molecular orbital method applied to Møller-Plesset perturbation theory

The three‐body energy expansion in the fragment molecular orbital method (FMO) was applied to the 2nd order Møller–Plesset theory (MP2). The accuracy of both the two and three‐body expansions was determined for water clusters, alanine n‐mers (α‐helices and β‐strands) and one synthetic protein, using the 6‐31G* and 6‐311G* basis sets. At the best level of theory (three‐body, two molecules/residues per fragment), the absolute errors in energy relative to ab initio MP2 were at most 1.2 and 5.0 mhartree, for the 6‐31G* and 6‐311G* basis sets, respectively. The relative accuracy was at worst 99.996% and 99.96%, for 6‐31G* and 6‐311G*, respectively. A three‐body approximation was introduced and the optimum threshold value was determined. The protein calculation (6‐31G*) at the production level (FMO2/2) took 3 h on 36 3.2‐GHz Pentium 4 nodes and had the absolute error in the MP2 correlation energy of only 2 kcal/mol. © 2007 Wiley Periodicals, Inc. J Comput Chem, 2007     

Comparison of levels of electronic structure theory in direct dynamics simulations of C2H5F --> HF + C2H4 product energy partitioning

Direct dynamics simulations at the MP2/6-311++G** level of theory were performed to study C(2)H(5)F --> HF + C(2)H(4) product energy partitioning. The simulation results are compared with experiment and a previous MP2/6-31G* simulation. The current simulation with the larger basis set releases more energy to HF vibration and less to HF + C(2)H(4) relative translation as compared to the previous simulation with the 6-31G* basis set. The HF rotation and vibration energy distributions determined from the current simulation are in overall very good agreement with previous experimental studies of C(2)H(5)F dissociation by chemical activation and IRMPA. A comparison of the simulations with experiments suggests there may be important mass effects for energy partitioning in HX elimination from haloalkanes. The transition state (TS) structures and energies calculated with MP2 and the 6-31G* and 6-311++G** basis sets are compared with those calculated using CCD, CCSD, CCSD(T), and the 6-311++G** basis set.     

An Ab Initio Molecular Orbital Study of Isophosphinines

Ab initio molecular orbital calculation at HF/6-31G*, HF/6-31G**, HF/6-311G**, HF/6-311++G**, RMP2-FC/6-31G*, and B3LYP/6-31G* levels of theory for geometry optimization and MP4(SDQ)/6-31G* for a single-point total energy calculation are reported for phosphinine and 13 isophosphinines 7-19 . Isomers 7-11 with an allenic system are calculated to be 8-18 kcal mol m 1 more stable than structures 12-17 with an acetylenic moiety. The calculated energy difference (66.19 kcal mol m 1 ) between phosphinine and the most stable isophosphinine (1-phospha-1,2,4-cyclohexatriene, 10 ) is smaller than the difference (78.96 kcal mol m 1 ) between benzene and the most stable isobenzene (cyclohexa-1,2,4-triene, 2 ). The isophosphinines 18 and 19 , with a butatriene moiety, are calculated to be the least stable isomers.     

1. 5-环辛二烯-3, 7-二炔(C~8H~4)结构和光谱性质的理论研究

采用量子化学从头算方法研究了1,5-环辛二烯-3,7-二炔(C~8H~4)的结构和光谱性质,根据等键反应分析和自然键轨道方法研究了它的稳定性、成键情况和共轭性。结果表明1,5-环辛二烯-3,7-二炔(C~8H~4)分子为平面刚性结构,可能稳定存在。分子中C≡C键与C=C键存在一定程度的共轭,可能具有芳香性。     

Calculated structures and relative stabilities of furoxan,some 1,2-dinitrosoethylenes and other isomers

The structures and relative stabilities of furoxan and some of its isomers, e.g., the 1,2-dinitrosoethylenes, have been determined by means of ab initio Hartee–Fock and Møller–Plesset calculations. Geometries were optimized at the HF/3-21G, HF/6-31G* and MP2/6-31G* levels, and subsequently used for computing MP2/6-31G*, MP3/6-31G*, and MP4/6-31G* energies. The results are markedly affected by the inclusion of electronic correlation, which renders three of the isomers unstable. It also emphasizes the importance of a zwitterionic contribution to the structure of furoxan, which promotes ring-opening through a cis 1,2-dinitrosoethylene intermediate/transition state that has an MP4/6-31G*//MP2/6-31G* energy that is 31.6 kcal/mol above furoxan.     

Rearrangement of azirine intermediates to nitriles: Theoretical study of cleavage of 3,4-dihydro-1aH-azirine[2,3-c]pyrrol-2-one to cyanoketene–formaldimine complex

A detailed investigation of the reaction path for the thermal rearrangement of 3,4-dihydro-1aH-azirine[2,3-c]pyrrol-2-one to yield a cyanoketene–formaldimine complex is carried out at the MP2/6-31G* and B3LYP/6-31G* levels of theory. The ring opening of the five-membered pyrrolinone ring and the formation of the nitrile group takes place in a concerted manner, presenting a significant strain energy release and allowing for an electronic stabilization by coarctate conjugation of the transition structure (TS). These two factors make possible a moderate energy barrier. Although the structural features B3LYP/6-31G* theoretical levels, it is found that the MP2 energy barrier (28.8) CCSD(T)/6-31G*//MP2/6-31G* value (17.1 kcal/mol). The complex electronic rearrangement can be rationalized using the theory of coarctate transition structures developed by Herges as the evolution of an azirine structure without referring to a hypothetical vinyl nitrene intermediate. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 912–922, 1998