Open Access
Issue
A&A
Volume 672, April 2023
Article Number A29
Number of page(s) 23
Section Galactic structure, stellar clusters and populations
DOI https://doi.org/10.1051/0004-6361/202245457
Published online 27 March 2023

© The Authors 2023

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

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

1. Introduction

Characterization of binary stellar systems in star clusters is a crucial step in shedding light on various fields of stellar astrophysics, including the dynamic evolution of stellar systems, star formation, and stellar evolution. For example, a robust determination of the physical properties of a star cluster, including the mass function (MF) and the total mass, would require significant knowledge of its populations of binary stars. Furthermore, the stellar evolution of a binary system can strongly differ from that of single stars, and depends on the physical properties of the system, including the mass ratio, binding energy, and orbital period.

Various approaches can be used to identify and characterize binary stars in stellar clusters. For example, binaries can be detected from radial-velocity variation or from photometric variability. While these methods have the advantage of constraining each binary system, they are either limited to bright stars or biased toward binaries with short periods.

In this work, we follow an alternative approach based on the fact that binary stars formed by couples of main sequence (MS) stars exhibit redder colors with respect to single MS stars. Hence, binary stars can be identified as stars lying on the red side of the MS fiducial line in the color–magnitude diagram (CMD; e.g., Romani & Weinberg 1991; Bolte 1992; Rubenstein & Bailyn 1997; Bellazzini et al. 2002; Clark et al. 2004; Richer et al. 2004; Zhao & Bailyn 2005; Milone et al. 2009). The main advantages of this approach are that (i) it requires observations in only two different filters, hence requiring a small amount of telescope time; (ii) it allows us to simultaneously investigate large numbers of stars in the CMD; and (iii) the detection efficiency does not depend on binary properties such as period and inclination.

Clearly, high-precision photometry is needed to disentangle binaries from single stars. Moreover, the correction for differential reddening and the identification of field stars that contaminate the cluster CMD are crucial ingredients to infer the fraction of binaries from photometry. On the other hand, the use of this method comes with some caveats, as follows.

This approach has been used to investigate the binaries of a large sample of 67 Galactic GCs (see e.g., Sollima et al. 2007; Milone et al. 2012, 2016; Ji & Bregman 2015) using homogeneous photometry from images collected with the Hubble Space Telescope, and for Galactic open clusters with the combination of multiple surveys (Malofeeva et al. 2022, 2023). Photometric binaries have been widely investigated in in Galactic open clusters (OCs), although most studies published so far are limited to a smaller number of clusters (e.g., Sharma et al. 2008; Sollima et al. 2010; Cordoni et al. 2018; Ebrahimi et al. 2022, and references therein). Historically, the main challenge in constraining the properties of photometric binaries in open clusters was provided by the contamination from field stars, which is critical in most Galactic open clusters.

In the past few years, the Gaia mission has provided high-precision photometry, proper motions, and parallaxes of nearly 2 billion stars in the Milky Way (Gaia Collaboration 2018, 2021, Gaia DR2 and DR3). This exquisite dataset has allowed accurate separations to be obtained between the bulk of the OC members and the background and foreground stars (e.g., Cantat-Gaudin et al. 2018, 2020; Bossini et al. 2019; Dias et al. 2022).

In this work, we use Gaia DR3 data to investigate a sample of 78 Galactic OCs, constrain their structure parameters, and characterize their photometric binaries. The paper is structured as follows: Sect. 2 describes the clusters sample choice, the cluster member selection, and differential reddening corrections. In Sect. 3, we derive the physical parameters of clusters, while Sect. 4 is dedicated to investigation of the binaries. Finally, we present and discuss the overall results in Sect. 5.

2. Data and star sample selection

To identify unresolved binary MS-MS stars in the CMD, we need high-precision photometry, low field contamination, negligible reddening effect, and a sufficiently large sample of stars. Hence, we downloaded the Gaia DR3 astro-photometric catalog of all the 269 open clusters analyzed in Bossini et al. (2019), and selected only those clusters that satisfied our requirements. Specifically, we included only clusters with a large number of stars, that is, more than 150 cluster members, low total reddening, namely E(B − V)≤0.5, and with reasonable differential reddening variation and field-star contamination. Specifically, as the determination of accurate binary star fractions relies on precise estimates of cluster age, distance, reddening, and metallicity, we opted to use the clusters discussed in Bossini et al. (2019) in place of larger cluster catalogs (see e.g., Cantat-Gaudin et al. 2018, 2020). Indeed, Bossini and collaborators restricted their analysis to those clusters for which the isochrone-fitting procedure would produce the most reliable age and distance estimates, thus selecting 269 out of the 1229 clusters analyzed in Cantat-Gaudin et al. (2018). Including clusters with less accurate estimates could result in additional sources of uncertainty. The same clusters are also included in the more recent analysis of Dias et al. (2021). Specifically, Dias et al. (2021) provide clusters’ proper motion estimates as well as isochrones fit parameters, such as age, reddening, distance modulus and metallicity, determined through direct isochrone fitting. Hence, in the following we will make use of the cluster’s properties inferred by Dias and collaborators. Moreover, the total number of clusters analyzed in this work is 78, and their main properties are summarized in Table C.1. The distributions of age and the total number of stars of the selected clusters are shown in Fig. 1. Clearly, the selected clusters span a wide age range, from a few million years (e.g., ∼10 Myr Pozzo 1) to almost 7.5 Gyr for NGC 6791. Similarly, the number of cluster members varies from a minimum of 130 for NGC 2632 to more than 5000 stars for NGC 6705. Therefore, the large variety of parameters allows us to explore the behavior of binary stars in relation to different cluster structural parameters; for example, mass1, age, dimension, and other physical parameters.

thumbnail Fig. 1.

Properties of analyzed clusters. Top: histogram distribution of the ages (in logarithm scale) of the 78 open clusters analyzed here. Bottom: histogram distribution of the number of cluster members. See the text for details.

The identification of unresolved binary stars relies on the separation between single MS stars and MS-MS binary stars in the CMDs and requires high-precision photometry. To select a sample of well-measured cluster stars, we exploited the exquisite photometry and astrometry provided by Gaia DR3. We first used some Gaia DR3 quality parameters that are indicative of the quality of the photometric and astrometric measurements, including the renormalized unit weight error (RUWE) and the photometric excess noise. These quantities are used to select a sample of stars with high-quality astrometry and photometry, as in Cordoni et al. (2018).

To select cluster members, we adopted the approach described in Cordoni et al. (2018) and illustrated in Fig. 2 for NGC 2099. The main steps of this method can be summarized as follows.

thumbnail Fig. 2.

Cluster members selection. Procedure carried out to select the probable NGC 2099 cluster members located within ∼30 arcmin from the cluster center. Panel a: stellar proper motions relative to the average cluster motion against the G magnitude. The blue line separates stars with cluster-like proper motions from the remaining stars. Panel b: stellar G magnitude as a function of parallax. Stars with cluster-like parallaxes are located between the two blue lines. The vector point diagram of proper motions is plotted in panel c, whereas panel d shows the spatial stellar distribution. Finally, panel e illustrates the G vs. GBP − GRP observed CMD, together with the best-fit isochrone. Likely cluster members are represented with azure crosses, while brown crosses indicate stars with cluster-like proper motions and parallaxes in the reference field.

We first analyzed the proper-motion diagram and computed the proper motion of each star relative to the cluster: μ R = ( μ α cos δ μ α cos δ ) 2 + ( μ δ μ δ ) 2 $ \mu_{\mathrm{R}} = \sqrt{(\mu_\alpha \cos\delta - \left\langle \mu_\alpha \cos\delta\right\rangle)^2 + (\mu_\delta - \left\langle \mu_\delta\right\rangle)^2} $. Here, ⟨μαcos δ⟩ and ⟨μδ⟩ are the average cluster proper motions from Dias et al. (2021). As an example, Fig. 2a shows μR as a function of the G magnitude for stars in the field of view of NGC 2099. We divided the interval between G = 9.0 and 20.0 mag into bins of 0.5 mag and calculated for each bin the median value of μR and the corresponding root-square mean (σ). To derive the blue line of Fig. 2a, we first shifted the median value of μR of each bin by 3σ and associated this quantity with the average magnitude of the star in the bin. The blue line is derived by linearly interpolating these points. We considered the stars on the left side of the blue line to be candidate cluster members (azure crosses in Fig. 2a). For completeness, we show the proper-motion diagram in Fig. 2c.

We then improved the selection of cluster members using stellar parallaxes (see Fig. 2b for NGC 2099). Similarly to what we did in panel a, we divided the diagram into magnitude bins and derived the median parallax of the proper-motion-selected stars and the corresponding rms (σ) of each bin. We shifted the median parallax value of each bin by a ±3σ region and associated the resulting values with the average magnitude of the stars in each of those bins. We derived the left and right blue lines plotted in panel d of Fig. 2 by linearly interpolating these points. We considered the proper-motion-selected stars between the two lines as probable cluster members and represented them with azure crosses.

Clearly, the sample of selected candidate cluster members could include residual field stars with proper motions and parallaxes similar to those of the cluster. To estimate the contribution of contaminating field stars, we identified a reference field located well outside the selected cluster field, that is, at a radial distance ten times larger than the radius containing 50% of the total number of stars determined in Dias et al. (2021). We selected the field stars that populate the same regions of the G vs. μR and G vs. ω diagrams populated by cluster members and that pass the selection criteria described above. These field stars are marked with brown crosses in Fig. 2.

Finally, we defined the outer cluster radius as the radius where the density of cluster stars is equal to the density of residual field stars. Figure 2d shows the spatial distribution of the selected NGC 2099 stars. The CMD of the 2031 cluster members and the 112 selected reference field stars is plotted in Fig. 2e.

2.1. Differential reddening

As our goal is to identify unresolved binary stars in the CMD, we need to minimize the broadening of stellar colors and magnitudes due to spatial reddening variations. To correct the photometry for differential reddening, we followed the method described in Milone et al. (2012). Briefly, we started by selecting a sample of reference stars, which is composed of bright MS cluster members, and calculated the corresponding fiducial line in the CMD. We used the absorption coefficients by Casagrande & VandenBerg (2018) to determine the direction of the reddening and calculate the color and magnitude distance of each reference star from the fiducial along the reddening direction.

To calculate the differential reddening associated with each star, we selected the 25 neighboring reference stars and calculated the median distance along the reddening line. We excluded each reference star from the determination of its own differential reddening. The corresponding error is calculated as the root mean scatter of the distance values divided by the square root of N − 1. We refer to Milone et al. (2012) for a more detailed description of the procedure.

The results of the differential reddening correction are illustrated in Fig. 3 for the open cluster NGC 6819. Specifically, panels a and b show the original and corrected photometry, respectively, while panel c shows the differential reddening map where blueish colors indicate negative DR and greenish colors indicate positive variations. Results on NGC 6819 highlight the effectiveness of the differential-reddening-correction procedure. Indeed, all the photometric sequences in the corrected CMDs are much narrower when compared to the original CMD.

thumbnail Fig. 3.

Differential reddening correction. Comparison of the GRP vs. GBP − GRP CMD before (panel a) and after the correction for differential reddening (panel b). The differential reddening map is shown in panel c, color-coded according to the top color bar.

We applied the above procedure to all the clusters of our sample. In 44 out of 78 clusters, the inferred reddening variations are significantly larger than the corresponding errors. In the following analysis, for these clusters marked with the flag DR in Table C.1 we make use of the photometry corrected for differential reddening. Their differential-reddening maps are available as supplementary material, together with the individual differential-reddening-correction catalogs.

2.2. The color–magnitude diagrams

The differential-reddening-corrected CMDs of the selected cluster members from the 78 clusters are shown in Figs. A.1A.4. Clusters are sorted by age, as determined in Dias et al. (2021), from the youngest, namely Pozzo 1 (10 Myr), to the oldest, namely NGC 6791 (7.5 Gyr). Visual inspection of these CMDs shows that at least 30−35 clusters exhibit the extended MS turn-off (eMSTO). Such clusters are flagged with eMSTO in Table C.1. They comprise about ∼45% of the 72 studied clusters younger than 2 Gyr, which is the limiting age for the occurrence of the eMSTO in star clusters (e.g., Milone et al. 2009; Goudfrooij et al. 2014; Niederhofer et al. 2015; Cordoni et al. 2018).

Our results corroborate the evidence that the eMSTO is a typical feature of young clusters. However, we note that at least 35 clusters with ages between 0.01 and 2 Gyr are consistent with simple populations. Hence, the eMSTO is not a ubiquitous property of these clusters. We refer to Cordoni et al. (in prep.) for further investigation of the eMSTO phenomenon.

3. Physical parameters of open clusters

Accurate estimates of the structural parameters of the main clusters are required to characterize their binary stars. In the following, we use Gaia DR3 astrometry and photometry to determine the core radius, density, and mass of each cluster.

3.1. Density profile

To constrain the density profile of each cluster, we divided the field of view into overlapping annular bins of constant width, defined as the difference between the external and internal radius. The width has been chosen so as to maximize the spatial information while ensuring good statistics. We calculated the stellar density in each annulus, (∑*), as the number of cluster stars brighter than GRP = 18.5 mag per square arcmin. We assumed that the Gaia DR3 catalog is complete for magnitudes brighter than GRP = 18.5 as discussed in Boubert & Everall (2020) and Everall & Boubert (2022). To account for residual field contamination, we determined the density of field stars, selected as described in Sect. 2, and subtracted such value from the observed cluster field density. The uncertainties that we associated with stellar density are determined as Poisson errors.

We plotted ∑* as a function of the average radius of each annulus (r) and performed a least-square fit of the observed density profile with the function discussed in Elson et al. (1987, EFF profile):

μ ( r ) = μ 0 ( 1 + r 2 a 2 ) γ / 2 , $$ \begin{aligned} \mu (r) = \mu _0 \left(1+\frac{r^2}{a^2}\right)^{-\gamma /2}, \end{aligned} $$(1)

where μ0 represents the central density in stars arcmin−2, a a scale radius, and γ is the index of the power law at large radial distances from the cluster center. The cluster core radius, rc, is linked to a, γ by the relation:

r c = a ( 2 2 / γ 1 ) 1 / 2 . $$ \begin{aligned} r_{\rm c} = a \left(2^{2/\gamma }-1 \right)^{1/2}. \end{aligned} $$(2)

The EFF profile differs from the classical King profile (King 1962) for the absence of tidal truncation. As an example, Fig. 4 compares the best-fit EFF and King profiles (black solid and dashed lines, respectively) for the test case cluster NGC 2099. We verified that the conclusions are not significantly affected by the adopted bin size. To this end, we repeated the analysis for bins of different sizes and used annuli with different radii that include the same number of cluster stars.

thumbnail Fig. 4.

Density profile of NGC 2099. Azure points with error bars represent the observed stellar density profile. The black solid line indicates the best-fitting EEF (Elson et al. 1987) profile, while the corresponding uncertainty is indicated by the gray shaded area. The inferred core radius is indicated with the black solid vertical line, while the shaded region is its uncertainty. The best-fitting parameters are quoted in the bottom-left inset. The dashed line shows the best-fitting King (King 1962) profile.

Moreover, as we limit our analysis to stars in a specific magnitude range, we lose the information about all stars that lie outside the analyzed range. Hence, the recovered central stellar density may be underestimated. To deal with this bias, we exploited the recovered MFs discussed in Sect. 3.2 to compute the total expected number of stars in each cluster. Specifically, we multiplied the fitted central density by the ratio between the expected total number of stars inferred from the MF, and the number of stars used to compute the density profile.

3.2. Cluster mass and mass function

To estimate the total stellar mass of each cluster, we first determined the MF of MS stars brighter than GRP = 18.0/18.5 according to each specific CMD quality. In particular, we adopted GRP = 18.5 in 46 clusters and GRP = 18.0 in the remaining 32 clusters. We first divided the MS into ten magnitude intervals, each containing 10% of the total number of MS stars. Stellar magnitudes are then converted into masses using the relations by Marigo et al. (2017).

As a by-product, we determined the half-mass radius (rh) interpolating the cumulative mass distribution.

As an example, the different mass intervals used to infer the MF of NGC 2099 are highlighted with different colors in the CMD of Fig. 5a, where we also indicate the mean mass of the MS stars in each bin (M). In Fig. 5b, we plot the logarithm of N normalized to the mass interval of each bin ΔM as a function of the logarithm of M.

thumbnail Fig. 5.

Mass function determination. Left panel: GRP vs. GBP − GRP CMD of NGC 2099. The magnitude intervals used to derive the MF are highlighted with different colors, while the median mass of each bin is indicated on the right. Right panel: MF for the stars plotted in the left panel, i.e., log(NM) vs. stellar mass. The blue and azure straight lines represent the best-fit straight lines of low- (M ≤ 1 M) and high-mass stars (M > 1 M), respectively. The red line is obtained by fitting the whole MF. The slopes are quoted in the bottom-left legend, while the total stellar mass of the cluster is indicated in the top-right inset.

The present-day stellar MF is described by a power law with index α:

d N / d M = k · m α , $$ \begin{aligned} \mathrm{d}N/\mathrm{d}M = k\cdot m^{-\alpha }, \end{aligned} $$(3)

with k a normalization constant. Taking the logarithm of Eq. (3), we obtain the equation of a straight line with slope α.

log ( d N / d M ) = log ( k ) α log ( m ) . $$ \begin{aligned} \log (\mathrm{d}N/\mathrm{d}M) = \log (k)-\alpha \log (m). \end{aligned} $$(4)

Hence, the observed MF can be interpolated with a straight line in log(ΔNM) vs. log(M) space by means of least squares, as in the right panel of Fig. 5 for NGC 2099. Visual inspection of the MFs of several clusters, including NGC 2099, reveals a change in the MF slope around M = 1 M.

To account for this feature of the MF, we divided each cluster stellar catalog into a low-mass region with M ≤ 1 M, and a high-mass region with M > 1 M, and fitted the function described by Eq. (4) to each of the two regions individually. Results of the three interpolations, that is, over the whole MS, just the low-mass region, and just the high-mass region, are represented with red, blue, and azure colors, respectively. We find that the change of slope is significant beyond the 3σ level in 30, or 38%, of the 78 analyzed clusters.

The curves that provide the best fit to all stars, low-mass stars, and high-mass stars of NGC 2099 are represented with red, blue, and azure lines, respectively, in Fig. 5, where we also quote the corresponding values of α, indicated on left. Specifically, we find an overall slope of α = 1.70 ± 0.17, while the best-fit slopes of low- and high mass-stars MF are 1.18 ± 0.04 and 2.08 ± 0.12, respectively.

Finally, we then used the derived overall MF slopes to infer cluster masses, which are listed in Table C.2. Speficically, we integrated the derived MF beyond the range of observed stellar masses, including evolved stars. Implicitly, we made the reasonable assumption that the MFs do not vary beyond the stellar mass interval used to derive the MF itself. However, as evolutionary timescales after the turn-off are very short, the mass range of stars beyond the turn-off would be negligible even in the case of a varying MF. We derived the half-mass relaxation time (trh) as in Spitzer (1987)

t rh = 0.138 M tot 1 / 2 r h 3 / 2 G 1 / 2 m ¯ ln ( γ M tot / m ¯ ) , $$ \begin{aligned} t_{\rm rh} = 0.138\frac{M_{\rm tot}^{1/2}r_{\rm h}^{3/2}}{G^{1/2}\bar{m}\ln (\gamma M_{\rm tot}/\bar{m})} , \end{aligned} $$(5)

where γ = 0.11 and m ¯ $ \bar{m} $ is the mean stellar mass. The inferred cluster parameters are listed in Table C.2.

The results for the MF fits are summarized in Table C.2 for the 78 analyzed clusters. As shown in Fig. 6a, we find that the MF slope distribution peaks around α ∼ 2.11 and ranges from ∼ − 0.5 to ∼3.5. Moreover, 61 out of 78 clusters (78%) are consistent with a Salpeter MF (α = 2.35) at the 1σ level. Remarkably, we find that the MF slope depends on the dynamical age of the cluster, which is defined as the ratio between the age and the half-mass relaxation time, as shown in Fig. 6b. Specifically, α is consistent with the Salpeter index in young and intermediate-age clusters, while it starts to decrease for clusters older than ∼500 Myr. The same trend is also found by Ebrahimi and collaborators (see e.g. figure of Ebrahimi et al. 2022).

thumbnail Fig. 6.

Mass functions properties. Panel a: distribution of the overall MF slope inferred for the 78 analyzed clusters. Panel b: MF slope α vs. cluster age. Panel c: cluster age vs. explored stellar mass range in MS stars.

To test the effect introduced by considering different stellar mass ranges in the determination of the cluster MFs, we adopted the following approach. Briefly, we simulated 1000 CMDs of clusters with the same age, distance, and metallicity as the observed clusters, and following a Salpeter MF. We the added errors in color and magnitude similar to the observed ones, and finally re-computed the MF of each simulation by following the same procedure and adopting the same mass interval adopted for the real clusters. Moreover, we randomly varied the observed stellar mass ranges by adding random Gaussian fluctuations. We find that the recovered MF slopes are consistent with the Salpeter input one to the level of 0.01. We can therefore safely assume that the bias introduced by analyzing different stellar mass intervals does not significantly affect our results. Moreover, in the case of a double-slope input MF, the recovered single slope is characterized by a larger spread with respect to the single-slope case. Specifically, in the case of a 1.35/2.35 input double slope MF, we find that the overall spread of the recovered single slope is of order 0.4.

To compare the MFs of clusters with different ages, we normalized each MF to the number of stars with −0.2 ≤ log(M/M)≤0.05, which is the mass interval common to all the analyzed clusters (see Fig. 6c). Figure 7 shows the stacked MF for dynamically young, intermediate, and old clusters (panels a–c, respectively), and for all clusters (panel d). Specifically, we define dynamically young clusters as those with log(age/trh) < 0.2, intermediate those with 0.2 ≤ log(age/trh) < 0.75, and old as those with log(age/trh)≥0.75. The adopted boundaries are defined so that approximately the same number of clusters fall into each group.

thumbnail Fig. 7.

A common mass function. Panels a–c: normalized stacked MFs for young (log(Age/trh)≤0.2), intermediate (0.2 < log(Age/trh) < 0.75), and old clusters (log(Age/trh)≥0.75). The normalization choice is described in Sect. 3.2. Blue and azure solid lines represent the straight line best-fit of the low-mass (M ≤ 1 M) and high-mass (M > 1 M) stars, respectively. The derived slope is indicated in the bottom-left legends. Panel d: same as panels a–c, but for the whole sample of analyzed clusters.

These figures corroborate the evidence for a change in the MF slope around 1 M. We derived the curves that provide the least-squares best fit of the stacked MFs by using the stars more massive and less massive than 1 M separately. The results are provided by the azure and blue straight lines shown in Fig. 7. The analysis shows that the MFs of the dynamically young, intermediate, and old clusters follow approximately the same pattern, with the slopes of low- and high-mass stars being consistent within 1σ. Finally, we repeated the procedure including all clusters, regardless of their age. The result is shown in Fig. 7d. Remarkably, we find that the general trend is, again, well approximated by a double power-law MF, with slopes consistent with those found in Sollima (2019) for ∼120 000 stars in the solar neighborhood; that is, 1.34 and 2.41/2.68 in the sub- and supersolar regimes, respectively.

4. The binary fraction

The unresolved binary systems composed of two MS stars (MS-MS binaries) exhibit redder colors and brighter magnitudes than both stellar components. Hence, they lie on the bright and red side of the MS fiducial line in most CMDs constructed with the Gaia filters. The exact location of MS-MS binaries in the CMD depends on the mass of the primary stars (Mprim) and on the ratio between the mass of the secondary and primary star, q = Msec/Mprim. The equal-mass binary systems (i.e., q = 1) form a sequence parallel to the MS but ∼0.75 mag brighter. On the other hand, binaries with q close to zero would approach the MS fiducial line. For intermediate mass-ratio values, we observe that the smaller the value of q, the closer the binary system will lie to the MS fiducial.

To estimate the fraction of unresolved binary systems, we follow the approach of Milone et al. (2012) and Cordoni et al. (2018), which is illustrated in Fig. 8 for NGC 2099. For each cluster, we limited the analysis to binaries with a mass ratio greater than a fixed threshold value, qlim, ranging from 0.6 to 0.7 according to the overall quality of the cluster CMD. Indeed, binaries with small mass-ratio values, for example q < qlim, lie too close to the MS and are almost indistinguishable from single stars.

thumbnail Fig. 8.

Binary stars selection. Panel a: illustration of the binary-selection procedure. The region enclosed by the azure solid line, dubbed region A, contains all MS stars and binary stars with a primary mass between 0.80 and 1.31 M. The gray-shaded region, i.e., region B, defines the locus of binary-star candidates with a mass ratio greater than 0.6, respectively. Finally, the dashed black line and red solid lines represent the fiducial line of MS stars and binary stars with mass ratio q = 0.6, as indicated in the top-right legend. Panel b: MS and binary stars identified in the cluster field of NGC 2099, marked with red and blue crosses, respectively. Panel c: same as panel b but for stars located in the reference field, well outside the cluster radius. Stars that fall into region A are marked with black crosses.

First, we defined region A, which is enclosed by the azure line in Fig. 8, and includes all single and binary cluster stars with G RP bright G RP G RP faint $ G_{\mathrm{RP}}^{\mathrm{bright}}\leq G_{\mathrm{RP}}\leq G_{\mathrm{RP}}^{\mathrm{faint}} $. Region B is instead defined as the portion of region A containing only candidate MS-MS binary stars with q > qlim, and is colored in gray in Figs. 8a–c.

To define these two regions, we first derived the fiducial line of MS stars. To do this, we divided the MS into magnitude bins and computed the median and the 68.27th percentiles of the color and magnitude distributions of the stars of each bin (hereafter σGBP − GRP and σGRP). The MS fiducial line (dotted line in Fig. 8a) is derived by linearly interpolating these median points.

Specifically, the boundaries of region A, which is enclosed by the azure line in Fig. 8, are defined as follows: (i) the blue border is obtained by shifting the fiducial line four times σGBP − GRP to the left; (ii) the bright and fainter boundaries are the loci of binary systems where the primary stars have magnitudes G RP bright $ G_{\mathrm{RP}}^{\mathrm{bright}} $ and G RP faint $ G_{\mathrm{RP}}^{\mathrm{faint}} $, respectively, and mass ratios ranging from 0 to 1; (iii) the reddest border is defined by shifting the fiducial line of equal mass binaries to the red by four times σGBP − GRP.

Region B, shaded in gray in Fig. 8, is the portion of region A that includes binaries with q > qlim. It is located on the red side of the fiducial line of binary stars with q = qlim. To avoid significant contamination from single stars, we have chosen the value of qlim of each cluster in such a way that the fiducial line of binaries with q = qlim is redder than the MS fiducial shifted by four times σGBP − GRP. We adopted qlim = 0.6 or qlim = 0.7 depending on the uncertainties of the photometry of each cluster.

Stars that fall into region B in NGC 2099 are colored in red in Fig. 8b, while the remaining stars that lie in region A are marked with blue crosses. To account for the residual contamination from field stars in the cluster CMD, we derived the numbers of field stars with cluster-like proper motions, and the parallaxes in the CMD of stars in the reference field normalized to the same area as the cluster field (see Fig. 8c for NGC 2099). The fraction of binaries is estimated as:

f bin q q lim = N cl B N fl B N cl A N fl A , $$ \begin{aligned} f_{\rm bin}^{q\ge q_{\rm lim}} = \frac{N_{\rm cl}^\mathrm{B}-N_{\rm fl}^\mathrm{B}}{N_{\rm cl}^\mathrm{A}-N_{\rm fl}^\mathrm{A}} , \end{aligned} $$(6)

where the subscripts “cl” and “fl” refer to the cluster and the reference field, respectively. We used a similar formula to calculate the fraction of binaries within the core.

The fractions of binaries range from less than 0.2 to ∼0.7. As shown in Fig. 9, there is no evidence for significant correlations between the fractions of binaries and or between cluster age and mass. However, Fig. 9c reveals a mild correlation between the total binary fraction and the central density derived in Sect. 4.4. Moreover, to account for the different distances of the analyzed clusters, we converted the central density from stars arcmin−2 to stars pc−2. In this section, we discuss only the relevant correlations between the total and core binary fractions and the physical parameters of the clusters; however, the plot in the appendix shows all the correlations (see e.g. Figs. C.1C.2). Remarkably, while there is no correlation between the total binary fraction and the age or the dynamical age of a cluster (ρ = 0.20, 0.27 respectively, bottom row of Fig. C.1), there are hints of a mild correlation between the core binary fraction and cluster age and dynamical age (ρ = 0.54, 0.49 respectively, bottom row of Fig. C.2). We investigate these relationships in the following sections.

thumbnail Fig. 9.

Relation between binary fractions and host clusters. Panels a–c: total binary fraction as a function of cluster dynamical age, i.e., age normalized over the half-mass relaxation time, total cluster mass, and central density in stars pc−2. Azure lines and shaded regions indicate the straight line that best fits the data and its 1σ dispersion. The values of Spearman’s rank coefficient (ρ) and the corresponding p-values are indicated in the top-left corner of each plot. Panels b1–b3: same as panel b, but with clusters grouped into three dynamical age groups, defined as young, intermediate, and old clusters. Panels c1–c3: same as panel c, for the three groups of clusters.

To further explore possible connections between binary fractions and the physical parameters of the clusters, we divided the clusters into the three “age” groups defined in Fig. 7 and repeat the analysis. Remarkably, the total binary fraction mildly correlates with cluster mass for the dynamically young cluster (Fig. 9b1), with log(age/trh)≤0.2, while it weakly anti-correlates in old clusters (Fig. 9b3). Concerning the central stellar density, Figs. 9c1–c3 reveal that the binary fraction is linked with the central density in young and intermediate age clusters, while these two quantities seem to be uncorrelated in older clusters.

4.1. Binaries and blue stragglers

To investigate the relations between binaries and blue stragglers (BSSs), we followed the procedure illustrated in Fig. 10 for NGC 188. We first defined the regions C and D in the CMD of stars in the cluster field. Region C, indicated by the azure shade, is populated by candidate BSSs. The left border is defined by the zero-age MS (ZAMS) isochrone (dashed line), while the right boundary is provided by the combination of the MS and Sub Giant Branc (SGB) fiducial line blue-shifted by four times the color error (grey line), and the fiducial line of equal-mass binaries, shown in red.

thumbnail Fig. 10.

Procedure used to estimate the fraction of BSSs in NGC 188. Cluster and reference-field stars are marked with gray circles and brown crosses, respectively. The gray dashed line represents the ZAMS. The gray solid line is the fiducial line of the MS and the SGB shifted to the blue by four times the observational color error (σ). The fiducial line of MS-MS and SGB-SGB binaries, shifted by four times σ to the red is represented by the red line. Region C of the CMD is indicated as an azure-shaded area, while region D, marked with the gray shade, is enclosed by gray and red solid lines and by the two horizontal black dashed lines. Cluster field stars in regions C and D of the CMD are marked with azure solid and black-shaded circles, respectively.

Region D, shaded in grey, is mostly populated by MS stars within 1 mag of the MS turn-off, that is, with GRP in the range ( G RP Toff , G RP Toff + 1 ) $ (G_{\mathrm{RP}}^{\mathrm{Toff}}, G_{\mathrm{RP}}^{\mathrm{Toff}}+1) $. The fraction of BSSs is estimated as

f BSS = N cl C N ref C N cl D N ref D , $$ \begin{aligned} f_{\rm BSS} = \frac{N_{\rm cl}^\mathrm{C}-N_{\rm ref}^\mathrm{C}}{N_{\rm cl}^\mathrm{D}-N_{\rm ref}^\mathrm{D}}, \end{aligned} $$(7)

where N cl C $ N_{\mathrm{cl}}^{\mathrm{C}} $ and N cl D $ N_{\mathrm{cl}}^{\mathrm{D}} $ are the numbers of stars in the regions C and D, respectively, in the cluster-field CMD, while N ref C $ N_{\mathrm{ref}}^{\mathrm{C}} $ and N ref D $ N_{\mathrm{ref}}^{\mathrm{D}} $ are the corresponding number of stars in the reference field.

We confirm the previous conclusion that clusters younger than ∼300 Myr do not host BSSs (Vaidya et al. 2022; Rain et al. 2021). As shown in Fig. 11, the fraction of BSSs shows no evidence of a correlation with cluster mass or metallicity, but there is a mild correlation with cluster age (see also Rain et al. 2021).

thumbnail Fig. 11.

Blue stragglers properties. Fraction of BSSs against cluster age (panel a), cluster mass (panel b), and cluster metallicity.

Figure 12 compares the fraction of BSSs and the binary fraction and shows no evidence of a significant correlation. Nevertheless, when we consider two groups of clusters with different central density values, expressed in stars pc−2 (see Sect. 3.1). We note that low-density clusters (μ0 < 1 stars pc−2) define a sequence that extends toward large values of fbin but small fractions of BSSs, whereas the remaining clusters (blue dots) have, on average, larger values of fBSS and define a steeper sequence in the fBSS vs. fbin plane. The high-density cluster Tombaugh 1 is a possible exception.

thumbnail Fig. 12.

Fraction of blue stragglers vs. binary fraction. Red colors indicate the clusters with central density μ0 < 10 stars pc−2, whereas the remaining clusters are colored blue.

4.2. Mass-ratio distribution

To investigate binaries with different mass ratios, we adapted the method of Milone et al. (2012) to open clusters. Briefly, we divided region B of the CMD into three regions, B1, B2, and B3 as shown in Fig. 13a for NGC 2099, where the candidate binaries are marked with white crosses. Each region encloses binaries with different mass ratios and the values of q corresponding to its left and right boundaries are derived with the criterion that each region covers similar areas in the CMD. Each subregion comprises binaries in a different mass-ratio interval Δqi. To compare the binary fractions inferred for the three mass-ratio intervals, we estimated the normalized binary fraction: νi = fbin, ii, where the index i ranges from 1 to 3, and the fbin, i is the fraction of binaries in the regions B1, B2, and B3 of the CMD.

thumbnail Fig. 13.

Binary stars mass-ratio as a function of mass-ratio in NGC 2099. Reproduction of the CMD of NGC 2099 plotted in Fig. 8b (panel a). The shaded regions (B1, B2, and B3) of the CMD comprise binaries with different mass ratios. The candidate binaries are marked with white crosses. Panel b shows the mass-ratio distribution of binaries in NGC 2099. The χ2 and the probability that the observed binary fractions come from a flat mass-ratio distribution are quoted in panel b.

Results are shown in Fig. 13b for NGC 2099, where we plot νi as a function of q. The horizontal gray line and shaded region indicate the total normalized binary fraction with q > qlim and its uncertainty, respectively.

We performed a χ2 test to investigate whether or not the binaries with q > qlim of each cluster are consistent with a flat mass-ratio distribution. We first computed the χ2 between the observed νi values and the flat mass-ratio distribution corresponding to the horizontal gray line of Fig. 13a, χ OBS 2 $ \chi^2_{\rm OBS} $. We then generated 1000 CMDs with the same number of stars and the same fraction of binaries with q > qlim as the observed one, but with a flat mass-ratio distribution. We calculate the normalized binary fractions, μi, corresponding to each simulation and derived the corresponding χ2 values with respect to the flat mass-ratio distribution, in close analogy with what we did for the observed CMD, χ SIM 2 $ \chi^2_{\rm SIM} $.

The probability that the observed distribution is consistent with a flat trend is estimated as the fraction of simulations with χ SIM 2 $ \chi^2_{\rm SIM} $ >  χ OBS 2 $ \chi^2_{\rm OBS} $ and corresponds to one in the case of a flat distribution. As an example, the observed mass-ratio distribution of binaries in NGC 2099 has a 9% probability of being flat.

Results for all clusters are shown in Fig. B.1, while the values of the χ2 and Pflat are provided in Table C.2. We find that in 66 out of 78 clusters, the mass-ratio distribution has a probability higher than 10% of being flat.

To further investigate the mass-ratio distribution of binary systems, we combine the results from all clusters as shown in Fig. 14. To do this, we divided the normalized binary fractions by the total fraction of binaries. The latter is estimated as f BIN TOT = f BIN q > q lim / ( 1 q lim ) $ f_{\mathrm{BIN}}^{\mathrm{TOT}} = f_{\mathrm{BIN}}^{q > q_{\mathrm{lim}}}/(1-q_{\mathrm{lim}}) $. We adopted a 2D-histogram representation to show the distribution of the normalized fraction of binary stars in each mass-ratio bin. The color code is indicated by the top-right color bar. When results from all clusters are combined together, we find a flat mass-ratio distribution for binary systems with q ≥ 0.6, as can be inferred by the position of the darkest cells in Fig. 14.

thumbnail Fig. 14.

Binary fraction mass-ratio distribution. The binary fraction in each mass-ratio bin has been normalized over the width of the bin, in q, and the total fraction of binaries extrapolated assuming a flat mass ratio distribution. Each cell is color coded according to the number of observations in that particular bin, as specified in the top right color bar.

4.3. Binary fraction as a function of the mass of the primary star

To determine the fraction of binary systems as a function of primary star mass, we divided region A of the CMD into three different subregions, as indicated in Fig. 15, each containing approximately the same number of stars. We computed the fraction of binary stars in each region by means of Eq. (6). The resulting distribution is shown in panel b2, where the horizontal gray line and gray shaded region indicate the binary fraction of NGC 2099. Visual inspection of the mass distribution reveals a weakly increasing or flat trend, as indicated by the Pflat value defined in the previous section. Specifically, the value of the χ2 and its probability prevent us from making further conclusions. Overall, we find that 72 out of 78 clusters are consistent with a flat mass distribution of binary stars, with Pflat > 10% (see Fig. 16).

thumbnail Fig. 15.

Binary stars as a function of primary stars Mprim in NGC 2099. The adopted regions are colored differently, as defined in the main text. The mass distribution of candidate binary stars with q ≥ _lim is shown in panel b.

thumbnail Fig. 16.

Binary fraction primary star distribution. The azure line and the shaded region represent the median trend and the 1σ confidence interval.

4.4. Radial distribution

To investigate the radial distribution of binary systems, we divided the field of view into three circular regions containing approximately the same number of stars. For each radial bin, we estimated the binary fraction with the procedure of Sect. 4.

The results are illustrated in Fig. B.3, where we plot the fraction of binaries as a function of the radial distance from the cluster center. We find that 15 clusters out of the 78 analyzed here exhibit a low probability of having a flat radial distribution, that is, Pflat < 5%.

To further investigate the radial distribution of binary stars, we considered the three groups of clusters with different ages – defined in Sect. 3.2 – separately. For each age group, we plotted the fraction of binaries normalized to the total binary fraction against the radial distance from the cluster center normalized to the core radius. As shown in Figs. 17a–c, we find that the binaries of both young and intermediate-age clusters (i.e., clusters younger than ∼800 Myr) exhibit a nearly flat radial distribution, with some hints of two peaks corresponding to the cluster center and to a radius of about four times the core radius. The binaries of the old clusters are more centrally concentrated than single stars. The binary fraction is maximum in the core and declines by a factor of two at the radial distance of four core radii. Such behavior is qualitatively consistent with the mild correlations between core binary fractions and cluster age. For completeness, we show in Fig. 17d the radial distribution of binaries for the entire sample of studied clusters.

thumbnail Fig. 17.

Binary stars radial distribution. Binary fraction with q ≥ qlim normalized to the total binary fraction, as a function radial distance from the cluster center in core radius units. The azure lines mark the median trends, while the shaded areas correspond to the 1σ dispersion. Panels a–c include clusters in the three different dynamical age intervals defined in the main text (see Sect. 3.2), while panel d refers to all clusters. Color coding and age bins are the same as in Fig. 6.

4.5. Comparison with Galactic globular clusters

In this subsection, we compare our findings on Galactic open clusters with those of Milone et al. (2012, 2016) for Galactic globular clusters (GCs). As the works by Milone and collaborators are limited to the central regions of GCs, we only compare the binaries in the core regions of globular and open clusters. Moreover, to minimize the effects of low statics on our conclusions, we restrict the investigation to the open clusters with more than 500 stars.

Results are shown in Fig. 18, where we plot the total fraction of binaries in the core of Galactic open clusters (aqua-blue circles) and in GCs (red crosses) as a function of the logarithm of cluster age (panel a), dynamical age (panel b), and cluster mass (panel c), as derived in Baumgardt & Hilker (2018).

thumbnail Fig. 18.

Comparison between OCs and GCs. Total fraction of binaries in the core of open clusters (aqua circles) and GCs (red crosses) as a function of the logarithm of cluster age (panel a), dynamical age (panel b), and cluster mass (panel c). The GC binary fractions are taken from Milone et al. (2012, 2016), while GC masses are from Baumgardt & Hilker (2018).

Intriguingly, Fig. 18a reveals that the observed core binary fractions with q > 0.6 of OCs is larger than that of GCs. However, panel b of the same figure may suggest that some GCs, namely NGC 6652, NGC 6838, and Ruprecht 106, characterized by a higher core binary fraction, are consistent with the trend defined by OCs. However, we observe that both GCs and OCs exhibit a qualitatively similar increasing pattern in the fbin vs. log(age/trh) plane. On the other hand, we observe opposite relations between core binary fractions and cluster mass. Indeed, OCs exhibit a correlation, while GCs clearly define an anti-correlation. Possible exceptions may be represented once again, by NGC 6652, NGC 6838, Rup 106, and NGC 6835, which seem to follow the relation defined by OCs. However, a more detailed analysis is required to draw firm conclusions.

While the fraction of binaries in the GC core anti-correlates with the mass of the host cluster, there is no evidence of a similar trend between the fraction of binaries in the core of open clusters and their masses.

5. Summary and conclusion

We used the Gaia DR3 data to investigate the population of unresolved binaries along the MS of 78 Galactic open clusters. To this end, we first carefully selected the probable cluster members based on their parallaxes and proper motions. Moreover, we corrected the photometry for effects of differential reddening when needed and accounted for residual field star contamination. We also exploited the isochrones from the Padova database (Marigo et al. 2017) that, according to Dias et al. (2021), provide the best match with the observed CMDs.

For each cluster, we measured the fraction of binaries with mass ratios larger than a fixed threshold, q > qlim, where qlim ranges from 0.6 to 0.7, depending on the photometric errors. We also extrapolated the total binary fraction by assuming a flat mass-ratio distribution. We derived the binary fractions in the entire cluster field of view and in the region within the cluster core. In addition, we investigated the fraction of binaries in three different mass-ratio bins and in different intervals of primary-star mass. Finally, we analyzed the radial distribution of binaries.

To properly investigate the relations between binary stars and the host star cluster, we derived various cluster properties:

  • We estimated the EFF profile that provides the best fit of the observed cluster-density profile. We therefore provide the central density and the core radius for each cluster.

  • We derived the MF of each cluster and constrained the MF slope, αMF, using a single power-law function. We find that the value of αMF depends on cluster age and ranges from ∼2.5 in ∼100 to 200 Myr-old clusters to less than 1.0 in clusters older than ∼2 Gyr. This finding is consistent with the recent conclusion by Ebrahimi et al. (2022) based on 15 nearby open clusters. The MFs are used to estimate the total cluster mass.

  • Visual inspection of the MFs of several clusters shows that stars in different mass intervals define MFs with different slopes. When we combine the MFs of all clusters together, we find that stars less massive than one solar mass are well fitted with a power-law function with α ∼ 1.5, whereas the MF of the other stars is consistent with α ∼ 2.5.

  • We derived the fraction of BSSs with respect to the MS stars in the luminosity interval between GRP = 0.0 and 1.4 mag from the turn-off. We confirm that clusters younger than ∼300 Myr show no evidence for BSSs (see Rain et al. 2021). In the clusters with BSSs, we find a mild correlation between the fraction of BSSs and cluster age and no correction with cluster mass.

Our main results on binaries can be summarized as follows:

  • When we analyze the entire cluster field, the fraction of binaries ranges from ∼15% in NGC 2669 to more than 60% in Haffner 26. There is no significant correlation between the binary fractions and either cluster age or cluster mass. A similar conclusion can be extended to the binaries in the core. Our results on open clusters differ from what was found in GCs, where there is a strong anti-correlation between the fraction of MS binaries in the core and the cluster mass (Milone et al. 2012, 2016).

  • There is no correlation between the fractions of binaries and BSSs. However, when we consider eight of nine star clusters with high central density (μ0 > 1 star pc−2), the fraction of BSSs correlates with the fraction of binaries.

  • There is no evidence for relations between the fraction of binaries of most clusters and either the mass of the primary star or the mass ratio.

  • The radial distribution of binaries significantly changes from one cluster to another and depends on cluster age. The binaries of the dynamically young clusters, that is, those with log(Age/trh) < 0.2, show a nearly flat radial distribution, and the distribution of binaries in intermediate-age clusters, 0.2 < log(Age/trh) < 0.75, exhibits hints of a double peak. On the contrary, the binaries of dynamically older clusters are more centrally concentrated than single stars, which is similar to what is observed in most GCs (Milone et al. 2012, 2016). This trend is also corroborated by the correlations between core binary fraction and both cluster age and dynamical age, and the absence of any correlation between total binary fraction and either cluster age or dynamical age. Our findings are qualitatively consistent with the theoretical predictions by Geller et al. (2013) and Sollima (2008). These authors used N-body simulations to model star clusters and binary evolution, and predict different radial distributions for binary stars of clusters with different ages. The fraction of binaries in young clusters (ages smaller than ∼300 Myr) would increase towards the external regions, while the binaries stars clusters older than ∼1 Gyr would exhibit an opposite trend. The clusters with intermediate ages would exhibit binaries with a double-peaked radial distribution.

  • We find signs of a mild correlation between total binary fraction and central stellar density (ρ = 0.54). Dividing clusters into three dynamical age groups, we find that this correlation is more evident in younger clusters, where the Spearman’s rank correlation coefficient reaches 0.63 for young clusters and 0.59 for intermediate-age clusters. We find no sign of correlation in dynamically old clusters. Moreover, while there is no correlation between cluster mass and binary fraction when considering the whole sample of clusters, dynamically young clusters exhibit hints of a mild correlation, that is, ρ = 0.45.

The results discussed in this work provide important constraints to the properties and evolution of Galactic OCs. In particular, we derived the properties of stellar binary systems and blue stragglers and connect them to the structural and physical parameters of the host cluster.


1

The procedure adopted to infer cluster radius and cluster mass is discussed in Sects. 3.1 and 3.2.

Acknowledgments

We thank the anonymous referee for various suggestions that improved the quality of the manuscript. This work has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research innovation programme (Grant Agreement ERC-StG 2016, No. 716082 ‘GALFOR’, PI: Milone, http://progetti.dfa.unipd.it/GALFOR) and from the European Union’s Horizon 2020 research and innovation programme under the Marie SklodowskaCurie Grant Agreement No. 101034319 and from the European Union – NextGenerationEU, beneficiary: Ziliotto. A.P.M., M.T., and E.D. acknowledge support from MIUR through the FARE project R164RM93XW SEMPLICE (PI: Milone). A.P.M. and E.D. have been supported by MIUR under PRIN program 2017Z2HSMF (PI: Bedin). A.F.M. acknowledges support from the project “Understanding the formation of globular clusters with their multiple stellar generations”. E.V. acknowledges support from NSF grant AST-2009193. S.J. acknowledges support from the NRF of Korea (2022R1A2C3002992, 2022R1A6A1A03053472).

References

  1. Baumgardt, H., & Hilker, M. 2018, MNRAS, 478, 1520 [Google Scholar]
  2. Bellazzini, M., Fusi Pecci, F., Messineo, M., Monaco, L., & Rood, R. T. 2002, AJ, 123, 1509 [NASA ADS] [CrossRef] [Google Scholar]
  3. Bolte, M. 1992, ApJS, 82, 145 [NASA ADS] [CrossRef] [Google Scholar]
  4. Bossini, D., Vallenari, A., Bragaglia, A., et al. 2019, A&A, 623, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Boubert, D., & Everall, A. 2020, MNRAS, 497, 4246 [NASA ADS] [CrossRef] [Google Scholar]
  6. Cantat-Gaudin, T., Jordi, C., Vallenari, A., et al. 2018, A&A, 618, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Cantat-Gaudin, T., Anders, F., Castro-Ginard, A., et al. 2020, A&A, 640, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Casagrande, L., & VandenBerg, D. A. 2018, MNRAS, 479, L102 [NASA ADS] [CrossRef] [Google Scholar]
  9. Clark, L. L., Sandquist, E. L., & Bolte, M. 2004, AJ, 128, 3019 [Google Scholar]
  10. Cordoni, G., Milone, A. P., Marino, A. F., et al. 2018, ApJ, 869, 139 [CrossRef] [Google Scholar]
  11. Dias, W. S., Monteiro, H., Moitinho, A., et al. 2021, MNRAS, 504, 356 [NASA ADS] [CrossRef] [Google Scholar]
  12. Dias, W. S., Monteiro, H., Caetano, T., et al. 2022, MNRAS, 512, 4464 [CrossRef] [Google Scholar]
  13. Ebrahimi, H., Sollima, A., & Haghi, H. 2022, MNRAS, 516, 5637 [NASA ADS] [CrossRef] [Google Scholar]
  14. Elson, R. A. W., Fall, S. M., & Freeman, K. C. 1987, ApJ, 323, 54 [Google Scholar]
  15. Everall, A., & Boubert, D. 2022, MNRAS, 509, 6205 [Google Scholar]
  16. Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Gaia Collaboration (Brown, A. G. A., et al.) 2021, A&A, 649, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Geller, A. M., de Grijs, R., Li, C., & Hurley, J. R. 2013, ApJ, 779, 30 [NASA ADS] [CrossRef] [Google Scholar]
  19. Goudfrooij, P., Girardi, L., Kozhurina-Platais, V., et al. 2014, ApJ, 797, 35 [CrossRef] [Google Scholar]
  20. Ji, J., & Bregman, J. N. 2015, ApJ, 807, 32 [Google Scholar]
  21. King, I. 1962, AJ, 67, 471 [Google Scholar]
  22. Malofeeva, A. A., Seleznev, A. F., & Carraro, G. 2022, AJ, 163, 113 [NASA ADS] [CrossRef] [Google Scholar]
  23. Malofeeva, A. A., Mikhnevich, V. O., Carraro, G., & Seleznev, A. F. 2023, AJ, 165, 45 [NASA ADS] [CrossRef] [Google Scholar]
  24. Marigo, P., Girardi, L., Bressan, A., et al. 2017, ApJ, 835, 77 [Google Scholar]
  25. Milone, A. P., Bedin, L. R., Piotto, G., & Anderson, J. 2009, A&A, 497, 755 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Milone, A. P., Piotto, G., Bedin, L. R., et al. 2012, A&A, 540, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Milone, A. P., Marino, A. F., Bedin, L. R., et al. 2016, MNRAS, 455, 3009 [CrossRef] [Google Scholar]
  28. Niederhofer, F., Georgy, C., Bastian, N., & Ekström, S. 2015, MNRAS, 453, 2070 [NASA ADS] [CrossRef] [Google Scholar]
  29. Rain, M. J., Ahumada, J. A., & Carraro, G. 2021, A&A, 650, A67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Richer, H. B., Fahlman, G. G., Brewer, J., et al. 2004, AJ, 127, 2771 [NASA ADS] [CrossRef] [Google Scholar]
  31. Romani, R. W., & Weinberg, M. D. 1991, ApJ, 372, 487 [NASA ADS] [CrossRef] [Google Scholar]
  32. Rubenstein, E. P., & Bailyn, C. D. 1997, ApJ, 474, 701 [NASA ADS] [CrossRef] [Google Scholar]
  33. Sharma, S., Pandey, A. K., Ogura, K., et al. 2008, AJ, 135, 1934 [NASA ADS] [CrossRef] [Google Scholar]
  34. Sollima, A. 2008, MNRAS, 388, 307 [NASA ADS] [CrossRef] [Google Scholar]
  35. Sollima, A. 2019, MNRAS, 489, 2377 [NASA ADS] [CrossRef] [Google Scholar]
  36. Sollima, A., Beccari, G., Ferraro, F. R., Fusi Pecci, F., & Sarajedini, A. 2007, MNRAS, 380, 781 [NASA ADS] [CrossRef] [Google Scholar]
  37. Sollima, A., Carballo-Bello, J. A., Beccari, G., et al. 2010, MNRAS, 401, 577 [NASA ADS] [CrossRef] [Google Scholar]
  38. Spitzer, L. 1987, Dynamical Evolution of Globular Clusters (Princeton: Princeton University Press) [Google Scholar]
  39. Vaidya, K., Panthi, A., Agarwal, M., et al. 2022, MNRAS, 511, 2274 [NASA ADS] [CrossRef] [Google Scholar]
  40. Zhao, B., & Bailyn, C. D. 2005, AJ, 129, 1934 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Open cluster CMDs

In the following, we show the differential reddening free GRP vs. GBP − GRP CMDs of the 78 analyzed open clusters. The CMDs include only cluster members, and are sorted according to the age of the cluster, indicated in the top inset in gigayears.

thumbnail Fig. A.1.

GRP vs. GBP − GRP CMDs corrected for differential reddening whenever needed. The clusters are sorted according to age (from Dias et al. 2021): Pozzo 1, ASCC 32, NGC 2547, Haffner 13, Trumpler 21, Trumpler 29, NGC 4609, NGC 6242, Collinder 394, NGC 4852, NGC 6087, NGC 2302, NGC 5617, Pismis 4, NGC 7243, Melotte 22, NGC 2353, Collinder 258, NGC 2448, and NGC 2546.

thumbnail Fig. A.2.

Same as Fig. A.1, but for the clusters NGC 5316, NGC 7063, NGC 2910, NGC 2422, Alessi 12, NGC 5168, NGC 6709, Lynga 2, IC 4725, NGC 2669, NGC 2396, NGC 6025, Collinder 463, Stock 12, IC 2488, NGC 1750, NGC 3114, NGC 2925, NGC 6152, and NGC 1912.

thumbnail Fig. A.3.

Same as Figure A.1, but for the clusters NGC 2516, NGC 6705, NGC 2287, NGC 6494, NGC 2374, NGC 1528, NGC 3532, NGC 2323, NGC 2482, IC 2714, NGC 7209, NGC 2567, NGC 2658, NGC 2358, Haffner 26, NGC 2281, NGC 2099, NGC 1664, NGC 2447, and NGC 2632.

thumbnail Fig. A.4.

Same as Figure A.1, but for the clusters NGC 2527, NGC 6811, NGC 2423, NGC 5822, NGC 2818, NGC 6940, NGC 1245, NGC 1817, NGC 2360, Tombaugh 1, Melotte 71, NGC 6208, NGC 2420, Haffner 22, NGC 6819, NGC 2682, NGC 188, and NGC 6791.

Appendix B: Individual binary-star distributions

Here, we show the individual binary distributions as a function of mass ratio (q, Fig. B.1), primary star mass (Mprim, Fig. B.2), and distance from the cluster center (B.3).

thumbnail Fig. B.1.

Mass-ratio distribution of binary stars for the 78 open clusters studied here. Clusters are sorted by name.

thumbnail Fig. B.2.

Mass distribution of binary stars for the 78 open clusters studied here.

thumbnail Fig. B.3.

Radial distribution of binary stars for the 78 open clusters studied here.

Appendix C: Correlations with clusters’ structural and physical properties

Figures C.1 and C.2 show the correlations between total and core binary fractions, respectively. We reiterate here that core binary fractions are computed only for clusters with more than 500 stars.

thumbnail Fig. C.1.

Correlations between total binary fraction and the physical parameters of clusters. The marker color and line color of each panel are indicative of the correlations as shown in the right color bar. The top insets indicate the considered physical parameters, i.e., the x-axis, while the values of the Spearman rank coefficients and p-values are indicated in the bottom inset. Variables μ0, Mcl, trh, Age/trh, Age are in logarithm.

thumbnail Fig. C.2.

Correlations between binary fraction in the core and the physical parameters of clusters. As discussed in the text, we limit the analysis of core binary fractions to clusters with more than 400 stars. The marker color and line color of each panel are indicative of the correlations as shown in the right color bar. The top insets indicate the considered physical parameters, i.e., the x-axis, while the values of the Spearman rank coefficients and p-values are indicated in the bottom inset. Variables μ0, Mcl, trh, Age/trh, Age are in logarithm.

Table C.1.

Cluster properties and density profiles. ID, right ascension, declination, proper motions, parallax, distance modulus, log(age), reddening, number of stars, differential reddening flag, core radius (arcmin) and associated error, core radius (parsecs) and associated error, central density (stars/arcmin2) and (stars/parsec2) and associated errors.

Table C.2.

Cluster structural parameters and binary fractions. IF, MF slope, low-mass MF slope, high-mass MF slope with relative uncertainties, cluster mass, half-mass radius, hal-fmass relaxation time minimum binaries mass-ratio q, binary fractions and errors, blue-straggler fraction, and number.

All Tables

Table C.1.

Cluster properties and density profiles. ID, right ascension, declination, proper motions, parallax, distance modulus, log(age), reddening, number of stars, differential reddening flag, core radius (arcmin) and associated error, core radius (parsecs) and associated error, central density (stars/arcmin2) and (stars/parsec2) and associated errors.

Table C.2.

Cluster structural parameters and binary fractions. IF, MF slope, low-mass MF slope, high-mass MF slope with relative uncertainties, cluster mass, half-mass radius, hal-fmass relaxation time minimum binaries mass-ratio q, binary fractions and errors, blue-straggler fraction, and number.

All Figures

thumbnail Fig. 1.

Properties of analyzed clusters. Top: histogram distribution of the ages (in logarithm scale) of the 78 open clusters analyzed here. Bottom: histogram distribution of the number of cluster members. See the text for details.

In the text
thumbnail Fig. 2.

Cluster members selection. Procedure carried out to select the probable NGC 2099 cluster members located within ∼30 arcmin from the cluster center. Panel a: stellar proper motions relative to the average cluster motion against the G magnitude. The blue line separates stars with cluster-like proper motions from the remaining stars. Panel b: stellar G magnitude as a function of parallax. Stars with cluster-like parallaxes are located between the two blue lines. The vector point diagram of proper motions is plotted in panel c, whereas panel d shows the spatial stellar distribution. Finally, panel e illustrates the G vs. GBP − GRP observed CMD, together with the best-fit isochrone. Likely cluster members are represented with azure crosses, while brown crosses indicate stars with cluster-like proper motions and parallaxes in the reference field.

In the text
thumbnail Fig. 3.

Differential reddening correction. Comparison of the GRP vs. GBP − GRP CMD before (panel a) and after the correction for differential reddening (panel b). The differential reddening map is shown in panel c, color-coded according to the top color bar.

In the text
thumbnail Fig. 4.

Density profile of NGC 2099. Azure points with error bars represent the observed stellar density profile. The black solid line indicates the best-fitting EEF (Elson et al. 1987) profile, while the corresponding uncertainty is indicated by the gray shaded area. The inferred core radius is indicated with the black solid vertical line, while the shaded region is its uncertainty. The best-fitting parameters are quoted in the bottom-left inset. The dashed line shows the best-fitting King (King 1962) profile.

In the text
thumbnail Fig. 5.

Mass function determination. Left panel: GRP vs. GBP − GRP CMD of NGC 2099. The magnitude intervals used to derive the MF are highlighted with different colors, while the median mass of each bin is indicated on the right. Right panel: MF for the stars plotted in the left panel, i.e., log(NM) vs. stellar mass. The blue and azure straight lines represent the best-fit straight lines of low- (M ≤ 1 M) and high-mass stars (M > 1 M), respectively. The red line is obtained by fitting the whole MF. The slopes are quoted in the bottom-left legend, while the total stellar mass of the cluster is indicated in the top-right inset.

In the text
thumbnail Fig. 6.

Mass functions properties. Panel a: distribution of the overall MF slope inferred for the 78 analyzed clusters. Panel b: MF slope α vs. cluster age. Panel c: cluster age vs. explored stellar mass range in MS stars.

In the text
thumbnail Fig. 7.

A common mass function. Panels a–c: normalized stacked MFs for young (log(Age/trh)≤0.2), intermediate (0.2 < log(Age/trh) < 0.75), and old clusters (log(Age/trh)≥0.75). The normalization choice is described in Sect. 3.2. Blue and azure solid lines represent the straight line best-fit of the low-mass (M ≤ 1 M) and high-mass (M > 1 M) stars, respectively. The derived slope is indicated in the bottom-left legends. Panel d: same as panels a–c, but for the whole sample of analyzed clusters.

In the text
thumbnail Fig. 8.

Binary stars selection. Panel a: illustration of the binary-selection procedure. The region enclosed by the azure solid line, dubbed region A, contains all MS stars and binary stars with a primary mass between 0.80 and 1.31 M. The gray-shaded region, i.e., region B, defines the locus of binary-star candidates with a mass ratio greater than 0.6, respectively. Finally, the dashed black line and red solid lines represent the fiducial line of MS stars and binary stars with mass ratio q = 0.6, as indicated in the top-right legend. Panel b: MS and binary stars identified in the cluster field of NGC 2099, marked with red and blue crosses, respectively. Panel c: same as panel b but for stars located in the reference field, well outside the cluster radius. Stars that fall into region A are marked with black crosses.

In the text
thumbnail Fig. 9.

Relation between binary fractions and host clusters. Panels a–c: total binary fraction as a function of cluster dynamical age, i.e., age normalized over the half-mass relaxation time, total cluster mass, and central density in stars pc−2. Azure lines and shaded regions indicate the straight line that best fits the data and its 1σ dispersion. The values of Spearman’s rank coefficient (ρ) and the corresponding p-values are indicated in the top-left corner of each plot. Panels b1–b3: same as panel b, but with clusters grouped into three dynamical age groups, defined as young, intermediate, and old clusters. Panels c1–c3: same as panel c, for the three groups of clusters.

In the text
thumbnail Fig. 10.

Procedure used to estimate the fraction of BSSs in NGC 188. Cluster and reference-field stars are marked with gray circles and brown crosses, respectively. The gray dashed line represents the ZAMS. The gray solid line is the fiducial line of the MS and the SGB shifted to the blue by four times the observational color error (σ). The fiducial line of MS-MS and SGB-SGB binaries, shifted by four times σ to the red is represented by the red line. Region C of the CMD is indicated as an azure-shaded area, while region D, marked with the gray shade, is enclosed by gray and red solid lines and by the two horizontal black dashed lines. Cluster field stars in regions C and D of the CMD are marked with azure solid and black-shaded circles, respectively.

In the text
thumbnail Fig. 11.

Blue stragglers properties. Fraction of BSSs against cluster age (panel a), cluster mass (panel b), and cluster metallicity.

In the text
thumbnail Fig. 12.

Fraction of blue stragglers vs. binary fraction. Red colors indicate the clusters with central density μ0 < 10 stars pc−2, whereas the remaining clusters are colored blue.

In the text
thumbnail Fig. 13.

Binary stars mass-ratio as a function of mass-ratio in NGC 2099. Reproduction of the CMD of NGC 2099 plotted in Fig. 8b (panel a). The shaded regions (B1, B2, and B3) of the CMD comprise binaries with different mass ratios. The candidate binaries are marked with white crosses. Panel b shows the mass-ratio distribution of binaries in NGC 2099. The χ2 and the probability that the observed binary fractions come from a flat mass-ratio distribution are quoted in panel b.

In the text
thumbnail Fig. 14.

Binary fraction mass-ratio distribution. The binary fraction in each mass-ratio bin has been normalized over the width of the bin, in q, and the total fraction of binaries extrapolated assuming a flat mass ratio distribution. Each cell is color coded according to the number of observations in that particular bin, as specified in the top right color bar.

In the text
thumbnail Fig. 15.

Binary stars as a function of primary stars Mprim in NGC 2099. The adopted regions are colored differently, as defined in the main text. The mass distribution of candidate binary stars with q ≥ _lim is shown in panel b.

In the text
thumbnail Fig. 16.

Binary fraction primary star distribution. The azure line and the shaded region represent the median trend and the 1σ confidence interval.

In the text
thumbnail Fig. 17.

Binary stars radial distribution. Binary fraction with q ≥ qlim normalized to the total binary fraction, as a function radial distance from the cluster center in core radius units. The azure lines mark the median trends, while the shaded areas correspond to the 1σ dispersion. Panels a–c include clusters in the three different dynamical age intervals defined in the main text (see Sect. 3.2), while panel d refers to all clusters. Color coding and age bins are the same as in Fig. 6.

In the text
thumbnail Fig. 18.

Comparison between OCs and GCs. Total fraction of binaries in the core of open clusters (aqua circles) and GCs (red crosses) as a function of the logarithm of cluster age (panel a), dynamical age (panel b), and cluster mass (panel c). The GC binary fractions are taken from Milone et al. (2012, 2016), while GC masses are from Baumgardt & Hilker (2018).

In the text
thumbnail Fig. A.1.

GRP vs. GBP − GRP CMDs corrected for differential reddening whenever needed. The clusters are sorted according to age (from Dias et al. 2021): Pozzo 1, ASCC 32, NGC 2547, Haffner 13, Trumpler 21, Trumpler 29, NGC 4609, NGC 6242, Collinder 394, NGC 4852, NGC 6087, NGC 2302, NGC 5617, Pismis 4, NGC 7243, Melotte 22, NGC 2353, Collinder 258, NGC 2448, and NGC 2546.

In the text
thumbnail Fig. A.2.

Same as Fig. A.1, but for the clusters NGC 5316, NGC 7063, NGC 2910, NGC 2422, Alessi 12, NGC 5168, NGC 6709, Lynga 2, IC 4725, NGC 2669, NGC 2396, NGC 6025, Collinder 463, Stock 12, IC 2488, NGC 1750, NGC 3114, NGC 2925, NGC 6152, and NGC 1912.

In the text
thumbnail Fig. A.3.

Same as Figure A.1, but for the clusters NGC 2516, NGC 6705, NGC 2287, NGC 6494, NGC 2374, NGC 1528, NGC 3532, NGC 2323, NGC 2482, IC 2714, NGC 7209, NGC 2567, NGC 2658, NGC 2358, Haffner 26, NGC 2281, NGC 2099, NGC 1664, NGC 2447, and NGC 2632.

In the text
thumbnail Fig. A.4.

Same as Figure A.1, but for the clusters NGC 2527, NGC 6811, NGC 2423, NGC 5822, NGC 2818, NGC 6940, NGC 1245, NGC 1817, NGC 2360, Tombaugh 1, Melotte 71, NGC 6208, NGC 2420, Haffner 22, NGC 6819, NGC 2682, NGC 188, and NGC 6791.

In the text
thumbnail Fig. B.1.

Mass-ratio distribution of binary stars for the 78 open clusters studied here. Clusters are sorted by name.

In the text
thumbnail Fig. B.2.

Mass distribution of binary stars for the 78 open clusters studied here.

In the text
thumbnail Fig. B.3.

Radial distribution of binary stars for the 78 open clusters studied here.

In the text
thumbnail Fig. C.1.

Correlations between total binary fraction and the physical parameters of clusters. The marker color and line color of each panel are indicative of the correlations as shown in the right color bar. The top insets indicate the considered physical parameters, i.e., the x-axis, while the values of the Spearman rank coefficients and p-values are indicated in the bottom inset. Variables μ0, Mcl, trh, Age/trh, Age are in logarithm.

In the text
thumbnail Fig. C.2.

Correlations between binary fraction in the core and the physical parameters of clusters. As discussed in the text, we limit the analysis of core binary fractions to clusters with more than 400 stars. The marker color and line color of each panel are indicative of the correlations as shown in the right color bar. The top insets indicate the considered physical parameters, i.e., the x-axis, while the values of the Spearman rank coefficients and p-values are indicated in the bottom inset. Variables μ0, Mcl, trh, Age/trh, Age are in logarithm.

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.