In this paper, the thermal decomposition physical models of different CL-20 polymorph crystals and their polymer bonded explosives (PBXs) bonded by polymeric matrices using polyisobutylene (PIB), acrylonitrile butadiene rubber (NBR), styrene butadiene rubber (SBR), Viton A, and Fluorel binders are obtained and used to predict the temperature profiles of constant rate decomposition. The physical models are further supported by the detailed decomposition pathways simulated by a reactive molecular dynamics (ReaxFF-lg) code. It has been shown that both ε-CL-20 and α-CL-20 decompose in the form of γ-CL-20, resulting in close activation energy (169 kJ mol −1 ) and physical model (first-order autoaccelerated model, AC1). Fluoropolymers could change the decomposition mechanism of ε-CL-20 from the "first-order autocatalytic" model to a "three-dimensional nucleation and growth" model (A3), while the polymer matrices of Formex P1, Semtex, and C4 could change ε-CL-20 decomposition from a single-step process to a multistep one with different activation energies and physical models. Compared to fluoropolymers, PIB, SBR and NBR may make ε-CL-20 undergo more complete N−NO 2 scission before collapse of the cage structure. This is likely the main reason why those polymer bases could greatly mitigate the decomposition process of ε-CL-20 from a single step to a multistep, resulting in lower impact sensitivity, whereas fluoropolymers have only a little effect on that. For ε-CL-20 and its PBXs, the impact sensitivity depends not only on the heat built-up period of their decomposition, but also on the probability of hotspot generation (defects in solid crystals and interfaces) especially when it decomposes in a solid state.