Origin of open clusters revealed by the evolution of the m_max$-$M_ecl relation

By analyzing Gaia DR3 data and comparing it with N-body simulations of gas expulsion and subcluster coalescence, the study concludes that the observed evolution of the most massive star versus cluster mass relation supports subcluster coalescence as the dominant formation pathway for open clusters.

The Gist

Here is an explanation of the paper using simple language and creative analogies.

The Big Question: How Do Star Clusters Grow Up?

Imagine a star cluster as a giant, bustling school. When the school first opens (the "birth" of the cluster), it has a specific number of students (stars) and a specific mix of sizes (masses).

A long-standing debate in astronomy is about how these "schools" are built:

1. The Random Lottery: Do stars form randomly, like drawing names out of a hat? Maybe you get a few giants and many small kids just by chance.
2. The Organized System: Is there a strict rule? Does the size of the school determine exactly how big the biggest student must be? (e.g., "If the school has 1,000 students, the biggest one must be 50kg.")

This paper investigates a specific rule called the $m_{max} - M_{ecl}$ relation. In plain English, this is the relationship between the mass of the biggest star in a cluster and the total mass of the cluster when it was born.

The Mystery: The "Missing" Giants

The authors looked at a massive database of star clusters (using data from the Gaia space telescope). They found something strange:

• Young Clusters (Newborns): The biggest stars fit the "Organized System" rule perfectly. The total mass predicts the biggest star's size accurately.
• Old Clusters (Teenagers/Adults): As the clusters get older (older than 5 million years), the rule breaks down. The biggest stars are often smaller than the rule predicted.

It's as if you walked into a 10-year-old school and found that the "biggest kid" was actually quite small, even though the school was huge. Where did the giants go?

The Investigation: Two Theories

To solve this mystery, the scientists ran supercomputer simulations (N-body simulations) to see how clusters evolve over time. They tested two main theories on how a cluster forms and changes:

Theory A: The "Single Giant" (Monolithic Formation)

Imagine a single, massive cloud of gas collapsing all at once to form one big school.

• The Process: The gas is expelled (blown away) by the stars' winds and radiation. This happens at different speeds (Fast, Moderate, or Slow).
• The Result: When gas leaves, the school loses students. The biggest stars are often the first to be kicked out or die because they are so heavy and unstable.
• The Problem: Even with different speeds of gas expulsion, these simulations didn't quite match the real data. The "biggest stars" in the simulations were still too heavy compared to what we see in the real universe.

Theory B: The "Merge" (Subcluster Coalescence)

Imagine a different scenario. Instead of one big school forming at once, imagine several small schools forming in the same neighborhood. Over time, they drift together and merge into one giant school.

• The Process: You have Cluster A, Cluster B, and Cluster C. Each has its own "biggest star." When they merge, the new "Big School" has a total mass equal to A+B+C.
• The Twist: The biggest star in the merged school is just the biggest star from the original small schools. It is not a new, super-massive star created by the merger.
• The Result: Because the total mass of the new school is huge (A+B+C), but the biggest star is only as big as the winner of the small schools, the ratio looks "off." The biggest star seems too small for the total mass.

The Verdict: The "Merge" Wins

The authors compared their simulations to the real observations and found a clear winner:

1. The "Merge" Scenario Matches Reality: The data shows that older clusters have "too small" biggest stars for their total mass. This happens perfectly in the Subcluster Coalescence model. When small clusters merge, the math works out exactly like the observations.
2. The "Single Giant" Scenario Fails: The single-cluster models (even with slow gas loss) kept predicting stars that were too massive compared to what we actually see.

The Analogy: The "Potluck Dinner"

Think of it like a Potluck Dinner:

• The "Single Giant" Theory: Imagine one person brings a huge platter of food. The size of the biggest dish on the table is directly related to how much food that one person brought.
• The "Merge" Theory: Imagine 10 different people bring small dishes. They all dump them onto one giant table.
  • The Total Mass of the table is huge (10 people worth of food).
  • But the Biggest Dish is still just the size of the biggest dish one person brought.
  • If you look at the table and try to guess the size of the biggest dish based on the total amount of food, you would be wrong. You'd expect a giant dish, but you only see a medium one.

The paper concludes that star clusters are like the Potluck Dinner. They don't form as one giant blob. Instead, they form as many small groups (subclusters) that eventually crash into each other and merge.

Why Does This Matter?

This discovery changes how we understand the birth of stars. It suggests that the universe is a bit more chaotic and collaborative than we thought. Stars aren't just born in isolation or in one giant explosion; they are born in small families that eventually grow up to become the massive clusters we see today.

In short: The "missing giants" aren't missing; they were never there to begin with. The clusters grew big by merging, not by growing a single giant star.

Technical Summary

Here is a detailed technical summary of the paper "Origin of open clusters revealed by the evolution of the $m_{\text{max}} - M_{\text{ecl}}$ relation".

1. Problem Statement

The paper addresses a fundamental question in astrophysics regarding the formation and evolution of star clusters: Do stars in a cluster form via stochastic (random) sampling from the Initial Mass Function (IMF), or is there a deterministic relationship between the total mass of the embedded cluster ($M_{\text{ecl}}$) and the mass of its most massive star ($m_{\text{max}}$)?

While a deterministic $m_{\text{max}} - M_{\text{ecl}}$ relation is well-established for young, embedded clusters ($<5$ Myr), it is unclear how this relation evolves as clusters age, lose gas, and disperse into open clusters (OCs). Specifically, the authors investigate whether the observed deviations in older clusters are due to standard dynamical evolution (mass loss) or if they imply a specific formation pathway, such as the coalescence of subclusters.

2. Methodology

The study employs a dual approach combining large-scale observational data analysis with extensive N-body simulations.

A. Observational Data

• Source: The Gaia DR3 open cluster catalog by Hunt & Reffert (2023, 2024).
• Sample Selection: The authors focused on high-quality bound open clusters (OCs) and unbound moving groups (MGs).
  • Criteria: Median color-magnitude diagram (CMD) classification $>0.5$ and Signal-to-Noise Ratio (S/N) $>5\sigma$.
  • Age Filter: Clusters with a 50th percentile age $>5$ Myr.
  • Mass Constraint: The uncertainty in the most massive star's mass must be $<50\%$.
• Data Processing:
  • Photometric masses were derived using PARSEC isochrones, corrected for selection effects and unresolved binaries.
  • To ensure consistency with simulations, the authors recalculated the "observed" cluster mass by restricting member stars to those within the theoretical tidal radius ($r_t$), rather than the catalog's observational radius.

B. N-Body Simulations

The authors utilized the PeTar code for dynamical evolution and McLuster for initial conditions. The simulations were divided into two main categories:

1. Individual Cluster Evolution:

   • Simulated clusters with initial stellar masses ranging from $300$to$10,000 M_\odot$.
   • Gas Expulsion Modes: Four modes were tested to model the removal of residual gas by stellar feedback:
     • Fast: Timescale $\tau_g$
     • Moderate 2: Timescale $2\tau_g$
     • Moderate 1: Timescale $5\tau_g$
     • Slow: Timescale $10\tau_g$
   • Initial Conditions: Plummer density profiles, fully mass-segregated ($S=1$), virial equilibrium ($Q=0.5$), and a canonical IMF. The initial $m_{\text{max}}$ was set according to the deterministic relation by Weidner et al. (2013).

2. Cluster Coalescence Simulations:

   • Modeled the merger of multiple subclusters formed within a single parental molecular cloud.
   • Scenarios:
     • "-fast": Subclusters initially at rest on the same plane.
     • "-fast-vd": Subclusters with spatial separation and velocity differences (following Larson's relation, $v \propto L^{0.5}$).
   • Setup: Subclusters were evolved individually for 1 Myr (expansion phase) before being merged to simulate the coalescence process.

3. Key Contributions

• Systematic Comparison: This is the first study to systematically compare the evolution of the $m_{\text{max}} - M_{\text{cluster}}$ relation in observations against both monolithic (single) cluster evolution and subcluster coalescence scenarios.
• Refined Tidal Radius Analysis: The authors introduced a rigorous method to align observational data with simulations by restricting observed cluster members to the theoretical tidal radius, ensuring a fair comparison of mass and $m_{\text{max}}$.
• Degeneracy Resolution: The study addresses the degeneracy between cluster mass, gas expulsion timescale, and formation channel, demonstrating that coalescence simulations produce distinct evolutionary tracks that better match observations than standard monolithic models.

4. Results

A. Evolution of the $m_{\text{max}} - M_{\text{cluster}}$ Relation

• Observations: As clusters age, the observed relation gradually deviates from the initial theoretical $m_{\text{max}} - M_{\text{ecl}}$ relation. For clusters older than 5 Myr, the most massive stars are significantly less massive than predicted by the initial relation for a given cluster mass. This manifests as age stratification.
• Individual Cluster Simulations:
  • The evolution of the relation follows a similar trajectory across all four gas expulsion modes, differing primarily in the speed of evolution.
  • Faster gas expulsion leads to rapid mass loss and a quicker deviation from the initial relation.
  • However, even with slow gas expulsion, individual cluster simulations tend to retain higher $m_{\text{max}}$ values for a given cluster mass compared to observations.
• Coalescence Simulations:
  • These simulations show a systematically lower $m_{\text{max}} - M_{\text{cluster}}$ relation compared to individual cluster simulations.
  • The merged cluster behaves dynamically like a single cluster undergoing slow/moderate gas expulsion, resulting in slower mass loss.
  • Crucially, the most massive star in a coalesced system is less massive than the most massive star in a monolithic cluster of the same total mass (because the total mass is distributed among multiple subclusters, each with its own $m_{\text{max}}$ limit).

B. Comparison of Observations and Simulations

• Best Fit: The coalescence simulations provide a significantly better match to the observational data than the individual cluster simulations.
  • Both observations and coalescence simulations exhibit lower $m_{\text{max}}$ values for a given cluster mass.
  • Observations show that for clusters $>5$ Myr, the most massive stars have already deviated from the initial relation. Individual cluster simulations often retain stars in the "initial relation" region for longer, whereas coalescence simulations show the deviation earlier and more consistently.
• Mass Loss Rates: The lower observed relation implies that open clusters lose mass more slowly than predicted by rapid gas expulsion models of monolithic clusters. Coalescence naturally reproduces this slower mass loss.

5. Significance and Conclusion

The paper concludes that the evolution of the $m_{\text{max}} - M_{\text{ecl}}$ relation strongly supports subcluster coalescence as the dominant pathway for open cluster formation.

• Implication for Star Formation: The results suggest that open clusters are not necessarily born as single, monolithic entities. Instead, they likely form through the hierarchical merging of smaller subclusters within a parent molecular cloud.
• Resolution of Discrepancies: This formation channel explains why observed open clusters have lower maximum stellar masses for their total mass compared to theoretical predictions based on random sampling or monolithic collapse.
• Consistency: These findings align with previous work by the authors (Zhou et al. 2024c, 2025) regarding the parameter space (mass, radius, age) of Milky Way open clusters.

In summary, the study provides robust evidence that the dynamical history of open clusters is dominated by the coalescence of substructures, fundamentally shaping the relationship between a cluster's total mass and its most massive stellar member.