Molecular design and theoretical study of oxadiazole-bifurazan derivatives

The design and synthesis of new high energy density materials is an important part of the research in the field of high energy materials. However, the synthesis of high-energy materials is very difficult and dangerous. Therefore, it is necessary to design the compounds in advance and evaluate the performance of the designed compounds, so as to screen the high-energy candidate compounds with excellent performance and provide reference for future synthesis and application. 1,2,5-oxadiazole (furazan) and 1,2,4-oxadiazole are five-membered nitrogen-oxygen heterocycles. Because their structures contain high-energy N-O, C=N bonds, they can effectively improve the energy density and oxygen balance of compounds, which has attracted widespread attention. In this paper, 42 kinds of oxadiazole-bifurazan energetic derivatives were designed by inserting different functional groups and changing the parent bridging groups with 1,2,4-oxadiazole and furazan as the basic structural units. Their electronic structures, aromaticity, heats of formation (HOFs), detonation properties, thermodynamic properties and electrostatic potential were systematically studied by density functional than theory (DFT). The results show that -C (NO₂)₃ has the greatest improvement effect on HOFs among all the substituent groups. The detonation performance of -N=N- bridged oxadiazole-bifurazan derivatives is better than that of -NH-NH- bridged derivatives. And -C(NO₂)₃ is the most effective group to improve the detonation performance and density of compounds. Compared with the parent compounds, when a -C(NO₂)₃ was introduced, the density increased by about 5.5%. A6 (D = 10.30 km·s^-1, P = 48.86 GPa) and D6 (D = 9.57 km·s^-1, P = 42.31 GPa) are the compounds with the best D and P among the designed compounds, which are higher than RDX and HMX, and are potential candidates for new high-energy materials.

With the help of Gaussian16 software and Multiwfn 3.8 package, the B3LYP method in density functional theory was selected. The 6-311G (d, p) basis set was used to optimize the structure of the 42 derivatives, and the high-precision def2-TZVPP basis set was used to calculate the energy.