A Monte Carlo simulation model for the kinetics of emulsion polymerization is proposed. In the present model, the formation of each polymer molecule is simulated by the use of only a couple of probability functions; therefore, the calculation can be handled well even on personal computers. It is straightforward to account for virtually any kinetic event, such as the desorption of oligomeric radicals and chain length dependence of kinetic parameters, and as a consequence very detailed information such as the full distributions of the dead polymer molecular weights and the macroradicals among various polymer particles can be obtained. When bimolecular terminations are the dominant chain stoppage mechanism, the instantaneous molecular weight distribution (produced in a very small time interval) becomes broader than that for homogeneous polymerizations due to a higher possibility that short and long polymer radicals react with each other if bimolecular reactions are fast enough. The increase in the polydispersity of the MWD is fairly large, especially when bimolecular termination by disproportionation is significant; however, the gel permeation chromatography (GPC) may not be a suitable analytical technique to detect such broadening since oligomeric peaks may not be observed in the elution curve. The present simulation method provides greater insight into the complicated phenomena of emulsion polymerizations.