Iron oxide (Fe3O4) nanoparticle was synthesized by coprecipitation and was modeled and solved using a hybrid (discrete-continuous) model, based on a kinetic Monte Carlo (kMC) simulation scheme. The latter was combined with the constant number MC method, to improve both speed and accuracy of the simulation. Complete particle size distribution (PSD) from simulation matches very well with PSD of both uncoated and coated (with either polyacrylic acid or dextran) Fe3O4 nanoparticles, obtained from our experiments. The model is general, as the time scales of various processes (nucleation, diffusion-growth and coagulation-growth) are incorporated in rate equations, while, input simulation parameters are experimentally measured quantities. With the help of the validated model, effect of coating agent on coagulation-growth was estimated by a single, fitted, coagulation-efficiency parameter. Our simulation shows that, logarithm of coagulation-efficiency scales linearly with logarithm of inverse of the molecular weight of the coating agent. With this scaling law, our model is able to a priori predict the experimental PSD of Fe3O4 nanoparticles, synthesized with an even higher molecular weight of dextran. © 2015 Elsevier Inc..