White dwarf planets in star clusters: gravitational scattering versus mass-loss effects

Resumo Técnico: Planetas de anãs brancas em aglomerados estelares: espalhamento gravitacional versus efeitos de perda de massa

Definição do Problema
As anãs brancas (WDs) servem como laboratórios únicos para compreender a sobrevivência e a evolução a longo prazo de sistemas planetários. Embora a perda de massa pós-sequência principal seja conhecida por alterar as propriedades orbitais planetárias e agitar discos de detritos (levando à poluição atmosférica), a influência do ambiente de nascimento estelar sobre esses sistemas permanece mal definida. Especificamente, não está claro como as interações dinâmicas dentro de regiões densas de formação estelar (como o espalhamento gravitacional e passagens próximas de estrelas/flybys) competem com ou modificam a evolução orbital impulsionada pela perda de massa estelar. Estudos anteriores modelaram interações dinâmicas em aglomerados ou a evolução estelar isoladamente, mas poucos acoplaram de forma autossuficiente a evolução dinâmica N-body com a evolução estelar desde a formação da região de formação estelar até que as estrelas progenitoras se tornem anãs brancas.

Metodologia
Os autores realizaram simulações N-body de regiões de formação estelar usando o integrador kira dentro do ambiente Starlab. As simulações evoluíram por 1 Gyr, um intervalo de tempo suficiente para que as estrelas mais massivas ($>2 M_\odot$) deixem a sequência principal e se tornem anãs brancas.

• Condições Iniciais: As simulações utilizaram $N_\star = 1000$ estrelas extraídas de uma Função de Massa Inicial (IMF) de Maschberger (2013). Cada estrela com massa inicial $M_\star \leq 2.5 M_\odot$ recebeu um único planeta de massa de Júpiter. Os semieixos maiores iniciais ($a_p$) foram amostrados a partir de uma distribuição uniforme (0,1–50 au) ou funções delta em valores específicos (0,1, 1, 5 e 30 au). As excentricidades e inclinações iniciais foram definidas como zero.
• Dinâmica de Aglomerados: Dois tipos de configurações iniciais de aglomerados foram testados:
  1. Distribuições fractais: Geradas usando um método de fractal de caixa com uma dimensão fractal $D=1,6$ para representar regiões de formação estelar altamente subestruturadas e ligadas. Os raios iniciais foram definidos como 1 pc (alta densidade, $\sim 10^4 M_\odot \text{pc}^{-3}$) e 5 pc (baixa densidade, $\sim 10^2 M_\odot \text{pc}^{-3}$).
  2. Esferas de Plummer: Distribuições suaves e centralmente concentradas com um raio de semimassa de 0,1 pc.
• Evolução Estelar: A evolução estelar e binária foi implementada usando tabelas de consulta SeBa, sincronizadas com o passo de tempo (timestep) do N-body. O código rastreou a perda de massa durante as fases de ramo gigante. Os planetas foram considerados engolfados se o raio estelar excedesse o semieixo maior do planeta. As simulações não incluíram binárias estelares primordiais ou campos de maré galácticos externos.
• Conjuntos de Simulação: Dez realizações foram executadas para cada conjunto de condições iniciais para avaliar a estocasticidade. O estudo focou no destino dos planetas orbitando progenitores que evoluem para anãs brancas (massas iniciais de 1–2,5 $M_\odot$).

Resultos Principais

1. Dominância da Perda de Massa: O principal motor da evolução orbital para os planetas sobreviventes é a perda de massa estelar, não o espalhamento dinâmico. Embora ambientes densos destruam ou alterem significativamente as órbitas de até 50% dos sistemas planetários via espalhamento, o efeito da perda de massa domina a dinâmica para os sistemas restantes. Este comportamento é amplamente independente da densidade estelar inicial da região de nascimento.
2. Expansão do Semieixo Maior: Para os planetas que sobrevivem ao engolfamento e permanecem ligados, o semieixo maior aumenta significativamente devido à perda de massa. No limite adiabático, o semieixo maior final escala inversamente com a massa estelar final ($a_f \propto M_{inicial}/M_{final}$). Como os progenitores perdem $\sim$50–80% de sua massa, as órbitas planetárias expandem-se por fatores de 2–5. Esta expansão cria uma característica distinta na distribuição do semieixo maior final (com pico em torno de 40–200 au para distribuições iniciais de 0,1–50 au) que não está presente na distribuição de todos os planetas (incluindo aqueles ao redor de estrelas de sequência principal) aos 100 Myr.
3. Mudanças de Excentricidade: Aproximadamente 40% das anãs brancas experimentam um aumento na excentricidade, comparado a 30% de todos os planetas. Crucialmente, estas mudanças de excentricidade são quase inteiramente impulsionadas por encontros dinâmicos ocorrendo antes da fase de anã branca. Uma vez que a estrela se torna uma anã branca, a perda de massa sozinha induz mudanças de excentricidade negligenciáveis para planetas dentro de $\sim$100 au.
4. Captura e Sistemas Triplos:
   • Planetas Capturados: Aproximadamente 10% dos sistemas planetários de anãs brancas finais consistem de planetas capturados (planetas de massa livre capturados por uma estrela). Estes sistemas tipicamente exibem grandes semieixos maiores ($>100$ au) e altas excentricidades ($e > 0,5$), assemelhando-se a uma distribuição térmica.
   • Sistemas Triplos: As simulações produziram sistemas triplos de longa duração contendo uma anã branca, um planeta e um terceiro corpo (seja uma estrela ou um planeta capturado). Estas configurações são potenciais análogos a sistemas como PSR B1620-26 (AB) b.
5. Taxas de Destruição: Em regiões densas, a taxa de destruição de sistemas planetários é alta. Para planetas orbitando progenitores de anãs brancas, a fração de planetas ligados sobreviventes cai para $\sim$52% após 1 Gyr, comparado a $\sim$71% para todas as estrelas. Isto é atribuído às taxas de interação mais elevadas experimentadas por estrelas mais massivas.

Significância e Conclusões
O artigo conclui que o ambiente dinâmico do aglomerado de nascimento tem um impacto modesto nas propriedades orbitais da maioria das anãs brancas sobreviventes. Especificamente:

• Para $\sim$80% dos sistemas sobreviventes, a configuração orbital atual é determinada pela perda de massa pós-sequência principal, em vez do espalhamento do aglomerado de nascimento.
• Para os $\sim$20% restantes, encontros dinâmicos alteram significativamente o semieixo maior.
• A presença de planetas de separação larga e alta excentricidade (ex: WD 0806-661 b) pode ser indicativa de uma origem por captura dentro de um aglomerado, um processo mais eficiente em aglomerados do que no campo Galáctico.
• Por outro lado, planetas muito próximos (ex: WD 1856+534 b) provavelmente chegaram às suas órbitas atuais via mecanismos de espalhamento (potencialmente envolvendo outros planetas ou estrelas) que ocorreram cedo na vida do aglomerado, seguidos de truncamento de maré, em vez de apenas a expansão por perda de massa.

Os autores fornecem uma síntese de população de planetas gigantes (1–100 au) ao redor de anãs brancas emergentes de aglomerados de formação, oferecendo restrições para observações futuras por instalações como o Telescópio Espacial Nancy Grace Roman, Gaia e JWST. O estudo enfatiza que, embora os ambientes de nascimento moldem a população através de destruição e captura, a expansão orbital fundamental dos sistemas sobreviventes é uma consequência robusta da evolução estelar.