On the fidelity of the core mass functions derived from dust column density data

Aims: We examine the recoverability and completeness limits of the dense core mass functions (CMFs) derived for a molecular cloud using extinction data and a core identification scheme based on two-dimensional thresholding. Methods: We performed simulations where a population of artificial cores was embedded into the variable background extinction field of the Pipe nebula. We extracted the cores from the simulated extinction maps, constructed the CMFs, and compared them to the input CMFs. The simulations were repeated using a variety of extraction parameters and several core populations with differing input mass functions and differing degrees of crowding. Results: The fidelity of the observed CMF depends on the parameters selected for the core extraction algorithm for our background. More importantly, it depends on how crowded the core population is. We find that the observed CMF recovers the true CMF reliably when the mean separation of cores is larger than their mean diameter (f>1). If this condition holds, the derived CMF is accurate and complete above M > 0.8-1.5 Msun, depending on the parameters used for the core extraction. In the simulations, the best fidelity was achieved with the detection threshold of 1 or 2 times the rms-noise of the extinction data, and with the contour level spacings of 3 times the rms-noise. Choosing larger threshold and wider level spacings increases the limiting mass. The simulations show that when f>1.5, the masses of individual cores are recovered with a typical uncertainty of 25-30 %. When f=1 the uncertainty is 60 %. In very crowded cases where f<1 the core identification algorithm is unable to recover the masses of the cores adequately. For the cores of the Pipe nebula f2.0 and therefore the use of the method in that region is justified.

💡 Research Summary

This paper investigates how reliably dense core mass functions (CMFs) can be derived from dust extinction maps when a two‑dimensional threshold‑based core identification algorithm is applied. The authors embed synthetic cores into the real, spatially variable extinction background of the Pipe Nebula and then extract those cores using a range of algorithmic parameters. By comparing the recovered CMFs with the known input CMFs, they quantify both the fidelity (how closely the shape matches) and the completeness (the lowest mass at which the CMF is reliably sampled).

A key variable introduced is the crowding factor f, defined as the ratio of the mean separation between cores to their mean diameter. When f > 1, cores are on average more isolated than they are large, and the algorithm can separate them cleanly. In this regime the recovered CMF reproduces the true CMF above a limiting mass of roughly 0.8–1.5 M⊙, depending on the detection threshold and contour spacing used. The simulations show that the best performance occurs with a detection threshold set to 1–2 times the rms‑noise of the extinction map and contour level spacings of 3 times the rms‑noise. Raising the threshold or widening the level spacing pushes the completeness limit to higher masses and reduces the number of low‑mass cores that are detected.

When f > 1.5, individual core masses are recovered with a typical uncertainty of 25–30 %. At f ≈ 1 the uncertainty grows to about 60 %, and for f < 1 (highly crowded fields) the algorithm fails to retrieve reliable masses altogether because cores blend into each other. The Pipe Nebula itself has f ≈ 2.0, placing it comfortably in the regime where the method is justified.

The study underscores that the choice of extraction parameters cannot be made in isolation; it must be matched to the underlying crowding of the core population. In regions where cores are densely packed, more sophisticated three‑dimensional techniques or advanced segmentation methods (e.g., machine‑learning approaches) will be required to avoid systematic biases in the CMF. Moreover, the completeness limit is directly tied to the noise level of the extinction data, suggesting that deeper, lower‑noise observations will extend the reliable mass range downward.

Overall, the paper provides a practical framework for assessing the reliability of CMFs derived from extinction data, offering clear guidelines on parameter selection and highlighting the critical role of core crowding. These insights are valuable for any study that uses CMFs to infer star‑formation efficiencies, initial mass functions, or the physical conditions within molecular clouds.