A stochastic and analytical model of hierarchical fragmentation: The fragmentation of gas structures into young stellar objects in the interstellar medium

Context. Molecular clouds are the most important incubators of young stars clustered in various stellar structures whose spatial extension can vary from a few AU to several thousand AU. Although the reality of these stellar systems has been established, the physical origin of their multiplicity remains an open question. Aims. Our aim was to characterise these stellar groups at the onset of their formation by quantifying both the number of stars they contain and their mass using a hierarchical fragmentation model of the natal molecular cloud. Methods. We developed a stochastic and predictive model that reconciles the continuous multi-scale structure of a fragmenting molecular cloud with the discrete nature of the stars that are the products of this fragmentation. In this model a gas structure is defined as a multi-scale object associated with a subregion of a cloud. Such a structure undergoes quasi-static subfragmentation until star formation. This model was implemented within a gravo-turbulent fragmentation framework to analytically follow the fragmentation properties along spatial scales using an isothermal and adiabatic equations of state (EOSs). Results. We highlighted three fragmentation modes depending on the amount of fragments produced by a collapsing gas structure, namely a hierarchical mode, a monolithic mode, and a mass dispersal mode. Using an adiabatic EOS we determined a characteristic spatial scale where further fragmentation is prevented, around a few tens of AU. We show that fragmentation is a self-regulated process as fragments tend to become marginally unstable following a M proportional to R Bonnor-Ebert-like mass-size profile. Supersonic turbulent fragmentation structures the cloud down to R approximate to 0.1 pc, and gradually turns into a less productive Jeans-type fragmentation under subsonic conditions so hierarchical fragmentation is a scale dependant process. Conclusions. Our work suggests that pre-stellar objects resulting from gas fragmentation, have to progressively increase their accretion rate in order to form stars. A hierarchical fragmentation scenario is compatible with both the multiplicity of stellar systems identified in Taurus and the multi-scale structure extracted within NGC 2264 molecular cloud. This work suggests that hierarchical fragmentation is one of the main mechanisms explaining the presence of primordial structures of stellar clusters in molecular clouds.