Apparent criticality from common-drive aliasing in self-exciting networks

Résumé Technique : Criticité Apparente par Repliement de Commande Commune dans les Réseaux Auto-excitateurs

Énoncé du Problème
Dans la physique statistique des systèmes en cascade, la criticité est souvent diagnostiquée via un paramètre de contrôle scalaire : le nombre de branchement ou de reproduction effectif. Dans les réseaux auto-excitateurs modélisés par des processus de Hawkes multivariés, ce paramètre correspond au rayon de Perron-Frobenius ($\rho$) de la matrice de reproduction. Si $\rho < 1$, le système est sous-critique (l'activité s'éteint) ; si $\rho > 1$, il est supercritique (les cascades se propagent indéfiniment).

Le problème central abordé est que les cascades empiriques sont rarement isolées. Elles sont fréquemment pilotées par des champs externes partagés et non stationnaires (par exemple, la productivité de type loi d'Omori pour les répliques, les changements de régime en finance, ou les annonces publiques sur les réseaux sociaux). Lorsque ce pilotage commun est omis du modèle et remplacé par une ligne de base constante, le champ externe n'est pas supprimé mais est plutôt replié (aliased) dans la matrice d'interaction estimée. Ce document examine si cette omission peut rendre un réseau physiquement sous-critique apparent comme critique, conduisant ainsi à de faux diagnostics de transitions de phase.

Méthodologie
Le document emploie une combinaison de théorie de la perturbation spectrale, d'analogies de physique statistique et d'audits empiriques sur des ensembles de données synthétiques et réels.

1. Cadre Théorique :

   • Modèle : Les auteurs considèrent un processus de comptage multivarié dont l'intensité conditionnelle dépend d'un fond $\mu(t)$ et de l'historique d'événements passés. La véritable matrice de branchement est $B_0$.
   • Mécanisme de Repliement (Aliasing) : Lorsque le véritable fond est un champ de rang faible et non stationnaire $\mu(t) = Cf(t)$ (où $C$ est une matrice de charge et $f(t)$ un signal temporel partagé) mais que le modèle ajusté suppose une ligne de base constante, la matrice de reproduction estimée devient $\hat{B} = B_0 + \Delta$.
   • Décalage Spectral : En utilisant la théorie de la perturbation du premier ordre, le document démontre que le décalage du paramètre critique apparent est régi par la projection du terme de repliement $\Delta$ sur le mode dominant (mode de Perron) de la matrice réelle $B_0$. Spécifiquement, pour une petite perturbation, $\rho(\hat{B}) - \rho(B_0) \approx \frac{u_0^\top \Delta v_0}{u_0^\top v_0}$, où $u_0$ et $v_0$ sont les vecteurs de Perron gauche et droit.
   • Contraintes de Rang : Le Théorème 1 établit que le terme de repliement $\Delta$ a un rang borné par le rang de la commande commune. Le Théorème 3 définit un « intervalle de déflation de rang un », fournissant un intervalle certifié $[\rho_{lo}, \rho_{up}]$ pour le véritable rayon de branchement étant donné un $\hat{B}$ observé et une direction de charge fixe.

2. Simulations Contrôlées :

   • Les auteurs simulent des réseaux de branchement sous-critiques ($\rho(B_0) < 1$) pilotés par une enveloppe partagée non stationnaire.
   • Ils comparent des ajustements à ligne de base constante avec des ajustements corrigés par un oracle pour démontrer l'émergence d'un « coin de fausse criticité » : une région où l'augmentation de l'amplitude de la commande commune pousse l'estimation de $\rho(\hat{B})$ au-dessus de 1, bien que la dynamique sous-jacente reste sous-critique.

3. Audits Empiriques :

   • Catalogues de Répliques Sismiques : Six catalogues sismiques spatialement résolus (Maule, Iquique, L'Aquila, Canterbury, Tocopilla, Illapel) sont analysés. Les auteurs partitionnent les événements en zones spatiales et ajustent des réseaux de Hawkes multivariés. Ils comparent les ajustements naïfs à ligne de base constante avec des méthodes de déclustrement espace-temps qui tiennent explicitement compte de la décroissance partagée de type Omori.
   • Flux Twitter du Higgs : Un ensemble de données de mentions, réponses et retweets entourant l'annonce du boson de Higgs en 2012 est analysé. L'annonce agit comme un champ externe global connu. Les auteurs comparent un ajustement à fond constant avec une régression de fond temporel flexible utilisant des nœuds de base temporels lissés.

Résultats Clés

• Théoriques : Le document prouve qu'une commande externe de faible dimension, lorsqu'elle est omise, agit comme une perturbation de rang faible sur la matrice d'interaction. La criticité apparente n'est pas déterminée par la positivité élémentaire du biais, mais par la projection de ce biais sur le mode de Perron du réseau. Même si le biais contient des entrées négatives, une projection positive sur le mode dominant est suffisante pour gonfler le rayon spectral.
• Simulations : Des expériences contrôlées confirment l'existence d'un coin de fausse criticité. À mesure que l'amplitude de la commande commune augmente, une matrice sous-critique fixe est inférée comme étant supercritique. Les simulations de Monte Carlo montrent que les estimations naïves s'accumulent près ou au-dessus du seuil critique ($\rho=1$), tandis que les estimations corrigées par un oracle restent centrées sur la valeur sous-critique réelle.
• Analyse des Répliques : Dans trois des six catalogues (Iquique, L'Aquila, Tocopilla), le rayon du réseau naïf à fond constant dépasse 1, suggérant une supercriticalité. Cependant, après application du déclustrement espace-temps pour supprimer le champ partagé de type Omori, le rayon du réseau corrigé retombe en dessous de 1 dans tous les cas. La supercriticalité naïve s'avère sensible au partitionnement spatial et aux seuils de magnitude, tandis que le certificat sous-critique demeure robuste.
• Analyse Twitter du Higgs : L'ajustement naïf produit un rayon quasi-critique ($\rho_{naïve} = 0,982$). L'ajustement à commande commune flexible réduit ce dernier de manière significative ($\rho_{flexible} = 0,156$). La matrice de différence est essentiellement de rang un et présente une forte projection de Perron (0,826), correspondant au déflation spectrale observée. Cela confirme que le signal quasi-critique a été fabriqué par l'explosion d'activité partagée faisant suite à l'annonce du CERN.

Signification et Revendications
Le document affirme identifier un mécanisme général — le repliement de commande commune — par lequel des réseaux auto-excitateurs sous-critiques peuvent apparaître critiques. La principale signification réside dans le recadrage de l'interprétation des rapports de branchement élevés dans les systèmes complexes :

1. Prudence Diagnostique : Un rayon de reproduction $\rho > 1$ n'est pas une preuve suffisante d'une transition physique vers un branchement supercritique si le système est observé sous un fond non stationnaire partagé.
2. Mécanisme Spectral : Le décalage du paramètre de contrôle est un phénomène spectral piloté par la projection du champ omis sur le mode de Perron dominant du réseau, plutôt que par un simple biais additif.
3. Flux de Travail Méthodologique : Les auteurs proposent un flux de travail diagnostique : estimer l'opérateur apparent, tester la présence d'une structure de commande commune de rang faible dans les résidus, calculer la projection de Perron de la composante supprimée, et rapporter un intervalle de déflation ou une correction structurelle plutôt que d'accepter le $\rho > 1$ naïf comme une transition critique.

Le document conclut que dans les domaines avec des champs externes connus (comme la sismologie), le déclustrement structurel peut récupérer la véritable nature sous-critique du processus de branchement. Dans les domaines avec une connaissance partielle, l'intervalle de déflation de rang un sert à identifier l'ambiguïté de la « fausse criticité », empêchant la mauvaise interprétation d'un forçage externe comme une criticité endogène.