Kick Velocities and Mass Function of Free-Floating Planets from Dynamical Ejection in Hierarchical Three-Body Systems

Resumo Técnico: Velocidades de Chute e Função de Massa de Planetas Livres de Flutuação a partir de Expulsão Dinâmica em Sistemas de Três Corpos Hierárquicos

Enunciado do Problema
Planetas livres de flutuação (FFPs), ou planetas errantes, são objetos subestelares não ligados a nenhuma estrela hospedeira. Sua distribuição de velocidade e massa servem como fósseis cinemáticos de sua formação e história de expulsão. Embora o microlente gravitacional tenha detectado uma população de FFPs variando de massa de Terra a massa de Júpiter, os mecanismos dinâmicos que governam suas velocidades de expulsão e a função de massa resultante permanecem incompletamente caracterizados. Especificamente, há uma necessidade de mapear os resultados de expulsão de sistemas de três corpos hierárquicos (uma estrela hospedeira, um perturbador gigante massivo e um planeta interno mais leve) para as propriedades observáveis da população galáctica de FFPs, incluindo suas velocidades assintóticas ($v_\infty$) e as escalas de tempo de expulsão ($t_{ejec}$).

Metodologia
Os autores apresentam um conjunto de simulações diretas de $N$-corpos usando o pacote de integração rebound com o integrador IAS15 (passo de tempo adaptativo de 15ª ordem) para garantir alta precisão durante encontros próximos. O estudo foca no problema de três corpos restrito, tratando o planeta interno como uma partícula de teste em relação a uma estrela hospedeira de massa solar ($M_\star = 1 M_\odot$) e um perturbador gigante massivo.

A grade de simulação é controlada e sistemática:

• Parâmetros Fixos: Massa da hospedeira ($1 M_\odot$), semieixo maior do planeta gigante ($a_J = 5,2$ UA, um análogo de Júpiter) e inclinações iniciais ($0^\circ$–$3^\circ$).
• Parâmetros Variados:
  • Massa do planeta expelido ($m_p$): $10^{-5}$ a $10^{-1} M_{Júpiter}$.
  • Massa do planeta gigante ($M_J$): $1$a$75 M_{Júpiter}$.
  • Semieixo maior do planeta expelido ($a_p$): $0,7$a$1,5 a_J$.
  • Excentricidades de ambos os planetas ($e_p$ e $e_J$): $0,0$a$0,9$.
• Realizações: Cada configuração é simulada com $N=10.000$ realizações independentes usando sementes aleatórias compartilhadas para as anomalias médias, a fim de isolar os efeitos dos parâmetros.
• Métricas: O estudo mede a velocidade de expulsão assintótica ($v_\infty$), a velocidade de chute de raio finito ($v_{kick}$) no momento da confirmação da expulsão e a escala de tempo de expulsão ($t_{ejec}$). A expulsão é definida por energia específica positiva, distância da hospedeira excedendo $10 \times$ o semieixo maior máximo do sistema e velocidade radial externa persistente por dez períodos orbitais.

Estrutura Analítica
Os resultados numéricos são interpretados através de uma estrutura semianalítica combinando:

1. Orçamento de Energia: O incremento de energia necessário para a expulsão é derivado da diferença entre a energia orbital ligada e a energia cinética assintótica.
2. Estilingue Gravitacional: A velocidade de chute máxima é limitada pelo mecanismo de estilingue no referencial do gigante, dependente da velocidade de encontro $u$ e da velocidade orbital do gigante.
3. Espalhamento na Escala de Hill: A velocidade mediana de expulsão é governada por encontros típicos dentro da esfera de Hill do gigante, escalando como $v_H \propto M_J^{1/3}$.
4. Parâmetro de Tisserand: Usado como um diagnóstico conservado para relacionar elementos orbitais a velocidades de encontro.
5. Evolução Secular: As escalas de tempo para instabilidade estão ligadas ao bombeamento secular e à sobreposição de ressonância de movimento médio (MMR).

Resultados Principais

1. Limite de Partícula de Teste: A massa do planeta expelido ($m_p$) tem impacto negligenciável nas estatísticas de expulsão ao longo de quatro ordens de magnitude ($10^{-5}$ a $10^{-1} M_{Júpiter}$). A fração de expulsão varia em menos de 2,3%, e as velocidades médias permanecem estáveis. Isso confirma que, para $m_p/M_J \lesssim 10^{-2}$, a dinâmica é dominada pelo perturbador gigante, e o corpo expelido se comporta como uma partícula de teste.
2. Escalonamento da Massa do Planeta Gigante: A massa do planeta gigante ($M_J$) define a escala e o tempo de expulsão.
   • A velocidade média de expulsão escala como $\langle v_\infty \rangle \propto M_J^{0.42}$.
   • A escala de tempo de expulsão escala como $\langle t_{ejec} \rangle \propto M_J^{-1,09}$, consistente com a previsão secular $t_{sec} \propto M_J^{-1}$.
   • A fração de expulsão é não monotônica, atingindo o pico próximo a $10 M_{Júpiter}$ e declinando em massas muito baixas e muito altas devido às limitações de escala de tempo secular e captura na esfera de Hill, respectivamente.
3. Dependência do Semieixo Maior: Para configurações dentro da zona de sobreposição de ressonância ($|a_p - a_J|/a_J \lesssim 0,40$ para um perturbador de $10 M_{Júpiter}$), os tempos de expulsão são constantes ($\sim 1,9 \times 10^4$ anos), governados pela escala de tempo de sobreposição caótica rápida em vez da evolução secular. As velocidades médias atingem o pico perto de $a_p = 0,9 a_J$.
4. Efeitos de Excentricidade (Resposta de Duas Velocidades):
   • Excentricidade do Planeta Interno ($e_p$): Afeta primariamente a cauda de alta velocidade. Aumentar $e_p$ de 0 para 0,9 eleva a velocidade de chute máxima por um fator de $\sim 3$ (de $\sim 7,8$ para $\sim 23$ km/s), mas aumenta a velocidade média por apenas um fator de $\sim 1,7$.
   • Excentricidade do Planeta Gigante ($e_J$): Tem um efeito substancialmente mais forte na cauda. Aumentar $e_J$ de 0 para 0,9 eleva a velocidade de chute máxima por um fator de $\sim 10$ (de $\sim 7,8$ para $\sim 80$ km/s). Isso é impulsionado pelo crescimento não linear da velocidade do periastro do gigante ($v_{peri} \propto \sqrt{(1+e_J)/(1-e_J)}$), o que permite encontros de estilingue de alta energia raros.
5. Envelope Superior Analítico: Os autores derivam uma fórmula analítica para a velocidade máxima de expulsão baseada no mecanismo de estilingue e no parâmetro de Tisserand. Este envelope prevê com precisão o teto das distribuições de velocidade simuladas, particularmente para $e_p \gtrsim 0,5$.

Significância e Implicações
O artigo afirma que estes resultados fornecem uma ligação direta entre a arquitetura de sistemas planetários e as propriedades observáveis da população galáctica de FFPs:

• Função de Massa: Como a fração de expulsão é independente da massa do planeta, a função de massa atual de FFPs rastreia diretamente a taxa de ocorrência de planetas subjacente ($dN_{FFP}/dm_p \approx \bar{f}_{eject} \cdot dN_{planet}/dm_p$). Isso implica que, se a taxa de ocorrência aumenta em direção a massas menores (como sugerido pelos dados do Kepler), a população de FFPs sub-Terras será numericamente dominante.
• Assinaturas de Microlente: A distribuição transformada do tempo de escala de microlente ($dN/d\log t_E$) é prevista como plana em $\log t_E$ se a taxa de ocorrência de planetas for plana, ou aumentada em escalas de tempo curtas se a taxa de ocorrência aumentar em direção a massas menores. Isso oferece uma previsão testável para o Roman Wide Field Survey e o Euclid.
• Cinemática: A maioria dos planetas expulsos possui velocidades assintóticas de apenas alguns km/s, tornando-os cinematicamente indistinguíveis de sua população estelar progenitora (alterando a dispersão de velocidade 3D em $\sim 1\%$). No entanto, sistemas com perturbadores gigantes altamente excêntricos ($e_J \gtrsim 0,5$) podem produzir uma cauda de alta velocidade ($v_\infty \sim 80$ km/s), potencialmente criando discrepantes de movimento próprio detectáveis pelo Gaia ou Roman.
• Escalas de Tempo de Expulsão: A maioria das expulsões de sistemas massivos ou excêntricos ocorre rapidamente ($10^3$–$10^4$ anos), frequentemente durante a fase pré-sequência principal. Consequentemente, a distribuição espacial da população geral de FFPs é efetivamente indistinguível do campo, com apenas objetos jovens, recentemente expulsos, retendo a memória de seu ambiente de nascimento.

Limitações
Os autores observam que a restrição a exatamente três corpos é a principal limitação, visto que sistemas reais podem envolver múltiplos planetas levando a instabilidades em cascata. Além disso, as simulações não incluem discos protoplanetários, perda de massa estelar ou uma gama completa de inclinações (especificamente o mecanismo de Kozai-Lidov), o que poderia alterar as estatísticas de expulsão. O trabalho futuro é identificado como necessário para incorporar esses fatores em um modelo completo de síntese de população.