A Test of Substellar Evolutionary Models with High-Precision Ages from Asteroseismology and Gyrochronology for the Benchmark System HR 7672AB

Título: Una prueba de modelos evolutivos subestelares con edades de alta precisión obtenidas mediante asterosismología y girocronología para el sistema de referencia HR 7672AB

1. El Problema

Los enanas marrones son objetos con masas entre los planetas gigantes gaseosos y las estrellas de menor masa (aprox. 13–75 $M_{Jup}$). A diferencia de las estrellas, no pueden sostener la fusión estable de hidrógeno y, por lo tanto, se enfrían y desvanecen con el tiempo. Los modelos evolutivos predicen su luminosidad, radio y temperatura efectiva en función de la masa y la edad. Sin embargo, validar estos modelos es extremadamente difícil porque:

• La mayoría de las enanas marrones son objetos flotantes libres, por lo que su edad y masa no son observables directamente.
• Los sistemas de referencia ("benchmarks") que tienen masa, edad y luminosidad medidas simultáneamente son escasos.
• Las mediciones de edad en sistemas existentes a menudo tienen incertidumbres grandes, lo que impide distinguir entre diferentes predicciones de modelos, especialmente cerca del límite estelar/subestelar.

El sistema HR 7672AB (una estrella tipo G0V y una compañera enana L) es un candidato ideal, pero carecía de una edad precisa y una masa dinámica refinada para realizar pruebas rigurosas.

2. Metodología

Los autores combinaron múltiples técnicas de observación de alta precisión para caracterizar el sistema:

• Asterosismología (Keck Planet Finder - KPF): Se obtuvieron tres noches de datos de velocidad radial (RV) de ultra alta precisión (~70 cm/s) utilizando el instrumento KPF en el telescopio Keck I. Se detectaron oscilaciones de 5 minutos en la estrella primaria (HR 7672A). Se aplicó deconvolución de Gold y ajuste de modelos de ondas senoidales para extraer frecuencias de modos acústicos ($\nu_{n,\ell}$).
• Girocronología: Se determinó el período de rotación de la estrella primaria mediante el análisis de series temporales del índice de actividad S (Ca II H&K) de múltiples instrumentos (Mount Wilson, Lick/APF, Keck/HIRES) y se modeló utilizando códigos de pérdida de momento angular (rotevol) y relaciones empíricas.
• Órbita Dinámica: Se combinaron nuevas mediciones de astrometría directa (imagen con NIRC2/Keck II) con datos históricos de imágenes y una larga serie temporal de velocidades radiales (incluyendo datos de California Planet Search y Gaia/Hipparcos). Se utilizó el código orvara (MCMC) para ajustar la órbita binaria y derivar las masas dinámicas.
• Fotometría TESS: Se analizaron datos de TESS casi simultáneos con las observaciones de KPF para detectar señales sísmicas débiles y comparar la relación de amplitud fotometría-RV con el Sol.
• Pruebas de Modelos: Se compararon las mediciones observacionales (masa, edad, luminosidad) con seis conjuntos de modelos de enfriamiento de enanas marrones: Burrows1997, SM2008, BHAC2015, ATMO2020, CBPD2023 y Diamondback. Se realizaron pruebas de consistencia (rejection sampling) y pruebas de bondad de ajuste (distancia de Mahalanobis).

3. Contribuciones Clave

• Edad de Alta Precisión: Se obtuvo una edad estelar combinada de 2.26 ± 0.40 Gyr (incertidumbre del 18%), derivada de la combinación de la edad sísmica (1.87 ± 0.65 Gyr) y la edad por girocronología (2.58 ± 0.47 Gyr). Esta es una de las edades más precisas para una estrella de tipo solar joven.
• Masas Dinámicas Refinadas: Se determinó una masa dinámica para la estrella primaria de 1.111 ± 0.017 $M_{\odot}$ y para la compañera subestelar de 75.39 ± 0.67 $M_{Jup}$. La incertidumbre en la masa de la compañera es del 0.9%.
• Sistema de Referencia Excepcional: La combinación de una incertidumbre de masa del 0.9% y una edad del 18% convierte a HR 7672B en un "benchmark" sin precedentes para probar la física de enfriamiento cerca del límite de fusión de hidrógeno.
• Validación de Amplitudes Sísmicas: Se detectó una señal sísmica débil en TESS y se confirmó que la relación de amplitud fotometría-RV es consistente con los valores solares y las simulaciones hidrodinámicas 3D.

4. Resultados Principales

• Comparación con Modelos:
  • El modelo CBPD2023 (Chabrier et al. 2023), que incorpora una nueva ecuación de estado (EOS) para el interior, mostró el mejor acuerdo con las observaciones, concordando dentro de < 0.3$\sigma$ en todas las pruebas (masa, edad, luminosidad).
  • Los modelos Diamondback, SM2008 y BHAC2015 (que incluyen tratamiento de nubes) predicen luminosidades ligeramente más altas que las observadas para la masa medida. Esto sugiere que, según estos modelos, el objeto podría estar en el límite inferior de la secuencia principal (estelar) en lugar de ser subestelar, o que su límite estelar/subestelar predicho es demasiado bajo en masa.
  • El modelo Burrows1997 predice las luminosidades más bajas y confirma claramente la naturaleza subestelar, pero con un desajuste mayor que CBPD2023.
• Límite Estelar/Subestelar: La masa de HR 7672B (75.39 $M_{Jup}$) cae justo en la zona donde los modelos difieren sobre si el objeto alcanzará la fusión de hidrógeno. Los modelos CBPD2023 y ATMO2020 predicen un límite más alto (>77 $M_{Jup}$), mientras que Diamondback y SM2008 predicen un límite más bajo (~69-72 $M_{Jup}$), lo que implicaría que HR 7672B es una estrella de muy baja masa en estos últimos modelos.
• Inconsistencias: Tres de los modelos (Diamondback, SM2008, BHAC2015) luchan por predecir una luminosidad lo suficientemente baja para HR 7672B dada su masa dinámica, independientemente de la edad, lo que indica una posible deficiencia en la física de enfriamiento o en la ecuación de estado utilizada en esos modelos.

5. Significado

Este trabajo representa un avance significativo en la astrofísica estelar y subestelar:

• Validación de Física Interna: La superioridad del modelo CBPD2023 respalda la importancia de las nuevas ecuaciones de estado (EOS) para el hidrógeno y helio en el interior de objetos de baja masa.
• Límite Crítico: HR 7672B se sitúa en el "desierto de enanas marrones" y en el límite estelar/subestelar. Su caracterización precisa permite descartar modelos que no pueden reproducir el enfriamiento correcto en este régimen de masa.
• Futuro: El sistema servirá como referencia fundamental para futuras misiones. Mejoras en la luminosidad (con JWST) y en la masa (con Gaia DR4) permitirán discriminar aún más entre los modelos, especialmente para objetos más viejos donde las diferencias en las predicciones de enfriamiento son más pronunciadas.

En resumen, el estudio demuestra que la combinación de asterosismología de alta precisión y dinámica orbital puede proporcionar pruebas rigurosas de los modelos evolutivos, revelando que las actualizaciones recientes en la ecuación de estado son cruciales para describir correctamente el enfriamiento de las enanas marrones masivas.