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

Technische Samenvatting: Schijnbare Kritikaliteit door Common-Drive Aliasing in Zelf-Exciterende Netwerken

Probleemstelling
In de statistische fysica van cascade-systemen wordt kritikaliteit vaak gediagnosticeerd via een scalaire controleparameter: het effectieve vertakkings- of reproductienummer. In zelf-exciterende netwerken die worden gemodelleerd door multivariate Hawkes-processen, komt deze parameter overeen met de Perron–Frobenius-radius ($\rho$) van de reproductiematrix. Als $\rho < 1$, is het systeem subkritisch (activiteit sterft uit); als $\rho > 1$, is het superkritisch (cascades propageren onbeperkt).

Het centrale probleem dat hier wordt behandeld, is dat empirische cascades zelden geïsoleerd zijn. Ze worden frequent gedreven door gedeelde, niet-stationaire externe velden (bijv. Omori-wet na-schok productiviteit, regimeveranderingen in de financiële sector, of publieke aankondigingen op sociale media). Wanneer een dergelijke gemeenschappelijke drive (common drive) wordt weggelaten uit het model en vervangen door een constante baseline, wordt het externe veld niet verwijderd maar in plaats daarvan gealiaseerd in de geschatte interactiematrix. Dit onderzoek onderzoekt of het weglaten hiervan een fysiek subkritisch netwerk kritiek kan doen lijken, wat leidt tot foutieve diagnoses van faseovergangen.

Methodologie
Het artikel maakt gebruik van een combinatie van spectrale perturbatietheorie, analogieën uit de statistische fysica en empirische audits over zowel synthetische als real-world datasets.

1. Theoretisch Kader:

   • Model: De auteurs beschouwen een multivariate telproces waarbij de conditionele intensiteit afhankelijk is van een achtergrond $\mu(t)$ en een historie van voorgaande gebeurtenissen. De ware vertakkingsmatrix is $B_0$.
   • Aliasing-mechanisme: Wanneer de ware achtergrond een laag-rang, niet-stationair veld is $\mu(t) = Cf(t)$ (waarbij $C$ een laadmatrix is en $f(t)$ een gedeeld temporeel signaal), maar het passende model een constante baseline veronderstelt, wordt de geschatte reproductiematrix $\hat{B} = B_0 + \Delta$.
   • Spectrale Verschuiving: Met behulp van eerste-orde perturbatietheorie leiden de auteurs af dat de verschuiving in de schijnbare kritieke parameter wordt bepaald door de projectie van de aliasing-term $\Delta$ op de dominante eigenmodus (Perron-modus) van de ware matrix $B_0$. Specifiek, voor een kleine perturbatie, geldt $\rho(\hat{B}) - \rho(B_0) \approx \frac{u_0^\top \Delta v_0}{u_0^\top v_0}$, waarbij $u_0$ en $v_0$ de linker- en rechter Perron-vectoren zijn.
   • Rang-beperkingen: Stelling 1 stelt vast dat de aliasing-term $\Delta$ een rang heeft die begrensd wordt door de rang van de gemeenschappelijke drive. Stelling 3 definieert een "rank-one deflation interval", wat een gecertificeerd bereik $[\rho_{lo}, \rho_{up}]$ biedt voor de ware vertakkingsradius, gegeven een geobserveerde $\hat{B}$ en een vaste laadrichting.

2. Gecontroleerde Simulaties:

   • De auteurs simuleren subkritische vertakkingsnetwerken ($\rho(B_0) < 1$) die worden gedreven door een gedeelde niet-stationaire envelop.
   • Zij vergelijken fits met een constante baseline met oracle-gecorrigeerde fits om de opkomst van een "false-critical wedge" aan te tonen: een regio waar het verhogen van de amplitude van de gemeenschappelijke drive de geschatte $\rho(\hat{B})$ boven 1 duwt, ondanks dat de onderliggende dynamica subkritisch blijft.

3. Empirische Audits:

   • Naspiek-catalogi (Aftershock Catalogs): Zes ruimtelijk opgeloste aardbevingscatalogi (Maule, Iquique, L'Aquila, Canterbury, Tocopilla, Illapel) worden geanalyseerd. De auteurs verdelen gebeurtenissen in ruimtelijke zones en passen multivariate Hawkes-netwerken toe. Zij vergelijken naïeve constant-baseline fits met ruimte-tijd declustering-methoden die expliciet rekening houden met de gedeelde Omori-type verval.
   • Higgs Twitter Stream: Een dataset van vermeldingen, reacties en retweets rond de aankondiging van het Higgs-boson in 2012 wordt geanalyseerd. De aankondiging fungeert als een bekende globale externe veld. De auteurs vergelijken een constant-achtergrond fit met een flexibele tijd-achtergrond regressie met gladde tijd-basis knopen (smooth time-basis knots).

Belangrijkste Resultaten

• Theoretisch: Het artikel bewijst dat een laag-dimensionale externe drive, wanneer deze wordt weggelaten, fungeert als een laag-rang perturbatie op de interactiematrix. De schijnbare kritikaliteit wordt niet bepaald door de entrywise positiviteit van de bias, maar door de uitlijning van deze bias met de Perron-modus van het netwerk. Zelfs als de bias negatieve termen bevat, is een positieve projectie op de dominante modus voldoende om de spectrale radius op te blazen.
• Simulaties: Gecontroleerde experimenten bevestigen het bestaan van een "false-critical wedge". Naarmate de amplitude van de gemeenschappelijke drive toeneemt, wordt een vaste subkritische matrix als superkritisch geïnterpreteerd. Monte Carlo-simulaties laten zien dat naïeve schattingen zich ophopen nabij of boven de kritieke drempel ($\rho=1$), terwijl oracle-gecorrigeerde schattingen gecentreerd blijven rond de ware subkritische waarde.
• Aftershock Analyse: In drie van de zes catalogi (Iquique, L'Aquila, Tocopilla) overschrijdt de naïeve constant-achtergrond netwerkradius 1, wat wijst op superkritikaliteit. Echter, na het toepassen van ruimte-tijd declustering om het gedeelde Omori-veld te verwijderen, valt de gecorrigeerde netwerkradius in alle gevallen onder 1. De naïeve superkritikaliteit blijkt gevoelig voor ruimtelijke partitionering en magnitude-drempels, terwijl het subkritische certificaat robuust blijft.
• Higgs Twitter Analyse: De naïeve fit levert een bijna-kritieke radius op ($\rho_{naive} = 0,982$). De flexibele common-drive fit vermindert dit aanzienlijk ($\rho_{flex} = 0,156$). De verschilmatrix is in essentie rank-één en vertoont een sterke positieve Perron-projectie (0,826), wat overeenkomt met de geobserveerde spectrale deflatie. Dit bevestigt dat het bijna-kritieke signaal werd gefabriceerd door de gedeelde burst van activiteit volgend op de CERN-aankondiging.

Betekenis en Claims
Het artikel claimt een algemeen mechanisme te hebben geïdentificeerd — common-drive aliasing — waardoor subkritische zelf-exciterende netwerken kritiek kunnen lijken. De primaire betekenis ligt in het herformuleren van de interpretatie van hoge vertakkingsratio's in complexe systemen:

1. Diagnostische Voorzichtigheid: Een reproductieradius $\rho > 1$ is niet voldoende bewijs voor een fysieke transitie naar superkritische vertakking als het systeem wordt geobserveerd onder een gedeelde niet-stationaire achtergrond.
2. Spectraal Mechanisme: De verschuiving in de controleparameter is een spectraal fenomeen dat wordt gedreven door de projectie van het weggelaten veld op de dominante eigenmodus van het netwerk, in plaats van een eenvoudige additieve bias.
3. Methodologische Workflow: De auteurs stellen een diagnostische workflow voor: schat de schijnbare operator, test op laag-rang common-drive structuren in de residuen, bereken de Perron-projectie van het verwijderde component, en rapporteer een deflatie-interval of een structurele correctie in plaats van de naïeve $\rho > 1$ te accepteren als een kritieke transitie.

Het artikel concludeert dat in domeinen met bekende externe velden (zoals seismologie), structurele declustering de ware subkritische aard van het vertakkingsproces kan herstellen. In domeinen met gedeeltelijke kennis dient het rank-one deflation interval om de "false-critical" ambiguïteit te identificeren, wat misinterpretatie van externe forcering als endogene kritikaliteit voorkomt.