Ce travail se situe en marge de ses voisins : à proximité, on trouve des études sur la structure à grande échelle de l'Univers, l'apprentissage automatique et les méthodes d'analyse de séries temporelles, mais aucun équivalent direct. Les plus proches voisins s'intéressent aux cartes de matière noire, à la prévision des éruptions solaires et aux horloges quantiques, c'est-à-dire qu'ils utilisent des approches par corrélation dans d'autres contextes. Cet article comble une niche méthodologique : il fournit un appareil statistique rigoureux pour l'analyse de corrélation de séries temporelles irrégulières, qui faisait jusqu'ici défaut.
L'idée clé est de remplacer les moyennes simples par des moyennes pondérées dans le temps. Pour la variance, cela semble naturel : l'intégrale sur le temps devient une somme où chaque terme est multiplié par la durée de l'intervalle correspondant. Pour l'autocorrélation, les auteurs introduisent un axe temporel centré par paires comme artifice auxiliaire, définissent le délai d'essai optimal comme la distance temporelle moyenne entre les paires de points, et ajoutent une pénalité gaussienne pour le désalignement résiduel au sein de la paire. Les poids issus de la règle trapézoïdale tiennent compte de la densité locale d'échantillonnage. Les intervalles de confiance sont construits par la méthode de Monte-Carlo : pour des temps d'observation fixes, les fluctuations de flux sont modélisées comme un bruit blanc, et sur de nombreuses réalisations, on calcule les seuils de significativité. L'erreur totale sur le délai temporel se compose de deux parties : l'erreur standard de la moyenne sur les réalisations de Monte-Carlo et l'erreur quadratique moyenne de l'irrégularité temporelle, calculée analytiquement à partir de la dispersion des différences par paires.
La méthode NUACF détecte avec succès des pics et des creux significatifs dans les courbes de lumière de la céphéide CEP J052707.40-651812.6 avec une période d'environ 10 jours, de la variable de type δ Scuti DSCT J061806.84-754759.0 avec huit composantes fréquentielles, et même de l'étoile AP2082775 avec de petites fluctuations non périodiques. Les méthodes traditionnelles (ACF rééchantillonnée et interpolée) donnent soit des profils déformés, soit ne disposent pas de critère de significativité. Pour la corrélation croisée sur les données du noyau actif de galaxie Fairall 9, la méthode NUCCF est en accord avec la CCF interpolée et l'ajustement de modèle CREAM dans 7 des 13 bandes, et dans la bande W1, elle donne un délai de signe opposé avec une significativité élevée. Pour le quasar à lentille HE 0435-1223, NUCCF donne un pic large autour d'un délai nul ; la largeur du pic significatif détermine la plage dans laquelle se situent les estimations des méthodes d'ajustement spécialisées, ce qui sert de validation indépendante à un niveau >4σ.
Contrairement aux travaux voisins qui appliquent des méthodes de corrélation à des problèmes astrophysiques spécifiques, cet article propose un outil universel, applicable à toute série temporelle irrégulière. Il fournit pour la première fois une formule analytique de la fonction de corrélation qui se réduit à la formule standard pour un échantillonnage régulier, et l'accompagne d'intervalles de confiance et d'une erreur totale sur le délai temporel. Cela comble une lacune de longue date : auparavant, les chercheurs étaient contraints soit d'interpoler les données, introduisant des biais incontrôlés, soit d'utiliser un regroupement en classes, perdant en résolution et sans estimation de la significativité.
La méthode pourrait être étendue aux corrélations multidimensionnelles, à l'analyse des sursauts radio rapides (FRB) pour étudier le lien entre la luminosité de la source radio persistante et l'activité des sursauts répétitifs, ainsi qu'à l'analyse conjointe de courbes de lumière dans différentes bandes de longueurs d'onde. Une adaptation est possible pour des données avec des lacunes et pour un échantillonnage irrégulier en coordonnées spatiales, par exemple pour les fonctions de corrélation de galaxies.
La méthode aura un impact sur l'analyse des données des grands relevés du ciel, tels que LSST, ZTF, ASAS-SN, Gaia, où l'irrégularité de l'échantillonnage est inévitable. Elle est également utile pour le reverberation mapping des noyaux actifs de galaxies et la mesure des délais temporels des quasars à lentille gravitationnelle.
Les auteurs prévoient d'appliquer la méthode aux données des sursauts radio rapides et à une classe plus large d'étoiles variables. Ils envisagent également de créer un logiciel public pour une analyse automatique.
Cet article résout le problème de l'absence d'un appareil statistique rigoureux pour l'analyse de corrélation de séries temporelles irrégulières en astrophysique. Cela est directement lié à des questions non résolues : comment distinguer une périodicité réelle des artefacts d'échantillonnage dans les étoiles variables, comment mesurer les délais temporels en reverberation mapping sans erreurs systématiques, et comment valider les délais des quasars à lentille sans hypothèses de modèle.
La méthode NUACF, appliquée de manière itérative à elle-même, agit comme un filtre passe-bas : dès la troisième itération, même un bruit pur acquiert une structure régulière. C'est pourquoi les auteurs recommandent de se limiter au deuxième ordre d'itérations pour ne pas créer de fausses périodicités.