Disk-Regulated Mass Transfer Between Rotating Non-Degenerate Stars: Insights from Be and sdOB Binaries

Technische Samenvatting: Schijf-gereguleerde Massatransfer tussen Roterende Niet-degeneratieve Sterren

Probleemstelling
Massatransfer tussen niet-degeneratieve sterren is een cruciaal maar nog onvoldoende begrepen proces in binaire evolutie. Een centrale onzekerheid is de accretie-efficiëntie—het deel van de getransfereerde massa dat door de accretor wordt behouden. Standaard binaire populatiesynthese (BPS) en gedetailleerde stellaire evolutiesimulaties (bijv. met behulp van MESA) hanteren doorgaans een "rotationeel beperkt" accretie-voorschrift. In dit kader wordt accretie beperkt zodra de accretor zijn kritische afbreekrotatiesnelheid nadert, wat vaak leidt tot zeer lage accretie-efficiënties.

Bestendige observaties van Be-sterren-binairiteit, met name Be-ster + subdwerg O/B-type (sdOB) systemen, suggereren echter aanzienlijk hogere accretie-efficiënties (vaak $\gtrsim 30\%$ tot $>50\%$). Bestaande rotationeel beperkte modellen onderschatten systematisch de massa's van Be-sterren in deze systemen en slagen er niet in het geobserveerde parameterrumte van Be+sdOB-binairiteit te reproduceren. Voorts wijzen eerdere studies erop dat noch rotationeel beperkte, noch thermisch beperkte voorschriften alleen voldoende zijn om aan de efficiëntie-eisen te voldoen die afgeleid zijn van geobserveerde Be+sdOB-populaties.

Methodologie
De auteurs stellen een fysisch gemotiveerd "schijf-ster koppeling" massaa-ccretievoorschrift voor om deze discrepanties aan te pakken. Het kernmechanisme is dat een accretie- of decretieschijf de naar de accretor getransporteerde impulsmoment reguleert, waardoor massaa-ccretie kan doorgaan zelfs wanneer de ster nabij zijn kritische rotatiesnelheid roteert.

1. Regulatie van het impulsmoment: Gebaseerd op de stationaire dunne schijfoplossingen van Popham & Narayan (1991), modelleren de auteurs het specifieke impulsmoment dat door de ster wordt geaccreteerd ($j_{acc}$) als een functie van de stellaire rotatiesnelheid ($\Omega/\Omega_{crit}$).

   • Bij lage rotatiesnelheden is $j_{acc}$ ongeveer de Keplerse waarde aan het stellaire oppervlak.
   • Wanneer de ster een drempelrotatiesnelheid nadert (geadopteerd als $\Omega/\Omega_{crit} \approx 0,9$), zorgt de afschuifspanning van de schijf ervoor dat het impulsmoment naar buiten toe van de ster naar de schijf wordt getransporteerd. Bijgevolg neemt $j_{acc}$ scherp af, mogelijk zelfs negatief te worden, wat voortgezette massaa-ccretie mogelijk maakt zonder dat de ster zijn kritische rotatie overschrijdt.
   • De auteurs implementeren een gladgestreken functionele vorm voor $j_{acc}$ om numerieke stabiliteit in 1D-simulaties te waarborgen, waarbij discontinuïteiten die zouden ontstaan door een strikte drempelwaarde worden vermeden.

2. Behoudswetten en Efficiëntieberekening: Het model gaat ervan uit dat een overtollig impulsmoment dat uit de ster wordt geëxtraheerd, door de schijf wordt getransporteerd en wordt verwijderd via massaverlies aan de buitenrand van de schijf ($R_{out}$). Door de behoud van massa en impulsmoment toe te passen, wordt de accretie-efficiëntie ($\eta = \dot{M}_{acc}/\dot{M}_{tr}$) zelfconsistent berekend. De efficiëntie hangt af van de ratio van het specifieke impulsmoment bij de schijfrand tot dat bij de accretieradius.

3. Simulatiekader: De auteurs hebben dit voorschrift geïmplementeerd binnen het POSYDON-framework, waarbij MESA (Modules for Experiments in Stellar Astrophysics) werd gebruikt om binaire modellen te evolueren.

   • Zij schakelden het impliciete windschema uit dat wordt gebruikt om superkritische rotatie te voorkomen, en vertrouwden in plaats daarvan op expliciete rotationeel versterkte winden en de gereduceerde accretiesnelheid in het superkritische regime om subkritische rotatie te herstellen.
   • Een grid van binaire modellen werd gesimuleerd met een zonale metaliciteit, waarbij de primaire massa's ($2-9 M_\odot$), massaverhoudingen ($0,4-0,95$) en orbitale perioden werden gevarieerd.
   • De studie onderzocht ook de impact van kern-overshoot behandelingen (waarbij de standaard POSYDON-instellingen werden vergeleken met een gereduceerd overshoot-voorschrift) op de uiteindelijke componentmassa's.

Belangrijkste Bijdragen en Resultaten

• Analytische Schattingen: De auteurs hebben analytische expressies voor de accretie-efficiëntie afgeleid. Zij tonen aan dat in het schijf-accretieregime de efficiëntie afhangt van de stellaire rotatie, de massaverhouding en de ratio van de stellaire straal tot de Roche-lobstraal. In het directe impact-regime hangt het af van rotatie en de massaverhouding. De algemene massa-gewogen accretie-efficiënties blijken dicht bij de waarden te liggen die verwacht worden nabij de drempelrotatiesnelheid, waar het specifieke geaccreteerde impulsmoment scherp daalt.
• Vergelijking met Binaire Evolutie:
  • Rotationeel Beperkt Model: In dit standaardmodel bereiken accretoren kritische rotatie na het accreteren van slechts een klein deel van de getransfereerde massa ($\sim 0,2 M_\odot$ in een testgeval), wat leidt tot lage efficiënties ($<10\%$ voor de meeste systemen). Dit resulteert in Be-sterren die systematisch minder massief zijn dan geobserveerd.
  • Schijf-Ster Koppeling Model: Dit model staat de accretor toe aanzienlijk meer massa te accreteren (bijv. $\sim 1,6 M_\odot$ in de snelle fase van het testgeval) met efficiënties die $\sim 40-60\%$ bereiken. De schijf extraheert impulsmoment, waardoor de ster kan uitzetten en accreteren terwijl hij nabij de kritische rotatiedrempel blijft.
• Reproduceren van Observaties:
  • Wanneer toegepast op Be+sdOB-systemen, faalt het rotationeel beperkte model om Be-sterren te produceren met massa's die consistent zijn met observaties, met name voor systemen met een brede baan (Geval B).
  • Het schijf-ster koppeling model produceert Be-sterren met massa's die een significant deel van het geobserveerde bereik beslaan.
  • De auteurs stelden vast dat de uiteindelijke componentmassa's gevoelig zijn voor de gehanteerde overshoot. Een model dat het schijf-ster koppeling voorschrift combineert met gereduceerde overshoot levert componentmassa's voor Be+sdOB-systemen die in nauwere overeenstemming zijn met observaties dan de standaard hoge-overshoot modellen.
• Orbitale Perioden: Hoewel het schijf-ster koppeling model de massaretentie verhoogt, zijn de resulterende orbitale perioden vergelijkbaar met of iets langer dan die in het rotationeel beperkte model. Dit komt doordat de meer conservatieve massatransfer de totale massaverlies uit het systeem vermindert, wat leidt tot extremere massaverhoudingen die de orbitale expansie versnellen, wat het verlies van impulsmoment via de schijfrand gedeeltelijk tegenwerkt.

Betekenis en Claims
Het artikel claimt dat de discrepantie tussen de geobserveerde eigenschappen van Be+sdOB-binairiteit en de huidige simulatiesignificante resultaten kan worden opgelost door een schijf-gereguleerd massatransfermodel te adopteren. De auteurs betogen dat:

1. Rotationeel beperkte accretievoorschriften de massa's van Be-sterren in Be+sdOB-systemen systematisch onderschatten.
2. Een schijf-ster koppeling mechanisme, dat een regulatie van de transport van impulsmoment mogelijk maakt, hogere accretie-efficiënties mogelijk maakt die consistent zijn met observationele beperkingen (specifiek de vereiste van efficiënties $>50\%$ in veel systemen).
3. De combinatie van dit nieuwe accretievoorschrift met gereduceerde kern-overshoot een meer fysisch consistente verklaring biedt voor de geobserveerde massa's en orbitale perioden van Be+sdOB-binairiteit.

De auteurs erkennen dat de werkelijke efficiëntie afhangt van onzekere factoren zoals schijfstructuur, magnetische velden en de specifieke condities van schijfvorming en massaverlies (bijv. via L2-punten of 'hot spots'). Zij concluderen echter dat het schijf-ster koppeling model een noodzakelijk fysisch mechanisme biedt om het bestaan van massieve Be-sterren in brede binaire systemen te verklaren waar deeb depletie door getijdenkrachten (tidal spin-down) inefficiënt is.