Comparison study of carbon clusters formation during thermal decomposition of 1,3,5-triamino-2,4,6-trinitrobenzene and benzotrifuroxan: a ReaxFF based sequential molecular dynamics simulation

Carbon rich clusters are usually found after the detonation of explosives, which greatly hinder their further decomposition into small molecules. A comparison study of thermal decomposition and clusters formation between 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) and benzotrifuroxan (BTF) crystals was conducted to uncover the mechanisms behind their distinct differences in sensitivity and reaction violence, which has not been investigated in detail. The simulations of heating at 3500 K, then expansion and cooling were conducted through reactive molecular dynamics using the ReaxFF-lg force field. As a result, the initial low decay rate indicates that TATB is more stable than BTF under high temperatures, while once ignited it decays faster than BTF. Nevertheless, BTF decomposes more completely with a higher potential energy release, a greater amount of final products, and higher reaction frequencies, and shows higher reaction violence than TATB. More and heavier clusters occur in TATB crystals compared with those in BTF. Large clusters form during the heating process and then partly dissociate during expansion and cooling. A faster cooling rate facilitates larger clusters formation. Graphitic geometries as well as carbon rings and carbon chains are common in the stable clusters. Besides, further simulations show that a lower heating temperature facilitates larger clusters formation both in TATB and BTF. Our results are expected to deepen the insight into the mechanisms of carbon clusters formation and the different performances of TATB and BTF.