Interaction between the ejecta, the accretion disk, and the secondary star in the recurrent nova system U Sco

Technische Samenvatting: Interactie tussen het ejecta, de accretieschijf en de secundaire ster in het recurrente nova-systeem U Sco

Probleemstelling
Hoewel er aanzienlijke inspanningen zijn geleverd voor het modelleren van de initiële thermonucleaire runaway en de ejectiefasen van recurrente nova's (RN's), blijft de daaropvolgende hydrodynamische interactie tussen het nova-ejecta, de accretieschijf die rond de witte dwerg draait, en de secundaire stellaire begeleider onderbelicht. Dit gat is cruciaal voor het begrijpen van de langetermijn evolutie van RN's, specifiek met betrekking tot de integriteit van de accretieschijf (die de recurrente perioden bepaalt), de massabudget van het binaire systeem (ontsnapping versus behoud), en de potentiële chemische contaminatie van de secundaire ster. Het artikel richt zich op U Sco, een prototype recurrente nova met een massieve witte dwerg, een subreus-begeleider en een korte orbitale periode, om deze dynamica aan te pakken.

Methodologie
De studie maakt gebruik van een hybride modelleringsbenadering die 1D en 3D hydrodynamica combineert:

1. Initiële condities: De explosie en de initiële ejectie worden gemodelleerd met behulp van de bol symmetrische Lagrangiaanse code SHIVA. De resulterende ejectastructuur (massa, straal, dichtheid, temperatuur en snelheidsprofielen) wordt vervolgens in een 3D-domein gemapt.
2. 3D-simulatie: De daaropvolgende evolutie wordt gesimuleerd met behulp van de GADGET-2 smoothed-particle hydrodynamics (SPH) code. Het systeem bevat:
   • Een puntvormige witte dwerg ($1,38 M_\odot$).
   • Een subreus-secundaire ster ($0,88 M_\odot$), waarbij alleen de buitenste lagen ($0,0054 M_\odot$) met SPH worden opgelost om de computationele belasting te beheersen, waarbij de kern als een puntmassa wordt behandeld.
   • Een accretieschijf die rond de witte dwerg draait, getest in verschillende geometrieën (uitgedijnde/flared versus V-vormig) en dichtheidsprofielen (afgeleid van Shakura-Sunyaev versus constant).
   • Nova-ejecta met massa's variërend van $1,1 \times 10^{-6}$ tot $2,1 \times 10^{-6} M_\odot$ en snelheden tot $10.000$ km s$^{-1}$.
3. Experimenteel ontwerp: Een reeks van tien simulaties (Modellen A–H) is uitgevoerd om de effecten van specifieke parameters te isoleren: ejectamassa, ejectasnelheid, ejectadichtheidsprofiel, schijfmassa en schijfgeometrie. Cruciaal is dat deze simulaties orbitale rotatie rond het zwaartepunt van het systeem bevatten, een kenmerk dat vaak wordt weggelaten in eerdere studies van wijde binaire systemen.

Belangrijkste Resultaten
De simulaties onthullen uiteenlopende uitkomsten op basis van de verhouding tussen de ejectamassa en de schijfmassa ($M_{ejecta}/M_{disk}$) en de schijfgeometrie:

• Schijfdisruptie: Voor modellen met uitgedijnde schijven waar $M_{ejecta}/M_{disk} \geq 1$ (Modellen A, B, C, D, G, H), resulteert de nova-explosie in de volledige disruptie en het wegvegen van de accretieschijf. De kinetische energie van het ejecta overstijgt de gravitationele bindingsenergie van de schijf aanzienlijk.
• Schijfoverleving: In contrast hiermee vertonen modellen met een massieve schijf ($M_{ejecta}/M_{disk} = 0,48$, Model E) of V-vormige schijven met een constante hoge initiële dichtheid (Model F) een gedeeltelijk voortbestaan van het schijfmateriaal na de impact.
• Massaverlies en Behoud:
  • Hoge ejectasnelheden (kenmerkend voor RN's, bijv. $10.000$ km s$^{-1}$) drijven een aanzienlijk deel van de gecombineerde schijf- en ejectamassa tot boven de ontsnappingssnelheid van het systeem. In het fiduciaire Model A verlaat $\sim97\%$ van het materiaal het binaire systeem, terwijl slechts $\sim3\%$ gravitationeel gebonden blijft aan de witte dwerg.
  • Dit staat in contrast met eerdere bevindingen voor klassieke nova's (Paper I), waar lagere snelheden leidden tot een hoger behoud.
  • Het verhogen van de ejectamassa of de snelheid verhoogt de totale massa die het systeem verlaat en verlaagt de massa die door de witte dwerg wordt geaccreteerd.
  • Omgekeerd verhogen dichter ejecta (constant versus profiel) en massievere schijven het fractie van het materiaal dat door de witte dwerg wordt behouden, door een verbeterde deceleratie tijdens de botsing.
• Interactie met de secundaire ster:
  • De impact van ejecta en schijfmateriaal op de subreus-metgezel is minimaal. Slechts een zeer klein deel van de deeltjes ($\sim0,07\%$ van de totale ejecta/schijfmassa in Model A) wordt opgenomen in de buitenste lagen van de secundaire ster, wat suggereert dat chemische contaminatie verwaarloosbaar is.
  • Een geringe ejectie van massa uit de buitenste lagen van de subreus wordt waargenomen tijdens de botsing in de latere fase. Een deel van deze secundaire massa wordt uit het systeem geslingerd, terwijl een klein deel wordt geaccreteerd door de witte dwerg.
• Rol van Rotatie: In de wijde baan van U Sco speelt rotatie een marginale rol in het totale massaverlies of behoud van het ejecta/schijf-mengsel ($<1\%$ verschil). Echter, rotatie vermindert de hoeveelheid massa die door de subreus wordt uitgestoten en vervolgens door de witte dwerg wordt geaccreteerd met ongeveer $20\%$.

Betekenis en Claims
De auteurs positioneren dit werk als een noodzakelijke kwalitatieve stap om het gat in het begrip van de post-uitbarsting evolutie van recurrente nova's te vullen. De studie claimt de eerste systematische 3D SPH-analyse te zijn van de interactie tussen nova-ejecta, accretieschijven en begeleiders in een roterend binair systeem met parameters specifiek voor U Sco.

De auteurs concluderen bescheiden dat hun bevindingen kwalitatief zijn vanwege resolutiebeperkingen (specifiek het lage aantal deeltjes voor laag-dens ejecta en de schijf). Desalniettemin stellen zij dat deze resultaten waardevolle richting geven aan toekomstig onderzoek, met name voor:

1. Het identificeren van de condities (massaverhoudingen, geometrie, dichtheid) die leiden tot schijfoverleving versus totale disruptie.
2. Het benadrukken van de gevoeligheid van massaretentie voor ejectasnelheid, een cruciaal onderscheidend kenmerk tussen recurrente en klassieke nova's.
3. Het benadrukken van de noodzaak voor gedetailleerde schijfgeometrie-modellen (uitgediend versus V-vormig) om de langetermijnevolutie van deze systemen nauwkeurig te simuleren.

De auteurs merken op dat hoewel hun simulaties in veel standaardscenario's wijzen op volledige schijfdisruptie, observationele claims van schijfoverleving (bijv. in de 2022 U Sco uitbarsting) kunnen duiden op de aanwezigheid van uitgebreide schijven of andere fysieke configuraties die niet volledig door de huidige parameterruimte worden gevangen.