Effects of physical conditions on the stellar initial mass function: The low-metallicity star-forming region Sh 2-209

Recent work suggested that the variation of the initial mass function (IMF) of stars depends on the physical conditions, notably, the metallicity and gas density. We investigated the properties of two clusters, namely the main cluster (MC) and the subcluster (SC), in the low-metallicity HII region Sh 2-209 (S209) based on recently derived IMFs. We tested three previously published correlations using previous observations: the top-heaviness of the IMF in dependence on metallicity, the half-mass radius, and the most massive star in dependence on the stellar mass of the embedded clusters. For this region, two different galactocentric distances, namely 10.5 kpc and 18 kpc, were considered, where an age-distance-degeneracy was found for the previously determined IMF to be consistent with other formulated metallicity and density dependent IMFs. The determined half-mass radius (0.080 +/- 0.005) pc and the embedded cluster density (0.2 +/- 0.1) 10^6 Msun pc^-3 for the MC with an age of 0.5 Myr in S209 assuming a galactocentric distance of 18 kpc support the assumption that a low-metallicity environment results in a denser cluster, which leads to a top-heavy IMF. Thus, all three tests are consistent with the previously published correlations. The results for S209 are placed in the context with the IMF determination within the metal-poor cluster in the star-forming region NGC 346 in the Small Magellanic Cloud.


Introduction
The initial mass function (IMF) gives information about the initial mass distribution of stars in a stellar system (e.g.Bastian et al. 2010;Kroupa et al. 2013;Hopkins 2018).Together with previously formulated relations, it describes the mass distribution of a population of stars at birth of the embedded cluster.These previously formulated relations are the m max −M ecl relation (Yan et al. 2023), with the maximum stellar mass m max and the stellar mass of the embedded cluster M ecl , and the variation of the IMF with metallicity and density of a molecular cloud core in which the population forms as an embedded star cluster (Marks et al. 2012;Jeřábková et al. 2018;Yan et al. 2021).Another relation is the r h −M ecl relation by Marks & Kroupa (2012), in which r h is the half-mass radius of the embedded cluster at birth, just prior to gas expulsion.The concepts and relations are explained in more detail in Sect. 2.
With this contribution, we test whether the previously determined constraints on the variation of the IMF with metallicity and density of a star-forming gas cloud are consistent with the recently published observations of two young star clusters by Yasui et al. (2023).These clusters lie in the low-metallicity starforming region Sh 2-209, which is located in the outer region of the Milky Way.Furthermore, we compare the observed clusters to the relation between r h and M ecl , obtained by Marks & Kroupa (2012).
The IMFs of solar metallicity regions have been wellstudied.One example is the Taurus region (e.g.Luhman 2000; Briceño et al. 2002;Thies & Kroupa 2008) with a solar metallicity [Fe/H] = −0.01 ± 0.05 (D'Orazi et al. 2011) or NGC 2264 by Sung & Bessell (2010) with [Fe/H] ≈ −0.13 (Netopil et al. 2016).Another solar metallicity region is the Trapezium cluster in the Orion nebular cluster (ONC) with [Fe/H] = −0.01 ± 0.04 (D'Orazi et al. 2009), for which Muench et al. (2002) derived the mass function from B stars to the deuterium-burning limit.On the other hand, the IMF of an entire low-metallicity region is usually only known over a limited stellar-mass range.For example, Harayama et al. (2008) determined the IMF of the NGC 3603 star cluster within the NGC 3603 HII region in the Milky Way (Goss & Radhakrishnan 1969) within a stellar mass range of 0.4−20 M .Kuncarayakti et al. (2016) found it to have a metallicity between Z = 0.004 and 0.008.Because the two-body relaxation time for high-mass stars is similar to the cluster age, Harayama et al. (2008) suggested that a composition of primordial and dynamical effects leads to the observed top-heaviness of the IMF.Another example is NGC 346.This star-forming region in the Small Magellanic Cloud (SMC) has a metallicity of [Fe/H] ≈ −0.72 (Rochau et al. 2007): the IMF was determined within a stellar mass range of 0.8−60 M by Sabbi et al. (2008) over a region spanning ≈40 pc.The authors reported that they did not find any environmental effects on the IMF.This holds for the global IMF derived by Sabbi et al. (2008), who noted its flattening at the centre and attributed it to mass segregation, which must be primordial because NGC 346 is younger than its mass-segregation timescale (see Table 1 in Sabbi et al. 2008).This means that more massive stars formed in the centre, as underlined by Fig. 5 in Sabbi et al. (2008), which supports the idea of variations of the IMF.We discuss NGC 346 in Sect.6.2.
S209 is the first low-metallicity ([O/H] = −0.5 dex) starforming region whose IMF was determined over the wide mass range of 0.1−20 M (Yasui et al. 2023).Our approach is to assume the previously constrained variation of the IMF (see Sect. 2.3) and apply it to the two clusters, the main cluster (MC) and the subcluster (SC), given the available information on their metallicities.With this, we calculate the radii of the two clusters and check them for consistency with the variable IMF and the independently derived relation of r h and M ecl by Marks & Kroupa (2012, see Sect. 2.4).This places constraints on the distances of both clusters.We also check for consistency in the relation of m max versus the M ecl (see Sect. 2.2) published by Yan et al. (2023).Section 2 documents the previously constrained birth relations, Sect. 3 introduces the observational data, Sect. 4 explains the methods, and Sect. 5 reports the results.The discussion and conclusion are provided in Sects.6 and 7, respectively.

Previous results on the birth relations
In this section, we give an overview of the results from previous research on the properties of newly formed embedded clusters that we applied to the data of S209.

The canonical initial mass function
In a stellar system, the initial mass distribution of the stars is of special interest to resolve the formation process and the stellar system's properties.The IMF gives the number of stars with a certain mass.While Salpeter introduced a one-power-law IMF (Salpeter 1955), Kroupa (2001) found that the mass distribution among stars in the solar neighbourhood follows a two-part power-law formulation, the canonical IMF.It is defined as the number of stars in the mass interval m to m + dm, dN = ξ(m) dm, where ξ(m) = k • k i m −α i , k and k i are constants, and the exponents are α 1 = 1.3 for 0.08 ≤ m/M < 0.5 and α 2 = 2.3 for 0.5 ≤ m/M < m max , m is the stellar mass, and m max is the most massive star in the embedded cluster (Sect.2.2).

The m max −M ecl relation
Based on the IMF being a probability density distribution function (PDF), it might be thought that the IMF of one massive star cluster is similar to that of multiple small star clusters.This contradicts the fact that in a massive embedded star cluster, more massive stars can form, which is not possible in low-mass clusters because the mass of a very massive star would exceed the mass of the small cluster.This concept is documented by the m max −M ecl relation (Weidner & Kroupa 2006).A comparison of the results of the two sampling methods, namely the random and optimal sampling introduced by Kroupa et al. (2013), shows that the m max −M ecl relation has a physical origin and is not a statistical effect (Yan et al. 2023).With this relation, an estimate of M ecl is obtained by only having m max .Yan et al. (2023) quantified the relation between m max and M ecl (see Figs. 5 and 6 in Yan et al. 2023).

The varying initial mass function
While many Milky Way and Large Magellanic Cloud clusters follow the canoncial IMF (Kroupa 2002), subsequent data on the mass-to-light ratios (Dabringhausen et al. 2009) and the occurrence of low-mass X-ray bright binaries (Dabringhausen et al. 2012) and of the stellar content of low-concentration globular clusters led to the suggestion that the IMF becomes top-heavy at low metallicity and high gas density (Marks et al. 2012).Star counts in nearby low-metallicity star-forming regions support this result (Schneider et al. 2018;Kalari et al. 2018 at the 2σ confidence level).Assuming this variation of the IMF, Marks et al. (2012), Jeřábková et al. (2018), andYan et al. (2021) introduced a metallicity-dependent relation between the initial gas density ρ gas and the exponent α 3 .
Thus, the exponent α 3 is related to the initial average cloud core density ρ gas and the initial gas metallicity of the cluster [Z/H] (Jeřábková et al. 2018;Yan et al. 2021), where ρ gas is the initial average density of the cloud core in units of 10 6 M pc −3 , and ρ ecl is the stellar mass density of the embedded cluster.These two are connected with the relation ρ ecl = • ρ gas , where the canonical value of = 0.3 is adopted for the star formation efficiency (Lada & Lada 2003;André et al. 2014;Megeath et al. 2016;Banerjee & Kroupa 2018;Wirth et al. 2022).Furthermore, is assumed to be constant throughout the cluster, such that the gas exactly tracks the stars.We note that ρ ecl constitutes the mathematically idealised moment of the highest effective density of the embedded cluster at the idealised time when all stars are born instantly, such that the binary population subsequently evolves to the observed populations (Marks & Kroupa 2012).The real physical embedded cluster has a more complex formation history in which the stars form over one to a few crossing times of the contracting molecular cloud core before it expands due to the expulsion of residual gas (e.g.Kroupa et al. 2001;Della Croce et al. 2023).

The r h −M ecl relation
Based on a study of binary populations, Marks & Kroupa (2012) introduced the canonical relation of r h and M ecl This is the half-mass radius of the embedded cluster in the most compact configuration, that is, when it is deeply embedded, and prior to any gas expulsion.It is required in order for the emanating open star cluster to have the observed population of binary stars, given that star formation provides a binary fraction of almost 100% (Kroupa 1995a ,b; Marks & Kroupa 2011).A108, page 2 of 9

Sh 2-209
S209 is an HII region located at a Galactic longitude 151 • .6062and Galactic latitude −0 • .24with right ascension RA = 4 11 06.7 and declination Dec = +51 09 44 (Wenger et al. 2000).With data from Gaia Collaboration (2021), Yasui et al. (2023) determined an astrometric distance of ≈2.5 kpc.Assuming the solar galactocentric distance to be R = 8.0 kpc, S209 has a galactocentric distance of about 10.5 kpc (Yasui et al. 2023).On the other hand, Foster & Brunt (2015) determined a heliocentric distance of ≈10.58 kpc from radial velocity measurements of the Canadian Galactic Plane Survey in the radio regime.This result is consistent with the value determined by Chini & Wink (1984), using photometric and spectroscopic data from exciting stars.The galactocentric distance for this case, again assuming R = 8.0 kpc, is therefore approximately 18 kpc.
S209 is a star-forming environment consisting of a MC for which a sample of 1500 objects is identified and a SC with 350 members (Yasui et al. 2023).The observed mass range is 0.1−20 M .Depending on the age and distance, the observed mass range differs.For instance, for younger ages, masses down to 0.02 M are considered.When the first break mass is lower than the minimum mass of 0.1 M , the first break mass is considered as the minimum mass.Following Yasui et al. (2023), we adopt the term "break mass" as the mass at which the slope of the IMF changes.The observed mass range was determined from the colour-magnitude diagram.However, some stars that are below this limit but have a large K excess were considered in the derivation of the K-band luminosity functions that were used to obtain the IMFs (Yasui et al. 2023).To identify the objects in the clusters, near-infrared image data were used, which also show that the cluster is surrounded by residual gas.This means that gas expulsion has already occurred.Using the metallicity map from Eilers et al. (2022), we estimated the metallicity of the environment to be [Fe/H] ≈ −0.25 dex for the smaller distance and [Fe/H] ≈ −0.5 dex for the outer position (see Fig. 1).Based on previous studies, the result by Yasui et al. (2023) has an oxygen abundance of [O/H] = −0.5 dex.For a more detailed discussion of the metallicities, we refer to Sect. 6.

Methods
In order to test for the relations we introduced in Sect.2, M ecl needs to be calculated.Ensuring that the IMF is smooth, we determined the constant k using the number of objects N in the star cluster and the equation N = m max m min ξ(m) dm.After constraining the IMF, the mass of the embedded cluster can be computed, with Here, m 1 (=0.5 M for the canonical IMF) is the first break mass, and m 2 (=1 M in Eq. ( 1)) is the second break mass, while m min = 0.08 M is approximately the hydrogen-burning limit.
In order to determine r h , ρ ecl needs to be computed.For this, Eqs. ( 2) and (3) were applied with the known metallicity.Using the relation ρ ecl = • ρ gas with = 0.3 as explained in Sect.2.3, ρ ecl was obtained.After determining M ecl and ρ ecl , we calculated the pre-gasexpulsion r h , assuming spherical symmetry, The results of these calculations are listed in Sect. 5.The uncertainties were calculated by Gaussian error propa- , where ∆ f is the uncertainty of the quantity f , x i is a variable, and ∆x i is its corresponding error).Because the constant k was calculated via integration and was included in the determination of M ecl , where a second integration was applied, a large error can be obtained for M ecl .Furthermore, the errors for the exponents Γ i contribute strongly.This can cause large uncertainties for the M ecl values, which lead to large uncertainties for some other results.Additional uncertainties such as extinction due to the interstellar medium or the intracluster medium were not accounted for here because a correction for them was already applied in the analysis by Yasui et al. (2023).
After computing M ecl and r h with Eqs. ( 5) and ( 6) based on the ρ gas obtained with Eqs. ( 2) and (3), we compared the results to the r h −M ecl relation described in Sect.2.4.This is displayed in Figs. 3 and 4 in Sect.6.

Results
In Yasui et al. (2023), the parameters given for the IMF varied for the different galactocentric distances assumed (10.5 kpc or 18 kpc), for the ages from 0.5 to 10 Myr, and for the considered cluster (MC or SC).In order to give an overview, different cases were considered in this paper.In the main part, the parameters of the best-fit IMF are evaluated, and in the Appendix A the results for ages of 1, 5, and 10 Myr are listed.
Using the best-fit results that Yasui et al. (2023) obtained, displayed in Table 1 for the MC and in Table 2 for the SC, we plot the IMFs in Fig. 2. The uncertainties of the parameters were adopted from Table 9 and Sect   Table 1.Best-fit values for the MC (Yasui et al. 2023).
Parameter Distance of 10.5 kpc Distance of 18 kpc Age/Myr 5 0.5 log 10 m 1 /M −1.2 −0.9 log 10 m 2 /M −0.9 +0.1 We proceeded as follows in order to test whether the observational constraints on α 3 by Yasui et al. (2023) are consistent with the formulation of the variation of the IMF given by Eqs. ( 1) and (2) (Sect.2.3).By adopting the reported values of α 3 , Eq. ( 3) was used to calculate ρ gas and Eq. ( 5) to calculate M ecl .The  The resulting densities, radii, and masses are displayed in Table 3 for the MC and in Table 4 for the SC.The provided results are formal solutions with their 1σ uncertainties.The 1σ upper limit would be physically relevant.
Ideally, all stellar masses would be summed from star counts as a constraint on the stellar mass of the cluster, but these data were not available to us.

Sh 2-209
Considering the m max −M ecl relation that we explained in Sect.2.2 and the results for the MC, m max = 20 M (Yasui et al. 2023), we obtain M ecl ≈ 800 M , which is compatible with the derived M ecl for the best-fit parameters of MC at a galactocentric distance of 18 kpc (see Table 3).For older ages, it is more consistent with the MC being at a galactocentric distance of 10.5 kpc (see Table A .11).Based on radio emission data, Richards et al. (2012) found the ionised gas mass to be greater than 10 3 M .This is more consistent with the values obtained for a galactocentric distance of 10.5 kpc (the authors even suggest a distance of 4.9 kpc).(Vilchez & Esteban 1996), which is consistent with the value derived here.Rudolph et al. (2006) used the data of Vilchez & Esteban (1996) and derived [O/H] ≈ −0.29 dex, which is 30% lower than the value for the inner position and about 50% higher than the value for the larger distance.Using the work of Caplan et al. (2000), Rudolph et al. (2006) derived an oxygen abundance of ≈−0.58 dex, which is inconsistent with the metallicity determined using the map by Eilers et al. (2022) for the inner position and is about 35% lower than the value for the outer position.Finally, Yasui et al. (2023) suggested [O/H] ≈ −0.5 dex, which is consistent with Eilers et al. (2022) for the outer position of S209.This maintains the preference that S209 has a galactocentric distance of 18 kpc, as derived by Foster & Brunt (2015) and Chini & Wink (1984) using spectroscopic and photometric data.
In the determination above, no binary stars are assumed.In the binary-star theorem by Kroupa (2008), the majority of star formation results in binary systems.Therefore, a binary fraction of 95% is adopted in the following, which corresponds to the updated pre-main-sequence eigenevolution model by Belloni et al. (2017).In order to give a maximum correction, we assumed the binaries to be composed of equal-mass stars, but in reality, the mass of the secondaries is somewhat lower than that of the primaries (Moe & Di Stefano 2013;Sana et al. 2012).The value of α 3 is not affected by a realistic binary, triple, and quadrupole population (Weidner et al. 2009;Kroupa & Jerabkova 2018).Hence, the mass M ecl is nearly doubled.Because the density is determined by the metallicity and α 3 (see Eqs. ( 1) and ( 3)), which means that it is independent of M ecl , the half-mass radius in Eq. ( 6) increases accordingly.
The canonical relation between r h and M ecl (Eq.( 4)) and the corresponding results are plotted in Fig. 3 for the MC and in Fig. 4 for the SC.The plots for the other ages are displayed in Appendix A. In Fig. 3, the value for a galactocentric distance of 18 kpc, which respects the binary-star theorem, lies in the error band and is therefore consistent within 1σ confidence, while the data points of the galactocentric distance of 10.5 kpc deviate more from the canonical Eq. ( 4), but still lie in the 1σ regime.The deviation of the result for a galactocentric distance of 10.5 kpc decreases with the binary assumption.In the plots in which the other ages are younger than 10 Myr, which are displayed in Appendix A, the data points for the MC at a galactocentric distance of 18 kpc fit the relation better, while the points for a galactocentric distance of 10.5 kpc do not deviate by more than 2σ.This changes for an older age (see Fig. A.8), where a galactocentric distance of 10.5 kpc corresponds better to the relation by Marks & Kroupa (2012).
All the half-mass radii are smaller than expected from the canonical Eq. ( 4), except for the age of 10 Myr (see Fig. A.8).This underlines the possibility that a low-metallicity starforming region might be denser because the mass is included in a smaller volume than is determined for other very young star clusters of higher metallicity.This would be reminiscent of the smaller radii of low-metallicity stars compared to metal-rich stars of the same mass, suggesting that self-regulation may play an important role for the constitution of an embedded cluster (see also Yan et al. 2023).
In Fig. 4, the values for a galactocentric distance of 18 kpc of the SC again agree with the relation by Marks & Kroupa (2012) within the 1σ confidence region.The results for the 10.5 kpc distance are only consistent with Eq. ( 4) within the 3σ confidence.Hence, the prediction by Marks & Kroupa (2012), that is, Eq. ( 4), agrees with the determined values for the MC and SC for both distances, while the results for a galactocentric distance of 18 kpc represent the observational data better.However, given the large uncertainties, no highly significant conclusions can be drawn.

NGC 346
NGC 346 is the largest active star-forming region in the Small Magellanic Cloud (Rickard et al. 2022), for which Sabbi et al. (2008) argued that the IMF is the same as in the Galaxy based on their survey, which spanned ≈40 pc around the central embedded cluster.Figure 5 in Sabbi et al. (2008) showed, however, that the centremost 4 pc have a top-heavy IMF with α 3 = (2.03 ± 0.14), which is consistent with a star cluster density of ρ ecl ≈ (0.25 ± 0.20) × 10 6 M pc −3 , according to Eq. ( 3), assuming a metallicity [Fe/H] ≈ −0.72 (Rochau et al. 2007), which we adopted for [Z/H] assuming an error of 0.3 to account for [α/H].With this, the total cluster mass is M tot ≈ 3.9 × 10 5 M , and with Eq. ( 6), we can give an upper limit on the half-mass radius of the central region of r h < (0.7 ± 0.2) pc, which is consistent with the foregoing calculation when the uncertainties caused by gas expulsion that lead to expansion are included.At larger distances between 4 and 9 pc, Sabbi et al. (2008) reported a canonical IMF (α 3 ≈ 2.36 ± 0.09), which is likely a mixture of on-site lower-mass embedded clusters and ejected massive stars from the starburst cluster (Oh & Kroupa 2016;Oh et al. 2015).At larger distances, α 3 ≈ (2.43 ± 0.09) and for distances larger than 14 pc, α 3 ≈ (3.08 ± 0.14), which reflects a population of low-mass embedded clusters lacking massive stars (Weidner & Kroupa 2006;Yan et al. 2023).Although the ejected massive stars induce a steeper slope of the MF in the inner part, mass segregation is more present in this case, which leads to the flatter slope in the central region.The centre of this whole star-forming system on a scale of ≈40 pc across is reminiscent of the Orion-South star-forming cloud where only the ONC formed massive stars, while much of the southern part of the cloud formed many low-mass embedded clusters lacking massive stars (Hsu et al. 2012).In contrast to the assertion by Sabbi et al. (2008), NGC 346 thus appears to be consistent with the variation of the IMF given by Eqs. ( 1) and ( 2), but more work is needed to quantify details.For example, it will be necessary to quantify whether many low-mass embedded clusters formed with IMFs lacking massive stars compared to the centred starburst cluster at larger distances from the centre of NGC 346.

Conclusions
We analysed two clusters of the star-forming region S209 with respect to their initial conditions.In particular, the low metallicity was taken into account, which allowed us to determine ρ ecl (Eq.( 3)), by using the determination of the IMF power-law index by Yasui et al. (2023).We constrained r h and compared it A108, page 5 of 9 to the independently developed canonical relation of r h and M ecl (Eq.( 4)).
Thus, the IMF data from Yasui et al. (2023) for the lowmetallicity region S209 result in a top-heavy IMF.The determined value of M ecl for ages of about 0.5 Myr for the MC at a galactocentric distance of 18 kpc is most consistent with the value M ecl ≈ 800 M deduced with the relation of m max and M ecl derived by Yan et al. (2023).This changes for older ages, where the resulting mass, M ecl , for the smaller galactocentric distance of 10.5 kpc is more consistent with the m max −M ecl relation.
The comparison of the determined metallicity with the map by Eilers et al. (2022) for the outer position and the independently derived metallicities by Yasui et al. (2023) and Vilchez & Esteban (1996) indicate a preferred galactocentric distance of 18 kpc.The value for the inner position (10.5 kpc) is not consistent with the metallicity of S209.
The calculated value for the r h of the MC at a galactocentric distance of 18 kpc is consistent with the canonical r h −M ecl relation for young ages, namely being within the 1σ confidence region.The binary theorem introduced by Kroupa (2008) tightens these results.The values of r h for the MC and SC at a galactocentric distance of 10.5 kpc we found are only consistent within 2σ or 3σ confidence for young ages.On the other hand, they replicate the m max −M ecl relation (Yan et al. 2021) better for an age of 10 Myr.This confirms the previously mentioned agedistance degeneracy.
In addition, we showed that the variation of the IMF (Eq.( 2)) is consistent with NGC 346, which shows a top-heavy IMF in the innermost region and a lack of massive stars with increasing distance to the centre.This is reminiscent of the ONC.
In conclusion, S209 is consistent with a low-metallicity region resulting in a denser star-forming region with a topheavy IMF, as quantified by Marks et al. (2012) and Yan et al. (2021).The findings for the MC with a galactocentric distance of 18 kpc are consistent with the relation of m max and M ecl by Yan et al. (2023) for young ages, while a galactocentric distance of 10.5 kpc is more consistent for an older age.When the binary-star theorem by Kroupa (2008) is included, the values determined for M ecl and r h in this paper are consistent with the r h −M ecl relation derived by Marks & Kroupa (2012).To quantify this result further, more low-metallicity environments have to be examined in the future, and the heliocentric distance measurement of S209 needs to be improved.

Fig. 1 .
Fig. 1.Metallicity map from Eilers et al. (2022), where the Galactic centre at (0, 0) is indicated by the cross, the Sun is located at (−8.0 kpc, 0) and is indicated by the circled dot, and the two possible positions of S209 are indicated by an asterisk.The arrows show the rotation of the stars and are coloured according to their mean self-calibrated metallicity.

Fig. 2 .
Fig. 2. Reconstructed IMF with values from Yasui et al. (2023).The dashed lines represent the canonical IMF (see Sect. 2.1).The upper line shows the MC, and the lower line shows the SC.

Fig. 3 .
Fig. 3. Canonical relation of M ecl and r h (Eq.(4)) shown by the blue line.The data points are the calculated values for the MC in S209.The caps indicate the 1σ and 2σ uncertainty boundaries.
we estimated the iron abundance [Fe/H] at the two possible positions of S209.Eilers et al. (2022) also provided a map for [O/Fe] (see their Fig.6), from which we estimated [O/Fe] ≈ 0.03 dex for the inner and [O/Fe] ≈ 0.07 dex for the outer position.With the relation [O/H] = [O/Fe] + [Fe/H], we obtained an oxygen abundance [O/H] of ≈−0.22 dex for the smaller and ≈−0.43 dex for the larger distance.Because the maps by Eilers et al. (2022) do not account for small-scale variations in the metallicity, the individually observed metallicities may differ from those implied by the Eilers et al. (2022) map.Yasui et al. (2023) discussed different values for [O/H], namely ≈−0.43 dex
Fig. A.5.As Fig. 4, but for the MC with an age of 5 Myr.

Table 3 .
Results for the MC.

Table 4 .
Results for the SC.

Table A .
(Yasui et al. 2023)for the SC with an age of 5 Myr(Yasui et al. 2023).
Table A.7. Results for the MC with an age of 5 Myr.Fig.A.6.As Fig. 4, but for the SC with an age of 5 Myr.Table A.8. Results for the SC with an age of 5 Myr.Table A.9. Best-fit values for the MC with an age of 10 Myr (Yasui et al. 2023).

Table A .
(Yasui et al. 2023s for the SC with an age of 10 Myr(Yasui et al. 2023).