A Metropolis scheme that can be used by various Monte Carlo methods for solving particle population balance equations has been proposed and verified. A detailed derivation of the scheme is presented. The computational accuracy of the proposed scheme is tested under various coagulation kernels. The computational performance of the scheme is also measured and compared to that of Inverse and acceptance-rejection scheme. (C) 2015 Published by Elsevier Ltd.