The pyrolysis kinetics of hexamethyldisiloxane (HMDSO) at various temperatures was studied using the reactive force field (ReaxFF) molecular dynamics simulations. Reaction rate constants and the main pyrolysis pathways were explored at the initial decomposition stage of HMDSO and intermediates decomposition stage. The activation energy and pre-exponential factor describing the reaction rate constants were obtained and further validated by experimental data and DFT theoretical calculations. The formation of C 5 H 15 OSi 2 fragment by Si-C bond dissociation was dominant at the initial decomposition stage of HMDSO at the simulation temperatures of 2500-4000 K. The subsequent reaction pathways involved the formation of C 5 H 14 OSi 2 and C 4 H 11 OSi 2 . After that, the pathways were different for 2500 and 3000 K. At 4000 K, small silicon-containing fragments were formed, including CH 3 Si, CH 4 Si, and C 3 H 9 Si, etc. The simulations also revealed that the major hydrocarbons generated during HMDSO pyrolysis were CH 3 and CH 4 . Also, CH 4 formation was more important in the end of HMDSO pyrolysis when simulation temperature was over 3500 K.