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 N2. We also observe the formation of large nitrogen- and oxygen-rich carbon aggregates, which delays the release of final reaction products.