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.     

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.     

Theoretical studies on nitrogen‐rich furoxan‐based heterocyclic derivatives

The density functional theory method was used to study the heats of formation (HOFs), energetic properties, electronic structure of a series of 4,4″‐dinitro(3,3′:4′,3′′)tris([1,2,5]oxadiazole)‐2′‐oxide (3,4‐bis[4′‐nitrofurazan‐3′‐yl]furoxan) derivatives. The results show that the substitution of the nitro group is very useful for improving their HOFs and detonation performances. The HOFs of the title compounds are all positive and larger than those of 1,3,5‐trinitro‐1,3,5‐triazinane and 1,3,5,7‐tetranitro‐1,3,5,7‐tetrazocane. The analysis of oxygen balance shows that the studied compounds need the oxygen in the explosive. Compound A1 has larger detonation velocity and detonation pressure than those of 1,3,5,7‐tetranitro‐1,3,5,7‐tetrazocane and can be regarded as a potential candidate for high‐energy compounds because of the moderate heat of detonation, high density, and high N. In addition, the energy gaps between the highest occupied molecular orbital and lowest unoccupied molecular orbital of the studied compounds are further investigated.     

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.     

Theoretical investigation of two‐photon absorption and fluorescence properties of cypridina luciferin‐based derivatives: 2,3,5‐trisubstituted pyrazine compounds

Two‐photon fluorescent probe materials are significant for achieving observation of living phenomena in entire organs and tissues. To explore new materials with high fluorescence and large two‐photon absorption (TPA) cross section, a series of 2,3,5‐trisubstituted pyrazine derivatives were designed. Their equilibrium geometries, one‐photon absorption, TPA, and luminescence properties have been studied by using density functional theory (DFT), time‐dependent DFT, and Zerner's intermediate neglect of differential overlap program. The results show that the introduction of styrene groups to 2,3,5‐trisubstituted pyrazine derivatives can efficiently increase the conjugated effect and enhance the TPA activity. Moreover, the luminescence properties of 2,3,5‐trisubstituted pyrazine derivatives were compared, and the effect of three substituents on the fluorescence of trisubstituted pyrazine derivatives was analyzed by means of different contribution of the basis functions localized on pyrazine fragment into the highest occupied molecular orbital and lowest unoccupied molecular orbital. The oscillator strengths in the excited state (?_em) for the pyrazine derivatives substituted by styryl are larger than that of other derivatives with acetylamino and indole groups substituted at 2‐site and 5‐site of the pyrazine core, and the ?_em of 3‐indolyl pyrazine derivatives is larger than that of 3‐styrene pyrazine derivatives. It suggests that the styrene group has a great influence on the luminescence property. In addition, the indole group substituted at 3‐site of the pyrazine derivatives can also promote the fluorescence property. Copyright © 2013 John Wiley & Sons, Ltd.     

Structure and detonation performance of a novel HMX/LLM‐105 cocrystal explosive

Intermolecular interactions and properties of octahydro‐1,3,5,7‐tetranitro‐1,3,5,7‐ tetrazocine (HMX) / 2,6‐diamino‐3,5‐dinitropyrazine‐1‐oxide (LLM‐105) cocrystal were studied by using the dispersion‐corrected density functionals (ωB97XD, B97D) and meta‐hybrid functional (M062x) methods. Binding energies, heats of formation, thermodynamic properties, atoms in molecules, and natural bond orbital analysis were performed to investigate HMX/LLM‐105 complexes. Results show that the main intermolecular interactions between HMX and LLM‐105 are CH…O, NH…O, N…O, and O…O interactions. In addition, Monte Carlo simulation was employed to predict the crystal structure of HMX/LLM‐105 cocrystal. The HMX/LLM‐105 cocrystal is most likely to crystallize in C2/c space group, and its corresponding cell parameters are Z = 8, a = 41.63 Å, b = 6.77 Å, c = 45.63 Å, ß = 164.55°, and ρ = 1.99 g/cm³. Detonation velocity and pressure of HMX/LLM‐105 cocrystal are 8.95 km/s, 37.69GPa, a little lower than those of HMX (9.10 km/s, 37.76GPa). However, according to the net charges of nitro group, HMX/LLM‐105 cocrystal exhibits less sensitive than HMX. Finally, bond dissociation energy calculation shows that HMX/LLM‐105 complexes are thermally stable. Considering thermal stability, sensitivity, and detonation performance, HMX/LLM‐105 cocrystal meets the requirements of insensitive high energy density materials. Copyright © 2013 John Wiley & Sons, Ltd.     

几种咪唑-二氧杂环的含能取代衍生物的理论研究

使用DFT的B3LYP方法对几种咪唑二氧杂环化合物的分子结构、红外光谱、生成焓、爆轰性能和化学/热稳定性进行了研究.四种不同含能基团-NO₂,-NH₂,-N₃和-ONO₂对该化合物各项性能的影响进行了比较.结果表明-NO₂和-ONO₂基团有效地增加了化合物的密度,而-N³基团极大地增加了化合物的生成焓.其中-NO₂取代物爆轰性能接近1,3,5,7-tetranitro-1,3,5,7-tetraazacyclooctane,-ONO₂取代物的爆轰性能低于hexahydro-1,3,5-trinitro-1,3,5-triazine.自然键轨道分析表明,这些化合物中相对较弱的键为取代基和咪唑环之间的键,以及二氧杂环中的C-O键.吸电子基团(-NO₂, -N³和-ONO₂)在分子中产生了明显的诱导效应,减弱了基团与咪唑环之间的链接,降低了对应键的键裂解能.研究表明基团的电负性与化合物的稳定性有关.考虑到化合物的爆轰性能和热稳定性,DNTNDI化合物满足高能量密度材料的要求.     

DFT study on the effect of relative positions of methyl-, nitro- and N→oxide groups on the molecular structure,thermal/kinetic stability,crystal density,heat of decomposition and performance characteristics of triazolones

Methyl-, nitro- and N→oxide substituted triazolones are of interest in the contest of high-energy density compounds and have been found to have true local energy minima at the B3LYP/aug-cc-pVDZ level. The optimised structures, harmonic frequencies and thermodynamic values for all the model molecules have been obtained in their ground state. The velocity of detonation (D) and detonation pressure (P) have been evaluated by the Kamlet–Jacob equations using the crystal density and the heat of explosion. The estimated performance properties are higher (D = 9.92–10.27 km/s, P = 48.10–52.52 GPa) compared with 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (D = 9.20 km/s, P = 42.0 Gpa). The higher densities are possibly due to the intramolecular hydrogen bonds and the layered structures in the crystal lattice. We speculate that the calculated heat of explosion and the density are for the gas phase compounds and in the reality they should be for the solid phase which would diminish the magnitude of the calculated values. The –N→O and –NO₂ group leads to the desirable consequences of higher heat of explosion and diminished sensitivities. The substituting of N–H hydrogen atom(s) of triazolones for a –CH₃ group decreases melting point, heat of formation and density; however, the methyl group increases the thermal stability.     

Theoretical studies on 2‐(5‐amino‐3‐nitro‐1,2,4‐triazolyl)‐3,5‐dinitropyridine (PRAN) and its derivatives

In this work, a set of derivatives of 2‐(5‐amino‐3‐nitro‐1,2,4‐triazolyl)‐3,5‐dinitropyridine (PRAN) with different energetic substituents (?N₃, –NO₂, –NH₂, –NF₂) have been studied at the Becke, three‐parameter, Lee–Yang–Parr/aug‐cc‐pvdz, Becke, three‐parameter, Lee–Yang–Parr/6‐31G(d), Becke, three‐parameter, Perdew 86/6‐31G(d), and Becke three‐parameter, Perdew–Wang 91/6‐31G(d,p) levels of density functional theory. The gas‐phase heats of formation were predicted with isodesmic reactions and the condensed‐phase HOFs were estimated with the Politzer approach. The effects of different functionals and basis sets were analyzed. –N₃ and –NO₂ greatly increase while –NH₂ and –NF₂ slightly decrease heats of formation. An analysis of the bond dissociation energies and impact sensitivity shows that all compounds have good stability. The crystal densities (1.82–2.00 g/cm³) computed from molecular packing calculations are big for all compounds and that of the –NF₂ derivative is the largest. All derivatives have higher detonation velocity and detonation pressure than PRAN. Compounds 3 and 4 (R = NO₂ and NF₂) have better performance than hexahydro‐1,3,5‐trinitro‐1,3,5‐trizine and the performance of 4 is quite close to that of 1,3,5,7‐tetranitro‐1,3,5,7‐tetraazacyclooctane, they are promising candidates of high energy compounds and worth further investigations. Copyright © 2012 John Wiley & Sons, Ltd.     

Theoretical investigations of pyridine derivatives as potential high energy density materials

Density function theory has been employed to study pyridine derivatives at the B3LYP/6‐31 G(d,p) and B3P86/6‐31 G(d,p) levels. The crystal structures were obtained by molecular mechanics methods. The heats of formation (HOFs) were predicted based on the isodesmic reactions. Detonation performance was evaluated by using the Kamlet–Jacobs equations based on the calculated densities and HOFs. The thermal stability of the title compounds was investigated by the bond dissociation energies and the energy gaps (ΔE_LUMO?HOMO) predicted. It is found that there are good linear relationships between detonation velocity, detonation pressure, and the number of nitro group. The simulation results reveal that molecule G performs similar to the famous explosive HMX and molecule D outperforms HMX. According to the quantitative standard of energetics and stability as high energy density materials, molecule D essentially satisfies this requirement. These results provide basic information for molecular design of novel high energetic density materials. Copyright © 2012 John Wiley & Sons, Ltd.     

Theoretical insights on a series of difluoramino group–based energetic molecules

A series of difluoramino group–based energetic molecules was designed and the relative properties were investigated by density functional theory. The results show that all the designed molecules have high positive heat of formation which ranges from 479.48 to 724.02 kJ/mol, detonation velocity ranges from 8.01 to 11.26 km/s, detonation pressure ranges from 28.03 to 63.46 GPa, and impact sensitivity ranges from 18.2 to 54.5 cm. Then, compounds D2, D3, D5, E4, E5, E6, and F2 were selected as the potential high energy density materials based on detonation properties and sensitivities. Natural bond orbital charges, electronic density, frontier molecular orbital, electrostatic potential on the surface, and thermal dynamic parameters of the screened molecules (compounds D2, D3, D5, E4, E5, E6, and F2) were also predicted at B3LYP/6‐31G(d,p) level to give a better understanding on the chemical and physical properties of them.     

Computational studies on two novel energetic nitrogen‐rich compounds based on tetrazolone

Two novel energetic nitrogen‐rich compounds 1,4‐diaminotetrazol‐5‐one ( DATO ) and 1,4‐dinitrotetrazol‐5‐one ( DNTO ) were proposed first and studied by quantum chemistry method with B3LYP/6‐31G* level of theory. The optimized geometry, IR predicted spectrum and thermochemical parameters, frontier molecular orbitals and molecular electrostatic potential were calculated for inspecting the electronic structure, molecular stability and chemical reactivity. The important macroscopic properties including density, enthalpy of formation, detonation parameters and impact sensitivity have been predicted as well. As a result, two designed compounds DATO and DNTO possess positive enthalpy of formation (395.79 and 342.77 kJ/mol), impressive detonation parameters (D = 8.80 km/s, P = 33.69 GPa; D = 8.89 km/s, P = 34.98 GPa) superior to the remarkable explosive RDX, acceptable sensitivities and might be promising candidates of energetic materials. Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

基于4，8-二氢二呋咱[3，}4-b，e]吡嗪高能材料的设计与筛选

为了寻找高能量密度的材料,本文设计了一系列基于4,8-二氢二呋咱[3,4-b,e]吡嗪的含能材料. 利用密度泛函理论研究了它们结构与性质之间的关系. 结果表明,这些设计化合物的性质受到含能基团和杂环取代基的影响. -N₃含能基团是提高设计化合物生成热的最有效取代官能团,而四唑环/-C(NO₂)₃基团对炸药的爆轰性能有较大贡献. 键解离能分析表明,引入-NHNH₂,-NHNO₂,-CH(NO₂)₃和-C(NO₂)₃基团会显著降低键解离能. 由于化合物A8,B8,C8,D8,E8和F8具有良好的爆轰性能和热稳定性,最终被筛选为潜在的高能密度材料. 此外,还计算了这些筛选化合物的电子结构.     

Computer simulations and analysis of structural and energetic features of crystalline cage energetic compound: 2, 4, 6, 8, 12‐pentanitro‐10‐(3, 5, 6‐trinitro (2‐pyridyl))‐2, 4, 6, 8, 12‐hexaazatetracyclo [5.5.0.03,11.05,9]dodecane

First principles molecular orbital and plane‐wave ab initio calculations have been used to investigate the structural and energetic properties of a new cage compound 2, 4, 6, 8, 12‐pentanitro‐10‐(3, 5, 6‐trinitro (2‐pyridyl))‐2, 4, 6, 8, 12‐hexaazatetracyclo [5.5.0.0^3,11.0^5,9]dodecane (PNTNPHATCD) in both the gas and solid phases. The molecular orbital calculations using the density functional theory methods at the B3LYP/6‐31G(d,p) level indicate that both the heat of formation and strain energy of PNTNPHATCD are larger than those of 2, 4, 6, 8, 10, 12‐hexanitro‐2, 4, 6, 8, 10, 12‐hexaazatetracyclo [5.5.0.0.0] dodecane (CL‐20). The infrared spectra and the thermodynamic property in gas phase were predicted and discussed. The calculated detonation characteristics of PNTNPHATCD estimated using the Kamlet–Jacobs equation equally matched with those of CL‐20. Bond‐breaking results on the basis of natural bond orbital analysis imply that C–C bond in cage skeleton, C–N bond in pyridine, and N–NO₂ bond in the side chain of cage may be the trigger bonds in the pyrolysis. The structural properties of PNTNPHATCD crystal have been studied by a plane‐wave density functional theory method in the framework of the generalized gradient approximation. The crystal packing predicted using the Condensed‐phase Optimized Molecular Potentials for Atomistic Simulation Studies (COMPASS) force fields belongs to the Pbca space group, with the lattice parameters a = 20.87 Å, b = 24.95 Å, c = 7.48 Å, and Z = 8, respectively. The results of the band gap and density of state suggest that the N–NO₂ bond in PNTNPHATCD may be the initial breaking bond in the pyrolysis step. As the temperature increases, the heat capacity, enthalpy, and entropy of PNTNPHATCD crystal all increase, whereas the free energy decreases. Considering that the cage compound has the better detonation performances and stability, it may be a superior high energy density compound. Copyright © 2015 John Wiley & Sons, Ltd.     

多硝基六氮杂金刚烷结构和相对稳定性的理论研究

基于密度泛函理论,在B3LYP/6-31G*水平下分别计算了潜在高能量密度化合物(HEDC)多硝基六氮杂金刚烷(PNHAAs)的结构,包括总能量、几何结构、氧平衡、偶极矩等．通过研究PNHAAs的结构参数,如N-NO2最大键长、N-NO2键的最小Mulliken集居数、-NO2所带最小负电荷、氧平衡,预测了10个标题物的相对稳定性或感度(起爆的难易程度．感度越高,化合物的稳定性越低)．结果发现,这四个参数基本上能得到一致的结论:随着硝基数目的增加,标题物的稳定性降低．氧平衡不能判别同分异构体的相对稳定性,但总能量和偶极矩可较好地判别同分异构体的相对稳定性．     

Synthesis,characterization, and electrochemical studies of new 5‐ and 6‐nitro N‐acyl‐1H‐indazoles

Two series of new N‐1 acylindazoles containing 5‐ or 6‐nitro groups were synthesized with moderate to good yields and characterized by IR and NMR spectroscopy. Cyclic voltammetry in aprotic media was utilized for the electrochemical characterization of the compounds. The calculated reduction potentials in physiological conditions are similar to those of known commercial antichagasic drugs. Therefore, the novel series reported herein are prospective candidates for antichagasic biological evaluation. Theoretical calculation results indicate that the studied dinitro derivatives undergo a single step reduction process because of the energy proximity of their radical and anionic state. Copyright © 2011 John Wiley & Sons, Ltd.     

Computational design and screening of promising energetic materials: Novel azobis(tetrazoles) with ten catenated nitrogen atoms chain

Density functional theory methods were used to study on 2 N10 compounds, 1,1′‐azobis(tetrazole) and 1,1′‐azobis(5‐methyltetrazole). We systematically investigated 10 novel substituted azobis(tetrazoles) with 10 catenated nitrogen atoms and various energetic groups (–CF₃ 1 , –C(NO₂)₃ 3 , –N₃ 5 , –NH₂ 6 , –NHNH₂ 7 , –NHNO₂ 8 , –NO₂ 9 , –OCH₃ 10 , –OH 11 , –ONO₂ 12 ). The optimized geometry, frontier molecular orbitals, electrostatic potential, Infrared and nuclear magnetic resonance spectrum were calculated for inspecting the molecular structure and stability as well as chemical reactivity. The effects of different substituents on the density, enthalpy of formation, heat of explosion, detonation velocity and pressure, and sensitivity of the azobis(tetrazole) derivatives have been investigated. Compound 9 with nitro was found to have remarkable detonation performances (D = 9.61 km/s, P = 42.14 GPa), which are close to the excellent explosive CL‐20. Results show that compounds 1 , 3 , 4 , 7 , 9 , 11, and 12 have high potential to replace RDX. It is surprising that compounds 1 , 3 , 9, and 12 possess better energetic properties than HMX. These novel substituted azobis(tetrazoles) with unique N10 structure may be promising candidates of HEDMs with outstanding performance and acceptable sensitivities.     

Theoretical study of one‐ and two‐photon absorption properties of expanded donor–acceptor calix[4]arenes

Two types of donor(D)–acceptor(A) calix[4]arenes have been theoretically studied using DFT//B3LYP/6‐31G(d) method and ZINDO/CISD method. The calculations show that the substitution of C? C by the conjugation bridge C?C and N?N plays an important part in altering one‐photon absorption (OPA) and two‐photon absorption (TPA) properties. The maximum OPA wavelengths of all studied compounds are less than 400 nm, which means high transparency. The geometry of the calixarenes strongly influences the TPA properties of the studied compounds. In addition, the nitro derivatives have a wider TPA response range than other non‐nitro derivatives. The tetrasubstituted calix[4]arenes (type B calixarenes) have a larger TPA cross‐section values than the bisubstituted calix[4]arenes (type A calixarenes). Copyright © 2009 John Wiley & Sons, Ltd.     

Molecular structure and vibrational spectra of 3 (or 4 or 6)‐methyl‐5‐nitro‐2‐pyridinethiones: FT‐IR,FT‐Raman and DFT quantum chemical calculations

IR and Raman spectra (RS) of polycrystalline 3‐(or 4 or 6)‐methyl‐5‐nitro‐2‐pyridinethione have been measured and analyzed by means of density functional theory (DFT) quantum chemical calculations. The B3LYP/6‐311G(2d,2p) approach has been applied for both the thiol and thione tautomers due to the possibility of the formation of these two thiole forms. Molecular structures of these compounds have been optimized starting from different molecular geometries of the thiol group and thione group. Two conformations of the 2‐mercaptopyridine, trans and cis, have been taken into account. It was shown that the studied compounds appear in the solid state in the thione form. The effect of the hydrogen‐bond formation in the studied compounds has been considered. Copyright © 2008 John Wiley & Sons, Ltd.