The reagents used in the study were provided by “Sigma-Aldrich, Fluka and Merck”. Melting temperature of the materials under study was measured by MELT-TEMP apparatus. The reaction steps and the final materials homogeneity were controlled by the thin-layer chromatography (PSL) method and silica gel on aluminum base plate (Merck).

The IR spectra were measured on KBr matrix by Perkin-Elmer spectrophotometer (FT-IR Spectrum BX II).

An LC-MS investigation was performed by Agilent technologies 6246 apparatus.

The general synthesis scheme is presented in Figure 1.

The structure of the molecules and their fragments under study were investigated by the density functional approach (DFT). We applied Becke’s three-parameter hybrid functional applying the non-local correlation provided by Lee, Yang, and Parr (B3LYP) [11] , a representative standard DFT method [12] . The cc-pVTZ basis set was also applied [13] . The approach was chosen taking into account

sufficient accuracy and comparably less computing time.

To obtain the most probable structure of the molecules, the Berny optimization was performed without any symmetry constraint of several different initial geometrical structures of the molecules with the same chemical compositions. The method we used enabled us to find the most stable geometrical structure. We also computed force constants and the resulting vibrational frequencies to be sure that equilibrium point was obtained.

To describe energetic properties of the compounds investigated, the detonation velocity and oxygen balance were calculated.

The detonation velocity was calculated as follows [14] :

where N is a number of the -NO₂ groups in the molecule. E is total energy of a molecule, a. u., and M is molar mass, g/mol. The regressed equations enabled us to get some fairly accurate ideas about detonation velocity values of explosives to compare their properties.

The oxygen balance was determined by using the following equation:

where M is molar mass, x―a number of C atoms, y―a number of H atoms, z―a number of O atoms [15] . Additionally, the toxicity and reactivity of the compounds investigated were obtained.

The transition states of the molecule under investigations were also obtained applying the Berny optimization procedure of the several different initial geometrical structures. The structures possessing the highest energy were selected for future investigations.

The IR spectra were obtained by using B3LYP/cc-pVTZ approach. The Gaussian program packages were applied [16] .

The views of the core molecule investigated are presented in Figure 2, the meaning of R1, R2 and R3 is presented in Table 1. Substitutions R1, R2 and R3 are -H or -NO₂. Here, the substitution ?NO₂ is introduced as ‘the additional NO₂ (nitro) group’ in order to separate it from other groups in the core molecule. The largest number of the additional nitro groups attached to the core molecule is two and their placement is different (Figure 2, Table 1). Hence, the molecules under the study are different due to different number of the nitro groups and their positions in the core molecules. We indicated the molecules investigated as HM-I - HM-VI to simplify discussion.

To observe the difference in the geometric structures of the molecules under study, the analysis of the bond length, angles and dihedral angles were performed. Generally, the bond length and angles of the core molecule

(HM-I) do not change significantly due to the additional nitro groups. However, the dihedral angles between planes that form (1, 2, 3), (5, 1, 2), (5, 1, 2) and (1, 2, 4) atoms (Figure 1) are this nitro group depended (Table 2).

Clearly the above dihedral angles are significantly different when an additional nitro group is in -R1 or/and -R2 position(s). The positions are rather close to the -NO₂ group of the core molecule (Figure 2). The difference between the geometrical structure of HM-I and HM-V, where the additional nitro group is in -R3 position, is insignificant (Figure 1). In this case, the additional nitro group is rather far from other nitro groups consisting of the molecules under investigations. It is obvious that the geometric structure of the molecules investigated is depended on the position of the additional group in the core molecule, i.e. depended on how far this group is from other nitro groups.

To compare the stability of the derivatives investigated, the binding energy per atom is calculated and presented in Table 4. It is evident, that the stability of the molecule HM-I, where R1, R2 and R3 are H atoms, is the lowest in comparison to that of other investigated compounds. However, the thermogravimetric analysis of HM -I was also performed to evaluate the thermal stability of this material. In the temperature range of 25˚C - 320˚C there was no observed mass change (Figure 3). That indicates that HM-I species are thermally stable, i.e. beyond ~320˚C temperature the material will begin to degrade.

It is evident that the stability of the other molecules should be higher in comparison to that of HM-I. However, we may not state that the stability of these molecules investigated increases when one of the above substitutions is -NO₂, because the binding energy per atom of molecule HM-II is only 0.013 eV, i.e. less than kT ~ 0.02 eV, and there is no possibility to recognize which of the above-mentioned molecules HM-I or HM-II is more stable (Table 3). Moreover, the comparison of the binding energy per atom of molecule HM-II and HM-VI (which are different because of the additional nitro group attached to different N atoms (Figure 2, Table 3), allows us to predict that even the placement of the additional -NO₂ group insignificantly influences the stability of N-(2,4,6-trinitrophenyl)-1H-1,2,4-triazol- 3-amine taking into account that in both cases of N-NO₂ similar additional chemical bonds are formed.

Molecules HM-V and HM-VI also have one additional -NO₂ group that is attached to the C and N atoms of the triazole ring, respectively. Referring to the obtained results of the binding energy per atoms, it is evident that the stability of molecule HM-V is greater than that of HM-IV. Recall that the stability of molecule HM-II and HM-VI, when only the N-NO₂ bond appeared, is not different. These observations allow us to conclude, that the stability of the molecules investigated could increase because of the presence of the additional -C-NO₂ bond. A similar phenomenon is obtained in the case of molecules HM-III and HM-IV. These molecules possess the two additional -NO₂ groups, and they are different due to the additional nitro group position in the triazole ring. In the case of HM-III, the N-NO₂ bond is formed, while C-NO₂ takes place in molecule HM-IV. Hence, the stability of N-(2,4,6-trinitrophenyl)-1H-1,2,4-triazol-3-amine could increase when the additional nitro group is attached to C atom. It is interesting that the number of the additional nitro groups does not have an influence on the stability of the investigated molecule: the binding energy per atom of molecule HM-IV and molecule HM-V is similar, although two and one additional nitro groups take place in these molecules respectively and one of them is attached to the C atom of the triazole ring.

To evaluate the chemical stability of the molecules under study, the HOMO, LUMO and the HOMO-LUMO gap are analyzed and presented in Table 4.

It is known that larger HOMO-LUMO gap indicates larger chemical stability of compounds. Hence, the sequence of the chemical stability of the molecules is as follows: HM-IV > HM-III > HM-II > HM-V > HM-VI > HM-I, i.e. HM-IV is chemically the most stable, while the chemical stability of the HM-I is the lowest. It must be emphasized, that the molecules consist of the two additional nitro groups are chemically more stable than other investigated molecules. On the other hand, the chemical stability of the HM-IV and HM-III also differs because one of the two additional nitro groups is joined with C, but not with N atom. Referring to the results, we may predict that an additional nitro group and its position could increase chemical stability of HM-I significantly.

To evaluate which of the molecules under study would be better to use as hazardous material, we calculated the activation energy―the minimum energy necessary to input to a system to cause a chemical reaction. The activation energy E_ac is evaluated as follows:

where E_gr is the total energy of the molecule at its equilibrium point, E_tr is that in the transition state. The activation energy obtained is presented in Table 5.

The results of the analysis of the activation energy exhibit, that the smallest energy must be input to the HM-IV to cause a chemical reaction. It is worth mentioning that the chemical stability of the molecule is the highest with respect to the other molecules under study. On the other hand, the activation energy of the molecules under study consisting of only one of the additional nitro groups is higher than that of HM-I. Hence, we may predict that the number of the additional nitro group have influences the amount of energy to be input to cause a chemical reaction.

There are two well-known and widely used indicators of the performance of the hazardous materials. These indicators are the detonation velocity and pressure. The detonation velocity, an indicator of strength, is calculated using Equation (1) and Equation (2). The results are presented in Table 6. The detonation velocity as an important property to be taken into consider when rating an explosive. The detonation velocities for high explosives range from 3300 fps to 29900 fps (1.01 km/s to 9.11 km/s). Hence, the results obtained indicate the molecules under investigations are highly explosive.

We would like to remind that the detonation velocity is calculated applying the regressed equations received on the basis of the quantum mechanical calculations [17] . Hence, it is not evident how precisely the above equations describe the molecules under investigation despite the fact that the Equation (1) and Equation (2) have been obtained on the basis of the research of the widely used compounds such as nitramines, aliphatic and aromatic-nitro. It implies, that the results obtained required an experimental confirmation. However, the detonation velocity of the molecules under investigation is obtained by using the same approach, i.e. the statistic errors of the investigations of the individual molecule could not be significantly different. It allows us to compare performance of the compounds investigated despite the model used and to rate the molecules under study. Referring to the results obtained the order of the strength of the materials is the following: HM-IV = HM-III > HM-V = HM-VI = HM-II > HM-I, i.e. the increasing the number of nitro groups in the molecule could strengthen the explosive performances.

We mention the detonation pressure as the other performance of the hazardous material to establish their ability to do the work and determine whether the explosive material is of high or low brisance.

Generally, the detonation velocity and detonation pressure are calculated by using Kamlet-Jacobs equation [18] :

where ρ is density of a compound (g/cm³), N is moles of gaseous detonation products, M is average molecular weight of gaseous products (g/mol), and Q is chemical energy of detonation (cal/g).

At present time we have no necessary data to evaluate detonation velocity by Kamlet-Jacobs Equation, but combining that, the following equations to calculate detonation pressure is

The detonation pressure obtained by Equation (7) is presented in Table 7.

It is necessary to mention the detonation pressure of TNT that is equal to 213 - 259 kbar and is used as a standard. So, the detonation pressure of the HM-I molecule is lower than that of TNT even then 10% - 20% deviation of the detonation velocity is accounted [19] . Referring to the results obtained, we could name the materials as high brisance (exception is HM-I) and rank them following: HM-I < TNT < HM-II < HM-V < HM-VI < HM-IV < HM-III. Moreover, the ability to do the work of the molecules consist of the two additional nitro groups is larger than that consisting only one or none of the above -NO₂ groups. Hence, the results obtained lead to the conclusion that the ability to do the work of the molecules under investigation depends on the number of the additional nitro groups.

Oxygen balance is used to indicate the degree to which an explosive can be oxidized. The sensitivity, strength, and brisance of hazardous materials are all somewhat dependent upon oxygen balance and tend to approach their maxima as oxygen balance approaches zero.

Referring to results presented in Table 8, we can conclude that the molecules under investigation contain less oxygen than needed and foresee that the combustion would be incomplete, and large amounts of toxic gases would be present. However, the negative oxygen balance could be used to indicate the shock sensitivity. It is expected that hazardous materials with the oxygen balance closer to zero will be more sensitive to a shock. Hence, the analysis of the oxygen balance indicates HM-I as less sensitive while HM-III or HM-IV as high sensitive. The results of the analysis indicate the relative sensitivity sequence of these molecules: HM-I > HM-III = HM-IV > HM-II = HM-V = HM-VI, i.e. the additional nitro groups could increase sensitivity of the materials.