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

技術要約：自己励起ネットワークにおける共通駆動によるエイリアシングに伴う見かけ上の臨界性

問題提起
カスケード・システムの統計物理学において、臨界性はしばしばスカラーの制御パラメータ、すなわち有効な分岐数または繁殖数によって診断される。マルチバリエート・ホークス過程（multivariate Hawkes processes）でモデル化される自己励起ネットワークでは、このパラメータは再生行列（reproduction matrix）のペロン・フロベニウス半径（$\rho$）に対応する。$\rho < 1$ であればシステムは亜臨界（活動が減衰する）であり、$\rho > 1$ であれば超臨界（カスケードが無限に伝播する）である。

本稿が扱う中心的な問題は、経験的なカスケードは決して孤立したものではないということである。それらは頻繁に、共有された非定常な外部場（例：オモリ則に従う余震の生産性、金融におけるレジーム変化、あるいはソーシャルメディアにおける公的な発表など）によって駆動されている。このような共通の駆動がモデルから除外され、定数ベースラインに置き換えられた場合、外部場は除去されるのではなく、推定された相互作用行列へとエイリアシング（折り返し雑音化）される。本論文は、この除外が、物理的には亜臨界であるネットワークを、あたかも臨界であるかのように見せかけ、結果として相転移の誤診を招く可能性があるかどうかを調査するものである。

手法
本稿では、スペクトル摂動論、統計物理学的なアナロジー、および合成データセットと実世界のデータセットを用いた経験的な監査を組み合わせて用いる。

1. 理論的枠組み:

   • モデル: 著者らは、条件付き強度が背景となる $\mu(t)$ と過去のイベント履歴に依存するマルチバリエート計数過程を検討する。真の分岐行列を $B_0$ とする。
   • エイリアシングのメカニズム: 真の背景が低ランクの非定常場 $\mu(t) = Cf(t)$ （ここで $C$ はローディング行列、$f(t)$ は共有された時間信号）であるにもかかわらず、適合させたモデルが定数ベースラインを仮定している場合、推定された再生行列は $\hat{B} = B_0 + \Delta$ となる。
   • スペクトル・シフト: 摂動論を用いて、著者らは、エイリアシング項 $\Delta$ が真の行列 $B_0$ の支配的な固有モード（ペロンモード）への射影によって決定されることを導出する。具体的には、微小な摂動に対して、$\rho(\hat{B}) - \rho(B_0) \approx \frac{u_0^\top \Delta v_0}{u_0^\top v_0}$ となる。ここで $u_0$ と $v_0$ は左および右のペロンベクトルである。
   • ランク制約: 定理1は、エイリアシング項 $\Delta$ のランクが共通駆動のランクによって制限されることを確立する。定理3は「ランク1のデフレーション区間」を定義し、観測された $\hat{B}$ と固定されたローディング方向から、真の分岐半径に対する証明された範囲 $[\rho_{lo}, \rho_{up}]$ を提供する。

2. 制御されたシミュレーション:

   • 著者らは、共有された非定常エンベロープによって駆動される亜臨界分岐ネットワークをシミュレートする。
   • 著者らは、定数ベースラインによる適合と、オラクル（正解）による補正適合を比較することで、「偽の臨界ウェッジ（false-critical wedge）」の出現を実証する。これは、共通の駆動の振幅を増加させることで、基礎となるダイナミクスが亜臨界のままであるにもかかわらず、推定された $\rho(\hat{B})$ が 1 を超えてしまう領域である。

3. 経験的監査:

   • 6つの空間的に解像された地震カタログ（Maule, Iquique, L'Aquila, Canterbury, Tocopilla, Illapel）を分析する。著者らはイベントを空間ゾーンに分割し、マルチバリエート・ホークス・ネットワークを適合させる。彼らは、素朴な定数ベースライン適合と、オモリ型の共有減衰を明示的に考慮した時空間デクラスタリング（脱クラスター化）手法を比較する。
   • ヒッグス粒子に関するTwitterのストリーム・データセット（メンション、リプライ、リツイートのイベント）を分析する。このデータセットは、CERNによる発表という既知のグローバルな外部場を持つ。著者らは、定数背景による適合と、滑らかな時間基底ノットを用いた柔軟な時間背景回帰を比較する。

主要な結果

• 理論的: 著者らは、低次元の外部駆動が、除外された場合に相互作用行列への低ランク摂動として作用することを証明している。見かけ上の臨界性は、バイアスの成分ごとの正負ではなく、このバイアスがネットワークの支配的な固有モード（ペロンモード）といかに整列しているかによって決定される。たとえバイアスに負の成分が含まれていても、支配モードへの正の射影があれば、スペクトル半径を膨張させるのに十分である。
• シミュレーション: 制御された実験は、偽の臨界ウェッジの存在を裏付けている。共通駆動の振幅が増加するにつれ、固定された亜臨界行列が超臨界として推論される。モンテカルロ・シミュレーションは、素朴な推定値が臨界閾値（$\rho=1$）付近またはその上に蓄積する一方で、オラクル補正された推定値は真の亜臨界値を中心にとどまることを示している。
• 余震分析: 6つのカタログのうち3つ（Iquique, L'Aquila, Tocopilla）において、素朴な定数背景ネットワークの半径が1を超え、超臨界性を示唆している。しかし、時空間デクラスタリングを適用して共有のオモリ場を除去した後、補正されたネットワーク半径はすべて1を下回った。素朴な超臨界性は空間分割やマグニチュードの閾値に対して敏感であるが、亜臨界の証明は堅牢であることが示された。
• ヒッグスTwitter分析: 素朴な適合は、臨界に近い半径（$\rho_{naive} = 0.982$）をもたらす。柔軟な共通駆動適合は、これを大幅に減少させる（$\rho_{flex} = 0.156$）。差分行列は本質的にランク1であり、強い正のペロン射影（0.826）を示しており、これは観測されたスペクトル・デフレーションと一致する。これは、CERNの発表に伴う活動の急増によって、臨界に近い信号が作り出されたことを裏付けている。

意義および主張
本論文は、亜臨界の自己励起ネットワークが臨界であるように見える原因となる一般的なメカニズム――共通駆動によるエイリアシング（common-drive aliasing）――を特定したと主張している。本研究の主要な意義は、複雑なシステムにおける高い分岐比の解釈を再定義することにある：

1. 診断上の注意: 共有された非定常な背景の下でシステムが観測されている場合、再生半径 $\rho > 1$ は、物理的な超臨界分岐への遷移を示す十分な証拠にはならない。
2. スペクトル的メカニズム: 制御パラメータのシフトは、単なる加法的なバイアスではなく、ネットワークの支配的な固有モードへの射影によって駆動されるスペクトル現象である。
3. 方法論的ワークフロー: 著者らは、診断ワークフローを提案している。すなわち、見かけ上のオペレーターを推定し、残差に低ランクの共通駆動構造があるかをテストし、除去された成分のペロン射影を計算し、素朴な $\rho > 1$ を臨界転移として受け入れるのではなく、デフレーション区間または構造的補正を報告することである。

本論文は、既知の外部場が存在する領域（地震学など）においては、構造的なデクラスタリングによって分岐過程の真の亜臨界性を回復できると結論付けている。部分的な知識しか得られない領域においては、ランク1のデフレーション区間が「偽の臨界」の曖昧さを特定する役割を果たし、外部からの強制力を内生的な臨界性と誤認することを防ぐ。