Abstract. The impact of stochastic fluctuations in cloud droplet growth is a matter of broad interest, since stochastic effects are one of the possible explanations of how cloud droplets cross the size gap and form the raindrop embryos that trigger warm rain development in cumulus clouds. Most theoretical studies on this topic rely on the use of the kinetic collection equation, or the Gillespie stochastic simulation algorithm. However, the kinetic collection equation is a deterministic equation with no stochastic fluctuations. Moreover, the traditional calculations using the kinetic collection equation are not valid when the system undergoes a transition from a continuous distribution to a distribution plus a runaway raindrop embryo (known as the sol-gel transition). On the other hand, the stochastic simulation algorithm, although intrinsically stochastic, fails to adequately reproduce the large end of the droplet size distribution due to the huge number of realizations required. Therefore, the full stochastic description of cloud droplet growth must be obtained from the solution of the master equation for stochastic coalescence.In this study the master equation is used to calculate the evolution of the droplet size distribution after the sol-gel transition. These calculations show that after the formation of the raindrop embryo, the expected droplet mass distribution strongly differs from the results obtained with the kinetic collection equation. Furthermore, the low-mass bins and bins from the gel fraction are strongly anticorrelated in the vicinity of the critical time, this being one of the possible explanations for the differences between the kinetic and stochastic approaches after the sol-gel transition. Calculations performed within the stochastic framework provide insight into the inability of explicit microphysics cloud models to explain the droplet spectral broadening observed in small, warm clouds.