The degradation of polystyrene was modeled at the mechanistic level using the method of moments to track structurally distinct polymer species. To keep the model size manageable, polymer species were lumped into groups, and within these groups, the necessary polymeric features for capturing the degradation chemistry were tracked. The pyrolysis reactions incorporated into the model included hydrogen abstraction, midchain β-scission, end-chain β-scission, 1,5-hydrogen transfer, 1,3-hydrogen transfer, radical addition, bond fission, radical recombination, and disproportionation. From the evolution of the zeroth, first, and second moments tracked for each dead species, polymer molecular weight distributions were constructed by summing the Schultz (