Numerical convergence of the collision–coalescence algorithm used in Lagrangian particle-based microphysics is studied in 2D simulations of an isolated cumulus congestus (CC) and in box and multi-box simulations of collision–coalescence. Parameters studied are the time step for coalescence and the number of super-droplets (SDs) per cell. A time step of 0.1 s gives converged droplet size distribution (DSD) in box simulations and converged mean precipitation in CC. Variances of the DSD and of precipitation are not sensitive to the time step. In box simulations, mean DSD converges for 10³ SDs per cell, but variance of the DSD does not converge as it decreases with an increasing number of SDs. Fewer SDs per cell are required for convergence of the mean DSD in multi-box simulations, probably thanks to mixing of SDs between cells. In CC simulations, more SDs are needed for convergence than in box or multi-box simulations. Mean precipitation converges for 5×10³ SDs, but only in a strongly precipitating cloud. In cases with little precipitation, mean precipitation does not converge even for 10⁵ SDs per cell. Variance in precipitation between independent CC runs is more sensitive to the resolved flow field than to the stochasticity in collision–coalescence of SDs, even when using as few as 50 SDs per cell.