Impact of Spacecraft Orbit Uncertainties and Velocity Mismodeling on the LISA Gravitational-Wave Response

Resumen Técnico: Impacto de las Incertidumbres en la Órbita de la Nave Espacial y el Error de Modelado de la Velocidad en la Respuesta de LISA a las Ondas Gravitacionales

Planteamiento del Problema
El Interferómetro de Ondas Gravitacionales por Láser y Espacio (LISA) está diseñado para detectar ondas gravitacionales (OG) en la banda de frecuencia de milihertz utilizando una constelación de tres naves espacales. El análisis preciso de los datos de OG depende críticamente de la función de respuesta del detector, que mapea la deformación de la OG entrante a la señal de fracción de frecuencia medida. Esta respuesta depende fundamentalmente de las posiciones, velocidades y longitudes de los brazos de las naves espaciales. Si bien estudios previos han abordado la calibración instrumental y la simulación de trayectorias, dos preguntas específicas sobre la fidelidad del modelo de respuesta permanecieron sin resolver:

1. ¿Cómo se propagan las incertidumbres en el conocimiento de las órbitas y las longitudes de los brazos de las naves espaciales en la respuesta de las OG?
2. ¿Cuál es el impacto de omitir los términos de velocidad de las naves espaciales (efecto Doppler/Doppler boosting) en el modelo de respuesta, particularmente para fuentes de binarias galácticas?

Metodología
Los autores cuantifican estos errores utilizando la métrica de "desajuste" (mismatch), la cual mide la pérdida fraccional en la relación señal-ruido (SNR) cuando se utiliza un modelo incorrecto. El estudio emplea dos escenarios orbitales complementarios:

1. Modelo de Juguete Estático: Se utiliza una configuración triangular simplificada e independiente del tiempo para aislar los efectos geométricos. Se aplican secuencialmente tres clases de perturbaciones:

   • Variación de la longitud del brazo: Perturbaciones gaussianas independientes ($\sigma_L = 1$ m) para simular un distanciamiento interferométrico a bordo de alta precisión.
   • Rotación de cuerpo rígido: Rotaciones aleatorias que inducen incertidumbres posicionales equivalentes a $\sigma_{rot} = 50$ km.
   • Traslación de cuerpo rígido: Traslaciones globales con $\sigma_{trans} = 50$ km para simular incertidumbres de posición absoluta provenientes del seguimiento radiométrico basado en tierra.

2. Escenario Evolutivo Realista: Una simulación de alta fidelidad utilizando los paquetes de software GODOT y MIDAS de la ESA. Esto implica:

   • Generar una trayectoria de "verdad" con perturbaciones realistas (autogravedad, fuerzas no gravitatorias, ruido de reloj).
   • Ejecutar un Filtro de Información de Raíz Cuadrada (SRIF) aumentado para estimar la trayectoria basada en observaciones simuladas (rango, Doppler, pseudorangos).
   • Muestrear a partir de la matriz de covarianza orbital resultante para generar un conjunto de estimaciones de trayectoria.
   • Comparar una respuesta de LISA calculada con estas trayectorias perturbadas frente a una respuesta nominal.

3. Error de Modelado de la Respuesta: El estudio compara la respuesta completa con efecto Doppler (incluyendo los términos de velocidad $\dot{r}/c$) frente a una aproximación de orden cero donde se omiten las velocidades de las naves espaciales. Esto se evalúa utilizando tanto una fuente monocromática de frecuencia única como un modelo de forma de onda de binaria galáctica rápida (jaxgb).

4. Estimación de Parámetros: Se realiza un análisis de Monte Carlo por Cadenas de Markov (MCMC) sobre una señal de binaria galáctica ($f = 10^{-4}$ Hz, SNR = 200) observada durante un año para evaluar los sesgos de parámetros resultantes del error de modelado de la velocidad.

Resultos Clave

• Incertidumbres de la Órbita:

  • En el escenario evolutivo realista, las incertidumbres en la determinación de la órbita producen desajustes insignificantes en toda la banda de LISA. El peor caso de desajuste es inferior a $10^{-7}$.
  • Los errores están dominados por efectos de alta frecuencia ($>10$ mHz) donde la respuesta es más sensible a las posiciones absolutas de las naves espaciales.
  • Comparado con el modelo de juguete estático, la simulación realista muestra que los errores de longitud de brazo son de 3 a 6 veces menores, los de posición de 8 a 20 veces menores, y los de ángulo aproximadamente 300 veces menores.

• Error de Modelado de la Velocidad:

  • Omitir las velocidades de las naves espaciales conduce a errores significativamente mayores, particularmente en bajas frecuencias.
  • A $\sim 10^{-4}$ Hz, los desajustes son del orden de $10^{-4}$, con errores de amplitud relativa alrededor de $10^{-2}$ y errores de fase que alcanzan un radián.
  • El error escala aproximadamente como $f^{-1}$ en bajas frecuencias, disminuyendo a $\sim 10^{-9}$ a 100 mHz.
  • El error de modelado de la velocidad domina el presupuesto de error por aproximadamente dos órdenes de magnitud a 10 mHz y nueve órdenes de magnitud a 0.1 mHz en comparación con las incertidumbres orbitales.

• Estimación de Parámetros:

  • Para una binaria galáctica con $f = 10^{-4}$ Hz y SNR = 200, omitir los términos de velocidad en la respuesta conduce a sesgos de parámetros menores a 1$\sigma$.
  • La diferencia en la media del log-verosimilitud entre la respuesta completa y el modelo que omite la velocidad es $\Delta \log \mathcal{L} \sim 10$.

Significancia y Reivindicaciones
Este trabajo proporciona la primera caracterización de cómo los errores en la respuesta de las ondas gravitacionales de LISA se propagan desde la deformación de la OG a través del detector hasta los parámetros estimados. Los autores concluyen que:

• Conocimiento de la Órbita: Las capacidades actuales de determinación de la órbita (con incertidumbres de $\sim 50$ km) son suficientes, ya que los desajustes resultantes son demasiado pequeños para afectar la explotación científica de LISA.
• Términos de Velocidad: Omitir las velocidades de las naves espaciales es la fuente principal de error de modelado para las binarias galácticas, particularmente en bajas frecuencias. Sin embargo, incluso en el peor escenario probado (SNR=200), los sesgos de parámetros resultantes permanecen por debajo de 1$\sigma$.
• Recomendación: Para la inferencia de alta fidelidad de binarias galácticas de baja frecuencia, se debe adoptar la respuesta completa con efecto Doppler, ya que no introduce un costo computacional adicional en comparación con la aproximación de orden cero.
• Trabajo Futuro: Los autores sugieren que los ceros de la respuesta dependientes de la frecuencia podrían servir como un diagnóstico de calibración y que estudios futuros deberían extender este análisis a Inspirales de Relación de Masa Extrema (EMRIs), binarias de origen estelar y fondos estocásticos.