Spectral origin of the topological gap exponent d + {\eta}: mechanism, kernel, decomposition, and scope

技術的要約：トポロジカル・ギャップ指数 $d + \eta$ のスペクトル起源

問題提起
臨界スピンモデル（例：イジング、ポッツ）に由来する点集合に適用されるパーシステント・ホモロジー（PH）は、「トポロジカル・ギャップ」$\Delta$ を明らかにする。これは、臨界構成における $H_1$ の全パーシステンスと、密度を一致させたシャッフルされた帰無モデルとの差として定義される。経験的な研究により、このギャップはシステムサイズ $L$ に対して $\Delta \sim L^{d+\eta}$ （ここで $d$ は空間次元、$\eta$ は異常次元）とスケーリングすることが確立されている。しかし、この特定の指数の物理的および解析的な起源は、第一原理からの説明がなされていない。さらに、このメカニズムの諸次元および普遍性クラスにわたる妥当性と、スケーリングを駆動するスペクトル・カーネルの性質には、厳密な導出と検証が必要である。

手法
本論文は、解析的な導出と広範な数値シミュレーションを組み合わせて用いている：

1. 解析的枠組み： 著者は、$\Delta$ がスペクトル積分 $I(\alpha) = \sum_{k \neq 0} S_{conn}(k) |k|^\alpha$ （ここで $S_{conn}(k)$ は連結構造因子）で表されると仮定している。臨界状態における $S_{conn}(k) \sim |k|^{-(2-\eta)}$ の赤外（IR）および紫外（UV）挙動を分析することで、$L^{d+\eta}$ 指数を回収するために必要なスケーリング条件を導出した。
2. 分解： スペクトル積分は、体積因子 $I_0$（全スペクトル重み）と形状因子 $I_{shape}$（スペクトル分布）に分解される。$I_0$ はパーセバルの等式を通じて磁化に関連付けられる。
3. 数値シミュレーション： スウェンズンド・ワング・クラスター更新法を用いて、2Dイジングモデル、2Dポッツモデル（$q=3, 4$）、および3Dイジングモデルに対してシミュレーションを行った。点集合は多数派スピンのサイトから生成された。
   • ツール： アルファ複体（alpha complex）の構築とPH計算にはGUDHIを使用；スペクトル積分にはFFTを使用。
   • テスト： 有限サイズスケーリング（FSS）解析（$I(\alpha)$）、比率の安定性テスト（$\Delta / I(\alpha)$）、構成ごとの相関解析（磁化を部分化）、およびポッツ $q=4$ モデルにおける最適なスペクトル重みを決定するための $\alpha$ スイープ。
4. ヘウリスティックな議論： 不均一なポアソン過程に対する生のPHカーネルの「平坦性」に関する推測が、なぜ実効的なカーネルが臨界状態で初めて $k$ 依存性を獲得するのかを説明するために提案されている。

主要な貢献と結果

• 指数の解析的導出： 本論文は、もしスペクトル積分がIR支配的であれば、$I(\alpha) \sim L^{2-\alpha-\eta}$ となることを証明した。$\alpha = -2\eta$ と設定すると、$d=2$ において特に $I(-2\eta) \sim L^{d+\eta}$ が得られる。これは、観測された指数を、臨界構造因子に直接結びつけることで、第一原理的な説明を提供している。
• 臨界次元の基準： IR支配メカニズムの必要条件として $d < 2 + \eta$ が導出された：
  • 2D では、この条件は任意の $\eta > 0$ に対して成立し、このメカニズムが本質的に二次元的であることを裏付けている。
  • 3D では、$\eta \ll 1$ であるため、$d < 2+\eta$ の条件が成立しない。したがって積分はUV支配となり（$I \sim L^d$）、先行研究で見られたように、なぜ $L^{d+\eta}$ スケーリングを回収するために密度正規化が必要なのかを説明している。
• 分散の分解： $\Delta$ の構成ごとの分散は、ほとんどすべてが磁化と強く相関する体積因子 $I_0 \propto (1-m^2)$ によって駆動されていることを示している。スペクトル形状因子 $I_{shape}$ は $L^\eta$ スケーリングに寄与するが、構成ごとの分散は低い。その結果、構成レベルでの $\Delta$ と $I(\alpha)$ の間の高い $R^2$ 相関は、スペクトルの形状メカニズムの直接的な証拠ではなく、主に共有された磁化依存性のアーティファクトである。
• クロスモデル検証：
  • 2Dイジングおよびポッツ $q=3$： 比 $\Delta / I(-2\eta)$ はシステムサイズ（$L=32$ から $512$）にわたって安定しており、比例関係と指数を確認している。
  • ポッツ $q=4$： $\alpha$ スイープ（$L=32$–$256$）により、最適なスペクトル重み$\alpha_{opt}$が$[-0.75, -0.5]$の範囲にあることが判明した。これはイジングの値（$-0.5$）とは一致するが、モデル固有の予測（$-2\eta_{q=4} = -1.0$）とは一致しない。著者は、これを$L=512$ での周辺ケースの自己相関問題や対数補正の可能性によるものとして、暫定的な結果であると指摘している。
• ハイパー・スケーリングと共分散： 積 $\langle m^2 \cdot I(-2\eta) \rangle$ がハイパー・スケーリングの診断として分析された。研究の結果、この積は $I_0$ を介したスペクトル積分と磁化の間の強い負の共分散（$r \approx -0.98$）によって支配されている。したがって、この積は独立したハイパー・スケーリングのテストではなく、$\langle m^2 \rangle$ と $\langle I \rangle$ の個別の測定が必要である。
• カーネルの平坦性： ヘウリスティックな議論は、不均一な強度の関数に対する生のPHカーネルは波数に依存せず平坦であることを示唆している。$|k|^{-2\eta}$ カーネルの出現は、臨界的な相関がアルファ複体の離散的なトポロジーと結合することから生じる、非摂動的な「ドレッシング（被覆）」効果として解釈される。

意義と主張
本論文は、トポロジカル・ギャップ指数 $d+\eta$ を、経験的な観察を超えて、臨界構造因子から初めて解析的に導出したと主張している。また、このメカニズムが本質的に二次元的であり、赤外（IR）支配に依存していること、および体積（磁化駆動）因子と形状（異常次元）因子の異なる役割を明確にしている。

著者は、スペクトル重みの普遍性については控えめな姿勢をとっている。メカニズムは2Dイジングおよびポッツ $q=3$ に対して確認されているが、ポッツ $q=4$ に関する結果は、周辺ケースの不安定性を解決するための大規模スケール（$L \ge 1024$）でのマルチシード熱平衡化を待つ必要がある「暫定的な観察」として提示されている。もし $\alpha_{opt} \approx -0.5$ という結果が確認されれば、それはスペクトル重みがモデル固有の異常次元ではなく、2Dにおけるアルファ複体の幾何学的な定数であることを示唆し、公式を $d + \eta_{eff}$ と再解釈することになる。

最後に、本論文はトポロジカルな観測量と磁化の関係を再構成し、構成ごとの相関がスペクトルの形状ではなく、多くの場合、共有された $I_0$ 振幅によって駆動されていることに警告を発している。また、ハイパー・スケーリングのテストは、この共分散を考慮する必要があるとしている。本研究は、非整数次のスペクトルモーメントを非ガウス的で密度変調されたスピン系に適用するという点で、ガウス・キネマティック公式（GKF）の系譜の中に、独自の立ち位置を確立している。