The quest for Magrathea planets. II. Orbital stability of exoplanets formed around double white dwarfs

Resumo Técnico: Estabilidade Orbital de Exoplanetas Formados em Torno de Anãs Brancas Duplas

Enunciado do Problema
Embora modelos teóricos sugiram que planetas de "segunda geração" podem se formar a partir de discos circumbinários criados pelo material ejetado durante a fase de envelope comum de binárias de anãs brancas duplas (DWD), a estabilidade dinâmica de longo prazo de tais sistemas permanece não testada. Estudos anteriores demonstraram o potencial de formação de planetas gigantes nesses ambientes (Ledda et al. 2023), contudo, não está claro se esses sistemas podem sobreviver às interações gravitacionais caóticas inerentes às arquiteturas circumbinárias de múltiplos planetas ao longo de escalas de tempo de milhões de anos. Além disso, a estabilidade desses sistemas é crítica para avaliar sua detectabilidade pela futura missão LISA (Laser Interferometer Space Antenna), que visa detectar exoplanetas orbitando DWDs via modulações de frequência de ondas gravitacionais (GW).

Metodologia
Os autores realizaram simulações de $N$-corpos de sistemas planetários circumbinários orbitando três configurações específicas de DWD (DWD2, DWD3 e DWD4) extraídas da população de DWD do LISA. O estudo empregou um integrador simplético híbrido (Chambers et al. 2002) especificamente adaptado para sistemas binários para lidar tanto com a evolução Kepleriana de longo prazo quanto com encontros próximos/colisões.

As simulações cobriram quatro populações distintas baseadas em condições iniciais:

• População A: Sistemas não ressonantes com planetas de menor massa ($M_p \leq 1.2 M_J$) em torno de binárias de menor massa ($M_{bin} < 0.6 M_\odot$).
• População B: Sistemas não ressonantes com planetas de maior massa ($M_p \geq 2.5 M_J$) em torno de binárias mais massivas ($M_{bin} \sim 1 M_\odot$).
• População C: Contrapartes ressonantes da População B, onde os planetas estão inicialmente travados em ressonâncias de movimento médio 2:1.
• População D: Sistemas de cinco planetas em uma ressonância específica de "NMS" (uma combinação de ressonâncias 4:3 e 5:3) em torno de DWD*4.

As condições iniciais incluíram massas planetárias extraídas de distribuições baseadas em Ledda et al. (2023), semieixos maiores entre o limite de estabilidade crítica ($a_{crit}$) e 1 au (ou escalonados para casos ressonantes), e excentricidades iniciais baixas ($e < 0.01$) e inclinações ($1^\circ-3^\circ$). As simulações rodaram por $3 \times 10^6$ a $10^7$ anos.

Para analisar os resultados, os autores utilizaram diversas métricas:

1. Déficit de Momento Angular Normalizado (NAMD): Para quantificar a excitação dinâmica (temperatura) do sistema.
2. Rastreamento de Eventos Catastróficos: Registrando ejeções (dispersão além de 100 au), colisões e aproximações próximas à binária.
3. Espaçamento Orbital e Concentração de Massa: Analisando mudanças no espaçamento orbital ($S_s$) e na fração de massa do planeta mais massivo ($S_m$).
4. Deslocamento do Centro de Massa (CoM): Calculando o deslocamento do baricentro do sistema planetário para inferir quais planetas foram perdidos.

Principais Resultados

• Estabilidade e Multiplicidade: Sistemas planetários orbitando DWDs podem permanecer estáveis ao longo dos períodos estudados, mas a estabilidade depende fortemente da multiplicidade inicial. Sistemas começando com dois planetas são os mais propensos a sobreviver inalterados. Por outro vez, sistemas com três, quatro ou cinco planetas frequentemente experimentam instabilidade catastrófica.
• Resultados das Populações:
  • População A (Baixa massa, não ressonante): 82,5% mantiveram sua contagem inicial de planetas; 17,5% sofreram eventos catastróficos.
  • População B (Alta massa, não ressonante): Apenas 49% mantiveram sua multiplicidade inicial; 51% sofreram perda de planetas.
  • População C (Alta massa, ressonante): 65,6% mantiveram a multiplicidade inicial. Embora a ressonância tenha aumentado a estabilidade para sistemas de baixa multiplicidade, sistemas ressonantes de quatro planetas foram mais propensos à disrupção do que seus contrapartes não ressonantes.
  • População D (Cinco planetas ressonantes): 100% dos sistemas sofreram eventos catastróficos. A maioria (58%) perdeu exatamente dois planetas, e 17,2% foram completamente disruptados (todos os planetas perdidos).
• Evolução Dinâmica: Sistemas que passaram por instabilidade tipicamente evoluíram para estados "dinamicamente quentes" (NAMD > $1.3 \times 10^{-3}$), caracterizados por maior espaçamento orbital e maior concentração de massa no planeta sobrevivente mais massivo. Uma forte anticorrelação foi observada entre a multiplicidade final e o NAMD: sistemas terminando com menos planetas geralmente exibiam maior excitação dinâmica.
• Efeitos de Ressonância: Embora as ressonâncias 2:1 iniciais tenham ajudado a preservar sistemas de baixa multiplicidade (2-3 planetas), elas não impediram a disrupção de cadeias de alta multiplicidade (4-5 planetas), que quase universalmente quebraram a ressonância e perderam planetas.
• Detectabilidade pelo LISA:
  • Sistemas de Planeta Único: Aproximadamente 7% dos sistemas totais simulados terminaram como sistemas de planeta único. A maioria dos sobreviventes possui massas e períodos que caem na janela de detecção do LISA, desde que haja localização de céu e distância favoráveis.
  • Sistemas de Múltiplos Planetas: Os autores derivaram uma fórmula (Eq. 15) para estimar a amplitude global do deslocamento de frequência de onda gravitacional ($\Delta f$) induzida por múltiplos planetas. Para a maioria dos sistemas simulados, este deslocamento está abaixo da resolução de frequência do LISA ($\sim 7.92 \times 10^{-9}$ Hz para uma missão de 4 anos). Apenas alguns sistemas ressonantes (População C) mostraram deslocamentos potencialmente detectáveis dentro de uma missão de 8 anos.
  • Aceleração do CoM: O estudo observou que ejeções planetárias podem induzir acelerações mensuráveis no centro de massa da binária (até $\sim 0.9$ au/ano$^2$ em casos específicos), o que poderia deixar uma marca no sinal de GW, embora um limiar de detecção específico para o LISA ainda não tenha sido definido.

Significância e Alegações
O artigo alega que sistemas planetários de segunda geração em torno de DWDs são candidatos dinamicamente viáveis para observação, mas suas arquiteturas finais são fortemente esculpidas por instabilidades dinâmicas precoces. Os autores enfatizam que:

1. A Sobrevivência é Comum, mas Modificada: Embora a disrupção completa seja rara (7,7% do total de sistemas), a perda de planetas é comum, levando a um aumento significativo (122%) na população de sistemas de dois planetas derivados de progenitores de alta multiplicidade.
2. NAMD como Diagnóstico: A métrica NAMD serve como uma ferramenta crucial para distinguir entre sistemas que se formaram com baixa multiplicidade e aqueles que foram dinamicamente excitados e perderam planetas. Valores altos de NAMD em sistemas de dois planetas observados podem indicar um histórico de instabilidade.
3. Potencial do LISA: O estudo apoia a viabilidade de detectar planetas "Magrathea" (planetas de segunda geração em torno de DWDs) com o LISA, particularmente sobreviventes de planeta único. No entanto, a detecção de sistemas de múltiplos planetas via modulação de frequência é improvável para a maioria das configurações devido aos pequenos deslocamentos de frequência, necessitando de estudos de forma de onda específicos e, potencialmente, dependendo de assinaturas de aceleração do CoM.

Os autores concluem que, embora suas simulações cubram a evolução dinâmica inicial (até $\sim 10$ Myr), os resultados fornecem uma estrutura fundamental para interpretar futuros dados do LISA e compreender o destino dos planetas em sistemas binários evoluídos. Eles declaram explicitamente que suas alegações de detectabilidade são ilustrativas e que análises de forma de onda dedicadas são necessárias para confirmar as perspectivas observacionais.