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

Resumo Técnico: Chronos: I. Um modelo hierárquico bayesiano de idade-lítio: Validação no aglomerado das Plêiades

Declaração do Problema
Estabelecer uma escala de idades estelares absoluta e autoconsistente permanece um desafio central na astrofísica, particularmente devido à falta de calibração cruzada sistemática entre diferentes técnicas de datação, que frequentemente produzem discrepâncias de até um fator de dois. Embora o limite de esgotamento do lítio (LDB) proporcione um cronômetro robusto e fundamentado fisicamente para populações estelares de baixa massa, sua aplicação é atualmente limitada por escalas de temperatura heterogêneas e dispersão astrofísica. Especificamente, a transição entre anãs ultra-frias (UCDs) totalmente convectivas e estrelas FGK parcialmente radiativas é frequentemente tratada com cortes de temperatura abruptos, falhando em propagar consistentemente as incertezas através da transição de regime. Além disso, a população FGK exibe uma dispersão significativa de lítio impulsionada pela rotação estelar, onde rotadores rápidos apresentam abundâncias de lítio aumentadas e maior dispersão em comparação com rotadores lentos. Métodos padrão frequentemente dependem de prescrições de valor único ou ajustes empíricos que não capturam plenamente essa diversidade intrínseca ou a física não linear do esgotamento do lítio.

Metodologia
Os autores apresentam o Chronos, um modelo hierárquico bayesiano (BHM) projetado para inferir idades de populações estelares ao modelar conjuntamente as temperaturas efetivas ($T_{\text{eff}}$), as abundâncias de lítio ($A_{\text{Li}}$) e a idade global do aglomerado. O framework integra três componentes principais:

1. Estrutura Hierárquica: O modelo trata os dados observados ($T_{\text{eff,obs}}, A_{\text{Li,obs}}$) como realizações ruidosas de parâmetros estelares verdadeiros latentes. Ele evita limites rígidos entre regimes estelares, utilizando em vez disso uma função logística suave para modelar a transição entre estrelas totalmente convectivas e parcialmente radiativas.
2. Surrogado de Rede Neural: Um perceptron multicamadas (MLP) treinado em modelos evolutivos BT-Settl serve como um surrogado para curvas teóricas de esgotamento de lítio. Esta rede prevê a abundância média teórica de lítio ($\mu_{A_{\text{Li}}}$) como uma função da idade e de $T_{\text{eff}}$, fornecendo uma aproximação rápida e contínua dos modelos físicos.
3. Modelo de Mistura Consciente da Rotação: Para contabilizar o aumento de lítio induzido pela rotação em estrelas FGK, o modelo emprega uma distribuição de mistura de dois componentes para a abundância verdadeira de lítio. Isto consiste em um componente de "rotador lento" (modelado como uma distribuição Normal centrada na previsão teórica) e um componente de "rotador rápido" (modelado como uma distribuição Normal Assimétrica com um deslocamento de amplitude $\Delta A_{\text{Li,rot}}$ e dispersão adicional). A probabilidade de uma estrela pertencer à população de rotadores rápidos é tratada como uma função de seu regime espectral e de um parâmetro global de fração de rotação ($p_{\text{rot}}$).

A inferência é realizada usando Cadeias de Markov Monte Carlo (MCMC) com o Amostrador No-U-Turn (NUTS) implementado no PyMC. O modelo foi validado usando conjuntos de dados sintéticos gerados a partir das mesmas suposições hierárquicas e aplicado a um conjunto de dados homogêneo de 122 membros do aglomerado das Plêiades.

Principais Contribuições

• Framework Unificado: O Chronos reformula a cronologia estelar baseada em lítio como um problema de inferência hierárquica coerente, permitindo a estimativa simultânea da idade do aglomerado, da dispersão intrínseca e dos parâmetros de rotação sem impor limites de temperatura arbitrários.
• Transição de Regime Probabilística: O modelo trata explicitamente a transição entre UCDs e estrelas FGK através de uma função probabilística contínua, garantindo uma propagação consistente de incertezas através do LDB.
• Modelagem de Rotação: Ao incorporar um modelo de mistura para efeitos de rotação, o framework contabiliza a dispersão assimétrica de lítio observada em estrelas FGK, distinguindo entre variância astrofísica intrínseca e aumento induzido pela rotação.
• Modelagem de Surrogado: O uso de uma rede neural treinada em modelos evolutivos permite a avaliação eficiente da física de esgotamento de lítio dentro do loop de inferência bayesiana, preenchendo a lacuna entre grades teóricas complexas e modelagem estatística.

Resultados
Aplicado ao aglomerado das Plêiades, o Chronos produz uma distribuição de idade posterior centrada em Idade = 124,53$^{+3,34}_{-2,70}$ Myr. Este resultado é consistente com estimativas clássicas de LDB e resultados recentes de ajuste de isócronas, mas apresenta uma restrição de incerteza mais apertada. O modelo recupera com sucesso a idade de entrada em testes sintéticos em uma faixa de 20–400 Myr, demonstrando robustez mesmo com números limitados de objetos no regime sensível ao LDB ($T_{\text{eff}} < 3700$ K).

Verificações de diagnóstico, incluindo estatísticas de convergência ($\hat{R} \approx 1$), análise de resíduos e verificações preditivas posteriores (PPCs), confirmam a estabilidade do modelo. A análise revela que a idade inferida é impulsionada principalmente por estrelas próximas à transição íngreme do LDB ($T_{\text{eff}} \sim 2800$ K) e permanece largely não correlacionada com os parâmetros de rotação, indicando que a determinação da idade é robusta contra a parametrização específica do aumento de lítio induzido pela rotação. O modelo também restringe com sucesso a fração de rotadores rápidos ($p_{\text{rot}} \approx 0,32$) e a magnitude do aumento de lítio ($\Delta A_{\text{Li,rot}} \approx 1,78$ dex).

Significado
O artigo posiciona o Chronos como um passo metodologicamente robusto em direção a uma escala de idades estelares absoluta e autoconsistente. Ao demonstrar que a cronologia baseada em lítio pode ser formulada como um problema de inferência hierárquica, o trabalho fornece um framework flexível e estatisticamente rigoroso para a determinação de idades em populações estelares jovens (1–600 Myr). Os autores enfatizam que esta abordagem permite a avaliação autoconsistente do poder de restrição relativo de múltiplos diagnósticos de idade. Embora a validação atual seja limitada às Plêiades e a modelos de metalicidade solar, o framework é projetado para ser extensível para incluir a metalicidade como uma dimensão e incorporar observáveis adicionais sensíveis à idade em iterações futuras. O trabalho não afirma resolver todas as discrepâncias de datação, mas estabelece uma base fisicamente interpretável e estatisticamente sólida para integrar diagnósticos de lítio em uma escala de idades mais ampla e unificada.