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.     

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 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     

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.     

多氰基立方烷生成热的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水平下，通过不破裂金刚烷分子骨架（即 选择金刚作为参考物）的等键反应设计，比较精确地计算了系列多硝基金金刚烷的 生成热。经验性基团加和法和半经验MO法（AM1，PM3，MNDO，MNDO/3）均不适用于 标题生成熟的估算。4种半经验MO方法中，以MP3计算结果略好些。探讨了生成热和 分子结构的关系，发现桥头C上硝基使生成热减小，而偕二硝基使生成热增大。运 用生成热、前沿轨道能级差和C-NO_2键级等计算结果，阐明了标题的相对稳定性， 为新一代高能量密度材料（HEDM）的分子设计提供了基础数据和规律性。     

Theoretical Study of Heterocyclic Nitro Compounds for Use as Novel Explosives and Propellants

The heats of formation (HOFs) of heterocyclic nitro compounds were obtained by using a density functional theory B3LYP method with 6‐31G* and 6‐311+G** basis sets. The isodesmic reactions designed for the evaluation of HOFs keep most of the basic ring structures of the title compounds and thus ensure the credibility of the results. The values of HOFs are 567.90, 874.29 and 975.83 kJ/mol at the B3LYP/6‐31G* level for hexanitrohexazaadamantane ( A ), nonanitrononaza‐tetracyclo[7.3.1.1^3,7.1^5,11] pentadecane ( B ) and tetranitrotetrazacubane ( C ) respectively. The predicted detonation velocities of the title compounds are larger than, and detonation pressures are much larger than that of the widely used 1,3,5,7‐tetranitro‐1,3,5,7‐tetraazacyclooctane (HMX). The dissociation energy for the weakest C‐N bonds in the cage skeleton of the title compounds are 137‐144 kJ/mol at the B3LYP/6‐31G* level.     

Computational DFT studies on a series of toluene derivatives as potential high energy density compounds

Based on the full optimized molecular geometric structures at B3LYP/6-31G**, B3LYP/6-31+G**, B3P86/6-31G**, and B3P86/6-31+G** levels, the densities (ρ), detonation velocities (D), and pressures (P) for a series of toluene derivatives, as well as their thermal stabilities, were investigated to look for high energy density compounds (HEDCs). The heats of formation (HOFs) are also calculated via designed isodesmic reactions. The calculations on the bond dissociation energies (BDEs) indicate that the BDEs of the initial scission step are between 48 and 59 kcal/mol, and pentanitrotoluene is the most reactive compound, while 2,4,6-trinitrotoluene is the least reactive compound for toluene derivatives studied. A good linear relationship between BDE/E and impact sensitivity is also obtained. The condensed phase HOFs are calculated for the title compounds. These results would provide basic information for the further studies of HEDCs. The detonation data of pentanitrotoluene show that it meets the requirement for HEDCs.     

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     

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

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

Structural changes, P-P bond energies, and homolytic dissociation enthalpies of substituted diphosphines from quantum mechanical calculations

The molecular structures of the diphosphines P(2)[CH(SiH(3))(2)](4), P(2)[C(SiH(3))(3)](4), P(2)[SiH(CH(3))(2)](4), and P(2)[Si(CH(3))(3)](4) and the corresponding radicals P[CH(SiH(3))(2)](2), P[C(SiH(3))(3)](2), P[SiH(CH(3))(2)](2), and P[Si(CH(3))(3)](2) were predicted by theoretical quantum chemical calculations at the HF/3-21G*, B3LYP/3-21G*, and MP2/6-31+G* levels. The conformational analyses of all structures found the gauche conformers of the diphosphines with C(2) symmetry to be the most stable. The most stable conformers of the phosphido radicals were also found to possess C(2) symmetry. The structural changes upon dissociation allow the release of some of the energy stored in the substituents and therefore contribute to the decrease of the P-P bond dissociation energy. The P-P bond dissociation enthalpies at 298 K in the compounds studied were calculated to vary from -11.4 kJ mol(-1) (P(2)[C(SiH(3))(3)](4)) to 179.0 kJ mol(-1) (P(2)[SiH(CH(3))(2)](4)) at the B3LYP/3-21G* level. The MP2/6-31+G* calculations predict them to be in the range of 52.8-207.9 kJ mol(-1). All the values are corrected for basis set superposition error. The P-P bond energy defined by applying a mechanical analogy of the flexible substituents connected by a spring shows less variation, between 191.3 and 222.6 kJ mol(-1) at the B3LYP/3-21G level and between 225.6 and 290.4 kJ mol(-1) at the MP2/6-31+G* level. Its average value can be used to estimate bond dissociation energies from the energetics of structural relaxation.     

Strain Energy Calculations of Caged Silanes

Abstract

Strain energies (SE) of typical caged silanes were evaluated at the B3LYP/6-31+G** and B3LYP/aug-cc-pVDZ levels, in combination with the isodesmic and homodesmotic reactions. The SE values of Si₄H₄ and Si₈H₈ are 512 kJ/mol and 336–338 kJ/mol, respectively, at the B3LYP/6-31+G** level, whereas the SE values of the caged silanes with five-number or larger rings are very small. In comparison, the SE value of Si₈H₈ is much smaller than 656–707 kJ/mol of cubane (C₈H₈), since the Si?Si bond length (2.274 Å) of Si₈H₈ is much larger than the corresponding C?C bond length of cubane. The large electron charge density at body-center of Si₄H₄ causes strong repulsion between the cage center and Si?Si bond, which leads to the Si?Si bond bending. The electron charge density at the face center of Si₄H₄ is more than two times of that of Si₈H₈, which also contributes significantly to the large SE difference between Si₄H₄ and Si₈H₈.     

Calorimetric and computational study of thiacyclohexane 1-oxide and thiacyclohexane 1,1-dioxide (thiane sulfoxide and thiane sulfone). Enthalpies of formation and the energy of the S=O bond

A rotating-bomb combustion calorimeter specifically designed for the study of sulfur-containing compounds [J. Chem. Thermodyn. 1999, 31, 635] has been used for the determination of the enthalpy of formation of thiane sulfone, 4, Delta(f)H(o) m(g) = -394.8 +/- 1.5 kJ x mol(-1). This value stands in stark contrast with the enthalpy of formation reported for thiane itself, Delta(f)H(o) m(g) = -63.5 +/- 1.0 kJ x mol(-1), and gives evidence of the increased electronegativity of the sulfur atom in the sulfonyl group, which leads to significantly stronger C-SO2 bonds. Given the known enthalpy of formation of atomic oxygen in the gas phase, Delta(f)H(o) m(O,g) = +249.18 kJ x mol(-1), and the reported bond dissociation energy for the S=O bond in alkyl sulfones, BDE(S=O) = +470.0 kJ x mol(-1), it was possible to estimate the enthalpy of formation of thiane sulfoxide, 5, a hygroscopic compound not easy to use in experimental calorimetric measurements, Delta(f)H(o) m(5) = -174.0 kJ x mol(-1). The experimental enthalpy of formation of both 4 and 5 were closely reproduced by theoretical calculations at the G2(MP2)+ level, Delta(f)H(o) m(4) = -395.0 kJ x mol(-1) and Delta(f)H(o) m(5) = -178.0 kJ x mol(-1). Finally, calculated G2(MP2)+ values for the bond dissociation energy of the S=O bond in cyclic sulfoxide 5 and sulfone 4 are +363.7 and +466.2 kJ x mol(-1), respectively.     

CX_2(X=H,F,Cl)与甲基异丙基醚C-H键插入反应的理论研究

用从头计算方法在MP2 /6 31G(d)水平上研究了CX2 (X =H ,F ,Cl)与甲基异丙基醚的C -H键插入反应。CCl2 与甲基异丙基醚两个不同的α C的C -H键插入势垒分别为 117.2kJ/mol (甲基 )和 2 0 .6kJ/mol (异丙基 )。CF2 与异丙基α C的C -H键上插入势垒为 12 0 .0kJ/mol,在插入甲基上C -H键时会引起C -O键的断裂。CH2 的插入反应则不需要势垒。对CX2 与二甲醚、甲乙醚、甲基异丙基醚、甲基苄基醚上各种不同的C -H键插入势垒进行了比较 ,甲基和苯基都促使其毗邻的C -H键更容易被CX2 所插入     

Computational evidence for methyl-donated hydrogen bonds and hydrogen-bond networking in 1,2-ethanediol-dimethyl sulfoxide

The 1:1 complex of 1,2-ethanediol with dimethyl sulfoxide was studied using density functional theory. A network of three hydrogen bonds holds the complex together, including two in which each methyl group donates to the same hydroxyl oxygen. Four lines of evidence support the existence of methyl-donated hydrogen bonds. The interaction energy is 36 +/- 5 kJ/mol using Becke's three parameter hybrid theory with the 1991 nonlocal correlation functional of Perdew and Wang, and a moderately large basis set (B3PW91/6-311++G**//B3PW91/6-31+G**). To determine the energy of each hydrogen bond, a relaxed potential energy scan was performed in a smaller basis set to break the weaker hydrogen bonds by forced systematic rotation of the methyl groups. Two cross-checking analyses show cooperative effects that cause individual hydrogen bond energies in the network to be nonadditive. When one methyl hydrogen bond is broken, the remaining interactions stabilize the complex by storing an additional 2-3 kJ/mol. With all hydrogen bonds intact, the O[bond]H...O[bond]S hydrogen bond contributes 26 +/- 2 kJ/mol stability, and each weak methyl bond stores 5 +/- 2 kJ/mol.     

Structural and spectral characteristics of the electrogenerated tetracyanoethylene dianion

Structural and spectral characteristics of the electrogenerated tetracyanoethylene dianion (TCNE2-) were experimentally and theoretically examined. Spectroelectrochemistry of TCNE gives the spectra of TCNE2- in CH3CN at 220 nm, and in CH2Cl2 at 300 nm. These spectral characteristics are well explained by CIS/6-31G(d) and semiempirical CNDO/S-CI calculations. The bands in CH3CN and in CH2Cl2 are assigned to the degenerate 1E<--1A1 transition at the D2d structure and the 1B2u<--1Ag transition at the D2h structure, respectively. The rotation barrier of the C=C bond in TCNE2- is estimated by Hartree-Fock (HF), second-order M?ller-Plesset perturbation (MP2) and fourth-order MP (MP4) calculations with 6-31G(d), 6-31+G(d) and 6-311+G(d) basis sets as 42-51 kJ mol(-1). The D2d structure is most stable, and the D2h structure represents the transition state of the internal rotation. The calculations reveal that the two-electron addition to the antibonding LUMO of TCNE causes an easy rotation around the C=C bond of TCNE2- characterized by the formal single bond. These results show that TCNE2- preferentially adopts D2d and D2h structures in solvents depending upon the solvent nature by virtue of the easy rotation around the C=C bond.     

Energetics of the N-O bonds in 2-hydroxyphenazine-di-N-oxide

The standard enthalpy of formation and the enthalpy of sublimation of crystalline 2-hydroxyphenazine-di-N-oxide, at T = 298.15 K, were determined from isoperibol static bomb combustion calorimetry and from Knudsen effusion experiments, as -76.7 +/- 4.2 kJ.mol(-1) and 197 +/- 5 kJ.mol(-1), respectively. The sum of these two quantities gives the standard enthalpy of formation in the gas-phase for this compound, delta(f)H(m)degrees(g) = 120 +/- 6 kJ.mol(-1). This value was combined with the gas-phase standard enthalpy of formation for 2-hydroxyphenazine retrieved from a group estimative method yielding the mean (N-O) bond dissociation enthalpy, in the gas-phase, for 2-hydroxyphenazine-di-N-oxide. The result obtained with this strategy is (DH(m)degrees (N - O)) = 263 +/- 4 kJ.mol(-1), which is in excellent agreement with the B3LYP/6-311+G(2d,2p)//B3LYP/6-31G(d) computed value, 265 kJ.mol(-1).     

炸药敏化剂的分子结构、电子结构和生成热的理论研究

在HF/6-31G*和B3LYP/6-31G*基组水平上, 对硝酸乙酯(EN)、硝酸正丙酯(NPN)、硝酸异丙酯(IPN)、硝酸异辛酯(EHN)和四甘醇二硝酸酯(TEGDN)五种炸药敏化剂进行理论计算, 研究了标题物的分子结构、电子结构和能量等方面的性质. 基于Mulliken布居和键长分析, 五种硝酸酯分子的热分解始于O2—N3键的断裂, 且由Mulliken电荷分布推知分子热解产生NO2气体. 在分析前线轨道能(EHOMO, ELUMO)和能量差(ΔE)的基础上对五种硝酸酯的相对热稳定性大小进行了评估. 由等键反应获得的EN、IPN、NPN、EHN和TEGDN的标准生成热分别是-155.972、-190.896、-175.279、-272.376和-790.733 kJ·mol-1.     

Clar goblet and related non-Kekulé benzenoid LPAHs. A theoretical study

The results of an ab initio and semiempirical study of Clar Goblet (1), a C(38)H(18) non-Kekulé diradical LPAH, and its constitutional isomers 4 and 5 are reported. Planar D(2)(h)-1 was only 87.4 (triplet) and 83.8 (singlet) kJ/mol less stable than its planar Kekulé isomer C(2)(v)-6 (at (U)B3LYP/6-31G). Planar C(s)-4 was 63.6 (triplet) and 76.5 (singlet) kJ/mol less stable than 6. Overcrowded C(1)-5 was 80.1 (triplet) and 98.1 (singlet) kJ/mol less stable than 6. In concealed non-Kekulé 1, the singlet was more stable then the triplet by 3.6 kJ/mol, while in nonconcealed non-Kekulé 4 and 5, the triplets were more stable than the corresponding singlets by 12.9 and 18.1 kJ/mol, respectively, in accordance with theory. The spin density in 1, 4, and 5 is delocalized throughout the positions corresponding to active peri-peri coupling positions of the radical anion of naphthanthrone (2). The bond lengths in 1, 4, and 5 are in the range expected for aromatic compounds, except for the central carbon-carbon bonds, which are considerably elongated. A certain stabilization is evident in the homodesmotic reaction singlet-1 + 10 + 10 --> 11 + 3 + 3, indicating a "communication" between the two benzo[cd]pyrenyl radical (3) units of diradical 1. The HOMA indices indicate that in both singlet 1 and triplet 1 all of the rings except the central one have a significant aromatic character. The central ring is essentially antiaromatic, having negative HOMA index (-0.140 at UB3LYP/6-31+G). The stabilities of 1(2)(-) and 1(2+) are decreased relative to 3(-) and 3(+), respectively.     

DFT Study of the Group Interactions in 1,3,5‐Triamino‐2,4,6‐Trinitrobenzene and 1,3,5‐Triamino‐2,4,6‐Tridifluoroaminobenzene

Density functional calculations on isodesmic disproportionation reactions of 1,3,5‐triamino‐2,4,6‐trinitrobenzene (TATB) and 1,3,5‐triamino‐2,4,6‐tridifluoroaminobenzene (TATDB) indicate that the interaction between nitro groups on meta carbons of TATB, which brings about unstability to the molecule, is surprisingly larger than that between difluroamino groups in TATDB. The electron‐withdrawing and electron‐donating groups generate large positive and very small negative values of E_{disproportion}, respectively. When both electron‐withdrawing and electron‐donating groups are attached to the benzene skeleton at the same time, large negative disproportionation energy is produced, which stabilizes the derivatives. The values of E_{disproportion} for TATB and TATDB are predicted to be ‐48.03 kJ/mol and ‐63.54 kJ/mol, respectively, indicating that the total interaction among groups with stabilization effects in TATDB is larger than that in TATB. The large difference of the E_{disproportion} values between TATB and TATDB is derived from the large difference between the interactions of the meta‐nitro group and those of meta‐difluoroamino groups. The energy barriers for the C‐N internal rotation of NO₂ group and NF₂ groups are 74.7 kJ/mol and 185.8 kJ/mol for TATB and TATDB, respectively. The large energy barrier for the rotation of the NF₂ group is caused by its stabilization interaction with neighbor amino groups, instead of steric effects. When the number of pairs of amino‐nitro or amino‐difluoroamino groups increases, there are more negative charges on the NO₂/NF₂ groups and on the O/F atoms.