Reassessing the Spin of Second-born Black Holes in Coalescing Binary Black Holes and Its Connection to the chi_eff-q Correlation

Using updated He-star wind prescriptions in isolated binary evolution models, this study finds that wind mass loss primarily determines the spin of second-born black holes and that neither stable mass transfer nor common-envelope channels produce a correlation between mass ratio and effective inspiral spin, suggesting alternative physical mechanisms are needed to explain the observed weak anti-correlation in GWTC-4.0.

The Gist

The Big Picture: A Cosmic Dance of Black Holes

Imagine the universe as a giant dance floor. On this floor, pairs of massive stars are born, dance together, and eventually die, turning into black holes. Sometimes, these black hole pairs get so close that they crash into each other, creating ripples in space-time called gravitational waves.

Scientists have been listening to these crashes (using detectors like LIGO and Virgo) and noticed a strange pattern. They thought they saw a rule: "If the two black holes are very different in size, they must be spinning fast." It was like saying, "If you see a giant and a dwarf dancing, the giant must be doing a pirouette."

However, with new data (called GWTC-4.0), this rule seems to be breaking down. The connection between their sizes and their spins is getting weaker. This paper asks: Why is this happening? And what actually controls how fast these black holes spin?

The Main Characters: The "First-Born" and the "Second-Born"

In a binary star system, one star dies first. Let's call it the First-Born Black Hole. It usually just sits there, spinning very slowly (or not at all).

The second star is still alive. It sheds its outer layers (like a snake shedding skin) and becomes a Helium Star. This star is the "Second-Born" candidate. Eventually, it dies and becomes the second black hole.

The big question is: How fast is this Second-Born Black Hole spinning when it is born?

The Two Forces in the Dance

The authors studied two main forces that decide the spin of this second black hole:

1. The Wind (Mass Loss): Massive stars blow off gas like a giant hair dryer. This wind carries away spin.

   • The Old Idea: Scientists used a standard "Dutch Wind" recipe that said these stars blow very hard, losing a lot of spin.
   • The New Idea: The authors used a new, more accurate recipe (SV2023+). They found that at lower metal levels (like in the early universe), these stars actually blow much weaker winds. They hold onto their spin better than we thought.

2. The Tides (The Partner's Pull): Just as the Moon pulls on Earth's oceans to create tides, the first black hole pulls on the second star. This "tidal force" tries to lock the star's rotation to the orbit, making them spin together.

The Surprising Findings

The authors ran thousands of computer simulations (like a cosmic flight simulator) to see what happens. Here is what they discovered:

• The "Wind" is the Boss: The most important factor isn't how fast the star started spinning or how big its partner is. It's how much gas the star loses to the wind.
  • Analogy: Imagine a figure skater spinning. If they are wearing a heavy, fluffy coat (strong wind), the friction slows them down. If they wear a tight suit (weak wind), they keep spinning fast. The authors found that for many stars, the "coat" is much lighter than we thought.
• Size Doesn't Matter (Much): Whether the partner black hole is small or huge doesn't change the final spin much.
• Timing Doesn't Matter: It doesn't matter if the tidal forces start when the star is young or old; the wind eventually washes out those differences.
• The "Internal Mixer": There is one tricky variable: how well the star mixes its insides. If the star acts like a solid ball (mixing everything well), it spins slowly. If it acts like a layered cake (keeping layers separate), it can spin faster. This is a major uncertainty.

The "Mass Ratio" Mystery Solved?

The paper tried to solve the mystery of the "Size vs. Spin" correlation (the rule that seemed to be breaking).

• The Scenario: In some cases, the stars swap roles. The smaller star steals mass from the bigger one, becomes the bigger one, and dies first. This is called Mass Ratio Reversal.
• The Result: The authors found that in the "Stable Mass Transfer" channel (a smooth dance), 86% of the pairs undergo this reversal. In the "Common Envelope" channel (a messy, chaotic dance where they get stuck in a gas cloud), only 3% do.
• The Big Conclusion: Despite all this swapping and dancing, they found NO correlation between the size difference and the spin.

In plain English: The universe is messier than we thought. You can have a giant and a dwarf dancing, and they might both be spinning slowly, or both spinning fast. The simple rule "different sizes = fast spin" doesn't hold up when you look at the physics of how these stars actually lose their gas.

Why This Matters

This paper is like updating the instruction manual for building black holes.

1. New Wind Recipe: We now know stars lose less gas than we thought, which changes how we predict their final spins.
2. No Simple Rule: We can't just look at the sizes of two black holes and guess how fast they are spinning. We need to know the complex history of how they lost their gas.
3. Future Work: The authors admit that if we change the "internal mixer" settings or how stars get kicked at birth, the results might change. They plan to keep testing these ideas as more gravitational wave data comes in.

The Takeaway: The universe is full of complex, chaotic dances. The "rules" we thought we knew about black hole spins are being rewritten, and it turns out the wind blowing off these stars is the most important choreographer of all.

Technical Summary

1. Problem Statement

The LIGO–Virgo–KAGRA (LVK) collaboration has reported a catalog of binary black hole (BBH) mergers (GWTC-3.0 and the newer GWTC-4.0). Previous studies identified an anti-correlation between the mass ratio ($q$) and the effective inspiral spin ($\chi_{\text{eff}}$), where systems with unequal masses tended to have higher effective spins. However, recent data from GWTC-4.0 suggests this anti-correlation is weaker or potentially non-existent.

The origin of this correlation remains debated. In the context of isolated binary evolution, the spin of the second-born black hole (formed from a helium star progenitor) is a critical determinant of $\chi_{\text{eff}}$. However, the physical mechanisms governing this spin—specifically the competition between stellar winds and tidal interactions during the helium-star phase—are subject to significant uncertainties. Key issues include:

• The strength of mass loss from helium stars (He-stars).
• The efficiency of angular momentum transport within the star.
• The impact of the companion mass and the evolutionary stage of the He-star at the onset of tidal locking.

2. Methodology

The authors employed a multi-stage approach combining detailed binary evolution and rapid population synthesis:

• Detailed Binary Evolution (MESA):

  • Used the MESA code (version r15140) to model the evolution of He-star + Black Hole (BH) binaries.
  • Wind Prescription: Adopted a recently proposed, theoretically motivated wind prescription for massive He-stars (SV2023+), which combines Sander & Vink (2020) and Sander et al. (2023) formulations. This was compared against the standard "Dutch" wind scheme (Nugis & Lamers 2000).
  • Tidal Interactions: Implemented a revised formulation for dynamical tides (correcting inconsistencies in previous implementations of Zahn's theory) to calculate synchronization timescales.
  • Angular Momentum Transport: Modeled internal transport via the Tayler–Spruit (TS) dynamo, which enforces quasi-solid-body rotation.
  • Parameters: Explored a wide range of initial conditions, including He-star masses (10–70 $M_\odot$), companion BH masses (3–40 $M_\odot$), initial orbital periods (0.3–1.5 days), metallicities ($1.0, 0.1, 0.01 Z_\odot$), and initial rotation rates.
  • Collapse Assumption: Assumed direct collapse to a BH at central carbon depletion with no mass loss or natal kicks.

• Population Synthesis (COMPAS):

  • Used the COMPAS code to generate a population of BH–He star systems.
  • Applied the fitting formula derived from the MESA models to predict the spin of the second-born BH ($\chi_2$) across the population.
  • Analyzed the resulting correlation between the final mass ratio ($q$) and effective spin ($\chi_{\text{eff}}$) for two formation channels: Stable Mass Transfer (SMT) and Common Envelope (CE).

3. Key Contributions

1. Wind Prescription Comparison: Demonstrated that the SV2023+ wind prescription is substantially weaker than the standard Dutch scheme, particularly at subsolar metallicities. This significantly alters mass-loss histories and angular momentum budgets.
2. Spin Determinants: Systematically isolated the factors determining the second-born BH spin, finding that wind mass loss is the dominant regulator, while the companion mass and the evolutionary stage of the He-star at the onset of tides have negligible effects.
3. Fitting Formula: Derived a new fitting formula for the spin magnitude of the second-born BH ($\chi_2$) as a function of initial He-star mass, orbital period, and metallicity, facilitating its use in large-scale population studies.
4. Re-evaluation of Correlations: Provided a rigorous test of the $q$–$\chi_{\text{eff}}$ anti-correlation using these updated physical prescriptions.

4. Key Results

A. Physics of the Second-born BH Spin

• Dominance of Winds: At solar metallicity, strong stellar winds efficiently remove angular momentum, resulting in low-spin BHs regardless of initial rotation. At low metallicity ($Z \le 0.1 Z_\odot$), winds are weak, and tidal interactions become dominant.
• Insensitivity to Companion Mass: The mass of the first-born BH (companion) has negligible impact on the final spin of the second-born BH, provided the system is tidally synchronized.
• Insensitivity to Evolutionary Stage: Whether tidal interactions begin when the He-star is at the Zero-Age Helium Main Sequence (ZAMS) or has already burned 30% of its core helium does not significantly alter the final spin.
• Role of Initial Rotation: Initial rotation only significantly affects the final spin at low metallicities where winds are weak and tidal coupling is inefficient. At high metallicities, winds wash out the memory of initial rotation.
• Angular Momentum Transport: The inclusion of the TS dynamo (efficient transport) results in significantly lower BH spins compared to models without it, as it prevents the core from spinning up independently of the envelope.

B. The $q$ – $\chi_{\text{eff}}$ Correlation

• Mass Ratio Reversal:
  • In the Stable Mass Transfer (SMT) channel, 85.8% of BBHs undergo mass-ratio reversal (the accretor becomes the more massive BH).
  • In the Common Envelope (CE) channel, only 2.8% exhibit mass-ratio reversal.
• Absence of Correlation: Crucially, the study found no correlation between the mass ratio $q$ and the effective spin $\chi_{\text{eff}}$ in either the SMT or CE channels when using the default fiducial model and the new wind prescriptions.
  • This contrasts with earlier findings that suggested a strong anti-correlation, suggesting that previous correlations may have been driven by different wind assumptions or population synthesis parameters.

5. Significance and Implications

• Resolving the GWTC-4.0 Discrepancy: The finding of no correlation between $q$ and $\chi_{\text{eff}}$ aligns with the weaker anti-correlation observed in the latest GWTC-4.0 data, suggesting that the "standard" isolated binary evolution models with updated wind physics may not naturally produce the strong anti-correlation seen in earlier catalogs.
• Refining Formation Models: The work highlights that the efficiency of angular momentum transport (TS dynamo) and the strength of He-star winds are the most critical uncertainties in predicting BBH spins.
• Upper Limits: The derived spins are treated as upper limits because the models assume direct collapse without mass loss or natal kicks. Future work must account for accretion feedback and potential natal kicks to fully explain the observed spin distribution, including events with negative $\chi_{\text{eff}}$.
• Future Directions: The authors plan to investigate how alternative physical prescriptions (e.g., different CE efficiencies, natal kick distributions, or accretion feedback) influence the $q$–$\chi_{\text{eff}}$ relationship to better match the growing LVK dataset.

In summary, this paper provides a critical reassessment of BBH spin formation physics, demonstrating that updated wind prescriptions and angular momentum transport mechanisms decouple the mass ratio from the effective spin in isolated binary evolution, challenging previous interpretations of the GWTC catalogs.