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

Sintesi Tecnica: Criticità Apparente da Aliasing di Drive Comune in Reti Auto-eccitanti

Problematica
Nella fisica statistica dei sistemi a cascata, la criticità viene spesso diagnosticata tramite un parametro di controllo scalare: il numero effettivo di ramificazione o di riproduzione. In reti auto-eccitanti modellate da processi di Hawkes multivariati, questo parametro corrisponde al raggio di Perron–Frobenius ($\rho$) della matrice di riproduzione. Se $\rho < 1$, il sistema è subcritico (l'attività si esaurisce); se $\rho > 1$, il sistema è supercritico (le cascate si propagano indefinitamente).

Il problema centrale affrontato è che le cascate empiriche sono raramente isolate. Esse sono frequentemente guidate da campi esterni comuni e non stazionari (ad esempio, la produttività delle scosse di assestamento secondo la legge di Omori, cambiamenti di regime in finanza o annunci pubblici sui social media). Quando tale drive comune viene omesso dal modello e sostituito da una linea di base costante, il campo esterno non viene rimosso, ma viene invece "aliassato" nella matrice di interazione stimata. Questo articolo investiga se tale omissione possa causare in una rete fisicamente subcritica l'apparenza di criticità, portando a false diagnosi di transizioni di fase.

Metodologia
L'articolo impiega una combinazione di teoria della perturbazione spettrale, analogie di fisica statistica e audit empirici su dataset sintetici e del mondo reale.

1. Quadro Teorico:

   • Modello: Gli autori considerano un processo di conteggio multivariato la cui intensità condizionata dipende da un background $\mu(t)$ e da una storia di eventi passati. La vera matrice di ramificazione è $B_0$.
   • Meccanismo di Aliasing: Quando il vero background è un campo a basso rango e non stazionario $\mu(t) = Cf(t)$ (dove $C$ è una matrice di carico e $f(t)$ è un segnale temporale condiviso), ma il modello stimato assume una linea di base costante, la matrice di riproduzione stimata diventa $\hat{B} = B_0 + \Delta$.
   • Spostamento Spettrale: Utilizzando la teoria della perturbazione del primo ordine, l'articolo deriva che lo spostamento nel parametro critico apparente è governato dalla proiezione del termine di aliasing $\Delta$ sul modo dominante (modo Perron) della matrice vera $B_0$. Nello specifico, per una piccola perturbazione, $\rho(\hat{B}) - \rho(B_0) \approx \frac{u_0^\top \Delta v_0}{u_0^\top v_0}$, dove $u_0$ e $v_0$ sono i vettori di Perron sinistro e destro.
   • Vincoli di Rango: Il Teorema 1 stabilisce che il termine di aliasing $\Delta$ ha un rango limitato dal rango del drive comune. Il Teorema 3 definisce un "intervallo di deflazione di rango uno", fornendo un intervallo certificato $[\rho_{lo}, \rho_{up}]$ per il vero raggio di ramificazione dato un $\hat{B}$ osservato e una direzione di carico fissata.

2. Simulazioni Controllate:

   • Gli autori simulano reti di ramificazione subcritiche ($\rho(B_0) < 1$) guidate da un inviluppo comune non stazionario.
   • Confrontano fit a baseline costante con fit corretti tramite "oracle" per dimostrare l'emergere di un "cuneo di falsa criticità": una regione in cui l'aumento dell'ampiezza del drive comune spinge il $\rho(\hat{B})$ stimato sopra 1, nonostante la dinamica sottostante rimanga subcritica.

3. Audit Empirici:

   • Cataloghi di Scosse di Assestamento: Vengono analizzati sei cataloghi di terremoti spazialmente risolti (Maule, Iquique, L'Aquila, Canterbury, Tocopilla, Illapel). Gli autori partizionano gli eventi in zone spaziali e applicano modelli di reti di Hawkes multivariati. Confrontano i fit ingenui a baseline costante con metodi di declustering spazio-temporale che tengono esplicitamente conto del decadimento condiviso di tipo Omori.
   • Flusso Twitter dell'Higgs: Viene analizzato un dataset di menzioni, risposte e retweet relativi all'annuncio della scoperta del bosone di Higgs nel 2012. L'annuncio funge da noto campo esterno globale. Gli autori confrontano un fit a background costante con un modello di regressione a background temporale flessibile con nodi di base temporale liscio.

Risultati Chiave

• Teorici: L'articolo dimostra che un drive esterno a basso rango, se omesso, agisce come una perturbazione a basso rango alla matrice di interazione. La criticità apparente è determinata non dalla positività elemento per elemento del bias, ma dall'allineamento di tale bias con il modo Perron della rete. Anche se il bias contiene voci negative, una proiezione positiva sul modo dominante è sufficiente per gonfiare il raggio spettrale.
• Simulazioni: Gli esperimenti controllati confermano l'esistenza di un cuneo di falsa criticità. All'aumentare dell'ampiezza del drive comune, una matrice fissa subcritica viene inferita come supercritica. Le simulazioni Monte Carlo mostrano che le stime naive si accumulano vicino o sopra la soglia critica ($\rho=1$), mentre le stime corrette tramite oracle rimangono centrate sul vero valore subcritico.
• Analisi delle Scosse di Assestamento: In tre dei sei cataloghi (Iquique, L'Aquila, Tocopilla), il raggio della rete con background costante ingenuo supera 1, suggerendo la supercriticità. Tuttavia, dopo l'applicazione del declustering spazio-temporale per rimuovere il campo Omori condiviso, il raggio della rete corretta scende sotto 1 in tutti i casi. La supercriticità ingenua si è dimostrata sensibile alla partizione spaziale e alle soglie di magnitudo, mentre il certificato di subcriticità rimane robusto.
• Analisi Twitter dell'Higgs: Il fit ingenuo produce un raggio quasi critico ($\rho_{naive} = 0.982$). Il fit con drive comune flessibile riduce significativamente questo valore ($\rho_{flex} = 0.156$). La matrice di differenza è essenzialmente di rango uno ed esibisce una forte proiezione di Perron positiva (0.826), corrispondendo alla deflazione spettrale osservata. Ciò conferma che il segnale di quasi-criticità era generato dal burst di attività condivisa seguito dall'annuncio del CERN.

Significatività e Rivendicazioni
L'articolo rivendica l'identificazione di un meccanismo generale — l'aliasing di drive comune — attraverso il quale reti auto-eccitanti subcritiche possono apparire critiche. La significatività primaria risiede nel riformulare l'interpretazione degli alti rapporti di ramificazione nei sistemi complessi:

1. Cautela Diagnostica: Un raggio di riproduzione $\rho > 1$ non è sufficiente prova di una transizione fisica verso una ramificazione supercritica se il sistema è osservato sotto un background comune non stazionario.
2. Meccanismo Spettrale: Lo spostamento nel parametro di controllo è un fenomeno spettrale guidato dalla proiezione del campo omesso sul modo dominante della rete, piuttosto che da un semplice bias additivo.
3. Workflow Metodologico: Gli autori propongono un workflow diagnostico: stimare l'operatore apparente, testare la presenza di una struttura di drive comune a basso rango nei residui, calcolare la proiezione di Perron della componente rimossa e riportare un intervallo di deflazione o una correzione strutturale, anziché accettare il $\rho > 1$ ingenuo come una transizione critica.

L'articolo conclude che in domini con campi esterni noti (come la sismologia), il declustering strutturale può recuperare la vera natura subcritica del processo di ramificazione. In domini con conoscenza parziale, l'intervallo di deflazione di rango uno serve a identificare l'ambiguità della "falsa criticità", prevenendo l'interpretazione errata del forcing esterno come criticità endogena.