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

Resumen Técnico: Chronos: I. Un modelo jerárquico bayesiano de edad-litio: Validación en el cúmulo de las Pléyades

Enunciado del Problema
Establecer una escala de edad estelar absoluta y autoconsistente sigue siendo un desafío central en astrofísica, particularmente debido a la falta de calibración cruzada sistemática entre diferentes técnicas de datación, que a menudo arrojan discrepancias de hasta un factor de dos. Si bien el límite de agotamiento del litio (LDB) proporciona un cronómetro robusto y fundamentado físicamente para poblaciones estelares de baja masa, su aplicación está actualmente limitada por escalas de temperatura heterogéneas y dispersión astrofísica. Específicamente, la transición entre enanas ultrafrías (UCDs) totalmente convectivas y estrellas FGK parcialmente radiativas a menudo se trata con cortes de temperatura abruptos, fallando en propagar las incertidumbres de manera consistente a través de la transición de régimen. Además, la población FGK exhibe una dispersión significativa del litio impulsada por la rotación estelar, donde los rotadores rápidos muestran abundancias de litio aumentadas y mayor dispersión en comparación con los rotadores lentos. Los métodos estándar a menudo se basan en prescripciones de valor único o ajustes empíricos que no capturan completamente esta diversidad intrínseca ni la física no lineal del agotamiento del litio.

Metodología
Los autores presentan Chronos, un modelo jerárquico bayesiano (BHM) diseñado para inferir edades de poblaciones estelares modelando conjuntamente las temperaturas efectivas ($T_{\text{eff}}$), las abundancias de litio ($A_{\text{Li}}$) y la edad global del cúmulo. El marco integra tres componentes principales:

1. Estructura Jerárquica: El modelo trata los datos observados ($T_{\text{eff,obs}}, A_{\text{Li,obs}}$) como realizaciones ruidosas de parámetros estelares verdaderos latentes. Evita límites rígidos entre regímenes estelares, utilizando en su lugar una función logística suave para modelar la transición entre estrellas totalmente convectivas y parcialmente radiativas.
2. Surrogado de Red Neuronal: Un perceptrón multicapa (MLP) entrenado en modelos evolutivos BT-Settl sirve como sustituto de las curvas teóricas de agotamiento del litio. Esta red predice la abundancia media teórica de litio ($\mu_{A_{\text{Li}}}$) en función de la edad y $T_{\text{eff}}$, proporcionando una aproximación rápida y continua de los modelos físicos.
3. Modelo de Mezcla Consciente de la Rotación: Para tener en cuenta la mejora del litio inducida por la rotación en estrellas FGK, el modelo emplea una distribución de mezcla de dos componentes para la abundancia verdadera de litio. Esto consiste en un componente de "rotador lento" (modelado como una distribución Normal centrada en la predicción teórica) y un componente de "rotador rápido" (modelado como una distribución Normal Sesgada con un desplazamiento de amplitud $\Delta A_{\text{Li,rot}}$ y dispersión adicional). La probabilidad de que una estrella pertenezca a la población de rotadores rápidos se trata como una función de su régimen espectral y un parámetro global de fracción de rotación ($p_{\text{rot}}$).

La inferencia se realiza utilizando Cadena de Markov Monte Carlo (MCMC) con el Muestreador No-U-Turn (NUTS) implementado en PyMC. El modelo fue validado utilizando conjuntos de datos sintéticos generados a partir de las mismas suposiciones jerárquicas y aplicado a un conjunto de datos homogéneo de 122 miembros del cúmulo de las Pléyades.

Contribuciones Clave

• Marco Unificado: Chronos reformula la cronología estelar basada en litio como un problema de inferencia jerárquico coherente, permitiendo la estimación simultánea de la edad del cúmulo, la dispersión intrínseca y los parámetros de rotación sin imponer umbrales de temperatura arbitrarios.
• Transición de Régimen Probabilística: El modelo maneja explícitamente la transición entre UCDs y estrellas FGK mediante una función probabilística continua, asegurando una propagación consistente de la incertidumbre a través del LDB.
• Modelado de Rotación: Al incorporar un modelo de mezcla para los efectos de la rotación, el marco da cuenta de la dispersión asimétrica del litio observada en estrellas FGK, distinguiendo entre la varianza astrofísica intrínseca y la mejora inducida por la rotación.
• Modelado de Surrogados: El uso de una red neuronal entrenada en modelos evolutivos permite la evaluación eficiente de la física del agotamiento del litio dentro del bucle de inferencia bayesiana, cerrando la brecha entre rejillas teóricas complejas y modelado estadístico.

Resultados
Aplicado al cúmulo de las Pléyades, Chronos produce una distribución de edad posterior centrada en Edad = 124.53$^{+3.34}_{-2.70}$ Myr. Este resultado es consistente con las estimaciones clásicas del LDB y los resultados recientes de ajuste de isocronas, pero presenta una restricción de incertidumbre más ajustada. El modelo recupera con éxito la edad de entrada en pruebas sintéticas en un rango de 20–400 Myr, demostrando robustez incluso con números limitados de objetos en el régimen sensible al LDB ($T_{\text{eff}} < 3700$ K).

Las verificaciones de diagnóstico, incluidas las estadísticas de convergencia ($\hat{R} \approx 1$), el análisis de residuos y las verificaciones predictivas posteriores (PPCs), confirman la estabilidad del modelo. El análisis revela que la edad inferida está impulsada principalmente por estrellas cerca de la transición LDB pronunciada ($T_{\text{eff}} \sim 2800$ K) y permanece en gran medida no correlacionada con los parámetros de rotación, lo que indica que la determinación de la edad es robusta frente a la parametrización específica de la mejora del litio inducida por la rotación. El modelo también restringe con éxito la fracción de rotadores rápidos ($p_{\text{rot}} \approx 0.32$) y la magnitud de la mejora del litio ($\Delta A_{\text{Li,rot}} \approx 1.78$ dex).

Significado
El artículo posiciona a Chronos como un paso metodológicamente robusto hacia una escala de edad estelar absoluta y autoconsistente. Al demostrar que la cronología basada en litio puede formularse como un problema de inferencia jerárquico, el trabajo proporciona un marco flexible y estadísticamente riguroso para la determinación de edades en poblaciones estelares jóvenes (1–600 Myr). Los autores enfatizan que este enfoque permite la evaluación autoconsistente del poder de restricción relativo de múltiples diagnósticos de edad. Si bien la validación actual se limita a las Pléyades y modelos de metalicidad solar, el marco está diseñado para ser extensible e incluir la metalicidad como una dimensión e incorporar observables adicionales sensibles a la edad en iteraciones futuras. El trabajo no afirma resolver todas las discrepancias de datación, pero establece una base físicamente interpretable y estadísticamente sólida para integrar diagnósticos de litio en una escala de edad más amplia y unificada.