P3.

- Ondas gravitacionales

El 11 de febrero de 2016 la colaboración aLIGO (advanced Laser Interferometer Gravitational-wave

Observatory) anunció al mundo la primera detección directa de ondas gravitacionales, predichas desde hace un
siglo por la Teoría de la Relatividad General de Einstein. La señal con el código GW150914 (figura 1) se había
detectado el 14 de septiembre de 2015 en los observatorios de Hanford y Livingston como una oscilación de
unos 40 Hz de frecuencia inicial, que aumentó hasta unos 250 Hz en menos de 20 ms. El análisis de la señal,
que llevó varios meses, permitió deducir que la onda se había producido a más de mil millones de años luz de la
Tierra durante el proceso de fusión de dos agujeros negros moviéndose en órbitas espirales, es decir de radio
decreciente, debido precisamente a la pérdida de energía por la emisión de una potente onda gravitacional
(figura 2), hasta que el sistema colapsó en un único agujero negro de masa apreciablemente inferior a la suma de
las masas de los dos agujeros negros originales.

Fig. 2 Impresión artística de las ondas

Fig. 1 Señales detectadas gravitacionales de dos agujeros negros en órbita.

Como sabrá, un agujero negro es un objeto de enorme masa que crea en su entorno un campo gravitatorio
tan intenso que ninguna partícula, incluido un fotón que viaja a la velocidad de la luz c, puede escapar a su
atracción si se encuentra a una distancia de su centro inferior a Rs (radio de Schwarzschild). En primera
aproximación, podemos imaginar el agujero negro como una esfera masiva de radio Rs .
a) Planteando que la velocidad de escape desde la distancia Rs es c, demuestre que el radio de
Schwarzschild de un agujero negro de masa M es R s  2 GM / c 2 , donde G es la constante de gravitación
universal.

Estudiaremos a continuación el movimiento orbital de dos agujeros negros
de igual masa M que interaccionan gravitatoriamente. Suponga que ambos
M O M
describen una trayectoria circular en torno al centro geométrico O del sistema
R Rs
(centro de masas), siendo R la distancia entre sus centros (figura 3).
b) Determine la velocidad angular  con que giran ambos cuerpos, en 
función de G, M y R. Fig. 3
Vamos a hacer algunos cálculos aproximados, cuando los agujeros negros están próximos a tocarse, es
decir cuando la distancia entre sus centros es un poco superior a 2Rs . Considere en concreto R  3Rs . La
frecuencia de la señal oscilante detectada en la Tierra1 en estas circunstancias fue fseñal  200 Hz (situación en
t ≈ 0,42 ms en la escala de la figura 1).
c) Obtenga una expresión para  en este caso, en función de G, M y c.
d) Calcule la masa M de cada agujero negro. Exprese su resultado en kg y en masas del Sol.

1
Ayuda: Tenga en cuenta que en cada revolución completa del sistema binario se emiten dos crestas de ondas gravitatorias,
de forma que la frecuencia de la señal detectada en la Tierra es el doble de la frecuencia de la órbita de los agujeros
negros.
 La Teoría de la Relatividad General permite determinar la potencia P emitida en forma de ondas
gravitacionales por un sistema de dos masas orbitando bajo la atracción gravitatoria mutua, como el que estamos
considerando. Suponiendo que las dos masas son iguales y que la órbita es circular, se obtiene

64 G 4 M 5
P
5 c5 R5
La órbita real es espiral, con radio decreciente, pero pueden hacerse cálculos estimativos considerando la
misma órbita circular con R  3Rs de los apartados anteriores. En la figura 1 se observa que la onda gravitatoria
se emitió principalmente durante un breve intervalo de tiempo, del orden de t  10 ms .

e) Suponiendo que toda la energía emitida durante el proceso de interacción y colapso de los agujeros
negros se traduce en una pérdida de masa del sistema, haga una estimación de esta pérdida, M . Exprese
su resultado en kg y en masas del Sol.
Vamos a hablar ahora del sistema de detección. Se ha empleado un
interferómetro de Michelson, que se esquematiza en la figura 4. Un haz de E2
luz láser, de longitud de onda , incide a 45º sobre una lámina
semiespejada (divisor de haz, D) donde se divide en dos haces, 1 y 2, que 2
viajan en direcciones perpendiculares. Cada haz se refleja normalmente en E1
un espejo plano, E1 y E2, y vuelve hacia D. Parte del haz 1 se refleja y parte  1
del 2 se transmite, de forma que en el haz 3 se superponen (interfieren) las
ondas luminosas 1 y 2 que han ido y vuelto por cada uno de los brazos del D L
interferómetro. La intensidad de la onda resultante depende de la diferencia 3
de fase, , entre estas dos ondas. Por ejemplo, si se superponen en F
contrafase () la intensidad resultante I 3 es nula (mínimo Fig. 4
interferencial). Cualquier pequeña variación en la longitud de los brazos
produce un cambio en el estado interferencial y por tanto en I 3 , que se mide con un fotodetector, F.
Siguiendo con el ejemplo anterior, si partiendo de un mínimo interferencial la longitud de uno de los
brazos aumenta en /4, el camino total recorrido por la luz en ese brazo (ida y vuelta) aumenta en /2, de forma
que las ondas pasan a interferir en fase y se tiene un máximo interferencial. Si el aumento de longitud fuese /2
se alcanzaría un nuevo mínimo nulo, en el orden interferencial siguiente al inicial.
En el interferómetro del observatorio de
Livingston los dos brazos tienen la misma Elongación
de la onda
longitud L0 = 4,0 km. En la figura 5 se
esquematiza, de forma muy exagerada, lo que t
ocurre cuando una onda gravitacional alcanza
la superficie de la Tierra. Estas ondas acortan
y alargan periódicamente la fábrica del L0A
espacio tiempo de forma que la longitud L de
los brazos del interferómetro oscila entre
L0  A y L0  A , con la particularidad de que L0+A
estas oscilaciones están en contrafase, es decir Fig. 5
cuando el brazo 1 alcanza su longitud máxima
L0  A , el brazo 2 tiene la mínima, L0  A , y viceversa.
Suele trabajarse en función de la deformación unitaria (“strain” en inglés) definida como h  L / L0 ,
donde  L es la diferencia entre las longitudes de los dos brazos. Esta magnitud adimensional h es la que aparece
en ordenadas de la gráfica de la figura 1, en la que se observan claramente las oscilaciones de h producidas por
la llegada de la onda gravitacional. La amplitud de esta oscilación llega a alcanzar el valor hmax  10 21 cuando
los agujeros negros empiezan a fusionarse.
f) Haga una estimación de la máxima amplitud de oscilación de los brazos del interferómetro de Livingston,
Amax , cuando recibió esta señal. Compare su resultado con el radio de un protón.
La deformación a detectar es en la práctica tan pequeña que han sido necesarios extraordinarios esfuerzos
técnicos para impedir que quede enmascarada por el “ruido” producido por microsismos o pequeñas variaciones
térmicas.
Para poder detectar ondas muy débiles conviene que L0 sea lo mayor E2
posible. Pero la longitud antes indicada, L0 = 4,0 km, no sería suficiente en
la práctica para poder detectar las ondas producidas por la mayoría de los 2
sucesos cósmicos previsibles. Para aumentar la longitud, se ha recurrido a E1
situar otro espejo dentro de cada brazo, entre el divisor de haz y el espejo  1
original del Michelson (figura 6). Estos espejos tiene una pequeña
transmitancia, de forma que, en promedio, la luz realiza unos N  280 D
3
viajes de ida y vuelta entre los dos espejos enfrentados antes de volver a
alcanzar el divisor de haz y contribuir a la interferencia en el haz 3. Con este F
procedimiento se consigue aumentar la longitud efectiva de los brazos del Fig. 6
Michelson hasta
Lef  N L0  1100 km .
Para terminar, vamos a caracterizar la extraordinaria sensibilidad del aparato. Observando la gráfica de la
figura 1, se deduce que el sistema es capaz de detectar deformaciones, no enmascaradas por el ruido, con una
amplitud mínima hmin  10 22 .
g) Exprese la sensibilidad del instrumento en fracciones de orden interferencial2,  min / 2 .

Datos: Constante de gravitación universal, G  6 , 6710 11 N m 2 / kg 2

Velocidad de la luz en el vacío, c  3, 0108 m / s
Masa del Sol, M Sol  2,0  1030 kg
Radio de un protón, r p  0 ,8810 15 m
Longitud de onda de la luz láser (Nd-YAG):   1, 06 μm

2
Ayuda: Cada nuevo orden interferencial corresponde a un cambio 2 en la diferencia de fase  entre las dos ondas que se
superponen en el haz 3.
 P3.- Solución

a) El radio de Schwarschild se obtiene planteando que la velocidad de escape (correspondiente a energía

mecánica nula) es c.
1 Mm 2 GM
E  0  mc 2  G → Rs 
2 Rs c2

1/ 2
M2 R  2 GM 
b) G  M 2 →    
R2 2  R3 

6GM c3
c) R  3Rs  → 
c2 108 GM

d) fseñal  200 Hz → f órbita 100 Hz ;   200  s 1

c3
M → M  6 , 2 10 31 kg  31 M Sol
108 G

Las masas deducidas del estudio detallado de la señal son M 1  ( 36  5) M Sol , M 2  ( 29  4) M Sol

2 GM
e) Rs   92 km ; R  3 R s  276 km
c2

64 G 4 M 5
P  6 , 0 10 49 W
5 c5R5

E  Pt  6,0  10 47 J  M c 2 → M  6,7  10 30 kg  3,3M Sol

La pérdida real estimada es M  ( 3 0,5) M Sol

f) Amplitud de la deformación máxima detectada en L0 = 4 km: hmax  10 21

L max 2 Amax h L
h max   → Amax  max 0
L0 L0 2

Amax  210 18 m r p / Amax  440

g) Si la diferencia de longitudes de los brazos es L, la diferencia de caminos ópticos (ida+vuelta) es 2L.
Por tanto, la diferencia de fase entre las dos ondas es
2
  k 2 L  2 h Lef

Donde k es el número de ondas.
Cada cambio de fase 2 corresponde a un nuevo orden iterferencial. La amplitud mínima detectable es
hmin  10 22 . Por tanto, la sensibilidad del instrumento, en fracciones de orden interferencial es
 min 2 Lef  min
 h → ~ 2 10 10
2  min 2
Nótese que, como se indica en el enunciado, para poder detectar señales muy débiles ( hmin lo menor
posible) interesa Lef lo mayor posible.
 P3.- Tabla de respuestas

Apartado Resultados analíticos Resultados numéricos Puntos

2 GM
a) Rs  1
c2

1/ 2
 2 GM 
b)     1,5
 R3 

c3
c)  1
108 GM

M  6,2  10 31 kg 0,8
d)
M  31M Sol 0,2

M  6,7  10 30 kg 1,8
e)
M  3, 3 M Sol 0,2

Amax  210 18 m 1,3

f)
r p / Amax  440 0,2

 min
g) ~ 2 10 10 2
2