An analytical approach to binary populations in globular clusters

This paper demonstrates through analytical calculations and N-body simulations that the low binary fractions observed in present-day globular clusters can be fully explained by the dynamical dissolution of primordial "soft" binaries, a process significantly influenced by stellar black holes, without requiring different initial star formation conditions compared to the solar neighborhood.

The Gist

The Big Mystery: Where Did All the Couples Go?

Imagine you walk into a crowded dance hall (our Solar Neighborhood). You look around and see that about half the people are dancing in pairs. It's a very social place; couples are everywhere.

Now, imagine walking into a different, much more chaotic dance hall called a Globular Cluster. These are ancient, super-dense balls of stars, packed tight like sardines in a can. When astronomers look inside these clusters, they are shocked: there are almost no couples. Only 1% to 10% of the stars are dancing in pairs. The rest are dancing solo.

The Question: Did these clusters start with fewer couples to begin with? Or did something happen to break them up?

The Paper's Answer: The "Soft" vs. "Hard" Dance

The authors of this paper propose a simple, elegant solution: The clusters started with the same number of couples as our neighborhood, but the chaotic environment broke up the "soft" ones.

To understand this, we need to define two types of dancing couples:

1. The "Hard" Couples (The Tight Embrace): These are pairs holding each other very tightly, spinning fast and close together. It takes a lot of energy to pull them apart. In the dance hall, they are stable and can survive the chaos.
2. The "Soft" Couples (The Loose Handshake): These are pairs holding hands loosely, standing far apart, and spinning slowly. They are "soft" because a gentle bump from a third dancer could easily knock them apart.

The Analogy of the Crowd:
Imagine a mosh pit (the dense cluster).

• If you are holding hands loosely with a friend (Soft), and a third person bumps into you, you and your friend will likely get separated. You become single dancers.
• If you are locked in a tight, spinning embrace (Hard), and someone bumps into you, you might spin faster, but you won't let go. In fact, the bump might make you hold on even tighter!

The Main Findings

The authors used math and computer simulations to prove three main things:

1. The "Soft" Couples Don't Last
In the beginning, the cluster was full of both tight and loose couples. Because the cluster is so crowded, the "Soft" couples (wide, loose orbits) were constantly getting bumped by other stars. Over millions of years, almost all of them were broken up. The "Hard" couples (tight orbits) survived. This explains why we see so few couples today: the loose ones were dissolved by the crowd.

2. The "Black Hole" Engine
You might ask: "Who is doing the bumping? Who has the energy to break up all those couples?"
The authors found that Stellar Black Holes are the secret agents.

• When stars die in these clusters, they often leave behind heavy black holes.
• Because they are heavy, they sink to the center of the cluster.
• These black holes act like a dynamical engine. They bounce around, colliding with the "Soft" couples and stealing their energy to break them apart.
• This is a good thing! If the black holes didn't do this, the cluster would collapse in on itself too quickly. The black holes are essentially "burning" the soft couples to keep the cluster stable for billions of years.

3. The "Birth Certificate" Check
The authors turned their logic around. They looked at the few couples that did survive in different clusters today. By counting how many survived, they could calculate how "crowded" the cluster must have been when it was born.

• The Result: They found that these ancient clusters were born with sizes and densities very similar to the massive, young star clusters we see forming today in the universe. This suggests that star formation rules are universal: whether a star is born in our quiet neighborhood or a chaotic ancient cluster, it starts with the same likelihood of finding a partner.

Why This Matters

For a long time, scientists debated whether the low number of couples in globular clusters meant they were born different (maybe low-metallicity stars don't like to pair up?) or if they were just broken up later.

This paper says: "It's the breakup."

• The Universal Rule: Stars everywhere start with the same "Initial Binary Distribution" (IBD).
• The Chaos Factor: The dense environment of the cluster acts like a giant pair-breaker, specifically targeting the loose couples.
• The Black Hole Hero: Black holes are the unsung heroes that manage this breakup process, preventing the cluster from collapsing and keeping the remaining "hard" couples safe.

The Bottom Line

Think of a Globular Cluster as a giant, ancient high school reunion. When the party started, everyone had a date (just like in our neighborhood). But over 12 billion years of dancing, bumping, and shoving, the couples who were holding hands loosely got separated. The ones who were locked in a tight embrace survived. And the heavy, mysterious "Black Holes" in the center of the room were the ones doing most of the shoving, keeping the party going without letting the whole building collapse.

The paper proves that if you know how many couples are left, you can figure out exactly how wild the party was when it started.

Technical Summary

1. Problem Statement

Globular clusters (GCs) exhibit significantly lower binary fractions (typically 1–10%) compared to main-sequence (MS) stars in the solar neighborhood (~50%). The physical origin of this discrepancy has been debated:

• Hypothesis A: Different star formation outcomes in the high-density, low-metallicity environments of GC progenitors.
• Hypothesis B: The dynamical processing (destruction) of primordial binaries over the cluster's lifetime.
• Hypothesis C: A combination of both.

Previous studies often relied on complex $N$-body simulations or assumed specific initial conditions that were difficult to verify. The authors aim to resolve this by testing the hypothesis of a Universal Initial Binary Distribution (IBD)—assuming GCs formed with the same binary demographics as the solar neighborhood—and determining if dynamical evolution alone can explain the observed paucity of binaries in present-day GCs.

2. Methodology

The authors employ a multi-faceted approach combining analytical formalism, energy budget calculations, and numerical simulations.

A. Analytical Formalism

• Cluster Model: Modeled as a Plummer sphere with mass $M_{cl}$ and scale length $b$.
• Binary Classification: Binaries are categorized as hard ($\eta > 1$) or soft ($\eta < 1$) based on the ratio of their orbital velocity to the cluster's velocity dispersion ($\sigma$).
  • Hard binaries: Tend to gain energy and harden (survive).
  • Soft binaries: Tend to lose energy and dissolve.
• Boundaries:
  • Hard/Soft Boundary ($a_{hs}$): Defined where the binary's binding energy equals the kinetic energy of a background star. This scales linearly with binary mass ($m_b$) and inversely with cluster density.
  • Tidal Limit ($a_{max}$): The maximum separation before the cluster's tidal field disrupts the binary.
• Initial Conditions: The authors adopt the observed binary fraction and separation distribution of the solar neighborhood (from Offner et al. 2023) as the universal IBD. They test two separation distributions: log-flat (Öpik's law) and log-normal.

B. Energy Budget Analysis

The authors calculate the energy required to dissolve the population of soft binaries ($|E_{b,s}|$) and compare it to the energy available from:

1. Primordial Hard Binaries: Burning of primordial hard binaries.
2. Stellar-Mass Black Holes (BHs): Dynamical burning of BH binaries (BBHs).
   They derive an analytical expression for the number of BHs required to supply the energy to dissolve soft binaries without causing the cluster to core-collapse prematurely.

C. Numerical Simulation

• Tool: The Cluster Monte Carlo (CMC) code.
• Setup: A simulation of a GC with $4 \times 10^5$ stars, initialized with a King profile ($W_0=5$), metallicity $Z=0.002$, and a primordial binary population matching the universal IBD (log-flat separation distribution).
• Physics: Includes direct integration of strong encounters, three-body interactions, tidal capture, and full stellar evolution (via COSMIC).
• Duration: Evolved to 13.8 Gyr.

3. Key Contributions

1. Universal IBD Validation: The paper provides a rigorous analytical argument that a universal IBD (identical to the solar neighborhood) is consistent with observed GC binary fractions, provided that soft binaries are dynamically dissolved.
2. Role of Black Holes: The study quantifies the "energy sink" problem of soft binary dissolution. It demonstrates that stellar-mass black holes are the primary energy source required to dissolve soft binaries in the early stages of cluster evolution, preventing the cluster from core-collapsing too early.
3. Birth Condition Constraints: By inverting the analytical model, the authors use observed binary fractions to constrain the birth radii of Milky Way GCs, linking present-day observations to progenitor properties.

4. Key Results

Analytical Predictions

• Soft Binary Dissolution: In dense clusters, the hard/soft boundary ($a_{hs}$) shifts to smaller separations for low-mass stars. Consequently, the vast majority of low-mass primordial binaries (which are typically wide) fall into the "soft" category and are rapidly dissolved.
• Mass Dependence: High-mass stars have a higher hard/soft boundary, meaning a larger fraction of their primordial binaries remain "hard" and survive. This predicts a rising binary fraction with stellar mass, consistent with observations of massive stars in Young Massive Clusters (YMCs).
• Quantitative Match: Using the solar neighborhood IBD, the analytical model predicts a surviving binary fraction of ~1–10% for low-mass stars in dense GCs, matching observational data perfectly.

Energy Budget & BH Burning

• Energy Cost: The energy required to dissolve the soft binary population ($|E_{b,s}|$) is significant, particularly in compact clusters.
• BH Sufficiency: The analysis shows that the energy released by the dynamical burning of a retained population of stellar-mass black holes is sufficient to dissolve the soft binaries.
• Constraint: If the primordial soft binary fraction were significantly higher than the solar neighborhood value, the required number of BHs to prevent core collapse would exceed plausible retention rates. This rules out extremely high primordial binary fractions for GC progenitors.

CMC Simulation Results

• Evolutionary Track: The simulation confirms the analytical picture.
  • Early Phase (< 3 Myr): Rapid dissolution of soft binaries causes a sharp drop in the total binary fraction (from ~30% to ~9%) and a core contraction.
  • BH Activation: Supernovae and BH formation occur. BH binaries begin to burn, injecting energy that halts core collapse and drives expansion.
  • Late Phase: The binary fraction continues to decline slowly on a relaxation timescale, powered by BH burning.
• Final State: At 13.8 Gyr, the global binary fraction settles at ~5%, and the core fraction at ~25%. These values align with both the analytical predictions and observational data for Milky Way GCs.

Application to Milky Way GCs

• Applying the model to 29 well-studied GCs, the authors infer their birth radii ($r_v$).
• The inferred birth radii span 0.5–20 pc, which is consistent with the observed sizes of Young Massive Clusters (YMCs) and Super Star Clusters (SSCs) in the local universe.
• This suggests that surviving Milky Way GCs formed with similar structural properties to modern YMCs, reinforcing the universality of the IBD.

5. Significance

• Resolving the Binary Fraction Discrepancy: The paper definitively shows that the low binary fractions in GCs are a result of dynamical processing (specifically the dissolution of soft binaries) rather than a fundamental difference in star formation physics between the solar neighborhood and GC progenitors.
• Importance of Black Holes: It reinforces the critical role of stellar-mass black holes as the "engine" regulating GC dynamics. Without BH burning, the energy sink of soft binary dissolution would likely drive most GCs into core collapse much earlier than observed.
• Modeling Implications: The authors argue that theoretical models of GC evolution must include a realistic, universal IBD (including soft binaries) to accurately reproduce the energy flow and structural evolution of clusters. Models ignoring soft binaries or assuming only hard primordial binaries may miss crucial early evolutionary phases.
• Observational Constraints: The method provides a new tool to infer the initial conditions (birth density and radius) of GCs based solely on their current binary fractions, bridging the gap between ancient GCs and modern YMCs.

In conclusion, the paper establishes that the dynamical dissolution of soft binaries, powered by black hole interactions, is the dominant mechanism shaping the binary populations of globular clusters, validating the hypothesis of a universal initial binary distribution across diverse galactic environments.