Molecular dynamics and kinetic study of carbon coagulation in the release wave of detonation products

We present a combined molecular dynamics and kinetic study of a carbon cluster aggregation process in thermodynamic conditions relevant for the detonation products of oxygen deficient explosives. Molecular dynamics simulations with the LCBOPII potential under gigapascal pressure and high temperatures indicate that (i) the cluster motion in the detonation gas is compatible with Brownian diffusion and (ii) the coalescence probability is 100% for two clusters entering the interaction cutoff distance. We used these results for a subsequent kinetic study with the Smoluchowski model, with realistic models applied for the physical parameters such as viscosity and cluster size. We found that purely aggregational kinetics yield too fast clustering, with moderate influence of the model parameters. In agreement with previous studies, the introduction of surface reactivity through a simple kinetic model is necessary to approach the clustering time scales suggested by experiments (1000 atoms after 100 ns, 10 000 atoms after 1 μs). However, these models fail to reach all experimental criteria simultaneously and more complex modelling of the surface process seems desirable to go beyond these current limitations.