where ri is the rate of reaction type i (i.e., each possible reaction shown in Scheme 1). A reaction is selected and executed stochastically based on these probabilities and the system is updated at each step; a monomer reacts and is removed from the pool in a propagation event or depropagated from a chain end and added back to the monomer pool. Subsequently, the composition of a polymer chain in the simulation is modified. This process is repeated until the monomer pool is exhausted or until the polymerization reaches equilibrium (i.e., the rates of propagation and depropagation are equal). The stochastic simulation method used in this work is similar to those reported previously following a general kinetic Monte Carlo (KMC) scheme but has been modified to include the reaction steps consistent with the general model shown in Scheme 1 (see the Stochastic Simulations section of the Supporting Information for details). For this method, we found that typically 1,000 chains were sufficient to model copolymerization (Figure ); we employed 5,000 to 10,000 for the comparisons shown in this work.