The winds of OBA hypergiants and luminous blue variables: Dynamically-consistent atmosphere models reveal multiple wind regimes

Probleemstelling

Massaverlies is een cruciale factor in de evolutie van massieve sterren ($M \gtrsim 8 M_\odot$). OBA-hypergiants en Luminous Blue Variables (LBVs) bevinden zich aan het koele einde van het regime van lijngedreven winden en staan dicht bij de Eddington-grens. Hoewel deze fase kort is, beïnvloedt hun hoge massaverlies de omgeving en de latere evolutie (bijv. naar Wolf-Rayet-sterren of type IIn-supernova's) aanzienlijk.

Bestaande theorieën, zoals de "bistability jumps" voorspeld door Vink et al. (1999), suggereren scherpe discontinuïteiten in het massaverlies ($\dot{M}$) en de terminale windsnelheid ($v_\infty$) veroorzaakt door recombinatie van ijzer (Fe) in de atmosfeer:

1. Eerste bistability jump (1-BiSJ): Bij $T_{\text{eff}} \approx 25$ kK (Fe IV $\to$ Fe III).
2. Tweede bistability jump (2-BiSJ): Bij $T_{\text{eff}} \approx 10$ kK (Fe III $\to$ Fe II).

Echter, recente waarnemingen van B-supergiants (BSGs) tonen geen duidelijke discontinuïteiten, wat suggereert dat de huidige massaverliesformules (zoals die van Vink et al. 2001) mogelijk niet volledig geldig zijn voor sterren die dichter bij de Eddington-grens staan. Er is een gebrek aan hydrodynamisch consistente modellen voor het specifieke parameterbereik van hypergiants en LBVs.

Methodologie

De auteurs gebruikten de hydrodynamisch consistente versie van de atmosfeercode PoWR (PoWRhd) om een reeks modellen te genereren die de fysica van deze sterren nauwkeurig weergeven.

• Modelparameters:
  • Massa: $20.6 M_\odot$.
  • Luminositeit: $\log(L/L_\odot) = 5.58$ (representatief voor hypergiants/LBVs).
  • Chemische samenstelling: Zonne-metalliciteit (met uitzondering van enkele tests bij lagere metalliciteit), licht verrijkt Helium ($X_{\text{He}} = 0.33$).
  • Turbulente druk: Een dieptegelijke turbulentiesnelheid van $v_{\text{turb}} = 20$ km/s (hoewel dit mogelijk een onderestimatie is voor heetere modellen).
  • Vullingfactor: $f_{\text{vol}} = 0.33$ (vastgehouden om clumping-effecten te isoleren).
• Variatie: De temperatuur aan de binnengrens ($T^* = T_{\text{eff}}(\tau=50)$) werd gevarieerd van $\sim 12.5$ kK tot $\sim 38.0$ kK.
• Benadering: Alle sterparameters werden vastgehouden, behalve de temperatuur, om het effect van temperatuurvariatie op de windstructuur te isoleren. De modellen zijn volledig hydrodynamisch consistent, wat betekent dat de windsnelheid en het massaverlies zelfoplossend worden bepaald in plaats van aangevoerd.

Belangrijkste Resultaten

1. Complexiteit van het Massaverlies en Windregimes

In plaats van een eenvoudige trend met twee scherpe sprongen, vonden de auteurs een complex patroon van massaverlies als functie van temperatuur, gekenmerkt door afwisselende "plateaus" (hoge $\dot{M}$) en "dalen" (lage $\dot{M}$). Ze onderscheiden twee hoofdtypen van windoplossingen:

• Dichte winden (Dense-wind solutions): Karakteristiek voor plateaus met hoge massaverliesraten en grote straaltransformatie ($R_t > 200 R_\odot$).
• Verdunde winden (Rarefied-wind solutions): Karakteristiek voor dalen met lagere massaverliesraten en kleinere $R_t$.

2. De Rol van Ijzer (Fe) en Bistabiliteit

• 1-BiSJ (Fe IV $\to$ Fe III): Rond $T^* \approx 27$ kK ($T_{2/3} \approx 22$ kK) treedt een scherpe toename in $\dot{M}$ op, veroorzaakt door de recombinatie van Fe IV naar Fe III. Dit komt overeen met de voorspelde eerste bistability jump, maar de auteurs vinden een plateau in plaats van een geleidelijke daling na de sprong.
• 2-BiSJ (Fe III $\to$ Fe II): Rond $T^* \approx 17$ kK ($T_{2/3} \approx 10$ kK) treedt een tweede grote toename in $\dot{M}$ op door de recombinatie van Fe III naar Fe II.
• Bistabiliteit: Bij beide overgangen vonden de auteurs een echte bistabiliteit van oplossingen (meerdere mogelijke windtoestanden voor dezelfde parameters), hoewel ze in de grafieken vaak één representatieve oplossing kozen voor duidelijkheid.

3. De "Dalen" en de Rol van Druk

Tussen de plateaus bevinden zich dalen waar het massaverlies drastisch daalt (met meer dan een orde van grootte).

• In deze dalen is de lijndekking (opacity) van Fe onvoldoende om een dichte wind te starten.
• Turbulente en gasdruk: In deze verdunde regimes wordt de windversnelling steeds meer ondersteund door gasdruk en turbulentiedruk. Voor de koudste modellen is turbulentiedruk essentieel om een wind te lanceren; zonder deze term zou er geen uitstroom ontstaan.
• De auteurs stellen dat de schakeling tussen dichte en verdunde regimes alleen optreedt binnen een bepaalde nabijheid van de Eddington-grens.

4. Snelheidsstructuur en Spectra

• Snelheid: De terminale snelheid ($v_\infty$) varieert sterk tussen de regimes. In verdunde winden is $v_\infty$ hoger (tot $\sim 1000$ km/s), terwijl deze in dichte winden lager is ($\sim 300$ km/s).
• Spectra: De synthetische spectra van de modellen komen overeen met waarnemingen van bekende hypergiants en LBVs in de respectievelijke temperatuurregimes.
  • Dichte winden vertonen sterke P-Cygni-profielen en emissielijnen (vergelijkbaar met HR Car, BP Cru).
  • Verdunde winden tonen zwakkere hypergiant-kenmerken (vergelijkbaar met late B-supergiants zoals HD 168625).
• Relatie $v_\infty$ - $T_{\text{eff}}$: De empirische relatie tussen terminale snelheid en effectieve temperatuur (Hawcroft et al. 2024) geldt goed voor verdunde winden, maar faalt voor dichte winden waar de relatie veel vlakker wordt.

5. Invloed van Metaliciteit

Tests bij lagere metaliciteit (0.2 - 0.7 $Z_\odot$) tonen aan dat bij zeer lage metaliciteit de lijndekking afneemt. Dit kan leiden tot een schakeling van een dichte naar een verdunde oplossing, zelfs bij temperaturen waar bij zonne-metaliciteit een dichte wind zou worden verwacht. Dit verklaart mogelijk waarom sommige LBVs in het LMC/SMC (Laag Metaliciteit) minder extreme winden vertonen.

Bijdrage en Betekenis

1. Bevestiging en Nuancering van Bistability Jumps: Het artikel bevestigt het bestaan van de eerste en tweede bistability jump in het hypergiant-regime, maar toont aan dat deze zijn ingebed in een complexer gedrag met afwisselende regimes, wat niet wordt vastgelegd door eenvoudige massaverliesformules.
2. Rol van Turbulente Druk: Een cruciale bevinding is dat turbulentiedruk een dominante rol speelt in de windversnelling van koele, massieve sterren, vooral in verdunde regimes en dicht bij de Eddington-grens. Dit onderstreept de noodzaak van geavanceerdere behandelingen van radiatief gedreven turbulentie in 1D-modellen.
3. Verbinding tussen Theorie en Observatie: De modellen sluiten aan bij de spectra van een breed scala aan objecten, van OBA-hypergiants tot LBVs en gele hypergiants (YHGs). Het suggereert dat objecten zoals V509 Cas (een YHG) en MWC 930 (een LBV) fysiek vergelijkbaar zijn, maar zich in verschillende windregimes bevinden afhankelijk van hun temperatuur en dichtheid.
4. Beperkingen van Bestaande Formules: De auteurs concluderen dat bestaande voorspellingen (zoals die van Vink et al. 2001) de complexiteit van de overgangen en het gedrag in de "dalen" niet correct voorspellen. Er is een nieuwe, meer dynamische benadering nodig voor massaverlies in het Eddington-dichtste regime.

Conclusie:
Dit onderzoek levert het eerste hydrodynamisch consistente beeld van de winden van OBA-hypergiants en LBVs. Het onthult dat het massaverlies niet lineair is, maar wordt gedicteerd door een complex samenspel tussen ijzer-recombinatie, Eddington-factoren en turbulentiedruk, wat leidt tot een bimodaal gedrag (dicht vs. verdund) dat sterk afhankelijk is van de temperatuur en metaliciteit.