Adiabatic tides in compact binaries on quasi-elliptic orbits: Radiation at the second-and-a-half relative post-Newtonian order

Resumen Técnico: Mareas Adiabáticas en Binarias Compactas con Órbitas Cuasi-elípticas

Título: Adiabatic tides in compact binaries on quasi-elliptic orbits: Radiation at the second-and-a-half relative post-Newtonian order
Autor: Quentin Henry (Universitat de les Illes Balears)

1. El Problema

La detección de ondas gravitacionales (OG) por el consorcio LIGO-Virgo-KAGRA (LVK) ha revelado eventos de coalescencia de sistemas binarios que contienen estrellas de neutrones (NS) y agujeros negros (BH). Un caso particularmente relevante es el evento GW200105, un sistema NS-BH que entró en la banda del detector con un movimiento excéntrico significativo (eccentricidad estimada $\sim 0.13$ a 20 Hz).

El problema central abordado en este trabajo es la incompletitud de los modelos de onda actuales para tales sistemas. La mayoría de los modelos de waveform utilizados para la estimación de parámetros (como Phenom o EOB) asumen órbitas cuasi-circulares o no incluyen efectos de tamaño finito (mareas) cuando la órbita es excéntrica. Aunque existen modelos que combinan mareas y excentricidad (como TEOBResumS-Dalí), carecen de información rigurosa de la teoría post-newtoniana (PN) sobre los términos de marea en la fuerza de reacción radiativa para órbitas no circulares. Esto puede introducir sesgos sistemáticos en la estimación de parámetros físicos y en la detección de desfasajes (dephasing) acumulados.

2. Metodología

El autor emplea el marco de la teoría post-newtoniana (PN) combinada con la acción efectiva de campo (EFT) para modelar las interacciones de marea. El trabajo se realiza en el orden relativo 2.5PN (que incluye efectos de disipación de energía y momento angular).

• Acción de la Materia: Se utiliza una acción efectiva que incluye interacciones de marea adiabática (cuadrupolo y octupolo de masa, y cuadrupolo de corriente) hasta el orden NNLO (Next-to-Next-to-Leading Order).
• Dinámica Conservativa: Se basa en los resultados del "Paper I" (artículo complementario), donde se resolvieron las ecuaciones de movimiento utilizando una parametrización cuasi-Kepleriana (QKP) para órbitas excéntricas, obteniendo las cantidades conservadas y la evolución secular de los elementos orbitales.
• Formalismo PN-MPM: Se utiliza el formalismo Post-Newtoniano-Multipolar-Post-Minkowskiano para derivar los momentos multipolares radiativos ($U_L, V_L$) a partir de los momentos de la fuente.
• Cálculo de Flujos: Se calculan los flujos de energía y momento angular radiados, separando las contribuciones en:
  • Instantáneas: Dependientes del estado actual del sistema.
  • Cola (Tail): Efectos hereditarios debidos a la dispersión de las ondas en el campo gravitatorio estático (integración sobre el pasado).
  • Memoria: Contribuciones no lineales (dejadas para trabajos futuros).
• Expansión y Resumación: Los resultados se expanden en serie de potencias de la excentricidad temporal ($e_t$) hasta el orden 12 ($O(e_t^{12})$) para las amplitudes y hasta el orden 14 para los flujos. Se aplica una técnica de resumación para los términos de cola, mejorando la precisión para excentricidades altas.
• Evolución Secular: Se integran numéricamente las ecuaciones de balance para obtener la evolución secular de los elementos orbitales ($x, e_t$) y se calcula la amplitud de la deformación (strain) descompuesta en armónicos esféricos de peso espín ($h_{\ell m}$).

3. Contribuciones Clave

• Derivación Rigurosa de Mareas Excéntricas: Es el primer trabajo que trata rigurosamente las correcciones de marea adiabática en el sector radiativo para órbitas binarias no cuasi-circulares hasta el orden 2.5PN.
• Flujos de Energía y Momento Angular: Se proporcionan expresiones analíticas completas para los flujos orbitales promediados, incluyendo términos instantáneos exactos en excentricidad y términos de cola resumidos.
• Modos de Onda Gravitacional: Se derivan los modos de amplitud ($h_{\ell m}$) descompuestos en armónicos esféricos, incluyendo correcciones post-adiabáticas (PA) necesarias para los modos $(2,2)$ y $(2,0)$.
• Expansión de Alto Orden: Se ofrecen resultados expandidos hasta $O(e_t^{12})$ para la amplitud y $O(e_t^{14})$ para los flujos, superando las limitaciones de modelos anteriores que solo llegaban a $O(e_t^6)$.
• Archivo de Datos: Todos los resultados analíticos complejos se han compilado en un archivo auxiliar (notebook de Mathematica) para su uso en modelos fenomenológicos (Phenom) y de cuerpo único efectivo (EOB).

4. Resultados Principales

• Evolución Orbital: Se obtuvieron las ecuaciones diferenciales acopladas para la evolución secular del parámetro de frecuencia ($x$) y la excentricidad ($e_t$) incluyendo correcciones de marea.
• Desfasaje (Dephasing): Mediante una integración numérica de sistemas binarios (incluyendo el caso GW200105 y binarias de estrellas de neutrones hipotéticas), se encontró que:
  • Las correcciones de marea excéntrica inducen un desfasaje en la fase de la onda gravitacional.
  • Para sistemas con alta excentricidad inicial o relaciones de masa altas (como binarias de estrellas de neutrones), este desfasaje puede ser detectable en los futuros observatorios de tercera generación (como el Einstein Telescope) o en las rondas de observación O5 de LIGO.
  • El efecto es pequeño para el caso específico de GW200105 debido a su excentricidad moderada, pero significativo en regiones específicas del espacio de parámetros.
• Precisión de los Modos: La expansión de alto orden en excentricidad evita características no físicas en la waveform para excentricidades $e > 0.3$, mejorando la fidelidad de los modelos para sistemas con órbitas muy excéntricas.

5. Significado e Impacto

• Mejora de Modelos de Waveform: Este trabajo proporciona la información teórica necesaria para actualizar los modelos de waveform (como Phenom y EOB) utilizados por el LVK, permitiendo una descripción más precisa de sistemas binarios con mareas y excentricidad.
• Estimación de Parámetros: La inclusión de estos términos es crucial para evitar sesgos sistemáticos en la inferencia de parámetros (masas, espines, parámetros de marea/Love numbers) de eventos futuros.
• Física de Estrellas de Neutrones: Al permitir una mejor caracterización de las mareas en órbitas excéntricas, se abre la puerta a medir con mayor precisión la ecuación de estado de la materia nuclear en estrellas de neutrones.
• Fundamento para Futuros Trabajos: Este estudio sienta las bases para el cálculo de efectos de marea dinámica y términos de memoria, esenciales para la modelización completa de sistemas compactos en órbitas excéntricas.

En resumen, el artículo cierra una brecha teórica importante al proporcionar las herramientas analíticas necesarias para modelar la radiación gravitacional de sistemas binarios compactos con mareas y órbitas excéntricas, mejorando la capacidad de los detectores actuales y futuros para explorar la física de la materia densa y la gravedad fuerte.