The accuracy of prediction against DFT is 1.15 and 1.51 kcal mol -1 for BDE and BDFE respectively. Previous data-driven models have focused on the prediction of C-F BDE values in PFAS molecules; 49 however, with the ability to predict across the breadth of bond types present in these compounds, we observe that C-S and C-C bonds (rather than C-F) are thermodynamically much more likely candidates for the primary site of homolysis. For 60% of the PFAS considered, the ML-predicted weakest bond matches DFT, while for the remaining 40% of cases, the weakest bond (from DFT) lies within 4 kcal mol -1 from the ML-minimum energy. This suggests one possible use of ML could be to quickly survey and rank possible homolytic cleavages, returning a focused set of candidate bonds to be investigated in greater depth with QM calculations.