Modeling of explosives: 1,4,2,3,5,6-dioxatetrazinane as a new green energetic material with enhanced performance

https://doi.org/10.1016/j.jpcs.2020.109458Get rights and content

Highlights

  • Two novel carbon-free crystalline energetic compounds are predicted and characterized.

  • One of them represents a powerful green explosive, yielding only N2 and H2O as detonation products.

  • Its detonation properties exceed those of most explosives hitherto synthesized or predicted.

  • The other compound is a very effective solid oxidant, which outperforms gaseous oxygen.

Abstract

Crystal structure prediction and characterization data have been computed for two novel carbon-free energetic compounds, namely 1,4,2,3,5,6-dioxatetrazinane (DOTZ) and its 2,5-dinitro derivative (DNDOTZ). Dynamic and mechanical stability of these crystalline materials at ambient pressure has been confirmed. Ab initio molecular dynamics simulations confirmed thermal stability of DOTZ at 298 K and of DNDOTZ at 77 K. A comprehensive spectral characterization (IR, Raman, UV/Vis, NMR) has been performed to provide information for future experimental identification. The parent molecule exhibits powerful detonation properties, which exceed those of most compounds hitherto obtained experimentally or predicted theoretically. Meanwhile, DNDOTZ is predicted to be an effective solid oxidant for common explosives with negative oxygen balance, such as TNT, TNB, DATB, and TATB. Such mixtures demonstrate better propulsive properties than corresponding mixtures of these explosives with gaseous oxygen. By virtue of its zero oxygen balance, the only detonation products of DOTZ are environmentally friendly molecular nitrogen and water. Therefore, this compound is of great potential interest for further energetic applications.

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-R3m , HNO3P21/m, H6N2O8-P1 , H8N2O-P3m1 , and H12N2O3Cm [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 P1). 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)

  • S.V. Bondarchuk et al.

    Super high-energy density single-bonded trigonal nitrogen allotrope—a chemical twin of the cubic gauche form of nitrogen

    Phys. Chem. Chem. Phys.

    (2017)
  • M.I. Eremets et al.

    Single-bonded cubic form of nitrogen

    Nat. Mater.

    (2004)
  • D. Tomasino et al.

    Pressure-induced symmetry-lowering transition in dense nitrogen to layered polymeric nitrogen (LP-N) with colossal Raman intensity

    Phys. Rev. Lett.

    (2014)
  • S. Duwal et al.

    Transformation of hydrazinium azide to molecular N8 at 40 GPa

    J. Chem. Phys.

    (2018)
  • G. Zhao et al.

    A safe and scaled up route to inert ammonia oxide hydroxylammonium azide (H7N5O2), hydrazinium azide (H5N5) and ammonium azide (H4N4)

    ACS Appl. Energy Mater.

    (2019)
  • Y. Xu et al.

    A series of energetic cyclo-pentazolate salts: rapid synthesis, characterization, and promising performance

    J. Mater. Chem. A

    (2019)
  • C. Yang et al.

    Synthesis and characterization of cyclo-pentazolate salts of NH4+, NH3OH+, N2H5+, C(NH2)3+, and N(CH3)4+

    J. Am. Chem. Soc.

    (2018)
  • R. Gilardi et al.

    A new class of flexible energetic salts, part 6: the structures of the hydrazinium and hydroxylammonium salts of dinitramide

    J. Chem. Crystallogr.

    (2000)
  • J.C. Bottaro et al.

    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.

    (1997)
  • K.O. Christe

    Recent advances in the chemistry of N5+, N5 and high-oxygen compounds

    Propellants, Explos. Pyrotech.

    (2007)
  • Ch Zhang et al.

    Synthesis and characterization of the pentazolate anion cyclo-N5 in (N5)6(H3O)3(NH4)4Cl

    Science

    (2017)
  • Y. Xu et al.

    A series of energetic metal pentazolate hydrates

    Nature

    (2017)
  • B.A. Steele et al.

    Pentazole and ammonium pentazolate: crystalline hydro-nitrogens at high pressure

    J. Phys. Chem. A

    (2017)
  • B.A. Steele

    Computational Discovery of Energetic Polynitrogen Compounds at High Pressure

    (2018)
  • B.A. Steele et al.

    Computational discovery of new high-nitrogen energetic materials

  • B.A. Steele et al.

    New phase of ammonium nitrate: a monoclinic distortion of AN-IV

    J. Chem. Phys.

    (2015)
  • B.A. Steele et al.

    Ternary inorganic compounds containing carbon, nitrogen, and oxygen at high pressures

    Inorg. Chem.

    (2017)
  • YuA. Bogdanova et al.

    Thermodynamic modelling of detonation H-N-O high explosives

    J. Phys.: Conf. Ser.

    (2016)
  • S.V. Bondarchuk, CarNit4 – A New Polymeric Energetic Material Based on Poly(1,5-Tetrazolediyl), Polymer, accepted...
  • T.S. Hermann

    Investigation of Oxygen- and Nitrogen-Rich Heterocyclic Compounds as Potential High-Energy Dense Oxidizers or Secondary Explosives

    (2018)
  • W. Chi et al.

    Density functional theory study on energetic nitro-triaziridine derivatives

    Struct. Chem.

    (2013)
  • R. Peverati et al.

    Ab initio quantum chemical computations of substituent effects on triaziridine strain energy and heat of formation

    Phys. Chem. Chem. Phys.

    (2009)
  • Materials Studio 7.0

    (2013)
  • S.J. Clark et al.

    First principles methods using CASTEP

    Z. für Kristallogr. - Cryst. Mater.

    (2005)
  • Cited by (13)

    • Modulation of ADN Crystal Surface Properties by Additives

      2024, Industrial and Engineering Chemistry Research
    View all citing articles on Scopus
    View full text