A Stellar Role Reversal: Multiple Features in the Mass and Mass Ratio Distributions of Merging Binary Black Holes from Stable Mass Transfer

Technische Zusammenfassung: Eine stellare Rollenumkehr: Mehrfache Merkmale in den Massen- und Massenverhältnisverteilungen verschmelzender binärer Schwarzer Löcher aus stabilem Massentransfer

Problemstellung
Die Ursprünge der Populationen verschmelzender binärer Schwarzer Löcher (BBH) bleiben trotz des wachsenden Katalogs an Gravitationswellen-Ereignissen ungewiss. Während Beobachtungen Merkmale in der BBH-Massenfunktion offenbart haben (z. B. Peaks nahe $10 M_\odot$ und $35 M_\odot$) sowie in den Massenverhältnisverteilungen, ist es schwierig, die Beiträge verschiedener Entstehungsmodelle zu entwirren, da dies aufgrund von Degenerationen in messbaren Eigenschaften und Unsicherheiten in der Sternentwicklung und den binären Interaktionen eine Herausforderung darstellt. Diese Arbeit untersucht spezifisch den Kanal des „stabilen Massentransfers“ (nur-SMT-Kanal), bei dem alle Phasen des Roche-Lobe-Überflusses (RLOF) stabil bleiben, unter Ausschluss des Common-Envelope-Kanals (CE). Die Autoren zielen darauf ab, zu bestimmen, wie Annahmen über die Akkretionseffizienz des Massentransfers ($\beta$) und Stabilitätskriterien die primären Massen- und Massenverhältnisverteilungen von verschmelzenden BBHs prägen.

Methodik
Die Studie verwendet einen dualen Ansatz, der die binäre Populationssynthese (BPS) mit einem analytischen Rahmen kombiniert:

1. Populationssynthese: Die Autoren nutzen den cosmic-Code (v4.1.0), um eine kosmologische Population verschmelzender BBHs zu simulieren.

   • Anfangsbedingungen: Die Simulationen decken ein Gitter von 30 Metallizitäten ($10^{-4}$ bis $0,03$) mit einer initialen Massenfunktion (Kroupa 2001) im Bereich von $0,08$ bis $100 M_\odot$ und einer Binärfraktion von 0,7 ab.
   • Physik: Das Modell beinhaltet metallizitätsabhängige Winde (Vink-Preskriptionen), tidegebundene Zirkularisierung und spezifische Stabilitätskriterien für den Massentransfer basierend auf dem Entwicklungszustand des Donors (Fall A, B, C und BC).
   • Bildung kompakter Objekte: Die Massen der Schwarzen Löcher (BH) werden mittels der „delayed“ Supernova-Preskription (Fryer et al. 2012) mit einer minimalen BH-Masse von $3 M_\odot$ zugewiesen. Paarinstabilitäts-Supernovae (PISN) werden modelliert, um die maximale BH-Masse bei $44 M_\odot$ zu begrenzen. Natal Kicks werden aus einer Lognormalverteilung gezogen und durch Fallback bei massereichen BHs unterdrückt.
   • Massentransfer: Das Referenzmodell nimmt einen konservativen Massentransfer ($\beta=1$) für akkretierende Sterne und eine Eddington-limitierte Akkretion ($\beta \approx 0$) für kompakte Objekt-Akkretoren an.
   • Kosmologie: Eine Sternentstehungsrate (SFH) basierend auf der IllustrisTNG-Simulation wird via cosmerge angewendet, um einen Katalog von $10^6$ Verschmelzungen zu generieren, der auf SMT-nur-Systeme gefiltert ist, die innerhalb einer Hubble-Zeit verschmelzen.

2. Analytischer Rahmen: Die Autoren passen ein analytisches Modell (van Son et al. 2022) an, um die Akkretionseffizienz ($\beta_1$), kritische Massenverhältnisse ($q_{crit}$) und initiale Massenverhältnisse ($q_{ZAMS}$) mit den resultierenden minimalen BH-Massen zu verknüpfen. Modifikationen umfassen einen variablen Kernmasseanteil ($f_{core}$) und spezifische Supernova-Massenverlustfunktionen, um die cosmic-Simulationen besser abzubilden.

Wesentliche Beiträge

• Identifizierung von Subpopulationen: Die Arbeit zeigt, dass der SMT-nur-Kanal natürlicherweise zwei beobachtbar unterscheidbare Subpopulationen erzeugt: eine Massenverhältnis-umgekehrte (MRR) Population und eine Nicht-MRR Population.
• Mechanismus der Trennung: Die Autoren ordnen die Grenzen dieser Subpopulationen direkt den Stabilitätskriterien des Massentransfers in verschiedenen Entwicklungsstadien (Fall A, B, C) zu. Sie zeigen, dass die Stabilitätskriterien für verschiedene Donortypen komplementäre Regionen in der Primärmassen–Massenverhältnis-Ebene ausschneiden.
• Parametersensitivität: Die Studie untersucht systematisch den Einfluss der Akkretionseffizienz ($\beta_1$) auf die Verschmelzungsrate und die Massenverteilung und bietet einen Rahmen, um zu verstehen, wie die Variation dieses unsicheren Parameters die resultierende BBH-Population verändert.

Ergebnisse

1. Distinkte Merkmale in Massenverteilungen:

   • MRR-Subpopulation: Charakterisiert durch hohe Primärmassen ($20-40 M_\odot$) und hohe Massenverhältnisse ($q_{BBH} > 0,75$, mit einem Peak nahe $0,9$). Diese Systeme durchlaufen vorwiegend Fall A Massentransfer während des ersten Massentransfer-Ereignisses (MT1).
   • Nicht-MRR-Subpopulation: Charakterisiert durch niedrigere Primärmassen (Peak nahe $10 M_\odot$) und eine breitere Massenverhältnisverteilung ($q_{BBH} \sim 0,3 - 0,7$). Diese Systeme entstehen primär durch Fall B oder C Massentransfer während MT1.
   • Trennung: Ein Primärmassen-Schnitt bei $\sim 20 M_\odot$ trennt diese beiden Subpopulationen im Referenzmodell effektiv.

2. Rolle der Akkretionseffizienz ($\beta_1$):

   • Eine Erhöhung von $\beta_1$ (konservativerer Massentransfer) reduziert die Gesamtverschmelzungsrate, indem sie den Parameterraum für stabile Entwicklung einschränkt.
   • Ein höheres $\beta_1$ erhöht den Anteil an MRR-Systemen und verschiebt die Primärmassenverteilung beider Subpopulationen zu höheren Massen.
   • Die MRR-Subpopulation wird prominenter und erstreckt sich zu höheren Massen, wenn $\beta_1$ steigt, während der Nicht-MRR-Peak nahe $10 M_\odot$ bleibt, aber schmaler wird.

3. Stabilitätskriterien als Grenzen:

   • Der Übergang zwischen MRR- und Nicht-MRR-Populationen wird durch das Zusammenspiel des initialen Massenverhältnisses ($q_{ZAMS}$) und der Stabilitätsgrenzen ($q_{crit}$) für spezifische Massentransferfälle gesteuert.
   • Fall A Massentransfer hat die strengsten Stabilitätskriterien, was hohe $q_{ZAMS}$ und hohe Initialmassen erfordert, um stabil zu bleiben, und somit den hochmassigen MRR-Peak besetzt.
   • Fall C Massentransfer erlaubt niedrigere $q_{ZAMS}$ und niedrigere Massen, was zum niedrigmassigen Nicht-MRR-Peak beiträgt.

Bedeutung und Behauptungen
Die Autoren behaupten, dass die Physik des stabilen Massentransfers distinkte, auflösbare Merkmale in Gravitationswellen-Populationen erzeugt. Speziell sagt der SMT-nur-Kanal eine Bimodalität in der Primärmassenverteilung und eine Dichotomie in den Massenverhältnissen voraus, die mit den aktuellen und zukünftigen Fähigkeiten von GW-Detektoren übereinstimmt. Die Studie legt nahe, dass die beobachteten Merkmale im GWTC-5 Katalog (z. B. der Niedrigmassen-Peak und die hochmassigen Hoch-Verhältnis-Systeme) Signaturen dieser distinkten SMT-Subpopulationen sein könnten.

Die Arbeit schließt bescheiden, dass aktuelle Modelle aufgrund von Degenerationen keine definitive Zuordnung beobachteter Ereignisse zu spezifischen Kanälen vornehmen können, aber die identifizierten Merkmale (der Split bei $\sim 20 M_\odot$ und die hoch-$q$ MRR-Population) als testbare Vorhersagen dienen. Zukünftige Beobachtungen der BBH-Massen- und Massenverhältnis-Ebene könnten helfen, die Physik des SMT, einschließlich der Akkretionseffizienz und der Stabilitätskriterien, einzugrenzen und somit unser Verständnis der Entwicklung massereicher binärer Sterne zu verbessern. Die Autoren räumen Einschränkungen hinsichtlich der Sternradius-Preskriptionen, des Wind-Massenverlusts und der Vereinfachung fester Stabilitätskriterien ein und merken an, dass diese Faktoren die exakten Positionen der vorhergesagten Merkmale verschieben könnten.