Diagnosing the Properties and Evolutionary Fates of Black Hole and Wolf-Rayet X-ray Binaries as Potential Gravitational Wave Sources for the LIGO-Virgo-KAGRA Network

This study utilizes revised accretion efficiency models and MESA binary evolution calculations to constrain the masses and spins of black holes in IC 10 X-1, NGC 300 X-1, and Cyg X-3, revealing that these Wolf-Rayet X-ray binaries are likely to evolve into binary black hole systems that will merge within a Hubble time and serve as detectable gravitational wave sources for the LVK network.

The Gist

Imagine the universe as a giant, chaotic dance floor. Most stars dance alone, but some pair up in tight duets. This paper is about three specific, very intense couples dancing on this cosmic floor: IC 10 X-1, NGC 300 X-1, and Cyg X-3.

Here is the simple story of what the scientists found, using everyday analogies.

The Characters: The "Heavyweight" and the "Wind-Blower"

In these three couples, one partner is a Black Hole (a super-dense, invisible object with gravity so strong it eats everything nearby). The other partner is a Wolf-Rayet star (a massive, aging star that is essentially a cosmic hair dryer, blasting out incredibly strong winds of gas).

• The Dance: The Black Hole is trying to "eat" the gas blowing off the Wolf-Rayet star.
• The Goal: The scientists wanted to figure out exactly how heavy these Black Holes are, how fast they are spinning, and what will happen to them in the future. Specifically, they wanted to know: Will these couples eventually crash into each other and create a "gravitational wave" (a ripple in space-time) that our detectors (LIGO/Virgo) can hear?

The Problem: The "Menu" Was Too Big

Before this study, scientists had a very loose idea of how heavy these Black Holes were. It was like trying to guess the weight of a person by looking at them from a mile away; the estimates ranged from "a bit heavy" to "absolutely massive."

The researchers used a super-computer simulation (called MESA) to act like a cosmic time machine. They ran thousands of simulations, tweaking the starting weights and speeds of these stars to see which ones could produce the exact systems we see today.

The Big Discoveries

1. The "Heavy" Black Holes Aren't as Heavy as We Thought

For the first two couples (IC 10 X-1 and NGC 300 X-1), the scientists realized the Black Holes are actually lighter than previously thought.

• The Analogy: Imagine you thought your neighbor was a 300-pound bodybuilder. After measuring their footprint and how much they eat, you realize they are actually a 150-pound athlete.
• The Result: The Black Hole in IC 10 X-1 is likely under 25 times the mass of our Sun (not 30+), and the one in NGC 300 X-1 is under 15 times the mass of our Sun. This is a big deal because it changes how we understand the "family tree" of these stars.

2. The "Spin" of the Cosmic Top

Black Holes spin like tops. The faster they spin, the more energy they have.

• The Discovery: For Cyg X-3 (the third couple), the scientists found that the Black Hole's spin is limited. It's spinning, but not as wildly as some thought. It's capped at about 60% of the maximum speed possible.
• Why it matters: If it were spinning faster, the physics of how it eats gas wouldn't match what we see in the X-ray light.

3. The "Lower-Mass Gap" Mystery

There is a weird gap in the universe where we rarely see Black Holes. They are usually either very light (like a neutron star) or very heavy. There is a "missing middle" (between 2 and 5 times the Sun's mass).

• The Surprise: The Wolf-Rayet star in Cyg X-3 is likely to die and turn into a Black Hole that falls right into this "missing middle" gap. It's like finding a rare coin in a pile of pennies and quarters.

The Future: The Grand Finale

The most exciting part of the paper is the prediction of the future.

• The Crash: Because these couples are dancing so close together, they are slowly losing energy. Eventually, they will spiral inward and smash into each other.
• The Ripple: When they smash, they will create a massive "tsunami" in space-time called a Gravitational Wave.
• The Verdict:
  • IC 10 X-1 and Cyg X-3 are almost guaranteed to crash and create a detectable wave within the next 10 billion years (a "Hubble time").
  • NGC 300 X-1 will likely crash too, unless its Black Hole is the lightest possible version (9 times the Sun's mass). If it's that light, they might drift apart forever and never meet.

The "Secret Sauce": A Better Recipe

The scientists improved their calculations by using a new, more accurate recipe for how the Black Hole "eats" the wind from the star.

• The Analogy: Imagine you are trying to catch rain in a bucket while running. The old recipe assumed you catch a lot of water. The new recipe realizes that if you run too fast or the wind is too strong, you actually catch much less.
• The Result: This small change in the "eating" recipe helped them rule out the super-heavy Black Hole theories and gave them the tighter, more accurate numbers mentioned above.

Summary

This paper is like a detective story where the scientists used a cosmic time machine to figure out the true identities of three mysterious couples. They found that the "villains" (the Black Holes) are slightly lighter than we thought, one of them is likely to create a rare "middle-weight" monster, and all of them are on a collision course that will eventually shake the fabric of the universe itself.

Technical Summary

1. Problem Statement

Wolf-Rayet (WR) X-ray binaries, consisting of a stellar-mass black hole (BH) accreting from a massive WR star, are considered immediate progenitors of double compact objects (specifically Binary Black Holes, BBHs). Three specific systems—IC 10 X-1, NGC 300 X-1, and Cyg X-3—are of particular interest due to their short orbital periods ($<1.5$ days), suggesting they will merge within a Hubble time and become detectable by the LIGO-Virgo-KAGRA (LVK) network.

However, previous estimates of their properties (such as BH masses and spins) and their evolutionary fates have been uncertain or based on population synthesis models (e.g., StarTrack, SEVN) with varying assumptions. There is a need for detailed, individual binary evolution calculations to:

• Constrain the current physical parameters (masses, spins, orbital periods) of these specific systems.
• Determine their ultimate evolutionary outcomes (merger times, final masses).
• Assess their detectability as gravitational wave (GW) sources.

2. Methodology

The authors utilized the Modules for Experiments in Stellar Astrophysics (MESA) code (release mesa-r15140) to perform detailed binary evolution calculations. Key methodological components include:

• Physics Implementation:

  • Convection & Mixing: Mixing-length theory ($\alpha_{mlt}=1.93$), Ledoux criterion for convective boundaries, step-overshooting, semiconvection, and thermohaline mixing.
  • Rotation & Angular Momentum: Modeled as diffusive processes accounting for Goldreich–Schubert–Fricke instability, Eddington–Sweet circulations, and shear mixing.
  • Mass Loss: Adopted updated WR wind prescriptions (Higgins et al. 2021) with rotationally enhanced mass loss.
  • Tides: Applied dynamical tides theory for helium-rich stars with radiative envelopes.
  • Supernova & Remnant Formation: Used the "delayed" supernova prescription (Fryer et al. 2012) to calculate remnant masses, assuming direct collapse for BH formation (no natal kicks).

• Accretion Modeling (Key Innovation):

  • The study adopted a revised Bondi-Hoyle-Lyttleton (BHL) accretion efficiency (Tejeda & Toalá 2025) to address nonphysical efficiencies when wind velocity is comparable to orbital velocity.
  • Wind velocity was modeled using a $\beta$-velocity law.
  • The formation of accretion disks was evaluated based on specific angular momentum; for disk-accreting systems, X-ray luminosity ($L_X$) was calculated using standard disk accretion prescriptions.

• Constraints & Grid Search:

  • The authors ran parameter grids tailored to each system, varying initial WR masses, orbital periods, and BH spins.
  • Models were filtered to match observed properties: WR mass, orbital period, and X-ray luminosity ($L_X$).
  • A chemical constraint was applied: C/He abundance ratio $< 0.01$ to ensure evolutionary consistency.

3. Key Contributions

• Refined Mass Constraints: Provided significantly tighter upper limits on BH masses for IC 10 X-1 and NGC 300 X-1 compared to previous literature.
• Spin Constraints: Established an upper limit on the BH spin for Cyg X-3 based on X-ray luminosity constraints.
• Evolutionary Fates: Predicted the specific merger timescales and final masses of the resulting compact binaries, identifying Cyg X-3 as a potential source of a "lower-mass-gap" black hole.
• GW Source Characterization: Calculated the chirp mass ($M_{chirp}$) and effective inspiral spin ($\chi_{eff}$) for the future BBH systems, providing specific targets for LVK detection.

4. Key Results

IC 10 X-1

• BH Mass: Constrained to $M_{BH} \lesssim 25 M_\odot$ (lower than previous estimates of $\sim 30 M_\odot$).
• WR Mass: Current mass $\sim 17–35 M_\odot$; initial mass $18–39 M_\odot$.
• Accretion: Accreted mass is low ($\sim 10^{-4} - 10^{-2} M_\odot$), insufficient to significantly spin up the BH.
• GW Fate: Will form a BBH merging within a Hubble time.
  • $M_{chirp} \approx 10.7 - 16.3 M_\odot$.
  • $\chi_{eff} \approx 0.3 - 0.6$.
  • Merger time: $\log(T_{merger}/\text{Myr}) \in [3.20, 4.06]$.

NGC 300 X-1

• BH Mass: Constrained to $M_{BH} \lesssim 15 M_\odot$ (significantly lower than previous estimates up to $28 M_\odot$).
• WR Mass: Current mass $\sim 23.95 - 25.0 M_\odot$.
• GW Fate: Generally expected to merge within a Hubble time.
  • Exception: If the BH mass is at the lower limit of previous studies ($9 M_\odot$), the system will not merge within a Hubble time.
  • $M_{chirp} \approx 10.4 - 13.0 M_\odot$.
  • $\chi_{eff} \approx 0.3 - 0.5$.
  • Merger time: $\log(T_{merger}/\text{Myr}) \in [3.48, 4.26]$.

Cyg X-3

• BH Mass: Well-constrained at $\approx 7.2 M_\odot$.
• BH Spin: Constrained to $\chi_1 \lesssim 0.6$ based on $L_X$ matching.
• WR Mass: Current mass $\approx 10.6 - 12.3 M_\odot$.
• GW Fate:
  • The WR star is predicted to collapse into a lower-mass-gap BH ($3.9 - 4.6 M_\odot$).
  • The resulting system may be a BBH or a BH-NS system (if the remnant is a neutron star, though the paper leans toward a low-mass BH).
  • The second-born BH is expected to be rapidly spinning ($\chi_2 \approx 0.55 - 0.75$) due to strong tidal effects.
  • $M_{chirp} \approx 4.6 - 5.0 M_\odot$.
  • $\chi_{eff} \approx 0.3 - 0.6$.
  • Merger time: $\log(T_{merger}/\text{Myr}) \in [2.17, 2.31]$ (significantly shorter than the other two).

5. Significance

• LVK Detection Prospects: The study confirms that these three systems are viable progenitors for GW events. Cyg X-3 is highlighted as a unique candidate for producing a merger involving a lower-mass-gap object, which is a critical area of study for understanding the BH mass distribution.
• Model Validation: The use of a revised accretion efficiency showed negligible impact on the specific systems studied, validating the robustness of standard BHL models for these parameters, while the detailed MESA modeling provided more precise constraints than population synthesis alone.
• Observational Guidance: The derived constraints on BH masses and spins provide specific targets for future X-ray spectroscopy and GW data analysis, helping to distinguish between different formation channels of compact binaries.
• Evolutionary Insight: The results reinforce that wind-fed accretion in these systems does not significantly spin up the BHs, implying that high spins observed in other HMXBs (like Cygnus X-1) likely require different accretion histories (e.g., conservative mass transfer or Roche-lobe overflow).