Modeling of explosives: 1,4,2,3,5,6-dioxatetrazinane as a new green energetic material with enhanced performance
Graphical abstract
Introduction
According to an analysis of the explosives market by Global Industry Analytics, Inc., the global market for explosives is projected to reach 30 million metric tons in 2024 [1]. Since most explosive formulations contain a significant amount of carbon, a huge mass of carbon oxides is released into the atmosphere, exacerbating the greenhouse effect. In this context, carbon-free nitrogen-rich explosives are highly desirable, since these would release environmentally friendly molecular dinitrogen as the primary detonation product. Ideally, pure nitrogen allotropes other than N2 [[2], [3], [4]] could replace conventional C–H–N–O explosives, but, to date, only three allotropes of nitrogen, both molecular and polymeric (cg-N, LP-N, and N8), have been characterized experimentally [[5], [6], [7]]. We should stress that these allotropes have been obtained only in very small quantities under extreme conditions; thus, their practical application is still unachievable.
Meanwhile, various H–N–O compounds have attracted close attention due to their superior performances as high-energy-density materials (HEDM). Among these compounds, the family of ammonium salts with various nitrogen-rich carbon-free anions occupies a special place [[8], [9], [10], [11], [12]]. For example, the burning rate of ammonium dinitramide is about 10 times faster than those of ammonium perchlorate, RDX, and HMX [13]. The most widely used cations and anions, as components of these salts, are illustrated in Chart 1. The latter not only have high nitrogen contents, but also display high enthalpies of formation, favorable detonation properties, and remarkable insensitivities [8]. The most promising species among those presented in Chart 1 is the long-sought pentazolate anion cyclo-N5– [14]. Though this anion was first observed in 2002, its synthesis in the form of an ammonium salt was only reported in 2017 [15,16]. A synthetic route to various salts of the cyclo-N5– anion has now been developed, making this hitherto elusive species much more readily available [9]. We should stress that the crystal structure of NH4+ cyclo-N5– was predicted earlier using an evolutionary algorithm, and it was calculated that this salt should be thermodynamically stable at pressures above 30 GPa [17]. Despite the wrong space group, the predicted crystalline environment is impressively close to that determined experimentally [9]. This clearly demonstrates the predictive force of modern theoretical methods.
Steele and Oleynik performed an extensive study of ternary C–N–O and H–N–O systems at 50 and 200 GPa using an evolutionary algorithm, and obtained the corresponding phase diagrams [[18], [19], [20], [21]]. A few interesting compounds, namely hydrazinium hydroxide (N2H5)(HO)–P21/m, diammonium oxide (NH4)2O-Cmcm, and nitric acid (HNO3)–P21/m, were found to be stable at 200 GPa [18]. At 50 GPa, six ternary compounds were located on the convex hull: H10N2O–C2/m, H8N4O- , HNO3–P21/m, H6N2O8- , H8N2O- , and H12N2O3–Cm [18]. Bogdanova et al. [22] calculated detonation parameters of different condensed high explosives of the general composition HxNyOz using a multiphase model of detonation products based on the equations of state.
The aforementioned studies inspired us to model a molecular form of an HxNyOz system having zero oxygen balance (H2O)xNy, thus releasing only H2O and N2, relatively high nitrogen content (>50 wt%), existing as relatively small molecules, having a high positive enthalpy of formation and crystal density, while maintaining dynamic, mechanical, and thermal stability. Moreover, as we have recently found, the novel energetic polymeric material CarNit4, poly(1,5-tetrazolediyl), demonstrates excellent potential as an explosive, but suffers from a lack of internal oxidant [23]. Thus, we have also modeled a positive oxygen balance material of the composition HxNyOz for potential application as a solid oxidant in explosive formulations or multi-component propellants, alternatives to various heterocyclic compounds [24], or strained nitro-triaziridine derivatives [[25], [26], [27]]. Herein, we present our crystal structure prediction and characterization with respect to dynamic, mechanical, and thermal stability of two crystalline materials, 1,4,2,3,5,6-dioxatetrazinane (DOTZ) and its 2,5-dinitro derivative (DNDOTZ), as well as information for their future spectral identification.
Section snippets
Computational details
The first-principles calculations employed in this work were performed in terms of Density Functional Theory (DFT) within the generalized gradient approximation (GGA) using the Materials Studio 7.0 suite of programs [28]. Geometry optimizations and calculations of band structure (BS), phonons, and vibrational spectra, as well as optical properties and electron excitation energies, were performed with the Cambridge Serial Total Energy Package (CASTEP) code [29]. Elastic constants and lattice
Structural features and topological analysis
The structures of the studied compounds (DOTZ and DNDOTZ) in molecular and crystalline forms are illustrated in Fig. 1a and b. In this work, we applied our modified eigenvector-following scheme to find crystal structures (Fig. S1 in the Supplementary data). This method avoids the optimization of a huge number of crystal structures, which are often very odd or improbable (as with purely automatic algorithms), and restrict the number of structures to just a few. Obviously, one can never be 100%
Conclusions
In summary, the presented results reveal that DOTZ and DNDOTZ are stable in the crystalline state (with space group ). This was corroborated by corresponding calculations of phonon dispersions and elastic constants as well as molecular dynamics simulations. DOTZ exhibits excellent detonation properties and outperforms all hitherto known explosives, both those obtained experimentally and those predicted theoretically. In comparison, DNDOTZ shows inferior detonation performance, mainly due to
CRediT authorship contribution statement
Sergey V. Bondarchuk: Conceptualization, Methodology, Software, Data curation, Writing - original draft, Visualization, Investigation, Writing - review & editing.
Declaration of competing interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Acknowledgments
This work was supported by the Ministry of Education and Science of Ukraine Research Fund (Grant No. 0118U003862).
References (68)
- et al.
Two-dimensional honeycomb (A7) and zigzag sheet (ZS) type nitrogen monolayers. A first principles study of structural, electronic, spectral, and mechanical properties, Comput
Mater. Sci.
(2017) - et al.
Modeling of ammonium dinitramide (ADN) monopropellant combustion with coupled condensed and gas phase kinetics
Combust. Flame
(2014) - et al.
Cyclic triamines as potential high energy materials. Thermochemical properties of triaziridine and triazirine
J. Mol. Struct.: THEOCHEM
(2005) - et al.
Time-dependent density functional theory within the Tamm-Dancoff approximation
Chem. Phys. Lett.
(1999) - et al.
Critic2: a program for real-space analysis of quantum chemical interactions in solids
Comput. Phys. Commun.
(2014) Impact sensitivity of aryl diazonium chlorides: limitations of molecular and solid-state approach
J. Mol. Graph. Model.
(2019)- et al.
High-throughput electronic band structure calculations: challenges and tools
Comput. Mater. Sci.
(2010) - et al.
Synthesis and properties of 2,3-dimethoxy-1,4,2,3-dioxadiazinane and dialkoxydiazene oxides
Mendeleev Commun.
(1992) Explosives – market analysis, trends and forecasts
Beyond molecular nitrogen: revelation of two ambient-pressure metastable single- and double-bonded nitrogen allotropes built from three-membered rings
Phys. Chem. Chem. Phys.
(2019)
Super high-energy density single-bonded trigonal nitrogen allotrope—a chemical twin of the cubic gauche form of nitrogen
Phys. Chem. Chem. Phys.
Single-bonded cubic form of nitrogen
Nat. Mater.
Pressure-induced symmetry-lowering transition in dense nitrogen to layered polymeric nitrogen (LP-N) with colossal Raman intensity
Phys. Rev. Lett.
Transformation of hydrazinium azide to molecular N8 at 40 GPa
J. Chem. Phys.
A safe and scaled up route to inert ammonia oxide hydroxylammonium azide (H7N5O2), hydrazinium azide (H5N5) and ammonium azide (H4N4)
ACS Appl. Energy Mater.
A series of energetic cyclo-pentazolate salts: rapid synthesis, characterization, and promising performance
J. Mater. Chem. A
Synthesis and characterization of cyclo-pentazolate salts of NH4+, NH3OH+, N2H5+, C(NH2)3+, and N(CH3)4+
J. Am. Chem. Soc.
A new class of flexible energetic salts, part 6: the structures of the hydrazinium and hydroxylammonium salts of dinitramide
J. Chem. Crystallogr.
1,1,3,3-Tetraoxo-1,2,3-triazapropene anion, a new oxy anion of nitrogen: the dinitramide anion and its salts
J. Am. Chem. Soc.
Recent advances in the chemistry of N5+, N5– and high-oxygen compounds
Propellants, Explos. Pyrotech.
Synthesis and characterization of the pentazolate anion cyclo-N5– in (N5)6(H3O)3(NH4)4Cl
Science
A series of energetic metal pentazolate hydrates
Nature
Pentazole and ammonium pentazolate: crystalline hydro-nitrogens at high pressure
J. Phys. Chem. A
Computational Discovery of Energetic Polynitrogen Compounds at High Pressure
Computational discovery of new high-nitrogen energetic materials
New phase of ammonium nitrate: a monoclinic distortion of AN-IV
J. Chem. Phys.
Ternary inorganic compounds containing carbon, nitrogen, and oxygen at high pressures
Inorg. Chem.
Thermodynamic modelling of detonation H-N-O high explosives
J. Phys.: Conf. Ser.
Investigation of Oxygen- and Nitrogen-Rich Heterocyclic Compounds as Potential High-Energy Dense Oxidizers or Secondary Explosives
Density functional theory study on energetic nitro-triaziridine derivatives
Struct. Chem.
Ab initio quantum chemical computations of substituent effects on triaziridine strain energy and heat of formation
Phys. Chem. Chem. Phys.
Materials Studio 7.0
First principles methods using CASTEP
Z. für Kristallogr. - Cryst. Mater.
Cited by (13)
Interpol review of the analysis and detection of explosives and explosives residues
2023, Forensic Science International: SynergyModulation of ADN Crystal Surface Properties by Additives
2024, Industrial and Engineering Chemistry Research