Cómo oír el ritmo de una estrella cuando toca fuera de compás ✎

Una orquesta sin director. Cada instrumento entra según su propio horario: el violín los lunes, la flauta los miércoles, los timbales los sábados. Los astrónomos no oyen la música, sino notas sueltas dispersas en el calendario. El telescopio ve la estrella solo cuando está sobre el horizonte y cuando el calendario le ha reservado tiempo. El resultado es una retícula temporal irregular en la que los métodos clásicos de correlación tropiezan.

La función de autocorrelación estándar exige que las mediciones se realicen a intervalos iguales. En la realidad eso no ocurre. Hay que elegir entre reconstruir los puntos faltantes (interpolación) o agrupar las observaciones en intervalos de tiempo (binning). La primera opción añade datos inventados; la segunda difumina los detalles. Y nadie puede decir qué significancia tiene el pico obtenido. Es como intentar oír el ritmo de una orquesta reordenando las notas en una fila uniforme.

Es como un director que no obliga a todos a tocar estrictamente con el metrónomo, sino que escucha cómo los músicos mantienen el tempo y adapta su batuta a sus retardos naturales. Para estimar los errores, los autores ejecutan el método de Monte Carlo: con los tiempos de observación fijos, las fluctuaciones del flujo se modelan como ruido blanco y, a partir de múltiples realizaciones, se calculan los límites de significancia. El error total del desfase temporal se compone de dos partes: el error estándar de la media de las realizaciones de Monte Carlo y el error cuadrático medio de la irregularidad temporal, calculado analíticamente a partir de la dispersión de las diferencias por pares.

El método se ha probado con curvas de luz reales. La cefeida CEP J052707.40-651812.6, con un período de unos 10 días, mostró picos y valles significativos. La variable de tipo δ Scuti DSCT J061806.84-754759.0, con ocho componentes de frecuencia, también reveló su estructura. Incluso en la estrella AP2082775, con pequeñas fluctuaciones no periódicas, el NUACF detectó regularidad donde los métodos tradicionales solo daban ruido. Para la correlación cruzada con los datos del núcleo activo de galaxia Fairall 9, el método NUCCF concuerda con la CCF interpolada y con el ajuste de modelo CREAM en 7 de las 13 bandas, y en la banda W1 produce un desfase de signo opuesto con alta significancia.

Para el cuásar con lente HE 0435-1223, el NUCCF da un pico ancho alrededor del desfase cero. La anchura del pico significativo determina el rango dentro del cual se encuentran las estimaciones de los métodos especializados de ajuste. Esto constituye una validación independiente a un nivel >4σ.

Antes de este estudio no existía una forma rigurosa de calcular funciones de correlación para observaciones con distribución temporal desigual. Ahora disponemos de un método analítico que incorpora directamente la irregularidad del muestreo, proporciona intervalos de confianza y un error completo de los desfases temporales, incluida la contribución de la irregularidad del calendario de observaciones. Se han sustituido los promedios simples por promedios ponderados en el tiempo, se ha introducido una penalización por desajuste entre pares y se han construido intervalos de confianza mediante el método de Monte Carlo con tiempos de observación fijos. El método puede aplicarse al análisis de curvas de luz de estrellas variables, al mapeo de reverberación de núcleos activos de galaxias, a la medición de desfases en cuásares con lente y al estudio de ráfagas rápidas de radio.

Por fin se ha oído a la orquesta sin director. El método NUACF/NUCCF no intenta obligar a los músicos a tocar sobre una retícula uniforme, sino que capta su ritmo natural. Y ese ritmo ahora no solo puede oírse, sino también medirse con un error completo que incluye la incertidumbre del propio calendario de observaciones.