JBCS



3:37, sex dez 6

Acesso Aberto/TP




Artigo

Theoretical studies on triacetonetriperoxide (tatp) derivatives to improve their performance

Zongcheng MiaoI; Faan LiII; Yi LuanIII,*

I. Key Laboratory of Organic Polymer Photoelectric Materials, School of Science, Xijing University, 710123 Xi'an, PR China
II. Department of Business Administration, School of Management and Economics, Beijing Institute of Technology, 100081 Beijing, PR China
III.Department of Materials Physics, School of Materials Science and Engineering, University of Science and Technology Beijing, 100083 Beijing, PR China

Recebido em: 14/05/2018
Aceito em: 05/09/2018
Publicado em: 27/09/2018

Endereço para correspondência

*e-mail: yiluan@ustb.edu.cn

RESUMO

Density functional theory (DFT) methods are used to study on TATP derivatives at DFT-B3LYP/6-311++G (d, p) level to improve TATP's performance. The derivatives are 3,6,9-trimethyl-3,6,9-tris(trifluoromethyl)-1,2,4,5,7,8-hexaoxonane (TATP3F), 3,6,9-trimethyl-3,6,9-tris(nitromethyl)-1,2,4,5,7,8-hexaoxonane (TATPNO2) and (3,6,9-trimethyl-1,2,4,5,7,8-hexaoxonane-3,6,9-triyl)trimethanamine (TATPNH2). Stability and detonation performance of TATP derivatives were investigated. Detonation performance was estimated using density, heat of formation and EXPLO 5 software. Stability was evaluated by geometrical geometry and electronic structure. The simulation results reveal that TATPNO2 has a best oxygen balance, detonation performance (Q=5424 kJ kg-1, P=23.7 Gpa, D=7699 m s-1) and insensitivity. Since TATPNO2 is more excellent than TATP and traditional explosive TNT in detonation performance, it is a candidate of novel high energy density materials.

Palavras-chave: TATP; density functional theory (DFT); detonation performance; sensitivities.

INTRODUCTION

High energy density materials (HEDMs) have attracted significant attention owing to their massive usage as explosive, propellants, pyrotechnic agents, gas generators and so on in a wide range of applications.1-4

A potentially initiating explosive peroxide triacetonetriperoxide (TATP), which is shown in Figure 1, is not viewed as HEDMs due to its poor performance. However, TATP has been extensively used by terrorists in terrorist attack because of its low-cost and simple production process.5,6 In the meantime, the obvious advantage of TATP is that it has an affluent oxygen content. There are two ways to improve the oxygen balance of HEDMs. One is to introduce nitro and difluoramino,7,8 and the other is to use oxygen-rich skeletons.9 TATP, with six oxygen atoms in a ring, is expected to be a promising energetic skeleton. The disadvantage of TATP is its poor sensitivity and detonation performance10 and modifications are needed to improve TATP's performance for meeting the requirements of HEDMs in military and civilian use.

 


Figure 1. Chemical structure of TATP

 

In this paper, we proposed to add three commonly energetic and stable substituents, -CF3, -NO2 and -NH2, separately to the rich oxygen skeleton of TATP in order to lower its sensitivity and enhance its detonation performance. The three new compounds are named 3,6,9-trimethyl-3,6,9-tris(trifluoromethyl)-1,2,4,5,7,8-hexaoxonane (TATP3F), 3,6,9-trimethyl-3,6,9-tris(nitromethyl)-1,2,4,5,7,8-hexaoxonane (TATPNO2) and (3,6,9-trimethyl-1,2,4,5,7,8-hexaoxonane-3,6,9-triyl)trimethanamine (TATPNH2), respectively (Figure 2).

 


Figure 2. Chemical structure of TATP derivatives

 

The properties of these three new compounds were studied by theoretical calculation. Theoretical calculations based on quantum chemistry have been accepted as a time-saving tool to select the expected candidates for effectively minimizing expensive and dangerous experimental tests.11,12 This paper presents crystal density (ρcry) and heat of formation (ΔHf) for calculating detonation performance of title compounds, such as explosion heat (Q), detonation velocity (D) and detonation pressure (P). Besides, this paper also presents electrostatic potential analysis (ESP) and energy gap (ΔELUMO-HOMO) of the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) for calculating sensitivity of title compounds. Their characteristics were used to compare with TATP's one to screen new energetic candidates with excellent performance and acceptable sensitivity. These results will provide theoretical support for molecular design of novel HEDMs.

 

COMPUTATIONAL METHODS

The calculations were performed with the Gaussian 09 package.13 DFT-B3LYP methods with the 6-311G++ (d, p) basis set14-16 were adopted to optimize the molecular structures and predict the heats of formation. Previous studies have shown that the 6-311G++ (d, p) basis set is able to precisely predict the molecular structures and energies of energetic organic compounds.17,18 All structures were identified to be local minima without imaginary frequencies.

Bond lengths, especially three O-O bond lengths, of title compounds were one of the major focuses of the calculation study. These bond lengths are used to compare to the ones of TATP. Generally speaking, the bond length is closely related to the bond stability: the longer the bond is, the less stable the bond is.19

The crystals density is an important parameter of energetic materials to predict detonation performance. It is predicted by using Equation (1) with considering the intermolecular interactions in crystals. In the Equation, M is the molecular weight, and Vm is the molecular volume defined as inside a contour of 0.001 au density that was evaluated using a Monte Carlo integration. Besides, σ2tot and v is total variance on the molecular surface and electrostatic balance parameter, respectively. The two parameters were obtained by Multiwfn program.20 Finally, in order to make a better account for intermolecular interactions and enhance the accuracy of predictions, Politzer et al. 21 proposed to use correction parameters α, β and γ. They compared experimental and calculated densities of 36 compounds and revealed that α, β and γ is 0.9183, 0.0028, 0.0443, respectively.

The heat of formation is another important parameter to predict detonation performance. It was predicted by the Equation (2)-(4) with isodesmic reactions (Schemes 1-3). Equation (2) was used to calculate the gas phase heats of formation of the title compounds where ΔfHP and ΔfHR are the heat of formation of the products and the reactants in isodesmic reactions, respectively. ΔE0 is the change in total energy between the products and the reactants at 0 K. ΔEZPE is the difference between the zero-point energy (EZPE) of the products and the reactants, and ΔHT is the thermal correction from 0 K to 298 K. ΔnRT is equal to zero here.22

 


Scheme 1. Isodesmic reactions of TATP3F

 

 


Scheme 2. Isodesmic reactions of TATPNO2

 

 


Scheme 3. Isodesmic reactions of TATPNH2

 

Equation (3) and (4) was used to calculate the condense phase heat of formation of the title compounds where ΔfH(c), ΔfH(g) and ΔHsub is the condense phase heat of formation, the gas phase heat of formation and enthalpy of sublimation, respectively. The meaning and value of n and σ2tot here is the same to Equation (1), and as is molecular surface area. In Equation (4), a, b, c developed by Politzer23 is 0.000267, 1.650087, 2.966078, respectively.

Based on the condense phase heat of formation and crystals density, explosion heat (Q), detonation velocity (D) and detonation pressure (P) of the title compounds can be calculated by EXPLO 5(v6.01), a program widely applied as an important tool to compute detonation performance.24-26

The energy gap ΔELUMO-HOMO between HOMO and LUMO is an important parameter of energetic materials to predict kinetic stability. It is worth noting that stability refers to the stability of the compound to photoinitiated electron-transfer processes.27 The energy gap shows that the intramolecular charges transfer from the electron donor to the electron acceptor. The higher the magnitudes are, the less stable the title compounds are.

The electrostatic potential (ESP) is another important parameter of energetic materials to predict stability, especially nucleophilic and electrophilic sites, as well as the molecular recognition mode for a long time. Their theoretical basis is that molecules always tends to approach each other.20,28,29 The ESP can be shown in a map in which the blue and red spheres correspond to electrophilic and nucleophilic ESP surfaces, respectively.

 

RESULTS AND DISCUSSION

Geometrical structures

The optimized geometries of the most stable conformers were calculated with B3LYP/6-311G++ (d, p) method. The optimized structure of TATP and its derivatives compounds were shown in Figure 3. All harmonic frequencies and summaries of their output parameters were shown in Supplementary Material. It is reported30 that due to the weak O-O bond, peroxides undergo facile thermal decomposition. Since the bond length is closely related to the bond stability, the O-O bond length is very important to title compounds.

 


Figure 3. The optimized structure of TATP and its derivatives

 

Three O-O bond lengths of these compounds were used to compare to the counterpart of TATP, shown in Table 1. Generally speaking, the bond length is closely related to the bond stability and the longer the bond is, the less stable the bond is. From the Figure 3 and Table 1, it can be seen that O-O bond lengths of TATP3F and TATPNO2 are shorter while the one of TATPNH2 is longer, since nitro and trifluoromethyl are electron-withdrawing groups which can decrease electron cloud densities of O-O bonds. Considering structure, TATP3F and TATPNO2 have an excellent O-O bond length and a stable structure. More information on standard orientation of every atom is presented in Supplementary Material in detail.

 

 

Detonation performance

In general, a high density causes better performance. For example, the detonation velocity of an explosive is proportional to its density, and the detonation pressure is proportional to its square.1,11 The crystal densities of title compounds were presented in Table 2. It is reported31 that experimental density of TATP is 1.22 g cm-3, consistent with calculated value 1.28 g cm-3. This results suggested that computational methods are effective. In Table 2, their densities order is TATP3F > TATPNO2 > TATPNH2 > TATP. Due to N, O and F atoms, density of TATP derivatives are higher than TATP. Fluorine- and oxygen-rich groups have long been recognized as useful building blocks for the design and synthesis of superior energetic compounds and these building blocks provide energy density from molecules containing them.8,32 Calculation results showed that TATP3F's density is higher than HMX's (1.90 g cm-3).33

 

 

The value of the heat of formation must be obtained to estimate the performances of an explosive molecule through theoretical calculations. Their condensed phase heats of formation were shown in Table 3. It is reported34 that heat of formation of TATP is -664.5 kJ mol-1, lower than the one of TATPNH2 and TATPNO2 with substituted energetic groups -NH2 and -NO2, respectively. The incorporation of three -NO2 groups into a molecule improves the heat of formation, which is consistent with a paper reported by Fei.35 It has also been found that compounds containing the halogen fluorine lower heats of formation.35,36 Calculation show that the value of TATP3F is least among them because of F atom replacing of H atom.

 

 

The detonation performances, such as detonation heat (Q), detonation velocity (D) and detonation pressure (P), of title compounds were summarized in Table 4. Detonation performance of TATP is also obtained by the same method to compare. It reveals that detonation pressure and detonation velocity of title compounds with substituent groups range from 14.0 to 18.1 Gpa and 6458 to 7699 m s-1 respectively, higher than that of TATP. Among the four compounds, TATPNO2 has relatively high values in both density and heat of formation, which contributes to its best performance. Detonation performance of TATPNO2 (Q=5424 kJ kg-1, P=23.7 Gpa, D=7699 m s-1) is higher than that of TNT (Q=5418 kJ kg-1, P=19.0 Gpa, D=6950 m s-1).38

 

 

Electronic structure

Electronic structure of title compounds can reflect photochemistry stability. In this paper, electronic structure includes electrostatic potential and energy gap ΔE between HOMO and LUMO. Electrostatic potential of TATP and title compounds are shown in Figure 4, in which the legend for the range of the ESPs is -0.03 (left) to +0.03 (right). While the HOMO and LUMO orbitals were depicted in Figure 5, in which red color represents positive phase and the blue color represents negative phase. As can be seen from the figure, the LUMO orbitals mainly locate on substituted groups but the HOMO orbitals distribute approximately on the rings. The intramolecular charges transfer from the electron donor to the electron acceptor through the rings and substituted groups.

 


Figure 4. Electrostatic potential of TATP and title compounds

 

 


Figure 5. Highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) of the four compounds

 

The values of energy gap ΔE of TATP and its derivatives were shown in Table 5. The electron density of TATPNH2 seems to be distributed more evenly over the surface of the molecule and that would account for a lower impact sensitivity.29 However, the sequence of energy gap is TATP3F > TATP > TATPNH2 > TATPNO2, which is the reverse of chemical activity in order. In conclusion, TATPNO2 has a good photochemistry stability according to electronic structure, and this result is consistent with geometrical structures in which O-O bond length is shorter.

 

 

CONCLUSION

In recent years TATP has been used as an improvised explosive because its precursor chemicals are readily obtained and its synthesis is straightforward. However, it is not the high energy density material because of its lower detonation performance and higher sensitivity. This paper is expected to find some alternative materials and TATP derivatives were studied for this reason. They are TATP3F, TATPNH2 and TATPNO2, which were introduced with trifluoromethyl, amino and nitro groups to the TATP backbone, respectively.

The TATP derivatives were studied to improve TATP's performance on detonation performance or stability. TATP3F, TATPNH2 and TATPNO2 were calculated by DFT methods based on 6-311G++ (d, p). Their geometrical structures, detonation performance and stability were investigated in detail.

Calculation suggests that in terms of both the detonation performance and stability, TATPNO2 is better than others, especially TATP. As far as stability is concerned, three nitros shorten O-O bond length and decrease energy gap ΔE between HOMO and LUMO. At the same time, TATPNO2 has a relatively high density and heat of formation which are a contribution to power. Three nitros enclosing to rich-oxygen skeleton enhance oxygen content which helps to burn. Besides, detonation heat, detonation velocity and detonation pressure of TATPNO2 (Q=5424 kJ kg-1, P=23.7 Gpa, D=7699 m s-1) is higher than that of traditional explosive TNT (Q=5418 kJ kg-1, P=19.0 Gpa, D=6950 m s-1). Therefore, TATPNO2 has a potential application as high energy density material because of its power, good stability, availability and rich-oxygen.

 

SUPPLEMENTARY MATERIAL

Standard orientations for atoms, harmonic frequencies of TATP and derivatives and calculation summaries of the Gaussian output parameters can be freely accessed at http://quimicanova.sbq.org.br, in pdf format.

 

ACKNOWLEDGEMENTS

This work is supported by Beijing Natural Science Foundation (2172037) and National Natural Science Foundation of China (No. 51673157 and 51503016). Z. C. M. also thanks the Science Research Foundation of Xijing University (Grant No. XJ16T02).

 

REFERENCES

1. He, P.; Zhang, J. G.; Wu, L.; Wu, J. T.; Zhang, T. L.; J. Phys. Org. Chem. 2017, 30, 1.

2. Zhang, W.; Wang, K.; Li, J.; Lin, Z.; Song, S.; Huang, S.; Liu, Y.; Nie, F. ; Zhang, Q.; Angew. Chem. 2018, 57, 1.

3. Chen, J.; Yu, Y.; Li, Y. ; Pang, S.; J. Fluorine Chem. 2018, 205, 35.

4. Zhao, G.; Lu, M.; Quim. Nova 2013, 36, 513.

5. Oxley, J. C.; Smith, J. L.; Bowden, P. R.; Rettinger, R. C.; Propellants, Explos., Pyrotech. 2013, 38, 244.

6. Oxley, J. C.; Smith, J. L.; Chen, H.; Propellants, Explos., Pyrotech. 2002, 27, 209.

7. Dalinger, I. L.; Kormanov, A. V.; Suponitsky, K. Y.; Muravyev, N. V.; Sheremetev, A. B.; Chem. - Asian J. 2018, 13, 1165.

8. Semenov, V. V.; Shevelev, S. A.; Bruskin, A. B.; Shakhnes, A. K.; Kuz'min, V. S.; Chem. Heterocycl. Compd. 2017, 53, 728.

9. Tsyshevsky, R.; Pagoria, P.; Zhang, M.; Racoveanu, A.; Parrish, D. A.; Smirnov, A. S.; Kuklja, M. M.; J. Phys. Chem. C 2017, 121, 23853.

10. Contini, A. E.; Bellamy, A. J.; Ahad, L. N.; Propellants, Explos., Pyrotech. 2012, 37, 320.

11. Zhao, G. Z.; Lu, M.; Bull. Korean Chem. Soc. 2012, 33, 1913.

12. He, P.; Zhang, J. G.;Wang, K.; Yin, X.; Zhang, T. L.; J. Org. Chem. 2015, 80, 5643.

13. Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Montgomery Jr., J. A.; Vreven, T.; Kudin, K. N.; Burant, J. C.; Millam, J. M.; Iyengar, S. S.; Tomasi, J.; Barone, V.; Mennucci, B.; Cossi, M.; Scalmani, G.; Rega, N.; Petersson, G. A.; Nakatsuji, H.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Klene, M.; Li, X.; Knox, J. E.; Hratchian, H. P.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Ayala, P. Y.; Morokuma, K.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Zakrzewski, V. G.; Dapprich, S.; Daniels, A. D.; Strain, M. C.; Farkas, O.; Malick, D. K.; Rabuck, A. D.; Raghavachari, K.; Foresman, J. B.; Ortiz, J. V.; Cui, Q.; Baboul, A. G.; Clifford, S.; Cioslowski, J.; Stefanov, B. B.; Liu, G.; Iiskorz, P.; Komaromi, I.; Martin, R. L.; Fox, D. J.; Keith, T.; Al-Laham, M. A.; Peng, C. Y.; Nanayakkara, A.; Challacombe, M.; Gill, P. M. W.; Johnson, B.; Chen, W.; Wong, M. W.; Gonzalez, C.; Pople, J. A.; Gaussian 09(Rev. A.01), Gaussian Inc., Wallingford CT., 2009.

14. Calais, J.-L.; Int. J. Quantum Chem. 1993, 47, 101.

15. Abashkin, Y.; Russo, N.; Sicilia, E.; Toscano, M.; Theor. Comput. Chem. 1995, 2, 255.

16. Krishnan, R.; Binkley, J. S.; Seeger, R.; Pople, J. A.; J. Chem. Phys. 1980, 72, 650.

17. Gao, H. X.; Ye, C. F.; Piekarski, C. M.; Shreeve, J. M.; J. Phys. Chem. C 2007, 111, 10718.

18. Fan, X. W.; Ju, X. H.; J. Comput. Chem. 2008, 29, 505.

19. Jin, X. H.; Zhou, J. H.; Wang, S. J.; Hu, B. C.; Quim. Nova 2016, 4, 467.

20. Lu, T.; Chen, F.; J. Comput. Chem. 2012, 33, 580.

21. Politzer, P.; Martinez, J.; Murray, J. S.; Concha, M. C.; Toro-Labbe, A.; Mol. Phys. 2009, 107, 2095.

22. Zhang, X.; Gong, X.; J. Chem. Eng. Data 2015, 60, 2869.

23. Byrd, E. F.; Rice, B. M.; J. Phys. Chem. A 2006, 110, 1005.

24. Suceska, M.; EXPLO5, version 6.01; Brodarski Institute, Zagreb, Croatia, 2013.

25. Suceska, M.; Propellants, Explos., Pyrotech. 1999, 24, 280.

26. Yu, Y.; Chen, J.; Zhang, R.; Li, Y.; Pang, S.; RSC Adv. 2017, 7, 23709.

27. Chi, W. J.; Li, L. L.; Li, B. T.; Wu, H. S.; J. Mol. Model. 2012, 18, 3695.

28. Murray, J. S.; Politzer, P.; Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2011, 1, 153.

29. Hammerl, A.; Klapotke, T. M.; Noth, H.; Warchhold, M.; Holl, G.; Propellants, Explos., Pyrotech. 2003, 28, 165.

30. Oxley, J. C.; Smith, J. L.; Chen, H.; Propellants, Explos., Pyrotech. 2002, 27, 209.

31. Williams, M.; Drug Dev. Res. 2013, 74, 339.

32. Lempert, D. B.; Nechiporenko, G. N.; Soglasnova, S. I.; Dorofeenko, E. M.; J. Phys. Chem. B 2017, 11, 443.

33. Talawar, M. B.; Sivabalan, R.; Mukundan, T.; Muthurajan, H.; Sikder, A. K.; Gandhe, B. R.; Rao, A. S.; J. Hazard. Mater. 2009, 161, 589.

34. Sinditskii, V. P.; Kolesov, V. I.; Egorshev, V. Y.; Patrikeev, D. I.; Dorofeeva, O. V.; Thermochim. Acta 2014, 585, 10.

35. Fei, T.; Du, Y.; Pang, S.; RSC Adv. 2018, 8, 10215.

36. Kettner, M. A.; Karaghiosoff, K.; Klapotke, T. M.; Suceska, M.; Wunder, S.; Chem. - Eur. J. 2014, 20, 7622.

37. https://webbook.nist.gov/chemistry/, accessed in September 2018.

38. Politzer, P.; Murray, J. S.; Cent. Eur. J. Energ. Mater. 2011, 8, 209.

On-line version ISSN 1678-7064 Printed version ISSN 0100-4042
Qu�mica Nova
Publica��es da Sociedade Brasileira de Qu�mica
Caixa Postal: 26037 05513-970 S�o Paulo - SP
Tel/Fax: +55.11.3032.2299/+55.11.3814.3602
Free access

GN1