Precision Solar System Dynamics for Ultralight Dark Matter Search

Resumen Técnico: Dinámica de Precisión del Sistema Solar para la Búsqueda de Materia Oscura Ultraligera

Planteamiento del Problema
La materia oscura ultraligera (ULDM, por sus siglas en inglés), un candidato bosónico con masa $m \lesssim 10^{-20}$ eV, exhibe fluctuaciones de densidad de orden uno en la escala de su longitud de de Broglie. Estas fluctuaciones estocásticas pueden perturbar el movimiento de cuerpos gravitantes, incluyendo estrellas y planetas. Si bien búsquedas previas han utilizado matrices de temporización de púlsares (PTA) y astrometría para restringir la ULDM, estos métodos dependen de efemérides planetarias precisas para modelar el baricentro del sistema solar. Esto crea un bucle lógico de retroalimentación: la precisión requerida para modelar el sistema solar para las PTA sugiere que las mediciones de rango subyacentes utilizadas para construir estas efemérides podrían ser ellas mismas sensibles a las perturbaciones inducidas por la ULDM. Los autores investigan si las mediciones de rango de radio interplanetario, actualmente utilizadas para la navegación y la construcción de efemérides, pueden servir como una sonda directa para las fluctuaciones estocásticas de densidad de la ULDM.

Metodología
Los autores modelan el campo de la ULDM como un campo aleatorio gaussiano clásico con un componente de densidad estocástica. Tratan estas fluctuaciones como "cuasipartículas" transitorias de masa efectiva $M_{\text{eff}}$ y tamaño $\lambda_{\text{DM}} \sim 1/(m\sigma)$, donde $\sigma$ es la dispersión de velocidad.

1. Modelado de la Señal: Los autores analizan el observable de rango interplanetario de doble vía, $\Delta \tau$, descomponiéndolo en perturbaciones orbitales, retraso de Shapiro y retraso temporal de Einstein. Mediante estimaciones de orden de magnitud y cálculos detallados del espectro de potencia (Apéndice A), demuestran que la señal de perturbación orbital domina sobre los efectos de retraso temporal relativista para los parámetros de interés. La señal se modela como una fuerza estocástica actuando sobre masas de prueba (planetas) en relación con el baricentro del sistema solar.
2. Marco Analítico: Asumiendo que las masas de prueba se comportan como partículas libres en escalas de tiempo cortas, los autores derivan el espectro de potencia de la señal y la densidad espectral cruzada entre dos mediciones de rango (por ejemplo, Tierra-Marte y Tierra-Venus). Definen un coeficiente de correlación, $\Gamma_{ab}$, que codifica la correlación espacial de la señal de la ULDM basada en la separación de los planetas y la longitud de onda de la ULDM.
3. Validación Numérica: Para abordar la no linealidad introducida por la aceleración gravitacional del Sol en escalas de tiempo largas, los autores realizan simulaciones numéricas del sistema solar (incluyendo los seis planetas interiores con correcciones post-newtonianas). Segmentan la serie temporal en intervalos más cortos que el periodo orbital ($T_{\text{seg}} < 1$ año) para aproximar el sistema como lineal. Los resultados numéricos para las autocorrelaciones y correlaciones cruzadas muestran concordancia con las predicciones analíticas derivadas bajo la aproximación de partícula libre.
4. Proyección de Sensibilidad: Utilizando un estadístico de correlación cruzada óptimo (adaptado de las búsquedas de ondas gravitacionales estocásticas), los autores proyectan la sensibilidad de una red de mediciones de rango. Asumen un conjunto de datos de referencia de $N=5$ mediciones de rango independientes (de la Tierra a Mercurio, Venus, Marte, Saturno y Júpiter) sobre una línea de base de 20 años, con características de ruido blanco que reflejan las incertidumbres de medición actuales ($\sigma_L \sim 1$ m).

Resultados Clave

• Pico de Sensibilidad: La sensibilidad proyectada alcanza su máximo en una masa de ULDM de $m \simeq 10^{-15}$ eV. En esta masa, la longitud de onda de la materia oscura ($\lambda_{\text{DM}} \approx 2.5$ UA) es comparable a las distancias interplanetarias, maximizando la respuesta.
• Restricciones de Densidad: Para los niveles de ruido de referencia, los autores muestran que las mediciones de rango actuales podrían sondear una densidad de materia oscura local de $\bar{\rho} \sim 10^5 \rho_0$ (donde $\rho_0 \approx 0.4$ GeV/cm$^3$) en la masa pico.
• Verificación de Correlación: Las simulaciones numéricas confirman que los coeficientes de correlación analíticos describen con precisión la señal para segmentos de tiempo más cortos de un año. Se encontró que la correlación entre los rangos Tierra-Marte y Tierra-Venus es consistente con las predicciones teóricas.
• Análisis de Ruido: Los autores identifican fuentes potenciales de ruido correlacionado, incluyendo la estabilidad de los relojes atómicos, las incertidumbres en el etiquetado de tiempo y los retrasos del plasma solar. Concluyen que los errores de reloj y de etiquetado de tiempo son insignificantes (escala sub-milimétrica) en comparación con el ruido de referencia. Aunque el plasma solar y las perturbaciones del cinturón de asteroides son relevantes, los autores argumentan que el piso de ruido actual está dominado por otros factores, y que la morfología espacial específica de la ULDM difiere del ruido inducido por asteroides, lo que sugiere que los límites siguen siendo robustos.

Significancia y Reivindicaciones
El artículo afirma establecer un nuevo canal para la búsqueda de materia oscura ultraligera utilizando la dinámica de precisión del sistema solar actual y futura.

• Complementariedad: Los resultados complementan las restricciones existentes de la temporización de púlsares (que son sensibles a masas más bajas, $m \sim 10^{-17}$ eV) y de los futuros detectores de ondas gravitacionales basados en el espacio como LISA. El enfoque de medición de rango sonda un rango de masa superior ($m \sim 10^{-15}$ eV) donde la longitud de onda de la materia oscura es comparable a la escala del sistema solar.
• Innovación Metodológica: El trabajo demuestra que los "residuales" en los datos de rango planetario, típicamente tratados como ruido o errores sistemáticos a minimizar para la navegación, contienen una estructura de señal estocástica que puede ser explotada mediante estadísticas de correlación cruzada.
• Limitaciones: Los autores señalan modestamente que su sensibilidad está limitada por la necesidad de segmentar los datos en intervalos cortos ($<1$ año) para evitar la dinámica orbital no lineal, lo que establece un límite inferior para la masa detectable ($m \gtrsim \sigma^2/T_{\text{seg}}$). Reconocen que los límites actuales ($\rho \sim 10^5 \rho_0$) están muy por encima de la densidad local esperada, pero las mejoras futuras en la precisión de la medición de rango (por ejemplo, BepiColombo apuntando a una precisión de nivel centimétrico) podrían mejorar significativamente estos límites.

El artículo concluye que, si bien los datos actuales aún no restringen la densidad local estándar de materia oscura, la metodología proporciona una vía viable para búsquedas futuras a medida que la precisión de la medición de rango mejore, potencialmente descubriendo nuevo espacio de parámetros para la materia oscura ultraligera.