Dynamic Monte Carlo simulation of atom-transfer radical polymerization

A dynamic Monte Carlo model was developed to simulate atom-transfer radical polymerization (ATRP). The algorithm used to describe the polymerization includes activation, deactivation, propagation, chain transfer, and termination by combination and disproportionation reactions. Model probabilities are calculated from polymerization kinetic parameters and reactor conditions. The model was used to predict monomer conversion, average molecular weight, poly-dispersity and the complete molecular weight distribution at any polymerization time or monomer conversion. The model was validated with experimental results for styrene polymerization and compared with simulation results from a mathematical model that uses population balances and the method of moments. The simulations agree well with experimental and theoretical results reported in the literature. We also investigated the control volume size and number of iterations to reduce computation time while keeping an acceptable noise level in the Monte Carlo results.