多异氰基立方烷生成热的理论研究——高能量密度材料(HEDM)的寻求

运用密度泛函理论(DFT)B3LYP方法和半经验MO(MINDO/3, MNDO, AM1和PM3)方法系统计算了全部21种多异氰基立方烷的生成热. 首先, 在DFT-B3LYP/6- 31G*水平下通过不破裂立方烷笼状骨架(亦即选择立方烷为参考物)的等键反应设计, 精确计算了8种多异氰基立方烷的生成热; 发现B3LYP/6-31G*结果分别地均与上述4种半经验MO方法求得的生成热之间存在良好的线性关系(相关系数均在0.9971以上), 且以PM3法得到的生成热结果最好. 其次, 其他13种多异氰基立方烷的精确生成热借助上述线性关系通过校正对应的PM3结果而获得. 多异氰基立方烷的生成热很高, 且随—CN基数目的增加而线性地增大, 表明它们属于极具潜力的“新一代高能量密度材料(HEDM)”而具开发价值.     

多氰基立方烷生成热的DFT-B3LYP和半经验MO研究

运用密度泛函理论(DFT)B3LYP方法和半经验MO(MINDO/3,MNDO,AM1和PM3)方法系统计算了全部21种多氰基立方烷的生成热,首先,在DFT-B3LYP/6-31G^*水平下通过不破裂立方烷笼状骨架(亦即选择立方烷为参考物)的等键反应设计,精确计算了9种多氰基立方烷的生成热;发现B3LYP/6-31G^*结果分别地均与上述四种半经验MO方法求得的生成热之间存在良好的线性关系(相关系数均在0.9994以上),且以AM1生成热与B3LYP/6-31G^*计算值最为接近,其次,其它12种多氰基立方烷的精确生成热借助上述线性关系通过校正对应的AM1结果而获得,多氰基立方烷的生成热很高,且随-CN基数目的增加而线性地增大,表明它们属于极具潜力的“新一低高能炸药”而具开发价值。     

多氰基立方烷生成热的DFT-B3LYP和半经验MO研究

运用密度泛函理论(DFT)B3LYP方法和半经验MO(MINDO/3,MNDO,AM1和PM3)方法系统计算了全部21种多氰基立方烷的生成热,首先,在DFT-B3LYP/6-31G^*水平下通过不破裂立方烷笼状骨架(亦即选择立方烷为参考物)的等键反应设计,精确计算了9种多氰基立方烷的生成热;发现B3LYP/6-31G^*结果分别地均与上述四种半经验MO方法求得的生成热之间存在良好的线性关系(相关系数均在0.9994以上),且以AM1生成热与B3LYP/6-31G^*计算值最为接近,其次,其它12种多氰基立方烷的精确生成热借助上述线性关系通过校正对应的AM1结果而获得,多氰基立方烷的生成热很高,且随-CN基数目的增加而线性地增大,表明它们属于极具潜力的“新一低高能炸药”而具开发价值。     

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

在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程序关联静态和动态理论指标,表明它们均可以用作预示标题物的热解引发机理和撞击感度.     

多硝基甲烷Mannich碱稳定性的理论研究

用MNDO SCF-MO方法全优化了系列多硝基甲烷Mannich碱的几何构型, 计算了它们的电子结构。在胺、醛组分相同时, 标题物>N-CH2-Y的特征键CH2-Y的键级随酸组分亲核性的增强而增大。该CH2-Y键级较小是造成Mannich碱在溶液中不稳定的重要原因。C-NO2键的键级在分子中较小, 可能是热或光解等受激分解的引发键, 从电子结构特征上阐明了以硝仿为酸组分的Mannich碱稳定性较差的原因。     

硝酸酯分子几何构型的量子化学研究

运用MINDO / 3、MNDO 和AM1 三种半经验分子轨道(MO)方法, 通过SCF计算, 首次系统地获得了32个硝酸酯化合物分子的全优化几何构型。三种方法的计算结果与已报道的四个化合物(硝酸甲酯、吉纳、硝化甘油和太安)的实验结果相比, AM1法较好。所有硝酸酯的酯基(-ONO~2)具有近似不变的几何参数。直链烷基硝酸酯的键长和键角极为相近, 全部重原子均共平面。二元直链和四元硝酸酯具有对称的分子构型。     

双（2,2-二硝基丙基）甲缩醛的结构和爆轰性能的量子化学研究

利用B3LYP/6-311＋G（2d,p）方法对一种新型含能增塑剂双（2,2-二硝基丙基）甲缩醛进行几何优化,计算了其红外光谱、生成焓和爆轰特性. 分析了最弱键的键离解能和键级并预测了目标化合物的热稳定性. 结果表明双（2,2-二硝基丙基）甲缩醛中的四个N-NO2键的键离解能都为164.38 kJ/mol. 表明目标化合物是一个热力学性能稳定的化合物. 以凝聚相生成焓和分子密度为基础,采用Kamlet-Jacobs方法预测其爆速和爆压. 目标化合物的晶体结构属于P21空间群.     

菁染料光解机理的前线分子轨道研究

通过SCF PM3和AM1 MO计算研究，揭示了正离子型菁染料与相应中性化合物在分子几何、电荷分布、前线分子轨道能级和组成等方面的不同特点。特别是根据前线轨道理论和成键三原则，通过比较氧的激发单态（^1O2^*)或超氧负离子自由基（O2^-.）与标题物之间的前线轨道相互作用，阐明了标题物的光解机理和正离子较中性化合物光稳定的根由。结论与实验事实良好相符。     

多硝基尿酸衍生物含能性质的理论研究

为了寻找兼具优异爆轰性能和良好热力学及动力学稳定性的高能材料, 本文设计了15个硝基尿酸化合物, 运用密度泛函理论, 对其性质进行了研究. 通过半经验的K-J方程和比冲量预测了其爆炸性能, 结果表明, 所设计分子的爆热、 分子密度、 爆炸速率和爆炸压强同硝基取代基数目之间存在较强的线性关系. 三硝基尿酸和四硝基尿酸衍生物的爆炸速率超过了8.0 km/s, 爆炸压强超过了30 GPa, 并且大多数衍生物的比冲量要高于目前经常使用的炸药黑索金. 通过计算N—NO₂键的解离能、 特征落高、 分子的自由空间预判了衍生物的稳定性和撞击感度, 结果显示, 绝大多数分子有大于80 kJ/mol的键解离能. 本文的理论结果可以为实验上设计合成新的高能材料提供一些有用的信息.     

Theoretical prediction on heats of formation for polyisocyanocubanes

The heats of formation (HOF) for all the 21 polyisocyanocubanes are calculated systematically with density functional theory (DFT) B3LYP and semiempirical MO(MINDO/3, MNDO, AM1 and PM3) methods. First, the accurate HOFs for the 8 title compounds are obtained by means of designed isodesmic reactions at DFT-B3LYP/6-31G* level, and the cubane cage skeleton has not been broken (i.e. choosing cubane as a reference compound) to produce more accurate and reliable results. It is found that there are good linear relationships between the HOFs calculated using the B3LYP/6-31G* and four semiempirical MO methods, respectively, and all of the linear correlation coefficients are more than 0.9971. The HOFs obtained from PM3 calculation are the best among the four semiempirical MO methods. Then, the accurate HOFs at B3LYP/6-31G* level of other 13 polyisocyanocubanes are obtained by systematically correcting their PM3-calculated HOFs. Polyisocyanocubanes have very high HOFs, and the HOFs increase linearly with the increasing of the number of isocyano groups in a molecule. The results show that polyisocyanocubanes are the new generation explosives with highly potential and exploitable value.     

Theoretical studies on heats of formation for cubylnitrates using density functional theory B3LYP method and semiempirical MO methods

The heats of formation (HOF) have been calculated for all the 21 cubylnitrate compounds using the semiemprical molecular orbital (MO) methods (MINDO/3, MNDO, AM1, and PM3) and for 8 of 21 cubylnitrates containing 1–4 ? ONO₂ groups using the density functional theory (DFT) method at the B3LYP/6‐31G* level by means of designed isodesmic reactions. The cubane cage skeletons in cubylnitrate molecules have been kept in setting up isodesmic reactions to produce more accurate and reliable results. It is found that there are good linear relationships between the HOFs of the 8 cubylnitrates calculated using B3LYP/6‐31G* and two semiempirical MO (PM3 and AM1) methods, and the linear correlation coefficients of PM3 and AM1 methods are 0.9901 and 0.9826, respectively. Subsequently, the accurate HOFs at B3LYP/6‐31G* level of other 13 cubylnitrates containing 4–8 ? ONO₂ groups are obtained by systematically correcting their PM3‐calculated HOFs. Compared with noncaged nitrates, all the 21 cubylnitrates have high heats of formation implying that they may be very powerful energetic materials and have highly exploitable value. The relationship between the HOFs and the molecular structures of cubylnitrates has been discussed. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Theoretical studies on the heats of formation and the interactions among the difluoroamino groups in polydifluoroaminocubanes

The heats of formation (HOFs) were calculated for a series of polydifluoroaminocubanes by using density functional theory (DFT), Hartree-Fock, and MP2 method with 6-31G basis set as well as semiempirical methods. The cubane skeleton was not broken in the process of designing isodesmic reactions; i.e., the cubane skeleton was chosen for a reference compound. The contribution of difluoroamino group to the heat of formation deviates from group additivity. The semiempirical MO (MNDO, AM1, and PM3) methods did not produce accurate and reliable results for the HOFs of the title compounds. The relationship between HOFs and molecular structures was discussed. It was found that the HOFs decreased dramatically initially and then gradually with each difluoroamino group attached to the cubane skeleton. The distance between difluoroamino groups influences the values of HOFs. The interacting energies of polydifluoroaminocubanes are in the range 14-20 kJ/mol. The interaction of neighbor difluoroamino groups discords with the group additivity. The average interaction energy between the nearest-neighbor NF(2) group in the most stable conformer of octadifluoroaminocubane is 13.94 kJ/mol at the B3LYP/6-31G level. The NF(2) group can rotate freely around the C-N bond. The relative stability of the title compounds was accessed on the basis of the calculated HOFs, the energy gaps between the frontier orbitals, and the bond order of C-NF(2). These results provide basic information for the molecular design of novel high energetic density materials.     

Theoretical prediction on heats of formation for polyisocyanocubanes——Looking for typical high energetic density material(HEDM)

The heats of formation (HOP) for all the 21 polyisocyanocubanes are calculated systematically with density functional theory (DFT) B3LYP and semiempirical MO(MINDO/3, MNDO, AM1 and PM3) methods. First, the accurate HOFs for the 8 title compounds are obtained by means of designed isodesmic reactions at DFT-B3LYP/6-31G* level, and the cubane cage skeleton has not been broken (i.e. choosing cubane as a reference compound) to produce more accurate and reliable results. It is found that there are good linear relationships between the HOFs calculated using the B3LYP/6-31G* and four semiempirical MO methods, respectively, and all of the linear correlation coefficients are more than 0.9971. The HOFs obtained from PM3 calculation are the best among the four semiempirical MO methods. Then, the accurate HOFs at B3LYP/6-31G* level of other 13 polyisocyanocubanes are obtained by systematically correcting their PM3-calculated HOFs. Polyisocyanocubanes have very high HOFs, and the HOFs increase linearly with the increa     

Substituent effects on heats of formation, group interactions, and detonation properties of polyazidocubanes

We have calculated the heats of formation (HOFs) for a series of polyazidocubanes by using the density functional theory (DFT), Hartree-Fock, and MP2 methods with 6-31G* basis set as well as semiempirical methods. The cubane skeleton was chosen for a reference compound, that is, the cubane skeleton was not broken in the process of designing isodesmic reactions. There exists group additivity for the HOF with respect to the azido group. The semiempirical AM1 method also produced reliable results for the HOFs of the title compounds, but the semiempirical MINDO3 did not. The relationship between HOFs and molecular structures was discussed. It was found that the HOF increases 330-360 kJ/mol for each additional number of the azido group being added to the cubane skeleton. The distance between azido groups slightly influences the values of HOFs. The interacting energies of neighbor azido groups in polyazidocubanes are in the range of 2.3 approximately 6.6 kJ/mol, which are so small and less related to the substituent numbers. The average interaction energy between nearest neighbor --N3 groups in the most stable conformer of octaazidocubane is 2.29 kJ/mol at the B3LYP/6-31G* level. The relative stability related to the number of azido groups of the title compounds was assessed based on the calculated HOFs, the energy gaps between the frontier orbitals, and the bond orders of the C--N3 and C--C bonds. The predicted detonation velocity of hepta- and octa-derivatives is over 9 km/s, and the detonation pressure of them is ca. 40 GPa or over.     

Computational investigation of the heat of formation, detonation properties of furazan-based energetic materials

The heats of formation (HOFs) for a series of furazan-based energetic materials were calculated by density functional theory. The isodesmic reaction method was employed to estimate the HOFs. The result shows that the introductions of azo and azoxy groups can increase the HOF, but the introduction of azo group can increase the more HOF, when compared with azoxy group. The detonation velocities and detonation pressures of the furazan-based energetic materials are further evaluated at B3LYP/6-31G* level. Dioxoazotetrafurazan and azoxytetrafurazan may be regarded as the potential candidates of high-energy density materials because of good detonation performance. In addition, there are good linear correlations between OB and detonation velocities, and OB and detonation pressures. The energy gaps between the HOMO and LUMO of the studied compounds are also investigated. These results provide basic information for the molecular design of novel high-energy density materials.     

叠氮化合物C6H6-n(N3)n(n=1~6)的密度泛函理论研究

运用密度泛函理论, 在B3LYP/6-31G*水平上, 对叠氮化合物C6H6-n(N3)n(n=1~6)进行理论计算, 并对所得到的几何结构进行了振动频率分析. 计算结果表明, 这些化合物是热力学稳定的. 基于自然键轨道理论, 分析了稳定结构的电荷分布及成键情况. 在不破坏苯环和叠氮基的原则下, 设计等键反应计算了这些化合物的生成热, 结果表明, 标题化合物的生成热都很高, 且随着叠氮基数目的增加而线性增大.     

Theoretical studies on four-membered ring compounds with NF2, ONO2, N3, and NO2 groups

Density functional theory (DFT) method has been employed to study the geometric and electronic structures of a series of four-membered ring compounds at the B3LYP/6-311G** and the B3P86/6-311G** levels. In the isodesmic reactions designed for the computation of heats of formation (HOFs), 3,3-dimethyl-oxetane, azetidine, and cyclobutane were chosen as reference compounds. The HOFs for N(3) substituted derivations are larger than those of oxetane compounds with --ONO2 and/or --NF2 substituent groups. The HOFs for oxetane with --ONO2 and/or --NF2 substituent groups are negative, while the HOFs for N3 substituted derivations are positive. For azetidine compounds, the substituent groups within the azetidine ring affect the HOFs, which increase as the difluoroamino group being replaced by the nitro group. The magnitudes of intramolecular group interactions were predicted through the disproportionation energies. The strain energy (SE) for the title compounds has been calculated using homodesmotic reactions. For azetidine compounds, the NF2 group connecting N atom in the ring decrease the SE of title compounds. Thermal stability were evaluated via bond dissociation energies (BDE) at the UB3LYP/6-311G** level. For the oxetane compounds, the O--NO2 bond is easier to break than that of the ring C--C bond. For the azetidine and cyclobutane compounds, the homolyses of C--NX2 and/or N--NX2 (X = O, F) bonds are primary step for bond dissociation. Detonation properties of the title compounds were evaluated by using the Kamlet-Jacobs equation based on the calculated densities and HOFs. It is found that 1,1-dinitro-3,3-bis(difluoroamino)-cyclobutane, with predicted density of ca. 1.9 g/cm(3), detonation velocity (D) over 9 km/s, and detonation pressure (P) of 41 GPa that are lager than those of TNAZ, is expected to be a novel candidate of high energy density materials (HEDMs). The detonation data of nitro-BDFAA and TNCB are also close to the requirements for HEDMs.     

Theoretical Studies on Energetic Property and Stability of Pyrazine and Pyridine Derivatives

Density functional calculations at the B3LYP level with 6‐311G** and aug‐cc‐pVDZ basis sets were performed to predict the heats of formation (HOFs) for two pyrazine derivatives and eight pyridine derivatives. In the isodesmic reactions designed for the computation of heats of formation (HOFs), pyrazine and pyridine were chosen as reference compounds. The N‐oxidations for the ring nitrogen of pyrazine and pyridine derivatives decrease the HOF values when N‐oxide oxygen is neighboring with amino groups, but increase when it neighbors with nitro groups. Thermal stability was evaluated via bond dissociation energies (BDE) at the UB3LYP/6‐311G** level. As a whole, the homolysis of C–NO₂ bonds is the main step for bond dissociation of the title compounds. The BDE values of title compounds are influenced by intramolecular hydrogen bonds. The hydrogen bond effects associated with the length of the H···O bonds were analyzed by the electron density at the critical points and natural bond orbital.