Designing compounds with high energetic performance and low sensitivity at the molecular level is the central challenge in molecular design for energetic materials. Molecular frameworks integrating 1,2,4-triazolo [4,3-b] [1,2,4,5] tetrazine (known as high-energy-density and low-sensitivity module) and monocyclic tetrazoles are constructed to design energetic compounds. A series of energetic compounds are constructed. Density functional theory (DFT) has been used to investigate geometries, frontier molecular orbital energy, heats of formation (HOFs), densities of the title compounds at B3PW91/6-31G (d, f) level. Heats of formation were calculated via isodesmic reactions. Crystal densities were predicted using Politzer's method. Detonation velocity (D) and detonation pressure (P) of the title compounds have been determined based on HOFs and densities through the Kamlet-Jacobs (K-J) equation. The effects of substituents on above properties are presented. Substitution with -NO2, -ONO2, -NH2, -NHNO2, and -N(NO2)2 groups can increase the heats of formation, densities, detonation velocity and detonation pressure of the compounds. Specifically, the compound substituted with the -N(NO2)2 group exhibits higher detonation performance than the high-energy explosive RDX, suggesting its potential as a promising high energy material. The results are of significant value, providing theoretical guidance for the molecular design of novel high-energy-density compounds and the optimization of established ones.
This is an Open Access article, distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution and reproduction in any medium or format, provided the original work is properly cited.
1,2,4-triazolo [4,3-b][1,2,4,5] Tetrazine, Density Functional Theory, Heats of Formation, Detonation Velocity, Detonation Pressure
1. Introduction
The development of high-energy, insensitive and green energetic compounds has always been a frontier issue in energetic materials research. How to balance the contradiction between high energy and low sensitivity at the micro level of molecular structure remains the core challenge in energetic material molecular design.
[1]
Tang, J., Yang, P. J., Yang, H. W., Xiong, H. L., Hu, W., Cheng, G. B. A Simple and Efficient Method to Synthesize High-nitrogen Compounds: Incorporation of Tetrazole Derivatives with N5 chains, Chemical Engineering Journal. 2020, 386(124027) 1-13.
Fan, H. H., Tang, J., Hu, W., Zheng, X. X., Yang, P. J., Cheng, G. B., Xiao, C., Yang, H. W., Combination of Tetrazole and 4-Azido-pyrazolotriazine Oxide: Balance of High Nitrogen, Energy, and Safety, Organic Letter, 2025, 27(3) 846-850.
Chen, B. H, Lu, H., Chen, J. Y., Chen, Z. X., Yin, S. F., Peng, L. F., Qiu, R. H. Recent Progress on Nitrogen-Rich Energetic Materials Based on Tetrazole Skeleton. Topics in Current Chemistry, 2023, 381(5), 2365-0869.
Cui, Z. Y., Wang, J. H., Wu, L. L., Yuan, X., Yu, Q., Zhang, J. Lin, K. F., Yang, Y. L., Xia, D. B., Strategy for Balance Energy and Safety: Salt Formation of Nitrogen-Rich Bicyclic Compounds Based on 1,2,4-Triazole, Crystal Growth & Design, 2025, 25(1), 88-100.
Nitrogen-rich fused-ring energetic compounds exhibit superior structural and performance characteristics compared to monocyclic and linked-ring nitrogen-rich heterocycles
[5]
Hu, L., He, C, L., Zhao, G., Imler, G. H., Parrish, D. A., Shreeve, J. M. Selecting Suitable Substituents for Energetic Materials Based on a Fused Triazolo-[1,2,4,5]tetrazine Ring, ACS Applied Energy Materials. 2020, 3(6), 5510-5516.
Liu, Y., Zhao, G., Tang, Y., Zhang, J., Hu, L., Imler, G. H., Parrish, D. A., Shreeve, J. M. Multipurpose [1,2,4] Triazolo[4,3-b] [1,2,4,5] Tetrazine-based Energetic Materials, Journal of Materials Chemistry A. 2019, 7(13), 7875-7884.
. Structurally, fused-ring systems are characterized by more energetic bonds (e.g., N=N, N–N, C–N, and N–O) and greater ring strain, which collectively enhance their thermochemical and detonation properties. These features result in significantly higher heats of formation (HOFs), greater crystal densities, and improved detonation velocities/pressures relative to monocyclic or linked-ring analogs. Furthermore, the coplanar structure of fused-ring systems facilitates extensive delocalized π-electron delocalization across the conjugated framework and promotes intermolecular π-π stacking interactions. These structural attributes not only reduce mechanical sensitivity (e.g., impact sensitivity, friction sensitivity) but also confer enhanced thermal stability (high decomposition temperatures). Collectively, these advantages establish nitrogen-rich fused-ring compounds as competitive candidates for next-generation insensitive high-energy materials (IHEMs).
Experimental results demonstrate that the fused 1,2,4-triazolo[4,3-b][1,2,4,5]tetrazine ring system exhibits superior insensitivity, with an impact sensitivity (IS) value exceeding 40 J and a friction sensitivity (FS) value surpassing 360 N
[1]
Tang, J., Yang, P. J., Yang, H. W., Xiong, H. L., Hu, W., Cheng, G. B. A Simple and Efficient Method to Synthesize High-nitrogen Compounds: Incorporation of Tetrazole Derivatives with N5 chains, Chemical Engineering Journal. 2020, 386(124027) 1-13.
. These values indicate markedly lower sensitivity compared to conventional explosives such as 1,3,5-trinitro-1,3,5-triazacyclohexane (RDX) (IS = 7.5 J, FS = 120 N) and 2,4,6-trinitrotoluene (TNT) (IS = 15 J, FS = 353 N). Thereby the fused 1,2,4-triazolo [4,3-b] [1,2,4,5]tetrazine ring system combines high energy with insensitivity.
Scheme 1. N-(2H-Tetrazol-5-yl)-[1,2,4]triazolo[4,3-b][1,2,4,5] Tetrazin-6-amine and their derivatives considered in this study.
T0: R=H T1: R=NO2 T2: R=ONO2 T3: R=CH3
T4: R= OCH3 T5: R=CN T6: R=N3 T7: R=NH2
T8: R= NHNO2 T9: R= N(NO2)2
In this project, we employ the 1,2,4-triazolo[4,3-b][1,2,4,5] tetrazine fused ring system coupled with nitrogen-rich monocyclic tetrazole units as fundamental building blocks to design derivatives that balance high energy and insensitivity. Density functional theory (DFT) are used to investigate the energy characteristics and stability of these compounds, aiming to screen high-energy and insensitive energetic materials. This research will provide theoretical guidance for the design and performance optimization of energetic materials.
2. Theory and Computational Details
The designed 10 molecules are shown in Scheme 1. Calculations were carried out with the Gaussian 09 program package
[7]
Frisch, M. J., Trucks, G. W., Schlegel, H. B., et al. Gaussian 09, Gaussian, Inc., Wallingford CT, 2009.
[7]
at B3PW91/6-31G (d, f) level. The geometries of molecules, including the title compounds and the reference species mentioned below to calculate HOFs for the title compounds are optimized without any constraints. Vibrational analysis has been performed for each stationary point to verify a minimum energy structure and to provide zero-point energies (ZPEs) and thermal correction (HT).
2.1. Heats of Formation
The standard heats of formation (HOFs) of a compound at 298.15 K is defined as the enthalpy change when one mole of the substance is formed from its constituent elements in their standard states under standard thermodynamic conditions (1 bar pressure). To compute HOFs at 298.15 K, isodesmic reaction methods have been used. Isodesmic reactions maintain the types and numbers of chemical bonds unchanged during the reaction process. Typically, they allow systematic errors arising from factors such as spin contamination and basis set superposition error (BSSE) in the system's bonding environment to cancel each other out, thus enabling the accurate determination of the standard heat of formation.
[8]
Hahre, W. J., Radom, L., Schleyer, P. V. R. Ab initio molecular orbital theory, New York: Wiley; 1986.
[8]
.
2.2. Density
Politzer's crystal electron structure theory, which analyzes the molecular surface electrostatic potential distribution and van der Waals interaction parameters, enhances the accuracy of density predictions for energetic crystals (mean absolute error < 0.03 g cm-3)
[9]
Politzer, P., Martinez, J., Murray, J. S., Concha, M. C., Toro-Labbé, A. An Electrostatic Interaction Correction for Improved crystal Density Prediction, Molecular Physics, 2009, 107(19) 2095–2101.
. In this study, we employ this approach shown below to predict the densities (g cm-3) of our study derivatives with consideration of the intermolecular interactions.
(1)
where M is the molecular mass (g mol-1) and Vm represents the molecular volume (cm3 mol-1) of an isolated gas-phase molecule, defined as the spatial volume enclosed within the 0.001 electrons bohr-3 electron density isosurface. ν characterizes the equilibrium ratio between positive and negative electrostatic potential surfaces, represents the variability measure of the electrostatic potential. α, β and γ are the coefficients assigned through fitting eq (1) to the experimental densities of a range of 36 energetic compounds. They are 0.9183, 0.0028, and 0.0443, respectively
[9]
Politzer, P., Martinez, J., Murray, J. S., Concha, M. C., Toro-Labbé, A. An Electrostatic Interaction Correction for Improved crystal Density Prediction, Molecular Physics, 2009, 107(19) 2095–2101.
was used to calculate surface electrostatic potential quantities.
2.3. Detonation Parameters
Detonation velocity (D) and detonation pressure (P) are two critical parameters for assessing the explosive performance of energetic materials. The Kamlet-Jacobs (K-J) equation is a well-established computational method originally developed by Kamlet M. J. and Jacobs S. J. to predict energetic material performance
[11]
Kamlet, M. J., Jacobs, S. J. Chemistry of detonations. I. A Simple Method for Calculating Detonation Properties of C–H–N–O Explosives, Journal of Chemical Physics. 1968, 48(1), 23–35.
. In this study, the K-J equations were employed to predict the detonation velocity and detonation pressure of the derivatives. The calculation formulas are as follows:
(2)
(3)
where D is detonation velocity in km s−1, P is detonation pressure in GPa, also known as Chapman-Jouguet (C-J) pressure, N is molar quantity of gaseous detonation products per unit mass in mol g−1, M is the average molar mass of detonation gas products in g mol−1, Q is chemical energy released per unit mass during detonation, ρ is the crystal density in g cm−3.
2.4. Frontier Molecular Orbital Energy
According to the Frontier Molecular Orbital (FMO) theory, the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) are key determinants of a system's chemical reactivity. The molecular orbital energy gap ΔE (ΔE = ELUMO - EHOMO) between the HOMO and the LUMO represents the energy required for electron transition. The higher the energy of HOMO and the lower the energy of LUMO of a compound, the smaller the energy gap ΔE, and the higher the chemical reactivity of the compound.
3. Results and Discussion
3.1. Heats of Formation
In designing isodesmic reactions, [1,2,4] triazolo [4,3-b] [1,2,4,5] Tetrazine and 2H-tetrazole backbone were selected as the reference framework to preserve bond types and quantities, thereby reducing computational errors. The experimental and calculated gas phase HOFs of all reference compounds are listed in Table 1. The isodesmic reactions designed are displayed in Table 2.
Table 1. The gas phase heat of formation () of the reference compounds.
Compd.
Compd.
[1,2,4] triazolo[4,3-b]-[1,2,4,5] Tetrazine
683.16a
2H-tetrazole
327.35a
CH4
-74.6 ± 0.3b
NH2ONO2
29.90a
C6H6
82.9 ± 0.9b
NH2OCH3
-31.33a
NH3
-45.9b
NH2CN
138.41a
CH3NHCH3
-19.0 ± 2.0b
NH2N3
421.12a
NH2NO2
-2.38a
NH2NHNO2
99.11a
NH2CH3
-23.5a
NH2N(NO2)2
181.80a
NH2NH2
95.35a
a The values are calculated at the G2 level from the atomization reaction.
b The experimental HOFs were taken from reference
[12]
NIST Chemistry WebBook; Linstrom, P. J., Mallard, W. G., Eds.; NIST Standard Reference Database Number 69; National Institute Standards and Technology: Gaithersburg MD, 2005.
c are the gas phase heat of formation in kJ mol-1.
For the isodesmic reactions, the heats of reaction at 298.15K can be calculated through the reaction enthalpies:
(4)
where and are the sums of the heats of formation for products and reactants in gas at 298.15K, respectively. Since the HOFs of reference compounds are available as shown in Table 1, the HOFs of the title compounds can be obtained if the heats of reaction are known. The can be calculated from the following equations:
(5)
where ΔEe and ΔZPE are the electronic energy difference and the zero-point energy difference between products and reactants at 0K, respectively; ΔHT is the changes in thermal correction to enthalpies between products and reactants; Δ(PV) equals ΔnRT for reaction in gas phase. For isodesmic reactions, Δn=0.
Table 2. Designed isodesmic work reactions for calculating gas heats of formation of the title compounds.
The solid phase heat of formation is required for determining the energetic properties of energetic materials. According to Hess's Law, the solid phase heat of formation can be predicted from the gas-phase heat of formation and the heat of sublimation enthalpy as,
(6)
(7)
where A is the isosurface area corresponding to an electron density of 0.001 electrons/Bohr3, ν and are defined and provided in the density calculation methodology (Section 2.2). a, b, and c are the coefficients optimized via iterative least-squares minimization. They were determined by the Byrd and Rice method
[13]
Byrd, E. F. C., Rice, B. M. Improved Prediction of Heats of For- mation of Energetic Materials Using Quantum Mechanical Calculations, Journal of Physical Chemistry A. 2006, 110(3), 1005–1013.
and are 2.670×10-4 kcal mol-1 Å-4, 1.650 kcal mol-1, and 2.966 kcal mol-1, respectively.
Calculated total energies (E0), Zero point Energies (ZPE), thermal corrections (HT), heats of sublimation ( ) and gas/solid phase heats of formation ( and ) of the title compoundsare listed in Table 3.
Table 3. Calculated total energies (E0), zero point energies (ZPE), thermal corrections (HT), heats of sublimation () and heats of formation (and ) of the title compounds*
Compd.
E0
ZPE
HT
T0
-756.0581
0.1151
29.31
994.45
27.22
880.55
T1
-960.4292
0.1159
36.33
1119.90
28.84
999.24
T2
-1035.5703
0.1198
39.10
1089.05
31.02
959.26
T3
-795.3635
0.1429
33.75
976.92
27.74
860.87
T4
-870.4824
0.1465
36.53
998.76
29.10
877.02
T5
-848.2194
0.1122
34.01
1249.21
28.36
1130.55
T6
-919.5298
0.1167
36.67
1473.05
29.52
1349.54
T7
-811.3607
0.1315
33.44
1091.91
29.00
970.57
T8
-1015.7535
0.1337
38.89
1173.36
31.52
1041.49
T9
-1220.1198
0.1339
46.50
1283.89
34.89
1137.92
*E0 and ZPE are in a.u., HT, and HOF are in kJ mol-1.
The solid phase heats of formation of different substituted compounds are shown in Figure 1. As shown in Figure 1, the solid phase HOF of compounds T1, T2, T5, T6, T7, T8, and T9 are all greater than that of T0. This indicates that substitutions with -NO2, -ONO2, -CN, -N₃, -NH2, -NHNO2, and -N(NO2)2 groups increase the standard heat of formation compared to the parent compound, except for CH₃ and OCH₃ groups. Notably, the -N₃ groups significantly enhance the HOF. Furthermore, comparing the HOF values of T0, T1, and T7–T9 reveals that the heat of formation increases with the number of nitro (-NO2) groups. Furthermore, all compounds exhibit heats of formation exceeding that of HMX (105 kJ mol-1), indicating that the fused-ring framework elevates the heat of formation, and consequently enhances the performance of energetic materials.
Figure 1. Comparison of solid phase heats of formation of different substituted compounds.
3.2. Density
The calculated crystal densities of the title compounds are listed in Table 4. To visually analyze substituent effects on density, Figure 2 was plotted.
It can be seen that substituted compounds exhibit significant density variations. Specifically, the densities of -NO2, -ONO2, -NHNO2, and -N(NO2)2 substituted compounds are 1.840 g cm⁻3, 1.867 g cm⁻3, 1.838 g cm⁻3, and 1.901 g cm⁻3, respectively. All values exceed that of unsubstituted compound T0 (1.799 g cm⁻3) and surpass the density of explosive RDX (1.82 g cm⁻3), demonstrating significant density enhancement. This phenomenon is attributed to intra- or intermolecular O···H hydrogen bonding between hydrogen atom and nitro-group oxygen atom. Critically, the -N(NO2)2 functionalized compound achieves a density of 1.901 g cm⁻3, surpassing HMX (1.900 g cm⁻3). In contrast, -CH₃, -OCH₃, -CN, -N₃, and -NH2 substituted compounds exhibit reduced densities.
3.3. Detonation Performance and HOMO-LUMO Energy Gap
Detonation velocity and pressure are critical parameters for evaluating the performance of energetic compounds. The values of the title compounds calculated via the K-J equations are tabulated in Table 4. The HOMO-LUMO energy gap is correlated with the sensitivity of energetic materials
[14]
Xu, X. J., Zhu, W. H., Xiao, H. M. DFT Studies on the Four Polymorphs of Crystalline CL-20 and the Influences of Hydrostatic Pressure on Epsilon-CL-20 Crystal, Journal of Physical Chemistry B. 2007, 111(8) 2090-2097.
, and the calculated values are also listed in Table 4. Substituents differentially affect detonation performance: -CH₃, -OCH₃ and -CN groups decrease D and P values, whereas -NO2, -ONO2, -N₃, -NH2, -NHNO2, and -N(NO2)2 groups enhance these detonation parameters. Although the energetic properties of -NO2、-ONO2、-N3、NHNO2 and N(NO2)2 substituted compounds remain inferior to HMX benchmarks (ρ = 1.90 g cm⁻3, D = 9.370 km s⁻¹, P = 40.20 GPa), they exhibit considerably higher ρ, D, and P values, which are compared to TNT (ρ = 1.65 g cm⁻3, D = 7.129 km s⁻¹, P = 21.37 GPa), demonstrating potential as novel energetic materials. Notably, the dinitramino-substituted compound T9 exhibits higher ρ, D, and P values than high-energy explosive RDX (ρ = 1.77 g cm⁻3, D = 8.882 km s⁻¹, P = 34.64 GPa), along with a larger HOMO-LUMO gap compared to its unsubstituted compound T0. These characteristics fulfill dual requirements for high-performance and low-sensitivity energetic materials, positioning it as a promising superior energetic material.
Figure 3. Detonation velocity and detonation pressure of different substituted compounds.
Table 4. Calculated density, detonation performance and molecular orbital energy gap ΔE for designed title compounds*
Compd.
D
P
Is
ΔE
T0
1.799
7.490
24.875
0.868
3.82
T1
1.840
8.291
30.896
0.947
3.95
T2
1.867
8.552
33.147
0.960
3.86
T3
1.677
6.919
20.335
0.852
3.73
T4
1.703
7.351
23.173
0.897
3.76
T5
1.751
7.371
23.702
0.909
3.97
T6
1.779
7.973
27.999
0.957
3.64
T7
1.781
7.610
25.522
0.886
3.65
T8
1.838
8.328
31.149
0.943
3.92
T9
1.901
8.927
36.504
0.984
3.96
*ρ is crystal density in g cm-3, D is detonation velocity in km s-1, P is detonation pressure in GPa and ΔE is molecular orbital energy gap in eV.
4. Conclusions
Ten nitrogen-rich fused-ring energetic compounds were constructed based on N-(2H-Tetrazol-5-yl)-[1,2,4]triazolo [4,3-b] [1,2,4,5] tetrazin-6-amine. The DFT-B3PW91 method was used to study the HOMO-LUMO energy gaps, heats of formation, densities and detonation performance.
Heats of formation were calculated via isodesmic reactions, crystal densities were predicted using Politzer's method, and detonation velocities/pressures were calculated using K-J equations. Key findings are summarized as follows:
1. Exceptfor -CH₃ and -OCH₃ groups, substitution with -NO2, -ONO2, -CN, -N₃, -NH2, -NHNO2, and -N(NO2)2 groups can increase the heats of formation of the compounds.
2. The densities of compounds substituted with -NO2, -ONO2, -NHNO2, and -N(NO2)2 groups are 1.840 g·cm⁻3, 1.867 g·cm⁻3, 1.838 g·cm⁻3, and 1.901 g·cm⁻3, respectively. These values are all higher than that of the unsubstituted compound T₀ (1.799 g·cm⁻3) and exceed the density of the explosive RDX (1.82 g·cm⁻3), indicating a significant improvement in compound density.
3. Detonation performance results show that substitution with -NO2, -ONO2, -N₃, -NH2, -NHNO2, and -N(NO2)2 groups improves the detonation performance. Specifically, the compound substituted with the -N(NO2)2 group exhibits higher values of density, detonation velocity, and detonation pressure than the high-energy explosive RDX. Additionally, its HOMO-LUMO energy gap is larger than that of the unsubstituted compound, suggesting its potential as a promising high-energy material.
Abbreviations
DFT
Density Functional Theory
HOF
Heat of Formation
D
Detonation Velocity
P
Detonation Pressure
HOMO
Highest Occupied Molecular Orbital
LUMO
Lowest Unoccupied Molecular Orbital
Author Contributions
Xinfang Su is the sole author. The author read and approved the final manuscript.
Data Availability Statement
The data supporting the outcome of this research work has been reported in this manuscript.
Funding
This work is supported by Science Foundation of Beijing University of Civil Engineering and Architecture (ZF17073, ZD202506).
Conflicts of Interest
The authors declare no conflicts of interest.
References
[1]
Tang, J., Yang, P. J., Yang, H. W., Xiong, H. L., Hu, W., Cheng, G. B. A Simple and Efficient Method to Synthesize High-nitrogen Compounds: Incorporation of Tetrazole Derivatives with N5 chains, Chemical Engineering Journal. 2020, 386(124027) 1-13.
Fan, H. H., Tang, J., Hu, W., Zheng, X. X., Yang, P. J., Cheng, G. B., Xiao, C., Yang, H. W., Combination of Tetrazole and 4-Azido-pyrazolotriazine Oxide: Balance of High Nitrogen, Energy, and Safety, Organic Letter, 2025, 27(3) 846-850.
Chen, B. H, Lu, H., Chen, J. Y., Chen, Z. X., Yin, S. F., Peng, L. F., Qiu, R. H. Recent Progress on Nitrogen-Rich Energetic Materials Based on Tetrazole Skeleton. Topics in Current Chemistry, 2023, 381(5), 2365-0869.
Cui, Z. Y., Wang, J. H., Wu, L. L., Yuan, X., Yu, Q., Zhang, J. Lin, K. F., Yang, Y. L., Xia, D. B., Strategy for Balance Energy and Safety: Salt Formation of Nitrogen-Rich Bicyclic Compounds Based on 1,2,4-Triazole, Crystal Growth & Design, 2025, 25(1), 88-100.
Hu, L., He, C, L., Zhao, G., Imler, G. H., Parrish, D. A., Shreeve, J. M. Selecting Suitable Substituents for Energetic Materials Based on a Fused Triazolo-[1,2,4,5]tetrazine Ring, ACS Applied Energy Materials. 2020, 3(6), 5510-5516.
Liu, Y., Zhao, G., Tang, Y., Zhang, J., Hu, L., Imler, G. H., Parrish, D. A., Shreeve, J. M. Multipurpose [1,2,4] Triazolo[4,3-b] [1,2,4,5] Tetrazine-based Energetic Materials, Journal of Materials Chemistry A. 2019, 7(13), 7875-7884.
Frisch, M. J., Trucks, G. W., Schlegel, H. B., et al. Gaussian 09, Gaussian, Inc., Wallingford CT, 2009.
[8]
Hahre, W. J., Radom, L., Schleyer, P. V. R. Ab initio molecular orbital theory, New York: Wiley; 1986.
[9]
Politzer, P., Martinez, J., Murray, J. S., Concha, M. C., Toro-Labbé, A. An Electrostatic Interaction Correction for Improved crystal Density Prediction, Molecular Physics, 2009, 107(19) 2095–2101.
Kamlet, M. J., Jacobs, S. J. Chemistry of detonations. I. A Simple Method for Calculating Detonation Properties of C–H–N–O Explosives, Journal of Chemical Physics. 1968, 48(1), 23–35.
NIST Chemistry WebBook; Linstrom, P. J., Mallard, W. G., Eds.; NIST Standard Reference Database Number 69; National Institute Standards and Technology: Gaithersburg MD, 2005.
Byrd, E. F. C., Rice, B. M. Improved Prediction of Heats of For- mation of Energetic Materials Using Quantum Mechanical Calculations, Journal of Physical Chemistry A. 2006, 110(3), 1005–1013.
Xu, X. J., Zhu, W. H., Xiao, H. M. DFT Studies on the Four Polymorphs of Crystalline CL-20 and the Influences of Hydrostatic Pressure on Epsilon-CL-20 Crystal, Journal of Physical Chemistry B. 2007, 111(8) 2090-2097.
Su, X. (2025). Theoretical Investigation of Compounds Based on 1, 2, 4-Triazolo [4, 3-b] [1,2,4,5] Tetrazine as High Energy and Low Sensitivity Energetic Materials. Modern Chemistry, 13(2), 40-47. https://doi.org/10.11648/j.mc.20251302.12
Su, X. Theoretical Investigation of Compounds Based on 1, 2, 4-Triazolo [4, 3-b] [1,2,4,5] Tetrazine as High Energy and Low Sensitivity Energetic Materials. Mod. Chem.2025, 13(2), 40-47. doi: 10.11648/j.mc.20251302.12
Su X. Theoretical Investigation of Compounds Based on 1, 2, 4-Triazolo [4, 3-b] [1,2,4,5] Tetrazine as High Energy and Low Sensitivity Energetic Materials. Mod Chem. 2025;13(2):40-47. doi: 10.11648/j.mc.20251302.12
@article{10.11648/j.mc.20251302.12,
author = {Xinfang Su},
title = {Theoretical Investigation of Compounds Based on 1, 2, 4-Triazolo [4, 3-b] [1,2,4,5] Tetrazine as High Energy and Low Sensitivity Energetic Materials
},
journal = {Modern Chemistry},
volume = {13},
number = {2},
pages = {40-47},
doi = {10.11648/j.mc.20251302.12},
url = {https://doi.org/10.11648/j.mc.20251302.12},
eprint = {https://article.sciencepublishinggroup.com/pdf/10.11648.j.mc.20251302.12},
abstract = {Designing compounds with high energetic performance and low sensitivity at the molecular level is the central challenge in molecular design for energetic materials. Molecular frameworks integrating 1,2,4-triazolo [4,3-b] [1,2,4,5] tetrazine (known as high-energy-density and low-sensitivity module) and monocyclic tetrazoles are constructed to design energetic compounds. A series of energetic compounds are constructed. Density functional theory (DFT) has been used to investigate geometries, frontier molecular orbital energy, heats of formation (HOFs), densities of the title compounds at B3PW91/6-31G (d, f) level. Heats of formation were calculated via isodesmic reactions. Crystal densities were predicted using Politzer's method. Detonation velocity (D) and detonation pressure (P) of the title compounds have been determined based on HOFs and densities through the Kamlet-Jacobs (K-J) equation. The effects of substituents on above properties are presented. Substitution with -NO2, -ONO2, -NH2, -NHNO2, and -N(NO2)2 groups can increase the heats of formation, densities, detonation velocity and detonation pressure of the compounds. Specifically, the compound substituted with the -N(NO2)2 group exhibits higher detonation performance than the high-energy explosive RDX, suggesting its potential as a promising high energy material. The results are of significant value, providing theoretical guidance for the molecular design of novel high-energy-density compounds and the optimization of established ones.
},
year = {2025}
}
TY - JOUR
T1 - Theoretical Investigation of Compounds Based on 1, 2, 4-Triazolo [4, 3-b] [1,2,4,5] Tetrazine as High Energy and Low Sensitivity Energetic Materials
AU - Xinfang Su
Y1 - 2025/06/26
PY - 2025
N1 - https://doi.org/10.11648/j.mc.20251302.12
DO - 10.11648/j.mc.20251302.12
T2 - Modern Chemistry
JF - Modern Chemistry
JO - Modern Chemistry
SP - 40
EP - 47
PB - Science Publishing Group
SN - 2329-180X
UR - https://doi.org/10.11648/j.mc.20251302.12
AB - Designing compounds with high energetic performance and low sensitivity at the molecular level is the central challenge in molecular design for energetic materials. Molecular frameworks integrating 1,2,4-triazolo [4,3-b] [1,2,4,5] tetrazine (known as high-energy-density and low-sensitivity module) and monocyclic tetrazoles are constructed to design energetic compounds. A series of energetic compounds are constructed. Density functional theory (DFT) has been used to investigate geometries, frontier molecular orbital energy, heats of formation (HOFs), densities of the title compounds at B3PW91/6-31G (d, f) level. Heats of formation were calculated via isodesmic reactions. Crystal densities were predicted using Politzer's method. Detonation velocity (D) and detonation pressure (P) of the title compounds have been determined based on HOFs and densities through the Kamlet-Jacobs (K-J) equation. The effects of substituents on above properties are presented. Substitution with -NO2, -ONO2, -NH2, -NHNO2, and -N(NO2)2 groups can increase the heats of formation, densities, detonation velocity and detonation pressure of the compounds. Specifically, the compound substituted with the -N(NO2)2 group exhibits higher detonation performance than the high-energy explosive RDX, suggesting its potential as a promising high energy material. The results are of significant value, providing theoretical guidance for the molecular design of novel high-energy-density compounds and the optimization of established ones.
VL - 13
IS - 2
ER -
Su, X. (2025). Theoretical Investigation of Compounds Based on 1, 2, 4-Triazolo [4, 3-b] [1,2,4,5] Tetrazine as High Energy and Low Sensitivity Energetic Materials. Modern Chemistry, 13(2), 40-47. https://doi.org/10.11648/j.mc.20251302.12
Su, X. Theoretical Investigation of Compounds Based on 1, 2, 4-Triazolo [4, 3-b] [1,2,4,5] Tetrazine as High Energy and Low Sensitivity Energetic Materials. Mod. Chem.2025, 13(2), 40-47. doi: 10.11648/j.mc.20251302.12
Su X. Theoretical Investigation of Compounds Based on 1, 2, 4-Triazolo [4, 3-b] [1,2,4,5] Tetrazine as High Energy and Low Sensitivity Energetic Materials. Mod Chem. 2025;13(2):40-47. doi: 10.11648/j.mc.20251302.12
@article{10.11648/j.mc.20251302.12,
author = {Xinfang Su},
title = {Theoretical Investigation of Compounds Based on 1, 2, 4-Triazolo [4, 3-b] [1,2,4,5] Tetrazine as High Energy and Low Sensitivity Energetic Materials
},
journal = {Modern Chemistry},
volume = {13},
number = {2},
pages = {40-47},
doi = {10.11648/j.mc.20251302.12},
url = {https://doi.org/10.11648/j.mc.20251302.12},
eprint = {https://article.sciencepublishinggroup.com/pdf/10.11648.j.mc.20251302.12},
abstract = {Designing compounds with high energetic performance and low sensitivity at the molecular level is the central challenge in molecular design for energetic materials. Molecular frameworks integrating 1,2,4-triazolo [4,3-b] [1,2,4,5] tetrazine (known as high-energy-density and low-sensitivity module) and monocyclic tetrazoles are constructed to design energetic compounds. A series of energetic compounds are constructed. Density functional theory (DFT) has been used to investigate geometries, frontier molecular orbital energy, heats of formation (HOFs), densities of the title compounds at B3PW91/6-31G (d, f) level. Heats of formation were calculated via isodesmic reactions. Crystal densities were predicted using Politzer's method. Detonation velocity (D) and detonation pressure (P) of the title compounds have been determined based on HOFs and densities through the Kamlet-Jacobs (K-J) equation. The effects of substituents on above properties are presented. Substitution with -NO2, -ONO2, -NH2, -NHNO2, and -N(NO2)2 groups can increase the heats of formation, densities, detonation velocity and detonation pressure of the compounds. Specifically, the compound substituted with the -N(NO2)2 group exhibits higher detonation performance than the high-energy explosive RDX, suggesting its potential as a promising high energy material. The results are of significant value, providing theoretical guidance for the molecular design of novel high-energy-density compounds and the optimization of established ones.
},
year = {2025}
}
TY - JOUR
T1 - Theoretical Investigation of Compounds Based on 1, 2, 4-Triazolo [4, 3-b] [1,2,4,5] Tetrazine as High Energy and Low Sensitivity Energetic Materials
AU - Xinfang Su
Y1 - 2025/06/26
PY - 2025
N1 - https://doi.org/10.11648/j.mc.20251302.12
DO - 10.11648/j.mc.20251302.12
T2 - Modern Chemistry
JF - Modern Chemistry
JO - Modern Chemistry
SP - 40
EP - 47
PB - Science Publishing Group
SN - 2329-180X
UR - https://doi.org/10.11648/j.mc.20251302.12
AB - Designing compounds with high energetic performance and low sensitivity at the molecular level is the central challenge in molecular design for energetic materials. Molecular frameworks integrating 1,2,4-triazolo [4,3-b] [1,2,4,5] tetrazine (known as high-energy-density and low-sensitivity module) and monocyclic tetrazoles are constructed to design energetic compounds. A series of energetic compounds are constructed. Density functional theory (DFT) has been used to investigate geometries, frontier molecular orbital energy, heats of formation (HOFs), densities of the title compounds at B3PW91/6-31G (d, f) level. Heats of formation were calculated via isodesmic reactions. Crystal densities were predicted using Politzer's method. Detonation velocity (D) and detonation pressure (P) of the title compounds have been determined based on HOFs and densities through the Kamlet-Jacobs (K-J) equation. The effects of substituents on above properties are presented. Substitution with -NO2, -ONO2, -NH2, -NHNO2, and -N(NO2)2 groups can increase the heats of formation, densities, detonation velocity and detonation pressure of the compounds. Specifically, the compound substituted with the -N(NO2)2 group exhibits higher detonation performance than the high-energy explosive RDX, suggesting its potential as a promising high energy material. The results are of significant value, providing theoretical guidance for the molecular design of novel high-energy-density compounds and the optimization of established ones.
VL - 13
IS - 2
ER -