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

Resumen Técnico: Reevaluación del Control de Calidad de Datos y el Modelado de Sistemas Múltiples en Pruebas de Gravedad de Binarias Anchas

Planteamiento del Problema
Las estrellas binarias anchas sirven como laboratorios naturales para probar la gravedad en el régimen de baja aceleración ($g \lesssim 10^{-9} \text{ m s}^{-2}$), donde los efectos de la materia oscura son despreciables. Análisis estadísticos recientes de los datos de Gaia DR3 han reportado una anomalía gravitatoria consistente con la Dinámica Newtoniana Modificada (MOND), caracterizada por un factor de aumento $\gamma \approx 1.3-1.6$. Sin embargo, estudios contradictorios (p. ej., Banik et al. 2024; Pittordis et al. 2025; Cookson et al. 2026) afirman que se prefiere la gravedad Newtoniana estándar o que no existe desviación alguna. Estos resultados conflictivos dependen de dos problemas metodológicos primarios: (1) la implementación del control de calidad de los datos, específicamente el "corte de Banik" sobre el error de velocidad normalizado ($\sigma_{\tilde{v}}$), y (2) el modelado realista de sistemas de estrellas múltiples jerárquicos (tríos y múltiples de orden superior) que pueden imitar aumentos gravitatorios. Este artículo reexamina estos problemas para determinar si la anomalía de tipo MOND reportada es robusta frente a estas críticas específicas.

Metodología
Los autores llevan a cabo extensas pruebas estadísticas de gravedad utilizando velocidades proyectadas en el cielo en 2D ($v_p$) de tres muestras distintas: K.-H. Chae (2023), I. Banik et al. (2024) y C. Pittordis et al. (2025) (PSS). El estudio emplea tres marcos de prueba primarios:

1. Prueba del Plano de Aceleración: Un método de deproyección de Monte Carlo que analiza la desviación ortogonal ($\Delta_\perp$) de la aceleración observada respecto a las predicciones Newtonianas en el plano de aceleración logarítmico.
2. Prueba de la Distribución de $\tilde{v}$: Comparación del histograma de la velocidad 2D normalizada ($\tilde{v} = v_p / v_c(r_p)$) contra las predicciones teóricas para binarias puras y aquellas con compañeros ocultos.
3. Prueba del Perfil de la Mediana de $\tilde{v}$: Análisis de la escala de la mediana de $\tilde{v}$ como función de la separación normalizada ($r_p/r_M$).

Pasos metodológicos clave incluyen:

• Análisis de Sesgo de los Cortes de Calidad: Utilizando binarias Newtonianas ficticias con errores controlados, los autores simulan los efectos del corte de Banik ($\sigma_{\tilde{v}} < 0.1 \max(1, \tilde{v}/2)$) y otros cortes (p. ej., $\sigma_{v_p} < 50 \text{ m s}^{-1}$).
• Calibración del Modelado de Estrellas Múltiples: Los autores implementan el algoritmo de modelado de sistemas triples propuesto por Pittordis et al. (2025) para muestras que satisfacen $ruwe < 1.2$. Crucialmente, calibran la fracción de tríos ($f_{trip}$) utilizando datos del régimen de alta aceleración (Newtoniano), donde las alineaciones fortuitas son insignificantes, en lugar de ajustar $f_{trip}$ simultáneamente con los modelos de gravedad en el régimen de transición/MOND.
• Comparación con Soluciones Numéricas: Las predicciones teóricas se comparan tanto con modelos analíticos aproximados como con soluciones numéricas realistas recientes de QUMOND de J. Pflamm-Altenburg (2025) que consideran la dinámica de dos cuerpos y el efecto de campo externo (EFE) galáctico.
• Análisis de Submuestra de Alta Calidad: Se analiza una submuestra de binarias cercanas ($d < 150 \text{ pc}$) que cumple con el estricto marco de calidad de Cookson et al. (2026) para abordar las críticas de insuficiencia estadística y contaminación de datos.

Contribuciones Clave y Resultados

1. El Corte de Banik Introduce Sesgo:
   El estudio demuestra que el corte de Banik no es una medida de control de calidad neutral. Debido a que el error relativo $\sigma_{\tilde{v}}$ aumenta con la separación ($r_p$), el corte elimina selectivamente puntos de datos de alta precisión con valores de $\tilde{v}$ más altos en el régimen de baja aceleración. Esto suprime artificialmente la tendencia ascendente de $\tilde{v}$, sesgando los datos hacia las predicciones Newtonianas. En contraste, un corte basado en el error de velocidad absoluta ($\sigma_{v_p} < 50 \text{ m s}^{-1}$) elimina los datos de manera uniforme sin introducir tal sesgo.

2. Calibración de las Fracciones de Estrellas Múltiples:
   Cuando la fracción de tríos ($f_{trip}$) se calibra utilizando datos del régimen Newtoniano (donde las binarias deben obedecer la gravedad Newtoniana), los valores resultantes de $f_{trip}$ son significativamente menores que los derivados por Pittordis et al. (2025), quienes ajustaron $f_{trip}$ en los regímenes de transición/MOND.

   • Para la muestra de Chae (2023) con $ruwe < 1.2$, $f_{trip} \approx 0.10 \pm 0.03$.
   • Para la muestra PSS, $f_{trip} \approx 0.12 \pm 0.03$ (o $0.07$ dependiendo de las suposiciones de sobrevuelo/flyby).
     Utilizando estas fracciones calibradas, la prueba del plano de aceleración revela una anomalía gravitatoria significativa ($\delta_{obs-newt} > 0$) con una significancia estadística $> 5\sigma$ en el régimen de baja aceleración. La preferencia previamente reclamada por la gravedad Newtoniana por parte de Pittordis et al. (2025) se atribuye al uso de un $f_{trip}$ no calibrado y sobreestimado, así como a un rango dinámico limitado que excluyó el régimen Newtoniano puro.

3. Robustez de la Anomalía con Datos de Alta Calidad:
   Al analizar una gran muestra ($N \approx 480$ en el régimen MOND) de binarias cercanas ($d < 150 \text{ pc}$) que satisface los estrictos cortes de calidad de Cookson et al. (2026), los autores encuentran:

   • El perfil de la mediana de $\tilde{v}$ aumenta con la separación, desviándose de la predicción Newtoniana plana.
   • La afirmación de "no hay evidencia de MOND" de Cookson et al. (2026) resulta ser una consecuencia de la baja estadística numérica ($N=61$ en el régimen MOND), lo que permite que las fluctuaciones aleatorias imiten un perfil plano, combinado con el uso de masas estelares sesgadas y un referente Newtoniano incorrecto ($\tilde{v} \approx 0.49$).
   • Al corregir los sesgos de masa y utilizar una estadística suficiente, los datos muestran una desviación $> 3\sigma$ de la gravedad Newtoniana, consistente con un factor de aumento $\gamma \approx 1.4-1.5$.

4. Acuerdo con Modelos Numéricos Realistas:
   El comportamiento observado de la anomalía gravitatoria en el régimen de transición ($10^{-10} \lesssim g_N \lesssim 10^{-9} \text{ m s}^{-2}$) se alinea más estrechamente con las soluciones numéricas realistas de QUMOND (J. Pflamm-Altenburg 2025) que con los modelos de partículas de prueba aproximados. La distinción entre la gravedad Newtoniana y MOND se confirma como robusta cuando se aplican la calibración adecuada y el modelado realista.

Significancia y Reivindicaciones
El artículo concluye que la anomalía gravitatoria de baja aceleración reportada en los estudios de binarias anchas es una característica robusta de los datos, no un artefacto de la mala calidad de los datos o de estrellas múltiples no modeladas. Los autores afirman que:

• El "corte de Banik" es un criterio de selección sesgado que debilita artificialmente la señal de MOND.
• Las afirmaciones de preferencia Newtoniana en la literatura reciente derivan de fallos metodológicos, específicamente la falta de calibración de las fracciones de sistemas jerárquicos en el régimen Newtoniano y el uso de tamaños de muestra insuficientes en el régimen de baja aceleración.
• Cuando se corrigen estos problemas, las pruebas estadísticas sobre velocidades 2D confirman consistentemente una ruptura de la gravedad estándar a bajas aceleraciones, favoreciendo teorías de tipo MOND con un factor de aumento de $\gamma \approx 1.3-1.6$.
• Los resultados son consistentes con las soluciones numéricas realistas recientes para órbitas de binarias anchas en QUMOND, lo que sugiere que los análisis estadísticos refinados están convergiendo con los resultados de modelado orbital en 3D.

Los autores enfatizan que, si bien el modelado orbital en 3D jugará un papel cada vez más importante con las futuras liberaciones de datos de Gaia, las pruebas estadísticas en 2D actuales, cuando se controlan adecuadamente para el sesgo y la contaminación jerárquica, proporcionan una confirmación fuerte e independiente de la anomalía gravitatoria.