Me das vueltas(cómo una estrella): En búsqueda de ondas gravitacionales continuas de estrellas de neutrones binarias desconocidas

Las estrellas de neutrones (NS, del inglés Neutron Stars) son el resultado de explosiones de supernova, el remanente de una estrella primitiva con una masa de entre 10 y 25 veces la masa del Sol. Estos objetos compactos tienen un radio típico de 10km y una masa similar a la de nuestro Sol, haciendo de las NS uno de los entornos más extremos donde se ha observado materia. Actualmente, la estructura y composición de estos objetos compactos es un campo de investigación activo e interdisciplinar, que une trabajos de las comunidades de física de partículas y astrofísica. La extrema compacidad (o densidad) de estos objetos hace de ellos un interesante laboratorio para poner a prueba efectos relativistas como la emisión de radiación gravitacional. De hecho, los detectores Advanced LIGO y Advanced Virgo han detectado numerosos eventos de ondas gravitacionales relacionados con la fusión de dos NS. Una de esas detecciones, GW170817, incluso estuvo acompañada por un homólogo electromagnético, generando la primera detección multimensajera de un evento astrofísico mediante ondas gravitacionales y luz.

Las ondas gravitacionales continuas (CW, del inglés Continuous Waves) son otro canal mediante el cual se puede investigar la estructura interna de NS en rotación rápida. Según la teoría, las NS pueden sostener desviaciones respecto a una forma simétrica, bien sea por las imperfecciones en la corteza exterior, una perturbación oscilante en la estructura interna o precesiones libres debidas a desalineaciones entre la simetría y los ejes de rotación. A medida que la NS gira rápidamente, dichas deformaciones producirían radiación gravitacional en forma de CW. Esta forma de radiación gravitacional es varios órdenes de magnitud más débil que las producidas durante la fusión de objetos compactos. Aun así, duran largos periodos de tiempo (de meses a años), permitiendo la integración de largosflujos de datos para acumular una alta relación señal/ruido.

Las búsquedas para este tipo de señal están divididas acorde con lo que se conoce sobre las potenciales fuentes. Por ejemplo, las búsquedas específicas (targeted searches en inglés) se centran en NS para las cuales la posición en el cielo y la frecuencia de rotación son conocidas por medios electromagnéticos, mientras que las búsquedas dirigidas apuntan a ciertas regiones del cielo donde pueden estar situadas NS desconocidas. Nosotros presentamos los resultados de una búsqueda del cielo completo, que significa que buscamos señales de CW provenientes de NS completamente desconocidas desde cualquier dirección del cielo. Particularmente, nos centramos en NS desconocidas en sistemas binarios. La señal esperada de NS aisladas es una onda de larga duración, cuya frecuencia se reduce lentamente debido a la emisión de energía a través de diferentes medios como la radiación electromagnética o gravitatoria. Este efecto de ralentización de la rotación (spindown en inglés) es suficientemente lento como para ser ignorado dada la población de NS consideradas por esta búsqueda. Desde el punto de vista del detector en la Tierra, la señal presenta una modulación de frecuencia debido a la rotación diaria y los movimientos orbitales de nuestro planeta. Esta señal es mucho más compleja para NS en un sistema binario, ya que se tiene que considerar una segunda modulación de Doppler a causa del movimiento relativo de la fuente alrededor de su cuerpo acompañante.

En esta búsqueda, aplicamos una variación de la transformada de Hough llamada BinarySkyHough, que opera utilizando un espectrograma de los datos. Este espectrograma aporta información sobre qué frecuencias son las más destacadas en los datos a medida que la observación progresa. La idea básica es que una señal CW podría aparecer como un exceso de potencia en el espectrograma que modula de una forma muy particular, como se ejemplifica en la Figura 1. Podemos describir trazas significativas usando un conjunto de parámetros relacionados con las propiedades físicas de la fuente considerada, como su frecuencia de rotación o su posición en el cielo; por tanto, buscar una CW es cuestión de identificar trazassignificativas en los datos. Debido a la inclusión de parámetros orbitales binarios, que describen la forma de una órbita binaria, tenemos que buscar más parámetros que cuando exploramos NS aisladas. Esto plantea un problema para los algoritmos de búsqueda típicos, a medida que el coste computacional, fácilmente, se vuelve insoportable. Una solución sencilla es usar Unidades de Procesamiento Gráfico (GPU) para realizar este tipo de búsquedas eficientemente, ya que son capaces de agilizar cálculossimplesllevando a cabo una gran cantidad de ellos en paralelo.

Usamos datos del principio del tercer periodo de observación de los detectores Advanced LIGO y Advanced Virgo, abarcando seis meses de datos desde abril hasta septiembre del 2019, para buscar CW de NS desconocidas en sistemas binarios desconocidos en la banda más sensible de los detectores. No encontramos evidencias de señales de CW. Como resultado, estimamos la sensibilidad de nuestra búsqueda mediante el estudio de la recuperación de una población de señales simuladas. Inicialmente esta sensibilidad está expresada en términos de la amplitud de CW más débil detectable por nuestra tubería. Conseguimos la mejor sensibilidad hasta la fecha para el espacio de parámetros analizado, con la amplitud detectable más débil un 60% menor que las estimaciones anteriores. Interpretamos estos resultados en términos de responder a dos preguntas astrofísicamente interesantes, concretamente, qué tan lejos de la Tierra fue capaz de sondear nuestra búsqueda y cuál es la deformación máxima permitida para una fuente de CW dentro del rango probado por nuestra búsqueda. Las figuras 2 y 3 resumen el alcance máximo de nuestra búsqueda y la máxima deformación permitida para una fuente de CW como una función de la frecuencia de la onda gravitacional. Como se ha comentado anteriormente, consideramos fuentes de CW con spindown negligible. Esta suposición establece un límite en la máxima deformación sondeada para una NS, representado por las regiones sombreadas en dichas figuras. En otras palabras, si no existen mecanismos de equilibrio, las regiones sombreadas se excluyen de nuestros resultados, ya que implicarían tasas de spindown más altas que las cubiertas por esta búsqueda. Un posible mecanismo de equilibrio que contrarresta el spindown de la NS podría ser la acreción de materia del cuerpo acompañante. Estos resultados muestran la mejora creciente de los detectores Advanced LIGO y Advanced Virgo a medida que alcanzan su sensibilidad de diseño, permitiéndonos acercarnos a los límites superiores de las propiedades de la NS cada vez más cerca de los resultados estimados proporcionados por las simulaciones numéricas.

Link para más información