EDP Sciences logo
Subscriber Authentication Point
Search
Menu
Home All issues Volume 693 (January 2025) A&A, 693 (2025) A51 Full HTML
Open Access
Issue
A&A
Volume 693, January 2025
Article Number A51
Number of page(s) 12
Section Galactic structure, stellar clusters and populations
DOI https://doi.org/10.1051/0004-6361/202450434
Published online 03 January 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.

Open Access funding provided by Max Planck Society.

1 Introduction and motivation

Understanding the rate at which molecular clouds convert their gas reservoirs into stars is of particular importance for many branches of astrophysics. The star formation rate (SFR) determines a wide range of the observed properties of galaxies and affects their chemical and dynamical evolution (Helou et al. 2000; Dib et al. 2009; Maraston et al. 2010; Dib 2011; Dariush et al. 2016; Pacifici et al. 2016; Barrera-Ballesteros et al. 2021; Grisdale et al. 2022; Garduno et al. 2023; Egorov et al. 2023). Within individual molecular clouds, the SFR is regulated by a number of physical processes that include the clouds’ self-gravity, supersonic turbulence and the way it is driven, their chemical composition, the incident radiation field, and tidal fields (Krumholz & McKee 2005; Dib et al. 2007, 2008, 2010, 2020; Padoan & Nordlund 2011; Glover & Clark 2012; Federrath & Klessen 2012; Schneider et al. 2022; Rani et al. 2022; Dib 2023a, Li 2024; Zhou et al. 2024a,b). On large scales, the SFR can be regulated by other processes, such as the gravity due to the existing stellar populations (Dib et al. 2017a; Shi et al. 2018; Marchuk 2018; Pessa et al. 2022), galactic shear (Seigar 2005; Dib et al. 2012; Aouad et al. 2020), and the overall galactic potential (Jog 2014; Meidt et al. 2018).

The SFR can be measured using a variety of star-formation indicators. Some of these methods rely on the stellar light that is directly emitted by young stars, such as in the ultraviolet, or by dust-processed stellar light in the infrared. For a detail review of these methods, we refer the reader to Calzetti (2013). While it should be noted that these methods are mostly employed for external galaxies where stellar populations are not well resolved, they make assumptions about the shape of the stellar IMF, and commonly assume that it is universal within and across galaxies. Recent work has shown that this is not the case, either on individual molecular clouds scales (Dib 2014; Dib et al. 2017b; Dib 2023b) or for entire galaxies (Dib & Basu 2018; Dib 2022). In nearby star forming regions in the Milky Way, a different method is usually used to infer the SFR. This method is based on the direct counting of young stellar objects (YSOs)1, which are optimally identified in the infrared (Forbrich 2009; Enoch et al. 2009; Evans et al. 2009; Heiderman et al. 2010; Gutermuth et al. 2011; Ybarra et al. 2013; Jose et al. 2013; Strafella et al. 2015; Stutz & Kainulainen 2015; Jose et al. 2016; Fischer et al. 2016; Hasan et al. 2023; Tobin & Sheehan 2024). The number of YSOs in a given evolutionary stage in a star-forming region, Ni, is converted into a SFRi using SFRi=Ni⟨M⟩τi,$\[\mathrm{SFR}_{i}=N_{i} \frac{\langle M\rangle}{\tau_{i}},\]$(1)

where the index i runs over the different protostellar classes (i=0, I, F, II, and III for objects in Class 0, Class I, Class F, Class II, and Class III, respectively), ⟨M⟩ is the mean protostellar mass present in the region, and τi is the typical lifetime of the YSOs spent in each class. The accuracy of the SFR measurements using Eq. (1) depends on the validity of some of the assumptions that are adopted when calculating the average mass ⟨M⟩ and protostellar lifetimes. For example, for Class II YSOs, τII is often taken to be ≈2 Myr. The mean mass is usually assumed to be ⟨M⟩ = 0.5 M⊙, which corresponds to the mean mass in a Milky Way-like IMF. These types of simplifications could potentially generate a bias in the estimates of the SFR. Another source of uncertainty that can affect the SFR measurements using YSO counts is the incompleteness of the YSO censuses and the level of incompleteness can vary with the evolutionary stage of the YSOs as well as with the specific details of how the observational data were obtained and reduced (Gutermuth et al. 2009; Winston et al. 2010; Broekhoven-Fiene et al. 2014; Dunham et al. 2015; Gutermuth & Heyer 2015). In this work, we attempt to estimate whether or not such biases exist and to determine their amplitudes. We approach this problem by generating synthetic populations of N* YSOs that form in a molecular cloud with a prescribed star formation history (SFH). The YSOs have masses that are drawn from an IMF and are assigned protostellar lifetimes that span the realistic range of each class. We then measure the real SFR and compare it to the SFR derived using Eq. (1) for YSOs in different classes. In Sect. 2, we describe the elements of the model, namely the protostellar mass function from which stellar masses are drawn, the function that describes the SFH of the YSOs, and their lifetimes spent within each class. In Sect. 3, we present the results of the time evolution of the SFRs derived from the various protostellar classes and compare them to the real SFR, exploring both the effects of the SFH and of the IMF. In Sect. 4, we measure the offsets between the real SFRs and those derived using the Class I and Class II populations. We explore the correlation between these offsets and the ratios between the numbers of Class I and II YSOs and those of Class III. Finally, in Sect. 5, we present our conclusions.

2 The model

2.1 The protostellar mass function

As we do not account for the mass evolution of YSOs, which is minimal for low-mass stars, we also refer to the protostellar mass function as the stellar IMF and assume that it can be described with a tapered power-law function (TPL; Dib 2014; Dib et al. 2017b; Dib 2023b). The TPL is characterized by two power laws in the low- and high-mass regimes, and a characteristic mass (i.e., peak mass) and is given by ϕi(M)=ϕ0M−Γ−1{1−exp⁡[−(MMch)γ+Γ]},$\[\phi_{\mathrm{i}}(M)=\phi_{0} M^{-\Gamma-1}\left\{1-\exp \left[-\left(\frac{M}{M_{\mathrm{ch}}}\right)^{\gamma+\Gamma}\right]\right\},\]$(2)

where ϕ0 is a normalization term, Γ is the exponent of the power law in the high-mass end, γ is the exponent in the low-mass regime, and Mch is the characteristic mass. The Salpeter (1955) value for the slope in the high-mass regime corresponds to Γ = 1.35. The Milky Way values for the two other variables for a single-star IMF (i.e., corrected for binarity) are γ = 0.51 and Mch = 0.35 M⊙ (Parravano et al. 2011). In this work, we consider cases where the parameters of the clusters’ IMF are assigned Milky Way-like values, and others where we vary the three parameters within the ranges that are permitted by the observations of young Galactic stellar clusters. The minimum and maximum stellar masses we consider are 0.02 M⊙ and 50 M⊙, respectively.

2.2 Protostellar classes and associated lifetimes

The classification into Classes I, II, and III was initially proposed by Lada & Wilking (1984) and is based on the value of the slope of the spectral energy distribution (SED) in the wavelength range of 2 to 20 μm. This classification system was subsequently used in numerous studies (Andre & Montmerle 1994; Evans et al. 2009; Rebull et al. 2014; Furlan et al. 2016; Sharma et al. 2017; Pokhrel et al. 2020; Kuhn et al. 2021; Sun et al. 2022). The Class 0 classification was later introduced by André et al. (2000) and is based on the indirect inference of the presence of a YSO based on the detection of a compact submillimeter source and a collimated CO outflow. These protostellar classes are associated with different phases of the formation of stars and each phase is characterized by a specific shape of the SED. The star-formation processes taking place in these phases are as follows: In Class 0, a YSO is formed in the central region of a protostellar core with an envelope mass that is far in excess of the YSO mass; in Class I, there is collapse of the envelope onto the central object with the transition between Class 0 and Class I being the point in time at which the envelope mass and the mass of the protostar are almost equal; in Class II, a disk emerges around the central star; and in Class III, the disk dissipates by various processes, such as the formation of planets, photo-evaporation, and tidal stripping. An intermediate class between Class 0 and Class I has been proposed by Greene et al. (1994) and was labelled the Flat class (i.e., a relatively flat SED in the wavelength range of 2 to 20 μm) but it was shown that this class is close to Class I in terms of the evolutionary status of the YSOs as it is associated with signs of a collapsing envelope (Calvet et al. 1994). Recent work shows that genuine flat-spectrum YSOs are rare, once contamination by reddened Class II YSOs and dusty star-forming galaxies is corrected for (Pokhrel et al. 2023)2.

The exact duration of each evolutionary stage for YSO classes is relatively uncertain. Estimates of the mean duration of each of these stages depend on the number of objects found in each class in star-forming regions and this may be affected by misclassifications due for example to YSOs being seen edge-on (Gutermuth et al. 2009). The usual approach is to use Eq. (1) and the numbers of YSOs found in different classes in order to derive relative ages for each class. The assumption is made that YSOs are formed at a constant rate over time and that their formation occurs over timescales that are typically longer than the longest duration of any class. By constraining the age of one of the classes, such as that of the Class II T Tauri stars, which is assumed to be of the order of their disk lifetime3, one can then derive the corresponding duration of the other classes for each YSO. A number of studies measured these timescales in nearby star-forming clouds and found them to fall in the range of 0.005 to 0.17 Myr for Class 0 YSOs (Andre et al. 1993; Froebrich et al. 2006; Enoch et al. 2009; Fischer et al. 2017; Kristensen & Dunham 2018), between 0.1 and 0.7 for Class I YSOs (Wilking et al. 1989; Greene et al. 1994; Kenyon & Hartmann 1995; Hatchell et al. 2007), and between 0.4 and 5 Myr for Class II YSOs (Wilking et al. 1989; Kenyon & Hartmann 1995; Hernandez et al. 2008), and found ages in excess of 2 Myr for Class III YSOs (Wilking et al. 2005).

In the present work, we adopt the values of the lifetimes measured by Evans et al. (2009). These were derived from observations of several nearby molecular clouds in the Spitzer legacy project, “From Molecular Cores to Planet-forming Disks” (c2d). Evans et al. (2009) found that the mean lifetimes of Class 0, Class I, and the Flat Class are τ¯0=0.1 Myr,τ¯I=0.44 Myr$\[\bar{\tau}_{0}=0.1 ~\mathrm{Myr}, \bar{\tau}_{\mathrm{I}}=0.44 ~\mathrm{Myr}\]$, and τ¯F=0.35$\[\bar{\tau}_{\mathrm{F}}=0.35\]$ Myr. These were derived assuming a constant SFR and a lifetime for Class II of τ¯II=2 Myr$\[\bar{\tau}_{\mathrm{II}}=2 ~\mathrm{Myr}\]$ and include a correction of the photometry for extinction effects. The uncertainties on these lifetimes are quite large. For Class II, the 1σ uncertainty is about στII = 1 Myr. Given all the assumptions, such as that of a constant SFR, the potential confusion of a fraction of the YSOs with extragalactic sources, and differences in the relative numbers of YSOs found in each region of the c2d cloud sample, we adopt uncertainties on the lifetimes that are 50% of their adopted values, such that στ0 = 0.05 Myr, στI = 0.22 Myr, and στF = 0.175 Myr. These results were corroborated by the study of Hsieh & Lai (2013), who reanalyzed the c2d data and found similar results for the lifetimes to those obtained by Evans et al. (2009). For the Class III YSOs, their census in the c2d is incomplete, and therefore their associated lifetimes were not estimated. In general, the lifetime associated with Class III YSOs is uncertain. An approximation of this duration is the Kelvin-Helmotz timescale, which is the typical time during which a young star contracts and radiates its kinetic energy and descends the Hayashi track to reach the main sequence (Hayashi 1961). The Kelvin-Helmotz timescale is given by τKH ≈ (GM2)/(2RL), where G is the gravitational constant and M, R, and L are the star’s mass, radius, and stellar luminosity, respectively. For solar-mass stars, τKH ≈ 30 Myr, and this value is lower for higher-mass stars and is substantially higher for low-mass stars. However, as we follow the evolution of the young stellar populations for the first 10 Myr, the choice of any value of τKH that is larger than 30 Myr will make little difference to our calculations. Here, we consider a single value across the entire stellar mass range of τ¯III =40 Myr$\[\bar{\tau}_{\text {III }}=40 ~\mathrm{Myr}\]$ with στIII = 20 Myr.

In the present work, we draw the lifetimes of the population of YSOs from Gaussian distribution functions given by P(τi)=1στi2πexp⁡(−12(τi−τ¯iστi)2),$\[P(\tau_{i})=\frac{1}{\sigma_{\tau_{i}} \sqrt{2 \pi}} \exp \left(-\frac{1}{2}\left(\frac{\tau_{i}-\bar{\tau}_{i}}{\sigma_{\tau_{i}}}\right)^{2}\right),\]$(3)

where i = 0, I, F, II, and III correspond to Class 0, Class I, the F Class, Class II, and Class III, respectively, and τ¯i$\[\bar{\tau}_{i}\]$ and στi are the characteristic lifetimes and the corresponding 1σ uncertainties for each class taken from Evans et al. (2009).

2.3 The star formation history

When deriving the SFR using Eq. (1), and in almost all observational studies, the assumption is made that the SFH in a star-forming region has been constant over the entire lifetime of the regions (Evans et al. 2009 and references therein). This assumption may be true in some molecular clouds, possibly in low-mass clouds where star formation proceeds at a slow pace, but is unlikely to be true for all clouds. The age distribution of stars in young clusters is not always a flat distribution as exemplified by the case of the Orion Nebula Cluster (ONC, Dib et al. 2013). In the ONC, the age distribution of stars is well reproduced by the sum of a linearly increasing SFH at old ages (up to 6 Myr) and an acceleration (i.e., a burst) of star formation at younger ages (see Figure 13 in Dib et al. 2013). Similar SFHs, consisting of a slow increase when clouds are in the process of assembling followed by a burst of star formation at later times as most of the gas sinks to the central regions of the clouds, are also observed in numerical simulations (e.g., Vazquez-Semadeni et al. 2017, Guszejnov et al. 2022). It has also been shown that such burst-like SFHs better describe the SFH of clusters such as the Orion Nebular cluster (Kounkel 2020; see also Megeath et al. 2022), and are necessary in order to reproduce the mass–luminosity relation of massive star-forming regions (Zhou et al. 2024c,d). In the present work, we test both assumptions. We consider cases where the SFH is constant and others where it varies over time. For cases with a constant SFH (labelled as constant or flat), we consider a fiducial case with a birthrate of 5 × 10−5 stars yr−1 (i.e., 500 stars form over a period of 107 yr). For the time-dependent SFH, a variety of functional forms can be used to describe a slowly increasing SFH followed by a burst. For simplicity, we use a Gaussian function, for which we are able to vary the position of the peak and the width of the distribution. The SFH is described by a probability for the time of birth of the stars, tb, and in the case of a time-varying SFH this is given by P(tb)=1σt2πexp⁡(−12(tb−tb,pσt)2),$\[P(t_{b})=\frac{1}{\sigma_{\mathrm{t}} \sqrt{2 \pi}} \exp \left(-\frac{1}{2}\left(\frac{t_{\mathrm{b}}-t_{\mathrm{b}, \mathrm{p}}}{\sigma_{t}}\right)^{2}\right),\]$(4)

where tb,p is the position of the peak and σt the standard deviation. For a cluster whose stars form at a constant rate, their tb are randomly drawn from a flat probability between t = 0 and tend, where tend is the time at which star formation ends. In the case of a time-dependent SFH, the tb of stars are randomly drawn from the P(tb) function given in Eq. (4). For the case with a time-dependent SFH, the fiducial parameters are tb,p = 2 Myr and σt = 1 Myr. We take a tend that is larger than any of the lifetimes associated with the protostellar classes and assume that tend = 10 Myr.

thumbnail Fig. 1

Four realizations of the Milky Way-like protostellar mass function. The masses of N* = 500 protostars are randomly drawn from a Milky Way-like stellar mass function in the stellar mass range of 0.02–50 M⊙.

3 Results

3.1 The effect of the SFH

First, we present results from two fiducial cases, one with a constant SFH and another with a Gaussian-like SFH. For these models, we consider a cluster that will end up containing N* = 500 stars sampled from a Milky Way-like single-star IMF (i.e., Γ = 1.35, γ = 0.51, and Mch = 0.35 M⊙). For each case, we generate 250 models. Figure 1 displays four realizations of the Milky Way-like randomly sampled IMF. Figures 2 and 3 each display four realizations of the birth times of protostars. For the case of the flat SFH displayed in Fig. 2, the birth times of protostars are randomly selected with a uniform probability in the time span between t = 0 and t = tend, whereas for the cases of the Gaussian-like SFH displayed in Fig. 3, the birth times of protostars are randomly sampled using Eq. (4). The next step is to assign to each protostar lifetimes that are associated with each protostellar class. As described above, these lifetimes are randomly drawn using the probability distribution functions described by Eq. (3). As fiducial values, we adopt here those listed in Sect. 2.2, namely τ¯0=0.1 Myr$\[\bar{\tau}_{0}=0.1 ~\mathrm{Myr}\]$ and στ0 = 0.05 Myr for Class 0, τ¯I=0.44 Myr$\[\bar{\tau}_{\mathrm{I}}=0.44 ~\mathrm{Myr}\]$ and στI = 0.22 Myr for Class I, τ¯F=0.35 Myr$\[\bar{\tau}_{\mathrm{F}}=0.35 ~\mathrm{Myr}\]$ and στF = 0.175 Myr for the F Class, τ¯II=2 Myr$\[\bar{\tau}_{\mathrm{II}}=2 ~\mathrm{Myr}\]$ and στII = 1 Myr for Class II, and τ¯III=4 Myr$\[\bar{\tau}_{\mathrm{III}}=4 ~\mathrm{Myr}\]$ and στIII = 2 Myr for Class III. As an example, Fig. 3 displays the distribution functions of the protostellar lifetimes for all YSOs that form in one cluster. When randomly sampling each YSO lifetime from the Gaussian function shown in Eq. (3), the randomly sampled values are drawn in the range of the ±1σ uncertainty of each quantity.

With the above input quantities, we track the number of YSOs found in each class as a function of time. Figure 4 (left panel) displays four examples (out of 250 generated models) for the case with a flat SFH, a Milky-like IMF, and N* = 500 in each cluster. Figure 4 indicates that within the temporal fluctuations due to the random sampling, the total number of Class 0 protostars is constant over time. The number of Class I and Class F protostars starts to increase after ≈0.1 and 0.4 Myr, respectively, and both reach a plateau after 1 Myr. The numbers of Class II and Class III PMS stars display the same behavior but at later times (i.e., after between two and several million years). The right panel of Fig. 4 displays the time evolution of the number of YSOs in each class for four cases with a Gaussian-like SFH (with the fiducial parameters). This figure clearly shows the peak in Class 0 protostars at t = tb,p = 2 Myr and this peak shifts to later epochs for more evolved protostars. WE note that in these cases, the number of Class II and Class III PMS stars also declines at later times because most of the protostars were formed early at t = 2 Myr and almost no new protostars were formed after ≈4 Myr.

Using the number of YSOs found in each class and at each epoch, as well as the mean values of their lifetimes (i.e., those suggested by Evans et al. 2009), we can now apply Eq. (1) and calculate the corresponding SFRs. We also measure the true SFR in each model. For this, we take a time step of Δt = 104 yr. This guaranties that the smallest protostellar lifetimes are well resolved. We then simply calculated the true SFR as the total mass of protostars that form within each time step: SFRreal⁡(t+Δt/2)=∑jM∗,jΔt,$\[\operatorname{SFR}_{\text {real}}(t+\Delta t / 2)=\frac{\sum_{j} M_{*, j}}{\Delta t},\]$(5)

where j runs over all stars that form between t and t + Δt. Figure 5 (left) displays the time evolution of the SFR for four cases where the SFH is flat. We recall that this is usually the assumption that is made in the observations when deriving the SFR. Figure 6 (left panel) displays the true SFR measured directly from the models (dashed lines) and an SFR estimate using the total number of YSOs (full line). For the latter, a characteristic YSO lifetime of 2.5 Myr was used (Pokhrel et al. 2020). The left panels of Figs. 5 and 6 show that, while the SFR estimates measured using a single class converge to the real value, they do so with a delay that increases for more evolved YSOs. This implies that for cases where the SFH is constant, measurements of the SFR will only be reliable if only Class 0 and I protostars are considered in young star forming regions and Class II PMS stars in older regions (i.e., ages ≳2 Myr). Figure 6 (left panel) additionally shows that using the total number of YSOs provides a very poor approximation of the real SFR; it underestimates the true SFR by a factor of between ≈10 and 100 in young star-forming regions and overestimates the true SFR by a factor of a few to several for older regions. The right panels in Figs. 5 and 6 display the same SFRs calculated for cases where the SFH is a Gaussian-like function with the fiducial values of the parameters, tp,b = 2 Myr and σt = 1 Myr. The same patterns that are observed with a constant SFH are also prevalent here, with the additional effect that at advanced epochs (i.e., ages ≳3.5 Myr), the estimates of the SFR that are either based on the Class II and Class III PMS stars or on the total number of YSOs overestimate the true SFRs by several orders of magnitude. This is due to the fact that while star formation has ceased, Class II and Class III PMS stars are long lived and their numbers bias the measurement of the SFR.

In order to get a clearer picture of the offset between the SFRs estimated using the YSO populations and the true SFR, we show the ratio between these quantities in Fig. 7. In both panels of this figure, the ratios were obtained by averaging over the 250 realizations (i.e., 250 clusters), each with N* = 500 and a Milky Way-like IMF. The top panel corresponds to cases with a flat SFH and the bottom panel to cases with the Gaussian-like SFH, with all clusters in the latter case having tb,p = 2 Myr and σt = 1 Myr. The light dashed line in both panels corresponds to the 1σ Poisson uncertainty on each measurement. Ideally, for these measurements to be accurate, the offset parameter ηi=0,I,F,II,II = SFRi/SFRreal would have to be close to the unity. However, this is far from being the case, with the worst departure from unity corresponding to the time-dependent, Gaussian-like SFH, especially around and after the peak of star formation at ≈2 Myr. For these time-dependent models of the SFH, we expect the level of departure from unity to depend on the duration of the burst (i.e., the width of the Gaussian). We therefore ran additional models, each with 250 clusters with Gaussian-like SFHs. As our fiducial case has a value of σt = 1 Myr, we considered additional cases where σt = 0.5 Myr (a narrow, high-amplitude burst) and σt = 2 Myr (an extended, low-amplitude burst). The results for these models are displayed in Fig. 8 for the narrow burst (top panel) and the extended burst (bottom panel). Qualitatively, both cases display the same features and are similar to the fiducial case displayed in Fig. 7 (top panel). However, a noticeable difference for the narrow-burst case (i.e., case with σt = 0.5 Myr) is that all SFR measurements using the different classes of YSOs or their total number, with the exception of the one based on the Class 0 population, fail to reproduce the true SFR.

The conclusion that can be drawn from Figs. 7 and 8 is that the only reliable estimate of the SFR at all epochs is the one based on the population of Class 0 protostars. The latter measurement is always accurate, irrespective of the SFH. Measurements of the SFR based on the population of more evolved YSOs can provide reliable estimates of the SFR in cases where the SFH is constant or displays a weak burst, and this is only valid when the star-forming regions are themselves more evolved (i.e., with ages of 2 Myr or older, and harboring a significant fraction of evolved YSOs). In all circumstances, estimates of the SFR that rely on the total number of YSOs and a characteristic YSO lifetime (in this work, 2.5 Myr) are not accurate and this is true irrespective of the SFH and the considered age of the region. Finally, we should mention that we ran models with lower and higher numbers of YSOs, namely N* = 250 and 1000, respectively, while keeping tend = 10 Myr in all models. Models with lower or higher numbers of protostars display the same effects as those observed in models with N* = 500 YSOs, but with an increased and decreased level of temporal fluctuation, respectively (see Appendix A).

thumbnail Fig. 2

Top: four realizations of the birth time of protostars with a flat, constant SFH. Bottom: four realizations of the birth time of protostars with a Gaussian SFH. The peak and standard deviation of the Gaussian are set to tb,p = 2 Myr and σt = 1 Myr. For improved visibility, the values of tb are binned with a bin size of 2 × 105 yr.
thumbnail Fig. 3

One realization of the distribution functions of the lifetimes associated with each protostellar class. The lifetimes are drawn using Eq. (4) with the mean value and standard deviation associated with each lifetime (see text for details). The bin size for each protostellar lifetime distribution is adjusted to the mean value of that distribution, namely the bin sizes are 5×103 yr, 2.2 × 104 yr, 1.75 × 104 yr, 105 yr, and 2 × 106 yr for the Class 0, Class I, the F Class, Class II, and Class III, respectively.
thumbnail Fig. 4

Left: four realizations of the time evolution of the number of YSOs in different protostellar classes. Each cluster has N* = 500, a Milky Way-like IMF and a constant SFH. Right: similar to left panel, but here each cluster has N* = 500, a Milky Way-like IMF, and a Gaussian-like SFH. The Gaussian’s peak is at 2 Myr and the standard deviation is 1 Myr. The bin-size on the x-axis is 105 yr and the same value is used in subsequent plots.
thumbnail Fig. 5

Left: four realizations of the time evolution of the SFR measured using the populations of YSOs found in different classes and their respective lifetimes. The dashed line displays the time evolution of the true SFR. Each cluster has N* = 500, a Milky Way-like IMF, and a constant SFH. Right: similar to the left panel, but here each cluster has N* = 500, a Milky Way-like IMF, and a Gaussian-like SFH. The Gaussian’s peak is at 2 Myr and the standard deviation is 1 Myr.
thumbnail Fig. 6

Left: four realizations of the time evolution of the true SFR (SFRreal, dashed line) and the SFR measured using the total number of YSOs and an average lifetime for all Classes. Each clusters has N* = 500, a Milky Way-like IMF and a constant SFH. Right: similar to the left panel, but here, each cluster has N* = 500, a Milky Way-like IMF, and a Gaussian-like SFH. The Gaussian’s peak is at 2 Myr and the standard deviation is 1 Myr.

3.2 The effect of the protostellar mass function

Recent studies suggest a non-negligible degree of variation in the parameters that characterize the shape of the IMF of young clusters in the Milky Way and in M31 (Dib 2014; Weisz et al. 2015; Lim et al. 2015; Dib et al. 2017a, Dib 2023b). Dib et al. (2017a) found that an intrinsic scatter in the IMF parameters is necessary in order to match the observed fraction of single O stars in young Galactic clusters (cluster ages ≲12 Myr). For a TPL function representing the system IMF, Dib et al. (2017a) found that the standard deviation in Γ, γ, and Mch is ≈0.6, 0.25, and 0.27 M⊙. Here, as we are dealing with a single-star protostellar mass function, we model the distribution of each of the mass function parameters with a Gaussian centered around the Galactic values (Γ = 1.35, γ = 0.51, and Mch = 0.35 M⊙) and with a standard deviation on each parameter that is 50% of its mean value, namely σΓ = 0.6, σγ = 0.25, and σMch = 0.17 M⊙. Similar to the previous calculations, we generate a sample of 250 clusters, each with N* = 500, and where the parameters of the IMF of each cluster are randomly drawn from these Gaussian distributions. The IMF parameters are randomly drawn from their Gaussian distribution within the ±1σ uncertainty range. Figure 9 displays the cases for four clusters with such a variable IMF. We then follow the same procedure as above and derive the SFRs based on the different YSO classes and an SFR based on the total YSO population. Adopting the case of a flat SFH, Fig. 10 displays the time evolution of these SFRs normalized by the true SFR. Comparing this figure with Fig. 7 (top panel), we conclude that a variable IMF has little effect on the derived SFRs. The same features observed in the top panel of Fig. 7 and where the IMF of clusters are randomly drawn from a Milky Way-like IMF are seen in Fig. 10; that is, we see a faithful reproduction of the true SFR (SFRreal) by the SFR based on the Class 0 population, and a reproduction of the SFRreal by the SFRs based on the more evolved populations (Class I and beyond) at more advanced epochs. The SFR derived using the total YSO population remains a poor approximation of the true SFR, as it severely underestimates SFRreal at early epochs and overestimates it at later epochs.

thumbnail Fig. 7

Top: time evolution of the ratio of the SFR measured using the populations of YSOs found in different classes (full lines) and the total population (dot-dashed line) and the real SFR that is generated in the models. Each line is an average over 250 clusters. Each cluster has N* = 500, a Milky Way-like IMF, and a flat SFH. Bottom: similar to the top panel, but here each clusters has N* = 500, a Milky Way-like IMF, and a Gaussian-like SFH. The Gaussian’s peak is at 2 Myr and the standard deviation is 1 Myr. In both panels, the dashed lines mark the 1σ uncertainty levels.
thumbnail Fig. 8

Top: time evolution of the ratio of the SFR measured using the populations of YSOs found in different classes (full lines) and the total population (dashed line) and the real SFR that is generated in the models. Each line is an average over 250 clusters. Each cluster has N* = 500, a Milky Way-like IMF, and a Gaussian-like SFH with the peak of the Gaussian located at tb,p = 2 Myr and a standard deviation of σt = 0.5 Myr. Bottom: similar to the top panel, but with tb,p = 2 Myr and σt = 2 Myr. For improved clarity, the 1σ uncertainties are not shown.
thumbnail Fig. 9

Four realizations of the variable protostellar mass function. The masses of N* = 500 protostars are randomly drawn from a TPL function in the stellar mass range of 0.02 to 50 M⊙. For each realization, the parameters that characterize the shape of the TPL mass function are randomly drawn from Gaussian distribution functions (see text for more details).
thumbnail Fig. 10

Same as the top panel of Fig. 7 but in this case the parameters that describe the shape of the IMF for each of the 250 clusters are randomly drawn from Gaussian distribution functions (see text for more detail). For more clarity, the 1σ uncertainties are not shown.

4 Constraining the SFH using protostellar counts and SFR estimates

In Sect. 3, we argue that the SFR derived using the Class 0 YSO count provides a reliable estimate of the true SFR in star-forming regions, irrespective of the SFH. However, as exemplified by the models with N* = 500 stars, the number of Class 0 YSOs remains quite modest in most, if not all, epochs. Furthermore, Class 0 YSOs are sometimes difficult to detect as they are deeply embedded in the clouds (Evans et al. 2009). In contrast, the numbers of evolved YSOs is much larger. There are approximately 20 times more Class II YSOs than Class 0 YSOs. For this reason, using Class II objects to derive the SFR can be an appealing option in order to make use of the larger number of Class II objects. A caveat here is that the offset between the Class II-derived SFR (SFRII) and the true SFR depends on the SFH. As shown in Sect. 3, for the case of a flat SFH, the offset between SFRII and SFRreal exists only at early epochs (≲2 Myr from the start of star formation). For the case of a burst-like SFH, the situation is more complex, and an offset between SFRII and SFRreal exists at all times. It is possible to reverse the arguments and assumes that the observed SFRII corresponds to the SFRreal and verify whether this assumption holds for the predictions of a constant SFH or a burst-like one. Figure 11 displays the Class I and Class II SFR offsets plotted against the ratios (NI/NIII) and (NII/NIII), respectively, for the case with a flat SFH. Figure 12 is similar to Fig. 11 but for the fiducial case with a Gaussian-like SFH. The dependencies between the offsets ηI and ηII and the number ratios are noticeably different between these cases with a different SFH.

We can test the usefulness of these relations in discriminating between SFHs by including observational data. We use the data presented in Tables 3 and 5 of Evans et al. (2009) for five nearby molecular clouds, namely Cha II, Lupus, Perseus, Serpens, and Ophiucus. Table 5 of Evans et al. (2009) lists the numbers of YSOs in Class I, the F Class, Class II, and Class III. As the numbers of Class II YSOs are larger than those in all other classes, this indicates that star formation started at least ≈1 Myr ago, which is roughly equal to τ¯0+τ¯I+τ¯F$\[\bar{\tau}_{0}+\bar{\tau}_{\mathrm{I}}+\bar{\tau}_{F}\]$ plus a fraction of the Class II lifetime. Table 3 of Evans et al. (2009) also lists SFR estimates for these clouds. These SFR measurements are derived using the Class II YSO populations and were presented in separate works, but are all based on c2d observations. Using the number counts of Class I and Class II YSOs, we use Eq. (1) to calculate SFRI and SFRII and consequently ηI and ηII. If the SFR values quoted for these clouds are the true values, the ratio ηII for these clouds should coincide with the prediction of the model with a constant SFR. The results shown in Fig. 11 (bottom panel) confirm that this is indeed the case. Furthermore, the upper panel of Fig. 11 shows that the ηI versus (NI/NIII) relation for these five clouds is also reasonably well matched by the model with a constant SFH. In contrast, the model with a burst-like SFH does not overlap with the observational data, as can be seen in the lower and upper panels of Fig. 12. The conclusion that can be drawn from Figs. 11 and 12 is that the SFH in these five nearby and low- to intermediate-mass clouds (from ≈400 to ≈5000 M⊙; Table 1 in Evans et al. 2009 and references therein) is relatively constant. It would be interesting to include data points for more massive star-forming regions and regions that are located further away in the Galaxy.

Below, we list the fits to the ηI − (NI/NIII) and ηII − (NII/NIII) relations. For the flat SFH case, we find: ηI=SFRISFRreal =1.66(±0.017)(NINIII)0.026±0.005ηII=SFRIISFRreal =1.59(±0.011)(NIINIII)−0.063±0.006,$\[\begin{aligned}\eta_{\mathrm{I}} & =\frac{\mathrm{SFR}_{\mathrm{I}}}{\mathrm{SFR}_{\text {real }}}=1.66( \pm 0.017)\left(\frac{N_{\mathrm{I}}}{N_{\mathrm{III}}}\right)^{0.026 \pm 0.005} \\\eta_{\mathrm{II}} & =\frac{\mathrm{SFR}_{\mathrm{II}}}{\mathrm{SFR}_{\text {real }}}=1.59( \pm 0.011)\left(\frac{N_{\mathrm{II}}}{N_{\mathrm{III}}}\right)^{-0.063 \pm 0.006},\end{aligned}\]$(6)

and for the time-dependent, Gaussian-like SFH, we find: ηI=SFRISFRreal =2.33(±0.029)(NINIII)−0.13±0.005ηII=SFRIISFRreal =27.25(±0.68)(NIINIII)−0.98±0.01.$\[\begin{aligned}\eta_{\mathrm{I}} & =\frac{\mathrm{SFR}_{\mathrm{I}}}{\mathrm{SFR}_{\text {real }}}=2.33( \pm 0.029)\left(\frac{N_{\mathrm{I}}}{N_{\mathrm{III}}}\right)^{-0.13 \pm 0.005} \\\eta_{\mathrm{II}} & =\frac{\mathrm{SFR}_{\mathrm{II}}}{\mathrm{SFR}_{\text {real }}}=27.25( \pm 0.68)\left(\frac{N_{\mathrm{II}}}{N_{\mathrm{III}}}\right)^{-0.98 \pm 0.01} .\end{aligned}\]$(7)

As pointed out in Sect. 2.2, the Class III YSOs observed by Spitzer for these nearby star-forming regions are incomplete and Eqs. (6) and (7) should be revised with more complete class III censuses. More complete censuses of Class III YSOs could be provided by X-ray observations, as shown by Preibisch et al. (2005) in the Orion Nebular Cluster. However, we note that more accurately accounting for Class III YSOs will reinforce our conclusions concerning the comparison of our models with the observations. Accounting for a higher number of missing Class III YSOs would cause the observational number ratios (NI/NIII) and (NII/NIII) in Figs. 11 and 12 to move to the left. This would not affect the agreement between the flat SFH models and the observations, but it would cause burst-like models to deviate further from the observations. While it is not yet possible to include a more massive region, such as the entire Orion cloud, because the census of Class II YSOs is incomplete (Furlan et al. 20164). Furlan et al. (2016) found that the number of Class 0 YSOs in Orion (N0 = 92). Using Eq. (1), we conclude that the true SFR for Orion is 4.6 × 10−4 M⊙ yr−1. This value sits in the middle between the SFRs measured by Lada et al. (2010) for the Orion A and Orion B clouds.

thumbnail Fig. 11

Top panel: offset between the SFR derived from the Class I population and the true SFR as a function of the ratio between the numbers of Class I and Class III PMS stars. The figure contains data from 250 clusters, each containing N* = 500, and for each cluster, it includes measurements at 100 epochs. The SFH is flat and the IMF resembles that of the Milky Way field. Bottom panel: similar to the top panel but showing the offset between the SFR derived from the Class II population and the true SFR as a function of the ratio of the numbers of Class II and Class III PMS stars.
thumbnail Fig. 12

Same as Fig. 11 but in this case the SFH is a time-dependent Gaussian-like function with the peak of the Gaussian located at 2 Myr and a standard deviation of σt = 1 Myr.

5 Conclusions

In this work, we explore the circumstances under which counting YSOs in star-forming regions can yield reliable estimates of the SFR. To this end, we developed a Monte Carlo model in which the masses, the birth times of the YSOs, and the lifetimes associated with the different YSO classes (Classes 0, I, F, II, and III) are all randomly drawn from distribution functions that describe each of these quantities. The masses of the protostars are randomly drawn from an IMF that can be either similar to that of the Milky Way field or variable in the range of variations observed among young clusters in the Galaxy. The birth times of protostars are either drawn from a flat distribution function (i.e., constant SFH), or from a time-dependent, burst-like function. Finally, the lifetimes of the YSOs associated with each protostellar class are randomly drawn from Gaussian distribution functions centered around those suggested by the observations (Evans et al. 2009). Using these prescriptions, we follow the time evolution of the number of YSOs in different evolutionary classes and calculate the corresponding SFR at every epoch. We find that the SFRs derived using the Class 0 population reproduce the true SFR, and this conclusion is valid irrespective of the shape of the SFH. For a constant SFH, the SFRs derived using the more evolved populations of YSOs (Class I, Class F, Class II, and Class III) reproduce the real SFR at later epochs. For example, for an SFR estimate that is based on the Class II population and lifetimes, the real SFR is reproduced after ≳2 Myr from the beginning of star formation. For a time-dependent burst-like SFH, all SFR estimates based on the number counts of the more evolved populations (all classes but Class 0) fail to reproduce the true SFR. Also, we find that SFR estimates that are based on the total number of YSOs associated with a characteristic lifetime for all populations fail to reproduce the real SFR, irrespective of the shape of the SFH. Furthermore, we show that all of these conclusions are independent of the shape of the stellar IMF.

The synthetic models presented in this work can help shed light on the SFH of observed star-forming regions. We show that the models with a constant SFH make different predictions for the SFR offsets versus the number ratios of YSOs from those with a burst-like SFH. For five, nearby, low-mass starforming regions that constitute the Evans et al. (2009) sample (i.e, the Cha II, Lupus, Perseus, Serpens, and Ophiucus clouds), we show that the relations between the SFR offset parameters ηI = SFRI and ηII = SFRII of the clouds and the number ratios NI/NIII and NII/NIII, respectively, are matched by models with a constant SFR. A comparison between our models and a larger number of star-formation regions, such as those that would be observed by the ABYSS project (Kounkel et al. 2023) and with the James Webb Science Telescope (e.g., Lenkic et al. 2024; Peltonen et al. 2024) will enable us to distinguish between the SFHs of star-forming clouds as a function of some of their fundamental properties, such as mass and surface density. In this work, we have not accounted for the effects of mass accretion, which induces variability in the luminosity of YSOs and may lead to their misclassification. We believe however that such effects will not change our conclusion, as only a small fraction of the YSO population shows signs of variability (Mairs et al. 2017).

Deriving accurate SFRs is crucial for the interpretation of the scaling laws of star formation on individual molecular cloud scales and for testing both theories and numerical simulations of star formation (e.g., Hennebelle & Chabrier 2011; Dib 2011; Padoan & Nordlund 2011; Federrath & Klessen 2012; Burkhart 2018; Eden & Teyssier 2024). In particular, the SFR–surface gas density (Σg) relation on local scales is observed to be steeper than the one for entire galaxies or the one derived on approximately kiloparsec scales within galaxies (e.g., Heiderman et al. 2010; Hony et al. 2015; Pokhrel et al. 2020, 2021). One interpretation of this steep relation between Σg and the SFR is that the higher SFR observed in individual regions is due to the selection bias introduced by studying a region that is actively forming stars, while the galaxy-averaged values include regions that are forming stars at different rates (Kruijssen & Longmore 2014). Another interpretation is that the steep Σg–SFR relation for individual clouds is the reflection of an evolutionary sequence in which the SFRs of contracting clouds increase much faster than their mean densities or mean surface densities and this leads to a steep supra-linear Σg–SFR relation (Zamora-Avilés et al. 2012). Further work is needed in order to disentangle these effects from the offsets in the SFR estimates uncovered in the present work, which are solely due to biases introduced by the star counting method.

Acknowledgements

We thank the anonymous referee for useful comments and suggestions that helped clarify aspects of this work. SC acknowledges funding from the State Research Agency (AEI-MCINN) of the Spanish Ministry of Science and Innovation under the grant ‘Thick discs, relics of the infancy of galaxies’ with reference PID2020-113213GA-I00. M.A.L.L. acknowledge support from the Ramón y Cajal program funded by the Spanish Government (reference RYC2020-029354-I), and from the Spanish grant PID2021-123417OB-I00.

Appendix A The effect of the number of stars

In the paper, we have chosen N* = 500 as the fiducial number of stars that form in each cluster, both for cases with a constant SFH and a time-dependent, burst-like SFH. Here we show examples of the calculated SFR in each protostellar Class for clusters with N* = 250 and N* = 1000. Fig. A.1 displays four realizations with N* = 250 and N* = 1000 for cases with a constant SFR and a Milky Way-like IMF. Decreasing (or increasing) the numbers of YSOs that form in the cluster by a certain factor (here a factor of 2) decreases (increases) the SFR by the same factor but a smaller N* (or larger N*) increases (decreases) the temporal fluctuations in the SFR.

thumbnail Fig. A.1

top: Four realizations of the time evolution of the SFR measured using the populations of YSOs found in different classes and their respective lifetimes. The dashed line displays the time evolution of the true SFR. Each cluster has N* = 250, a Milky Way-like IMF and a constant SFH. bottom: Similar to the top panel, but here, each cluster has N* = 1000.

Appendix B The effect of the YSOs’ lifetime distribution functions

In the main text, we have adopted a width of the Gaussian distributions for each YSO lifetime that is 0.5 times the mean lifetime value. Here, we explore how the width of such distribution affects our results. We consider two cases where the width of the lifetimes’ Gaussian distributions is 0.25 and .0.75 their mean values. In both cases, we use a flat SFH, N* = 500, and adopt a Milky Way-like IMF. The results for the offset parameters are displayed in the top row. These figures show that when the distribution of YSO lifetimes are narrower, the saturated level of the offset ratio is reached earlier. This is simply due to the fact that narrower (wider) distributions imply that there are less (more) YSO that transition from a given class to the next one with lifetimes that are larger than the mean value of the corresponding lifetime. The effects of the YSOs lifetime distribution is reflected in the variations that can be observed in the offset vs. number ratios diagram. A comparison of the lower panels in Fig. A.2 along with Fig. 11 shows that these diagrams are affected by the level of variations in YSO lifetimes for each class, even though at this stage, it is difficult to distinguish which set of models fits the data best given the limited number of observational data.

thumbnail Fig. A.2

Top row: same as Fig. 7 (top panel) but for cases where the width of the Gaussian distribution of the YSO lifetimes are taken as 0.25 the mean value (left) and 0.75 the mean value (right). Bottom row: Similar to Fig. 11 but for cases where the width of the Gaussian distribution of the YSO lifetimes are taken as 0.25 the mean value (left) and 0.75 the mean value (right).

References

  1. André, P., & Montmerle, T. 1994, ApJ, 420, 837 [Google Scholar]
  2. André, P., Ward-Thompson, D., & Barsony, M. 1993, ApJ, 406, 122 [Google Scholar]
  3. André, P., Ward-Thompson, D., & Barsony, M. 2000, in Protostars and Planets IV (Tucson: University of Arizona Press), 59 [Google Scholar]
  4. Aouad, C. J., James, P. A., & Chilingarian, I. V. 2020, MNRAS, 496, 5211 [CrossRef] [Google Scholar]
  5. Barrera-Ballesteros, J. K., Heckman, T., Sánchez, S. F., et al. 2021, ApJ, 909, 131 [NASA ADS] [CrossRef] [Google Scholar]
  6. Broekhoven-Fiene, H., Matthews, B. C., Harvey, P. M. et al. 2014, ApJ, 786, 37 [NASA ADS] [CrossRef] [Google Scholar]
  7. Burkhart, B. 2018, ApJ, 863, 118 [CrossRef] [Google Scholar]
  8. Calvet, N., Hartmann, L., Kenyon, S. J., & Whitney, B. A. 1994, ApJ, 434, 330 [NASA ADS] [CrossRef] [Google Scholar]
  9. Calzetti, D. 2013, in Secular Evolution of Galaxies, eds. Falcón-Barroso J., Knapen J. H., 419 [CrossRef] [Google Scholar]
  10. Dariush, A., Dib, S., Hony, S. et al. 2016, MNRAS, 456, 2221 [NASA ADS] [CrossRef] [Google Scholar]
  11. Dib, S. 2011, ApJ, 737, L20 [NASA ADS] [CrossRef] [Google Scholar]
  12. Dib, S. 2014, MNRAS, 444, 1957 [NASA ADS] [CrossRef] [Google Scholar]
  13. Dib, S. 2022, A&A, 666, A113 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Dib, S. 2023a, MNRAS, 524, 1625 [NASA ADS] [CrossRef] [Google Scholar]
  15. Dib, S. 2023b, ApJ, 959, 88 [NASA ADS] [CrossRef] [Google Scholar]
  16. Dib, S., & Basu, S. 2018, A&A, 614, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Dib, S., Kim, J., Vázquez-Semadeni, E., et al. 2007, ApJ, 661, 262 [NASA ADS] [CrossRef] [Google Scholar]
  18. Dib, S., Brandenburg, A., Kim, J., et al. 2008, ApJ, 678, 105 [Google Scholar]
  19. Dib, S., Walcher, C. J., Heyer, M., et al. 2009, MNRAS, 398, 1201 [NASA ADS] [CrossRef] [Google Scholar]
  20. Dib, S., Hennebelle, P., Pineda, J. E., et al. 2010, ApJ, 723, 425 [NASA ADS] [CrossRef] [Google Scholar]
  21. Dib, S., Helou, G., Moore, T. J. T., et al. 2012, ApJ, 758, 125 [NASA ADS] [CrossRef] [Google Scholar]
  22. Dib, S., Gutkin, J., Brandner, W., & Basu, S. 2013, MNRAS, 436, 3727 [NASA ADS] [CrossRef] [Google Scholar]
  23. Dib, S., Hony, S., & Blanc, G. 2017a, MNRAS, 469, 1521 [NASA ADS] [CrossRef] [Google Scholar]
  24. Dib, S., Schmeja, S., & Hony, S. 2017b, MNRAS, 464, 1738 [NASA ADS] [CrossRef] [Google Scholar]
  25. Dib, S., Bontemps, S., Schneider, N., et al. 2020, A&A, 642, A177 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Dunham, M. M., Allen, L. E., Evans, N. J., II, et al. 2015, ApJS, 220, 11 [Google Scholar]
  27. Eden, G., & Teyssier, R. 2024, MNRAS, 527, 6779 [NASA ADS] [Google Scholar]
  28. Egorov, O. V., Kreckel, K., Glover, S. C. O., et al. 2023, A&A, 678, A153 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Enoch, M. L., Evans, N. J., II, Sargent, A. I., & Glenn, J. 2009, ApJ, 692, 973 [Google Scholar]
  30. Evans, N. J. II, Dunham, M. M., Jørgensen, J. K. et al. 2009, ApJS, 181, 321 [NASA ADS] [CrossRef] [Google Scholar]
  31. Evans, N. J., II., Heiderman, A., & Vutisalchavakul, N. 2014, ApJ, 782, 114 [CrossRef] [Google Scholar]
  32. Federrath, C., & Klessen, R. S. 2012, ApJ, 761, 156 [Google Scholar]
  33. Federman, S., Megeath, S. T., Tobin, J. J., et al. 2023, ApJ, 944, 49 [NASA ADS] [CrossRef] [Google Scholar]
  34. Fischer, W. J., Padgett, D. L., Stapelfeldt, K. L., & Sewilo, M. 2016, ApJ, 827, 96 [NASA ADS] [CrossRef] [Google Scholar]
  35. Fischer, W. J., Megeath, S. T., Furlan, E., et al. 2017, ApJ, 840, 69 [NASA ADS] [CrossRef] [Google Scholar]
  36. Forbrich, J., Lada, C. J., Muench, A. A., et al. 2009, ApJ, 704, 292 [NASA ADS] [CrossRef] [Google Scholar]
  37. Froebrich, D., Schmeja, S., Smith, M. D., & Klessen, R. S. 2006, MNRAS, 368, 435 [NASA ADS] [CrossRef] [Google Scholar]
  38. Furlan, E., Fischer, W. J., Ali, B., et al. 2016, ApJS, 224, 5 [Google Scholar]
  39. Garduño, L. E., Zaragoza-Cardiel, J., Lara-López, M. A., et al. 2023, MNRAS, 526, 2479 [CrossRef] [Google Scholar]
  40. Glover, S. C. O., & Clark, P. C. 2012, MNRAS, 426, 377 [Google Scholar]
  41. Greene, T. P., Wilking, B. A., André, P., et al. 1994, ApJ, 434, 614 [NASA ADS] [CrossRef] [Google Scholar]
  42. Grisdale, K., Hogan, L., Rigopoulou, D., et al. 2022, MNRAS, 5133906 [Google Scholar]
  43. Guszejov, D., Markey, C., Offner, S. S. R., et al. 2022, MNRAS, 515, 167 [NASA ADS] [CrossRef] [Google Scholar]
  44. Gutermuth, R. A., & Heyer, M. 2015, AJ, 149, 64 [Google Scholar]
  45. Gutermuth, R. A., Megeath, S. T., Myers, P. C., et al. 2009, ApJS, 184, 18 [Google Scholar]
  46. Gutermuth, R. A., Pipher, J. L., Megeath, S. T., et al. 2011, ApJ, 739, 84 [Google Scholar]
  47. Haisch, K. E. Jr., Lada, E. A., & Lada, C. J. 2001, ApJ, 553, 153 [Google Scholar]
  48. Hasan, P., Mudasir, R., Saifuddin, Md., & Hasan, S. N. 2023, JAA, 44, 41 [NASA ADS] [Google Scholar]
  49. Hatchell, J., Fuller, G. A., Richer, J. S., et al. 2007, A&A, 468, 1009 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Hayashi, C. 1961, PASJ, 13, 450 [NASA ADS] [Google Scholar]
  51. Heiderman, A., Evans, N. J., II, Allen, L. E., et al. 2010, ApJ, 723, 1019 [NASA ADS] [CrossRef] [Google Scholar]
  52. Helou, G., Lu, N. Y., Werner, M. W., et al. 2000, ApJ, 532, 21 [Google Scholar]
  53. Hennebelle, P., & Chabrier, G. 2011, ApJ, 743, 29 [Google Scholar]
  54. Hernández, J., Hartmann, L., Megeath, T., et al. 2007, ApJ, 662, 1067 [Google Scholar]
  55. Hernández, J., Hartmann, L., Calvet, N., et al. 2008, ApJ, 686, 1195 [Google Scholar]
  56. Heyer, M., Gutermuth, R., Urquhart, J. S., et al. 2016, A&A, 588, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Hony, S., Gouliermis, D. A., Galliano, F., et al. 2015, MNRAS, 448, 1847 [NASA ADS] [CrossRef] [Google Scholar]
  58. Hsieh, T.-H., & Lai, S.-P. 2013, ApJS, 205, 5 [NASA ADS] [CrossRef] [Google Scholar]
  59. Jog, C. J. 2014, AJ, 147, 132 [NASA ADS] [CrossRef] [Google Scholar]
  60. Jose, J., Pandey, A. K., Samal, M. R. et al. 2013, MNRAS, 432, 3445 [NASA ADS] [CrossRef] [Google Scholar]
  61. Jose, J., Kim, J. S., Herczeg, G. J. et al. 2016, ApJ, 822, 49 [NASA ADS] [CrossRef] [Google Scholar]
  62. Kenyon, S. J., & Hartmann, L. 1995, ApJS, 101, 117 [Google Scholar]
  63. Kounkel, M. 2020, ApJ, 902, 122 [CrossRef] [Google Scholar]
  64. Kounkel, M., Zari, E., Covey, K., et al. 2023, ApJS, 266, 10 [NASA ADS] [CrossRef] [Google Scholar]
  65. Kristensen, L. E., & Dunham, M. M. 2018, A&A, 618, A158 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Kruijssen, J. M. D., & Longmore, S. N. 2014, MNRAS, 439, 3239 [NASA ADS] [CrossRef] [Google Scholar]
  67. Krumholz, M. R., & McKee, C. F. 2005, ApJ, 630, 250 [Google Scholar]
  68. Kuhn, M. A., de Souza, R. S, Krone-Martins, A., et al. 2021, ApJS, 254, 33 [NASA ADS] [CrossRef] [Google Scholar]
  69. Lada, C. J., & Wilking, B. A. 1984, ApJ, 287, 610 [NASA ADS] [CrossRef] [Google Scholar]
  70. Lada, C. J., Lombardi, M., & Alves, J. F. 2010, ApJ, 724, 687 [Google Scholar]
  71. Lenkić, L., Nally, C., Jones, O. C., et al. 2024, ApJ, 967, 110 [CrossRef] [Google Scholar]
  72. Li, G.-X. 2024, MNRAS, 528, 52 [Google Scholar]
  73. Lim, B., et al. 2015, in IAU General Assembly, Meeting 29, 2246322 [Google Scholar]
  74. Mairs, S., Johnstone, D., Kirk, H., et al. 2017, ApJ, 849, 107 [NASA ADS] [CrossRef] [Google Scholar]
  75. Mamajek, E. 2009, AIP Conf. Ser. 1158, Exoplanets and Disks: Their Formation and Diversity, eds. T. Usuda, M. Tamura, & M. Ishii (Melville, NY: AIP), 3 [CrossRef] [Google Scholar]
  76. Maraston, C., Pforr, J., Renzini, A., et al. 2010, MNRAS, 407, 830 [NASA ADS] [CrossRef] [Google Scholar]
  77. Marchuk, A. A. 2018, MNRAS, 476, 3591 [NASA ADS] [CrossRef] [Google Scholar]
  78. Megeath, S. T., Gutermuth, R., Muzerolle, J., et al. 2012, AJ, 144, 192 [NASA ADS] [CrossRef] [Google Scholar]
  79. Megeath, S. T., Gutermuth, R., Muzerolle, J., et al. 2016, AJ, 151, 5 [Google Scholar]
  80. Megeath, S. T., Gutermuth, R. A., & Kounkel, M. A. 2022, PASP, 134, 042001 [CrossRef] [Google Scholar]
  81. Meidt, S. E., Leroy, A. K., Rosolowsky, E., et al. 2018, ApJ, 854, 100 [NASA ADS] [CrossRef] [Google Scholar]
  82. Pacifici, C., Kassin, S. A., Weiner, B. J., et al. 2016, ApJ, 832, 79 [Google Scholar]
  83. Padoan, P., & Nordlund, Å. 2011, ApJ, 730, 40 [Google Scholar]
  84. Parravano, A., McKee, C. F., & Hollenbach, D. J. 2011, ApJ, 726, 27 [Google Scholar]
  85. Peltonen, J., Rosolowsky, E., Williams, T. G., et al. 2024, MNRAS, 527, 10668 [Google Scholar]
  86. Pessa, I., Schinnerer, E., Leroy, A. K. et al. 2022, A&A, 663, A61 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  87. Pokhrel, R., Guthermuth, R. A., Betti, S. K., et al. 2020, ApJ, 896, 60 [NASA ADS] [CrossRef] [Google Scholar]
  88. Pokhrel, R., Guthermuth, R. A., Krumholz, M. R., et al. 2021, ApJ, 912, L19 [CrossRef] [Google Scholar]
  89. Pokhrel, R., Megeath, S. T., Gutermuth, R. A., et al. 2023, ApJS, 266, 32 [NASA ADS] [CrossRef] [Google Scholar]
  90. Preibisch, T., McCaughrean, M. J., Grosso, N., et al. 2005, ApJS, 160, 582 [NASA ADS] [CrossRef] [Google Scholar]
  91. Rani, R., Moore, T. J. R., Eden, D. J., & Rigby, A. J., 2022, MNRAS, 515, 271 [CrossRef] [Google Scholar]
  92. Rebull, L. M., Cody, A. M., Covey, K. R., et al. 2014, AJ, 148, 92 [NASA ADS] [CrossRef] [Google Scholar]
  93. Robitaille, T. P., & Whitney, B. A. 2010, ApJ, 710, 11 [Google Scholar]
  94. Salpeter, E. E. 1955, ApJ, 121, 161 [Google Scholar]
  95. Schneider, N., Ossenkopf-Okada, V., Clarke, S., et al. 2022, A&A, 666, A165 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  96. Seigar, M. S. 2005, MNRAS, 36, 20 [Google Scholar]
  97. Sharma, S., Pandey, A. K., Ojha, D. K., et al. 2017, MNRAS, 467, 2943 [NASA ADS] [CrossRef] [Google Scholar]
  98. Shi, Y., Yan, L., Armus, L. et al. 2018, ApJ, 853, 149 [NASA ADS] [CrossRef] [Google Scholar]
  99. Strafella, F., Lorenzetti, D., Giannini, T., et al. 2015, ApJ, 798, 104 [NASA ADS] [CrossRef] [Google Scholar]
  100. Stutz, A. M., & Kainulainen, J. 2015, A&A, 577, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  101. Sun, J., Gutermuth, R. A., Wang, H., et al. 2022, MNRAS, 516, 5244 [Google Scholar]
  102. Tobin, J. J., & Sheehan, P. D. 2024, ARA&A, 62, 203 [NASA ADS] [CrossRef] [Google Scholar]
  103. Vázquez-Semadeni, E., González-Samaniego, A., & Colín, P. 2017, MNRAS, 467, 1313 [NASA ADS] [Google Scholar]
  104. Weisz, D. R., Johnson, L. C., Foreman-Mackey, D., et al. 2015, ApJ, 806, 198 [NASA ADS] [CrossRef] [Google Scholar]
  105. Wilking, B. A., Lada, C. J., & Young, E. T. 1989, ApJ, 340, 823 [NASA ADS] [CrossRef] [Google Scholar]
  106. Wilking, B. A., Meyer, M. R., Robinson, J. G., & Greene, T. P. 2005, AJ, 130, 1733 [NASA ADS] [CrossRef] [Google Scholar]
  107. Winston, E., Megeath, S. T., Wolk, S. J., et al. 2010, AJ, 140, 266 [NASA ADS] [CrossRef] [Google Scholar]
  108. Ybarra, J. E., Lada, E. A., Román-Zúñiga, C. G., et al. 2013, ApJ, 769, 140 [NASA ADS] [CrossRef] [Google Scholar]
  109. Zamora-Avilés, M., Vázquez-Semadeni, E., & Colín, P. 2012, ApJ, 751, 77 [CrossRef] [Google Scholar]
  110. Zhou, J. W., Dib, S., Wyrowski, F., et al. 2024a, A&A, 682, A173 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  111. Zhou, J. W., Dib, S., Juvela, M., et al. 2024b, A&A, 686, A146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  112. Zhou, J. W., Kroupa, P., & Dib, S. 2024c, PASP, 136, 094301 [NASA ADS] [CrossRef] [Google Scholar]
  113. Zhou, J. W., Kroupa, P., & Dib, S. 2024d, PASP, 136, 094302 [NASA ADS] [CrossRef] [Google Scholar]
1

We use the acronym YSO to designate the young stellar population in its ensemble. When referring specifically to young YSOs (Class 0, Class I, and the F Class), we employ the terms protostars. The older YSOs (Class II and Class III) are referred to as pre-main sequence (PMS) stars.

2

The exact nature of the flat-spectrum YSOs is still under investigation. A recent analysis using observations from the ATACAMA Compact Array (ACA) suggests that these objects are more likely disk dominated than envelope dominated and, as such, are more likely to be late-stage protostars or edge-on Class II YSOs (Federman et al. 2023).

3

The disk lifetime of a class II YSO is estimated as being the timescale over which the fraction of YSOs with a disk decays as a function of the mean isochronal ages of star forming regions of a certain range of ages (Haisch et al. 2001; Hernandez et al. 2007, 2008; Mamajek 2009).

4

Megeath et al. (2012, 2016) performed a comprehensive survey of the YSO population in the Orion clouds. However the YSOs in these works were categorized in the broad category of protostars and PMS stars and for this reason, the Orion data points could not be included in Figs. 11 and 12.

All Figures

thumbnail Fig. 1

Four realizations of the Milky Way-like protostellar mass function. The masses of N* = 500 protostars are randomly drawn from a Milky Way-like stellar mass function in the stellar mass range of 0.02–50 M⊙.
In the text
thumbnail Fig. 2

Top: four realizations of the birth time of protostars with a flat, constant SFH. Bottom: four realizations of the birth time of protostars with a Gaussian SFH. The peak and standard deviation of the Gaussian are set to tb,p = 2 Myr and σt = 1 Myr. For improved visibility, the values of tb are binned with a bin size of 2 × 105 yr.
In the text
thumbnail Fig. 3

One realization of the distribution functions of the lifetimes associated with each protostellar class. The lifetimes are drawn using Eq. (4) with the mean value and standard deviation associated with each lifetime (see text for details). The bin size for each protostellar lifetime distribution is adjusted to the mean value of that distribution, namely the bin sizes are 5×103 yr, 2.2 × 104 yr, 1.75 × 104 yr, 105 yr, and 2 × 106 yr for the Class 0, Class I, the F Class, Class II, and Class III, respectively.
In the text
thumbnail Fig. 4

Left: four realizations of the time evolution of the number of YSOs in different protostellar classes. Each cluster has N* = 500, a Milky Way-like IMF and a constant SFH. Right: similar to left panel, but here each cluster has N* = 500, a Milky Way-like IMF, and a Gaussian-like SFH. The Gaussian’s peak is at 2 Myr and the standard deviation is 1 Myr. The bin-size on the x-axis is 105 yr and the same value is used in subsequent plots.
In the text
thumbnail Fig. 5

Left: four realizations of the time evolution of the SFR measured using the populations of YSOs found in different classes and their respective lifetimes. The dashed line displays the time evolution of the true SFR. Each cluster has N* = 500, a Milky Way-like IMF, and a constant SFH. Right: similar to the left panel, but here each cluster has N* = 500, a Milky Way-like IMF, and a Gaussian-like SFH. The Gaussian’s peak is at 2 Myr and the standard deviation is 1 Myr.
In the text
thumbnail Fig. 6

Left: four realizations of the time evolution of the true SFR (SFRreal, dashed line) and the SFR measured using the total number of YSOs and an average lifetime for all Classes. Each clusters has N* = 500, a Milky Way-like IMF and a constant SFH. Right: similar to the left panel, but here, each cluster has N* = 500, a Milky Way-like IMF, and a Gaussian-like SFH. The Gaussian’s peak is at 2 Myr and the standard deviation is 1 Myr.
In the text
thumbnail Fig. 7

Top: time evolution of the ratio of the SFR measured using the populations of YSOs found in different classes (full lines) and the total population (dot-dashed line) and the real SFR that is generated in the models. Each line is an average over 250 clusters. Each cluster has N* = 500, a Milky Way-like IMF, and a flat SFH. Bottom: similar to the top panel, but here each clusters has N* = 500, a Milky Way-like IMF, and a Gaussian-like SFH. The Gaussian’s peak is at 2 Myr and the standard deviation is 1 Myr. In both panels, the dashed lines mark the 1σ uncertainty levels.
In the text
thumbnail Fig. 8

Top: time evolution of the ratio of the SFR measured using the populations of YSOs found in different classes (full lines) and the total population (dashed line) and the real SFR that is generated in the models. Each line is an average over 250 clusters. Each cluster has N* = 500, a Milky Way-like IMF, and a Gaussian-like SFH with the peak of the Gaussian located at tb,p = 2 Myr and a standard deviation of σt = 0.5 Myr. Bottom: similar to the top panel, but with tb,p = 2 Myr and σt = 2 Myr. For improved clarity, the 1σ uncertainties are not shown.
In the text
thumbnail Fig. 9

Four realizations of the variable protostellar mass function. The masses of N* = 500 protostars are randomly drawn from a TPL function in the stellar mass range of 0.02 to 50 M⊙. For each realization, the parameters that characterize the shape of the TPL mass function are randomly drawn from Gaussian distribution functions (see text for more details).
In the text
thumbnail Fig. 10

Same as the top panel of Fig. 7 but in this case the parameters that describe the shape of the IMF for each of the 250 clusters are randomly drawn from Gaussian distribution functions (see text for more detail). For more clarity, the 1σ uncertainties are not shown.
In the text
thumbnail Fig. 11

Top panel: offset between the SFR derived from the Class I population and the true SFR as a function of the ratio between the numbers of Class I and Class III PMS stars. The figure contains data from 250 clusters, each containing N* = 500, and for each cluster, it includes measurements at 100 epochs. The SFH is flat and the IMF resembles that of the Milky Way field. Bottom panel: similar to the top panel but showing the offset between the SFR derived from the Class II population and the true SFR as a function of the ratio of the numbers of Class II and Class III PMS stars.
In the text
thumbnail Fig. 12

Same as Fig. 11 but in this case the SFH is a time-dependent Gaussian-like function with the peak of the Gaussian located at 2 Myr and a standard deviation of σt = 1 Myr.
In the text
thumbnail Fig. A.1

top: Four realizations of the time evolution of the SFR measured using the populations of YSOs found in different classes and their respective lifetimes. The dashed line displays the time evolution of the true SFR. Each cluster has N* = 250, a Milky Way-like IMF and a constant SFH. bottom: Similar to the top panel, but here, each cluster has N* = 1000.
In the text
thumbnail Fig. A.2

Top row: same as Fig. 7 (top panel) but for cases where the width of the Gaussian distribution of the YSO lifetimes are taken as 0.25 the mean value (left) and 0.75 the mean value (right). Bottom row: Similar to Fig. 11 but for cases where the width of the Gaussian distribution of the YSO lifetimes are taken as 0.25 the mean value (left) and 0.75 the mean value (right).
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.