The age of the Universe from a large sample of the oldest Galactic stars

Resumen Técnico: La edad del Universo a partir de una gran muestra de las estrellas galácticas más antiguas

Planteamiento del problema
El artículo aborda la determinación de la edad del Universo ($A_U$) mediante la estimación de la edad latente de las estrellas más antiguas de la Vía Láctea ($A_\star$). Esta medición sirve como un límite inferior crítico para $A_U$ (donde $A_U = A_\star + t_f$, siendo $t_f \approx 0.2$ Gyr el tiempo de formación de las primeras estrellas de larga vida). El estudio está motivado por la "tensión de Hubble", una discrepancia estadísticamente significativa entre la constante de Hubble ($H_0$) derivada del Fondo Cósmico de Microondas (CMB) bajo el modelo estándar $\Lambda$CDM ($H_0^{CMB} \approx 67.4$ km/s/Mpc) y las mediciones locales de la escala de distancia ($H_0^{local} \approx 73.5$ km/s/Mpc).

Las soluciones a esta tensión que involucran nueva física antes de la recombinación generalmente predicen un Universo más joven ($A_U \approx 12.9 \pm 0.2$ Gyr) para acomodar el mayor $H_0 manteniendo la consistencia con los datos del CMB. Por el contrario, el modelo$\Lambda$CDM estándar calibrado por el CMB predice $A_U \approx 13.8$ Gyr. Los autores pretenden probar estas predicciones contrapuestas derivando $A_\star$ a partir de una muestra grande y de alta calidad de estrellas galácticas sin imponer fuertes priors cosmológicos que sesguen el resultado hacia cualquiera de los dos modelos.

Metodología
Los autores utilizan una muestra de 247,103 estrellas subgigantes del catálogo de Xiang & Rix (2022, en adelante X22), que combina espectroscopía de LAMOST DR7 con astrometría de Gaia eDR3. Las edades estelares se estimaron inicialmente utilizando isócronas de Yonsei-Yale (YY). Para asegurar la fiabilidad de las estimaciones de edad para la inferencia cosmológica, los autores aplicaron una serie de rigurosos cortes de calidad a nivel de población, diseñados para eliminar estrellas con determinaciones de edad problemáticas (por ejemplo, aquellas con incertidumbres irrealmente bajas a edades altas) sin introducir sesgos cosmológicos circulares.

Los pasos clave de control de calidad incluyeron:

1. Cortes astrométricos: Requerir una incertidumbre de paralaje $\sigma_\varpi/\varpi < 0.1$ y $\varpi > 0.2$ mas (limitando la muestra a dentro de 5 kpc).
2. Contaminación por binarias: Excluir estrellas con valores altos de ipd_frac_multi_peak en los datos de Gaia para eliminar binarias marginalmente resueltas.
3. Consistencia química: Aplicar cortes dependientes de la edad en la metalicidad ([Fe/H]) y la abundancia de elementos $\alpha$ ([$\alpha$/Fe]). Los autores exigen que las estrellas viejas sean pobres en metales y ricas en $\alpha$, de acuerdo con los modelos de evolución química galáctica. Las estrellas que parecen viejas pero son ricas en metales o pobres en $\alpha$ fueron rechazadas como probables errores de medición.
4. Validación entre catálogos: Comparar las edades de YY con las edades de FLAME derivadas únicamente de los datos de Gaia DR3. Se ajustó una línea de tendencia lineal a la relación edad-masa para imputar las edades de FLAME faltantes, seguido de un rechazo de valores atípicos de $3\sigma$ entre las edades de YY y FLAME.

Tras estos cortes, la muestra nominal final contenía 155,600 estrellas. Los autores emplearon dos enfoques estadísticos primarios para reconstruir la distribución de la edad latente $P(A)$ e inferir $A_\star$:

• Reconstrucción iterativa: Un método no paramétrico donde la distribución de la población sirve como un prior para refinar los posteriors individuales, iterado hasta la convergencia.
• Reconstrucción MCMC: Un enfoque de Monte Carlo por Cadenas de Markov para inferir $P(A)$ en 100 bins de edad (0.1–20.1 Gyr). Los autores modelaron el log-posterior de la distribución, incluyendo una penalización de gradiente para suavizar el ruido. Para inferir $A_\star$, analizaron la relación de la densidad de probabilidad reconstruida con su incertidumbre ($P/\sigma$), identificando un punto de "cruce" agudo donde la distribución transiciona de una tendencia decreciente a un "suelo de sesgo" plano (donde el algoritmo fuerza $P > 0$).

Resultados Clave
El análisis arroja una estimación robusta para la edad de la estrella más antigua de la muestra:
$$A_\star = 13.73^{+0.18}_{-0.15} \text{ Gyr}$$
Este resultado se deriva de la muestra nominal y el método de análisis. Los autores probaron la robustez de este valor variando los cortes de calidad:

• Cortes de metalicidad: Relajar el techo de [Fe/H] (permitiendo más estrellas ricas en metales) aumentó $A_\star$ ligeramente (hasta $14.02^{+0.18}_{-0.15}$ Gyr), mientras que endurecerlo la disminuyó (hasta $13.31^{+0.21}_{-0.18}$ Gyr). Sin embargo, los autores argumentan que endurecer el corte en 0.2 dex es acientífico, ya que elimina estrellas consistentes con la tendencia edad-metalicidad galáctica y excluye poblaciones conocidas como la fusión Gaia-Sausage/Enceladus.
• Enriquecimiento de $\alpha$ y cortes de validación entre catálogos: Las variaciones en el suelo de [$\alpha$/Fe] y la rigurosidad del acuerdo YY-FLAME ($2.5\sigma$ vs $3\sigma$) resultaron en valores de $A_\star$ consistentes con el resultado nominal dentro de las incertidumbres.

Un análisis secundario centrado en una submuestra de "química temprana" (1,556 estrellas con bajo [Fe/H] y alto [$\alpha$/Fe]) utilizando un modelo paramétrico de $\tau$-retrasado para $P(A)$ produjo un resultado consistente de $A_\star = 13.85^{+0.05}_{-0.04}$ Gyr, apoyando aún más el hallazgo principal.

Significancia y Reivindicaciones
El artículo afirma que su $A_\star$ inferido es consistente con la predicción estándar de $\Lambda$CDM calibrada por el CMB ($A_\star^{CMB} \approx 13.6$ Gyr, asumiendo $t_f = 0.2$ Gyr) pero en tensión significativa ($>2\sigma$) con las predicciones de las soluciones de "tiempo temprano" para la tensión de Hubble ($A_\star^{ETS} \approx 12.7$ Gyr).

Los autores argumentan que:

1. Restricciones sobre la nueva física: El acuerdo con la edad de $\Lambda$CDM estándar pone en duda las soluciones a la tensión de Hubble que dependen únicamente de nueva física antes de la recombinación, ya que estas generalmente requieren una edad cósmica de $12.9 \pm 0.2$ Gyr. Reconciliar una edad tan baja con el resultado de los autores requeriría que las incertidumbres en el modelado estelar fueran significativamente mayores de lo que se entiende actualmente, lo cual es improbable dada la baja metalicidad de la muestra y las restricciones asterosísmicas independientes.
2. Robustez Metodológica: El estudio demuestra que los cortes de calidad a nivel de población y de gran muestra (basados en la química y la consistencia entre catálogos) son esenciales para eliminar valores atípicos que de otro modo sesgarían las estimaciones de edad. El uso de la reconstrucción MCMC permite un manejo riguroso de la distribución de la edad latente y sus incertidumbres.
3. Consistencia con trabajos previos: El resultado se alinea con estimaciones previas de los cúmulos globulares más antiguos y estrellas individuales (por ejemplo, HD 140283), reforzando la visión de que el Universo tiene aproximadamente 13.8 mil millones de años.

El artículo concluye que, si bien existen incertidumbres en el modelado estelar, estas son insuficientes para reconciliar la alta edad inferida de las estrellas más antiguas con las bajas edades predichas por las soluciones de tiempo temprano de la tensión de Hubble, apoyando así la edad estándar del Universo según el modelo $\Lambda$CDM.