Open Access
Issue
A&A
Volume 695, March 2025
Article Number L25
Number of page(s) 7
Section Letters to the Editor
DOI https://doi.org/10.1051/0004-6361/202452427
Published online 25 March 2025

© The Authors 2025

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

Molecular clouds fragment into dense, compact cores, with gravitationally bound cores being widely recognized as the primary sites of protostellar formation. The process of gravitational fragmentation is considered key to core creation in molecular clouds (e.g. Larson 1985; André et al. 2010). Recent observations from the Atacama Large Millimeter/submillimeter Array (ALMA) appear to have provided empirical support for this concept. These observations show that the core separations within high-mass, cluster-forming molecular clumps are in reasonable agreement with the predictions of thermal Jeans fragmentation (Klaassen et al. 2018; Beuther et al. 2022; Lu et al. 2020; Zhang et al. 2022; Ishihara et al. 2024; Sanhueza et al. 2019; Morii et al. 2024). This suggests that thermal gravitational, or thermal Jeans, fragmentation governs core formation in high-mass regions.

In particular, Ishihara et al. (2024) used 1.3mm dust continuum images of 30 high-mass cluster-forming clumps observed as part of the DIHCA (Digging into the Interior of Hot Cores with ALMA) survey, identifying around 103 dense cores. Their results show that the typical core separation is about 7800 au, which is consistent with the thermal Jeans length of the parent clumps: L J = c s π / G ρ $ L_J= c_{\mathrm{s}} \sqrt{\pi/G\rho} $, where cs is the sound speed, G is the gravitational constant, and ρ is the density. The corresponding core mass, given by the Jeans mass, is MJ = 4πρ(LJ/2)3/3.

These results suggest that thermal Jeans fragmentation plays a significant role in core formation within dense, cluster-forming regions. However, it remains unclear whether the same dynamics apply in other environments. In several nearby molecular clouds, such as Orion A and ρ-Ophiuchus, most dense cores have been found to be gravitationally unbound and pressure-confined (Maruta et al. 2010; Kirk et al. 2017; Takemura et al. 2023). The mechanisms responsible for the formation of these unbound cores remain uncertain.

This paper focuses on nearby low-mass clouds, specifically the Polaris Flare (hereafter Polaris) and Lupus I, where the role of cloud self-gravity is relatively diminished compared to the influence of supersonic turbulence, which may even dominate over self-gravity. To achieve a deeper understanding of core formation in various cloud environments, we apply methods similar to those used by Ishihara et al. (2024) to investigate the physical processes governing core formation in these less dense, low-mass molecular clouds.

We consider three key fragmentation mechanisms: thermal Jeans fragmentation, turbulent Jeans fragmentation, and turbulent fragmentation. In thermal Jeans fragmentation, the balance between thermal pressure and self-gravity determines how gravitationally bound structures form. The characteristic mass and length scales in this scenario are the Jeans mass (MJ) and Jeans length (LJ), which depend on the cloud’s temperature and density (Jeans 1928; Larson 1985).

Turbulent Jeans fragmentation is a variation of thermal Jeans fragmentation, whereby turbulent pressure supports the cloud, leading to larger characteristic mass and length scales than in the thermal case. However, this model does not account for the full properties of cloud turbulence and is only applicable to larger structures for which turbulent motions can be approximated as isotropic (e.g. Mac Low & Klessen 2004; Palau et al. 2015).

In turbulent fragmentation, dense cores form through the dynamical compression of local turbulence (Padoan & Nordlund 2002). The density probability distribution function (PDF) in turbulent clouds tends to follow a log-normal distribution due to repeated compression (Padoan et al. 2020). In a turbulent medium with a power-law energy spectrum, there are no inherent characteristic mass and length scales. However, in the interstellar medium, the power law breaks at the sonic scale, lS, where turbulent velocity dispersions become comparable to the sound speed (Federrath & Klessen 2012). For nearby molecular clouds, this scale is typically measured to be lS ∼ 0.1 pc (≈20 000 au). Therefore, the sonic scale, lS, may represent the characteristic length for turbulent fragmentation, with the corresponding mass (MS) estimated as MS ≈ 4πρ⟩(lS/2)3/3, respectively, where ⟨ρ⟩ is the average cloud density. Since these three fragmentation mechanisms predict different characteristic mass and length scales, we aim to constrain the dominant fragmentation process in these nearby clouds by deriving core masses and separations from observational data.

The paper is structured as follows: Sect. 2 details the data used in this study. In Sect. 3, we present the results of core identification and their physical properties. We then compare these core properties with predictions from gravitational fragmentation and turbulent fragmentation models in Sects. 4 and 5, respectively, and discuss the physical processes involved in core formation in these molecular clouds in Sect. 6.

2. Targets, data, methods, and core identification

From the Herschel Gould Belt Survey (André et al. 2010), we have chosen two clouds detailed in Table 1: Polaris and Lupus I. High-resolution (18″) column density maps of these clouds were obtained from the Herschel Gould Belt Survey’s website1. At the distances to the clouds (355 pc for Polaris, and 182 pc for Lupus I), the 18″ resolution corresponds to 6390 au and 2700 au, respectively (see Fig. A.1). The average density was computed by dividing the total cloud mass by the volume of a sphere with a radius of R cl = Area / π $ R_{\mathrm{cl}} = \sqrt{\mathrm{Area}/\pi} $, where the cloud area was determined by that of the trunk identified by astrodendro. As is shown below, these clouds exhibit a sufficient number (several hundred or more) of identified dense cores, using astrodendro.

Table 1.

Physical properties of target clouds.

Polaris and Lupus I are categorized as a quiescent region and a low-mass star-forming region, respectively, according to Schneider et al. (2022). In Lupus I, nine protostars are identified (Benedettini et al. 2018), whereas Polaris has not formed protostars yet. Both clouds share similar average column densities and temperatures, with 〈Ncl〉 ≈ 8 − 12 × 1020 cm−2 and 〈Tcl〉 ≈ 15 K. As a result, the Jeans mass of the clouds is estimated to be around 70–80 M, with a corresponding Jeans length of about 3–4 pc. However, the degree of density concentration is different (see Fig. A.1). The molecular gas in Polaris appears to be more uniformly and spatially extended. The gravitational boundness of the cloud is likely weaker, given that the virial parameters for the Polaris and Lupus I clouds are estimated to be about 3.1 and 2.1, respectively, based on an observed velocity dispersion of 1.51 km s−1 and 1.46 km s−1 (Spilker et al. 2022).

We applied astrodendro (Rosolowsky et al. 2008) to the column density maps with the threshold value of min_value (Smin = 5σ), the minimum step of min_delta (δmin = 3σ), and the minimum pixel number of min_npix (θmin was set to the number of pixels equivalent to the beam area). Here, we adopted 1-σ noise levels of 7.67 × 1019 cm−2 for Polaris and 8.02 × 1019 cm−2 for Lupus I, respectively (for this choice of parameters, see Appendix E). The astrodendro identifies three structures: a leaf, branch, and trunk. We defined the minimum structure of the ‘leaf’ as a dense core and derived the core mass by subtracting the background gas components, which are defined as the column density of the corresponding trunk component. Following Takemura et al. (2021), we used additional condition to select cores: the peak column density of the core should be more than 2× min_value (see also Shimajiri et al. 2015). The positions of the identified cores are indicated in Fig. A.1. Our result show that in Lupus I, nine protostars are identified (compared with the catalogue by Benedettini et al. 2018), whereas Polaris has not formed protostars yet. This is consistent with previous studies (Benedettini et al. 2018; Ward-Thompson et al. 2010). A more detailed comparison is given in Appendix C. The core identification also depends on the parameters chosen. In Appendix E, we examine the dependence of astrodendro parameters. In Appendix C, we also compare the core properties with the ones of Benedettini et al. (2018), who used getsources.

Then, we applied the nearest neighbour search (NNS) method to derive the core separations on the 2D H2 column density image (Knuth 1973). The core separations (S3D) were obtained by multiplying a correction factor of 4/π (see Appendix F).

3. Physical properties of dense cores

We identified the dense cores; their physical properties are summarized in Table 2. Here, the core radius, Rcore, was calculated as half of the geometric mean of the FWHMmajor and FWHMminor provided by astrodendro; specifically, R core = 0.5 FWHM major × FWHM minor $ R_{\mathrm{core}} = 0.5 \sqrt{\mathrm{FWHM}_{\mathrm{major}} \times \mathrm{FWHM}_{\mathrm{minor}}} $. The core radius was then beam-deconvolved. The core density, ncore, was estimated under the assumption of spherical symmetry using ncore = Mcore/(4πRcore3/3). The identified cores were classified into two groups: gravitationally bound cores and unbound cores, based on the Bonnor-Ebert ratio ( α BE = M BE crit / M core $ \alpha_{\mathrm{BE}} = M_{\mathrm{BE}}^{\mathrm{crit}}/M_{\mathrm{core}} $), which compares the core mass (Mcore) to the critical Bonnor-Ebert mass (see Sect. 4.1 of Könyves et al. (2010) for the definition and discussion of M BE crit $ M_{\mathrm{BE}}^{\mathrm{crit}} $.). The number of identified cores, 𝒩core, in each cloud significantly exceeds the cloud Jeans number (NJ, cl = Mcl/MJ, cl), suggesting that cloud structures are hierarchical and that cores primarily form in denser regions. This implies that core formation may not be predominantly controlled by gravitational fragmentation, or may be a combination of both hierarchical structures and other processes.

Table 2.

Physical properties of identified cores.

In Polaris, the cores are somewhat larger than the ones in Lupus I (see also Fig. 1). The average core radii are estimated to be 0.077 pc for Polaris and 0.037 pc for Lupus I, while the average core masses are 0.28 M and 0.24 M, respectively. The core density tends to be higher in Lupus I (see also the left panel of Fig. 1).

thumbnail Fig. 1.

Core mass versus radius diagram (left) and histograms of core mass (middle) and core radius (right) for Polaris and Lupus I by astrodendro. The red circle and blue cross symbols represent the cores in Polaris and Lupus I, respectively. The error bar was calculated by Poisson statistics.

To judge the gravitational boundness of cores, we derived αBE (André et al. 2014). In Polaris, only seven gravitationally bound cores were found, whereas about 18% of the cores (≈80) in Lupus I are gravitationally bound. This is consistent with the fact that Lupus I exhibits star formation activity in its denser regions.

4. Comparison with gravitational fragmentation

The upper panel of Fig. 2 shows the histogram of core separations for the two clouds. Interestingly, the distribution of core separation is quite similar to each other, peaking at ∼1 − 2 × 104 au (≈0.1 − 0.2 pc). Given that the Jeans length for these clouds is about 4 pc (=8 × 105 au), the peak separation is about 20–40 times smaller than the cloud Jeans length. The density at which the Jeans length coincides with the peak separation is estimated to be of the order of 105 cm−3. This density is significantly higher than what is observed in these clouds.

thumbnail Fig. 2.

Histogram of core separation (left) and core mass versus separation diagrams in physical scale (middle) and normalized scale (right). The red circle and blue cross symbols represent the cores in Polaris and Lupus I, respectively. For the normalized scale, the core mass was divided by the local thermal Jeans mass (mJ) and the core separation was divided by the local thermal Jeans length (lJ), where mJ and lJ are estimated in the trunks that include the corresponding cores (leaves). The error bar was calculated by Poisson statistics.

To further compare with the thermal Jeans fragmentation model, we show in Fig. 2 the relationship between the core mass and the separation. For comparison, we plot the core mass and separation normalized to the local Jeans mass (mJ) and Jeans length (lJ), where mJ and lJ are estimated in the trunks that include the corresponding cores (leaves). If the cores form by thermal Jeans fragmentation, the cores should be preferentially distributed in the upper right part of Mcore/mJ ≥ 1 and S3D/lJ ≥ 1. There are no such cores found in that part. Therefore, thermal Jeans fragmentation seems to play a minor role in the core formation in these low-mass clouds.

For turbulent Jeans fragmentation, both the characteristic length scale and the mass scale become larger than those of thermal Jeans fragmentation, since the internal cloud turbulence contribute to additional dynamical support. The turbulent Mach numbers for both clouds are measured to be ∼10 (Spilker et al. 2022; Dame et al. 2001), resulting in the turbulent Jeans length being about three times greater than the thermal Jeans length. Therefore, the turbulent Jeans fragmentation is also in disagreement with the observed core properties.

Gravitational fragmentation is strongly influenced by cloud geometry. In the case of filamentary clouds, dense cores typically form with separations around four times the filament width, which is generally several times longer than the thermal Jeans length (see Ishihara et al. 2024, and references therein). Given the widely accepted filament width of 0.1 pc, this would predict core separations of about 0.4 pc–larger than the observed values in our analysis. Therefore, our findings contradict gravitational fragmentation scenarios and suggest that this mechanism is unlikely to explain the core formation in these clouds.

5. Comparison with turbulent fragmentation

5.1. Core separation and mass

In a turbulent medium, the energy power spectrum can be characterized by a power-law distribution. In such scale-free distribution, there are no characteristic length scales. However, in actual cloud turbulence, the thermal width becomes dominant over the nonthermal width at smaller scales. The scale below which the nonthermal width coincides with the thermal width is called a sonic scale, lS, (Goodman et al. 1998; Federrath & Klessen 2012; Federrath et al. 2021), and below lS, the thermal pressure becomes larger than the turbulent pressure.

Thus, the characteristic length scale of turbulent fragmentation corresponds to the cloud’s sonic scale, which has been estimated to be about 0.24 pc (∼47 000 au) based on CO observations for the Polaris (Ossenkopf & Mac Low 2002). This agrees well with the observed peak core separation of the Polaris. For Lupus I, there are no precise measurements from current observations, but the sonic scales in nearby clouds with similar environments (e.g. Taurus, Ophiuchus, and Orion) have been measured to be about 0.1 pc (Brunt 2010; Yun et al. 2021a,b). Therefore, the sonic scale of Lupus I is likely close to 0.1 pc (see also Eq. (13) of Federrath & Klessen 2012).

If the gas is accumulated from the surrounding area within lS/2, the characteristic mass scale is given as

M S = 4 π ρ ( l S / 2 ) 3 / 3 , $$ \begin{aligned} M_S = 4\pi \left\langle \rho \right\rangle (l_S/2) ^3/3 \ ,\end{aligned} $$(1)

where ⟨ρ⟩ (≈3.2 × 102 cm−3 for Polaris and ≈3.0 × 103 cm−3 for Lupus I) is estimated in the main branches identified by astrodendro. This mass scale is estimated to be 0.16 M for Polaris and 0.1 M for Lupus I, which agrees with those of the observations by a factor of a few. Therefore, we conclude that the core separation and mass are in good agreement with the prediction of turbulent fragmentation.

5.2. Core density PDFs

Another important prediction of turbulent fragmentation is its density PDF, which follows a log-normal form. The dispersion of the PDF was determined by the Mach number of cloud turbulence (ℳ) as:

σ ρ 2 = ln ( 1 + M 2 b 2 ) , $$ \begin{aligned} \sigma _{\rho }^2 = \ln (1 + \mathcal{M}^2 b^2), \end{aligned} $$(2)

where b is a constant representing the ratio between compressible and solenoidal components, typically ranging from b ≈ 1/3 − 1 (Federrath & Klessen 2012).

To investigate the effect of cloud turbulence, we calculated the density PDFs for both clouds. The density PDFs are sometimes modelled as a combination of a log-normal distribution and a power-law function at higher densities (Spilker et al. 2022; Schneider et al. 2022), with the power-law tail presumably influenced by cloud self-gravity. Here, we applied the fitting method performed by Schneider et al. (2015). In Schneider et al. (2015), the deviation point where the log-normal distribution transitions into a power-law distribution is defined as the residual between the fitted log-normal distribution and the measured PDF being greater than three times the statistical noise in each bin of the PDF. However, the sample size in this study was small and the statistical noise was large, so the conditions were relaxed from three times to one.

Core densities were calculated under the assumption of spherical symmetry. The resulting density PDFs are plotted in Fig. 3. The peak densities were estimated to be 1.7 × 103 cm−3 and 1.9 × 104 cm−3 for Polaris and Lupus I, respectively. The density PDFs of both clouds seem to be well fitted by single log-normal functions, having high-density tails. The dispersions of these density PDFs are measured as σρ = 1.09 for Polaris, and 1.64 for Lupus I, respectively. Using Eq. (2), we can deduce the expected turbulent Mach number for these clouds. The expected Mach numbers calculated are ℳ = 3.8 − 5.1 and 9.3 − 12.4 with b = 1/3 − 0.4 (solenoidal dominant forcing – mixing of solenoidal and compressive forcing; Federrath & Klessen 2012) for Polaris and Lupus I, respectively. For Polaris, Federrath et al. (2010) analyzed the density PDF, suggesting that the solenoidal component is dominant for the Polaris cloud, and this implies that b ≈ 1/3 may be reasonable. The turbulent Mach numbers derived from molecular line observations are in reasonable agreement with the deduced Mach numbers. Spilker et al. (2022) derived the Mach number using the CO(J = 1 − 0) data by Dame et al. (2001): ℳ = 12 for Polaris and ℳ = 11 for Lupus I.

thumbnail Fig. 3.

Histograms of mean volume density of identified core for Polaris (top) and Lupus I (bottom). The white and diagonal bins represent the sample of all cores and only bound cores, respectively. The solid red curves and dotted red line indicate the results of the log-normal fit and power-law fit. The error bar was calculated by Poisson statistics.

6. Conclusion

To understand the physics behind this process, we identified dense cores and derived their masses and separations in the two nearby low-mass clouds, Polaris and Lupus I. These clouds are less dense, and thus the effect of self-gravity is less important than in high-mass star-forming clumps.

The most prevalent cloud fragmentation mechanism is gravitational fragmentation, which can be divided into two types: thermal Jeans fragmentation and turbulent Jeans fragmentation. In thermal Jeans fragmentation, the characteristic length and mass scales are the Jeans length (lJ) and Jeans mass (mJ). In turbulent Jeans fragmentation, internal nonthermal motion’s support is included in the effective sound speed as c s eff = c s 1 + M 2 $ c_{\mathrm{s}} ^{\mathrm{eff}} = c_{\mathrm{s}} \sqrt{ 1+ {\cal M}^2} $, leading to larger characteristic length and mass scales than those of thermal Jeans fragmentation.

Recent studies have emphasized gravitational fragmentation of filaments as a primary mode of core formation (André et al. 2010). The fragmentation scale or separation for filament fragmentation is about four times the filament width, generally a few times longer than the thermal Jeans length (Ishihara et al. 2024). Therefore, the gravitational fragmentation scenarios – turbulent Jeans fragmentation and filament fragmentation – predict core masses and separations larger than those of thermal Jeans fragmentation.

Contrary to these predictions, we found that the core masses and separations in both Polaris and Lupus I are smaller than those expected from thermal Jeans fragmentation. The peak of core separation appears to be comparable to the sonic scale of ∼0.1 pc, below which thermal support becomes dominant, in both clouds. In addition, the density PDFs constructed from the core’s volume densities are fitted well, with log-normal distributions having dispersions estimated from the turbulent Mach numbers. These observational facts strongly indicate that the cores preferentially form by turbulent fragmentation.

In high-mass star-forming clumps, Ishihara et al. (2024) found that core separation is comparable to the thermal Jeans length, concluding that core formation is likely controlled by thermal Jeans fragmentation. This apparent difference between high-mass clumps and nearby low-mass clouds may suggest that different physical processes are responsible for core formation in different environments. Alternatively, the higher densities of high-mass star-forming clumps could mean that the characteristic mass scale at the sonic scale, MS, may already exceed or be at least comparable to the thermal Jeans mass. If this were the case, turbulent fragmentation might still be the primary driver of core formation, even in high-mass clumps. The sonic scales in some high-mass star-forming clumps are measured to be ∼0.01 pc (≈2 × 103 au), about ten times smaller than the ones of nearby molecular clouds (Ohashi et al. 2016; Li et al. 2023). These smaller sonic scales appear to be consistent with the core separation measured in high-mass star-forming clumps within a factor of a few (Ishihara et al. 2024). More precise measurements of the sonic scale in these regions could provide further insights into understanding the cloud fragmentation process in various environments.


Acknowledgments

This work was supported in part by The Graduate University for Advanced Studies, SOKENDAI. This work was financially supported by JSPS KAKENHI Grant Numbers JP22H01271 (PS), JP23H01218 (FN), and JP23H01221 (PS).

References

  1. André, P., Men’shchikov, A., Bontemps, S., et al. 2010, A&A, 518, L102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. André, P., Di Francesco, J., Ward-Thompson, D., et al. 2014, in Protostars and Planets VI, eds. H. Beuther, R. S. Klessen, C. P. Dullemond, & T. Henning, 27 [Google Scholar]
  3. Benedettini, M., Pezzuto, S., Schisano, E., et al. 2018, A&A, 619, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Beuther, H., Schneider, N., Simon, R., et al. 2022, A&A, 659, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Brunt, C. M. 2010, A&A, 513, A67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Dame, T. M., Hartmann, D., & Thaddeus, P. 2001, ApJ, 547, 792 [Google Scholar]
  7. Federrath, C., & Klessen, R. S. 2012, ApJ, 761, 156 [Google Scholar]
  8. Federrath, C., Roman-Duval, J., Klessen, R. S., Schmidt, W., & Mac Low, M. M. 2010, A&A, 512, A81 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Federrath, C., Klessen, R. S., Iapichino, L., & Beattie, J. R. 2021, Nat. Astron., 5, 365 [NASA ADS] [CrossRef] [Google Scholar]
  10. Goodman, A. A., Barranco, J. A., Wilner, D. J., & Heyer, M. H. 1998, ApJ, 504, 223 [Google Scholar]
  11. Ishihara, K., Sanhueza, P., Nakamura, F., et al. 2024, ApJ, 974, 95 [NASA ADS] [CrossRef] [Google Scholar]
  12. Jeans, J. H. 1928, Astronomy and Cosmogony (Cambridge University Press) [Google Scholar]
  13. Kirk, H., Friesen, R. K., Pineda, J. E., et al. 2017, ApJ, 846, 144 [Google Scholar]
  14. Klaassen, P. D., Johnston, K. G., Urquhart, J. S., et al. 2018, A&A, 611, A99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Knuth, D. E. 1973, The Art of Computer Programming. Vol. 3: Sorting and Searching, (Addison-Wesley Professional) [Google Scholar]
  16. Könyves, V., André, P., Men’shchikov, A., et al. 2010, A&A, 518, L106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Larson, R. B. 1985, MNRAS, 214, 379 [NASA ADS] [Google Scholar]
  18. Li, S., Sanhueza, P., Zhang, Q., et al. 2023, ApJ, 949, 109 [NASA ADS] [CrossRef] [Google Scholar]
  19. Lu, X., Cheng, Y., Ginsburg, A., et al. 2020, ApJ, 894, L14 [NASA ADS] [CrossRef] [Google Scholar]
  20. Mac Low, M.-M., & Klessen, R. S. 2004, Revi. Mod. Phys., 76, 125 [Google Scholar]
  21. Maruta, H., Nakamura, F., Nishi, R., Ikeda, N., & Kitamura, Y. 2010, ApJ, 714, 680 [NASA ADS] [CrossRef] [Google Scholar]
  22. Men’shchikov, A., André, P., Didelon, P., et al. 2010, A&A, 518, L103 [Google Scholar]
  23. Miville-Deschênes, M. A., Martin, P. G., Abergel, A., et al. 2010, A&A, 518, L104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Morii, K., Sanhueza, P., Zhang, Q., et al. 2024, ApJ, 966, 171 [Google Scholar]
  25. Ohashi, S., Sanhueza, P., Chen, H.-R. V., et al. 2016, ApJ, 833, 209 [NASA ADS] [CrossRef] [Google Scholar]
  26. Ossenkopf, V., & Mac Low, M. M. 2002, A&A, 390, 307 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Padoan, P., & Nordlund, Å. 2002, ApJ, 576, 870 [NASA ADS] [CrossRef] [Google Scholar]
  28. Padoan, P., Pan, L., Juvela, M., Haugbølle, T., & Nordlund, Å. 2020, ApJ, 900, 82 [NASA ADS] [CrossRef] [Google Scholar]
  29. Palau, A., Ballesteros-Paredes, J., Vázquez-Semadeni, E., et al. 2015, MNRAS, 453, 3785 [Google Scholar]
  30. Panopoulou, G. V., Clark, S. E., Hacar, A., et al. 2022, A&A, 657, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Robitaille, J. F., Motte, F., Schneider, N., Elia, D., & Bontemps, S. 2019, A&A, 628, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Rosolowsky, E. W., Pineda, J. E., Kauffmann, J., & Goodman, A. A. 2008, ApJ, 679, 1338 [Google Scholar]
  33. Sanhueza, P., Contreras, Y., Wu, B., et al. 2019, ApJ, 886, 102 [Google Scholar]
  34. Schneider, N., Ossenkopf, V., Csengeri, T., et al. 2015, A&A, 575, A79 [CrossRef] [EDP Sciences] [Google Scholar]
  35. Schneider, N., Ossenkopf-Okada, V., Clarke, S., et al. 2022, A&A, 666, A165 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Shimajiri, Y., Kitamura, Y., Nakamura, F., et al. 2015, ApJS, 217, 7 [NASA ADS] [CrossRef] [Google Scholar]
  37. Spilker, A., Kainulainen, J., & Orkisz, J. 2022, A&A, 667, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Takemura, H., Nakamura, F., Ishii, S., et al. 2021, PASJ, 73, 487 (Paper I) [NASA ADS] [CrossRef] [Google Scholar]
  39. Takemura, H., Nakamura, F., Arce, H. G., et al. 2023, ApJS, 264, 35 [Google Scholar]
  40. Ward-Thompson, D., Kirk, J. M., André, P., et al. 2010, A&A, 518, L92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Yun, H.-S., Lee, J.-E., Choi, Y., et al. 2021a, ApJS, 256, 16 [NASA ADS] [CrossRef] [Google Scholar]
  42. Yun, H.-S., Lee, J.-E., Evans, N. J., et al. 2021b, ApJ, 921, 31 [NASA ADS] [CrossRef] [Google Scholar]
  43. Zhang, C., Evans, N. J., Liu, T., et al. 2022, MNRAS, 510, 4998 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Core distribution in the clouds

Figs. A.1a and A.1b show the column density maps of Polaris and Lupus I, respectively. For the Polaris, we trimmed the original map to account for variations in background levels and noise across different regions. In the trimmed map, we used astrodendro to identify the cores. The identified cores are marked with crosses in each panel. For Lupus I, we excluded ∼ 24′ wide parts from the edges of the map due to the prominence of (artificial) scan patterns.

thumbnail Fig. A.1.

Column density distributions of target clouds. The H2 column density maps of Polaris (top) and Lupus I (bottom). Red star and black cross symbols are the positions of the protostellar cores identified by (Benedettini et al. 2018) and those of the dense cores identified by astrodendro (this work).

Appendix B: Distances to the targets

The target distances are recently updated by the Gaia data. In this paper, we adopted the distances, 355 pc and 182 pc, calculated by Panopoulou et al. (2022) for Polaris and Benedettini et al. (2018) for Lupus I, respectively. These distances are based on the Gaia data. For the Polaris distance, one of the positions in the Appendix of Panopoulou et al. (2022) (l, b) = (123.7, 24.8), overlaps with the Herschel map at (RA, Dec) = (01h58m11s,87d34m17s). Therefore, we adopted 355 pc.

Appendix C: Calculation of map noise levels

In our analysis, we calculated the background emission to be 3.80 × 1020 cm−2 ± 7.67 ×1019 cm−2 for Polaris and 2.16 × 1020 cm−2 ± 8.02 ×1019 cm−2 for Lupus I, respectively, based on the areas shown in Fig. A.1. For core identification, we subtracted these background components from the column density maps. The astrodendro algorithm was then applied to the processed maps, using the calculated standard deviations as the 1-σ noise levels. When estimating core masses, we added the subtracted background component back to the core mass. Notably, this background subtraction does not significantly influence the analysis presented in this paper.

The Herschel maps are also subject to cosmic infrared background (CIB) emission from unresolved high-redshift infrared galaxies (Miville-Deschênes et al. 2010; Robitaille et al. 2019). As a result, some of the identified cores could correspond to such sources. However, we consider the fraction of these "fake" cores to be negligible in our analysis.

It is worth mentioning that Ward-Thompson et al. (2010) identified only five cores in the Polaris. Their results agree with ours if we adopt their 1-σ threshold of ∼1.2 × 1021 cm−2. However, André et al. (2010) estimated the 1σ noise level of the Polaris image to be much lower, around 6 × 1019 cm−2, due to the significantly low background contamination in the Polaris region. This discrepancy in the number of identified cores is likely due to the different 1σ noise levels used in these studies, with many more cores being detected at the lower threshold. In fact, Men’shchikov et al. (2010) reported ∼300 starless cores in Polaris. If we adopt 1σ ∼ 1 × 1021 cm−2, our identified core number becomes ∼ 5, the same as that of Ward-Thompson et al. (2010). Many our cores in Polaris are distributed along the elongated structures seen in yellow–red in the right panel of Fig. A.1.

Appendix D: Comparison with the getsources core catalogue

In Figs. D.1 and D.2, we present the core masses, radii, and separations for the Lupus I cloud, comparing results derived using astrodendro (this paper) with those from the core catalogue of Benedettini et al. (2018), obtained via getsources. The mean values and overall distributions show good agreement, indicating that our conclusion regarding the importance of turbulent fragmentation is robust and independent of the core identification method. Furthermore, the typical cores identified by getsources exhibit smaller values of αBE (⟨αBE⟩≈ 0.25), which is consistent with the dominance of unbound cores (about 80%), similar to findings in Ophiuchus and Orion A.

thumbnail Fig. D.1.

Comparison of between astrodendro cores (this paper) and getsources cores (Benedettini et al. 2018). Core mass versus radius diagram (left) and histograms of core mass (middle) and core radius (right) identified by Benedettini et al. (2018) using getsources and this work using astrodendro. The red circle and blue cross symbols represent the cores by this work and Benedettini et al. (2018), respectively. The error bar was calculated by Poisson statistics.

thumbnail Fig. D.2.

Comparison of between astrodendro cores (this paper) and getsources cores (Benedettini et al. 2018). Histogram of core separation (left) and core mass versus separation diagrams in physical scale (right). The red circle and blue cross symbols represent the cores by this work and Benedettini et al. (2018), respectively. The error bar was calculated by Poisson statistics.

Appendix E: Dependence of astrodendro parameter

The astrodendro algorithm relies on three key parameters: min_value, min_delta, and min_npix. The third parameter, min_npix, is typically set to match the angular resolution of the map. To examine the influence of the first two parameters on core properties, we varied min_value and min_delta in the ranges of 3σ to 7σ, and 2σ to 4σ, respectively. The results are briefly summaries in Tab. E.1. Within these parameter ranges, we found that the total number of cores is sensitive to min_delta, while not to min_value. The mean core mass and core separation are also sensitive to min_delta, but there are no significant differences. Therefore, we adopted min_value = 5σ and min_delta = 3σ as the representative values in this paper.

Table E.1.

Physical properties for various dendrogram parameters.

Appendix F: Nearest neighbour search and the core separation

A NNS is a method of finding the closest point to a given query point from a set of data points. As the input data, we used the core center locations in the 2D plane of sky. Using the tool in the matlab library, we computed the distance from the query core center to the nearest neighbour core center. We used the distance as a core separation of the query core. Ishihara et al. (2024) derived the correction factor to convert the 2D core separation to 3D value. Assuming that the observed cores are uniformly distributed in a spherical space, the 3D core separation is computed as

S 3 D = S obs ( 4 π ) 1 π / 2 π / 2 sin θ × 2 π sin θ d θ = S 2 D × 4 π $$ \begin{aligned} S_{\rm 3D} = \frac{S_{\rm obs}}{(4\pi )^{-1} \int ^{\pi /2}_{-\pi /2}\sin \theta \times 2\pi \sin \theta d \theta } = S_{\rm 2D} \times \frac{4}{\pi } \end{aligned} $$(F.1)

See Ishihara et al. (2024) for more details.

All Tables

Table 1.

Physical properties of target clouds.

Table 2.

Physical properties of identified cores.

Table E.1.

Physical properties for various dendrogram parameters.

All Figures

thumbnail Fig. 1.

Core mass versus radius diagram (left) and histograms of core mass (middle) and core radius (right) for Polaris and Lupus I by astrodendro. The red circle and blue cross symbols represent the cores in Polaris and Lupus I, respectively. The error bar was calculated by Poisson statistics.

In the text
thumbnail Fig. 2.

Histogram of core separation (left) and core mass versus separation diagrams in physical scale (middle) and normalized scale (right). The red circle and blue cross symbols represent the cores in Polaris and Lupus I, respectively. For the normalized scale, the core mass was divided by the local thermal Jeans mass (mJ) and the core separation was divided by the local thermal Jeans length (lJ), where mJ and lJ are estimated in the trunks that include the corresponding cores (leaves). The error bar was calculated by Poisson statistics.

In the text
thumbnail Fig. 3.

Histograms of mean volume density of identified core for Polaris (top) and Lupus I (bottom). The white and diagonal bins represent the sample of all cores and only bound cores, respectively. The solid red curves and dotted red line indicate the results of the log-normal fit and power-law fit. The error bar was calculated by Poisson statistics.

In the text
thumbnail Fig. A.1.

Column density distributions of target clouds. The H2 column density maps of Polaris (top) and Lupus I (bottom). Red star and black cross symbols are the positions of the protostellar cores identified by (Benedettini et al. 2018) and those of the dense cores identified by astrodendro (this work).

In the text
thumbnail Fig. D.1.

Comparison of between astrodendro cores (this paper) and getsources cores (Benedettini et al. 2018). Core mass versus radius diagram (left) and histograms of core mass (middle) and core radius (right) identified by Benedettini et al. (2018) using getsources and this work using astrodendro. The red circle and blue cross symbols represent the cores by this work and Benedettini et al. (2018), respectively. The error bar was calculated by Poisson statistics.

In the text
thumbnail Fig. D.2.

Comparison of between astrodendro cores (this paper) and getsources cores (Benedettini et al. 2018). Histogram of core separation (left) and core mass versus separation diagrams in physical scale (right). The red circle and blue cross symbols represent the cores by this work and Benedettini et al. (2018), respectively. The error bar was calculated by Poisson statistics.

In the text

Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.