Accretion of Primordial Black Holes in Stellar Interiors

Technische Samenvatting: Accretie van Primordiale Zwarte Gaten in Stellaire Interieurs

Probleemstelling
Primordiale Zwarte Gaten (PBH's) blijven een levensvatbare kandidaat voor donkere materie, met name in het slecht geconstreinde sub-maan massa-bereik ($10^{-16} M_\odot \lesssim M_{BH} \lesssim 10^{-10} M_\odot$). Sterren bieden een uniek laboratorium om PBH's in dit regime te onderzoeken, aangezien een gevangen PBH ingebed in een stellaire kern de omringende gas kan accreteren, wat de evolutie van de ster potentieel kan veranderen. Eerdere studies (bijv. Bellinger et al. 2023) hebben dit scenario gemodelleerd met behulp van 1D stellaire evolutiecodes, maar vertrouwden op aannames afgeleid van dunne accretieschijftheorie. Specifiek namen zij een vaste radiatieve efficiëntie $\eta \approx 0,08$ aan en pasten zij de Eddington-luminositeitslimiet toe over alle massa's. Echter, voor PBH's in dit massabereik is de Bondi-straal ($r_B$) ordes van grootte kleiner dan de stellaire circularisatieradius, wat de vorming van een schijf voorkomt. In plaats daarvan is de accretie sferisch, waarbij het gas wordt gecomprimeerd en verhit tot relativistische temperaturen ($T \sim 10^{10}–10^{12}$ K) nabij de Schwarzschild-straal ($r_S$). In dit regime verschillen de microfysica (relativistische bremsstrahlung, paarprocessen, Klein–Nishina verstrooiing) en de optische diepte van de flow fundamenteel van het dunne-schijfparadigma, wat een zelfconsistente berekening van de accretiesnelheid ($\dot{M}$), de radiatieve efficiëntie ($\eta$) en de flowstructuur noodzakelijk maakt.

Methodologie
De auteurs lossen de tijdafhankelijke sferische Euler-vergelijkingen op met een impliciete koelingsterm om de accretieflow naar een PBH ingebed in een zonachtige stellaire kern te modelleren.

• Microfysica: De koelingsterm $\epsilon(\rho, T)$ bevat relativistische bremsstrahlung (elektron-ion en elektron-elektron) en elektron-positron annihilatie. De opaciteit houdt rekening met Thomson-verstrooiing, free-free absorptie en Klein–Nishina onderdrukking bij hoge fotonenergieën.
• Numerieke Aanpak: De auteurs maken gebruik van radbondi, een open-source Python-pakket, om de flow te evolueren van een initiële adiabatische Bondi-profiel naar een zelfconsistente stationaire toestand. Deze tijdafhankelijke aanpak is noodzakelijk omdat, in het koelingsgedomineerde regime, het geluidspunt een focus (complexe eigenwaarden) wordt in plaats van een zadelpunt, waardoor standaard stationaire ODE-schietmethoden ongeldig worden.
• Regime Afhandeling:
  • Optisch Dun ($M_{BH} \lesssim 5 \times 10^{-13} M_\odot$): De Euler-solver berekent direct de fysieke radiatieve efficiëntie en de accretiesnelheid.
  • Optisch Dik ($M_{BH} \gtrsim 5 \times 10^{-13} M_\odot$): De Bondi-sfeer wordt optisch dik en komt in een foton-trapping regime terecht. De auteurs hanteren het analytische voorschrift van Begelman (1979), waarbij lokaal geproduceerde fotonen binnen de trapping-straal ($r_{trap}$) terug naar het zwarte gat worden geadvecteerd, terwijl de ontsnappende luminositeit zelfregulerend is bij $L_{esc} \lesssim 0,6 L_{Edd}$.
• Feedback: De studie modelleert zelfconsistent de radiatieve feedback, waarbij de luminositeit die buiten de Bondi-sfeer wordt gedeponeerd, het omgevingplasma opwarmt. Deze verhitting wordt gemodereerd door convectief transport (Mixing Length Theory), wat extreme temperatuurverhogingen voorkomt die door pure radiatieve diffusie zouden worden voorspeld.

Belangrijkste Resultaten
De studie identificeert drie onderscheidende sferische accretieregimes en levert de volgende resultaten op:

1. Drie Accretieregimes:

   • Hot Bondi ($M_{BH} \lesssim 10^{-14} M_\odot$): Koeling is dynamisch verwaarloosbaar. De flow is bijna adiabatisch en bereikt temperaturen van $T \sim 10^{11}$ K nabij $r_S$. De binnenste flow wordt collisionless voordat deze de paarproductie-temperaturen bereikt, wat de paar-evenwicht onderdrukt. Bremsstrahlung domineert de koeling.
   • Bremsstrahlung Koeling ($10^{-14} \lesssim M_{BH} \lesssim 5 \times 10^{-13} M_\odot$): Koeling wordt efficiënt en voert thermische energie af voordat het gas de paarproductie-temperaturen bereikt. De flow nadert een isothermische toestand ($T \approx T_\infty$), en de accretiesnelheid wordt met een factor $\sim 2–7$ verhoogd ten opzichte van de adiabatische Bondi-rate.
   • Foton Trapping ($M_{BH} \gtrsim 5 \times 10^{-13} M_\odot$): De Bondi-sfeer is optisch dik. Hoewel de lokaal uitgezonden luminositeit snel stijgt met de massa, wordt de ontsnappende luminositeit begrensd door de Begelman-limiet ($L_{esc} \lesssim 0,6 L_{Edd}$). Cruciaal is dat de accretiesnelheid dicht bij de Bondi-waarde blijft ($\dot{M} \propto M_{BH}^2$) omdat er geen sferische Eddington-limiet is voor de instroom; stralingsdruk moduleert enkel de effectieve zwaartekracht.

2. Radiatieve Efficiëntie: De zelfconsistent berekende radiatieve efficiëntie is $\eta \approx 10^{-2}$ over het optisch dunne bereik, een orde van grootte lager dan de canonieke dunne-schijfwaarde ($\eta \approx 0,08$) die in eerdere literatuur werd aangenomen. De efficiëntie piekt nabij $10^{-15} M_\odot$ maar daalt in de koelings- en isotherme limieten.

3. Groeitrajecten: Omdat de sferische Eddington-limiet niet van toepassing is, blijft de groei van PBH's super-exponentieel ($\dot{M} \propto M_{BH}^2$) gedurende het foton-trapping regime. De groeibottleneck wordt volledig bepaal door de trage Hot Bondi-fase bij de laagste massa's.

4. Kritische Initiële Massa: De kritische initiële massa vereist voor een PBH om een zonachtige ster binnen een Hubble-tijd te consumeren is $M_{0,crit} \sim 10^{-16} M_\odot$. Dit is vijf ordes van grootte lager dan eerdere schattingen ($\sim 10^{-11} M_\odot$) die een hogere efficiëntie en een Eddington-bottleneck aannamen.

5. Magnetische Velden: De auteurs stellen vast dat de accumulatie van magnetische flux verwaarloosbaar is. Het omgevingsveld in de zonnekern is te zwak om de flow dynamisch te stoppen, en magnetische reconnectie beperkt elke fluxopbouw tot niveaus die ver onder de dynamische belangrijkheid liggen ($P_{mag}/P_{ram} \ll 1$).

Betekenis en Claims
Het artikel claimt de eerste zelfconsistente berekening te bieden van de radiatieve efficiëntie en accretiesnelheden voor sferische PBH-accretie over het sub-maan massabereik. Door de dunne-schijfaannames te vervangen door een microfysische behandeling van sferische flow, laten de auteurs zien dat:

• De radiatieve efficiëntie aanzienlijk lager is dan voorheen aangenomen wordt.
• Er geen Eddington-gelimiteerde fase is voor sferische accretie; de flow blijft super-Eddington in termen van massavloed zonder te worden begrensd door stralingsdruk.
• Bijgevolg zijn de restricties op PBH donkere materie uit stellaire observaties veel strenger dan voorheen gedacht, aangezien PBH's met initiële massa's zo laag als $10^{-16} M_\odot$ kunnen groeien tot ze hun gaststerren consumeren binnen een Hubble-tijd.

De auteurs merken op dat hun werk zich richt op de accretie-microfysica en de terugwerking op de globale stellaire structuur niet bevat (wat koppeling aan een stellaire evolutiecode zoals MESA zou vereisen), maar zij leveren wel de noodzakelijke $\eta(M_{BH})$ en $\dot{M}(M_{BH})$ relaties om dergelijke toekomstige studies mogelijk te maken. Zij stellen ook de discussie over schijfvorming (die voorkomt bij hogere massa's of sneller roterende sterren) uit naar een begeleidend artikel.