Оркестр играет без дирижёра. Каждый инструмент вступает по своему расписанию: скрипка в понедельник, флейта в среду, литавры в субботу. Астрономы слышат не музыку, а отдельные ноты, разбросанные по календарю. Телескоп видит звезду только когда она над горизонтом и когда расписание выделило ей время. Получается неровная временная сетка, на которой классические методы корреляции спотыкаются.
Стандартная автокорреляционная функция требует, чтобы измерения шли через равные промежутки. В реальности так не бывает. Приходится либо достраивать пропущенные точки (интерполяция), либо складывать наблюдения в корзины по времени (бинирование). Первое добавляет выдуманные данные, второе размывает детали. И никто не может сказать, насколько полученный пик значим. Это как пытаться услышать ритм оркестра, переставляя ноты в ровную строчку.
Это похоже на дирижёра, который не заставляет всех играть строго по метроному, а слушает, как музыканты сами держат темп, и подстраивает палочку под их естественные задержки. Для оценки ошибок авторы запускают метод Монте-Карло: при фиксированных временах наблюдений флуктуации потока моделируются как белый шум, и по множеству реализаций вычисляются границы значимости. Полная ошибка временной задержки складывается из двух частей: стандартной ошибки среднего по реализациям Монте-Карло и среднеквадратичной ошибки временной нерегулярности, вычисляемой аналитически из разброса парных разностей.
Метод проверен на реальных кривых блеска. Цефеида CEP J052707.40-651812.6 с периодом около 10 дней показала значимые пики и впадины. Переменная типа δ Щита DSCT J061806.84-754759.0 с восемью частотными компонентами тоже проявила свою структуру. Даже у звезды AP2082775 с малыми непериодическими флуктуациями NUACF выделил регулярность там, где традиционные методы давали лишь шум. Для кросс-корреляции на данных активного ядра галактики Fairall 9 метод NUCCF согласуется с интерполированной CCF и модельной подгонкой CREAM в 7 из 13 полос, а в полосе W1 даёт задержку противоположного знака с высокой значимостью.
Для линзированного квазара HE 0435-1223 NUCCF даёт широкий пик около нулевой задержки. Ширина значимого пика определяет диапазон, внутри которого лежат оценки специализированных методов подгонки. Это независимая валидация на уровне >4σ.
До этого исследования не было строгого способа вычислять корреляционные функции для неравномерно распределённых во времени наблюдений. Теперь появился аналитический метод, который напрямую учитывает неравномерность выборки, даёт доверительные интервалы и полную ошибку временных задержек, включая вклад нерегулярности расписания наблюдений. Заменили простые средние на взвешенные по времени, ввели штраф за рассогласование пар и построили доверительные интервалы методом Монте-Карло при фиксированных временах наблюдений. Метод может применяться для анализа кривых блеска переменных звёзд, реверберационного картирования активных ядер галактик, измерения задержек линзированных квазаров и изучения быстрых радиовсплесков.
Оркестр без дирижёра наконец-то услышан. Метод NUACF/NUCCF не пытается заставить музыкантов играть по ровной сетке, а ловит их естественный ритм. И этот ритм теперь можно не только услышать, но и измерить с полной ошибкой, включающей неопределённость самого расписания наблюдений.
Метод NUACF при итеративном применении к самому себе действует как фильтр нижних частот: уже на третьей итерации даже чистый шум приобретает регулярную структуру. Поэтому авторы рекомендуют ограничиваться вторым порядком итераций, чтобы не создавать ложные периодичности.