Moments of Random Permanents and Hafnians: A Replica Approach

技術要約：ランダムなパーマネントおよびハフニアンのモーメント：レプリカ法によるアプローチ

問題提起
本論文は、独立かつ同一分布（i.i.d.）に従う成分を持つ$N \times N$（または$2N \times 2N$）ランダム行列における、ランダム行列のパーマネントおよびハフニアンの統計力学、具体的にはそのモーメント $\langle | \text{per } M |^{2k} \rangle$ および $\langle | \text{hf } M |^{2k} \rangle$ の計算を取り扱う。行列式（determinant）のモーメントは行列のスペクトルとの関連からよく理解されているが、パーマネントはスペクトル表現を欠き、計算量的に#P完全である。この困難さゆえに、ランダムなパーマネントの高次モーメントはこれまでほとんど未開拓であったが、これらはフォトニック量子コンピューティング実験（単一光子のボソンサンプリングや、スクイーズド状態のガウス型ボソンサンプリング）の出力分布を定量化する上で極めて重要な役割を果たす。著者らは、実数および複素ガウスアンサンブル、ならびに一般的な対称分布の両方について、指数的な成長率、前係数、および有限サイズ（$1/N$）補正に関する閉形式の表現を導出することを目指している。

手法
著者らは、平均場スピングラスの分配関数のモーメント計算に類似した「レプリカ・アプローチ」を採用している。核心となる手法は以下の3つのステップからなる：

1. 複素平面における輪郭積分表示： パーマネントおよびハフニアンは、複素平面における多重輪郭積分として表現される。行列 $M$ に対し、パーマネントは以下のように表される：
   $$\text{per}(M) = \prod_{i=1}^N \oint \frac{du_i}{2\pi i u_i^2} \prod_{j=1}^N \oint \frac{dv_j}{2\pi i v_j^2} \exp\left( \sum_{i,j} u_i M_{ij} v_j \right)$$
   これにより、パーマネントの累乗の平均化という問題が、補助変数に関する積分へと変換される。

2. ハバード・ストラトノビッチ（HS）デカップリング： ランダム行列の成分 $M_{ij}$ についての平均化を行うために、著者らは複素（または実数）HS場を導入する。これにより、行列の行と列がデカップリングされ、問題はこれらの補助場に関する積分へと変換される。ガウス分布に従う成分の場合、これはHS場行列のハフニアン（実行列の場合）またはパーマネント（複素行列の場合）を含む有効作用をもたらす。

3. 鞍点解析： 大規模な $N$ の極限において、積分は鞍点法を用いて評価される。著者らは、鞍点の安定性を分析し、**レプリカ対称（RS）**解と潜在的なレプリカ対称性の破れ（RSB）を区別する。著者らは、正規化されたモーメントを $Y$ としたとき、レート関数 $\Lambda(k) = \lim_{N \to \infty} \frac{1}{N} \ln \langle Y^k \rangle$ を導出し、鞍点周りのガウス揺らぎを計算することで、前係数および最初の $1/N$ 補正項を得る。

主な貢献と結果

• $k > 2$ におけるレプリカ対称性： $2k$ 次のモーメントについて、支配的な鞍点はすべての $k > 2$ でレプリカ対称である。著者らは、正確な指数的成長率 $\Lambda(k)$、前係数 $C_k$、および最初の $1/N$ 補正項 $a_1/N$ を閉形式で導出した。

  • 実ガウス行列： レート関数は $\Lambda_R(k) = 2k \ln 2 + 2 \ln \Gamma(k + 1/2) - \ln \pi - k \ln(2k - 1)$ である。前係数および補正項は成分分布の累積モーメントに依存する。具体的には、第四累積モーメント $\kappa_4$ が前係数を修正し、高次の累積モーメントは $1/N$ 項のみに影響を与える。
  • 複素ガウス行列： レート関数は $\Lambda_C(k) = 2 \ln \Gamma(k+1) - k \ln k$ である。著者らは、正確なモーメント・サイズ双対性 $\langle |\text{per } M|^{2k} \rangle_N = \langle |\text{per } M|^{2N} \rangle_k$ を確立した。表現論を通じて以前に導出されたこの双対性が、ここでのHSデカップリングを通じて回収されている。

• $k=2$ での不安定性と大偏差： $k=2$ において、レプリカ対称解は不安定になる（ヘッセ行列がゼロまたは負の固有値を持つ）。これは、「分布のバルク（本体）」と「大偏差の裾（テール）」の間の転移を合図するものである。著者らは裾の挙動を計算し、転移領域において普遍的なべき乗則の減衰 $P(y) \sim y^{-3}$（正規化された変数に対して）を示し、続いて伸長指数関数的な裾を示すことを発見した。バルク分布（$k \le 2$）にはRSBの手法が必要であり、これは今後の課題として記されている。

• ハフニアンとガウス型ボソンサンプリング： 著者らは、実ガウス・パーマネントのモーメントと、対称複素ガウス行列のハフニアンのモーメントとの間の双対性を導出した。具体的には、$\langle (\text{per } M)_{real}^{2k} \rangle_N = \langle |\text{hf } M|^{2N} \rangle_{2k}$ である。この結果は、出力確率がハフニアンの二乗に比例するガウス型ボソンサンプリングに直接適用可能である。

• 普遍性： 指数的成長率 $\Lambda(k)$ は、すべてのモーメントを持つ任意の成分分布（例：ラダーマー $\pm 1$ 成分）に対して普遍的であり、成分が実数か複素数かにのみ依存することが示されている。非ガウス的な特徴（高次の累積モーメント）は、劣次的な前係数および有限サイズ補正にのみ影響を与える。

意義と主張
本論文は、無秩序系の手法を用いて、ランダムなパーマネントおよびハフニアンの高次モーメントに対する初の体系的な解析的取り扱いを提供したと主張している。著者らは以下の点を強調している：

1. 解析的アクセス： 組み合わせ論的な計数問題を解可能な統計力学の問題へと見事に変換し、小さな $N$（例：$N=8, 9$）に対する厳密な数値計算と驚くほどよく一致する結果を得た。
2. 実験的関連性： これらの結果は、フォトニック実験における出力確率分布の「ヘビーテイル（重い裾）」を定量化する。著者らは、高次の分散（モーメント $k \ge 3$ によって決定される）は、クロスエントロピー・スコアのような統計的ベンチマークが $N$ に対して指数関数的に増大するサンプル数を必要とすることを意味しており、標準的なガウス信頼区間は分布のヘビーテイルな性質のために不十分であることを指摘している。
3. サンプリングと評価の区別： 単一のパーマネントを計算することは困難であるが、ランダムなパーマネメントの統計的性質（サンプリングの困難さを決定するもの）は、これらのモーメントを通じて特徴付けられることを明確にしている。
4. 限界： 著者らは、自身のレプリカ対称アプローチが $k > 2$ で有効であることを控えめに述べている。分布のバルク（$k \le 2$）の挙動（反集中性（anti-concentration）など、困難さの証明に不可欠な性質を決定するもの）は、この特定の手法では到達できず、将来のRSB計算を必要とする。

本研究は、ランダム行列の理論とスピングラスの物理学を結びつけ、行列式の「探索されにくい対応物」（パーマネント）が、レプリカ対称性とその破れの観点から理解できることを示唆している。