Bayesian analysis of Gaia epoch astrometry and radial velocities with kima

Resumo Técnico: Análise Bayesiana de Astrometria de Época e Velocidades Radiais da Gaia com o kima

Problema e Contexto
A próxima liberação de dados (DR4) do telescópio espacial Gaia, prevista para dezembro de 2026, fornecerá astrometria de época individual para aproximadamente dois bilhões de fontes. Este conjunto de dados promete a detecção de dezenas de milhares de novos exoplanetas, bem como restrições orbitais 3D para estrelas binárias e buracos negros. No entanto, analisar este volume de dados requer ferramentas dedicadas capazes de lidar com ajustes orbitais complexos, comparação de modelos e o processamento simultâneo de dados astrométricos e de velocidade radial (RV). Ferramentas existentes como kepmodel, BINARYS e Octofitter podem ajustar a astrometria de época da Gaia, mas carecem das capacidades específicas para amostragem trans-dimensional eficiente e comparação robusta de modelos Bayesianos inerentes ao código-fonte do kima.

Metodologia
Os autores estendem o código de ajuste de órbita de código aberto kima (originalmente desenvolvido para análise de RV) através da implementação de dois novos modelos:

1. GAIAmodel: Projetado para analisar apenas a astrometria de época da Gaia. Ele modela a posição ao longo da varredura ($w$) como uma soma de um modelo de estrela única de linha de base (5 parâmetros: deslocamentos de posição, movimento próprio, paralaxe), termos de aceleração opcionais (7 ou 9 parâmetros), sinais dependentes do ângulo de varredura (harmônicos do ângulo de varredura $\psi$) e componentes orbitais Keplerianos.
2. RVGAIAmodel: Combina o GAIAmodel com dados de velocidade radial. Ele vincula a semi-amplitude de RV ($K$) e o semieixo maior do fotocentro ($a$) através de uma relação física, compartilhando parâmetros orbitais (período, excentricidade, argumento do periastro, etc.) entre os dois conjuntos de dados, permitindo parâmetros de jitter e tendência independentes.

Ambos os modelos utilizam o algoritmo de amostragem aninhada difusiva (via o pacote DNEST4) para explorar o espaço de parâmetros. Uma característica fundamental é que o número de sinais Keplerianos ($N_p$) é tratado como um parâmetro livre, permitindo a comparação de modelos Bayesianos eficiente (cálculo de evidências) e estimativa de parâmetros em uma única execução. Os modelos suportam as parametrizações de Thiele-Innes e geométrica, sendo esta última preferida para o RVGAIAmodel para lidar com degenerações de inclinação.

Principais Contribuições

• Implementação de Software: A adição de modelos personalizados ao kima para analisar a astrometria de época da Gaia, seja de forma isolada ou conjunta com RVs.
• Discriminação de Sinal de Ângulo de Varredura: A inclusão de um modelo específico para contabilizar sinais dependentes do ângulo de varredura (usando 3º, 5º e 7º harmônicos), permitindo a distinção entre órbitas Keplerianas reais e artefatos causados por binárias marginalmente resolvidas.
• Limites de Compatibilidade: Uma demonstração de como o kima pode gerar "limites de compatibilidade" (um equivalente Bayesiano aos limites de detecção) diretamente das amostras posteriores de sinais não detectados, aplicável tanto a dados simulados quanto reais.

Resultados
Os autores testaram os modelos em dados simulados (derivados de resíduos reais da Gaia com sinais injetados) e dados reais do sistema Gaia BH3:

• Validação: Os modelos recuperaram com sucesso sinais planetários injetados (alta e baixa relação sinal-ruído, períodos longos e curtos) com distribuições posteriores consistentes com os valores reais. Os resultados foram consistentes com os obtidos usando o kepmodel.
• Comparação de Modelos: Em testes envolvendo um sinal conhecido de ângulo de varredura (Target 4), o modelo com harmônicos de ângulo de varredura foi decisivamente favorecido sobre o modelo Kepleriano (Fator de Bayes $\approx 640$). Por outro lado, para um Kepleriano de baixa S/N (Target 6), o modelo Kepleriano foi favorecido sobre os modelos de ângulo de varredura. No entanto, para um Kepleriano de alta S/N (Target 5), o modelo de ângulo de varredura foi ligeiramente favorecido, ilustrando que distinguir esses sinais nem sempre é definitivo sem acompanhamento.
• Parametrização: Comparações entre as parametrizações de Thiele-Innes e geométrica mostraram distribuições posteriores consistentes, mas notaram que a escolha da parametrização afeta o Fator de Bayes devido ao "fator de Occam" (razão volume prior-posterior).
• Aplicação em Dados Reais: Os modelos foram aplicados ao Gaia BH3. Os autores observaram que para razões sinal-ruído muito altas ($S/N \gtrsim 3000$), o modo "Known-Object" com priors mais estreitos é recomendado em vez de buscas cegas para garantir que o amostrador aninhado localize o pico agudo de verossimilhança.
• Limites de Compatibilidade: O estudo gerou limites de compatibilidade para período vs. massa. Os resultados mostraram formas esperadas, incluindo um mínimo em torno do intervalo de tempo dos dados, e destacaram potenciais degenerações onde órbitas parcialmente cobertas ou múltiplos sinais podem mimetizar soluções Keplerianas únicas.

Significância
O artigo afirma que o kima fornece uma ferramenta eficiente e de código aberto para a análise da astrometria de época da Gaia, sendo único em sua capacidade de ajustar simultaneamente qualquer número de Keplerianos e realizar comparação de modelos rigorosa em uma única execução. Os autores afirmam que a ferramenta é validada para analisar dados da Gaia, capaz de distinguir órbitas reais de artefatos de ângulo de varredura em casos específicos, e adequada para calcular taxas de ocorrência e limites de compatibilidade. O trabalho posiciona o kima como uma ferramenta crítica para a análise futura do DR4 da Gaia, particularmente para estudos que exigem a análise conjunta de astrometria e velocidades radiais.