This study focuses on a statistical description of droplet sizes, created as a result of capillary-induced breakup of ligaments. Direct numerical simulations of air-water systems are employed by solving the two-phase Navier-Stokes equations on adaptive Octree grids (http://basilisk.fr/), using the VOF methodology coupled with height-function-based curvature modeling. Breakup of individual ligaments are triggered by initial surface corrugations, the dynamics of which are deterministic. Stochasticity is introduced in the mix by conducting an ensemble of simulations of slender corrugated ligaments, each realization corresponding to a random but unique initial configuration. Probability density functions of the droplet sizes are computed for different ensemble sizes. These results (fig. 1) combining the effects of stochasticity with the capillarity-driven nonlinear dynamics are compared to the predictions of a Gaussian random process theory for near-monochromatic waves, which facilitate our understanding of the nature of drop size distributions encountered in realistic and convoluted fluid fragmentation scenarios.