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

Probleemstelling

Deze studie richt zich op een unieke klasse van röntgendubbelsterren: systemen bestaande uit een sterrenmassa-zwart gat (BH) dat materiaal accreteert van een Wolf-Rayet (WR)-ster. Drie specifieke systemen worden onderzocht: IC 10 X-1, NGC 300 X-1 en Cyg X-3.
Deze systemen zijn van groot belang omdat ze met hun korte omlooptijden (< 1,5 dagen) directe voorlopers zijn van dubbele zwarte gaten (BBH) of zwarte gat-neutronenster (BHNS) systemen. Deze eindstadia zijn potentiële bronnen van zwaartekrachtsgolven (GW) die detecteerbaar zijn door het LIGO-Virgo-KAGRA (LVK) netwerk. Eerdere schattingen van de massa's en eigenschappen van deze systemen waren echter onzeker of gebaseerd op minder geavanceerde populatiesynthesemodellen. Het doel van dit onderzoek is om de eigenschappen van deze systemen nauwkeuriger te construeren en hun evolutionaire lot te voorspellen als GW-bronnen.

Methodologie

De auteurs hebben uitgebreide berekeningen van dubbelster-evolutie uitgevoerd met behulp van de MESA-code (Modules for Experiments in Stellar Astrophysics, versie r15140). De belangrijkste methodologische aspecten zijn:

• Fysische modellen: Er zijn modellen gebruikt voor convectie (mixing-length theory), rotatie-menging, en windverlies. Voor het massaverlies van WR-sterren is een recentere voorspelling gehanteerd die rekening houdt met rotatie-geïnduceerd massaverlies.
• Accretie-efficiëntie: Een cruciale update in deze studie is de toepassing van een herziene accretie-efficiëntie binnen het standaard Bondi-Hoyle-Lyttleton (BHL) raamwerk. De auteurs gebruiken een uitdrukking van Tejeda & Toalá (2025) die rekening houdt met situaties waar de windsnelheid vergelijkbaar is met of lager is dan de omloopsnelheid, wat voorkomt dat onfysische accretie-efficiënties ontstaan.
• Modelgrid: Voor elk van de drie systemen werden specifieke grids van modellen gelopen, variërend in initiële massa's van de WR-ster, initiële omlooptijden en BH-massa's.
• Beperkingen: De modellen werden gefilterd op basis van waargenomen eigenschappen: de massa van de WR-ster, de omlooptijd, de röntgenluminositeit ($L_X$) en de rochelobevullingsfactor ($f_{RL}$). Voor Cyg X-3 werd ook de BH-spin gebruikt als beperking.
• Toekomstvoorspelling: Na het construeren van de huidige eigenschappen, werd de evolutie verder doorgerekend tot de vorming van een tweede zwart gat en de daaropvolgende samensmeltingstijd ($T_{merger}$) berekend met de formule van Peters (1964).

Belangrijkste Bijdragen en Resultaten

1. IC 10 X-1:

• BH-massa: De bovenste limiet voor de BH-massa is geconstrueerd tot $\lesssim 25 M_\odot$, wat significant lager is dan eerdere schattingen (die tot 30 $M_\odot$ gingen). Dit verschil wordt toegeschreven aan de huidige beperkingen op de röntgenluminositeit.
• Accretie: De geaccreteerde massa is klein ($\sim 2 \times 10^{-4}$ tot $1.7 \times 10^{-2} M_\odot$), wat betekent dat de spin van het zwarte gat niet significant is toegenomen door accretie.
• Evolutie: Het systeem zal evolueren naar een BBH dat binnen de Hubble-tijd zal samensmelten. De chirp-massa ($M_{chirp}$) ligt tussen 10,7 en 16,3 $M_\odot$.

2. NGC 300 X-1:

• BH-massa: De bovenste limiet voor de BH-massa is geconstrueerd tot $\lesssim 15 M_\odot$. Dit is aanzienlijk lager dan eerdere schattingen (tot 28 $M_\odot$).
• Uitzondering: Als de BH-massa aan de ondergrens ligt (9 $M_\odot$, zoals gesuggereerd door eerdere studies met een continuüm-fit methode), zal het systeem niet binnen de Hubble-tijd samensmelten.
• Evolutie: Voor de meeste geaccepteerde massa's zal het systeem een BBH vormen met een chirp-massa van 10,4 tot 13,0 $M_\odot$ die binnen de Hubble-tijd samensmelt.

3. Cyg X-3:

• Spin-beperking: De spin van het huidige zwarte gat is geconstrueerd tot $\chi_1 \lesssim 0.6$.
• WR-ster: De WR-ster heeft een massa van ongeveer 10,6 - 12,3 $M_\odot$.
• Evolutionair Lot: De WR-ster zal waarschijnlijk evolueren naar een zwart gat in de zogenaamde "lower-mass gap" (tussen 3,9 en 4,6 $M_\odot$).
• GW-signatuur: Het resulterende BBH-systeem heeft een zeer korte samensmeltingstijd (log($T_{merger}$) $\approx$ 2,17 - 2,31 Myr) en een lage chirp-massa (4,6 - 5,0 $M_\odot$). Dit maakt het een interessante kandidaat voor detectie, mogelijk zelfs als een BHNS-systeem als de tweede component als neutronenster wordt geïnterpreteerd (hoewel de studie uitgaat van een BH-vorming).

Significantie

Deze studie levert een belangrijke bijdrage aan het begrip van de evolutie van hoge-massa dubbelsterren en de populatie van toekomstige zwaartekrachtgolfbronnen:

1. Nauwkeurigere Massa's: Door het toepassen van een herziene accretie-efficiëntie en het strikt construeren op waargenomen röntgenluminositeit, worden de massa's van de zwarte gaten in IC 10 X-1 en NGC 300 X-1 aanzienlijk lager geschat dan eerder gedacht. Dit heeft gevolgen voor de voorspellingen van de massa-verdelingen van BBH-systemen die door LVK worden gedetecteerd.
2. Spin-beperkingen: De studie suggereert dat de spins van deze zwarte gaten niet extreem hoog zijn (zoals soms wordt aangenomen voor wind-gevoede systemen), maar beperkt blijven tot $\chi \lesssim 0.6$.
3. Voorspelling van GW-bronnen: Het bevestigt dat deze drie systemen directe voorlopers zijn van samensmeltingen binnen de leeftijd van het heelal. Cyg X-3 is bijzonder omdat het een potentieel bron is voor een samensmelting in de "mass gap", een gebied dat van groot belang is voor het begrijpen van supernova-mechanismen en de vorming van compacte objecten.
4. Methodologische Vooruitgang: Het gebruik van gedetailleerde evolutiecodes (MESA) in plaats van snellere populatiesynthesecodes (zoals StarTrack of SEVN) biedt een dieper inzicht in de specifieke fysische processen die deze unieke systemen bepalen.

Kortom, dit werk verfijnt de parameters van drie sleutelsystemen en biedt een robuuste voorspelling voor hun rol als bronnen voor de volgende generatie zwaartekrachtgolfobservatoria.