The mass function of clumps observed in molecular clouds raises interesting theoretical issues, especially in its relation to the stellar initial mass function. We propose a statistical model of the mass function of prestellar cores (CMF), formed in self-gravitating isothermal clouds at a given stage of their evolution. The latter is characterized by the mass-density probability distribution function ($\rho$-PDF), which is a power-law with slope $q$. The variety of MCs is divided in ensembles according to the PDF slope and each ensemble is represented by a single spherical cloud. The cores are considered as elements of self-similar structure typical for fractal clouds and are modeled by spherical objects populating each cloud shell. Our model assumes relations between size, mass and density of the statistical cores. Out of them a core mass-density relationship $\rho\propto m^x$ is derived where $x=1/(1+q)$. We found that $q$ determines the existence or non-existence of a threshold density for core collapse. The derived general CMF is a power law of slope $-1$ while the CMF of gravitationally unstable cores has a slope $(-1 + x/2)$, comparable with the slopes of the high-mass part of the stellar initial mass function and of observational CMFs.

We present a new approach to extract the power-law part of a density/column-density probability density function (rho-pdf/N-pdf) in star-forming clouds. It is based on the mathematical method bPLFIT of Virkar & Clauset (2014) and assesses the power-law part of an arbitrary distribution, without any assumptions about the other part of this distribution. The slope and deviation point are derived as averaged values as the number of bins is varied. Neither parameter is sensitive to spikes and other local features of the tail. This adapted bPLFIT method is applied to two different sets of data from numerical simulations of star-forming clouds at scales 0.5 and 500 pc and displays rho-pdf and N-pdf evolution in agreement with a number of numerical and theoretical studies. Applied to Herschel data on the regions Aquila and Rosette, the method extracts pronounced power-law tails, consistent with those seen in simulations of evolved clouds.

In this paper we continue to investigate the energy conservation equation obtained in our previous work. We set ourselves three new goals. The first one is to rewrite the main equations in terms of density profile in order to give more physical insight. The second one is to investigate the significance of two new terms in the energy conservation equation. They originate from the gravity of the outer shells of cloud and the masses outer to the cloud, respectively. The third goal is to investigate the main equation in the case when the kinetic turbulent term scales according to Larson's law and it is independent, formally, of the accretion, in contrast to the previous work. The combination of supersonic turbulence and spherical symmetry raises a caveat which is commented in our conclusions. We obtained two solutions for the density profile. They scale with slopes -2 and -3/2, respectively. The energy balance for the second solution is the same as in the previous paper: this is a free-fall. For the first solution there are two cases. The first one: if the turbulent term does not scale, then it could be important for the energy balance of the cloud. The second one: if the turbulent term does scale, then it is not important for the energy balance of the cloud. The two new gravitational terms don't affect the existence of the two solutions, but the gravitation of the outer masses calibrate the energy balance for the first solution.

Spatial association of clumps from different tracers turns out to be a valuable tool to determine the physical properties of molecular clouds. It provides a reliable estimate for the $X$-factors, serves to trace the density of clumps seen in column densities only and allows to measure the velocity dispersion of clumps identified in dust emission. We study the spatial association between clump populations, extracted by use of the GAUSSCLUMPS technique from $^{12}$CO (1-0), $^{13}$CO (1-0) line maps and Herschel dust-emission maps of the star-forming region Rosette, and analyse their physical properties. All CO clumps that overlap with another CO or dust counterpart are found to be gravitationally bound and located in the massive star-forming filaments of the molecular cloud. They obey a single mass-size relation $M_{\rm cl}\propto R_{\rm cl}^\gamma$ with $\gamma\simeq3$ (implying constant mean density) and display virtually no velocity-size relation. We interpret their population as low-density structures formed through compression by converging flows and still not evolved under the influence of self-gravity. The high-mass parts of their clump mass functions are fitted by a power law ${\rm d}N_{\rm cl}/{\rm d}\,\log M_{\rm cl}\propto M_{\rm cl}^{\Gamma}$ and display a nearly Salpeter slope $\Gamma\sim-1.3$. On the other hand, clumps extracted from the dust-emission map exhibit a shallower mass-size relation with $\gamma=2.5$ and mass functions with very steep slopes $\Gamma\sim-2.3$ even if associated with CO clumps. They trace density peaks of the associated CO clumps at scales of a few tenths of pc where no single density scaling law should be expected.

We have set ourselves the task of obtaining the probability distribution function of the mass density of a self-gravitating isothermal compressible turbulent fluid from its physics. We have done this in the context of a new notion: the molecular clouds ensemble. We have applied a new approach that takes into account the fractal nature of the fluid. Using the medium equations, under the assumption of steady state, we show that the total energy per unit mass is an invariant with respect to the fractal scales. As a next step we obtain a nonlinear integral equation for the dimensionless scale Q which is the third root of the integral of the probability distribution function. It is solved approximately up to the leading-order term in the series expansion. We obtain two solutions. They are power-law distributions with different slopes: the first one is -1.5 at low densities, corresponding to a equilibrium between all energies at a given scale, and the second one is -2 at high densities, corresponding to a free fall at small scales.

We present a method to derive the density scaling relation $\langle n\rangle \propto L^{-\alpha}$ in regions of star formation or in their turbulent vicinities from straightforward binning of the column-density distribution ($N$-pdf). The outcome of the method is studied for three types of $N$-pdf: power law ($7/5\le\alpha\le5/3$), lognormal ($0.7\lesssim\alpha\lesssim1.4$) and combination of lognormals. In the last case, the method of Stanchev et al. (2015) was also applied for comparison and a very weak (or close to zero) correlation was found. We conclude that the considered `binning approach' reflects rather the local morphology of the $N$-pdf with no reference to the physical conditions in a considered region. The rough consistency of the derived slopes with the widely adopted Larson's (1981) value $\alpha\sim1.1$ is suggested to support claims that the density-size relation in molecular clouds is indeed an artifact of the observed $N$-pdf.

We introduce the concept of a class of equivalence of molecular clouds represented by an abstract spherically symmetric, isotropic object. This object is described by use of abstract scales in respect to a given mass density distribution. Mass and average density are ascribed to each scale and thus are linked to the density distribution: a power-law type and an arbitrary continuous one. In the latter case, we derive a differential relationship between the mean density at a given scale and the structure parameter which defines the mass-density relationship. The two-dimensional (2D) projection of the cloud along the line of sight is also investigated. Scaling relations of mass and mean density are derived in the considered cases of power-law and arbitrary continuous distributions. We obtain relations between scaling exponents in the 2D and 3D cases. The proposed classes of equivalence are representative for the general structure of real clouds with various types of column-density distributions: power law, lognormal or combination of both.

Recently Stanchev et al. (2015) proposed a technique to derive density scaling relations in a star-forming region from analysis of the probability distribution function of column density. We address the possible dependence of the outcome on the selection of probe zones, applying the method to Planck dust-opacity data on Orion A. The derived steep scaling relation of mean density with index -1.6 in the molecular cloud (so called `Central filament') points to its self-gravitating nature. The result is reproduced also for large parts of the clouds' vicinity which indicates major role of gravity in the energy balance of the entire star-forming region.

We present a model for describing the general structure of molecular clouds (MCs) at early evolutionary stages in terms of their mass-size relationship. Sizes are defined through threshold levels at which equipartitions between gravitational, turbulent and thermal energy $|W| \sim f(E_{\rm kin} + E_{\rm th})$ take place, adopting interdependent scaling relations of velocity dispersion and density and assuming a lognormal density distribution at each scale. Variations of the equipartition coefficient $1\le f\le 4$ allow for modelling of star-forming regions at scales within the size range of typical MCs ($\gtrsim$4 pc). Best fits are obtained for regions with low or no star formation (Pipe, Polaris) as well for such with star-forming activity but with nearly lognormal distribution of column density (Rosette). An additional numerical test of the model suggests its applicability to cloud evolutionary times prior to the formation of first stars.

We present an analysis of probability distribution functions (pdfs) of column density in different zones of the star-forming region Perseus and its diffuse environment based on the map of dust opacity at 353 GHz available from the Planck archive. The pdf shape can be fitted by a combination of a lognormal function and an extended power-law tail at high densities, in zones centred at the molecular cloud Perseus. A linear combination of several lognormals fits very well the pdf in rings surrounding the cloud or in zones of its diffuse neighbourhood. The slope of the mean density scaling law $\langle\rho\rangle_L \propto L^\alpha$ is steep ($\alpha=-1.93$) in the former case and rather shallow ($\alpha=-0.77\pm0.11$) in the rings delineated around the cloud. We interpret these findings as signatures of two distinct physical regimes: i) a gravoturbulent one which is characterized by nearly linear scaling of mass and practical lack of velocity scaling; and ii) a predominantly turbulent one which is best described by steep velocity scaling and by invariant for compressible turbulence $\langle\rho\rangle_L u_L^3/L$, describing a scale-independent flux of the kinetic energy per unit volume through turbulent cascade. The gravoturbulent spatial domain can be identified with the molecular cloud Perseus while a relatively sharp transition to predominantly turbulent regime occurs in its vicinity.

Volume and column density PDFs in molecular clouds are important diagnostics for understanding their general structure. We developed a novel approach to trace the cloud structure by varying the lower PDF cut-off and exploring a suggested mass-density relationship with a power-law index $x^\prime$. The correspondence of x' as a function of spatial scale to the slope of the high-density PDF tail is studied. To validate the proposed model, we use results from hydrodynamical simulations of a turbulent self-gravitating cloud and recent data on dust continuum emission from the Planck mission.

The statistical approach for derivation of the clump mass function (ClMF) developed by Donkov, Veltchev & Klessen is put to observational test through comparison with mass distributions of clumps from molecular emission and dust continuum maps of Galactic cloud complexes, obtained by various authors. The results indicate gravitational boundedness of the dominant clump population, with or without taking into account the contribution of their thermal and magnetic energy. The ClMF can be presented by combination of two power-law functions separated by a characteristic mass from about ten to hundreds solar masses. The slope of the intermediate-mass ClMF is shallow and nearly constant (-0.25 \gtrsim \Gamma_{IM} \gtrsim -0.55) while the high-mass part is fitted by models that imply gravitationally unstable clumps and exhibit slopes in a broader range (-0.9 \gtrsim \Gamma_{IM} \gtrsim -1.6), centered at the value of the stellar initial mass function (\Gamma_{HM} \gtreqless -1.3).

The statistical approach for description of molecular cloud substructure, proposed by Donkov, Veltchev and Klessen (2011, 2012), allows for alternative models, operating with different type of objects: an ensemble of clumps or a larger cloudlet. We demonstrate briefly the predictive power of both models, applied to molecular emission and dust extinction studies of Galactic clouds.

We develop a statistical approach for description of dense structures (cores) in molecular clouds that might be progenitors of stars. Our basic assumptions are a core mass-density relationship and a power-law density distribution of these objects as testified by numerical simulations and observations. The core mass function (CMF) was derived and its slope in the high-mass regime was obtained analytically. Comparisons with observational CMFs in several Galactic clouds are briefly presented.

We derive the mass function of condensations (clumps) which were formed through a turbulent cascade over a range of spatial scales $L\le20$ pc during early, predominantly turbulent evolution of a molecular cloud. The approach rests upon the assumption of a statistical clump mass-density relationship $n\propto m^x$ with a scale dependence of the exponent $x$ obtained from equipartition relations between various forms of energy of clumps. The derived clump mass function (ClMF) could be represented by series of 2 or 3 power laws, depending on the chosen equipartition relation, the velocity scaling index and the type of turbulent forcing. The high-mass ClMF exhibits an average slope $\Gamma\simeq-1$, typical for fractal clouds, whereas its intermediate-mass part is shallower or flattened, in agreement with some observational studies.

We study the mass-density relationship n ~ m^x in molecular cloud condensations (clumps), considering various equipartition relations between their gravitational, kinetic, internal and magnetic energies. Clumps are described statistically, with a density distribution that reflects a lognormal probability density function (pdf) in turbulent cold interstellar medium. The clump mass-density exponent $x$ derived at different scales $L$ varies in most of the cases within the range $-2.5\lesssim x \lesssim-0.2$, with a pronounced scale dependence and in consistency with observations. When derived from the global size-mass relationship m ~ l^{\gamma_{glob}} for set of clumps, generated at all scales, the clump mass-density exponent has typical values $-3.0\lesssim x(\gamma_{glob}) \lesssim -0.3$ that depend on the forms of energy, included in the equipartition relations and on the velocity scaling law whereas the description of clump geometry is important when magnetic energy is taken into account.