Chronos: Towards a self-consistent and absolute stellar age scale. I. A Bayesian hierarchical lithium-age model: Validation on the Pleiades cluster

Technische Zusammenfassung: Chronos: I. Ein bayessches hierarchisches Lithium-Altersmodell: Validierung am Plejaden-Sternhaufen

Problemstellung
Die Etablierung einer konsistenten und absoluten Sternalter-Skala bleibt eine zentrale Herausforderung in der Astrophysik, insbesondere aufgrund des Fehlens einer systematischen Kreuzkalibrierung zwischen verschiedenen Altersbestimmungsmethoden, die häufig Diskrepanzen von bis zu einem Faktor zwei aufweisen. Während die Lithium-Verarmungsgrenze (LDB) eine robuste, physikalisch fundierte Chronometrie für Populationen von massearmen Sternen bietet, ist ihre Anwendung derzeit durch heterogene Temperaturskalen und astrophysikalische Dispersion eingeschränkt. Insbesondere wird der Übergang zwischen vollständig konvektiven ultra-kühlen Zwergen (UCDs) und teilweise radiativen FGK-Sternen oft mit scharfen Temperaturgrenzen behandelt, was eine konsistente Fortpflanzung von Unsicherheiten über den Regime-Übergang hinweg verhindert. Darüber hinaus zeigt die FGK-Population eine signifikante Lithium-Dispersion, die durch die Sternrotation angetrieben wird, wobei schnelle Rotatoren im Vergleich zu langsamen Rotatoren erhöhte Lithium-Häufigkeiten und eine größere Streuung aufweisen. Standardmethoden verlassen sich häufig auf einwertige Vorschriften oder empirische Anpassungen, die diese inhärente Vielfalt oder die nichtlineare Physik der Lithium-Verarmung nicht vollständig erfassen.

Methodik
Die Autoren stellen Chronos vor, ein bayessches hierarchisches Modell (BHM), das entwickelt wurde, um Alter von Sternpopulationen abzuleiten, indem es gleichzeitig die effektiven Temperaturen ($T_{\text{eff}}$), die Lithium-Häufigkeiten ($A_{\text{Li}}$) und das globale Alter des Sternhaufens modelliert. Das Framework integriert drei Kernkomponenten:

1. Hierarchische Struktur: Das Modell behandelt beobachtete Daten ($T_{\text{eff,obs}}, A_{\text{Li,obs}}$) als verrauschte Realisierungen latenter wahrer Sternparameter. Es vermeidet harte Grenzen zwischen Sternregimen und verwendet stattdessen eine glatte logistische Funktion, um den Übergang zwischen vollständig konvektiven und teilweise radiativen Sternen zu modellieren.
2. Neuronales Netz als Surrogat: Ein auf BT-Settl-Entwicklungsmodellen trainiertes Multilayer-Perceptron (MLP) dient als Surrogat für theoretische Lithium-Verarmungskurven. Dieses Netzwerk sagt die mittlere theoretische Lithium-Häufigkeit ($\mu_{A_{\text{Li}}}$) als Funktion von Alter und $T_{\text{eff}}$ voraus und bietet eine schnelle und kontinuierliche Näherung der physikalischen Modelle.
3. Rotationsbewusstes Mischungsmodell: Um die rotationsinduzierte Lithium-Anreicherung bei FGK-Sternen zu berücksichtigen, verwendet das Modell eine zweikomponentige Mischungsverteilung für die wahre Lithium-Häufigkeit. Diese besteht aus einer Komponente für „langsame Rotatoren" (modelliert als Normalverteilung, zentriert auf die theoretische Vorhersage) und einer Komponente für „schnelle Rotatoren" (modelliert als Schief-Normalverteilung mit einer Amplitudenverschiebung $\Delta A_{\text{Li,rot}}$ und zusätzlicher Dispersion). Die Wahrscheinlichkeit, dass ein Stern zur Population der schnellen Rotatoren gehört, wird als Funktion seines Spektralregimes und eines globalen Rotationsanteil-Parameters ($p_{\text{rot}}$) behandelt.

Die Inferenz wird mittels Markov-Chain-Monte-Carlo (MCMC) mit dem No-U-Turn-Sampler (NUTS), implementiert in PyMC, durchgeführt. Das Modell wurde mit synthetischen Datensätzen validiert, die unter denselben hierarchischen Annahmen generiert wurden, und auf einen homogenen Datensatz von 122 Mitgliedern des Plejaden-Sternhaufens angewendet.

Hauptbeiträge

• Einheitliches Framework: Chronos fasst die auf Lithium basierende Sternchronologie als kohärentes hierarchisches Inferenzproblem neu, was die gleichzeitige Schätzung des Haufenalters, der intrinsischen Dispersion und der Rotationsparameter ohne die Auferlegung willkürlicher Temperaturschwellen ermöglicht.
• Probabilistischer Regime-Übergang: Das Modell behandelt den Übergang zwischen UCDs und FGK-Sternen explizit über eine kontinuierliche probabilistische Funktion und gewährleistet so eine konsistente Fortpflanzung von Unsicherheiten über die LDB hinweg.
• Rotationsmodellierung: Durch die Einbeziehung eines Mischungsmodells für Rotationseffekte berücksichtigt das Framework die beobachtete asymmetrische Lithium-Dispersion bei FGK-Sternen und unterscheidet zwischen intrinsischer astrophysikalischer Varianz und rotationsinduzierter Anreicherung.
• Surrogat-Modellierung: Die Verwendung eines auf Entwicklungsmodellen trainierten neuronalen Netzes ermöglicht eine effiziente Bewertung der Physik der Lithium-Verarmung innerhalb der bayesschen Inferenzschleife und schließt die Lücke zwischen komplexen theoretischen Gittern und statistischer Modellierung.

Ergebnisse
Angewendet auf den Plejaden-Sternhaufen liefert Chronos eine posteriori-Altersverteilung mit einem Zentrum bei Alter = 124,53$^{+3,34}_{-2,70}$ Myr. Dieses Ergebnis ist konsistent mit klassischen LDB-Schätzungen und jüngsten Ergebnissen der Isochronen-Anpassung, weist jedoch eine engere Unsicherheitsbeschränkung auf. Das Modell recovert das Eingabealter in synthetischen Tests erfolgreich über einen Bereich von 20–400 Myr und demonstriert Robustheit selbst bei begrenzten Objektzahlen im LDB-empfindlichen Regime ($T_{\text{eff}} < 3700$ K).

Diagnostische Checks, einschließlich Konvergenzstatistiken ($\hat{R} \approx 1$), Residuenanalyse und posteriori-prädiktive Checks (PPCs), bestätigen die Stabilität des Modells. Die Analyse zeigt, dass das abgeleitete Alter primär durch Sterne in der Nähe des steilen LDB-Übergangs ($T_{\text{eff}} \sim 2800$ K) getrieben wird und weitgehend unkorreliert mit den Rotationsparametern bleibt, was darauf hindeutet, dass die Altersbestimmung robust gegenüber der spezifischen Parametrisierung der rotationsinduzierten Lithium-Anreicherung ist. Das Modell schränkt zudem erfolgreich den Anteil der schnellen Rotatoren ($p_{\text{rot}} \approx 0,32$) und das Ausmaß der Lithium-Anreicherung ($\Delta A_{\text{Li,rot}} \approx 1,78$ dex) ein.

Bedeutung
Der Artikel positioniert Chronos als methodisch robusten Schritt hin zu einer konsistenten und absoluten Sternalter-Skala. Durch den Nachweis, dass die auf Lithium basierende Chronologie als hierarchisches Inferenzproblem formuliert werden kann, bietet die Arbeit ein flexibles und statistisch rigoroses Framework für die Altersbestimmung in jungen Sternpopulationen (1–600 Myr). Die Autoren betonen, dass dieser Ansatz eine konsistente Bewertung der relativen einschränkenden Kraft multipler Altersdiagnostiken ermöglicht. Während die aktuelle Validierung auf den Plejaden und Modelle mit solarer Metallizität beschränkt ist, ist das Framework so konzipiert, dass es erweiterbar ist, um Metallizität als Dimension einzubeziehen und in zukünftigen Iterationen zusätzliche alterssensitive Observablen zu integrieren. Die Arbeit beansprucht nicht, alle Diskrepanzen bei der Altersbestimmung zu lösen, sondern legt eine physikalisch interpretierbare, statistisch fundierte Grundlage für die Integration von Lithium-Diagnostiken in eine breitere, einheitliche Altersskala.