Nii-body: Bayesian Inference of Multiplanet Dynamics via N-body Simulations

1. O Problema

A maioria dos sistemas exoplanetários descobertos são configurações de múltiplos planetas cujas dinâmicas de longo prazo são governadas por interações gravitacionais de corpo N (N-body). Métodos tradicionais de ajuste orbital, que utilizam a superposição Kepleriana (somando órbitas de dois corpos independentes), são computacionalmente eficientes, mas falham catastróficamente em sistemas onde as interações planeta-planeta são fortes, como em ressonâncias de movimento médio ou arquiteturas compactas. Nessas situações, as perturbações mútuas acumulam-se ao longo do tempo, causando deriva nos observáveis (velocidade radial, astrometria, variações no tempo de trânsito) que o modelo Kepleriano não consegue reproduzir. Para caracterizar com precisão esses sistemas, é necessário um modelo de forward (avançado) de corpo N autoconsistente, integrado a métodos estatísticos Bayesianos. No entanto, a integração de N-body dentro de cadeias de Markov (MCMC) é computacionalmente custosa, tornando o ajuste de parâmetros um desafio significativo.

2. Metodologia

Os autores desenvolveram o Nii-body, um código de código aberto que implementa um framework de inferência Bayesiana acoplado a simulações de N-body. A metodologia baseia-se nos seguintes pilares:

• Integrador Numérico (RKF78): O código utiliza um solver Runge–Kutta–Fehlberg 7(8) adaptativo com controle de passo e tolerância de erro local. As equações de movimento newtonianas são resolvidas em coordenadas bariocêntricas. O código é escrito em C para garantir alta eficiência computacional e inclui bibliotecas de conversão entre coordenadas cartesianas e elementos orbitais Keplerianos.
• Algoritmo de Otimização (MCMC): O framework emprega um algoritmo MCMC de Temperamento Paralelo Automático (baseado na biblioteca Nii-C). Esta técnica utiliza múltiplas cadeias em diferentes "temperaturas" para escapar de mínimos locais e explorar o espaço de parâmetros globalmente. Um sistema de controle automático ajusta as distribuições de proposta durante uma fase de sintonia inicial, garantindo altas taxas de aceitação nas fases posteriores.
• Geração de Observáveis Sintéticos: O código gera sinais sintéticos de astrometria e velocidade radial (RV) diretamente a partir das integrações de N-body. O modelo considera o movimento de reflexo da estrela hospedeira em relação ao baricentro do sistema, projetando-o no plano do céu (para astrometria) e na linha de visada (para RV).
• Inferência Bayesiana: O código calcula a verossimilhança comparando os dados observados com os sinais sintéticos, utilizando priores uniformes para todos os parâmetros do modelo (massas, elementos orbitais, incerteza observacional).

3. Contribuições Principais

• Desenvolvimento do Nii-body: Apresentação de uma nova ferramenta de código aberto que combina um integrador RKF78 de alta precisão com um motor MCMC de temperamento paralelo automatizado, especificamente projetado para o ajuste de sistemas multiplanetários.
• Validação de Desempenho: Benchmarking do integrador RKF78 contra o pacote REBOUND (integrador IAS15) em um sistema de teste (Sol-Terra-Júpiter), demonstrando concordância numérica quase perfeita com alta eficiência computacional.
• Prova de Conceito em Astrometria: Demonstração da capacidade do código de recuperar parâmetros orbitais de um sistema de dois planetas (Kepler-9) em ressonância 2:1, utilizando apenas dados astrométricos sintéticos, sem necessidade de "chutes" iniciais (busca cega).
• Comparação Crítica: Evidência empírica de que modelos de superposição Kepleriana falham em recuperar parâmetros verdadeiros em sistemas ressonantes, enquanto o modelo de N-body succeeds.

4. Resultados

• Precisão do Integrador: Os testes de benchmark mostraram que o Nii-body/RKF78 produz trajetórias idênticas às do REBOUND/IAS15 (com erros quadráticos médios da ordem de $10^{-14}$ a $10^{-8}$ AU), com tempos de execução comparáveis ou ligeiramente inferiores.
• Recuperação de Parâmetros (Kepler-9): Em simulações com dados astrométricos sintéticos do sistema Kepler-9 (5 anos de dados, ruído de 1 µas):
  • O Nii-body conseguiu recuperar com sucesso os parâmetros planetários (massas, períodos, excentricidades, inclinações) em 6 de 8 execuções independentes (75% de taxa de sucesso) a partir de buscas cegas no espaço de 15 parâmetros.
  • O tempo de convergência foi de aproximadamente 24 horas em um laptop com processador Intel i7-12650H.
  • As execuções que não convergiram totalmente (Run 3 e Run 4) revelaram degenerescências entre parâmetros (ex: excentricidade vs. inclinação) ou a presença de ilhas de alta densidade no espaço de parâmetros, destacando a necessidade de múltiplas cadeias para confirmação.
• Falha do Modelo Kepleriano: Quando o mesmo conjunto de dados foi ajustado usando um modelo de superposição Kepleriana, nenhuma das 8 execuções conseguiu recuperar os parâmetros verdadeiros, resultando em intervalos de credibilidade extremamente amplos e falha de convergência.

5. Significado e Impacto

O artigo estabelece o Nii-body como uma ferramenta viável e necessária para a próxima geração de estudos de exoplanetas. À medida que missões futuras (como THEIA e CHES) e observatórios atuais (como o Gaia) alcançam precisão astrométrica na escala de microarcsegundos, a capacidade de modelar interações gravitacionais complexas em sistemas compactos torna-se crítica.

O trabalho demonstra que:

1. O ajuste de N-body é computacionalmente factível para sistemas de poucos planetas, mesmo com o custo inerente das integrações.
2. Para sistemas com interações dinâmicas fortes, o modelo Kepleriano é insuficiente e pode levar a conclusões errôneas sobre a arquitetura do sistema.
3. O framework é extensível para incluir dados combinados (RV, TTVs e Astrometria) e, futuramente, restrições de estabilidade dinâmica, tornando-se uma ferramenta versátil para a caracterização precisa de sistemas exoplanetários reais.