Detailed Investigation of a W UMa Contact Binary with an Ultralow Mass Ratio and a Third-Body as a Potential Merger Candidate

Resumo Técnico: Investigação Detalhada de um Binário de Contato W UMa com uma Razão de Massa Ultrabaixa e um Terceiro Corpo como um Potencial Candidato à Fusão

Problema
O limite inferior da razão de massa ($q = M_2/M_1$) em sistemas binários de contato W UMa permanece um objeto de incerteza, com evidências observacionais sugerindo a existência de sistemas abaixo das previsões de estabilidade teórica. A estabilidade dinâmica de tais sistemas de razão de massa extremamente baixa ($q < 0,1$) é debatida, pois são potenciais precursores de fusões estelares. Estudos anteriores de sistemas como VSX J082700.8+462850 e TYC 4002-2628-1 basearam-se fortemente em modelagem fotométrica, muitas vezes carecendo de confirmação espectroscópica do componente secundário devido à sua intrínseca baixa luminosidade. Um estudo recente por Li et al. (2024) analisou o alvo TYC 3801-1529-1, relatando uma razão de massa de $q \approx 0,0356$. No entanto, os autores deste estudo observaram discrepâncias no ajuste da curva de luz e buscaram reavaliar os parâmetros do sistema, as variações do período orbital e o status de estabilidade usando dados do TESS e técnicas de modelagem avançadas.

Metodologia
Os autores realizaram uma reanálise abrangente de TYC 3801-1529-1 (TIC 445065429) utilizando as seguintes abordagens:

1. Modelagem Fotométrica: Os dados da curva de luz do TESS dos setores 20, 47, 60 e 74 foram processados usando o Quick-Look Pipeline (QLP). A análise empregou o código PHOEBE 2.4.9 para validação inicial e a aplicação BSN (um acelerador MCMC baseado em Windows) para a otimização final. A aplicação BSN, configurada com 28 walkers e 1200 iterações, foi utilizada para derivar as incertezas e distribuições dos parâmetros, amostrando cinco parâmetros principais: temperaturas efetivas ($T_1, T_2$), razão de massa ($q$), fator de preenchimento (fillout factor, $f$) e inclinação orbital ($i$).
2. Estimativa de Parâmetros Absolutos: Utilizando a razão de massa derivada e assumindo a massa primária consistente com Li et al. (2024), os autores calcularam a massa secundária, o semieixo maior (via terceira lei de Kepler), os raios estelares e as luminosidades.
3. Análise de Variação do Período Orbital (O-C): Os autores compilaram 408 tempos de mínimo do TESS, um do ZTF e 17 de observações terrestres (incluindo novos dados do Observatório de San Pedro Mártir), cobrindo um período de seis anos. Eles aplicaram o método Observado menos Calculado (O-C) para investigar mudanças de período.
4. Análise de Terceiro Corpo e Aliasing: Para distinguir entre mudanças de período intrínsecas e influências externas, os autores modelaram os resíduos O-C usando uma função senoidal para testar o efeito de tempo de viagem de luz de um terceiro corpo. Eles também realizaram uma análise de periodograma de Lomb-Scargle e testes de prewhitening para descartar efeitos de aliasing causados por padrões de amostragem sazonal.
5. Análise de Estabilidade e de Derivada: O estudo calculou a razão entre o momento angular de rotação e o momento angular orbital ($J_{spin}/J_{orb}$) usando raios de griação ($k^2$) atualizados e apropriados para as estruturas estelares específicas. Além disso, o método da derivada de Kouzuma (2023, 2025) foi aplicado à curva de luz para estimar independentemente a razão de massa.

Resultados Principais

• Razão de Massa Ultrabaixa: A reanálise resultou em uma razão de massa de $q = 0,024^{+0,002}_{-0,001}$, significativamente inferior aos $0,0356$ relatados por Li et al. (2024). Esta representa a menor razão de massa determinada experimentalmente para um binário de contato até o momento.
• Propriedades dos Componentes: O componente secundário é estimado ter uma massa de $M_2 \approx 0,050 M_\odot$, posicionando-o dentro do regime de anãs marrons ($0,013–0,08 M_\odot$). A massa primária é $M_1 \approx 2,096 M_\odot$. O sistema é classificado como um binário de contato do tipo A com um fator de preenchimento médio ($f \approx 0,315$).
• Presença de Terceiro Corpo: O diagrama O-C exibe uma clara variação cíclica. Após a remoção de uma tendência linear, os resíduos foram melhor ajustados por um modelo senoidal com um período de $P_{mod} \approx 8,75$ anos. A função de massa derivada e simulações de Monte Carlo sugerem um terceiro corpo com uma massa de $M_3 \approx 0,49 M_\odot$ (provavelmente uma anã M de baixa massa) contribuindo com luminosidade desprezível ($l_3 \approx 0,001$).
• Variações de Período: Uma vez removida a contribuição do terceiro corpo, as variações de período intrínsecas do binário puderam ser descritas por uma tendência linear com inclinação negativa ou uma tendência quadrática com curvatura descendente. Métricas estatísticas (AIC/BIC) favoreceram o modelo linear mais simples, embora os autores notem que um ajuste quadrático implicaria em um período orbital decrescente.
• Instabilidade Dinâmica: A razão calculada entre o momento angular de rotação e o orbital é $J_{spin}/J_{orb} \approx 1,18$. Este valor excede significativamente o limiar de estabilidade de Darwin ($\sim 0,33$), indicando que o sistema é dinamicamente instável.
• Validação da Razão de Massa: A razão de massa derivada via método da derivada de Kouzuma ($q \approx 0,031$) foi encontrada consistente com os resultados de MCMC e com Li et al. (2024), reforçando a conclusão de que o sistema possui uma razão de massa extremamente baixa.

Significância e Alegações
O artigo afirma que o TYC 3801-1529-1 representa um candidato notável para uma fusão binária devido à sua configuração dinamicamente instável e razão de massa ultrabaixa. O estudo destaca os desafios na detecção de secundárias tênues em tais sistemas, que frequentemente levam a degenerações nas soluções fotométricas. Ao demonstrar que as variações de período do sistema são provavelmente impulsionadas por um companheiro terciário, em vez de atividade magnética (o mecanismo de Applegate foi descartado devido à insuficiência de variação do momento quadrupolar), os autores fornecem um quadro mais claro da arquitetura do sistema.

Os autores enfatizam que, embora o sistema seja um forte candidato à fusão, a direção definitiva da mudança do período orbital (aumento vs. diminuição) permanece ambígua sem observações de longo prazo adicionais. O estudo ressalta a necessidade de acompanhamento espectroscópico de alta resolução para confirmar a razão de massa e a natureza da secundária, que é atualmente inferida como sendo uma anã marrom. Os achados contribuem para o debate contínuo sobre a razão de massa mínima para binários de contato e as vias evolutivas que levam a fusões estelares.