Detonation processes probed with atomistic details have remained elusive due to highly complex reactions in heterogeneous shock structures. Here, we provide atomistic details of the initial reaction pathways during shock-induced decomposition of 2,4,6-triamino-1,3,5-trinitrobenzene (TATB) crystal using large reactive molecular dynamics simulations based on reactive force fields. Simulation results reveal the existence of three competing intermolecular pathways for the formation of N 2 . We also observe the formation of large nitrogen-and oxygen-rich carbon aggregates, which delays the release of final reaction products.