Comparative theoretical studies of differently bridged nitramino-substituted ditetrazole 2-<Emphasis Type="Italic">N</Emphasis>-oxides with high detonation performance and an oxygen balance of around zero

In this work, six (A–F) nitramino (–NHNO₂)-substituted ditetrazole 2-N-oxides with different bridging groups (–CH₂–, –CH₂–CH₂–, –NH–, –N=N–, and –NH–NH–) were designed. The six compounds were based on the parent compound tetrazole 2-N-oxide, which possesses a high oxygen balance and high density. The structure, heat of formation, density, detonation properties (detonation velocity D and detonation pressure P), and the sensitivity of each compound was investigated systematically via density functional theory, by studying the electrostatic potential, and using molecular mechanics. The results showed that compounds A–F all have outstanding energetic properties (D: 9.1–10.0 km/s; P: 38.0–46.7 GPa) and acceptable sensitivities (h ₅₀: 28–37 cm). The bridging group present was found to greatly affect the detonation performance of each ditetrazole 2-N-oxide, and the compound with the –NH–NH– bridging group yielded the best results. Indeed, this compound (F) was calculated to have comparable sensitivity to the famous and widely used high explosive 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX), but with values of D and P that were about 8.7% and 19.4% higher than those for HMX, respectively. The present study shows that tetrazole 2-N-oxide is a useful parent compound which could potentially be used in the design of new and improved high-energy compounds to replace existing energetic compounds such as HMX.     

Theoretical study of the effect of N-oxides on the performances of energetic compounds

In order to study the effects of N-oxide on structure and performance, six categories of energetic compounds were systemically investigated. The results indicated that the C–C bonds in the rings were shortened, and the C–N bonds close to the N?→?O bond were elongated when N atoms was oxidized to form N?→?O bonds. N?→?O bonds can increase the densities of most categories of compounds, and the increment will increase with the number of N?→?O bonds. As to their detonation performances, almost all categories of compounds had an increased trend, except for some NO₂-, NHNO₂- and ONO₂-substituted compounds. The contribution of 1,2,3,4-tetrazine and 1,2,4,5-tetrazine to performances was better than that of pyrazine and [1,2,5] oxadiazolo [3,4-b] pyrazine on the whole, and the groups, especially energetic groups, made a huge contribution to performance. When R was a NH₂ or ONO₂ group, all compounds had lower impact sensitivities, and thus represent candidates for novel energetic compounds. However, other than the sixth category of compounds, all compounds had higher impact sensitivities when R was a NO₂ or NHNO₂ group, and have little significance in application.

Open image in new window

Graphical abstract To study the effects of N-oxide on the structure and performance of energetic compounds, and to propose theoretical guidance for the design of novel compounds, the six categories (94 species) listed in the figure were investigated systemically by density functional theory methods and some empirical formulae

     

Computational study of energetic nitrogen-rich derivatives of 1,4-bis(1-azo-2,4-dinitrobenzene)-iminotetrazole

The heats of formation (HOFs), electronic structure, energetic properties, and thermal stabilities for a series of 1,4-bis(1-azo-2,4-dinitrobenzene)-iminotetrazole derivatives with different substituents and substitution positions and numbers of nitrogen atoms in the nitrobenzene rings were studied using the DFT-B3LYP method. All the substituted compounds have higher HOFs than their parent compounds. As the number of nitrogen atoms in the nitrobenzene ring increases, the HOFs of the derivatives with the same substituent rise gradually. Replacing carbon atoms in the nitrobenzene with nitrogen atoms to form N–N bonds is very helpful in improving their HOFs. Most of the substituted compounds have higher HOMO–LUMO gaps than the corresponding unsubstituted compounds. Substitution of the –NO₂, –NF₂, or –ONO₂ group and an increase in the number of nitrogen atoms in the nitrobenzene rings are useful for enhancing their detonation performance. The substituents’ substitution is not favorable for improving thermal stability. Considering detonation performance and thermal stability, five compounds may be considered potential candidates for high energy density compounds (HEDCs).     

A theoretical study on 1,5-diazido-3-nitrazapentane (DANP) and 1,7-diazido-2,4,6-trinitrazaheptane (DATNH): molecular and crystal structures, thermodynamic and detonation properties, and pyrolysis mechanism

1,5-Diazido-3-nitrazapentane (DANP) and 1,7-diazido-2,4,6-trinitrazaheptane (DATNH) are two energetic plasticizers. To better understand them, a detailed theoretical investigation was carried out using density functional theory and molecular mechanics methods. The crystal structures, spectra, thermodynamic properties, heats of formation, detonation velocity, detonation pressure, specific impulse and thermal stability were estimated. Possible initiation steps of pyrolysis were discussed by considering the bond breaking of N–NO₂, C–N₃, and N–N₂ (via hydrogen transfer) for both compounds and the cyclization of the adjacent nitro and azido groups for DATNH. Results show that the rupture of N–NO₂ and N–N₂ (via hydrogen transfer) may happen simultaneously as the initial step of pyrolysis. Both crystals have P-1 symmetry as was observed experimentally. DANP has higher stability than DATNH, while DATNH has better detonation performance than DANP. In addition, DANP has a lower while DATNH has a higher specific impulse than RDX, which shows their prospects as propellant components.     

Theoretical design of novel energetic salts derived from bicyclo-HMX

We designed three novel cage energetic anions by introducing ionic bridges containing NΘ, N(OΘ) and N(NΘNO₂) into cis-2,4,6,8-tetranitro-1H,5H-2,4,6,8- tetraazabicyclo[3.3.0] octane (bicyclo-HMX or BCMHX). The properties of 21 energetic salts, based on cage anions and ammonium-based cations, were studied by density functional theory (DFT) and volume-based thermodynamics (VBT) calculations. Compared to the parent nonionic BCHMX, most title salts have lower predicted impact sensitivities, higher predicted densities, larger predicted heats of formation (HOFs) and better predicted detonation properties. In particular, 11 energetic salts not only exhibit excellent predicted energetic properties, superior to 2,4,6,8,10,12-hexanitrohexaazaisowurtzitane (CL-20), but also have lower predicted sensitivity than CL-20. The best salt had a predicted detonation velocity of 10.06 km s^?1, a predicted detonation pressure of 48.54 GPa and a predicted sensitivity (h₅₀) of 23.99 cm. By introducing ionic bridges into highly nitrated rings, or modifying the original bridge with ionic bridges, some highly nitrated cage compounds with both excellent performance and low sensitivity can be developed strategically.

Open image in new window

Graphical abstract Heats of detonation, detonation velocities, and detonation pressures of salts derived from bicyclo-HMX

     

Probing the effect of nitro groups in nitramine based energetic molecules: a DFT and charge density study

The structure, electron density distribution, energetic and electrostatic properties of simple nitramine based energetic TMA, DMNA, MDA and TNA molecules were determined using density functional theory (B3LYP) with the 6-311G^** and aug-cc-pVDZ basis sets coupled with Bader's theory of atoms in molecules. In the NO₂ group substituted molecules, the N–N bond distance increases with the increase of NO₂ groups, whereas in C–N bonds, this effect is relatively less, and the distances are almost equal. The topological analysis of electron density reveals that the electron density ρ_bcp(r) of C–N and N–N bonds are significantly decreasing with the increase of NO₂ groups in the nitramine molecules. The Laplacian of electron density ▽²ρ_bcp(r) of N–NO₂ bonds [DMNA: ? 16.7 eÅ^? 5, MDA: ? 12.8 eÅ^? 5 and TNA: ? 7.9 eÅ^? 5] of the molecules are relatively less negative, and the values also decrease with the increase of NO₂ groups; this implies that the charge concentration decreases with the increase of NO₂ groups, which leads to weakening the N–N bonds of the molecules. The isosurface of molecular electrostatic potential displays high electronegative regions around the NO₂ groups. The oxygen balance OB₁₀₀ of the molecules increases as the number of NO₂ group increases in the molecules, in which, the TNA molecule having maximum OB₁₀₀ value [+7.89]. The band gap, heat of detonation, bond dissociation energy and charge imbalance are predominantly depends on the number of NO₂ group present in the molecule. The charge imbalance parameter (ν) has been calculated for all molecules, which reveals that TNA is a highly sensitive molecule, the corresponding ν value is 0.047.     

Theoretical design of energetic nitrogen-rich derivatives of 1,7-diamino-1,7-dinitrimino-2,4,6-trinitro-2,4,6-triazaheptane

The heats of formation (HOFs), energetic properties, and thermal stability of a series of 1,7-diamino-1,7-dinitrimino-2,4,6-trinitro-2,4,6-triazaheptane derivatives with different substituents, different numbers of substituents, and different original chains are found by using the DFT-B3LYP method. The results show that -NO₂ or -NH₂ is an effective substituent for increasing the gas-phase HOFs of the title compounds, especially -NO₂ group. As the numbers of substitutents increase, their HOFs enhance obviously. Increasing the length of original chain is helpful for improving their HOFs. The substitution of -NO₂ is useful for enhancing their detonation performances and the effects of the length of original chains on detonation properties are coupled with those of the substituents. An analysis of the BDE of the weakest bonds indicates that the substitution of the -NH₂ groups and replacing the -NO₂ groups of N-NO₂ by the -NH₂ groups are favorable for improving their thermal stability, while the substitution of -NO₂ and increasing the length of original chain decrease their thermal stability. Considering the detonation performance and thermal stability, seven compounds may be considered as the potential candidates of high energy density compounds.     

Is 1-nitro-1-triazene a high energy density material?

An azo bridge (–N?=?N–) can not only desensitize explosives but also dramatically increase their heats of formation and explosive properties. Amino and nitro are two important high energy density functional groups. Here, we present calculations on 1-nitro-1-triazene (NH₂–N?=?N–NO₂). Thermal stability and detonation parameters were predicted theoretically at CCSD(T)/6-311G* level, based on the geometries optimized at MP2/6-311G* level. It was found that the p?→?π conjugation interaction and the intramolecular hydrogen bonding that exist in the system together increase the thermal stability of the molecule. Moreover, the detonation parameters were evaluated to be better than those of the famous HMX and RDX. Finally, the compound was demonstrated to be a high energy density material.     

Theoretical studies on the crystal structure, thermodynamic properties, detonation performance and thermal stability of cage-tetranitrotetraazabicyclooctane as a novel high energy density compound

The B3LYP/6-31G (d) method of density functional theory (DFT) was used to study molecular geometry, electronic structure, infrared spectrum (IR) and thermodynamic properties. The heat of formation (HOF) and calculated density were estimated to evaluate the detonation properties using Kamlet–Jacobs equations. Thermal stability of 3,5,7,10,12,14,15,16-octanitro- 3,5,7,10,12,14,15,16-octaaza-heptacyclo[7.5.1.1^2,8.0^1,11.0^2,6.0^4,13.0^6,11]hexadecane (cage-tetranitrotetraazabicyclooctane) was investigated by calculating the bond dissociation energy (BDE) at unrestricted B3LYP/6-31G (d) level. The calculated results show that the N–NO₂ bond is a trigger bond during thermolysis initiation process. The crystal structure obtained by molecular mechanics (MM) methods belongs to Pna2₁ space group, with cell parameters a?=?12.840 Å, b?=?9.129 Å, c?=?14.346 Å, Z?=?6 and ρ?=?2.292 g·cm^?3. Both the detonation velocity of 9.96 km·s^?1 and the detonation pressure of 47.47 GPa are better than those of CL-20. According to the quantitative standard of energetics and stability, as a high energy density compound (HEDC), cage-tetranitrotetraazabicyclooctane essentially satisfies this requirement.     

Theoretical insights into the structures and mechanical properties of HMX/NQ cocrystal explosives and their complexes,and the influence of molecular ratios on their bonding energies

Molecular dynamics (MD) methods were employed to study the binding energies and mechanical properties of selected crystal planes of 1,3,5,7-tetranitro-1,3,5,7-tetrazacyclooctane (HMX)/nitroguanidine (NQ) cocrystals at different molecular molar ratios. The densities and detonation velocities of the cocrystals at different molar ratios were estimated. The intermolecular interaction and bond dissociation energy (BDE) of the N–NO₂ bond in the HMX:NQ (1:1) complex were calculated using the B3LYP, MP2(full) and M06-2X methods with the 6-311++G(d,p) and 6-311++G(2df,2p) basis sets. The results indicated that the HMX/NQ cocrystal prefers cocrystalizing in a 1:1 molar ratio, and the cocrystallization is dominated by the (0 2 0) and (1 0 0) facets. The K, G, and E values of the ratio of 1:1 are smaller than those of the other ratios, and the 1:1 cocrystal has the best ductility. The N–NO₂ bond becomes stronger upon the formation of the intermolecular H-bonding interaction and the sensitivity of HMX decreases in the cocrystal. This sensitivity change in the HMX/NQ cocrystal originates not only from the formation of the intermolecular interaction but also from the increment of the BDE of N–NO₂ bond in comparison with isolated HMX. The HMX/NQ (1:1) cocrystal exhibits good detonation performance. Reduced density gradient (RDG) reveals the nature of cocrystallization. Analysis of the surface electrostatic potential further confirmed that the sensitivity decreases in complex (or cocrystal) in comparison with that in isolated HMX.

Open image in new window

Graphical Abstract Binding energies and mechanical properties of HMX/NQ cocrystals in different molecular molar ratios were studied using molecular dynamics methods. The origin of the sensitivity change in the HMX/NQ cocrystal originates from formation of intermolecular interactions and the bond dissociation energy increment of the N–NO₂ bond

     

Does the environment around the H-cluster allow coordination of the pendant amine to the catalytic iron center in [FeFe] hydrogenases? Answers from theory

[FeFe] hydrogenases are H₂-evolving enzymes that feature a diiron cluster in their active site (the [2Fe]_H cluster). One of the iron atoms has a vacant coordination site that directly interacts with H₂, thus favoring its splitting in cooperation with the secondary amine group of a neighboring, flexible azadithiolate ligand. The vacant site is also the primary target of the inhibitor O₂. The [2Fe]_H cluster can span various redox states. The active-ready form (H_ox) attains the Fe^IIFe^I state. States more oxidized than H_ox were shown to be inactive and/or resistant to O₂. In this work, we used density functional theory to evaluate whether azadithiolate-to-iron coordination is involved in oxidative inhibition and protection against O₂, a hypothesis supported by recent results on biomimetic compounds. Our study shows that Fe–N(azadithiolate) bond formation is favored for an Fe^IIFe^II active-site model which disregards explicit treatment of the surrounding protein matrix, in line with the case of the corresponding Fe^IIFe^II synthetic system. However, the study of density functional theory models with explicit inclusion of the amino acid environment around the [2Fe]_H cluster indicates that the protein matrix prevents the formation of such a bond. Our results suggest that mechanisms other than the binding of the azadithiolate nitrogen protect the active site from oxygen in the so-called H _ox ^inact state.     

A comparative theoretical investigation into the strength of the trigger-bond in the Na+, Mg2+ and HF complexes involving the nitro group of R–NO2 (R = –CH3, –NH2 and –OCH3) or the C = C bond of (E)-O2N–CH = CH–NO2

A comparative theoretical investigation into the change in strength of the trigger-bond upon formation of the Na⁺, Mg²⁺ and HF complexes involving the nitro group of RNO₂ (R?=? –CH₃, –NH₂, –OCH₃) or the C?=?C bond of (E)-O₂N–CH?=?CH–NO₂ was carried out using the B3LYP and MP2(full) methods with the 6-311++G**, 6-311++G(2df,2p) and aug-cc-pVTZ basis sets. Except for the Mg²⁺?π system with (E)-O₂N–CH?=?CH–NO₂ (i.e., C₂H₂N₂O₄?Mg²⁺), the strength of the trigger-bond X–NO₂ (X?=?C, N or O) was enhanced upon complex formation. Furthermore, the increment of bond dissociation energy of the X–NO₂ bond in the Na⁺ complex was far greater than that in the corresponding HF system. Thus, the explosive sensitivity in the former might be lower than that in the latter. For C₂H₂N₂O₄?Mg²⁺, the explosive sensitivity might also be reduced. Therefore, it is possible that introducing cations into the structure of explosives might be more efficacious at reducing explosive sensitivity than the formation of an intermolecular hydrogen-bonded complex. AIM, NBO and electron density shifts analyses showed that the electron density shifted toward the X–NO₂ bond upon complex formation, leading to a strengthened X–NO₂ bond and possibly reduced explosive sensitivity.

Characterization of the hydrogen bond network in guanosine quartets by internucleotide 3hJNC′ and 2hJNN scalar couplings

Scalar coupling correlations across hydrogen bonds with carbonyl groups as acceptors have been observed in a variety of proteins, but not in nucleic acids. Here we present a pulse scheme that allows such an observation and quantification of trans-hydrogen bond ^3hJ_NC correlations in nucleic acid base pairs, between the imino nitrogen ¹⁵N1 and the carbonyl ¹³C6 nuclei within the guanine quartets of the Oxy-1.5 DNA-quadruplex. Intra- and internucleotide N-H···O=C connectivities can be traced around each guanine quartet, allowing the hydrogen bonding partners to be unambiguously assigned. Absolute values of the ^3hJ_NC couplings are approximately 0.2 Hz as quantified by a selective long-range H(N)CO experiment and are thus on average smaller than the analogous ^3hJ_NC couplings observed in proteins. In addition, an improved version of the pseudo-heteronuclear H(N)N-COSY [Majumdar et al. (1999) J. Biomol. NMR, 14, 67–70] is presented which allows simultaneous detection of the ¹⁵N-donor and ¹⁵N-acceptor resonances connected by ^2hJ_NN couplings in hydrogen bonds involving amino groups. Using this experiment, values ranging between 6 and 8 Hz are determined for the ^2hJ_NN couplings between ¹⁵N2 and ¹⁵N7 nuclei in the guanine quartet. These values are not strongly influenced by the presence of a significant amount of chemical exchange broadening due to amino group rotations.     

Computational studies on polynitropurines as potential high energy density materials

As part of a search for high energy density materials (HEDMs), a series of purine derivatives with nitro groups were designed computationally. The relationship between the structures and the performances of these polynitropurines was studied. Density functional theory (DFT) at the B3LYP/6-311G** level was employed to evaluate the heats of formation (HOFs) of the polynitropurines by designing an isodesmic reaction method. Results indicated that the HOFs were influenced by the number and positions of substituent groups. Detonation properties were evaluated using the Kamlet–Jacobs equations, based on the theoretical densities and heats of formation of the polynitropurines. The relative stabilities of the polynitropurines were studied via the pyrolysis mechanism and the UB3LYP/6-311G** method. Homolysis of the ring–NO₂ bond is predicted to be the initial step in the thermal decomposition of these purine derivatives. Considering their detonation properties and relative stabilities, the tetranitropurine (D1) derivatives may be regarded as potential candidates for practical HEDCs. These results may provide useful information for further investigations.     

Structural, molecular mechanics, and DFT study of cadmium(II) in its crown ether complexes with axially coordinated ligands, and of the binding of thiocyanate to cadmium(II)

The role of relativistic effects (RE) in the structures of Cd(II) complexes with crown ethers, and the reason the ‘soft’ Cd(II) strongly prefers to bind to SCN⁻ through N, are considered. The synthesis and structures of [Cd(18-crown-6)(thiourea)₂] (ClO₄)₂.18-crown-6 (1) and [Cd(Cy₂-18-crown-6)(NCS)₂] (2) are reported. (18-crown-6 = 1,4,7,10,13,16-hexaoxacyclooctadecane; Cy₂-18-crown-6 = cis-anti-cis-2,5,8,15,18,21-hexaoxatricylo[20.4.0.0(9,14)]hexacosane). In 1 Cd is coordinated in the plane of the crown which has close to D_3d symmetry, with long Cd-O bonds averaging 2.688 Å. The two thiourea molecules form relatively short Cd-S bonds that average 2.468 Å, with an S-Cd-S angle of 164.30°. This structure conforms with the idea that Cd(II) can adopt a near-linear structure involving two covalently-bound donor atoms (the S-donors) with short Cd-S bonds, which resembles gas-phase structures for species such as CdCl₂. The structure of 2 is similar, with the two SCN⁻ ligands N-bonded to Cd, with short Cd-N bonds of 2.106 Å, and N-Cd-N angle of 180°. The crown in 2 forms long Cd-O bonds that average 2.698 Å. Molecular mechanics calculations suggest that a main reason Cd(II) prefers to bind to SCN⁻ through N is that when bound through S, the small Cd-S-C angle, which is typically close to 100°, brings the ligand into close contact with other ligands present, and causes steric destabilization. In contrast, the Cd-N-C angles for SCN⁻ coordinated through N are much larger, being 171.4° in 2, which keeps the SCN⁻ groups well clear of the crown ether. DFT (density functional theory) calculations are used to generate the structures of [Cd(18-crown-6)(H₂O)₂]²⁺ (3) and [Cd(18-crown-6)Cl₂] (4). In 3, the Cd(II) is bound to only three O-donors of the macrocycle, with Cd-O bonds averaging 2.465 Å. The coordinated waters form an O-Cd-O angle of 139.47°, with Cd-O bonds of 2.295 Å. In contrast, for 4, the Cd is placed centrally in the cavity of the D_3d symmetry crown, with long Cd-O bonds averaging 2.906 Å. The Cl groups form a Cl-Cd-Cl angle of 180°, with short Cd-Cl bonds of 2.412 Å. With ionically bound groups on the axial sites of[Cd(18-crown-6)X₂] complexes, such as with X = H₂O in 3, the Cd(II) does not adopt linear geometry involving the two X groups, with long Cd-O bonds to the O-donors of the macrocycle. With covalently-bound X = Cl in 4, short Cd-Cl bonds and a linear [Cl-Cd-Cl] unit results, with long Cd-O bonds to the crown ether.     

Theoretical studies on the structure, thermochemical and detonation properties of amino and nitroso substituted 1,2,4-triazol-5-one-N-oxides

DFT calculations at the B3LYP/aug-cc-pVDZ level have been carried out to explore the structure, stability, electron density, heat of formation, detonation velocity and detonation pressure of substituted amino- and nitroso-1,2,4-triazol-5-one-N-oxides. Heats of formation of substituted triazol-5-one-N-oxides have been computed at the B3LYP/aug-cc-pVDZ level via isodesmic reaction procedure. Materials Studio 4.1 package was used to predict the crystal density of model compounds. Kamlet-Jacob equations were used to calculate detonation properties based on the calculated heat of explosion and crystal density. The designed compounds 4, 6, 7 and 8 have shown higher performance compared with those of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane and octanitrocubane. Atoms-in-molecule (AIM) analyses have also been carried out to understand the nature of intramolecular interactions in the designed molecules.     

Theoretical studies on the thermodynamic properties, densities, detonation properties, and pyrolysis mechanisms of trinitromethyl-substituted aminotetrazole compounds

Trinitromethyl-substituted aminotetrazoles with –NH₂, –NO₂, –N₃, and –NHC(NO₂)₃ groups were investigated at the B3LYP/6-31G(d) level of density functional theory. Their sublimation enthalpies, thermodynamic properties, and heats of formation were calculated. The thermodynamic properties of these compounds increase with temperature as well as with the number of nitro groups attached to the tetrazole ring. In addition, the detonation velocities and detonation pressures of these compounds were successfully predicted using the Kamlet–Jacobs equations. It was found that these compounds exhibit good detonation properties, and that compound G (D = 9.2 km/s, P = 38.8 GPa) has the most powerful detonation properties, which are similar to those of the well-known explosive HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazocine). Finally, the electronic structures and bond dissociation energies of these compounds were calculated. The BDEs of their C–NO₂ bonds were found to range from 101.9 to 125.8 kJ/mol^-1. All of these results should provide useful fundamental information for the design of novel HEDMs.     

Molecular design of new nitramine explosives: 1,3,5,7-tetranitro-8-(nitromethyl)-4-imidazolino[4,5-b]4-imidazolino-[4,5-e] pyridine and its N-oxide

Two new nitramine compounds containing pyridine, 1,3,5,7-tetranitro-8-(nitromethyl) -4-imidazolino[4,5-b]4-imidazolino-[4,5-e]pyridine and its N-oxide 1,3,5,7-tetranitro-8- (nitromethyl)-4-imidazolino[4,5-b]4-imidazolino-[4,5-e]pyridine-4-ol were proposed. Density functional theory (DFT) has been employed to study the molecular geometries, electronic structures, infrared spectra, and thermodynamic properties at the B3LYP/6-31G* level. Their detonation performances evaluated using the Kamlet-Jacobs equations with the calculated densities and heats of formation are superior to those of HMX. The predicted densities of them were ca. 2 g*cm^-3, detonation velocities were over 9 km*s^-1, and detonation pressures were about 40 GPa, showing that they may be potential candidates of high energy density materials (HEDMs). The natural bond orbital analysis indicated that N-NO₂ bond is the trigger bond during thermolysis process. The stability of the title compounds is slightly lower than that of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12- hexaazaisowurtzitane (CL-20). The results of this study may provide basic information for the molecular design of new HEDMs.     

Screening for potential energetic C–N cages with high energy and good stability: a theoretical comparative study

This paper presents a theoretical study on carbon-nitrogen cages as potentially high energy density materials (HEDMs) using density functional theory. The energetic properties, detonation performance, and stability of two C₆N₆H₁₂ cages were researched comparison with two similar common cage compounds hexaazaisowurtzitane and cubane. Results indicate that both of two C₆N₆H₁₂ cages have high positive heat of formation and good stability. Their densities and detonation characteristics are equivalent or slightly superior to those of hexaazaisowurtzitane and cubane, indicating that they have good detonation performance.     

Modulating weak intramolecular interactions through the formation of beryllium bonds: complexes between squaric acid and BeH2

The electronic structure of the two most stable isomers of squaric acid and their complexes with BeH₂ were investigated at the B3LYP/6-311?+?G(3df,2p)// B3LYP/6-31?+?G(d,p) level of theory. Squaric acid forms rather strong beryllium bonds with BeH₂, with binding energies of the order of 60 kJ?mol^?1. The preferential sites for BeH₂ attachment are the carbonyl oxygen atoms, but the global minima of the potential energy surfaces of both EZ and ZZ isomers are extra-stabilized through the formation of a BeH···HO dihydrogen bond. More importantly, analysis of the electron density of these complexes shows the existence of significant cooperative effects between the beryllium bond and the dihydrogen bond, with both becoming significantly reinforced. The charge transfer involved in the formation of the beryllium bond induces a significant electron density redistribution within the squaric acid subunit, affecting not only the carbonyl group interacting with the BeH₂ moiety but significantly increasing the electron delocalization within the four membered ring. Accordingly the intrinsic properties of squaric acid become perturbed, as reflected in its ability to self-associate.