🔭 astrophysics

Stellar Multiplicity in Open Clusters: Investigating Binary Fractions and Their Relationship with Cluster Properties

By analyzing 773 open clusters, this study reveals that while binary fractions are largely independent of cluster age and evolution, they exhibit an anticorrelation with metallicity and show that binary systems are preferentially concentrated in cluster centers with mass ratios decreasing toward the outskirts, suggesting that multiplicity is primarily established during early formation rather than shaped by subsequent dynamical evolution.

Original authors: Ruan P. Alves, Hektor Monteiro, Wilton S. Dias, Gabriel R. Hickel

Published 2026-07-17

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Ruan P. Alves, Hektor Monteiro, Wilton S. Dias, Gabriel R. Hickel

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the night sky not as a sea of lonely, solitary suns, but as a bustling cosmic dance floor where most stars are actually holding hands. In the grand theater of astronomy, a "binary system" is simply a pair of stars orbiting a common center, like a spinning couple in a ballroom, while a "multiple system" is a whole group of stars dancing together. These pairs aren't just romantic; they are the heavyweights of the universe's dynamics. When stars form, they often pop out in pairs, and these relationships dictate how stars evolve, how planets are born, and even how spectacular explosions happen later in life. But here's the big question: does the environment where these stars are born change how often they pair up? Is a star more likely to find a partner in a young, chaotic nursery or an old, settled neighborhood? To answer this, astronomers look at "open clusters"—groups of stars born from the same cloud of gas and dust, like a massive family reunion where everyone shares the same birthday and genetic makeup. By studying these clusters, scientists can see if the "rules of pairing" change based on the cluster's age, its chemical makeup, or its location in the galaxy.

This paper dives deep into that very question, acting like a cosmic detective examining 773 different open clusters to see how many stars are dancing in pairs. The researchers, using a massive dataset of star positions and brightness, built a detailed map to identify which stars are single and which are part of a binary system. They didn't just count the pairs; they looked at how heavy the partners are compared to each other (the "mass ratio") and where these pairs are hanging out within their clusters.

Here is what they found, and it's a bit more surprising than you might expect. First, they discovered that the "binary fraction"—the percentage of stars that are in pairs—doesn't seem to change much as a cluster gets older. Whether a cluster is a fresh 4-million-year-old baby or a 6-billion-year-old veteran, the number of pairs stays roughly the same. This suggests that the decision to form a pair happens very early in a star's life, and once the stars are born, the cluster's aging process doesn't really break them up or create new ones. The paper explicitly argues against the idea that the cluster's age or its evolutionary stage is a major factor in changing these numbers.

However, the chemical makeup of the cluster does seem to matter. The study found a weak but noticeable link: clusters with lower metallicity (which means they have fewer heavy elements like iron) tend to have slightly more binary stars. It's as if the "ingredients" of the birth cloud influence how likely stars are to pair up.

The researchers also looked at the "mass ratios," or how similar the two stars in a pair are in weight. They found that stars love to pair up with partners of almost equal weight (a ratio near 0.9), creating "twin" systems. Interestingly, this preference for twins is much stronger in young clusters. In older clusters, the number of these perfect twins drops off, suggesting that over billions of years, the chaotic motion of the cluster might be breaking up these delicate, equal-weight pairs, leaving behind more mismatched couples.

Another fascinating discovery is about where these pairs live. Binary stars seem to prefer the "VIP section" near the center of the cluster, while single stars are more likely to be found wandering in the outer edges. As you move further away from the center, the average weight difference between partners gets bigger. This concentration in the middle might actually be a safety mechanism, protecting the pairs from being flung out of the cluster by gravitational tugs.

Finally, the team looked at how likely a star is to be the "leader" (the primary) or the "follower" (the companion) in a pair. They found that the heavier a star is, the more likely it is to be the primary component. For lighter stars (less than the mass of our Sun), the chance of being in a binary system increases as the star gets heavier. But for the massive stars, the odds stay steady at a high 60% to 70%, meaning almost all the heavy hitters are part of a team.

In short, this study suggests that the "pairing up" of stars is mostly a decision made at birth, influenced by the chemical ingredients of the cloud they form in, rather than by the drama of their later years. While the cluster's age doesn't seem to change the total number of pairs, the passage of time does seem to shuffle the deck, breaking up some of the perfect twin pairs and leaving the survivors clustered tightly in the center of the group. The paper doesn't claim to have solved every mystery—there are still questions about how exactly these pairs form and how rotation might confuse the data—but it provides a massive, clear picture of how stellar relationships hold up across the galaxy.

Technical Summary: Stellar Multiplicity in Open Clusters

Problem Statement
Stellar multiplicity is a fundamental component of stellar formation and evolution, significantly influencing the dynamical history of open clusters (OCs) and the formation of planetary systems. While it is established that most stars form in multiple systems, the extent to which environmental properties—such as cluster age, metallicity, and dynamical state—influence the occurrence and characteristics of binary systems remains a subject of investigation. Previous studies have reported varying binary fractions ( $f_b$ ) across different OC samples, often limited by small sample sizes or heterogeneous methodologies. This work aims to investigate how binary fractions and mass ratios ( $q = M_2/M_1$ ) vary across a large, homogeneous sample of OCs and to determine their relationship with fundamental cluster parameters.

Methodology
The study utilizes a sample of 773 open clusters derived from the catalog of A. Almeida et al. (2023), which itself is based on Gaia DR3 data and the W. S. Dias et al. (2021) cluster catalog. The clusters span a wide parameter space, with distances ranging from 0.14 to 6.87 kpc, ages from 4 Myr to 6 Gyr, and metallicities between −0.37 and 0.43 dex.

Binary Identification and Mass Estimation: The authors adopted the methodology of A. Almeida et al. (2023), which employs a Monte Carlo approach to compare observed photometric data with synthetic clusters generated using Padova PARSEC isochrones. This process classifies member stars as single or binary and estimates the masses of primary and companion stars. Systems with companion mass uncertainties $\delta M_2/M_2 > 0.5$ were excluded.
Correction for Observational Bias: The authors identified systematic biases where binary fractions were artificially inflated in distant clusters (due to the preferential detection of massive stars) and in highly reddened clusters (due to CMD broadening). To mitigate this, three correction methods were tested: a mean-based correction, a moving-average (MA) correction, and a similarity-based correction. The mean-based correction was selected as the reference, as it most effectively removed the dependence of $f_b$ on extinction ( $A_V$ ) and distance.
Statistical Analysis: The study calculated the total binary fraction ( $f_b$ ) and the high mass-ratio binary fraction ( $f_{b, q>0.5}$ ). Probabilistic definitions were established for stars being in a binary system ( $P_{BS}$ ), acting as a primary ( $P_{Prim}$ ), or acting as a secondary ( $P_{Comp}$ ). Uncertainties were estimated by combining systematic errors (via Monte Carlo resampling) and statistical errors (binomial statistics).
Dynamical Context: The radial distribution of binaries was analyzed relative to the half-mass radius ( $R_h$ ), and the dynamical state of clusters was quantified using the parameter $\tau = \text{age}/t_{relax}$ .

Key Results

Binary Fractions: The estimated total binary fraction for the sample ranges from 0.10 to 0.85, with a median of 0.40. For high mass-ratio systems ( $q > 0.5$ ), the fraction ranges from 0.06 to 0.79, with a median of 0.36. The authors note these values are lower limits, as the analysis is restricted to systems with $q > 0.34$ .
Mass Ratios: The distribution of mass ratios peaks near $q \approx 0.9$ , indicating a preference for nearly equal-mass binaries. A notable feature is a significant decrease in the mean mass ratio ( $\langle q \rangle$ ) for primary masses between 1 and 3 $M_\odot$ , reaching a minimum of $\langle q \rangle = 0.65$ at $M_1 \approx 2.75 M_\odot$ . This dip is observed in both young and intermediate-age clusters, suggesting it is not solely an evolutionary artifact of very young populations.
Dependence on Cluster Properties:
- Age and Evolution: No strong dependence of the binary fraction on cluster age or evolutionary stage was found. However, the mass-ratio distributions of young and old clusters are statistically distinct. Younger clusters exhibit a more pronounced peak at $q \approx 0.9$ and a higher occurrence of "twin" binaries compared to older clusters, where the distribution for $q > 0.6$ is flatter.
- Metallicity: A weak anticorrelation was identified between metallicity and binary fraction, consistent with theoretical expectations that metal-poor environments favor molecular cloud fragmentation. However, the authors caution that this result is affected by small sample sizes at the extremes of the metallicity distribution and potential methodological biases.
- Mass Dependence: The probability of a star being a primary component increases with stellar mass. Conversely, the probability of a star belonging to a binary system increases with mass for $M < 1 M_\odot$ and stabilizes between 60% and 70% for higher masses.
Spatial Distribution: Binary systems are preferentially located in the inner regions of clusters compared to single stars. Furthermore, the mean mass ratio decreases with increasing distance from the cluster center.

Significance and Claims
The paper claims that the binary fraction in open clusters appears to be established at early times and remains largely unaffected by subsequent cluster dynamical evolution, as evidenced by the lack of strong correlation between $f_b$ and cluster age or dynamical state. The concentration of binaries in the inner regions suggests that these systems are less susceptible to dynamical escape than single stars. The distinct mass-ratio distributions between young and old clusters are attributed to dynamical interactions over time, where wide, low-mass-ratio binaries are more likely to be disrupted, while twin binaries remain stable.

The study provides a large-scale, homogeneous statistical analysis of stellar multiplicity in OCs, offering a robust dataset that refines previous estimates and highlights the complex interplay between initial formation conditions and dynamical processing. The authors emphasize that their results serve as a lower limit to the true binary fraction due to the exclusion of low-mass companions ( $q < 0.34$ ).

Technical Summary: Stellar Multiplicity in Open Clusters

More like this