Revisiting Data Quality Control and Multiple-star Modeling in Wide Binary Gravity Tests: Confirmation of MOND-type Gravitational Anomaly at Low Acceleration

Resumo Técnico: Revisitando o Controle de Qualidade de Dados e a Modelagem de Múltiplas Estrelas em Testes de Gravidade de Binárias Longejas

Definição do Problema
Estrelas binárias longejas servem como laboratórios naturais para testar a gravidade no regime de baixa aceleração ($g \lesssim 10^{-9} \text{ m s}^{-2}$), onde os efeitos da matéria escura são negligenciáveis. Análises estatísticas recentes dos dados do Gaia DR3 relataram uma anomalia gravitacional consistente com a Dinâmica Newtoniana Modificada (MOND), caracterizada por um fator de reforço $\gamma \approx 1,3-1,6$. No entanto, estudos conflitantes (ex: Banik et al. 2024; Pittordis et al. 2025; Cookson et al. 2026) afirmam que a gravidade Newtoniana padrão é preferida ou que nenhuma desviação existe. Esses resultados conflitantes baseiam-se em dois problemas metodológicos primários: (1) a implementação do controle de qualidade de dados, especificamente o "corte de Banik" na velocidade normalizada ($\sigma_{\tilde{v}}$), e (2) a modelagem realista de sistemas de múltiplas estrelas hierárquicas (tríplas e múltiplas de ordem superior) que podem mimetizar reforços gravitacionais. Este artigo reexamina essas questões para determinar se a anomalia do tipo MOND relatada é robusta contra essas críticas específicas.

Metodologia
Os autores conduzem extensos testes estatísticos de gravidade usando velocidades 2D projetadas no céu ($v_p$) de três amostras distintas: K.-H. Chae (2023), I. Banik et al. (2024) e C. Pittordis et al. (2025) (PSS). O estudo emprega três estruturas de teste primárias:

1. Teste do Plano de Aceleração: Um método de deprojeção de Monte Carlo analisando o desvio ortogonal ($\Delta_\perp$) da aceleração observada das previsões Newtonianas no plano logarítmico de aceleração.
2. Teste da Distribuição de $\tilde{v}$: Comparação do histograma da velocidade 2D normalizada ($\tilde{v} = v_p / v_c(r_p)$) contra previsões teóricas para binárias puras e aquelas com companheiros ocultos.
3. Teste do Perfil da Mediana de $\tilde{v}$: Análise da escala da mediana de $\tilde{v}$ como uma função da separação normalizada ($r_p/r_M$).

Principais etapas metodológicas incluem:

• Análise de Viés de Cortes de Qualidade: Usando binárias Newtonianas fictícias com erros controlados, os autores simulam os efeitos do corte de Banik ($\sigma_{\tilde{v}} < 0,1 \max(1, \tilde{v}/2)$) e outros cortes (ex: $\sigma_{v_p} < 50 \text{ m s}^{-1}$).
• Calibração de Modelagem de Triplas: Os autores implementam o algoritmo de modelagem de estrelas triplas proposto por Pittordis et al. (2025) para amostras que satisfazem $ruwe < 1,2$. Crucialmente, eles calam a fração de triplas ($f_{trip}$) usando dados do regime de alta aceleração (Newtoniano), onde alinhamentos casuais são negligenciáveis, em vez de ajustar $f_{trip}$ simultaneamente com modelos de gravidade no regime de transição/MOND.
• Comparação com Soluções Numéricas: As previsões teóricas são comparadas tanto com modelos analíticos aproximados quanto com soluções numéricas realistas recentes de QUMOND por J. Pflamm-Altenburg (2025) que consideram a dinâmica de dois corpos e o efeito de campo externo galáctico (EFE).
• Análise de Subamostra de Alta Qualidade: Uma subamostra de binárias próximas ($d < 150 \text{ pc}$) que satisfaz o rigoroso framework de qualidade de Cookson et al. (2026) é analisada para abordar as alegações de estatística insuficiente ou contaminação de dados.

Principais Contribuições e Resultados

1. O Corte de Banik Introduz Viés:
   O estudo demonstra que o corte de Banik não é uma medida de controle de qualidade neutra. Como o erro relativo $\sigma_{\tilde{v}}$ aumenta com a separação ($r_p$), o corte remove seletivamente pontos de dados de alta precisão com valores de $\tilde{v}$ mais altos no regime de baixa aceleração. Isso suprime artificialmente a tendência crescente de $\tilde{v}$, enviesando os dados para as previsões Newtonianas. Em contraste, um corte baseado no erro de velocidade absoluta ($\sigma_{v_p} < 50 \text{ m s}^{-1}$) remove dados uniformemente sem introduzir tal viés.

2. Calibração de Frações de Múltiplas Estrelas:
   Quando a fração de triplas ($f_{trip}$) é calibrada usando dados do regime Newtoniano (onde as binárias devem obedecer à gravidade Newtoniana), os valores de $f_{trip}$ resultantes são significativamente menores do que os derivados por Pittordis et al. (2025), que ajustaram $f_{trip}$ nos regimes de transição/MOND.

• Para a amostra de Chae (2023) com $ruwe < 1,2$, $f_{trip} \approx 0,10 \pm 0,03$.
• Para a amostra PSS, $f_{trip} \approx 0,12 \pm 0,03$ (ou $0,07$ dependendo das suposições de flyby).
  Usando essas frações calibradas, o teste do plano de aceleração revela uma anomalia gravitacional significativa ($\delta_{obs-newt} > 0$) com significância estatística $> 5\sigma$ no regime de baixa aceleração. A preferência anteriormente reivindicada pela gravidade Newtoniana por Pittordis et al. (2025) é atribuída ao uso de um $f_{trip}$ não calibrado e superestimado, além de um alcance dinâmico limitado que excluiu o regime Newtoniano puro.

3. Robustez da Anomalia com Dados de Alta Qualidade:
   Ao analisar uma grande amostra ($N \approx 480$ no regime MOND) de binárias próximas ($d < 150 \text{ pc}$) que satisfaz os cortes de qualidade rigorosos de Cookson et al. (2026), os autores encontram:

• O perfil da mediana de $\tilde{v}$ sobe com a separação, desviando-se da previsão Newtoniana plana.
• A afirmação de "nenhuma evidência para MOND" por Cookson et al. (2026) é mostrada como uma consequência de baixa estatística numérica ($N=61$ no regime MOND), o que permite que flutuações aleatórias mimetizem um perfil plano, combinada com o uso de massas estelares enviesadas e um benchmark Newtoniano incorreto ($\tilde{v} \approx 0,49$).
• Quando corrigido para vieses de massa e utilizando estatísticas suficientes, os dados mostram um desvio $> 3\sigma$ da gravidade Newtoniana, consistente com um fator de reforço $\gamma \approx 1,4-1,5$.

4. Concordância com Modelos Numéricos Realistas:
   O comportamento observado da anomalia gravitacional no regime de transição ($10^{-10} \lesssim g_N \lesssim 10^{-9} \text{ m s}^{-2}$) alinha-se mais proximamente com as soluções numéricas realistas de QUMOND (J. Pflamm-Altenburg 2025) do que com modelos de partícula de teste aproximados. A distinção entre a gravidade Newtoniana e a MOND é confirmada como sendo robusta quando a calibração adequada e a modelagem realista são aplicadas.

Significância e Alegações
O artigo conclui que a anomalia gravitacional de baixa aceleração relatada em estudos de binárias longejas é uma característica robusta dos dados, não um artefato de má qualidade de dados ou de múltiplas estrelas não modeladas. Os autores afirmam que:

• O "corte de Banik" é um critério de seleção enviesado que enfraquece artificialmente o sinal MOND.
• As alegações de preferência Newtoniana na literatura recente derivam de falhas metodológicas, especificamente a falta de calibração das frações de sistemas hierárquicos no regime Newtoniano e o uso de tamanhos de amostra insuficientes no regime de baixa aceleração.
• Quando esses problemas são corrigidos, testes estatísticos em velocidades 2D confirmam consistentemente uma quebra da gravidade padrão em baixas acelerações, favorecendo teorias do tipo MOND com um fator de reforço de $\gamma \approx 1,3-1,6$.
• Os resultados são consistentes com soluções numéricas realistas recentes para órbitas de binárias longejas em QUMOND, sugerindo que análises estatísticas refinadas estão convergindo com os resultados de modelagem orbital 3D direta.

Os autores enfatizam que, embora a modelagem orbital 3D desempenhará um papel cada vez mais importante com os futuros lançamentos de dados do Gaia, os testes estatísticos 2D atuais, quando devidamente controlados para viés e contaminação hierárquica, fornecem uma confirmação forte e independente da anomalia gravitacional.