A Bayesian Search for Planet Engulfment Signatures in Solar Analogs

Resumo Técnico: Uma Busca Bayesiana por Assinaturas de Engolfamento Planetário em Análogos Solares

Definição do Problema
Planetas podem ser levados para dentro de suas estrelas hospedeiras via decaimento de maré ou instabilidade dinâmica, potencialmente deixando impressões químicas mensuráveis na fotosfera estelar se o material acolhido se misturar em um envelope raso. Embora buscas baseadas em abundâncias anteriores tenham identificado com sucesso tais assinaturas em sistemas binários co-natais (onde uma estrela serve como um controle quimicamente pareado), permanece incerto se assinaturas comparáveis podem ser identificadas na população muito maior de estrelas de campo não co-natais. Métodos existentes para estrelas únicas frequentemente dependem da correlação entre a abundância elemental e a temperatura de condensação (inclinação $T_c$) ou razões entre voláteis e refratários. No entanto, essas abordagens comprimem padrões multi-elementares em números únicos, o que pode obscurecer estruturas detalhadas e falhar em distinguir o engolfamento planetário de outros processos astrofísicos, como a Evolução Química Galáctica (GCE), difusão atômica ou fatores estocásticos de formação.

Metodologia
Os autores apresentam uma estrutura Bayesiana sistemática para modelar padrões de abundância completos em 113 gêmeos solares e análogos, dos quais 45 abrigam planetas ou anãs marrons conhecidos ou candidatos. O conjunto de dados compreende medições de abundância diferencial linha a linha de alta precisão para 19 elementos (C até Zn) de duas fontes: a amostra de Bedell et al. (2018) e o levantamento Planets Around Solar Twins/Analogs (PASTA).

A metodologia central envolve a comparação de três classes de modelos para o vetor de abundância de cada estrela (relativo ao Sol ou à média da amostra):

1. Modelos Nulos:
   • Modelo Plano ($F$): Um deslocamento independente de elemento mais dispersão intrínseca, representando a hipótesia de ausência de estrutura resolvida.
   • Evolução Química Galáctica (GCE): Duas variantes foram testadas: $G_{Bedell}$ (relações [X/Fe]-idade lineares) e $G_{GALAH}$ (um modelo de dois processos combinando enriquecimento imediato e tardio).
2. Modelo de Engolfamento ($E$): Assume material acolhido bem misturado na zona de convecção estelar. O modelo prevê aumentos de abundância baseados na massa engolfada ($M_p$) e na composição do material acolhido, testado contra duas composições: Bulk Earth (Terra Primitiva) e condrito CM.

A análise utiliza uma verossimilhança Gaussiana incluindo incertezas relatadas e um parâmetro de dispersão intrínseca ($\sigma_{int}$). A evidência Bayesiana ($Z$) é calculada usando o algoritmo de amostragem aninhada dynesty. A preferência do modelo é quantificada via o logaritmo do fator de Bayes ($\Delta \ln Z$). Para estabelecer limiares de detecção robustos e evitar cortes fixos, os autores realizaram bootstrapping paramétrico em 121 vetores de abundância fictícios para determinar o $\Delta \ln Z$ máximo esperado sob hipóteses nulas, correspondendo a uma probabilidade de falso-positivo de $\sim 1/1210$.

A análise fiduciária adota o Sol como linha de base e utiliza razões [X/Fe] para mitigar assinaturas de difusão atômica. Um "teste espelho" (revertendo as diferenças de abundância) confirmou que a configuração Sol-[X/Fe] é a única configuração onde reverter os dados reduz significativamente os fatores de Bayes para o engolfamento, validando o Sol como uma linha de base fisicamente motivada e sem histórico recente de engolfamento.

Principais Resultados
Sob a configuração fiduciária Sol-[X/Fe], o estudo identificou três estrelas candidatas de alta significância cujos padrões de abundância favorecem fortemente o modelo de engolfamento planetário sobre as alternativas de GCE e o modelo plano:

• TOI-3342: Favorece a composição Bulk Earth; massa engolfada inferida de $16 \pm 4 \, M_{\oplus}$.
• HIP 101905: Favorece a composição Bulk Earth; massa engolfada inferida de $7.5 \pm 1.0 \, M_{\oplus}$.
• TOI-2426: Favorece a composição condrito CM; massa engolfada inferida de $33^{+9}_{-7} \, M_{\oplus}$.

Adicionalmente, HIP 30502 emergiu como um candidato sob a configuração média-[X/Fe]. O estudo também identificou 14 candidatos "otimistas" que superam o modelo plano, mas que ainda podem ser explicados por GCE.

A taxa de detecção bruta de assinaturas de engolfamento planetário nesta amostra é de aproximadamente 2,7% (3/113). Quando separada por amostra, a taxa é de 1,3% para a amostra genérica de Bedell e maior para a amostra de hospedeiros planetários PASTA, embora os autores alertem contra a interpretação excessiva dessa diferença devido à função de seleção e variações sistemáticas entre os conjuntos de dados.

Significância e Alegações
O artigo afirma estender as buscas de abundância por engolfamento além do contexto de sistemas binários, fornecendo um arcabouço para investigar a coevolução estrela-planeta em análogos solares individuais. As principais contribuições são:

• Avanço Metodológico: Demonstrar que a modelagem Bayesiana baseada em modelos de template de padrões de abundância completos é superior à análise de inclinação $T_c$. Os autores mostram que grandes inclinações $T_c$ não são condições necessárias nem suficientes para corresponder a um template de engolfamento, e que modelos de GCE podem mimetizar tendências que seriam interpretadas incorretamente por análises de apenas inclinação.
• Taxa de Detecção: Estabelecer uma taxa de detecção nominal de $\sim 1–3\%$ de assinaturas de engolfamento planetário entre análogos solares, o que é consistente com limites superiores derivados de estudos binários, mas derivado aqui para a população de campo.
• Identificação de Candidatos: Fornecer candidatos específicos de alta confiança (TOI-3342, HIP 101905, TOI-2426), com massas engolfadas inferidas variando de $\sim 7.5$ a $33 \, M_{\oplus}$, todas bem abaixo do total de sólidos em uma nebulosa de massa mínima solar.

Os autores mantêm uma postura modesta, reconhecendo que a taxa de detecção é uma estimativa bruta sem correções de sensibilidade e que outros efeitos (ex: processos de formação planetária, GCE, difusão atômica) podem introduzir assinaturas degeneradas. Eles enfatizam que trabalhos futuros exigem amostras maiores e homogêneas com arquiteturas planetárias bem caracterizadas para restringir precisamente as taxas de engolfamento e distinguir entre as competidoras origens físicas.