The Gillespie-based kinetic Monte Carlo (kMC) method is employed to simulate terpolymerization of which the main algorithm is included our recent review. In addition, molecular information on monomer sequences is stored via coupled matrices, making the overall method coupled matrix-based Monte Carlo (CMMC) of which the main algorithm embedding the Gillespie-based algorithm has been published in full detail in Figueira et al. and De Keer et al. The reader is referred to the Supporting Information for more details on the CMMC model.