A Bayesian Search for Planet Engulfment Signatures in Solar Analogs

Resumen Técnico: Una Búsqueda Bayesiana de Firmas de Engullimiento Planetario en Análogos Solares

Planteamiento del Problema
Los planetas pueden ser impulsados hacia sus estrellas anfitrionas mediante el decaimiento de marea o la inestabilidad dinámica, lo que potencialmente deja huellas químicas mensurables en la fotosfera estelar si el material acrecentado se mezcla en una envoltura superficial. Si bien las búsquedas previas basadas en abundancias han identificado con éxito tales firmas en sistemas binarios co-natales (donde una estrella sirve como control químicamente emparejado), sigue sin estar claro si es posible identificar firmas comparables en la población mucho más amplia de estrellas de campo no co-natales. Los métodos existentes para estrellas individuales suelen basarse en la correlación entre la abundancia elemental y la temperatura de condensación (pendiente de $T_c$) o en las relaciones entre volátiles y refractarios. Sin embargo, estos enfoques comprimen patrones multielementales en valores únicos, lo que puede oscurecer estructuras detalladas y no lograr distinguir el engullimiento planetario de otros procesos astrofísicos como la evolución química galáctica (GCE), la difusión atómica o factores estocásticos de formación.

Metodología
Los autores presentan un marco bayesiano sistemático para modelar patrones de abundancia completos en 113 gemelos solares y análogos, 45 de los cuales albergan planetas o enanas marrones conocidos o candidatos. El conjunto de datos comprende mediciones de abundancia diferencial de alta precisión, línea por línea, para 19 elementos (de C a Zn) provenientes de dos fuentes: la muestra de Bedell et al. (2018) y el sondeo Planets Around Solar Twins/Analogs (PASTA).

La metodología central consiste en comparar tres clases de modelos para el vector de abundancia de cada estrella (relativo al Sol o a la media de la muestra):

1. Modelos Nulos:
   • Modelo Plano ($F$): Un desplazamiento independiente de los elementos más dispersión intrínseca, que representa la hipótesis de ausencia de estructura resuelta.
   • Evolución Química Galáctica (GCE): Se probaron dos variantes: $G_{Bedell}$ (relaciones [X/Fe]-edad lineales) y $G_{GALAH}$ (un modelo de dos procesos que combina enriquecimiento rápido y retardado).
2. Modelo de Engullimiento ($E$): Supone material acrecentado bien mezclado en la zona de convección estelar. El modelo predice aumentos de abundancia basados en la masa engullida ($M_p$) y la composición del material acrecentado, probado contra dos composiciones: Bulk Earth (Tierra total) y condrito CM.

El análisis utiliza una verosimilitud gaussiana que incluye las incertidumbres reportadas y un parámetro de dispersión intrínseca ($\sigma_{int}$). La evidencia bayesiana ($Z$) se calcula utilizando el algoritmo de muestreo anidado dynesty. La preferencia del modelo se cuantifica mediante el logaritmo del factor de Bayes ($\Delta \ln Z$). Para establecer umbrales de detección robustos y evitar cortes fijos, los autores realizaron un bootstrapping paramétrico sobre 121 vectores de abundancia simulados para determinar el $\Delta \ln Z$ máximo esperado bajo hipótesis nulas, correspondiente a una probabilidad de falsos positivos de $\sim 1/1210$.

El análisis fiduciario adopta al Sol como línea de base y utiliza relaciones [X/Fe] para mitigar las firmas de difusión atómica. Una "prueba de espejo" (reversión de las diferencias de abundancia) confirmó que la configuración Sol-[X/Fe] es la única configuración donde revertir los datos reduce significativamente los factores de Bayes para el engullimiento, validando al Sol como una línea de base físicamente motivada y sin historia de engullimiento reciente.

Resultos Clave
Bajo la configuración fiduciaria Sol-[X/Fe], el estudio identificó tres estrellas candidatas de alta significación cuyas firmas de abundancia favorecen fuertemente el modelo de engullimiento planetario sobre las alternativas de GCE y el modelo plano:

• TOI-3342: Favorece la composición Bulk Earth; masa engullida inferida de $16 \pm 4 \, M_{\oplus}$.
• HIP 101905: Favorece la composición Bulk Earth; masa engullida inferida de $7.5 \pm 1.0 \, M_{\oplus}$.
• TOI-2426: Favorece la composición condrito CM; masa engullida inferida de $33^{+9}_{-7} \, M_{\oplus}$.

Además, HIP 30502 emergió como candidato bajo la configuración de la media-[X/Fe]. El estudio también identificó 14 candidatos "optimistas" que superan al modelo plano pero que aún podrían explicarse por la GCE.

La tasa de detección bruta de firmas de engullimiento planetario en esta muestra es de aproximadamente el 2.7% (3/113). Al separarse por muestra, la tasa es del 1.3% para la muestra genérica de Bedell y mayor para la muestra de anfitriones planetarios de PASTA, aunque los autores advierten contra la sobreinterpretación de esta diferencia debido a la función de selección y las variaciones sistemáticas entre los conjuntos de datos.

Significancia y Reivindicaciones
El artículo afirma extender las búsquedas de engullimiento basadas en abundancias más allá del contexto de estrellas binarias, proporcionando un marco para sondear la coevolución estrella-planeta en análogos solares individuales. Las contribuciones clave son:

• Avance Metodológico: Demostrar que el modelado bayesiano basado en plantillas de patrones de abundancia completos es superior al análisis de la pendiente de $T_c$. Los autores muestran que las pendientes grandes de $T_c$ no son condiciones necesarias ni suficientes para coincidir con una plantilla de engullimiento, y que los modelos de GCE pueden imitar tendencias que serían malinterpretadas por análisis de solo pendiente.
• Tasa de Detección: Establecer una tasa de detección nominal de $\sim 1–3\%$ de firmas de engullimiento planetario entre análogos solares, lo cual es consistente con los límites superiores derivados de estudios binarios, pero derivado aquí para la población de campo.
• Identificación de Candidatos: Proporcionar candidatos específicos de alta confianza (TOI-3342, HIP 101905, TOI-2426), con masas engullidas inferidas que oscilan entre $\sim 7.5$ y $33 \, M_{\oplus}$, todas muy por debajo de los sólidos totales en una nebulosa protosolar de masa mínima.

Los autores mantienen una postura modesta, reconociendo que la tasa de detección es una estimación bruta sin correcciones de sensibilidad y que otros efectos (por ejemplo, procesos de formación planetaria, GCE, difusión atómica) pueden introducir firmas degeneradas. Enfatizan que el trabajo futuro requiere muestras más grandes y homogéneas con arquitecturas planetarias bien caracterizadas para constreñir precisamente las tasas de engullimiento y distinguir entre las competencias de los orígenes físicos.