Comparative theoretical studies of high energetic cyclic nitramines

Density functional theory studies on cyclic nitramines were performed at B3LYP/6‐311G(d,p) level. The crystal structures were obtained by molecular packing calculations. Heats of formation (HOFs) were predicted through designed isodesmic reactions. Results indicate that the value of HOF relates to the number of =N–NO₂ group and aza nitrogen atom and increases with the augment of the number of =N–NO₂ group and aza nitrogen atom for cyclic nitramines. Detonation performance was evaluated by using the Kamlet–Jacobs equations based on the calculated densities and HOFs. All the cyclic nitramines exhibit better detonation performance than 1,3,5‐trinitro‐1,3,5‐triazacyclohexane and 1,3,5,7‐tetranitro‐1,3,5,7‐tetraazacyclooctane. The stability of cyclic nitramines was investigated by the bond dissociation energies. The result shows that the increase of =N?NO₂ group or aza nitrogen atom reduces the stability of the title compounds. These results provide basic information for molecular design of novel high energetic density materials. Copyright © 2013 John Wiley & Sons, Ltd.     

Computational study of pyrazine‐based derivatives and their N‐oxides as high energy materials

Gas‐phase heats of formation (HOF), solid‐phase HOF, detonation properties, electronic structure and thermal stability for a series of polynitro pyrazine derivatives containing three heterocycles have been investigated using density functional theory. It is found that the nitro group is an efficient tool to improve HOF of pyrazine derivatives. Furthermore, detonation velocities and detonation pressures of these compounds are evaluated using empirical Kamlet–Jacobs equations. As a result, it indicates that the nitro group is useful to enhance detonation properties. Detonation velocities of five compounds are 9.67, 9.20, 9.74, 9.76 and 9.87 km/s, respectively, which are significantly larger than that of HMX (9.10 km/s). Bond dissociation energy is also performed to investigate their thermal stability, showing that thermal stability of these compounds is little affected by nitro groups or the position of substituent groups. Considering solid‐phase HOF, detonation properties and thermal stability, some of pyrazine derivatives can be potential high energy density materials. Copyright © 2013 John Wiley & Sons, Ltd.     

Theoretical studies on a series of nitroaliphatic energetic compounds

Density functional theory calculations at the B3LYP/6-311G＊＊ level are performed to study the geometric and elec- tronic structures of a series of nitroaliphatic compounds. The heats of formation （HOF） are predicted through the designed isodesmic reactions. Thermal stabilities are evaluated via the homolytic bond dissociation energies （BDEs）. Further, the correlation is developed between impact sensitivity h50% and the ratio （BDE/E） of the weakest BDE to the total energy E containing zero point energy correction. In addition, the relative stability of the title compounds is evaluated based on the analysis of calculated Mulliken population and the energy gaps between the frontier orbitals. The calculated BDEs, HOFs, and energy gaps consistently indicate that compound 1,1,1,6,6,6-hexanitro-3-hexyne is the most unstable and the compound 3,3,4,4,-tetranitro-hexane is the most stable. These results provide basic information for the molecular design of novel high energetic density materials.     

Theoretical study of energetic trinitromethyl‐substituted tetrazole and tetrazine derivatives

Density functional theory method was used to study the heats of formation, energetic properties, and thermal stability for a series of trinitromethyl‐substituted tetrazole and tetrazine derivatives with different substituents. It is found that the group ―NO₂, ―NHNO₂, or ―NF₂ play a very important role in increasing the heats of formation of the derivatives. The calculated detonation velocities and pressures indicate that the group ―CF₂NF₂, ―NHNO₂, ―1H‐tetrazolyl, ―2H‐tetrazolyl, or ―1,2,4,5‐tetrazinyl is an effective structural unit for enhancing their detonation performance. An analysis of the bond dissociation energies for several relatively weak bonds indicates that incorporating the group ―NHNO₂ and ―NH₂ into parent ring decreases their thermal stability. Considering the detonation performance and thermal stability, 37 compounds may be considered as the potential high‐energy compounds. Their oxygen balances are close to zero. These results provide basic information for the molecular design of novel high‐energy compounds. Copyright © 2013 John Wiley & Sons, Ltd.     

聚对苯二甲酸丁二醇酯二聚体键离能的理论研究

采用密度泛函理论方法对聚对苯二甲酸丁二醇酯(PBT)二聚体的键离能进行了计算.为了选取较为精确的方法来计算PBT各个键的键离能,以与PBT具有相同的酯基官能团的乙酸乙酯为模型参照物.采用M062X, B3P86, M06, PBE0, wB97xD方法分别在基组6-31G(d), 6-311G(d), 6-311+G(d, p), 6-311++G(d, p), cc-pVDZ, cc-pVTZ水平下对乙酸乙酯的键离能进行计算.通过对比计算结果与iBonD数据库的乙酸乙酯实验测定值可知,M062X在基组6-311G(d)水平下计算结果与实验值最为接近.因此,本研究采用M062X方法在基组6-311G(d)水平下对聚对苯二甲酸丁二醇酯(PBT)二聚体的键离能进行计算.计算结果表明：在PBT的各键中C-C_arcmatic键的键离能最大,主链上的C-C键离能最小,为370.9 kJ/mol.其次就是C-O键,为404.6 kJ/mol.基于PBT键离能的计算结果,设计了3条PBT二聚体热降解过程可能形成的反应路径,分析了热解产物的形成机理.结果表明PBT二聚体热解过程可...     

含能双环HMX衍生物分子设计的密度泛函研究

运用DFT-B3LYP/6-311G(d,p)方法,计算了所设计的三种双环HMX（2,4,6,8-四硝基-2,4,6,8-四氮杂双环[3.3.0]辛烷）衍生物分子。基于理论晶体密度和固态生成热计算衍生物分子的爆轰性能;通过前线轨道能与特征高度（h50）评价衍生物分子的感度。结果表明,理论晶体密度均高于1.90 g.cm-3,爆速高于9.0 km.s-1,爆压约为40 GPa。三种双环HMX衍生物分子是潜在的高能量密度材料。     

Computational design of tetrazolone‐based high‐energy density energetic materials: Property prediction and decomposition mechanism

The density functional theory methods are used to design a series of new highly energetic tetrazolone‐based molecules by the combination of the linked tetrazolone framework and versatile substitutes. The molecular and electronic structures, physicochemical, and energetic properties were analyzed and predicted. The decomposition mechanisms were computationally simulated, and 3 potential decomposition channels were proposed. These newly designed tetrazolone‐based compounds show high densities (up to 2.08 g/cm³) and highly positive heats of formation (407.0‐1377.9 kJ/mol) due to all right content of nitrogen and oxygen. Most of them exhibit good detonation velocity (8.31‐9.62 km/s) and detonation pressure (32.40‐43.86 GPa), and some are comparative to excellent explosive CL‐20. Results show that compounds 6 , 10 , 11 , 12 , 15 , 16 , 17 , 22 , 23 , and 24 own superior detonation performance than widely used explosive HMX and may be promising candidates of green high‐performance energetic materials.     

DFT studies on a high-energy density cage compound 1, 3, 5, 7, 9, 11-hexo(N(CH3)NO2)-2, 4, 6, 8, 10, 12-hexaazatetracyclo[5, 5, 0, 0, 0] dodecane

The infrared and Raman spectra, heat of formation (HOF) and thermodynamic properties were investigated by B3LYP/6-31G^** method for a new designed polynitro cage compound 1,3,5,7,9,11-hexo(N(CH₃)NO₂)-2,4,6,8,10,12-hexaazatetracyclo[5,5,0,0,0]dodecane. The detonation velocity (D) and pressure (P) were predicted by the Kamlet–Jacobs equations based on the theoretical density and condensed HOF. The bond dissociation energies and bond orders for the weakest bonds were analysed to investigate the thermal stability of the title compound. The computational result shows that the detonation velocity and pressure of the title compound are superior to those of hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX), but inferior to those of 1,3,5,7-tetranitro-1,3,5,7-tetraazacyclooctane (HMX) and hexanitrohexaazaisowurtzitane (HNIW). And the analysis of thermal stability shows that the first step of pyrolysis is the rupture of the N₇–NO₂ bond. The crystal structure obtained by molecular mechanics belongs to the P2₁ space group, with the lattice parameters Z = 2, a = 11.8246 Å, b = 10.4632 Å, c = 15.9713 Å, ρ = 1.98 g cm^?3.     

三硝基甲烷键离解能和生成焓的理论计算

采用密度泛函(DFT)四种交换/相关函数(B3LYP、B3P86、B3PW91和PBE0)结合不同的基函数,求得了三硝基甲烷C-NO2键的离解能(BDE),并且通过合理选择参考物硝基甲烷,设计等键等电子对反应,计算了气相三硝基甲烷分子的生成焓(HOF).与实验数据进行比较,PBE0/6-31g*计算出的BDE值最好,误差为-2.1 kcal mol-1;PBE0密度泛函结合带极化函数的6-31g基组得到的HOF值与实验值吻合的最好(误差在0.1 kcal mol-1以内).     

Theoretical investigation of a series of bis(1H-tetrazol-5-yl)furazan and bis(1H-tetrazol) derivatives as high-energy-density materials

Using density functional theory (DFT), a series of bis(1H-tetrazol-5-yl)furazan and bis(1H-tetrazol) derivatives with different linkages and substituents are investigated theoretically as potential high-energy-density materials (HEDMs). The heat of formation (HOF), detonation properties, natural bond orbital (NBO) and thermal stabilities are calculated and reported. The introduction of a furazan ring, an –N=N– bridge group and an –N₃ substituent is beneficial to increase the HOF of the title compounds. NBO analysis shows that there are electronic delocalisation effects among the bridge groups, furazan and tetrazole rings, and substituted groups. The conjugation effects and electronic transitions are influenced by the different linkages and substituents. The estimated detonation velocities and pressures indicate that the –ONO₂ and –NO₂ groups and the –N=N– linkage play important roles in enhancing the detonation properties. The bond dissociation energy (BDE) calculations reveal that the –NO₂ group is the substituent group which causes the least thermal stability. The bond between the substituent group and the tetrazole ring is the weakest bond in the title molecules. Considering the detonation performance and the thermal stability, 17 compounds may be promising candidates for HEDMs with good performance. Eight of them (A3, A4, C3, C4, D3, F3, G1 and G3) have better detonation properties than HMX.     

硝化纤维素生成焓、键离解能及爆轰性能的理论研究

运用密度泛函理论方法结合等键反应,研究了硝化纤维素(NC)单体的生成热及其不同位置上氧硝基键(ONO2)的键离解能,基于半经验的Kamlet-Jacobs方程估算了其爆速与爆压。结果表明:B3LYP和B3P86两种方法得到的计算结果较为接近,生成热均为负值,其大小与氧硝基的数目有关,与氧硝基的取代位置无关。理论计算得到的氧硝基键的键离解能与含氮量无关,而是与实际断裂的氧硝基键的数目相关。爆速和爆压的计算表明,NC单体中所含氧硝基的数量增多,会降低其爆速和爆压。     

Theoretical studies on structure and performance of [1,2,5]‐oxadiazolo‐[3,4‐d]‐pyridazine‐based derivatives

Based on energetic compound [1,2,5]‐oxadiazolo‐[3,4‐d]‐pyridazine, a series of functionalized derivatives were designed and first reported. Afterwards, the relationship between their structure and performance was systematically explored by density functional theory at B3LYP/6‐311 g (d, p) level. Results show that the bond dissociation energies of the weakest bond (N–O bond) vary from 157.530 to 189.411 kJ · mol^?1. The bond dissociation energies of these compounds are superior to that of HMX (N–NO_2, 154.905 kJ · mol^?1). In addition, H1, H2, H4, I2, I3, C1, C2, and D1 possess high density (1.818–1.997 g · cm^?3) and good detonation performance (detonation velocities, 8.29–9.46 km · s^?1; detonation pressures, 30.87–42.12 GPa), which may be potential explosives compared with RDX (8.81 km · s^?1, 34.47 GPa ) and HMX (9.19 km · s^?1, 38.45 GPa). Finally, allowing for the explosive performance and molecular stability, three compounds may be suggested as good potential candidates for high‐energy density materials. Copyright © 2016 John Wiley & Sons, Ltd.     

两种太根(Tri-EGDN和Tetra-EGDN)的几何结构和键裂解能的理论研究

利用密度泛函理论(DFT)中的B3LYP、B3PW91、B3P86方法研究了二缩三乙二醇二硝酸酯(Tri-EGDN)和三缩四乙二醇二硝酸酯(Tetra-EGDN)的平衡分子构型、前线轨道能(EHOMO、ELUMO)、能差(ΔE=ELUMO-EHOMO)和O-NO2键的裂解能(BDE),并分析了两种硝酸酯分子和对应自由基的热稳定性.Tri-EGDN和Tetra-EGDN中的O-NO2键的裂解能用B3P86方法计算最接近真实值.由O-NO2键的裂解能很好地符合文献中对应物质的表观活化能,推知两种物质的热分解反应只是单分子O-NO2键的均裂反应.     

Quantum chemical calculations of bond dissociation energies for COOH scission and electronic structure in some acids

Quantum chemical calculations are performed to investigate the equilibrium C-COOH bond distances and the bond dissociation energies(BDEs) for 15 acids.These compounds are studied by utilizing the hybrid density functional theory(DFT)(B3LYP,B3PW91,B3P86,PBE1PBE) and the complete basis set(CBS-Q) method in conjunction with the 6311G** basis as DFT methods have been found to have low basis sets sensitivity for small and medium molecules in our previous work.Comparisons between the computational results and the experimental values reveal that CBS-Q method,which can produce reasonable BDEs for some systems in our previous work,seems unable to predict accurate BDEs here.However,the B3P86 calculated results accord very well with the experimental values,within an average absolute error of 2.3 kcal/mol.Thus,B3P86 method is suitable for computing the reliable BDEs of C-COOH bond for carboxylic acid compounds.In addition,the energy gaps between the highest occupied molecular orbital(HOMO) and the lowest unoccupied molecular orbital(LUMO) of studied compounds are estimated,based on which the relative thermal stabilities of the studied acids are also discussed.     

Theoretical study of the interaction between pyridine derivatives and atomic chlorine. Substituent effect and nature of the bonding

The interaction of pyridine derivatives (H, 4-NH₂, 4-CH₃, 4-F, 4-CN, 4-NO₂) with atomic chlorine is investigated theoretically by the density functional theory (DFT)-based LC-BLYP/aug-cc-pVDZ method. Pyridines and Cl° are held together by a (2c--3e) bond and the intermolecular distances range from 2.313 to 2.343 Å. The existence of a N…Cl bond is confirmed by the atom-in-molecule analysis of the systems. The binding energies of the adducts, ranging from ?42.08 to ?53.96 kJ mol^?1, are linearly correlated to the proton affinity of the pyridines. The charge transfer from pyridine to Cl° varies between 0.222 and 0.277 e. The spin density analysis shows that the strongest complex has the highest (2c–3e) character. The CH bonds are contracted and the ν(CH) vibrations are blueshifted owing to the decrease in σ*(CH) occupation. This decrease results not only from the classical anomeric effect but also from the σ(CC) and σ(CN) → σ*(CH) delocalisation in the heteroaromatic ring.     

Advances and challenges in DFT-based energy materials design

The growing worldwide energy needs call for developing novel materials for energy applications. Ab initio density functional theory (DFT) calculations allow the understanding and prediction of material properties at the atomic scale, thus, play an important role in energy materials design. Due to the fast progress of computer power and development of calculation methodologies, DFT-based calculations have greatly improved their predictive power, and are now leading to a paradigm shift towards theory-driven materials design. The aim of this perspective is to introduce the advances in DFT calculations which accelerate energy materials design. We first present state-of-the-art DFT methods for accurate simulation of various key properties of energy materials. Then we show examples of how these advances lead to the discovery of new energy materials for photovoltaic, photocatalytic, thermoelectric, and battery applications. The challenges and future research directions in computational design of energy materials are highlighted at the end.     

The polarization energy of normal liquid 3He

We present new quantum Monte Carlo calculations of the ground-state energy of liquid ³He using backflow and 3-body wavefunctions and twist-averaged boundary conditions. We estimate the effects of the 3-body potential and compare to experimental measurements of the unpolarized equation of state and the magnetic susceptibility. There is still a large discrepancy with respect to the magnetic susceptibility, requiring either improved nodal surfaces or going beyond the fixed-node method.     

Theoretical investigations on 4,4′,5,5′‐tetranitro‐2,2′‐1H,1′H‐2,2′‐biimidazole derivatives as potential nitrogen‐rich high energy materials

The ―NH₂, ―NO₂, ―NHNO₂, ―C(NO₂)₃ and ―CF(NO₂)₂ substitution derivatives of 4,4′,5,5′‐tetranitro‐2,2′‐1H,1′H‐2,2′‐biimidazole were studied at B3LYP/aug‐cc‐pVDZ level of density functional theory. The crystal structures were obtained by molecular mechanics (MM) methods. Detonation properties were evaluated using Kamlet–Jacobs equations based on the calculated density and heat of formation. The thermal stability of the title compounds was investigated via the energy gaps (?E_{LUMO ? HOMO}) predicted. Results show that molecules T5 (D = 10.85 km·s^?1, P = 57.94 GPa) and T6 (D = 9.22 km·s^?1, P = 39.21 GPa) with zero or positive oxygen balance are excellent candidates for high energy density oxidizers (HEDOs). All of them appear to be potential explosives compared with the famous ones, octahydro‐1,3,5,7‐tetranitro‐1,3,5,7‐tetraazocane (HMX, D = 8.96 km·s^?1, P = 35.96 GPa) and hexanitrohexaazaisowurtzitane (CL‐20, D = 9.38 km·s^?1, P = 42.00 GPa). In addition, bond dissociation energy calculation indicates that T5 and T6 are also the most thermally stable ones among the title compounds. Copyright © 2014 John Wiley & Sons, Ltd.     

GeHn±(n=0,1,2)的势能函数和光谱常数研究

用密度泛函B3lyp/6-311++g(d,p)方法对 分子离子进行理论研究。结果表明：GeH,GeH+,GeH-1均能稳定存在,GeH+2和GeH-2有亚稳定态存在,其基态电子态分别是: 2Π(GeH),1Σ(GeH+),3Σ(GeH-),2Σ(GeH+2)和2Π(GeH-2),GeH+2和GeH-2的势能函数呈明显的‘火山态’型,导出了相应的分子离子的解析势能函数、光谱数据和力常数,比较了四参数、八参数Murrell-Sorbie势和Zhu-Wang势对不稳定分子势能函数的拟合情况,指出了用八参数Murrell-Sorbie势对‘火山态’型势能函数的拟合也是合适的;同时讨论了电荷对势能函数和能级的影响。     

Mixed second and third energy derivatives from auxiliary density perturbation theory†

Abstract

The working equations for the calculation of mixed second- and third-order energy derivatives in the framework of auxiliary density functional theory are presented. The perturbations with respect to nuclear displacements and external homogeneous electric field components are calculated with auxiliary density perturbation theory. The presented energy derivative working equations were implemented in deMon2k and validated by vibrational spectra simulations within the double harmonic approximation. The effect of the auxiliary functions on the IR and Raman spectra simulation were analysed for the C₆₀ fullerene. As applications, vibrational spectra of icosahedral carbon fullerenes with up to 540 atoms are calculated without employing symmetry constraints.