The main advantage for applications of graphene and related 2D materials is that they can be produced on large scales by liquid phase exfoliation. The exfoliation process shall be considered as a particular fragmentation process, where the 2D character of the exfoliated objects will influence significantly fragmentation dynamics as compared to standard materials. Here, we used automatized image processing of Atomic Force Microscopy (AFM) data to measure, one by one, the exact shape and size of thousands of nanosheets obtained by exfoliation of an important 2D-material, boron nitride, and used different statistical functions to model the asymmetric distribution of nanosheet sizes typically obtained.Being the resolution of AFM much larger than the average sheet size, analysis could be performed directly at the nanoscale and at the single sheet level. We find that the size distribution of the sheets at a given time follows a log-normal distribution, indicating that the exfoliation process has a "typical" scale length that changes with time and that exfoliation proceeds through the formation of a distribution of random cracks that follow Poisson statistics. The validity of this model implies that the size distribution does not depend on the different preparation methods used, but is a common feature in the exfoliation of this material and thus probably for other 2D materials.The huge scientic and technological interest for graphene has triggered in the last few years the development of a wide range of techniques to produce and process nanosheets that, having nanometric thickness and mesoscopic lateral size, shall be considered as quasi 2-dimensional (2D) objects. Besides their novel properties, even the way these 2D sheets are produced in solution, by exfoliation, 1 is an original process, still not completely understood.The exfoliation of a 2D object from a 3D bulk material is a process spanning from the nano-to meso-scale due to bubble cavitation, intercalation and disruptive fragmentation, as we described in recent work. 2 Exfoliation always yields a polydispersed range of nanosheet thickness and lateral size. When characterizing these 2D sheet solutions, their average size and size standard deviation are commonly reported, oen assuming that their size follows a "Gaussian" (a.k.a. "normal") distribution. Conversely, the experimental data show that the size distribution of these materials is highly asymmetric and nonGaussian.It is noteworthy that this asymmetry in size distribution shall be observed in very different systems such as the distribution of chemical elements in rocks, the species abundance in biology, the lengths of latent periods of infectious diseases in medicine, and the distribution of galaxies in astronomy (Fig. 1). 3,4 A better modelling of the size distribution of 2D materials is needed both from a fundamental point of view (to understand the exfoliation mechanism) and from a technological point of view (to improve the metrology of 2D materials for applications and quality control).Here, we ...