Radiating Bondi Flows I: Dimensionless Framework and Constant Opacity Solutions

Titel: Radiating Bondi Flows I: Dimensionless Framework and Constant Opacity Solutions

Auteurs: Avery P. Bailey, Andrew N. Youdin, en Kaitlin M. Kratter
Datum: 24 maart 2026 (Conceptversie)

---

1. Probleemstelling

De klassieke theorie van Bondi-accretie (H. Bondi, 1952) beschrijft de sferisch symmetrische, stationaire accretie van gas op een zwaartekrachtsbron onder de aanname van adiabatische thermodynamica. Hoewel deze theorie fundamenteel is voor vele astrofysische scenario's, negeert ze stralingsfeedback.

Bestaande uitbreidingen naar stralingsgedreven accretie zijn voornamelijk gericht op omgevingen waar stralingsdruk dominant is (zoals bij zwarte gaten of in AGN-schijven). Er is echter een gebrek aan theoretisch werk voor gas-druk-gedomineerde omgevingen (waar thermische druk de stralingsdruk overtreft), zoals protostellair en protoplanetaire schijven. In deze koude omgevingen ($\sim 100$ K) kan radiatieve verwarming en afkoeling de accretie significant beïnvloeden, maar de dynamica is anders dan in stralings-gedomineerde regimes.

Het doel van dit werk is om de theorie van niet-adiabatische Bondi-accretie uit te breiden voor gas-druk-gedomineerde systemen, met name om te kwantificeren hoe radiatieve feedback de accretiesnelheid ($\dot{M}$) onderdrukt ten opzichte van de adiabatische snelheid ($\dot{M}_{ad}$).

2. Methodologie

De auteurs ontwikkelen een dimensioneel raamwerk en numerieke oplossingen voor stationaire, sferisch symmetrische accretie met radiatieve verwarming en afkoeling.

• Fysische Aannames:

  • Sferische symmetrie en stationaire toestand.
  • Ideaal gas met een adiabatische index $\gamma = 7/5$ (diatomisch gas).
  • Constante opaciteit: $\kappa_P = \kappa_R = \kappa$ (Planck- en Rosseland-middelpunt zijn gelijk en constant).
  • Gas-druk domineert stralingsdruk (geen stralingsdrukkrachten in de impulsvergelijking).
  • Stralingsvervoer wordt gemodelleerd via een diffusiebenadering die overgaat in vrije straling, met behoud van de energiebalans.

• Dimensionale Parameters:
  Het probleem wordt gereduceerd tot vier dimensieloze parameters:

  1. $\gamma$: Adiabatische index (vastgesteld op $7/5$).
  2. $\tau_B$: Optische diepte door de Bondi-radius ($\kappa \rho_\infty r_B$).
  3. $\tilde{L}_\infty$: Dimensieloze luminostiteit op oneindig ($L_\infty / L_B$).
  4. $\beta$: Dimensieloze afkoeltijd ($t_{cool} c_{s,\infty} / r_B$).

• Numerieke Oplossing:
  De auteurs integreren een stelsel gekoppelde differentiaalvergelijkingen voor dichtheid, Mach-getal, entropie, luminostiteit en stralingsenergie-dichtheid.

  • Sonic punt-probleem: Het sonic punt (waar de stroming supersonisch wordt) is niet a priori bekend in het radiatieve geval. De auteurs gebruiken een schietmethode (shooting method) waarbij ze de accretiesnelheid iteratief aanpassen totdat de oplossing het sonic punt correct passeert (Type I oplossing).
  • Randvoorwaarden: Integratie begint bij een grote straal ($r \approx 10^{12} r_B$) met voorwaarden op oneindig. Om numerieke instabiliteit bij grote stralen te voorkomen (waar radiatieve evenwicht geldt), wordt eerst een evenwichtsoplossing gebruikt en later overgeschakeld naar het volledige niet-evenwichtssysteem.

3. Belangrijkste Bijdragen en Resultaten

De studie levert een uitgebreide kaart van de accretie-efficiëntie $f_{acc} = \dot{M} / \dot{M}_{ad}$ over een parameter ruimte van $\tau_B, \tilde{L}_\infty, \beta \in [10^{-3}, 10^3]$.

A. Onderdrukking van Accretie

In tegenstelling tot stralings-gedomineerde omgevingen, waar straling de accretie kan versnellen, werkt radiatieve feedback in gas-gedomineerde omgevingen primair als een onderdrukkende factor.

• Verwarming van het invallende gas verhoogt de geluidssnelheid, wat de sonic radius naar binnen verplaatst en de effectieve accretie-oppervlakte verkleint.
• De onderdrukking is het sterkst bij hoge optische diepte ($\tau_B$), hoge luminostiteit ($\tilde{L}_\infty$) en lange afkoeltijden ($\beta$).

B. Analytische Schaalwetten

De auteurs leiden analytische schaalwetten af voor verschillende regimes (onder de aanname van radiatief evenwicht, geldig voor kleine $\beta$):

1. Isotherm Regime (Snelle afkoeling, lage luminostiteit):

   • Voor $\tilde{L}_\infty \lesssim \min(1, 1/\tau_B)$ en snelle afkoeling ($\beta$ klein).
   • De accretie nadert de isotherme limiet: $f_{acc} \approx 2.97$ (ongeveer 3 keer de adiabatische snelheid).

2. Dunne Optische Regime (Hoge luminostiteit, lage $\tau_B$):

   • Voor $\tilde{L}_\infty \gtrsim \max(1, \tau_B)$.
   • De accretie wordt onderdrukt volgens:
     $$f_{acc} \sim \tilde{L}_\infty^{-5/4}$$
   • Dit is onafhankelijk van de opaciteit en afkoeltijd.

3. Dikke Optische Regime (Lage luminostiteit, hoge $\tau_B$):

   • Voor $\tilde{L}_\infty \lesssim \tau_B$ en lage luminostiteit.
   • De accretie wordt onderdrukt volgens:
     $$f_{acc} \sim \tau_B^{-10/11} \beta^{-5/11}$$
   • Opmerking: Deze oplossingen met radiaal variërende luminostiteit zijn vaak niet astrofysisch relevant omdat ze geen zelf-consistente vrije-val-oplossingen toelaten in het binnenste gebied.

4. Intermediair Regime (Hoge $\tau_B$ én hoge $\tilde{L}_\infty$):

   • Voor $\tilde{L}_\infty \gtrsim \tau_B^{5/11} \beta^{8/11}$.
   • De onderdrukking hangt af van beide parameters:
     $$f_{acc} \sim (\tilde{L}_\infty \tau_B)^{-5/8}$$
   • Dit is onafhankelijk van de afkoeltijd $\beta$.

C. Gebied zonder Oplossing

Er bestaat een regio in de parameter ruimte (voornamelijk bij lage luminosties en specifieke $\tau_B$ waarden) waar geen transsonische oplossing bestaat. In dit gebied zou de luminostiteit negatief worden voordat het sonic punt wordt bereikt, wat fysisch onmogelijk is. Astrofysische systemen met realistische accretieluminostiteiten vermijden dit gebied van nature.

D. Convectie en Stabiliteit

De onderdrukte oplossingen vertonen een negatieve entropiegradiënt buiten het sonic punt, wat in statische systemen zou leiden tot convectieve instabiliteit (Schwarzschild-criterium).

• De auteurs analyseren echter dat de advectietijd (tijd voor het gas om naar binnen te stromen) en de afkoeltijd vaak korter zijn dan de convectieve groeitijd.
• Hierdoor is convectie vaak inefficiënt en kan de stralingsdynamica de entropieprofielen handhaven.
• Zelfs als convectie optreedt (bijvoorbeeld in modellen met lange afkoeltijden), heeft het slechts een beperkt effect op de totale accretiesnelheid, tenzij het convectieve gebied zeer diep in de Bondi-radius doordringt.

4. Significantie en Toepassingen

• Planetaire Vorming: Dit werk is direct relevant voor de "runaway" fase van de accretie van gasreuzen in protoplanetaire schijven. In dit stadium is de gasdruk dominant, maar kan de intrinsieke luminostiteit van de jonge planeet of de accretieluminostiteit de accretie significant remmen. Traditionele modellen die een adiabatische snelheid aannemen, kunnen de accretie dus overschatten.
• Analytisch Gereedschap: De auteurs bieden een gedetailleerde analytische formule (zie Appendix B) en een numerieke database die gebruikt kan worden in populatie-synthese modellen om realistischere accretie-snelheden te berekenen zonder dure hydrodynamische simulaties.
• Andere Toepassingen: Hoewel gericht op planeten, kunnen de modellen ook relevant zijn voor sterren in AGN-schijven (waar gasdruk dominant kan zijn) en wind-accretie in dubbelsterstelsels.

Conclusie

Dit artikel legt de theoretische basis voor radierende Bondi-stromen in gas-gedomineerde omgevingen. Het toont aan dat radiatieve feedback een krachtig mechanisme is om accretie te onderdrukken, met schaalwetten die sterk afhangen van de optische diepte en de luminostiteit. De resultaten bieden een cruciale correctie op bestaande adiabatische modellen, met name voor de vorming van gasreuzen.