Radiatively Cooled Binary Mass Transfer: Flow Structure, Luminosities, and L2 Outflows Across Mass Transfer Rates

Technische Samenvatting: Radiatief Gekoelde Massatransfer in Binair Stelsels

Probleemstelling
Massatransfer (MT) in binair stelsels is een kritiek proces in de sterrenontwikkeling, dat voorafgaat aan de vorming van bronnen van zwaartekrachtsgolven en kerninstortings-supernova's. Een grote onzekerheid bij het modelleren van deze systemen is de mate waarin MT conservatief is—specifiek, of de accretor al het overgedragen materiaal vasthoudt of dat massa het systeem verlaat, waarbij hoekmoment (AM) wordt meegenomen. Eerdere werken (bijv. Scherbak et al. 2025, "Paper I") maakten gebruik van hydrodynamische simulaties van snelle, stabiele MT onder adiabatische aannames, waarbij bleek dat materiaal uitstroomt via het buitenste Lagrangepunt ($L_2$), waardoor circumbinaire uitstromen (CBO's) ontstaan. Paper I negeerde echter radiatieve koeling, wat naar verwachting significant wordt bij gematigde MT-rates ($\lesssim 10^{-2} M_\odot \text{yr}^{-1}$). Het ontbreken van koeling in eerdere modellen liet vragen open over de overgang tussen conservatieve en niet-conservatieve MT-regimes, het specifieke AM dat door ontsnappende massa wordt meegenomen, en de observationele eigenschappen (helderheid en temperatuur) van deze transienten.

Methodologie
De auteurs voeren 3D Newtoniaanse hydrodynamische simulaties uit met de PLUTO-code om stabiele MT te modelleren tussen een donorster van $10 M_\odot$ en een accretor met een puntmassa van $5 M_\odot$ in een vaste cirkelvormige baan. In tegenstelling tot Paper I omvat deze studie een benaderend model voor radiatieve koeling.

• Koelvoorschrift: Het koelvermogen per volume wordt geschat onder de aanname dat fotonen loodrecht op het baanvlak ontsnappen. Het model interpoleert opaciteitstabellen (zonne-metalliciteit) uit MESA en lost de gastemperatuur op door gas- en stralingsdruk in evenwicht te brengen. De koelterm werkt om de inwendige energie te verminderen, onderworpen aan veiligheidsbodems om onfysische waarden te voorkomen (bijv. een temperatuurbodem van 2000 K in gebieden met hoge dichtheid en een drukbodem in gebieden met lage dichtheid).
• Simulatiegrid: De auteurs variëren de baanafstand $a$ om een reeks fysische MT-rates ($\dot{M}_{\text{phys}}$) te onderzoeken van $\sim 10^{-5}$ tot $\sim 10^{-1} M_\odot \text{yr}^{-1}$. Vier primaire simulaties worden uitgevoerd: lowest_mdot ($2.0 \times 10^{-5}$), mid_mdot ($2.2 \times 10^{-3}$), high_mdot ($2.3 \times 10^{-2}$), en highest_mdot ($1.5 \times 10^{-1}$).
• Beperkingen: De simulaties omvatten geen expliciete viscositeit, viskeuze verwarming of gedetailleerd stralingsvervoer. De accretor wordt weergegeven door een Plummer-verzachtingspotentiaal. De donorster wordt aangedreven om zijn Roche-lobe te vullen via warmte-injectie op een Kelvin-Helmholtz-tijdschaal.

Belangrijkste Resultaten

1. Stroomstructuur en Conservatisme:

   • Massa stroomt van de donor naar een schijf rond de accretor. Bij hoge MT-rates ($\gtrsim 10^{-3} M_\odot \text{yr}^{-1}$) treden significante, equatoriaal geconcentreerde uitstromen op via $L_2$.
   • Het fractie massa dat door de accretor wordt vastgehouden, $\beta$, neemt scherp af naarmate de MT-rate toeneemt. Voor $\dot{M} \gtrsim 10^{-2} M_\odot \text{yr}^{-1}$ is MT sterk niet-conservatief ($\beta \approx 0$). Voor lagere rates ($\dot{M} \lesssim 10^{-3} M_\odot \text{yr}^{-1}$) is MT grotendeels conservatief ($\beta \gtrsim 0.9$), hoewel er nog steeds een kleine hoeveelheid materiaal via $L_2$ ontsnapt.
   • De overgang naar sterk niet-conservatieve MT treedt ongeveer op wanneer de accretieluminositeit aan de buitenrand van de schijf van de accretor de Eddington-grens overschrijdt. De auteurs merken op dat hun drempel voor niet-conservatief gedrag hoger ligt (met een factor van $\sim 2-5$) dan semianalytische voorspellingen van Lu et al. (2023).

2. Hoekmoment en Energiebalans:

   • In alle gevallen asymptoteert het specifieke AM ($h_{\text{loss}}$) van het uitstromende gas naar een waarde die ongeveer gelijk is aan die van het $L_2$-punt ($h_{L2}$), ongeacht de MT-rate. Dit impliceert dat zelfs kleine uitstromen efficiënt AM uit het binair stelsel onttrekken.
   • De asymptotische radiale snelheid van de uitstroom neemt toe naarmate de MT-rate afneemt (van $\sim 10\%$ van de baansnelheid bij hoge rates tot $\sim 50\%$ bij lagere rates). Echter, in fysische eenheden blijven de snelheden voor de meeste simulaties in het bereik van $30-60 \text{ km s}^{-1}$.
   • Bernoulli-parameters geven aan dat uitstromen bij hoge MT-rates gebonden zijn, terwijl die bij lagere rates marginaal ongebonden zijn.

3. Luminositeiten en Observationeel Uiterlijk:

   • De koelluminositeit van het gas ($L$) en de temperatuur ($T$) nemen toe met de MT-rate.
   • Bij lage MT-rates wordt koeling gedomineerd door de schijf van de accretor. Bij hoge rates ($\sim 10^{-1} M_\odot \text{yr}^{-1}$) draagt de CBO significant bij, waardoor de totale luminositeit verdubbelt.
   • Piekluminositeiten voor koeling bereiken $\sim 10^{39} \text{ erg s}^{-1}$ voor de hoogste MT-rate en $\sim 10^{38} \text{ erg s}^{-1}$ voor de hoge MT-rate.
   • Effectieve temperaturen ($T_{\text{eff}}$) variëren van $\sim 5 \times 10^4 \text{ K}$ (UV) in de schijf van de accretor bij hoge rates tot $\sim 10^2 \text{ K}$ (middelinfrarood) in de CBO bij lagere rates. De auteurs waarschuwen dat stofvorming mogelijk is maar moeilijk definitief te voorspellen binnen hun simulatiedomein, aangezien temperaturen die nodig zijn voor condensatie mogelijk pas worden bereikt op stralen groter dan waar hun grid zich uitstrekt.

Betekenis en Beweringen
Het artikel beweert het regime te schetsen waarin radiatieve koeling de dynamiek van stabiele MT significant verandert ten opzichte van adiabatische modellen. Door systematisch MT-rates te variëren, identificeren de auteurs een kritieke drempel ($\sim 10^{-3}$ tot $10^{-2} M_\odot \text{yr}^{-1}$) waarbij het systeem overgaat van conservatief naar sterk niet-conservatief, gedreven door het super-Eddington-karakter van de buitenschijf.

De auteurs stellen dat hun resultaten noodzakelijke invoer leveren voor populatiesyntheseberekeningen die de stabiliteit van MT en de orbitale evolutie van binair compacte objecten voorafgaand aan de samensmelting onderzoeken. Zij benadrukken dat het specifieke AM dat door $L_2$-uitstromen wordt meegenomen, consistent hoog is, wat de langetermijnevolutie van de baan van het binair stelsel zal beïnvloeden. Bovendien bieden de berekende luminositeiten en temperaturen voorspellingen voor het observationeel uiterlijk van transienten geassocieerd met snelle MT, zoals die mogelijk worden waargenomen in systemen zoals SS 433, W Serpentis of V1309 Sco. De auteurs merken bescheiden op dat hun resultaten beperkt zijn tot een enkele massaverhouding en zonne-metalliciteit, en dat toekomstig werk vereist is om stofvorming op grotere stralen te modelleren en viskeuze accretiefysica te incorporeren.