La orquesta toca sin director. Cada instrumento entra según su propio calendario: el violín el lunes, la flauta el miércoles, los timbales el sábado. Los astrónomos no oyen 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. Resulta una cuadrí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 no es así. Hay que reconstruir los puntos faltantes (interpolación) o agrupar las observaciones en intervalos de tiempo (binning). Lo primero añade datos inventados; lo segundo difumina los detalles. Y nadie puede decir cuán significativo es el pico obtenido. Es como intentar oír el ritmo de una orquesta reordenando las notas en una línea uniforme.
Esto se parece a 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 por sí mismos y adapta su batuta a sus retardos naturales. Para estimar los errores, los autores ejecutan el método de Monte Carlo: con 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, NUACF destacó una 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 del modelo CREAM en 7 de las 13 bandas, y en la banda W1 da un desfase de signo opuesto con alta significancia.
Hasta esta investigación no existía una forma rigurosa de calcular funciones de correlación para observaciones con distribución temporal no uniforme. Ahora contamos con un método analítico que tiene en cuenta 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 de 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 por 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.
La orquesta sin director, por fin, se ha escuchado. El método NUACF/NUCCF no intenta obligar a los músicos a tocar sobre una cuadrí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, incluida la incertidumbre del propio calendario de observaciones.
El método NUACF, aplicado iterativamente sobre sí mismo, actúa como un filtro de paso bajo: ya en la tercera iteración, incluso el ruido puro adquiere una estructura regular. Por eso los autores recomiendan limitarse al segundo orden de iteraciones para no crear periodicidades falsas.