Lente Gravitacional de Ondas Gravitacionales
en un Sistema Triple Jerárquico

Tratamiento en óptica de ondas — resultados
Resumen

Calculamos el efecto de lente que un tercer cuerpo del mismo sistema triple jerárquico ejerce, en óptica de ondas, sobre la onda gravitacional (GW) de una binaria compacta — el mismo triple en dos épocas de su inspiral: un chirp tardío (lente congelado) y una fase temprana cuasi-monocromática (el lente se mueve, y produce pulsos de lensing repetidos).

Esta presentación es la versión visual; report.pdf tiene el mismo contenido con más detalle en texto, y theory/theory.pdf la derivación completa. Todo viene de cases/, verificado en tests/ — nada acá está sin chequear.

Contenidos
  1. 1El sistema Qué es un triple jerárquico, y por qué éste: el mismo inspiral mirado en dos épocas que caen a lados opuestos de la transición entre difracción y óptica geométrica.
  2. 2Caso A — el chirp, lente estático La binaria barre de 10 Hz al merger en 15 s. El lente no llega a moverse: franjas de interferencia en frecuencia, y un eco en el tiempo.
  3. 3Caso B — cuasi-monocromático, lente en movimiento La misma binaria, mucho antes, a 0.05 Hz durante 12 días. Ahora el que se mueve es el lente, y el patrón se recorre en el tiempo.
  4. 4Lo que muestran juntos Un único factor de amplificación F(w,y), recorrido por dos ejes distintos, y el puente entre óptica de ondas y óptica geométrica.
  5. 5Alcance y limitaciones Qué se modeló, qué no, y cuánto vale cada aproximación — declarado, no dejado implícito.
1El sistema

Un triple jerárquico, en dos épocas del mismo inspiral

Un triple jerárquico: una binaria interna (la fuente de GWs) orbita un tercer cuerpo mucho más masivo (el lente), lo bastante lejos como para que la marea sea despreciable. Seguimos el mismo sistema en dos momentos de la vida de la binaria interna: el chirp final (Caso A) y una etapa muy anterior, cuasi-monocromática (Caso B).

Triple jerárquico, visto de costado

Parámetros

Binaria interna
(la fuente)
m1 = 20 M☉ m2 = 15 M☉ ℳ = 15.05 M☉ (masa de chirp)
Lente ML = 5 × 104 M☉ DL = 5 kpc
Órbita externa aout = 1.82 UA Pout = 4 d iout = 87° (casi de canto, circular)
cargando…
2Caso A

El chirp, lente estático

La binaria interna hace inspiral de 10 Hz a 125.6 Hz en 15.2 s — muchísimo más rápido que la órbita externa (4 días), así que el lente queda congelado en un único parámetro de impacto (yA=1.589) durante toda la señal. La forma de onda es IMRPhenomD real (inspiral+merger+ringdown), no una aproximación truncada.

Ese chirp no cae en cualquier punto de la órbita externa: lo ubicamos en la conjunción — la fase 0.25 del período, cuando la fuente pasa exactamente por detrás del lente. De ahí sale yA=1.589, que es el mínimo de toda la órbita, el mismo valor al que llega el Caso B en el pico de cada pulso. La elección importa y es sostenible: los 15 s del chirp son 4×10−5 del período externo, así que la binaria coalesce entera sin salir de esa fase. En cualquier otra, el lente quedaría más lejos de la línea de visión y el efecto sería menor.

Caso A

La señal, con lente y sin lente

En frecuencia, lo que el lente hace es una multiplicación: hL(f)=F(f)·h(f). En el tiempo eso no es un reescaleo — es la señal original más una copia suya, retrasada y más débil. Acá están las dos curvas sobre el mismo reloj que el esquema de la coalescencia, y la diferencia entre ellas es del orden del 5%: chica, pero no constante, porque F oscila con f mientras el chirp barre la banda.

Los dos picos llegan en el mismo instante. No es casualidad: F está referida a la llegada de la imagen fuerte, así que no queda ningún corrimiento global arbitrario y cualquier desfasaje que se vea entre las curvas es físico.

La binaria interna, cayendo

h(t): ventana de 0.35 s, el mismo reloj que el esquema

cargando…
Caso A

La segunda imagen llega tarde y débil

En óptica geométrica un lente puntual da dos imágenes, en x±=12(y±y2+4): una es un mínimo del tiempo de llegada y la otra un punto de silla. La del silla recorre más camino, así que llega ΔT = 3.434 s más tarde, y está menos magnificada.

Eso es el eco: la misma forma de onda otra vez, con una amplitud de 0.2395 del pico sin lente — que es exactamente la magnificación |μ₋| = 0.2396 de esa imagen, al 0.05%. Un número que sale del cálculo en el dominio del tiempo y coincide con otro que sale de la óptica geométrica, sin nada ajustado en el medio. La segunda imagen entra en banda 11.7 s antes del merger, así que durante casi toda la señal las dos están presentes y se interfieren; el eco se ve separado recién después del merger, cuando la primera ya se apagó.

Dos caminos hasta el detector

Amplitud de la señal medida

cargando…
Caso A

Dos imágenes interfieren: F recorre una circunferencia

El lente está quieto, pero f no, y F depende de f. Al barrer la banda, F da una vuelta completa por cada franja — 397 vueltas entre 10 Hz y el merger.

Que el recorrido sea una circunferencia es la suma de las dos imágenes: un vector fijo de largo |μ₊| = 1.0283, la imagen fuerte, más otro de largo |μ₋| = 0.2396 que gira con la fase w·ΔT. Como w·ΔT es lineal en f, el período de franja es 1/ΔT = 0.2912 Hz, el mismo en toda la banda.

Y como el círculo no encierra al origen, arg F nunca completa una vuelta: las dos tangentes desde el origen la dejan acotada en |arg⁡F|≤arcsen⁡(|μ₋|/|μ₊|)=0.235 rad. Ahí está la asimetría del efecto. El lente mueve la amplitud entre 0.789 y 1.268 — ±0.24 alrededor de 1.028 — y la fase, apenas un cuarto de radián.

F en el plano complejo

La fase, arg F(f)

cargando…
Caso A

El patrón se aprieta mientras el chirp sube

La fuente está quieta en yA=1.589; lo que cambia es la frecuencia, y con ella el espaciado de los anillos — 12.5× más finos entre el inicio de la banda y el merger. La fuente los cruza cada vez más rápido, y eso es la modulación que el lente le imprime a la amplitud: sobre el final pasan decenas de franjas por segundo, así que en la señal del merger ya no se alcanza a ver.

Ampliación sobre la fuente

El patrón completo

La modulación que ve la fuente

cargando…
Caso A

Lo que el lente le hace a la onda

Amplificación de la muestra pico: 1.0557. Amplificación RMS: 1.0557, que por Parseval es A(yA). Coinciden: la amplitud sola no delata nada. Lo que hace que la forma de onda con lente sea una copia desfasada y no simplemente reescalada es arg F, que oscila a lo largo de la banda dentro de ±0.235 rad.

El pico con lente llega exactamente al mismo tiempo que el pico sin lente (F está referida a la llegada de la imagen fuerte, sin desfasaje absoluto arbitrario). Una segunda copia, ~4× más débil — la amplitud exacta coincide con la magnificación geométrica de la imagen débil, |μ₋|, a menos del 1% — de todo el merger reaparece más tarde, exactamente en el retraso fijo que predice la óptica geométrica entre las dos imágenes. Ambas construidas sólo del pasado de la fuente: la causalidad nunca estuvo en duda. Un eco genuino, visible en caseA_spectrogram.png como una segunda traza en el plano tiempo-frecuencia.

3Caso B

Cuasi-monocromático, lente en movimiento

Más temprano en el mismo inspiral: f=0.05 Hz, casi constante durante una observación de 12 días (3 períodos externos). Acá el lente sí se mueve — el parámetro de impacto barre el mismo tipo de patrón de difracción una vez por período orbital externo, y produce pulsos de amplificación repetidos.

Caso B

La órbita recorriendo el patrón

El mismo reloj para todo: dónde está la fuente, dónde cae eso en el patrón, qué mide el detector, y el retardo de Roemer que la misma órbita escribe encima — la fuente se aleja y se acerca, y su luz tarda hasta 905 s más o menos en llegar. Acá se ven dos de las tres vueltas de la observación, para poder mirarlas con detalle; las figuras estáticas muestran los 12 días enteros.

Dónde está la fuente

Dónde cae en el patrón (wB=0.31)

Patrón de difracción

Las dos huellas de la órbita

La onda, a la amplificación del instante

cargando…
Caso B

En media órbita no hay lente en absoluto

En exactamente la mitad de la órbita no hay lente: cuando la fuente está delante del lente y no detrás, no hay ninguna configuración de lente que resolver — no es una aproximación, es F=1 por construcción, y se ve en la diapositiva anterior como el marcador saliendo del patrón (cian/gris) durante la mitad "de adelante" de cada vuelta.

Esto es exactamente el origen de los pulsos periódicos: "lensing repetido", el análogo en óptica de ondas de lo que D'Orazio & Loeb (2020) encontraron en óptica geométrica.

Caso B

La órbita hace dos cosas: pulsos de lente y fase cinemática

Todo lo anterior es lo que la órbita externa hace ópticamente. Pero la misma órbita también arrastra a la fuente a lo largo de la línea de visión, a βlos=0.0165, y el retardo de Roemer que eso produce es de 1810 s pico a pico: a fB, unos 90.5 ciclos de GW de fase escritos sobre la forma de onda. La fase que escribe el lente mismo, arg F/2π, vale 0.017 ciclos. Un factor 5206×.

Retardo y corrimiento no son lo mismo, y no ocurren en la misma fase orbital. Son la misma cantidad separada por una derivada: el retardo va como la separación a lo largo de la línea de visión, el corrimiento como su velocidad. En la conjunción la velocidad radial se anula — la fuente está en ese instante lo más lejos que llega — así que ahí el retardo está en un extremo (+905 s de ±905) y el corrimiento pasa por cero (|vlos|/c = 1.3×10⁻⁶ contra un máximo orbital de 1.6×10⁻²). Los extremos del Doppler caen un cuarto de período antes y después, donde no hay lensing en absoluto. La cinemática domina la fase y deja casi intacta la envolvente; el lente produce el pulso de amplitud. Esa diferencia permite separar ambos efectos.

Roemer/Doppler vs lensing, en ciclos de GW
Las dos huellas de la misma órbita externa, una por panel y cada una a su propia escala: un cociente de 5206 no se puede dibujar en un eje común. Arriba, el retardo de Roemer en ciclos de GW junto al corrimiento Doppler en eje gemelo — se cruzan como debe ser, el verde pasa por cero justo donde el violeta pica, sobre la punteada de la conjunción, y tiene sus propios extremos un cuarto de período más allá. Abajo, la fase del lente sobre su propia grilla, con el sistemático de evaluar F a w fijo dibujado como banda: las dos excursiones grandes valen 21 veces esa banda y son un resultado, el ringing lejano cae adentro y no lo es.
Caso B

Y sin embargo no cambia ningún resultado de amplitud

No es un defecto del cálculo de lensing — es cómo se ve realmente un triple jerárquico, y es por eso que D'Orazio & Loeb (2020), el mismo trabajo que este caso sigue para el lensing repetido, modelan el Doppler junto con él (su Apéndice A).

El Doppler se aplica evaluando la forma de onda en el tiempo de emisión retardado, y por ser puro timing deja la envolvente intacta a 1.1×10−5: todos los resultados de amplitud quedan sin cambio, de manera demostrada.

4Lo que muestran juntos

Un solo F(w,y), recorrido por dos ejes distintos

5Alcance y limitaciones

Qué se modeló, qué no, y cuánto vale cada aproximación

  • Sólo lente de masa puntual, sin distribución extendida.
  • IMRPhenomD es un modelo calibrado (no re-derivado acá); el generador de respaldo TaylorF2 (si pycbc no está disponible) no tiene merger ni ringdown.
  • Sin calibración absoluta de strain — sólo amplitudes relativas y el factor de amplificación son afirmaciones físicas acá.
  • Detector idealizado: sin ruido, patrón de antena, ni filtrado adaptado.
  • La aproximación de lente delgado/paraxial se empuja cerca de cada cruce del Caso B — cuantificado (ventana de ~60 s de 4 días) y mostrado que no afecta nada reportado.
  • La cinemática de la fuente se trata a primer orden en βlos. Lo que queda más allá (corrimiento gravitacional + Doppler transverso) es exactamente constante para una órbita circular, 4.08×10−4, y por lo tanto exactamente degenerado con la masa de chirp: se reporta como sesgo sobre 𝓜, no se aplica.
  • El Doppler entra en la forma de onda pero no en F, que se evalúa al wB nominal fijo. Cuantificado: 1.0×10−4 de error en |F| en el pulso, hasta 3.5×10−3 en el peor punto de la mitad con lente.

Para seguir

Derivación completa: theory/theory.pdf
Versión en PDF de este resultado: report.pdf
Reproducir todo: README.md → ./reproduce.sh