agreement is also quantitative, for suitably chosen λ , for low chucking rates and when the solute particles are small compared to the colloid. This work is structured as follows. In the next Section, we describe in detail our simulation model and methods. In Section III we report results for a diffusing chucker, a dragged chucker and a chucker close to a surface. Finally, Section IV contains a discussion and our conclusions. We model the chucker as a hard sphere of radius R$_{1}$ , with long-time diffusion coefficient D$_{1}$ . With a rate k$_{c}$ , the sphere produces (“chucks”) smaller particles (solute) isotropically at its surface. The solute particles have radius R$_{2}$ and long-time diffusion coefficient D$_{2}$ , which is related to D$_{1}$ by Stokes’ law: D$_{2}$/D$_{1}$ = R$_{1}$/R$_{2}$ . In our simulations, the length unit is the chucker radius R$_{1}$ and we use two different values for the solute radius R$_{2}$ : R$_{2}$/R$_{1}$ = 0 . 1 ( R$_{1}$/R$_{2}$ = 10 ) or R$_{2}$/R$_{1}$ = 0 . 01 ( R$_{1}$/R$_{2}$ = 100 ). Our main simulation results are obtained using a dynamical Monte Carlo algorithm. In this algorithm, a particle is chosen at random and subjected to a trial displacement δ , in a randomly chosen direction and with magnitude uniformly distributed in the range 0 → 0 . 001 R$_{1}$ (for the chucker) and 0 → 0 . 001 R$_{1}$ √ R$_{1}$/R$_{2}$ for the solute particles (this ratio of displacements is to ensure the correct ratio between the diffusion coefficients D$_{1}$ and D$_{2}$ , as given by Stokes’ Law). If this trial displacement results in an overlap between chucker and solute, it is rejected; otherwise it is accepted. One MC cycle consists of N trial displacements, where N is the current number of particles in the system. It has been shown in previous literature [22–25, 28–30] that this algorithm is equivalent to Brownian Dynamics for small trial displacements. An MC cycle can be interpreted as the time unit in a dynamical trajectory, where the diffusion constant in reduced length units is given by D$_{1}$/R 2 1 = ( δR 1 /R 1 ) 2 6 [MC cycles] − $^{1}$. At the beginning of every MC cycle, a new solute particle is inserted with uniform probability of k$_{c}$ , randomly at the surface of the chucker. Note that because the solutes do not interact, all particle insertions are accepted. All values for II. MODEL AND METHODS A. Monte Carlo simulations 5