Ages and masses of asymptotic giant branch stars from the period--luminosity diagram

Titel: Leeftijden en massa's van asymptotische reuzentak-sterren (AGB) afgeleid uit het periode-luminositeit-diagram

1. Het Probleem

Asymptotische reuzentak-sterren (AGB) vertegenwoordigen de laatste kernbrandingsfase van sterren met lage en intermediaire massa. Ze spelen een cruciale rol in de chemische verrijking van het interstellair medium (ISM) door middel van massaverlies. Echter, het bepalen van de leeftijd en de massa (zowel de initiële als de huidige massa) van individuele AGB-sterren is extreem moeilijk.

• Bestaande methoden: Deze zijn vaak afhankelijk van het Hertzsprung-Russell-diagram en de relatie tussen kernmassa en luminositeit. De envelopmassa is echter sterk onzeker door de degeneratie met metaliciteit, waardoor nauwkeurige bepalingen van massa en leeftijd voor veldsterren vrijwel onmogelijk zijn.
• De uitdaging: De complexiteit van de AGB-fase, inclusief thermische pulsen, derden-dredge-up (het naar boven brengen van koolstof) en massaverlies, maakt modellering moeilijk. Er is behoefte aan een methode die gebruikmaakt van observabele grootheden om statistisch betrouwbare schattingen te geven voor populaties van sterren.

2. Methodologie

De auteur presenteert een probabilistische methode om de leeftijd en massa van AGB-sterren te bepalen door hun positie in het periode-luminositeit ($P-L$) diagram te vergelijken met theoretische evolutionaire modellen.

• Modellering: Er wordt gebruikgemaakt van een raster van evolutionaire tracks van het Padova CMD (Colibri S_37 AGB tracks, Bressan et al. 2012; Pastorelli et al. 2020), aangevuld met niet-lineaire pulsatiemodellen (Trabucchi et al. 2021). Deze modellen simuleren de pulsatieperioden (fundamentele modus en overtonen) voor sterren met verschillende initiële massa's ($0.8 - \sim6 M_\odot$) en metaliciteiten.
• Probabilistische Weegfactoren: Voor elk punt in het $P-L$-diagram wordt een waarschijnlijkheidswaarde ($w_{ip}$) berekend die rekening houdt met:
  • De Initiële Massa Functie (IMF) (Kroupa 2001).
  • De Stervormingsgeschiedenis (SFH) van de Melkweg (gebaseerd op Gaia-data binnen 100 pc van de Zon, Alzate et al. 2021).
  • De tijdsintervallen tussen isochronen en de dichtheid van sterren in het diagram.
• Berekening: Voor een waargenomen ster met een periode ($P_{obs}$) en een luminositeitproxy ($L_{obs}$, vaak een Wesenheit-index of magnitude), wordt de relatieve waarschijnlijkheid berekend dat de ster een bepaalde eigenschap (bijv. leeftijd of massa) heeft. Dit resulteert in een statistische verdeling in plaats van een enkele waarde.
• Selectiecriteria: De methode houdt rekening met observatiebias, zoals limieten in pulsatieperiode (bijv. $P > 60$ dagen voor LPV's), amplitude, en chemische samenstelling (O-rijk vs. C-rijk).

3. Belangrijkste Bijdragen

• Probabilistische Raamwerk: Dit is het eerste werk dat volledig probabilistische berekeningen toepast op de $P-L$-diagrammen van AGB-sterren, inclusief generieke code die door anderen kan worden gebruikt.
• Integratie van Pulsatie: Het gebruik van recente, niet-lineaire pulsatiemodellen in de Padova-isochronen maakt het mogelijk om de pulsatievolgorde (A, B, C, C', D) nauwkeurig te reproduceren, wat essentieel is voor het koppelen van observaties aan evolutionaire status.
• Kwantificering van Bias: De studie analyseert systematisch hoe verschillende surveys (Gaia, NESS, DEATHSTAR, ATOMIUM) de waargenomen populatie van AGB-sterren vertekenen ten opzichte van de theoretische verwachtingen.

4. Resultaten

A. Statistische Verdelingen (Melkweg)

• De typische AGB-ster die de Melkweg verlaat, heeft een initiële massa van ongeveer $1.1 M_\odot$.
• Ongeveer 95% van de huidige AGB-sterren heeft een initiële massa onder de $2 M_\odot$.
• Wanneer gewogen naar massaverlies (bijdrage aan het ISM), is de mediane initiële massa iets hoger, ongeveer $1.22 M_\odot$. De meeste massaverlies komt echter van sterren met een initiële massa onder de $2.8 M_\odot$.
• De huidige massa is vaak lager dan de initiële massa door massaverlies, maar voor de meeste sterren is dit verlies nog beperkt ten opzichte van de totale levensduur.

B. Nauwkeurigheid en Onzekerheden
Voor steekproeven van Melkweg-sterren worden de volgende foutmarges geschat (statistisch + systematisch):

• Leeftijd: $+29/-35 \pm 15\%$
• Initiële Massa: $+14/-7 \pm 7\%$
• Huidige Massa: $+17/-11 \pm 27\%$
• De methode werkt het beste voor statistische steekproeven. Voor individuele sterren, vooral die met overtone-pulsaties of afwijkende posities in het diagram, is de nauwkeurigheid lager.

C. Analyse van Survey-Steekproeven
De auteur toont aan dat bestaande surveys vaak vertekend zijn:

• Gaia LPV & NESS (300 pc): Deze steekproeven komen het dichtst bij de theoretische verwachtingen, maar tonen een lichte bias naar jongere sterren en hogere massa's (waarschijnlijk door Malmquist-bias of selectie-effecten).
• NESS (Tier 3 & 4), DEATHSTAR en ATOMIUM: Deze surveys zijn sterk vertrokken naar jongere sterren met hogere massa's en extreme massaverliesraten.
  • Oorzaak: Wetenschappelijke interesse in chemisch peculiar sterren (C-rijke, S-type) en sterren met hoge massaverliesraten, die per definitie zeldzamer en jonger zijn.
  • Conclusie: Deze surveys zijn niet representatief voor de algemene populatie van AGB-sterren, maar wel voor de "extreme" bijdragers aan het ISM.

5. Betekenis en Conclusie

De paper concludeert dat de $P-L$-diagrammethode, gekoppeld aan moderne evolutionaire modellen, een krachtig hulpmiddel is om de eigenschappen van AGB-sterrenpopulaties te karakteriseren.

• Beperkingen: De methode is niet geschikt voor het bepalen van de exacte leeftijd van een enkele ster met hoge precisie, vooral niet als de pulsatiemodus onzeker is of als de ster afwijkt van de standaardevolutie (bijv. door binair interactie). Chemische classificatie (C/O-ratio) wordt afgeraden als Bayesian evidence vanwege de onzekerheid in modellering van dredge-up en hot bottom burning.
• Toepassing: De methode is uiterst waardevol voor het vergelijken van verschillende steekproeven en het begrijpen van hoe selectie-effecten de waargenomen populatie beïnvloeden.
• Algemene Inzichten: De meeste bijdrage aan het interstellair medium in de Melkweg komt van sterren met een initiële massa van ongeveer $1.2 M_\odot$, en niet van de "extreme" zware sterren die vaak in de literatuur worden belicht. Dit corrigeert een veelvoorkomend misverstand over de dominantie van extreme massaverliesraten in de totale chemische verrijking van de Melkweg.

Kortom, dit werk biedt een robuust statistisch raamwerk om de "onzichtbare" parameters van AGB-sterren (leeftijd en massa) te ontrafelen op basis van hun pulsatie-eigenschappen, terwijl het tegelijkertijd waarschuwt voor de systematische vertekeningen in bestaande astronomische surveys.