Due to the prominent role of the propargyl radical for hydrocarbon growth within combustion environments, it is important to understand the kinetics of its formation and loss. The ab initio transition state theory-based master equation method is used to obtain theoretical kinetic predictions for the temperature and pressure dependence of the thermal decomposition of propargyl, which may be its primary loss channel under some conditions. The potential energy surface for the decomposition of propargyl is first mapped at a high level of theory with a combination of coupled cluster and multireference perturbation calculations. Variational transition state theory is then used to predict the microcanonical rate coefficients, which are subsequently implemented within the multiple-well multiple-channel master equation. A variety of energy transfer parameters are considered, and the sensitivity of the thermal rate predictions to these parameters is explored. The predictions for the thermal decomposition rate coefficient are found to be in good agreement with the limited experimental data. Modified Arrhenius representations of the rate constants are reported for utility in combustion modeling.