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

Résumé technique : Chronos : I. Un modèle hiérarchique bayésien de l'âge au lithium : Validation sur l'amas des Pléiades

Énoncé du problème
Établir une échelle d'âge stellaire absolue et auto-cohérente demeure un défi central en astrophysique, particulièrement en raison du manque d'étalonnage croisé systématique entre différentes techniques de datation, qui produisent souvent des écarts pouvant atteindre un facteur deux. Bien que la limite d'appauvrissement en lithium (LDB) fournisse un chronomètre robuste et physiquement fondé pour les populations stellaires de faible masse, son application est actuellement limitée par des échelles de température hétérogènes et une dispersion astrophysique. Plus précisément, la transition entre les naines ultra-froides (UCD) entièrement convectives et les étoiles FGK partiellement radiatives est souvent traitée par des coupes nettes en température, ne parvenant pas à propager les incertitudes de manière cohérente à travers la transition de régime. De plus, la population FGK présente une dispersion significative du lithium, pilotée par la rotation stellaire, où les rotateurs rapides affichent des abondances de lithium accrues et une dispersion plus importante par rapport aux rotateurs lents. Les méthodes standard reposent souvent sur des prescriptions à valeur unique ou des ajustements empiriques qui ne capturent pas pleinement cette diversité intrinsèque ni la physique non linéaire de l'appauvrissement en lithium.

Méthodologie
Les auteurs présentent Chronos, un modèle hiérarchique bayésien (BHM) conçu pour inférer les âges des populations stellaires en modélisant conjointement les températures effectives ($T_{\text{eff}}$), les abondances en lithium ($A_{\text{Li}}$) et l'âge global de l'amas. Le cadre intègre trois composants principaux :

1. Structure hiérarchique : Le modèle traite les données observées ($T_{\text{eff,obs}}, A_{\text{Li,obs}}$) comme des réalisations bruitées de paramètres stellaires latents véritables. Il évite les frontières rigides entre les régimes stellaires, utilisant à la place une fonction logistique lisse pour modéliser la transition entre les étoiles entièrement convectives et partiellement radiatives.
2. Surrogate de réseau neuronal : Un perceptron multicouche (MLP) entraîné sur des modèles d'évolution BT-Settl sert de substitut aux courbes théoriques d'appauvrissement en lithium. Ce réseau prédit l'abondance moyenne théorique en lithium ($\mu_{A_{\text{Li}}}$) en fonction de l'âge et de $T_{\text{eff}}$, fournissant une approximation rapide et continue des modèles physiques.
3. Modèle de mélange sensible à la rotation : Pour tenir compte de l'augmentation du lithium induite par la rotation chez les étoiles FGK, le modèle emploie une distribution de mélange à deux composantes pour l'abondance réelle en lithium. Celle-ci comprend une composante « rotateur lent » (modélisée comme une distribution Normale centrée sur la prédiction théorique) et une composante « rotateur rapide » (modélisée comme une distribution Skew-Normale avec un décalage d'amplitude $\Delta A_{\text{Li,rot}}$ et une dispersion supplémentaire). La probabilité qu'une étoile appartienne à la population des rotateurs rapides est traitée comme une fonction de son régime spectral et d'un paramètre global de fraction de rotation ($p_{\text{rot}}$).

L'inférence est réalisée à l'aide de la méthode de Monte Carlo par chaînes de Markov (MCMC) avec l'échantillonneur No-U-Turn (NUTS) implémenté dans PyMC. Le modèle a été validé à l'aide de jeux de données synthétiques générés à partir des mêmes hypothèses hiérarchiques et appliqué à un jeu de données homogène de 122 membres de l'amas des Pléiades.

Contributions clés

• Cadre unifié : Chronos reformule la chronologie stellaire basée sur le lithium comme un problème d'inférence hiérarchique cohérent, permettant l'estimation simultanée de l'âge de l'amas, de la dispersion intrinsèque et des paramètres de rotation sans imposer de seuils de température arbitraires.
• Transition de régime probabiliste : Le modèle gère explicitement la transition entre les UCD et les étoiles FGK via une fonction probabiliste continue, assurant une propagation cohérente des incertitudes à travers la LDB.
• Modélisation de la rotation : En incorporant un modèle de mélange pour les effets de rotation, le cadre prend en compte la dispersion asymétrique observée du lithium chez les étoiles FGK, distinguant la variance astrophysique intrinsèque de l'augmentation induite par la rotation.
• Modélisation par substitut : L'utilisation d'un réseau neuronal entraîné sur des modèles d'évolution permet une évaluation efficace de la physique de l'appauvrissement en lithium au sein de la boucle d'inférence bayésienne, comblant le fossé entre les grilles théoriques complexes et la modélisation statistique.

Résultats
Appliqué à l'amas des Pléiades, Chronos produit une distribution d'âge a posteriori centrée sur Âge = 124,53$^{+3,34}_{-2,70}$ Myr. Ce résultat est cohérent avec les estimations classiques de la LDB et les résultats récents d'ajustement d'isochrones, mais présente une contrainte d'incertitude plus serrée. Le modèle récupère avec succès l'âge d'entrée dans les tests synthétiques sur une plage de 20 à 400 Myr, démontrant sa robustesse même avec un nombre limité d'objets dans le régime sensible à la LDB ($T_{\text{eff}} < 3700$ K).

Les vérifications diagnostiques, incluant les statistiques de convergence ($\hat{R} \approx 1$), l'analyse des résidus et les vérifications prédictives a posteriori (PPC), confirment la stabilité du modèle. L'analyse révèle que l'âge inféré est principalement piloté par les étoiles proches de la transition raide de la LDB ($T_{\text{eff}} \sim 2800$ K) et reste largement non corrélé aux paramètres de rotation, indiquant que la détermination de l'âge est robuste face à la paramétrisation spécifique de l'augmentation du lithium induite par la rotation. Le modèle contraint également avec succès la fraction de rotateurs rapides ($p_{\text{rot}} \approx 0,32$) et l'amplitude de l'augmentation du lithium ($\Delta A_{\text{Li,rot}} \approx 1,78$ dex).

Signification
L'article positionne Chronos comme une étape méthodologiquement robuste vers une échelle d'âge stellaire absolue et auto-cohérente. En démontrant que la chronologie basée sur le lithium peut être formulée comme un problème d'inférence hiérarchique, ce travail fournit un cadre flexible et statistiquement rigoureux pour la détermination de l'âge dans les populations stellaires jeunes (1–600 Myr). Les auteurs soulignent que cette approche permet l'évaluation auto-cohérente du pouvoir contraignant relatif de multiples diagnostics d'âge. Bien que la validation actuelle se limite aux Pléiades et aux modèles à métallicité solaire, le cadre est conçu pour être extensible afin d'inclure la métallicité comme dimension et d'intégrer des observables supplémentaires sensibles à l'âge dans les itérations futures. Ce travail ne prétend pas résoudre tous les écarts de datation, mais établit une fondation physiquement interprétable et statistiquement solide pour intégrer les diagnostics au lithium dans une échelle d'âge plus large et unifiée.