Issue 
A&A
Volume 571, November 2014



Article Number  A95  
Number of page(s)  22  
Section  Interstellar and circumstellar matter  
DOI  https://doi.org/10.1051/00046361/201423647  
Published online  14 November 2014 
On the probability distribution function of the mass surface density of molecular clouds. II.^{⋆}
Canadian Institute for Theoretical Astrophysics, University of
Toronto, 60 St. George
Street, Toronto,
ON M5S3H8, Canada
email:
fischera@cita.utoronto.ca
Received: 14 February 2014
Accepted: 11 September 2014
The probability distribution function (PDF) of the mass surface density of molecular clouds provides essential information about the structure of molecular cloud gas and condensed structures out of which stars may form. In general, the PDF shows two basic components: a broad distribution around the maximum with resemblance to a lognormal function, and a tail at high mass surface densities attributed to turbulence and selfgravity. In a previous paper, the PDF of condensed structures has been analyzed and an analytical formula presented based on a truncated radial density profile, ρ(r) = ρ_{c}/ (1 + (r/r_{0})^{2})^{n/ 2} with central density ρ_{c} and inner radius r_{0}, widely used in astrophysics as a generalization of physical density profiles. In this paper, the results are applied to analyze the PDF of selfgravitating, isothermal, pressurized, spherical (BonnorEbert spheres) and cylindrical condensed structures with emphasis on the dependence of the PDF on the external pressure p_{ext} and on the overpressure q^{1} = p_{c}/p_{ext}, where p_{c} is the central pressure. Apart from individual clouds, we also consider ensembles of spheres or cylinders, where effects caused by a variation of pressure ratio, a distribution of condensed cores within a turbulent gas, and (in case of cylinders) a distribution of inclination angles on the mean PDF are analyzed. The probability distribution of pressure ratios q^{1} is assumed to be given by P(q^{1}) ∝ q^{−k1}/ (1 + (q_{0}/q)^{γ})^{(k1 + k2) /γ}, where k_{1}, γ, k_{2}, and q_{0} are fixed parameters. The PDF of individual spheres with overpressures below ~100 is well represented by the PDF of a sphere with an analytical density profile with n = 3. At higher pressure ratios, the PDF at mass surface densities Σ ≪ Σ(0), where Σ(0) is the central mass surface density, asymptotically approaches the PDF of a sphere with n = 2. Consequently, the powerlaw asymptote at mass surface densities above the peak steepens from P_{sph}(Σ) ∝ Σ^{2} to P_{sph}(Σ) ∝ Σ^{3}. The corresponding asymptote of the PDF of cylinders for the large q^{1} is approximately given by P_{cyl}(Σ) ∝ Σ^{−4/3}(1 − (Σ/Σ(0))^{2/3})^{−1/2}. The distribution of overpressures q^{1} produces a powerlaw asymptote at high mass surface densities given by ∝ Σ^{−2k2 − 1} (spheres) or ∝ Σ^{−2k2} (cylinders).
Key words: ISM: clouds / methods: statistical / ISM: structure / methods: analytical
Appendices are available in electronic form at http://www.aanda.org
© ESO, 2014
1. Introduction
The observed probability distribution function (PDF) of the mass surface density of molecular clouds suggests a combination of two separate components that produce a broad distribution around the peak, which is generally modeled using a lognormal function and a tail at high mass surface densities that often has a powerlaw form (e.g., Kainulainen et al. 2009; Kainulainen & Tan 2013; Schneider et al. 2012, 2013). The broad distribution around the maximum is typically attributed to turbulence of the molecular gas, the tail to selfgravitating structures. The relative amount of gas in gravitationally bound structures seems to be indicative of the star formation rate in the cloud. The PDF of starforming molecular clouds like Taurus or Ophiuchus, for example, is characterized by a strong tail that seems to be very low or absent in clouds with no apparent star formation, such as the “Coalsack” (Kainulainen et al. 2009).
Images of the column densities of molecular clouds (Schneider et al. 2012; Kainulainen & Tan 2013) show that the high column density values are not randomly distributed within the cloud, but are indeed related to small regions, as expected for cold condensed structures surrounded by warmer or more turbulent gas. In a previous paper (Fischera 2014a, Paper I) the PDFs of spherical and cylindrical condensed structures have been analyzed assuming a truncated smooth radial density profile used in astrophysics as a generalization of physical density profiles. It was found that the PDF of such a profile can be described by a simple implicit analytical function. In Paper I, we showed how the geometric shape, the radial density profile, and the ratio ρ_{c}/ρ(r_{cl}) of the central density and the density at the peripheral regions affect the functional form of the PDF. For example, we found that the PDF of a sphere is truncated, while the PDF of a cylinder has a pole at the highest mass surface density.
In this paper, the results of Paper I are applied to analyze the PDF of isothermal selfgravitating pressurized spheres and cylinders as an approximation of individual condensations in giant molecular clouds. In the case of the global statistical properties of molecular clouds, the condensed structures probably show a variety of different physical conditions. The consequences of a number of distributions on the mean PDF are therefore also addressed.
The paper is divided into two main sections. In the first part (Sect. 2), the properties of the PDFs of individual spheres and cylinders are discussed. We show how the highest position and the asymptotes at low and high mass surface densities are affected by the pressure ratio p_{c}/p_{ext} and the external pressure p_{ext} surrounding the structures. The second part (Sect. 3) is about the properties of a mean PDF for a variety of distributions as they may apply to molecular clouds. In Sect. 3.2, we analyze the mean PDF of an ensemble of spheres and cylinders for a distribution of gravitational states. In Sect. 3.3, the effect caused by a distribution of cylinders with a variety of inclination angles is addressed, and in Sect. 3.4, we consider a mean PDF of cores that are uniformly distributed within a turbulent medium. The results are discussed in Sect. 4 and compared with observed PDFs, as derived by Kainulainen et al. (2009). The paper is summarized in Sect. 5.
2. Probability distribution function of the mass surface density of condensed structures
2.1. Density profile of the isothermal condensed structures
In Paper I, we have analyzed the onepoint statistic for condensed spheres and cylinders assuming an analytical density profile given by (1)where ρ_{c} is the central density and r_{0} the inner radius. In the following we show how the density profile is related to the density profile of selfgravitating isothermal spheres, known as BonnorEbert spheres (Ebert 1955; Bonnor 1956), and cylinders.
In isothermal clouds, the gas pressure and density are proportional, p = Kρ, where the proportionality constant is given by K = kT/ (μm_{H}), where k, T, μ, and m_{H} are the Boltzmann constant, the effective temperature, the mean molecular weight, and the mass of a hydrogen atom. For an idealized gas we have , where c_{s} is the sound speed.
To analyze isothermal selfgravitating clouds, it is convenient to introduce a unitfree radius θ = rA. The constant is given by A^{2} = (4πGρ_{c}) /K, where G is the gravitational constant. The radius r_{0} is related to the effective temperature T and the central density through where ξ_{n} is a fixed appropriate constant for given n so that . The parameters for the considered density profiles are listed in Table 1. For BonnorEbert spheres the analytical density profiles are only representations of the correct density profile at low and and in the limit of high θ. The corresponding constants ξ_{n} are chosen, as described in Sect. 2.1.1, to match the correct density profiles where the approximations are valid.
The pressure ratio q^{1} = p_{c}/p_{ext} of the central pressure p_{c} and the external pressure p_{ext} are referred to as “overpressure”. As a reference pressure we assume a mean interstellar medium (ISM) pressure p_{ext}/k = 2 × 10^{4} K cm^{3}, as previously suggested by Curry & McKee (2000) based on the study of Boulares & Cox (1990), where the pressure is attributed to thermal, turbulent, and magnetic pressure components. The same value has recently been inferred from the physical parameters and profiles of filaments (Fischera & Martin 2012b,a).
Fig. 1 Normalized radial density profiles of some cases considered in this study given as a function of the unitfree radius θ = rA where A^{2} = (4πGρ_{c}) /K. The filled symbol marks the condition for a critically stable sphere where θ_{crit} ≈ 6.451 and ρ_{c}/ρ(θ_{crit}) ≈ 14.04. The curves of the analytical density profile given by Eq. (1) are labeled with their asymptotic behavior for . 

Open with DEXTER 
For completeness, we discuss in Appendix B the PDF of the local density for spheres and cylinders.
Parameters of analytical density profiles.
2.1.1. BonnorEbert sphere
The density profile of an isothermal selfgravitating sphere, or BonnorEbert sphere, is given by (2)where ω = φ/K is the unitfree potential of the gravitational potential φ determined by the LaneEmden equation, (3)Spheres have a wellknown critically stable configuration with an overpressure above which the sphere becomes gravitationally unstable under compression (Ebert 1955; Bonnor 1956; McCrea 1957; Curry & McKee 2000; Fischera & Dopita 2008; Fischera 2011; Fischera & Martin 2012b; see Appendix A.1.1). The critically stable configuration is characterized through an overpressure of p_{c}/p_{ext} ~ 14.04 at θ_{crit} = 6.451. However, theoretically, equilibrium solutions with higher overpressures can also be constructed. As in previous papers (Fischera 2011; Fischera & Martin 2012b), we refer to spheres with overpressures below the critical value as subcritical and spheres with higher overpressure as supercritical.The term “supercritical” is often used in the literature to characterize clouds with masses above the critically stable mass. The supercritical spheres in this paper have masses that are lower than the mass of a critically stable sphere if we let the temperature and external pressure remain constant, as discussed by (Fischera & Martin 2012b; see also Appendix A).
As shown in Fig. 1, the profile of the BonnerEbert sphere below θ = 10 or ρ(r) > ~ 0.01 ρ_{c} is well represented by an analytical density profile as given by Eq. (1) with n = 3. For the constant ξ_{3} we have chosen a value that is consistent with the overpressure and the size of a critically stable sphere: (4)This value is used throughout the paper.
At large sizes θ ≫ 10 the density profile of the BonnorEbert sphere approaches asymptotically a powerlaw profile ρ(θ) /ρ_{c} ∝ C/θ^{2} where C is a constant. Inserting the powerlaw profile in Eq. (3), we find C = 2. The same asymptotic behavior is found for an analytical density profile as given in Eq. (1) with n = 2 and ξ_{2} = 2 in the limit of large sizes . We consider this profile as the asymptotic profile of BonnorEbert spheres for high overpressures. Although stable configurations with overpressure above the critical value probably do not exist, the profiles of sub and supercritical solutions seem to reflect within the observed uncertainties the profiles not only of stable clouds, but also seem to mimic the density profile of collapsing clouds (Kandori et al. 2005; Keto & Caselli 2010; Keto et al. 2014).
Physical parameters of BonnorEbert spheres such as mass and radius are discussed in more detail in Appendix A in comparison with spheres with a density profile with n = 3 or n = 2. In Appendix B.1, we show the variation of the local density PDF with overpressure.
2.1.2. Selfgravitating cylinders
For n = 4 the density profile is identical to the profile of an isothermal selfgravitating cylinder with ξ_{4} = 8 (Stodólkiewicz 1963; Ostriker 1964). The physical parameters of pressurized cylinders have been studied in great detail by Fischera & Martin (2012a). Here, we wish to describe the most important properties for understanding the PDF of the mass surface density of cylinders.
As an important characteristic, a highest possible mass line density (M/l)_{max} = 2K/G exists for cylinders, which corresponds to a cylinder with infinite overpressure and infinitely small size. As shown in Fischera & Martin (2012a), for example, physical properties such as size, full width at half maximum (FWHM), and stability considerations can be expressed through the mass ratio f_{cyl} = (M/l)/(M/l)_{max} of the mass line density and the highest possible mass line density. It is related to the pressure ratio q by (5)Cylinders with high mass ratio f_{cyl} have a steep density profile ρ(r) ∝ r^{4} at radii . However, Fig. 1 shows that this asymptotic profile is only established at a relatively large radius θ ≫ ~ 10, which approximately corresponds to cylinders with an overpressure q^{1}> ~ 100 or equally to cylinders with mass ratio f_{cyl}> ~ 0.9. Measurements of the density profile of observed filaments in general indicate a flatter density profile more consistent with r^{2}, as discussed by Fiege & Pudritz (2000), suggesting significantly lower mass ratios. Indications for a filament with a high mass ratio and a possible steep density profile are found for a filament in the cloud IC 5146 (Fischera & Martin 2012a)^{1}.
2.2. Profile of the mass surface density
The profile of the mass surface density of a pressurized, isothermal, selfgravitating sphere or cylinder seen at inclination angle i is in general given by (6)where θ_{⊥} and θ_{cl} are the projected unitfree radius and cloud radius. λ is a power index with λ = 0 for spheres and λ = 1 for cylinders. In this convention the inclination angle is i = 0° for a cylinder seen from the side. Profiles for spheres and cylinders for a range of overpressures can be found in the work of Fischera & Martin (2012b). It is noticeable that the profile does not depend on the gas temperature. Consequently, isothermal clouds in a certain pressure region with the same overpressure will have the same profile Σ(θ) independently of their mass or mass line density, which vary for given overpressure with or (Appendix A.1). Their profiles only vary in terms of size as for fixed overpressure (Appendix A.2).
The profiles of the mass surface density for a sphere and a cylinder with a truncated radial density profile as given by Eq. (1) were presented in Paper I. It was found to be convenient to introduce a unitfree mass surface density X_{n} defined by (7)The unitfree mass surface density is given by (8)where y_{n} = (1 − q^{2 /n})(1 − x^{2}) and where x = r_{⊥}/r_{cl} is the ratio of the projected radius r_{⊥} and the cloud radius r_{cl}, which we refer to as the normalized impact parameter. The upper limit of the integral is . For density profiles with n = 2, 3, and 4 the profiles are simple analytical functions and can be found in Paper I.
For isothermal selfgravitating pressurized clouds, we can replace the inner radius using and the mass surface density becomes (9)We note that for n = 2 the relation is independent of the pressure ratio q.
The mass surface density profiles through BonnorEbert spheres is shown in Fig. 2. At small projected radius θ_{⊥} ≪ 10 the profile is dominated by the innermost part of the radial density profile and we can use the profile of the analytical density profile with n = 3 as approximation. At large projected radius θ_{⊥} the profile of a BonnorEbert sphere fluctuates asymptotically toward a power law Σ ∝ 1 /θ_{⊥}. With the exception of small θ_{⊥}, the mass surface densities of a BonnorEbert sphere is approximately described by a sphere with an analytical density profile with n = 2.
Fig. 2 Profiles of the mass surface density through BonnorEbert spheres and spheres with an analytical density profile as given in Eq. (1). The gray curves are profiles for different cloud radii θ_{cl}. For BonnorEbert spheres the cloud radius corresponds to q = q_{crit} (a), q = 10^{2}q_{crit} (b), and q = 10^{4}q_{crit} (c), where q_{crit} = 1/14.04. For each assumed cloud radius the central pressure p_{c} in the considered spheres is the same. The profiles are normalized by the corresponding central mass surface density through a BonnorEbert sphere. 

Open with DEXTER 
To compare the PDF of BonnorEbert spheres and the PDF of the approximations it is convenient to measure the mass surface density in units of (10)For the truncated analytical density profile with index n the mass surface density becomes (11)It is common to measure the projected density in visual extinction A_{V} assuming the optical properties derived for the diffuse interstellar medium. To provide an estimate of the mass surface density in units of extinction A_{V} we assumed a gastodust ratio of N_{H}/E(B − V) = 5.8 × 10^{21} cm^{2} mag^{1} (Bohlin et al. 1978) and an absolutetorelative extinction R_{V} of 3.1 (Fitzpatrick 1999) so that (12)where is the assumed mean atomic gas mass in units of a hydrogen atom with mass m_{H}.
The critically stable BonnorEbert sphere (q_{crit} = 1/14.04) has a normalized central mass surface density of which corresponds to a central visual extinction of .
As a special important case we give the behavior of the central mass surface density in the limit of high overpressures. In this regime the central mass surface density through BonnorEbert spheres is approximately given by (Eq. (C.4), Fischera & Martin 2012b) (13)For spheres and cylinders with analytical density profiles as given in Eq. (1), the upper limit of the integral in Eq. (8) behaves as u_{max} → ∞ and the central mass surface density becomes approximately (14)where (15)is the betafunction and where Γ(x) is the Γfunction. For a sphere with n = 3, we obtain with the constant ξ_{3} = 8.63 in the limit of high overpressures (16)which is only 3% lower than the correct asymptotic value of the BonnorEbert sphere.
2.3. Probability distribution function of BonnorEbert spheres
The probability of measuring a mass surface density Σ of an individual sphere in the range from Σ...Σ + dΣ is given by (17)where P(r_{⊥}) is the probability of measuring a projected radius r_{⊥}. For a sphere, we have .
The derivative of the mass surface density of a cloud with a radial density profile is in general given by (18)where r_{cl} is the cloud radius and r_{⊥} the projected radius. The integral only applies to the first term in brackets. After substituting the density with ρ(θ) = ρ_{c}e^{−ω(θ)}, we obtain (19)where we have replaced r, r_{⊥}, and r_{cl} through the corresponding unitfree radii θ = rA, θ_{⊥} = r_{⊥}A, and θ_{cl} = r_{cl}A with . Inserting the expression in Eq. (17) provides for the PDF of the BonnorEbert sphere (20)where e^{−ω(zcl)} = q, and where we have again replaced the radius and the projected radius by the corresponding unitfree expressions and Kρ_{c} with p_{ext}/q.
The PDF of BonnorEbert spheres is compared with spheres with an analytical density profile as given in Eq. (1) with n = 3 and n = 2. The PDF of spheres with such a profile with arbitrary n is given by the implicit function (Paper I) (21)where (22)It is straightforward to show that for the approximations of the radial density profile of BonnorEbert spheres the two different expressions 20 and 21 with 22 for the PDF are identical (Appendix D).
Fig. 3 Probabilities of a number of mass surface densities through BonnorEbert spheres shown as a function of the central mass surface density. The curves are labeled by the corresponding mass surface density given in units of . The probabilities are compared with those of analytical density profiles as given in Eq. (1) with n = 3 and n = 2. The probabilities for n = 3 are only shown for where the analytical density profile becomes a valid approximation of the density profile of a BonnorEbert sphere. The overpressures given in the upper axis are only valid for BonnorEbert spheres. The gray lines show the asymptotic probability for BonnorEbert spheres in the limit of infinite overpressure. The vertical dotted line marks the location of the critically stable BonnorEbert sphere. 

Open with DEXTER 
The behavior of the PDF of BonnorEbert spheres with increasing overpressure is shown in Fig. 3. At low overpressures the probabilities are well described by a sphere with an analytical density profile with n = 3. The probability for increasing overpressure fluctuates asymptotically toward the value of a sphere with an analytical density profile with n = 2 in the limit of infinite overpressure.
It is common in observational studies of molecular clouds to measure the logarithmic PDF P(lnΣ) of the mass surface density by estimating the number of mass surface densities Σ within a logarithmic mass surface density bin ΔlnΣ. This is motivated by the large scales in mass surface densities involved, but also by the fact that turbulent cloud structures should show a lognormal density distribution (see Sect. 3.4). As P(Σ) dΣ = ΣP(Σ) dΣ/Σ = ΣP(Σ) dlnΣ we have P(lnΣ) = ΣP(Σ). In Paper I, P(Σ) is therefore referred to as linear PDF and ΣP(Σ) as logarithmic PDF of the mass surface density Σ. It follows from Eqs. (6) and (20) that the amplitude of the logarithmic PDF of BonnorEbert spheres is independent of the external pressure. The same is true for clouds with analytical density profiles as X_{n}P(X_{n}) = Σ_{n}P(Σ_{n}). The statistical properties of the clouds are thus visualized using the logarithmic PDF.
Fig. 4 PDF of the mass surface density of BonnorEbert spheres embedded in a medium with p_{ext}/k = 2 × 10^{4} K cm^{3} for different gravitational states (solid black curves). The curves correspond to overpressures with i = −1,0,1 where is the overpressure of a critically stable BonnorEbert sphere. The PDFs are directly compared with those of spheres with a simple density profile (Eq. (1)) with n = 3 (dotted curves). The black dashed lines are asymptotes of the PDF Σ_{3}P_{sph}(Σ_{3}) in the limit of high mass surface densities (Eq. (39)). The PDF of a BonnorEbert sphere in the limit of infinite overpressure is shown as a gray curve and is derived assuming a sphere with an analytical density profile with n = 2. The gray dashed lines are the corresponding powerlaw asymptotes in the limit of low (Eq. (34)) and high (Eq. (40)) mass surface densities. 

Open with DEXTER 
Fig. 5 Similar to Fig. 4, but for a wider range of pressure ratios. We show the transition from a PDF of a critically stable BonnorEbert sphere to the PDFs of highly supercritical spheres. The overpressure is varied in units of the critical value . Also shown are the corresponding PDFs for a sphere with the analytical density profile with n = 2 and n = 3. The powerlaw asymptotes close to the truncation point (black dashed lines) are scaled to the PDF of the Bonnor Ebert spheres at central mass surface density. 

Open with DEXTER 
Fig. 6 Normalized impact parameter (lefthand panel) and mass surface density (righthand panel) at the maximum of the PDF of the mass surface densities of BonnorEbert spheres, given as a function of the central mass surface density. The gray and black curves correspond to the logarithmic and the linear pdf. The values are compared with the corresponding impact parameters and mass surface densities of pressurized spheres with a density profile as given in Eq. (1) with n = 2 and n = 3. At low central mass surface densities the maximum of the PDF is at impact parameter [ x ] _{max} = 0 and the curves in the righthand panel follow the long dashed line where [ Σ ] _{max} = Σ(0). The location of the critically stable BonnorEbert sphere is marked by a vertical dotted line and the corresponding normalized impact parameters and mass surface densities at the PDF maximum as filled circles. The horizontal dotted lines in the left panel are the asymptotic values of the impact parameters of the PDF maximum of BonnorEbert spheres in the limit of high overpressure. The overpressures shown in the upper axis are only valid for the BonnorEbert sphere. 

Open with DEXTER 
The PDFs of BonnorEbert spheres for a series of overpressures are shown in Figs. 4 and 5. We might interpret the series of curves as those of an initially stable cloud that becomes gravitationally unstable through cooling or mass accretion and finally collapses.
The PDF shows characteristics discussed in Paper I for spheres with a truncated analytical density profiles as given in Eq. (1), for example, the cutoff at the central mass surface density, the broad maximum above a certain overpressure, and the powerlaw asymptote in the limit of low mass surface densities where ΣP(Σ) ∝ Σ^{2}. For high overpressures, the PDF shows oscillations at high mass surface densities related to the oscillations of the radial density profile (Fig. 1) or the profile of the mass surface densities (Fig. 2). With decreasing mass surface density, the same PDFs fluctuate asymptotically toward the PDF of a sphere with an n = 2profile. A quantitative analysis of the highest position and the properties at low and high mass surface densities of the PDF of BonnorEbert spheres is given in Sects. 2.3.1−2.3.3.
2.3.1. Maxima position
In this section we quantitatively analyze the highest position of the PDF of BonnorEbert spheres. The normalized impact parameters [ x ] _{max} and the mass surface density [ Σ ] _{max} at the maxima of the linear and logarithmic PDF of BonnorEbert spheres are shown in Fig. 6. The values were derived by solving the conditional equations for maxima positions as given in Appendix C.
The behavior of the PDF maxima of BonnorEbert spheres are compared with the maxima positions of the corresponding PDFs of spheres with analytical density profiles as given in Eq. (1) with n = 3 and n = 2. As found in Paper I, above a certain overpressure the local maxima of the PDF of spheres with analytical density profiles is related to a fixed parameter [ y_{n} ] _{max}, and consequently, considering Eq. (8), to a fixed normalized mass surface density [ X_{n} ] _{max}. The corresponding values for n = 2 and n = 3 are given in Table 2. The mass surface density [ Σ_{n} ] _{max} at the PDF maximum is obtained using the relation given in Eq. (11).
At low overpressures, the maximum coincides in general with the central mass surface density, so that the corresponding impact parameter is [ x ] _{max} = 0. A local maximum with [ x ] _{max}> 0 occurs
Location of the PDF maxima.
at a sufficiently high overpressure [ q^{1} ] _{0,max}. For analytical density profiles the lowest overpressure for a local maximum is given by (23)From Eq. (1) we find the corresponding unitfree cloud radius (24)A local maximum with x> 0 is otherwise related to a central mass surface density higher than (25)Above the lowest overpressure the normalized impact parameter at the PDF maximum of spheres with an analytical density profile behaves as (26)The impact parameter strongly moves outward with increasing overpressure and asymptotically approaches a constant value at high overpressure given by (27)We see from Fig. 6 that at low overpressure q^{1}< 100, the values of the normalized impact parameter and the mass surface density at PDF maximum of the BonnorEbert sphere follow the values of a sphere with n = 3. In this regime, a higher overpressure is related to a lower mass surface density at the PDF maximum. For the logarithmic PDF we find for the maximum (28)Assuming diffuse dust properties, this relates to a visual extinction (29)At higher overpressure, the impact parameter and the mass surface density at the maximum of the PDF of the BonnorEbert sphere fluctuate asymptotically to the corresponding constant value of a sphere with an analytical density profile with n = 2. If we consider again the logarithmic PDF, the asymptotic position of the PDF maximum is given by (30)or (31)Characteristic values of the maxima of the PDF are summarized in Table 2. We also add the corresponding values of critically stable conditions both for the BonnorEbert sphere () and the sphere with n = 3 where (see Sect. A.1.1). Isothermal selfgravitating pressurized cylinders are known to be stable under compression for all overpressures as dp(r_{cl})/dr_{cl}< 0 (see Fischera & Martin 2012b and references therein). The same would apply for a sphere with a hypothetical analytical profile with n = 2.
2.3.2. Asymptotic behavior for Σ ≪ [ Σ ] _{max}
The behavior of the PDF at low mass surface densities is related to the truncation of the density profile and the limit Σ ≪ [ Σ ] _{max} is equivalent to θ_{⊥} → θ_{cl}. In this regime the integrand in Eq. (6) for the mass surface density of BonnorEbert spheres becomes approximately constant with ρ/ρ_{c} ≈ e^{−ω(θcl)} = q and it follows that (32)Applying Eq. (17) and using r_{⊥} = θ_{⊥}/A, we obtain for the PDF (33)The replacement of the unitfree radius with provides (34)In the limit of infinite overpressure we obtain with n = 2 and ξ_{2} = 2 the asymptote (35)
2.3.3. Asymptotic behavior for Σ ≫ [ Σ ] _{max}
In Paper I we showed that in the regime of high mass surface densities the PDF of spheres with a truncated radial density profile as given Eq. (1) asymptotically approaches a power law (36)where (37)where B(a,b) is the betafunction. Using the relation Eq. (9) for the unitfree mass surface density X_{n} and the relation Eq. (21) provides for the approximations of BonnorEbert spheres the asymptotes (38)For the logarithmic PDF, the asymptote becomes ΣP(Σ) ∝ Σ^{−β} with β = 2/(n − 1). The slope agrees with what is expected for spheres with simple powerlaw profiles ρ ∝ r^{−n} of the local density (Kritsuk et al. 2011; Federrath & Klessen 2013).
We showed (Fig. 2) that the mass surface density at low impact parameters θ_{⊥} is dominated by the innermost part of the radial density profile of BonnorEbert spheres approximately described by the n = 3profile, which explains the slope of the asymptote (β = 1) close to the highest mass surface density. For BonnorEbert spheres with low overpressure we find using Eq. (38) (39)Figure 4 shows that for critically stable spheres the asymptotic behavior is not established. The slope of the logarithmic PDF is more consistent with β = 0.7 instead of 1. As Fig. 5 shows, Eq. (39) is not applicable for BonnorEbert spheres with high overpressure.
To approximate a BonnorEbert sphere in the limit of high overpressure (n = 2, ξ_{2} = 2, q ≈ 0) we obtain for the asymptote using Eq. (38) (40)The fluctuation of the PDF at high mass surface densities of BonnorEbert spheres with high overpressure (Fig. 5) are also seen in collapsingcloud models (Kritsuk et al. 2011). The slope in these simulations is somewhat steeper than for a simple BonnorEbert sphere (β ≈ 2.5), suggesting a flatter density profile (n ≈ 1.8, for negligible background).
2.4. Probability distribution function of selfgravitating cylinders
Fig. 7 PDF of the mass surface density of isothermal selfgravitating pressurized cylinders for various overpressures. The vertical line marks the pole position of the PDFs and is labeled with the corresponding pressure ratio q^{1}. The asymptotes at high and low mass surface density are shown as black dashed lines. The gray curves are the corresponding asymptotic PDFs (Eq. (43)) where the pole is removed. The powerlaw asymptotes at high mass surface densities are shown as gray dashed lines. 

Open with DEXTER 
The PDF of the mass surface density of cylinders with an analytical density profile as given in Eq. (1) is given by the implicit function (Paper I) (41)where (42)P_{sph}(X_{n})is the PDF of spheres of the unitfree mass surface density given by Eq. (22) and x is the normalized impact parameter given by x = r_{⊥}/r_{cl}.
The PDFs of the mass surface density of isothermal selfgravitating cylinders where n = 4 with various pressure ratios are shown in Fig. 7. The pressure ratios and the external pressure are the same as assumed for spheres in Fig. 5. Compared with the PDFs of BonnorEbert spheres, the PDFs of cylinders have a pole at the highest mass surface density. The PDFs of cylinders with high overpressure have a local maximum at low mass surface densities that continues to shift to lower mass surface densities for higher overpressures. For the same overpressure the central mass surface density is slightly lower than for spheres (see also Fischera 2011). For comparison, the asymptotic PDFs of cylindrical clouds are shown given by (Paper I) (43)Figure 7 shows that the asymptotic PDF provides the correct functional behavior of the cylindrical PDF apart from the pole, which is removed. For example, the asymptotic PDF provides the correct powerlaw asymptotes at low and high mass surface densities as well as the correct position of the local maximum of the PDF maximum for high overpressures.
2.4.1. Highest position
Because the pole of the PDF of cylinders are close to the highest mass surface density, the analysis of the highest positions of the PDF for isothermal selfgravitating pressurized cylinders have been based upon the asymptotic PDF. Qualitatively, the situation for the asymptotic PDF of cylinders is similar to the PDF of spheres. At low overpressure the local maximum coincides with the central mass surface density. Above a characteristic overpressure the asymptotic PDF develops a local maximum with [ x ] _{max}> 0, which corresponds to a fixed [ y_{n} ] _{max} and [ X_{n} ] _{max}. The values listed in Table 2 for cylinders with an analytical density profile with n = 4 are taken from Paper I. The characteristic parameters can be derived as shown for spheres with an analytical density profile.
The impact parameters and the mass surface densities at the local maximum of the asymptotic PDFs with [ x ] _{max}> 0 are shown in Fig. 8. The normalized impact parameter of the local maximum follows Eq. (26) for n = 4. It follows from Eq. (11) that the local maxima [ Σ_{4} ] _{max} decrease proportional to q^{1/4} with increasing overpressure. The dependence is slightly stronger than for BonnorEbert spheres with a low overpressure (q^{1}< 100). The local maxima of the logarithmic asymptotic PDF are located at (44)or (45)For low overpressures, the maxima of the asymptotic logarithmic PDFs agree within a factor of two with the maxima of BonnorEbert spheres. Figure 7 shows that at low overpressures the correct PDF of cylinders is still strongly dominated by the pole. In this case, the highest position of the asymptotic PDF approximately corresponds to the location of a shoulder in the correct PDF.
2.4.2. Asymptotic behavior for Σ ≪ [ Σ ] _{max}
In Paper I, we found for low mass surface densities the asymptote (46)For isothermal cylinders we obtain (47)
Fig. 8 Same as Fig. 6, but for the asymptotic PDF of isothermal selfgravitating cylinders where the PDF has no pole at the highest mass surface density. 

Open with DEXTER 
2.4.3. Asymptotic behavior for Σ ≫ [ Σ ] _{max}
For high overpressure the PDF of a cylinder with the analytical density profile with n> 1 behaves at high mass surface densites Σ_{n} ≫ [ Σ_{n} ] _{max} as (Eq. (B.7), Paper I) (48)For selfgravitating isothermal cylinders we obtain using the defintion Eq. (7) for X_{n} and the relation Eq. (41) for the different PDFs for X_{n} and Σ_{n} the asymptote (49)where (50)is the mass surface density through the center of a cylinder in the limit of high overpressure (Eq. (14)). For sufficiently high overpressure the PDF shows a powerlaw slope for [ Σ ] _{max} ≪ Σ ≪ Σ(0) given by (51)Figure 7 shows that this powerlaw behavior is only established for cylinders with overpressures greater than ~10^{3}.
3. Averaged probability distribution function
The origin of the observed PDF of the mass surface density of giant molecular clouds is very complex. In general, observations indicate that the global PDF of starforming molecular clouds is a combination of different components of randomly moving gas, generally attributed to turbulence, and the condensed structures. But structures related to windblown bubbles of massive stars or star clusters might also contribute to the global PDF.
The parameters determining the condensed spherical or elongated structures in a molecular cloud are not constant, but show certain variations. The main parameters for the pressurized structures of a certain geometry and density profile (power index n) are the pressure ratio q, the external pressure p_{ext}, and the temperature T_{cl}. The condensed clouds are located on a certain background Σ_{b} that might also vary for different clouds. The bounding pressure of the condensed structures is probably related to a pressure profile in the larger molecular cloud. Moreover, the condensed structures might form in a hierarchical fashion so that the bounding pressure of the innermost structure is also the central pressure of the surrounding cloud.
Observations (Schneider et al. 2012, 2013) suggest that as a first approach the condensations can be treated as individual structures within the larger molecular cloud. In this case, we can obtain the mean PDF for the cloud through an average of the PDFs of individual condensations. We furthermore assume that the different parameters determining the individual structures are independent (uncorrelated) variables.
We refer to the PDF of a single condensation under the condition of parameters a as P(Σ,a), where a = (q,p_{ext},Σ_{b}). Under these assumptions, the averaged probability distribution over the distribution functions P(a) is then given by (52)where κ = 2 for spheres and κ = 1 for cylinders.
3.1. Temperature distribution
The mass of a core for a given external pressure and overpressure varies strongly with the kinetic temperature (M_{sph} ∝ T^{2}) (see Appendix A.1). But as we showed in the previous section, the temperature has no effect on the PDF of the column density if both the external pressure and the overpressure are fixed. All critically stable cores in a certain pressure region, for example, have the same PDF of their mass surface density independent of their temperature. However, because the total size scales with temperature (see Appendix A.2), the temperature distribution affects the covering factor and therefore the ratio between the PDFs of the cores and the PDF of the surrounding gas.
3.2. Distribution of the gravitational states
In general, the condensed cores in giant molecular clouds are probably in various different gravitational states that are characterized by different pressure ratios q. The situation of the condensed cores might be similar to the situation of Bok globules, which, according to the study of Kandori et al. (2005), show a variety of qvalues.
To demonstrate the effect of a distribution of the gravitational state on the average PDF of the column density we consider a large sample of cores that are pressurized by the same external pressure, which is taken to be p_{ext}/k = 2 × 10^{4} K cm^{3}. In addition, we consider the cores to have the same effective temperature and the same K. As shown in Appendix A.1, the masses of supercritical cores, for example, are smaller than the mass of the critically stable core by up to a factor 2. By fixing p_{ext} and K, the mass of the cores does not vary monotonically with overpressure.
It seems reasonable to assume that the distribution of gravitational states of subcritical and stable spheres () follows functional forms different from the distribution of clouds in a supercritical state. As indicated in the study of Bok globules by Kandori et al. (2005), most cores are probably close to the critical value. This might be related to the fact that clouds with a deeper gravitational potential become more stable with a deeper gravitational potential against external disturbances. It might also be a simple observational bias because the central extinction through BonnorEbert spheres strongly increases when the mass approaches the critical value (Fischera & Dopita 2008). The distribution of supercritical states is most likely related to gravitational collapse, as was also considered in the study by Kandori et al. (2005). Because a detailed analysis goes beyond the scope of this paper, we simply assumed that the probability of supercritical states decreases with increasing central density or overpressure q^{1}.
To be able to consider two different distributions for spheres with low and high overpressure we assume for the probability function of the overpressures the function (53)where k_{2}> 1. The constant C can be neglected in this study because of the additional renormalization of Eq. (52). The parameter is a characteristic overpressure determining the transition between the two regimes. The distribution at q ≫ q_{0} and q ≪ q_{0} are power laws given by P(q^{1}) ~ Cq^{−k1} and . The parameter γ determines the smoothness of the transition between the two power laws. In the limit γ → ∞ the distribution, for example, becomes a broken power law. In this study we assume γ = 2. For the distribution we assume a power index k_{1} = 2, which produces a decreasing probability at lower gravitational states (smaller q^{1}) for .
The distribution does not vanish for q^{1} → 1, which seems unrealistic. However, the product of the surface area with the probability distribution of the overpressure does indeed vanish as r_{sph}(q) ~ r_{n}(q) → 0 when q → 1 according to Eq. (A.8).
Fig. 9 Mean PDF of the mass surface density for an ensemble of cores assumed to be BonnorEbert spheres with a distribution of gravitational states or pressure ratios q^{1} as given in Eq. (53). For comparison the mean PDF is also derived for spheres with an analytical density profile as given in Eq. (1) with n = 3. The characteristic pressure ratio of the probability distribution of the gravitational states is either q_{0} = q_{crit} ~ 1/14.04 of a critically stable sphere (top panel) or q_{0} ≫ 1 (bottom panel). The various curves labeled a), b), and c) correspond to different power indices k_{2} assumed to be 2.0, 1.5, and 1.2. The power k_{1} = 2 and the smoothness parameter γ = 2 are the same in all calculations. The black dashed lines are powerlaw approximations of the PDFs in the limit of low and high mass surface densities for an analytical density profile with n = 3 (Appendix E.1). The external pressure is assumed to be p_{ext}/k = 2 × 10^{4} K cm^{3}. Also shown are the PDFs of a critically stable sphere and a sphere with infinite overpressure (gray lines and gray annotation). The corresponding powerlaw asymptotes of the PDFs of single cores in the limit of low (Eq. (34)) and high (Eq. (39) and Eq. (40)) mass surface densities are shown as gray dashed lines. 

Open with DEXTER 
To show the effect of the characteristic overpressure we considered two different cases. In a first set of calculations the characteristic overpressure is taken to be the overpressure of a critically stable sphere (). In a second set of calculations the power law of the high q^{1}states extends down to lowest gravitational states by choosing . For a direct comparison, the same assumptions were made for both spheres and cylinders. The power index k_{2} was varied for given geometry and chosen characteristic overpressure to visualize its dependence on the mean PDF. For spheres the assumed power indices k_{2} were 1.2, 1.5, and 2.0 and for cylinders 1.5, 2.0, and 3.0.
3.2.1. Mean PDF of spheres
The mean PDFs for distributions of BonnorEbert spheres are shown in Fig. 9. They are compared with the mean PDFs of spheres with analytical density profiles as given in Eq. (1) with n = 3.
Fig. 10 Mean PDF of the column density for an ensemble of isothermal selfgravitating cylinders with a variety of gravitational states or overpressures q^{1}, but fixed inclination angle i. The probability distribution of the overpressures is given by Eq. (53) with fixed q_{0}, k_{1} = 2, and γ = 2 where q_{0} is either 1/14.04 (lefthand figure) or q_{0} ≫ 1 (righthand figure). The power index k_{2} is assumed to be 1.5, 2, and 3. The cylinders are pressurized by a medium with p_{ext}/k = 2 × 10^{4} K cm^{3}. The black dashed lines are asymptotes at low and high mass surface densities (Appendix E.2). For comparison, the PDF of a single cylinder with overpressure q^{1} = 14.04 is shown (gray line). The gray dotted line is the asymptotic PDF as defined in Eq. (43), and the gray dashed lines are the corresponding asymptotes at low (Eq. (47)) and high (Eq. (51)) mass surface densities. 

Open with DEXTER 
The mean PDF is characterized by a broad peak at mass surface densities between the highest position of a critically stable sphere and a sphere with infinite overpressure. For the reference pressure the maxima lie at A_{V} ~ (1−2) mag. At low and high mass surface densities the mean PDF asymptotically approaches the power laws discussed quantitatively in Appendix E for analytical density profiles.
The power law at low mass surface densities is identical to the powerlaw asymptote at low mass surface densities of pressurized spheres (Sect. 2.3.2) and cylinders (Sect. 2.4.2) with . The absolute probabilities follow the probabilities of single highly supercritical BonnorEbert spheres within a factor of two.
The power law at high mass surface densities is independent of the power n of the radial density profile and is given by . The slope steepens as one would expect for larger k_{2} and lower probabilities for spheres with high overpressures. The values of k_{2} are limited to k_{2}> 1. Figure 9 shows that for k_{2} → 1 the mean PDF approaches the PDF of BonnorEbert spheres with infinite overpressure.
The mean PDF for a characteristic pressure ratio q_{0} = 1/14.04 deviates considerably from the PDF of a critically stable BonnorEbert sphere. At the central mass surface density of a critically stable sphere the mean PDF shows a knee with a flatter curvature at lower and steeper curvature at higher mass surface densities. The feature is more pronounced for higher k_{2}. The mean PDF for q_{0} ≫ 1 has no additional feature apart from the prominent peak.
As we see in the figure, for high k_{2} and the assumed characteristic pressures q_{0} the mean PDF of a distribution of BonnorEbert spheres is well represented by the corresponding mean PDF of a sphere with the analytical density profile with n = 3. The two curves only start to deviate for low k_{2} where spheres with higher overpressure than q^{1} ~ 100 attribute to the mean PDF where the analytical density profile with n = 3 is not a valid approximation of a BonnorEbert sphere. We have seen that the PDF of an individual sphere with an analytical density profile with q^{1}> 100 is higher than the PDF of a BonnorEbert sphere (Fig. 5) at both low and high mass surface densities. In this case, the mean PDF of the analytical density profile is therefore too high with respect to the correct mean PDF, as seen for k_{2} = 1.2 and q_{0} = q_{crit}.
3.2.2. Mean PDF of cylinders
The mean PDF of an ensemble of selfgravitating isothermal cylinders pressurized by the same external pressure, but with a distribution of pressure ratios q, is shown in Fig. 10. The mean PDF is compared with the PDF of a single cylinder with an overpressure q^{1} = 14.04.
One of the main characteristics is that the mean PDF of a distribution of cylinders has no pole. If we consider the mean PDF for the PDF has a broad feature located approximately at the pole position of the PDF of a single cylinder with a pressure ratio q = q_{0}. The feature broadens and its maximum shifts to higher mass surface densities for a flatter distribution of q^{1} or equivalently lower k_{2}. The mean PDF for q_{0} ≫ 1, as is also the case for BonnorEbert spheres, has no additional feature apart from the prominent peak at A_{V} ~ 1 mag.
At low and high mass surface densities the mean PDF has powerlaw asymptotes, similar to the mean PDF of spheres. The asymptotes shown in Fig. 10 are derived in Appendix E.2. At low mass surface densities the asymptote is again given by . At high mass surface densities the slope is flatter than that of a distribution of spheres with .
Fig. 11 Angleaveraged PDF of the column density for cylinders with a mean inclination angle of 0° (lefthand figure) and 60° (righthand figure). The cosine of the inclination angle is assumed to be Gaussian distributed. The curves are labeled with the assumed standard deviation σ_{μ}. The cylinders are pressurized by a medium with p_{ext}/k = 2 × 10^{4} K cm^{3} and have an overpressure p_{c}/p_{ext} = 14.04. The intrinsic PDF for the limit σ_{μ} → 0 is shown as dotted curves. The gray vertical lines mark the highest extinction values of a cylinder seen under the assumed inclination angle and the position of the peak of the averaged PDF of randomly oriented cylinders. The thick curves are the PDFs for randomly distributed cylinders. 

Open with DEXTER 
3.3. Angleaveraged PDF for cylinders
For filaments, we have the additional complication that the orientation most probably varies not only between different filaments, but may also vary along individual filamentary structures. To analyze its effect on the mean PDF we consider an ensemble of cylinders with a certain distribution of inclination angles i that are otherwise identical (same q).
For randomly distributed filaments, the probability of the cosine of the projection angle i is a constant. However, observations reveal that the orientation of dense massive structures, which probably dominate the PDF, is not completely random. The molecular cloud IC 5146 (Arzoumanian et al. 2011) or the Taurus complex (Kainulainen et al. 2009) are dominated by massive filamentary structures with a certain mean orientation. We ignore the complication that the orientation can change at 90° angles, as is the case for the Taurus filament. To show the effect of any distribution of orientation angles, we assume that the cosine μ = cosi of the inclination angle i has a Gaussian distribution P(μ) around a certain mean μ_{0} with a standard deviation σ_{μ}. The distribution is renormalized so that (54)The limit σ_{μ} ≫ 1 produces a flat distribution of randomly orientated cylinders with P(μ) = 1.
If P_{cyl}(Σ) is the PDF for cylinders seen edgeon, the PDF of a cylinder seen at an inclination angle i is given by (55)The average over all inclination angles is then given by (56)where (57)and where μ_{min} is related to the greatest length of the filament. For infinitely long filaments we have μ_{min} = 0. The mean PDF for Σ ≤ Σ(0) is an average over all inclination angles. For Σ > Σ(0) the mean PDF is related to cylinders seen at increasingly high inclination angles.
For randomly distributed cylinders we derive(58)For Σ′> Σ(0) the angleaveraged PDF of randomly distributed cylinders becomes a simple power law where (59)This powerlaw slope is the same as for spheres with a density profile with n = 3 in the limit of high overpressure (Figs. 4 and 5).
The effect of the angular distribution on the PDF for cylinders is shown in Fig. 11. A distribution of inclination angles reduces the peak at the highest column density and produces a tail beyond the peak that flattens for wider distributions up to P_{cyl}(Σ) ∝ Σ^{2} for randomly oriented cylinders. The effect on the peak is considerably stronger for cylinders with a mean inclination angle i> 0 (μ_{0}< 1). Figure 11 shows that for μ_{0}< 1 the angular distribution also broadens the feature and produces a shift to lower mass surface densities. The PDF of randomly distributed cylinders has a small feature at the central mass surface density for i = 0. Because the filamentary structure in the ISM is unlikely to have a welldefined overpressure (p_{c}/p_{ext}), this feature might be much less prominent than shown in Fig. 11, or might even be absent.
3.4. Backgroundaveraged PDF
In the previous sections we have considered the PDFs of the mass surface density of condensed structures for a negligible background. However, this is certainly no longer valid if the condensed structures are located in molecular clouds or if they are seen through the interstellar medium. As an example, we consider randomly distributed cores in an otherwise turbulent medium. As we show below, the combined PDF of the turbulent structure and the condensed structures resemble the main features of the observed PDF of typical starforming clouds.
3.4.1. Turbulent medium
As shown in hydrodynamical simulation of driven turbulence, turbulence produces a lognormal density distribution of the local density (VazquezSemadeni 1994; Padoan et al. 1997a; Passot & VázquezSemadeni 1998).
The lognormal function of a statistical variable s with a mean value ⟨s⟩ is given by (60)where ŝ = s/ ⟨s⟩ is the normalized variable and σ_{lnŝ} the standard deviation of the lognormal function, which is related to the standard deviation of the variable s through (61)According to simulations of driven turbulence (Padoan et al. 1997a; Nordlund & Padoan 1999), there exists a simple linear relation between the standard deviation of the local density and the Mach number M, where (62)The proportionality factor b and the applicability for the interstellar medium is still a debated question, as discussed in Paper III (Fischera 2014b) or by Kainulainen & Tan (2013). Observational studies of IC 5146 and the Taurus molecular cloud complex suggest a value of b ~ 0.5 (Padoan et al. 1997a; Brunt 2010). A lower value (b ~ 0.20) has been proposed by Kainulainen & Tan (2013) based on high dynamicrange extinction mapping of infrared dark clouds. Simulations indicate a value between 0.26 and 1 (Padoan et al. 1997b; Passot & VázquezSemadeni 1998; Kritsuk et al. 2007; Beetz et al. 2008; Federrath et al. 2008, 2010; Price et al. 2011).
The PDF of the column density through simulated turbulence was also found to be close to a lognormal function (Ostriker et al. 2001; VázquezSemadeni & García 2001; Brunt et al. 2010a). The results are supported by studies of the column density through nonstarforming molecular clouds that agree well with a simple lognormal function, as in the case of Lupus V and the Coalsack (Kainulainen et al. 2009).
The statistical properties of the local density are directly related to the properties of the column density. Their general functional dependence has been described by Fischera & Dopita (2004) and Brunt et al. (2010b,a). Fischera & Dopita (2004) have studied the relationship between the column density and the local density of an isothermal turbulent screen by examining the one and twopoint statistics. They assumed a simple lognormal function of the density distribution of the local density and a simple power law of the power spectrum, where is the Fourier coefficient of the local density ρ and k is the wavenumber.
It has been found that if the screen thickness Δ is greater than the highest turbulent scale L_{max}, the variance of the column density of thick screens varies as (Fischera & Dopita 2004) (63)where m ≠ −2 and m ≠ −3. For isolated molecular clouds the thickness should be Δ ≥ L_{max}. Relation 63 allows an estimate of the variance of the local density by measuring both the variance of the column density and the power m of the power spectrum. For a Kolmogorovlike power spectrum of the local density^{2}, where m = −11/3, and for a wide range of turbulent length scales (k_{max} ≫ k_{min}), we have, for example^{3}, (64)
3.4.2. Combined PDF of cores and turbulent gas
The cores are considered to be small relative to the molecular cloud so that the variations of the column density of the background of individual cores can be neglected. In addition, we consider the simplified case where the cores do not overlap and are homogeneously distributed within the turbulent molecular cloud. The backgroundaveraged PDF of the mass surface density of the cores is then given by the convolution (65)where P_{cl}(Σ) is the PDF of a single core and P_{turb}(Σ) the PDF of the turbulent medium given by the lognormal function discussed in the previous section. For a narrow distribution of the mass surface density caused by the turbulent medium, the PDF of the background becomes a delta function, and the mean is simply given by (66)The cores are assumed to cover a certain fraction f_{A} of the area of the cloud. The total PDF of the molecular cloud is then given by the combined PDF of the turbulent background and the backgroundaveraged PDF of the cores, so that (67)Examples of the combined PDF of cores and the turbulent medium are shown in Fig. 12 for two different gravitational states of the cores and three different assumptions of the turbulent medium. The cores are either critically stable or highly supercritical BonnorEbert spheres to mimic collapsing cores. The density profiles are approximated using the corresponding analytical density profiles for low (n = 3) and high (n = 2) overpressure. The pressure within the turbulent medium is again assumed to be p_{ext}/k = 2 × 10^{4} K cm^{3}. The mean background of the turbulent gas is related to the central mass surface density of a critically stable BonnorEbert sphere. The ratios ⟨Σ_{b}⟩/Σ(0) for the mean backgrounds ⟨Σ_{b}⟩ and central mass surface density Σ(0) through a critically stable sphere are 10^{1.0}, 10^{0.5}, and 10^{0.0} to cover the range from a low to a very high background level. The lowest background is close to the peak position for a number of PDFs of molecular clouds studied by Kainulainen et al. (2009).
Fig. 12 Combined PDF of a molecular cloud consisting of a turbulent medium (dotted curves) and uniformly distributed condensed cores (dasheddotted curves). The cores are either critically stable (left) or in a state of collapse (right) where the density profiles are described using the corresponding approximations for BonnorEbert spheres. The turbulent medium is modeled through a simple lognormal function and the pressure in the turbulent gas is assumed to be p_{ext}/k = 2 × 10^{4} K cm^{3}. The spheres do not overlap and cover f_{A} = 40% of the area of the molecular cloud. The distributions are labeled with the assumed ratio ⟨Σ_{b}⟩/Σ(0) of the mean background ⟨Σ_{b}⟩ and the central mass surface density Σ(0) through a critically stable BonnorEbert sphere. The scaled intrinsic nonbackgroundaveraged PDFs of BonnorEbert spheres (see Fig. 4) are shown as solid gray lines. Also given are the corresponding powerlaw asymptotes at low (Eq. (34)) and at high (Eq. (39) and Eq. (40)) mass surface densities. 

Open with DEXTER 
To assign a standard deviation of the mass surface density to the different mean backgrounds we considered an idealized turbulent slab with a clearly defined mean density ⟨ρ⟩ and a fixed greatest turbulent length scale L_{max} so that . With Eq. (63), a higher mean background is then related to a lower fluctuation of the normalized mass surface densities Σ/⟨Σ⟩. As a typical standard deviation we assumed for the lowest mean background σ_{Σ/⟨Σ⟩} = 0.40 (Table 1, Kainulainen et al. 2009), which provides a standard deviation of 0.22 and 0.13 for the next higher backgrounds. This implies, considering Eq. (64), a standard deviation of the local density of . Assuming for the correlation constant b = 0.5 between σ_{ρ/ ⟨ρ⟩} and M, this suggests a Mach number of . According to Kainulainen & Tan (2013), molecular clouds typically have a Mach number of about M ~ 10, which would imply Δ ≫ L_{max}, or, considering a thickness Δ ≈ L_{max}, a smaller correlation coefficient or possibly a flatter power spectrum.
The combined PDF resembles several features of the global PDF of starforming molecular clouds. The PDF is characterized by a broad peak and a tail at high mass surface densities. The mean PDF of critically stable cores embedded in the turbulent medium is furthermore characterized by a strong decrease of probabilities at high mass surface densities.
The background clearly modifies the shape of the PDF of individual cores. The backgroundaveraged PDF appears squeezed compared with the intrinsic PDF, where the probabilities are shifted to higher mass surface densities. This is visible, for example, in the shift of the highest position or in the shift of the knee in the case of critically stable cores. Naturally, the effect on the PDF of the cores increases with background level and is strongest at low mass surface densities.
As discussed in Sect. 2.3.3 and shown in Fig. 12, the high end tail of the PDF of individual cores without background only becomes a power law in the limit of high mass surface densities. However, a low background can reduce the difference between the PDF of the cores and the asymptote and may even produce a tail at high mass surface densities below a possible knee that approximately is a power law. On the other hand, for higher backgrounds with ⟨Σ_{b}⟩/Σ(0) ≫ 0.1, the tail at high mass surface densities develops a curved shape where the probabilities at the low end are considerably higher than the powerlaw asymptote.
As a result of the additive nature of two different components, the functional form of the mean PDF is not necessarily a simple composition of a turbulent part at low mass surface densities and a tail at high mass surface densities. The appearance will depend on the covering factor and the background level for the condensed cores. Indeed, the theoretical curves of the mean PDFs derived for the high backgrounds show a smooth transition of the lognormal part and a curved tail. For a low background, however, the tail at high mass surface densities shows a broad feature at the low end that is caused by the shifted peak of the backgroundaveraged PDF.
4. Discussion
We have analyzed the PDF of the projected density for several different aspects related to isothermal and selfgravitating structures that are pressurized by the ambient gas. In the following we qualitatively compare the theoretical predictions with the basic properties of observed PDFs of star and nonstarforming molecular cloud as derived by Kainulainen et al. (2009). Paper IV (Fischera 2014c) will present a quantitative analysis of the PDF at low mass surface densities, including the broad maximum referred to as the turbulent part and the tail at high mass surface densities.
4.1. Highend tail of the global PDF
The tail at high mass surface densities often has a powerlaw form. For the PDFs with a maximum at A_{V} ~ 1 mag, the index β of the power law approximation Σ^{−β} of the tail of the logarithmic PDF lies for most cases in between 1 and 2. An unusually steep slope with β> 4 can be found for the global PDF of the Pipe molecular cloud, which is also characterized by an unusually high maximum position around A_{V} ≈ 3 mag.
We showed that a distribution of inclination angles of not overlapping cylinders would produce a steep slope at high mass surface densities. The flattest tail with β = 1 is expected for randomly distributed cylinders. In principle, the slope would be steeper if the filaments had a certain alignment, but this would also be accompanied by a strong feature that is not detected in the observed PDFs. In this case, additional effects need to be considered that might reduce or even remove the feature as a distribution of gravitational states, as shown in this study, or distortions along filaments possibly associated to turbulence, which will be analyzed in Paper III. Still, the tail is expected to show features caused by overlapping filaments that are not accounted for.
We showed that a powerlaw distribution P(q^{1}) ∝ 1/(q^{1})^{k2} of the gravitational states of BonnorEbert spheres or cylinders would lead to a powerlaw asymptote in the limit of high mass surface densities with β = 2k_{2} for spheres or β = 2k_{2} − 1 for cylinders. A steeper distribution of q^{1} would be related to a steeper powerlaw asymptote.
The simplest model for the tail at high mass surface densities are single BonnorEbert spheres. We showed that spheres with low overpressure <~100 have a rather flat asymptote at high mass surface densities with β = 1 and that the slope above the PDF maximum steepens with overpressure. Clouds with regions essentially in gravitational collapse and in the process of forming new stars are expected to have a tail with a slope close to β = 2, as observed for the Ophiuchus cloud or Taurus. A flatter tail would be expected for clouds without any star formation activity. In this case, the tail is expected to have a knee or a break at high mass surface densities. These characteristics can be found in the tail of the PDF derived for the Musca molecular cloud, which is classified as nonstarforming. However, the agreement may be coincidental, considering its elongated structure.
A break can also be observed in some other PDFs such as Lupus III, Lupus I, LDN 1228, LDN 204, Ori A, the Perseus cloud, or Ori B. It occurs around A_{V} = 6−15 mag, which approximately coincides with the central extinction through critically stable spheres pressurized by a medium with p_{ext}/k ~ 2 × 10^{4} K cm^{3}. However, it cannot be excluded that the break may also be caused by the resolution of the map where the high mass surface densities in the central regions of collapsing clouds maybe smoothed over a larger area. These observational problems need to be addressed and can be solved using maps with higher resolution.
The shape of the tail suggests that the condensed structures are seen against a certain background and that the slope of the tail at high mass surface densities steepens with the highest position of the global PDF, as expected for cores that are randomly distributed within a turbulent medium (Sect. 3.4).
We showed that the tail would appear in the case of a low background close to the power law, which is otherwise only expected as an asymptote in the limit of high mass surface densities. For the PDF of the Taurus complex the power law continues almost down to the highest position without a clear separation of the two components.
The rather steep slope of the PDF tail (β> 4) of the Pipe molecular cloud might be caused by a rather high background, as suggested by the high maximum position at A_{V} ~ 3 mag. The same may apply to the global PDF of the Rosette molecular cloud. The logarithmic PDF derived by Schneider et al. (2012, 2013) using Herschel observations does not only have a relatively high maximum at A_{V} = 2 mag, but also an unusually steep tail with β ~ 3.1 (Fig. 6, Schneider et al. 2012). Relating this slope to a powerlaw density profile ρ(r) ∝ r^{−n} using the relation n = (2 /β) + 1 (Kritsuk et al. 2011; Federrath & Klessen 2013; Fischera 2014a) would indicate a rather flat density profile with power n = 1.65^{4}. However, the steeper slope may equally well be a direct consequence of condensed structures seen against a relatively high background and not directly related to the density profile.
The PDF of Ophiuchus has a broad peak or a shoulder at the low end of the tail around A_{V} ~ 3 mag, similar to the broad feature in the theoretical mean PDF of condensed cores located on a low background (Fig. 12). The feature at the low end of the PDF tail of the Pipe molecular cloud might have a similar explanation.
4.2. Turbulent part of the global PDF
Most of the derived PDFs of molecular clouds studied by Kainulainen et al. (2009) have a maximum around A_{V} = 1 mag, which is close to the maximum of the PDFs of BonnorEbert spheres pressurized by the assumed ISM pressure. Although it may also be coincidental, this would suggest that the molecular gas as a whole is not only pressurized, but also selfgravitating.
At mass surface densities below the maximum, the observed PDFs show an obvious deviation to a simple lognormal function. For a number of cases the PDF at low mass surface densities has a functional form close to a power law. The deviation may be related to an unsymmetrical nature of the local density PDF, as seen in more current simulations of forced turbulence (Federrath et al. 2008; Konstandin et al. 2012; Federrath & Klessen 2013). The behavior in the limit of low mass surface density might also be related to a truncated density profile of the larger cloud if it is pressurized by a low dense medium. We showed that in the case of a negligible background, the PDF of pressurized nonturbulent spheres and cylinders asymptotically approaches a power law ΣP(Σ) ∝ Σ^{2} at low mass surface densities. The actual shape may be flatter because of the turbulent nature of the gas, even if the PDF of the local density is lognormal, as will be shown in Paper III. It still needs to be considered that the shape in particular at low mass surface densities will be affected by additional material along the line of sight that is not related to the molecular cloud, which will shift the probabilities to higher mass surface densities.
A first estimate of the density fluctuation in the turbulent cloud surrounding the condensed structures can be derived using Eq. (63), which is based on the infinite turbulent slab approximation. The estimate of the density contrast would be uncertain by the unknown power spectrum of the density structure and the number of turbulent length scale. A statistical analysis of molecular clouds should therefore be combined with studies of the power spectrum of the projected density. A better estimate of the density contrast might be obtained by taking geometric effects and a possible background into account (Papers III and IV).
It appears to be reasonable to assume that most of the starforming gas lies in the central region of the cloud where the gas is shielded from the interstellar radiation field to allow the gas to cool. The gas also experiences the highest pressure in the center, so that in the case of dynamic equilibrium the gas becomes most easily gravitationally unstable because the critical mass varies as . This would suggest a threshold for star formation as indicated by observational studies (Kirk et al. 2006; Foster et al. 2009). The model described in this paper should allow an independent determination of the central pressure in the turbulent gas and might therefore provide a better understanding of the star formation process.
5. Summary and conclusion
We have analyzed the properties of the PDF of the mass surface density of pressurized isothermal selfgravitating spheres, known as BonnorEbert spheres, and cylinders where we applied the results obtained in Paper I for a simple analytical density profile as a generalization of physical density profiles given by ρ(r) = ρ_{c}/ (1 + (r/r_{0})^{2})^{n/ 2}. We showed that for ρ(r) /ρ_{c}> 0.01 the density profile of BonnorEbert spheres is well approximated by the analytical profile with power n = 3. At larger cloud radius the radial density profile fluctuates asymptotically to the analytical profile with n = 2. Emphasis was given on the dependence of the PDF of given geometry on the pressure ratio p_{ext}/p_{c} and the external pressure assumed to be p_{ext}/k = 2 × 10^{4} K cm^{3}.
The main properties of the PDFs of individual spheres and cylinders are found to be as follows:

1.
The PDF of critically stable spheres is truncated at the centralmass surface density or extinction value .

2.
Below an overpressure of ~100, the PDF of a BonnorEbert sphere can be well approximated by the PDF of a sphere with an analytical density profile with n = 3. At fixed mass surface density, the PDF of spheres fluctuates asymptotically with overpressure to the value of a sphere with an analytical sphere with n = 2 in the limit of infinite overpressure.

3.
For the reference pressure the maximum of the logarithmic PDF of a sphere lies between A_{V} ~ 1 mag and 3 mag. For spheres with an overpressure <~100, the highest position of the PDF shifts proportionally to q^{1/6} to lower mass surface densities. For highly supercritical clouds the highest position of the logarithmic PDF asymptotically approaches a constant given by mag.

4.
The flattening of the radial density profile of spheres with overpressure is related to a steepening of the PDF at mass surface densities above the peak maximum. The slope β of the powerlaw asymptote ΣP_{sph}(Σ) ∝ Σ^{−β} varies from β = 1 for spheres with overpressures below ~100 to β = 2 for highly supercritical spheres.

5.
The highest position of the logarithmic asymptotic PDF (underlying PDF without a pole) of an isothermal selfgravitating pressurized cylinder decreases with q^{1/4}. For overpressures below 100, the maxima of cylinders and spheres agree within a factor of two.

6.
For cylinders with high overpressure, the PDF at high mass surface densities is approximately given by P_{cyl}(Σ) ∝ Σ^{−4/3}(1 − (Σ/Σ(0))^{2/3})^{−1/2}.

7.
At low mass surface densities the PDF of both spheres and cylinders approaches a power law P_{cl}(Σ) ∝ Σ.
In addition to individual clouds, an ensemble of spheres or cylinders with a distribution of overpressures q^{1} = p_{c}/p_{ext} given by were considered. Studied were distributions with fixed parameters k_{1} = 2 and γ = 2 and different values for q_{0} and k_{2}. The corresponding mean PDF has the following properties:

1.
The distribution does not change the asymptotic behavior of thePDF at low mass surface density, which is, as in case of individualclouds, proportional to a power law ⟨P(Σ)⟩ ∝ Σ.

2.
At high mass surface densities, the PDF asymptotically approaches a power law that is independent of the radial density profile. This asymptote is either for spheres or for cylinders.

3.
For cylinders, the distribution effectively decreases the high probabilities at the highest mass surface density of single cylinders.

4.
For q_{0} ≫ 1 the mean PDF maximum of spheres and cylinders lies at A_{V} = 1−2 mag.
The mean PDF of an ensemble of cylinders with the same overpressure, but a distribution of inclination angles, decreases the high probabilities at the highest mass surface density of single cylinders. Randomly oriented and not overlapping cylinders with a narrow distribution of gravitational states will produce a PDF with a powerlaw tail ΣP(Σ) ∝ Σ^{1} at high mass surface densities.
A simple model of the global PDF of molecular clouds was presented based on a combination of a turbulent medium and embedded randomly distributed BonnorEbert spheres. The model apparently reproduces the basic features of many of the observed PDFs derived by Kainulainen et al. (2009).

1.
The combined model produces a lognormal function aroundthe peak and a tail at high mass surface densities where the relativeratio of the two components is related to the covering factor of thecondensed cores.

2.
At low covering factor, the combined PDF shows a break between the two components.

3.
The functional form of the tail at high mass surface densities is affected by the background level (PDF maximum of the turbulent medium) relative to the PDF maximum of the cores. a) A low background can produce a tail that across the whole range is approximately a power law. b) A high background will lead to a curved tail, where the probabilities at the low end are higher than a power law.
For the assumed reference pressure the background can be considered to be low for A_{V}< ~ 1 mag and high for A_{V}> ~ 1 mag. The curvature for a high background might not be easily detectable because of limiting resolution and noise in the observed data, so that the assumption of a simple power law would imply a steeper slope than expected for the radial density profile.
On the other hand, there are claims that starforming filaments in particular have mass ratios partly considerably higher than the theoretical highest value (Arzoumanian et al. 2011). However, this might be related to the low temperature considered in observational studies (Fischera & Martin 2012b).
This assumption is correct for a contaminant driven by a Kolmogorov velocity field as might apply to dust particles (Lazarian & Pogosyan 2000; Fischera & Dopita 2004).
Fischera & Dopita (2004) took the power index of the Kolmogorov velocity field to be m = −10/3, which leads to .
The even flatter power law given in the study of Schneider et al. (2012) is inaccurate and has been corrected (Schneider et al. 2013).
Acknowledgments
This work was supported by grants from the Natural Sciences and Engineering Research Council of Canada and the Canadian Space Agency. I would like to thank Richard Tuffs for his support and Jouni Kainulainen for helpful discussions. Furthermore, I would like to thank the referee for a careful reading and the suggestions that helped to improve the manuscript. Many thanks go to my parents who enabled me to finish the project.
References
 Arzoumanian, D., André, P., Didelon, P., et al. 2011, A&A, 529, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Beetz, C., Schwarz, C., Dreher, J., & Grauer, R. 2008, Phys. Lett. A, 372, 3037 [NASA ADS] [CrossRef] [Google Scholar]
 Bohlin, R. C., Savage, B. D., & Drake, J. F. 1978, ApJ, 224, 132 [NASA ADS] [CrossRef] [Google Scholar]
 Bonnor, W. B. 1956, MNRAS, 116, 351 [NASA ADS] [CrossRef] [Google Scholar]
 Boulares, A., & Cox, D. P. 1990, ApJ, 365, 544 [NASA ADS] [CrossRef] [Google Scholar]
 Brunt, C. M. 2010, A&A, 513, A67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Brunt, C. M., Federrath, C., & Price, D. J. 2010a, MNRAS, 405, L56 [NASA ADS] [CrossRef] [Google Scholar]
 Brunt, C. M., Federrath, C., & Price, D. J. 2010b, MNRAS, 403, 1507 [NASA ADS] [CrossRef] [Google Scholar]
 Curry, C. L., & McKee, C. F. 2000, ApJ, 528, 734 [NASA ADS] [CrossRef] [Google Scholar]
 Ebert, R. 1955, Zeitschrift fur Astrophysik, 37, 217 [NASA ADS] [Google Scholar]
 Federrath, C., & Klessen, R. S. 2013, ApJ, 763, 51 [NASA ADS] [CrossRef] [Google Scholar]
 Federrath, C., Klessen, R. S., & Schmidt, W. 2008, ApJ, 688, L79 [NASA ADS] [CrossRef] [Google Scholar]
 Federrath, C., RomanDuval, J., Klessen, R. S., Schmidt, W., & Mac Low, M.M. 2010, A&A, 512, A81 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Federrath, C., Sur, S., Schleicher, D. R. G., Banerjee, R., & Klessen, R. S. 2011, ApJ, 731, 62 [NASA ADS] [CrossRef] [Google Scholar]
 Fiege, J. D., & Pudritz, R. E. 2000, MNRAS, 311, 85 [NASA ADS] [CrossRef] [Google Scholar]
 Fischera, J. 2011, A&A, 526, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fischera, J. 2014a, A&A, 565, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fischera, J. 2014b, A&A, submitted [Google Scholar]
 Fischera, J. 2014c, A&A, submitted [Google Scholar]
 Fischera, J., & Dopita, M. 2004, ApJ, 611, 911 [NASA ADS] [CrossRef] [Google Scholar]
 Fischera, J., & Dopita, M. 2008, ApJS, 176, 164 [NASA ADS] [CrossRef] [Google Scholar]
 Fischera, J., & Martin, P. G. 2012a, A&A, 547, A86 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fischera, J., & Martin, P. G. 2012b, A&A, 542, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fitzpatrick, E. L. 1999, PASP, 111, 63 [NASA ADS] [CrossRef] [Google Scholar]
 Foster, J. B., Rosolowsky, E. W., Kauffmann, J., et al. 2009, ApJ, 696, 298 [NASA ADS] [CrossRef] [Google Scholar]
 Joos & Richter. 1978, Höhere Mathematik für den Praktiker (Verlag Harri Deutsch−Thun−Frankfurt) [Google Scholar]
 Kainulainen, J., & Tan, J. C. 2013, A&A, 549, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kainulainen, J., Beuther, H., Henning, T., & Plume, R. 2009, A&A, 508, L35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kandori, R., Nakajima, Y., Tamura, M., et al. 2005, AJ, 130, 2166 [NASA ADS] [CrossRef] [Google Scholar]
 Keto, E., & Caselli, P. 2010, MNRAS, 402, 1625 [NASA ADS] [CrossRef] [Google Scholar]
 Keto, E., Rawlings, J., & Caselli, P. 2014, MNRAS, 440, 2616 [NASA ADS] [CrossRef] [Google Scholar]
 Kirk, H., Johnstone, D., & Di Francesco, J. 2006, ApJ, 646, 1009 [NASA ADS] [CrossRef] [Google Scholar]
 Konstandin, L., Girichidis, P., Federrath, C., & Klessen, R. S. 2012, ApJ, 761, 149 [NASA ADS] [CrossRef] [Google Scholar]
 Kritsuk, A. G., Norman, M. L., Padoan, P., & Wagner, R. 2007, ApJ, 665, 416 [NASA ADS] [CrossRef] [Google Scholar]
 Kritsuk, A. G., Norman, M. L., & Wagner, R. 2011, ApJ, 727, L20 [NASA ADS] [CrossRef] [Google Scholar]
 Lazarian, A., & Pogosyan, D. 2000, ApJ, 537, 720 [NASA ADS] [CrossRef] [Google Scholar]
 McCrea, W. H. 1957, MNRAS, 117, 562 [NASA ADS] [CrossRef] [Google Scholar]
 Nordlund, Å. P., & Padoan, P. 1999, in Insterstellar Turbulence, eds. J. Franco, & A. Carramiñana (Cambridge University Press), 218 [Google Scholar]
 Ostriker, E. C., Stone, J. M., & Gammie, C. F. 2001, ApJ, 546, 980 [NASA ADS] [CrossRef] [Google Scholar]
 Ostriker, J. 1964, ApJ, 140, 1056 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Padoan, P., Jones, B. J. T., & Nordlund, A. P. 1997a, ApJ, 474, 730 [NASA ADS] [CrossRef] [Google Scholar]
 Padoan, P., Nordlund, A., & Jones, B. J. T. 1997b, MNRAS, 288, 145 [NASA ADS] [CrossRef] [Google Scholar]
 Passot, T., & VázquezSemadeni, E. 1998, Phys. Rev. E, 58, 4501 [NASA ADS] [CrossRef] [Google Scholar]
 Price, D. J., Federrath, C., & Brunt, C. M. 2011, ApJ, 727, L21 [NASA ADS] [CrossRef] [Google Scholar]
 Schneider, N., Csengeri, T., Hennemann, M., et al. 2012, A&A, 540, L11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Schneider, N., Csengeri, T., Hennemann, M., et al. 2013, A&A, 551, C1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Stodólkiewicz, J. S. 1963, Acta Astron., 13, 30 [NASA ADS] [Google Scholar]
 VazquezSemadeni, E. 1994, ApJ, 423, 681 [NASA ADS] [CrossRef] [Google Scholar]
 VázquezSemadeni, E., & García, N. 2001, ApJ, 557, 727 [NASA ADS] [CrossRef] [Google Scholar]
Online material
Appendix A: Approximations of a BonnorEbert sphere
Here, the mass, radius, and mean mass surface density of spheres with an analytical profile given by Eq. (1) are compared with the correct values of BonnorEbert spheres. They are shown in Fig. A.1. Characteristic parameters are listed in Table A.1.
Fig. A.1 Mass, radius, and mean mass surface density of spherical clouds with fixed bounding pressure p(z_{sph}) and temperature as a function of the pressure ratio p_{c}/p(z_{sph}). The values of a BonnorEbert sphere are compared with the corresponding values of spherical clouds with smooth density profiles as given in Eq. (1) with n = 2 and n = 3 (dasheddotted and dashed line, respectively). The dotted vertical and horizontal lines mark the values at maximum radius (open circle) and at critical stability (filled circle) of BonnorEbert spheres. The mass, radius, and mean mass surface density of a BonnorEbert sphere in the limit of infinite overpressure are shown as gray horizontal lines. 

Open with DEXTER 
Appendix A.1: Mass of BonnorEbert spheres
The mass of a BonnorEbert sphere is given by (Fischera & Dopita 2008) (A.1)where the size θ_{cl} is determined by pressure equilibrium p_{c}e^{−ω(θcl)} = p(θ_{cl}) = p_{ext}.
In Fig. A.1 the mass of BonnorEbert spheres is compared with the masses of spheres using the approximation n = 3 and n = 2, which are given by (A.2)and (A.3)At given external pressure and temperature, the mass of a BonnorEbert sphere increases with increasing overpressure up to the critical value where the sphere has the highest mass. For higher overpressure the mass is lower than the critical value and fluctuates to an asymptotic value at infinite overpressure given by (A.4)
Appendix A.1.1: Critical stability
BonnorEbert spheres have a wellknown critically stable configuration related to the response of the gas pressure at cloud radius in relation to compression. As long as the compression leads to a pressure increase (dp(r_{cl})/dr_{cl}< 0), a cloud is considered stable, otherwise as unstable. Critical stability is given for (A.5)As the derivative of the radius with respect to q is nonzero at the pressure maximum, the critical condition is (A.6)For BonnorEbert spheres, a critically stable sphere of the first pressure maximum, which is also the global maximum, is characterized through an overpressure of 14.04. When we apply the same criterion to spheres with an analytical profile as given by Eq. (1) with n = 3, we find from Eq. (A.2) for fixed mass and fixed K a critical overpressure , close to the critical value of a BonnorEbert sphere. The physical parameters corresponding to the critical value of a sphere with n = 3 are listed in Table A.1.
As pointed out in the paper of Fischera & Martin (2012b), the critical mass of a BonnorEbert sphere is also the highest possible mass for fixed K and external pressure p_{ext} , but varying pressure ratio q. The same applies for a sphere with n = 3. A cylinder, by comparison, has no critical stability with regard to compression. The highest mass line density is related to infinite overpressure. The situation for spheres with a n = 2profile would be the same as for cylinders.
Appendix A.2: Radius of BonnorEbert spheres
The radius of the BonnorEbert sphere is given by (A.7)For a cloud with a truncated analytical density profile we find from Eq. (1) and the inner radius where A can be expressed through A^{2} = 4πGp_{ext}/K^{2}/q,(A.8)For a given q and p_{ext} the radius increases proportionally to K. Clouds with the same q and K will be smaller in higher pressure regions. Clouds with n> 2 have the greatest extension at an overpressure (A.9)Figure A.1 shows that the radius of the BonnorEbert sphere behaves similarly as a function of overpressure as does the mass. However, the cloud with the greatest extension is subcritical with an overpressure of 4.990. The greatest extension of a sphere with n = 3 corresponds to an overpressure (A.10)which is close to the corresponding value of BonnorEbert spheres. For a higher overpressure, the radius of a BonnorEbert sphere is generally smaller than the radius at critical stability, but at a given q it is larger than a sphere with an analytical density profile with n = 3. In the limit of high overpressure the radius of the BonnorEbert sphere fluctuates asymptotically toward (A.11)which is identical with the radius of a sphere with n = 2 in the limit of infinite overpressure.
Appendix A.3: Mean mass surface density
It has been shown (Fischera & Martin 2012b) that for overpressures where the gas pressure at the cloud boundary of a sphere with fixed mass has local maxima for varying overpressure the mean mass surface density of BonnorEbert spheres is simply given by (A.12)The mean value applies, for example, to a critically stable sphere and for a sphere with infinite overpressure (Fig. A.1). For a sphere with infinite overpressure this can be easily verified by replacing the density ratio e^{−ω(θ)} by the asymptotic density profile at large scales θ, which is shown in Sect. 2.1.1 to be 2 /θ^{2}.
Characteristic cloud parameters.
Appendix B: Probability distribution function of the local density
In this section we provide for the sake of completeness the PDF of the local density of spheres and cylinders with a truncated density profile as given by Eq. (1).
The logarithmic PDF is given by (B.1)where V_{cl} is the cloud volume and where V(r) is the cloud volume within radius r.
For a truncated density profile with a pressure ratio p_{ext}/p_{c} = q we obtain for spheres (B.2)In the limit of high overpressure this becomes for (ρ/ρ_{c})^{2 /n} ≪ 1 a power law (B.3)or ρP_{sph}(ρ) ∝ ρ^{−3 /n}, as expected for simple radial density profiles ρ ∝ r^{−n} (Kritsuk et al. 2011; Federrath et al. 2011).
For cylinders the PDF of the local density is given by (B.4)The PDF for a given q over the full range is a power law. For isothermal selfgravitating cylinders we obtain .
Appendix B.1: Local density PDF of BonnorEbert spheres
Fig. B.1 Local density PDF of critical and supercritical BonnorEbert spheres for three different pressure ratios, given in units of the overpressure of a critically stable BonnorEbert sphere. They are directly compared with the PDFs of spheres with density profiles as given in Eq. (1) with n = 3 and n = 2. Their corresponding powerlaw asymptotes in the limit of infinite overpressure are shown as gray lines. 

Open with DEXTER 
Replacing in Eq. (B.1) the density by the solution of the LaneEmden equation ρ = ρ_{c}e^{−ω(θ)} and the radius by the unitfree size θ = rA provides for the local density PDF of a BonnorEbert sphere (B.5)The PDFs for critically stable and supercritical BonnorEbert spheres (here defined as a selfgravitating sphere with an overpressure p_{c}/p_{ext}> 14.04) are shown in Fig. B.1. We showed in Sect. 2.1 that the inner part up to p_{c}/p(r) < ~ 100 of the density profile of selfgravitating spheres can be approximated by an analytical profile as given by Eq. (1) with n = 3. The shape of the PDF of the BonnorEbert spheres at high density values is consequently determined by this profile. For example, the PDF of BonnorEbert spheres at densities close to the central density value is approximately a power law with ρ^{1}. For supercritical BonnorEbert spheres, the PDF of density values ρ ≪ ρ_{c} asypmtotically approaches the PDF of a sphere with an analytical density profile with n = 2. The PDF becomes close to a power law ρ^{1.5}.
The values are close to the values obtained with collapse models. Kritsuk et al. (2011) studied the theoretical density distribution in starforming interestellar clouds. The collapse model provided a PDF that can be approximated over a wide range by a power law ρ^{1.695}. At highest densities the PDF showed a flatter shape approximated by a power law ρ^{1} and therefore the same as for BonnorEbert spheres. However, the physical explanation given are rotationally supported cores.
Appendix C: Estimating the maxima of P_{BE}(Σ) and ΣP_{BE}(Σ)
The maxima positions for linear and logarithmic PDFs of BonnorEbert spheres were derived by estimating the zero points of the first derivatives. The highest positions of the linear and logarithmic PDF fulfill the conditions
where the linear PDF of the BonnorEbert sphere is given by (C.3)
Appendix C.1: Linear PDF
Inserting the expression C.3 into the condition Eq. (C.1) for the maxima, we obtain (C.4)The condition becomes (C.5)where the expression in brackets is referred to as g(θ_{⊥}). From the derivative of the mass surface density Eq. (19), it follows that (C.6)where the integration along the radius θ has been changed to the integration along the depths at impact radius θ_{⊥}. The change of g(θ_{⊥}) by an infinitely small increase dθ_{⊥} is obtained by expanding the upper limit of the integral and the functions depending on θ_{⊥} into Taylor series. To first order, we find (C.7)where (C.8)For the derivative we obtain (C.9)
Appendix C.2: Logarithmic PDF
The condition for the highest positions of the logarithmic PDF is derived in a similar manner. Equation (C.2) provides (C.10)
Appendix D: Approximation of the P_{BE}(Σ)
In this appendix we show for the sake of completeness that for the approximations of BonnorEbert spheres expression 20 of the PDF is identical to Eqs. (21) and (22). Replacing in Eq. (20) the density profile e^{−ω(θ)} = ρ(θ) /ρ_{c} by the analytical density profile as given in Eq. (1) and the potential by ω = −ln [ ρ(θ) /ρ_{c} ] , we obtain (D.1)where we have changed the integration along the radius θ to the integration along the depth . By introducing the integration variable (D.2)and by replacing the projected radius and the cloud radius by θ_{⊥} = xθ_{cl} and the expression (D.1) can be transformed into (D.3)where y_{n} = (1 − q^{2 /n})(1 − x^{2}) and where the upper limit is given by . The exponent of the integrand can be reduced by using (simplified version of Eqs. (6)−(57) of Joos & Richter 1978) (D.4)which provides (D.5)After applying Eq. (8) for the unitfree mass surface density, we obtain (D.6)which is identical with Eqs. (21) and (22).
Appendix E: Asymptotes at high and low mass surface densities for qaveraged PDFs
In this appendix we provide the asymptotes at high and low mass surface densities of the mean PDF for a considered distribution of the pressure ratios q of an ensemble of isothermal, selfgravitating, and pressurized spheres or cylinders as shown in Figs. 9 and 10. For the sake of simplicity, for BonnorEbert spheres the analytical approximations for low and high overpressure are considered. The external pressure p_{ext} is considered the same for all clouds in the sample.
For convenience, we introduce the unitfree radius defined by (E.1)The mean PDF as defined by Eq. (52) at a given mass surface density Σ_{n} is then given by the integral (E.2)where P(q^{1}) is the distribution function of the overpressures q^{1} given by Eq. (53), P_{cl}(Σ_{n},q) the conditional PDF under the pressure ratio q, and where is a normalization constant given by (E.3)where κ = 2 for spheres and κ = 1 for cylinders. The lowest overpressure is related to the central mass surface density with Σ_{n}(0) = Σ_{n}.
Appendix E.1: Sphere
Appendix E.1.1: Asymptote at high mass surface densities
As shown in Paper I and in Sect. 2.3.3, the PDF of spheres with high overpressure asymptotically approaches within the limit of high mass surface densities the power laws given by Eq. (38). For the product we derive(E.4)Because the probabilities at high mass surface densities are related to spheres with high overpressure (Sect. 2.3), we can simplify the probability distribution Eq. (53) of the overpressure q^{1} through the corresponding powerlaw asymptote in the limit , so that . For the integral (E.2) we obtain (E.5)where we have replaced the lowest overpressure by the central mass surface density using Eq. (14).
The estimate of the normalization constant (Eq. (E.3)) is straightforward and is(E.6)where (E.7)is the normalized incomplete betafunction with and where (E.8)
Appendix E.1.2: Asymptote at low mass surface densities
In a similar manner as in the previous section, we can derive an asymptote of the mean PDF for low mass surface densities. As discussed in Sect. 2.3.2, the PDF of individual clouds approaches in the limit of low mass surface densities a simple power law given by Eq. (34). We find for the product of the square of the unitfree radius and the PDF the approximation (E.9)The value does not depend on the overpressure q^{1}.
In the limit of low mass surface densities we can approximate . The mean asymptotic value of the PDF is then given by (E.10)with and (E.11)The slope of the asymptotic power law for low mass surface densities is the same as for individual clouds.
Appendix E.2: Cylinder
Appendix E.2.1: Asymptote at high mass surface densities
To estimate the mean value of the PDF for high mass surface densities we can apply the approximation Eq. (49). For the product of radius and PDF we obtain (E.12)where we have replaced the central mass surface density by the asymptotic value given by Eq. (50) valid in the limit of high overpressure. The same replacement has been made for the fixed mass surface density Σ_{4} in the denominator, which is related to the lowest overpressure . To estimate the integral (E.2) we can make the same simplification as in Appendix E.1.1 by replacing the probability distribution Eq. (53) by the corresponding power law . For the mean PDF of selfgravitating cylinders we obtain as the asymptote at high mass surface densities (E.13)To provide a simple analytical estimate for the normalization constant we can consider the two extreme cases we discussed here. For a distribution in q^{1} with the radius can be simplified to . The estimate of the integral (E.3) is again straightforward and the normalization constant becomes (E.14)where (E.15)For , the distribution of overpressure is a simple power law and the integral can be solved without further simplification. The normalization constant becomes (E.16)
Appendix E.2.2: Asymptote at low mass surface densities
In the limit of low mass surface densities it follows from Eqs. (34) and (41) with (42) that the asymptote of the PDF in the limit of low mass surface densities is given by (E.17)For low mass surface densities the mean PDF is dominated by x ≈ 1. For the product of radius and PDF we obtain (E.18)For large the contribution to the mean PDF is related to high overpressures so that we can simplify the expression by removing the square root using . The asymptote of the mean PDF at low mass surface densities becomes (E.19)where (E.20)For the limit , where the probability distribution P(q^{1}) becomes a simple power law, we find (E.21)
All Tables
All Figures
Fig. 1 Normalized radial density profiles of some cases considered in this study given as a function of the unitfree radius θ = rA where A^{2} = (4πGρ_{c}) /K. The filled symbol marks the condition for a critically stable sphere where θ_{crit} ≈ 6.451 and ρ_{c}/ρ(θ_{crit}) ≈ 14.04. The curves of the analytical density profile given by Eq. (1) are labeled with their asymptotic behavior for . 

Open with DEXTER  
In the text 
Fig. 2 Profiles of the mass surface density through BonnorEbert spheres and spheres with an analytical density profile as given in Eq. (1). The gray curves are profiles for different cloud radii θ_{cl}. For BonnorEbert spheres the cloud radius corresponds to q = q_{crit} (a), q = 10^{2}q_{crit} (b), and q = 10^{4}q_{crit} (c), where q_{crit} = 1/14.04. For each assumed cloud radius the central pressure p_{c} in the considered spheres is the same. The profiles are normalized by the corresponding central mass surface density through a BonnorEbert sphere. 

Open with DEXTER  
In the text 
Fig. 3 Probabilities of a number of mass surface densities through BonnorEbert spheres shown as a function of the central mass surface density. The curves are labeled by the corresponding mass surface density given in units of . The probabilities are compared with those of analytical density profiles as given in Eq. (1) with n = 3 and n = 2. The probabilities for n = 3 are only shown for where the analytical density profile becomes a valid approximation of the density profile of a BonnorEbert sphere. The overpressures given in the upper axis are only valid for BonnorEbert spheres. The gray lines show the asymptotic probability for BonnorEbert spheres in the limit of infinite overpressure. The vertical dotted line marks the location of the critically stable BonnorEbert sphere. 

Open with DEXTER  
In the text 
Fig. 4 PDF of the mass surface density of BonnorEbert spheres embedded in a medium with p_{ext}/k = 2 × 10^{4} K cm^{3} for different gravitational states (solid black curves). The curves correspond to overpressures with i = −1,0,1 where is the overpressure of a critically stable BonnorEbert sphere. The PDFs are directly compared with those of spheres with a simple density profile (Eq. (1)) with n = 3 (dotted curves). The black dashed lines are asymptotes of the PDF Σ_{3}P_{sph}(Σ_{3}) in the limit of high mass surface densities (Eq. (39)). The PDF of a BonnorEbert sphere in the limit of infinite overpressure is shown as a gray curve and is derived assuming a sphere with an analytical density profile with n = 2. The gray dashed lines are the corresponding powerlaw asymptotes in the limit of low (Eq. (34)) and high (Eq. (40)) mass surface densities. 

Open with DEXTER  
In the text 
Fig. 5 Similar to Fig. 4, but for a wider range of pressure ratios. We show the transition from a PDF of a critically stable BonnorEbert sphere to the PDFs of highly supercritical spheres. The overpressure is varied in units of the critical value . Also shown are the corresponding PDFs for a sphere with the analytical density profile with n = 2 and n = 3. The powerlaw asymptotes close to the truncation point (black dashed lines) are scaled to the PDF of the Bonnor Ebert spheres at central mass surface density. 

Open with DEXTER  
In the text 
Fig. 6 Normalized impact parameter (lefthand panel) and mass surface density (righthand panel) at the maximum of the PDF of the mass surface densities of BonnorEbert spheres, given as a function of the central mass surface density. The gray and black curves correspond to the logarithmic and the linear pdf. The values are compared with the corresponding impact parameters and mass surface densities of pressurized spheres with a density profile as given in Eq. (1) with n = 2 and n = 3. At low central mass surface densities the maximum of the PDF is at impact parameter [ x ] _{max} = 0 and the curves in the righthand panel follow the long dashed line where [ Σ ] _{max} = Σ(0). The location of the critically stable BonnorEbert sphere is marked by a vertical dotted line and the corresponding normalized impact parameters and mass surface densities at the PDF maximum as filled circles. The horizontal dotted lines in the left panel are the asymptotic values of the impact parameters of the PDF maximum of BonnorEbert spheres in the limit of high overpressure. The overpressures shown in the upper axis are only valid for the BonnorEbert sphere. 

Open with DEXTER  
In the text 
Fig. 7 PDF of the mass surface density of isothermal selfgravitating pressurized cylinders for various overpressures. The vertical line marks the pole position of the PDFs and is labeled with the corresponding pressure ratio q^{1}. The asymptotes at high and low mass surface density are shown as black dashed lines. The gray curves are the corresponding asymptotic PDFs (Eq. (43)) where the pole is removed. The powerlaw asymptotes at high mass surface densities are shown as gray dashed lines. 

Open with DEXTER  
In the text 
Fig. 8 Same as Fig. 6, but for the asymptotic PDF of isothermal selfgravitating cylinders where the PDF has no pole at the highest mass surface density. 

Open with DEXTER  
In the text 
Fig. 9 Mean PDF of the mass surface density for an ensemble of cores assumed to be BonnorEbert spheres with a distribution of gravitational states or pressure ratios q^{1} as given in Eq. (53). For comparison the mean PDF is also derived for spheres with an analytical density profile as given in Eq. (1) with n = 3. The characteristic pressure ratio of the probability distribution of the gravitational states is either q_{0} = q_{crit} ~ 1/14.04 of a critically stable sphere (top panel) or q_{0} ≫ 1 (bottom panel). The various curves labeled a), b), and c) correspond to different power indices k_{2} assumed to be 2.0, 1.5, and 1.2. The power k_{1} = 2 and the smoothness parameter γ = 2 are the same in all calculations. The black dashed lines are powerlaw approximations of the PDFs in the limit of low and high mass surface densities for an analytical density profile with n = 3 (Appendix E.1). The external pressure is assumed to be p_{ext}/k = 2 × 10^{4} K cm^{3}. Also shown are the PDFs of a critically stable sphere and a sphere with infinite overpressure (gray lines and gray annotation). The corresponding powerlaw asymptotes of the PDFs of single cores in the limit of low (Eq. (34)) and high (Eq. (39) and Eq. (40)) mass surface densities are shown as gray dashed lines. 

Open with DEXTER  
In the text 
Fig. 10 Mean PDF of the column density for an ensemble of isothermal selfgravitating cylinders with a variety of gravitational states or overpressures q^{1}, but fixed inclination angle i. The probability distribution of the overpressures is given by Eq. (53) with fixed q_{0}, k_{1} = 2, and γ = 2 where q_{0} is either 1/14.04 (lefthand figure) or q_{0} ≫ 1 (righthand figure). The power index k_{2} is assumed to be 1.5, 2, and 3. The cylinders are pressurized by a medium with p_{ext}/k = 2 × 10^{4} K cm^{3}. The black dashed lines are asymptotes at low and high mass surface densities (Appendix E.2). For comparison, the PDF of a single cylinder with overpressure q^{1} = 14.04 is shown (gray line). The gray dotted line is the asymptotic PDF as defined in Eq. (43), and the gray dashed lines are the corresponding asymptotes at low (Eq. (47)) and high (Eq. (51)) mass surface densities. 

Open with DEXTER  
In the text 
Fig. 11 Angleaveraged PDF of the column density for cylinders with a mean inclination angle of 0° (lefthand figure) and 60° (righthand figure). The cosine of the inclination angle is assumed to be Gaussian distributed. The curves are labeled with the assumed standard deviation σ_{μ}. The cylinders are pressurized by a medium with p_{ext}/k = 2 × 10^{4} K cm^{3} and have an overpressure p_{c}/p_{ext} = 14.04. The intrinsic PDF for the limit σ_{μ} → 0 is shown as dotted curves. The gray vertical lines mark the highest extinction values of a cylinder seen under the assumed inclination angle and the position of the peak of the averaged PDF of randomly oriented cylinders. The thick curves are the PDFs for randomly distributed cylinders. 

Open with DEXTER  
In the text 
Fig. 12 Combined PDF of a molecular cloud consisting of a turbulent medium (dotted curves) and uniformly distributed condensed cores (dasheddotted curves). The cores are either critically stable (left) or in a state of collapse (right) where the density profiles are described using the corresponding approximations for BonnorEbert spheres. The turbulent medium is modeled through a simple lognormal function and the pressure in the turbulent gas is assumed to be p_{ext}/k = 2 × 10^{4} K cm^{3}. The spheres do not overlap and cover f_{A} = 40% of the area of the molecular cloud. The distributions are labeled with the assumed ratio ⟨Σ_{b}⟩/Σ(0) of the mean background ⟨Σ_{b}⟩ and the central mass surface density Σ(0) through a critically stable BonnorEbert sphere. The scaled intrinsic nonbackgroundaveraged PDFs of BonnorEbert spheres (see Fig. 4) are shown as solid gray lines. Also given are the corresponding powerlaw asymptotes at low (Eq. (34)) and at high (Eq. (39) and Eq. (40)) mass surface densities. 

Open with DEXTER  
In the text 
Fig. A.1 Mass, radius, and mean mass surface density of spherical clouds with fixed bounding pressure p(z_{sph}) and temperature as a function of the pressure ratio p_{c}/p(z_{sph}). The values of a BonnorEbert sphere are compared with the corresponding values of spherical clouds with smooth density profiles as given in Eq. (1) with n = 2 and n = 3 (dasheddotted and dashed line, respectively). The dotted vertical and horizontal lines mark the values at maximum radius (open circle) and at critical stability (filled circle) of BonnorEbert spheres. The mass, radius, and mean mass surface density of a BonnorEbert sphere in the limit of infinite overpressure are shown as gray horizontal lines. 

Open with DEXTER  
In the text 
Fig. B.1 Local density PDF of critical and supercritical BonnorEbert spheres for three different pressure ratios, given in units of the overpressure of a critically stable BonnorEbert sphere. They are directly compared with the PDFs of spheres with density profiles as given in Eq. (1) with n = 3 and n = 2. Their corresponding powerlaw asymptotes in the limit of infinite overpressure are shown as gray lines. 

Open with DEXTER  
In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.