We study the moments of $\overline{|\det(H-E)|^q}$ and the associated large deviations of $\log |\det(H-E)|$ where $H$ are random matrix operators involving Laplace operators and random potentials. This includes as a special case Hessians of random elastic manifolds at a generic energy configuration. In one dimension $d=1$ these are $N \times N$ matrix valued random Schrödinger operators and $\log | \det(H-E) | $ is the sum of the $N$ associated Lyapunov exponents. Using a mapping to a stochastic matrix Ricatti equation we make a connection between the spectral properties of these operators and the total $N$ particle current of a Dyson Brownian motion (DBM) in a cubic potential. The latter model was studied by Allez and Dumaz [1] who showed that for $N=+\infty$ it exhibits a sharp transition between a phase with non-zero current and a confined (zero current) phase. We compute the barrier-crossing probability of the DBM at large but finite $N$, which gives an estimate of the exponential tail of the average density of states of a matrix Schrodinger operator below the edge of its spectrum. The barrier behaves as $\sim N (-E)^{3/2}$ at large negative energy and vanishes as $\sim N(E^*-E)^{5/4}$ near the edge. For $q=1$ the present work provides an independent derivation of the total complexity of stationary points for an elastic string embedded in $N$ dimension in presence of disorder.
論文ID : 2511.00954タイトル : Large deviations of spectral determinants of matrix-valued random Schrödinger operators and Dyson Brownian motion in cubic potentials著者 : Yan Fyodorov (King's College London), Pierre Le Doussal (ENS Paris), Alexander Ossipov (King's College London)分類 : math-ph, cond-mat.dis-nn, cond-mat.stat-mech, math.MP, math.PR提出日時 : 2025年11月2日論文リンク : https://arxiv.org/abs/2511.00954 本論文は、行列値ランダムシュレーディンガー演算子の谱行列式のモーメント ∣ det ( H − E ) ∣ q ‾ \overline{|\det(H-E)|^q} ∣ det ( H − E ) ∣ q およびその大偏差特性を研究する。ここで H H H はラプラシアンと乱数ポテンシャルを含む。一次元の場合、これらは N × N N \times N N × N 行列値ランダムシュレーディンガー演算子であり、log ∣ det ( H − E ) ∣ \log|\det(H-E)| log ∣ det ( H − E ) ∣ は N N N 個のリャプノフ指数の和に等しい。ランダム行列リッカチ方程式への写像を通じて、著者はこれらの演算子の谱特性とキューブポテンシャル中のダイソンブラウン運動(DBM)の N N N 粒子全流との関連性を確立した。大きいが有限の N N N に対するDBMのポテンシャル障壁貫通確率を計算し、行列シュレーディンガー演算子の谱端以下の平均状態密度の指数テールの推定値を与えた。障壁は大きな負のエネルギーで ∼ N ( − E ) 3 / 2 \sim N(-E)^{3/2} ∼ N ( − E ) 3/2 として振る舞い、端付近で ∼ N ( E ∗ − E ) 5 / 4 \sim N(E^*-E)^{5/4} ∼ N ( E ∗ − E ) 5/4 に消滅する。
本論文は行列値ランダムシュレーディンガー演算子の谱行列式統計特性を研究し、特に以下に焦点を当てる:
谱行列式 ∣ det ( H − E ) ∣ q |\det(H-E)|^q ∣ det ( H − E ) ∣ q のモーメントおよびその成長率 Σ q \Sigma_q Σ q 対数谱行列式 log ∣ det ( H − E ) ∣ \log|\det(H-E)| log ∣ det ( H − E ) ∣ の大偏差特性 一次元の場合とリャプノフ指数およびダイソンブラウン運動との関連性 物理的応用 :このような演算子は複数の物理的背景で自然に現れる:アンダーソン局在:多チャネル準一次元系(線または条)のモデル 無秩序弾性系のヘッシアン行列:ガラス状態の特性を記述 DMPK法:無秩序導体理論の数学的基礎 数学的意義 :ランダム行列理論における谱行列式モーメント研究の拡張 スカラー場合(N = 1 N=1 N = 1 )の一般化、一般化されたリャプノフ指数を含む Kac-Rice公式との関連性、エネルギー汎関数の臨界点総数の計算が可能 複雑性理論 :q = 1 q=1 q = 1 の場合、量 Σ 1 \Sigma_1 Σ 1 は臨界点の焼きなまし複雑度(annealed complexity)として解釈され、無秩序系における豊富な力学平衡状態を記述するスカラー場合(N = 1 N=1 N = 1 )は研究されているが、行列値場合(N > 1 N>1 N > 1 )の研究は少ない 先行研究は主に q = 1 q=1 q = 1 の場合に集中している 一般的な q q q 値と大偏差特性の系統的研究が欠けている ダイソンブラウン運動との深い関連性は十分に探索されていない 著者は以下を目指している:
Σ q \Sigma_q Σ q を計算するための2つの相補的方法の開発:行列リッカチ法と鞍点法行列シュレーディンガー演算子とキューブポテンシャル中のDBMの正確な対応の確立 谱端のリフシッツテールとDBMポテンシャル障壁貫通の関係の理解 無秩序弾性系に対する独立した複雑度計算の提供 二重方法論 :モーメント成長率 Σ q \Sigma_q Σ q を計算するための2つの独立で相補的な方法を開発:行列リッカチ方程式法(d = 1 d=1 d = 1 、任意 N N N ) 鞍点法(任意 d d d 、大 N N N 極限) DBM関連性 :行列シュレーディンガー演算子とキューブポテンシャル中のダイソンブラウン運動の深い関連性を確立:谱端 E ∗ E^* E ∗ がDBMの相転移点に対応することを証明 状態密度とDBM粒子流の対応関係 自己平均特性の独立検証 大偏差理論 :ルジャンドル変換を通じて率関数 Φ ( e ) \Phi(e) Φ ( e ) を取得。ここで e = 1 L d N log ∣ det ( H − E ) ∣ e = \frac{1}{LdN}\log|\det(H-E)| e = L d N 1 log ∣ det ( H − E ) ∣ :
Φ ( e ) = 1 2 ( f − 1 ( δ e + f ( μ / J ) ) − μ / J ) 2 \Phi(e) = \frac{1}{2}\left(f^{-1}(\delta e + f(\mu/J)) - \mu/J\right)^2 Φ ( e ) = 2 1 ( f − 1 ( δe + f ( μ / J )) − μ / J ) 2 ポテンシャル障壁貫通分析 :有限 N N N に対するDBMのポテンシャル障壁貫通確率を計算:端付近:U ≃ 4 5 2 × 3 1 / 4 β 1 / 6 ( E ∗ − E ) 5 / 4 U \simeq \frac{4}{5}\sqrt{2 \times 3^{1/4}}\beta^{1/6}(E^*-E)^{5/4} U ≃ 5 4 2 × 3 1/4 β 1/6 ( E ∗ − E ) 5/4 大きな負のエネルギー:U ≃ 4 3 ∣ E ∣ 3 / 2 U \simeq \frac{4}{3}|E|^{3/2} U ≃ 3 4 ∣ E ∣ 3/2 リフシッツテール推定:ρ K ( α ) ∼ exp ( − N U / J ~ 2 ) \rho_K(\alpha) \sim \exp(-NU/\tilde{J}^2) ρ K ( α ) ∼ exp ( − N U / J ~ 2 ) 相図の特性化 :2つの相を識別:単純相 (μ > μ b \mu > \mu_b μ > μ b ):ゼロ複雑度領域(q = 1 q=1 q = 1 時)複雑相 (μ < μ b \mu < \mu_b μ < μ b ):非ゼロ複雑度、谱内部に対応相境界:μ b = μ c − ( q − 1 ) J 2 ∫ k 1 μ c − t Δ ( k ) \mu_b = \mu_c - (q-1)J^2\int_k \frac{1}{\mu_c - t\Delta(k)} μ b = μ c − ( q − 1 ) J 2 ∫ k μ c − t Δ ( k ) 1 入力 :
行列値ランダム演算子 H = K + X + μ I H = K + X + \mu I H = K + X + μ I 、ここで:
K K K :ラプラシアンとGOEランダム行列を含む構造化演算子X X X :対角ガウスノイズパラメータ:N N N (行列次元)、L L L (システムサイズ)、J J J (無秩序強度)、c c c (ノイズ相関パラメータ) 出力 :
モーメント成長率:Σ q = lim L → ∞ 1 N L d log Y q ‾ \Sigma_q = \lim_{L\to\infty} \frac{1}{NL^d}\log\overline{Y_q} Σ q = lim L → ∞ N L d 1 log Y q 、ここで Y q = ∣ det ( H − E ) ∣ q Y_q = |\det(H-E)|^q Y q = ∣ det ( H − E ) ∣ q 大偏差率関数:Φ ( e ) \Phi(e) Φ ( e ) は P ( e ) ∼ e − N L d Φ ( e ) P(e) \sim e^{-NL^d\Phi(e)} P ( e ) ∼ e − N L d Φ ( e ) を満たす リャプノフ指数和(d = 1 d=1 d = 1 ):N e = ∑ j = 1 N γ j Ne = \sum_{j=1}^N \gamma_j N e = ∑ j = 1 N γ j 一次元連続モデルの関数行列式を行列リッカチ方程式解のトレース積分として表現:
det ( H − E ) = y ( L ) , ∣ y ( L ) ∣ = e ∫ 0 L d τ Tr Z ( τ ) \det(H-E) = y(L), \quad |y(L)| = e^{\int_0^L d\tau \text{Tr}Z(\tau)} det ( H − E ) = y ( L ) , ∣ y ( L ) ∣ = e ∫ 0 L d τ Tr Z ( τ )
ここで Z ( τ ) Z(\tau) Z ( τ ) はランダム行列リッカチ方程式を満たす:
∂ τ Z = − E − Z 2 + W ~ ( τ ) , Z ( 0 ) = + ∞ \partial_\tau Z = -E - Z^2 + \tilde{W}(\tau), \quad Z(0) = +\infty ∂ τ Z = − E − Z 2 + W ~ ( τ ) , Z ( 0 ) = + ∞
Z Z Z の固有値 λ i ( τ ) \lambda_i(\tau) λ i ( τ ) を導入すると、ランダム進化方程式を満たす:
d λ i ( τ ) = − ( E + λ i ( τ ) 2 ) d τ + J ~ 2 d τ N ∑ j ≠ i 1 λ i ( τ ) − λ j ( τ ) + 2 J ~ 2 N d B i ( τ ) + J ~ ξ ~ ( τ ) d τ d\lambda_i(\tau) = -(E + \lambda_i(\tau)^2)d\tau + \frac{\tilde{J}^2 d\tau}{N}\sum_{j\neq i}\frac{1}{\lambda_i(\tau)-\lambda_j(\tau)} + \sqrt{\frac{2\tilde{J}^2}{N}}dB_i(\tau) + \tilde{J}\tilde{\xi}(\tau)d\tau d λ i ( τ ) = − ( E + λ i ( τ ) 2 ) d τ + N J ~ 2 d τ ∑ j = i λ i ( τ ) − λ j ( τ ) 1 + N 2 J ~ 2 d B i ( τ ) + J ~ ξ ~ ( τ ) d τ
これはまさにキューブポテンシャル中のダイソンブラウン運動 方程式である。ここで:
キューブポテンシャル:V ( λ ) = E λ + λ 3 / 3 V(\lambda) = E\lambda + \lambda^3/3 V ( λ ) = E λ + λ 3 /3 粒子間対数反発相互作用 ブラウンノイズと共通ノイズ ξ ~ ( τ ) \tilde{\xi}(\tau) ξ ~ ( τ ) 経験密度のStieltjes変換を定義:
G ( z , τ ) = 1 N ∑ i 1 λ i ( τ ) − z G(z,\tau) = \frac{1}{N}\sum_i \frac{1}{\lambda_i(\tau)-z} G ( z , τ ) = N 1 ∑ i λ i ( τ ) − z 1
大 N N N 極限では、G G G は進化方程式を満たす(ξ ~ ( τ ) = ξ ˉ + η ( τ ) / N \tilde{\xi}(\tau) = \bar{\xi} + \eta(\tau)/\sqrt{N} ξ ~ ( τ ) = ξ ˉ + η ( τ ) / N に分解):
∂ τ G ( z , τ ) = ∂ z [ z + ( E + z 2 ) G ( z , τ ) ] + 1 2 J ~ 2 ∂ z G ( z , τ ) 2 \partial_\tau G(z,\tau) = \partial_z[z + (E+z^2)G(z,\tau)] + \frac{1}{2}\tilde{J}^2\partial_z G(z,\tau)^2 ∂ τ G ( z , τ ) = ∂ z [ z + ( E + z 2 ) G ( z , τ )] + 2 1 J ~ 2 ∂ z G ( z , τ ) 2
定常状態解は代数方程式を満たす:
z + ( E + z 2 ) G ( z ) + 1 2 J ~ 2 G ( z ) 2 = J z + (E+z^2)G(z) + \frac{1}{2}\tilde{J}^2 G(z)^2 = \mathcal{J} z + ( E + z 2 ) G ( z ) + 2 1 J ~ 2 G ( z ) 2 = J
ここで積分定数 J ( E ) \mathcal{J}(E) J ( E ) の虚部は粒子流を与える:
j = 1 π Im J ( E ) j = \frac{1}{\pi}\text{Im}\mathcal{J}(E) j = π 1 Im J ( E )
禁閉相 (E < E ∗ = − 3 4 ( 2 J ~ 2 ) 2 / 3 E < E^* = -\frac{3}{4}(2\tilde{J}^2)^{2/3} E < E ∗ = − 4 3 ( 2 J ~ 2 ) 2/3 ):J \mathcal{J} J は実数、ゼロ流、密度は有限サポート流動相 (E > E ∗ E > E^* E > E ∗ ):J \mathcal{J} J は虚部を持つ、非ゼロ流、密度サポートは全実軸重要な識別:E ∗ = − μ c E^* = -\mu_c E ∗ = − μ c はまさに行列シュレーディンガー演算子の谱端である
定常状態平均トレース:
⟨ Tr Z ⟩ s t = − Re J ( E − J ~ ξ ˉ ) \langle \text{Tr}Z \rangle_{st} = -\text{Re}\mathcal{J}(E-\tilde{J}\bar{\xi}) ⟨ Tr Z ⟩ s t = − Re J ( E − J ~ ξ ˉ )
Gelfand-Yaglom関係を利用:
lim N , L → ∞ 1 N L log ∣ y ( L ) ∣ = − Re J ( E − J ~ ξ ˉ ) \lim_{N,L\to\infty} \frac{1}{NL}\log|y(L)| = -\text{Re}\mathcal{J}(E-\tilde{J}\bar{\xi}) lim N , L → ∞ N L 1 log ∣ y ( L ) ∣ = − Re J ( E − J ~ ξ ˉ )
以下を証明した:
d d E Re J ( E − J ~ ξ ˉ ) = − PV ∫ d α ρ K ( α ) α − E + J ~ ξ ˉ \frac{d}{dE}\text{Re}\mathcal{J}(E-\tilde{J}\bar{\xi}) = -\text{PV}\int d\alpha \frac{\rho_K(\alpha)}{\alpha - E + \tilde{J}\bar{\xi}} d E d Re J ( E − J ~ ξ ˉ ) = − PV ∫ d α α − E + J ~ ξ ˉ ρ K ( α )
ここで ρ K ( α ) \rho_K(\alpha) ρ K ( α ) は演算子 K K K の平均谱密度である。
重要な観察:大 N N N 極限では、対数行列式は自己平均を示す:
⟨ e q log ∣ det [ K + X + ( μ − E ) I ] ∣ ⟩ GOE’s ≈ e q ⟨ Tr log ∣ K + X + ( μ − E ) I ∣ ⟩ GOE’s \langle e^{q\log|\det[K+X+(\mu-E)I]|}\rangle_{\text{GOE's}} \approx e^{q\langle\text{Tr}\log|K+X+(\mu-E)I|\rangle_{\text{GOE's}}} ⟨ e q l o g ∣ d e t [ K + X + ( μ − E ) I ] ∣ ⟩ GOE’s ≈ e q ⟨ Tr l o g ∣ K + X + ( μ − E ) I ∣ ⟩ GOE’s
これは問題を以下に簡約する:
Y q ∣ N ≫ 1 ∼ ∏ x ∫ R d ξ ( x ) 2 π / N e − N S [ ξ ] Y_q|_{N\gg 1} \sim \prod_x \int_{\mathbb{R}} \frac{d\xi(x)}{\sqrt{2\pi/N}} e^{-NS[\xi]} Y q ∣ N ≫ 1 ∼ ∏ x ∫ R 2 π / N d ξ ( x ) e − NS [ ξ ]
ここで作用量:
S [ ξ ] = ∑ x 1 2 c ξ ( x ) 2 − q N ⟨ Tr ( log ∣ K + X + ( μ − E ) I ∣ − log ∣ − Δ + μ ∣ ) ⟩ GOE’s S[\xi] = \sum_x \frac{1}{2c}\xi(x)^2 - \frac{q}{N}\langle\text{Tr}(\log|K+X+(\mu-E)I| - \log|-\Delta+\mu|)\rangle_{\text{GOE's}} S [ ξ ] = ∑ x 2 c 1 ξ ( x ) 2 − N q ⟨ Tr ( log ∣ K + X + ( μ − E ) I ∣ − log ∣ − Δ + μ ∣ ) ⟩ GOE’s
作用量の最小化は鞍点条件をもたらす:
ξ q ∗ = q f ′ ( ξ q ∗ + μ / J ) , f ( ξ ) : = ∫ d α ln ∣ α + J ξ ∣ ρ K ( α ) \xi^*_q = qf'(\xi^*_q + \mu/J), \quad f(\xi) := \int d\alpha \ln|\alpha+J\xi|\rho_K(\alpha) ξ q ∗ = q f ′ ( ξ q ∗ + μ / J ) , f ( ξ ) := ∫ d α ln ∣ α + J ξ ∣ ρ K ( α )
重要な性質:鞍点 ξ q ∗ \xi^*_q ξ q ∗ は位置 x x x に独立である (凸性論証により証明)
平均谱密度 ρ K ( α ) \rho_K(\alpha) ρ K ( α ) はresolvent i r λ ir_\lambda i r λ により決定される:
ρ K ( λ ) = 1 π Im ( i r λ ) ∣ Im λ = 0 − \rho_K(\lambda) = \frac{1}{\pi}\text{Im}(ir_\lambda)|_{\text{Im}\lambda=0^-} ρ K ( λ ) = π 1 Im ( i r λ ) ∣ Im λ = 0 −
ここで i r λ ir_\lambda i r λ は以下を満たす:
i r λ = ∫ k 1 λ + t Δ ( k ) − i r λ J 2 ir_\lambda = \int_k \frac{1}{\lambda + t\Delta(k) - ir_\lambda J^2} i r λ = ∫ k λ + t Δ ( k ) − i r λ J 2 1
実部と虚部 i r λ = x λ + i y λ ir_\lambda = x_\lambda + iy_\lambda i r λ = x λ + i y λ を導入し、鞍点 λ = − J ξ q ∗ − μ \lambda = -J\xi^*_q - \mu λ = − J ξ q ∗ − μ で:
単純相 (y = 0 y=0 y = 0 ):
ξ q ∗ = J q ∫ k 1 μ q − t Δ ( k ) , μ q = μ + ( 1 − 1 q ) J ξ q ∗ \xi^*_q = Jq\int_k \frac{1}{\mu_q - t\Delta(k)}, \quad \mu_q = \mu + (1-\frac{1}{q})J\xi^*_q ξ q ∗ = J q ∫ k μ q − t Δ ( k ) 1 , μ q = μ + ( 1 − q 1 ) J ξ q ∗
q = 1 q=1 q = 1 の場合:Σ 1 = 0 \Sigma_1 = 0 Σ 1 = 0 (ゼロ複雑度)
複雑相 (y > 0 y>0 y > 0 ):
1 = J 2 ∫ k 1 ( μ q − t Δ ( k ) ) 2 + J 4 y 2 1 = J^2\int_k \frac{1}{(\mu_q - t\Delta(k))^2 + J^4 y^2} 1 = J 2 ∫ k ( μ q − t Δ ( k ) ) 2 + J 4 y 2 1 ξ q ∗ = J q ∫ k μ q − t Δ ( k ) ( μ q − t Δ ( k ) ) 2 + J 4 y 2 \xi^*_q = Jq\int_k \frac{\mu_q - t\Delta(k)}{(\mu_q - t\Delta(k))^2 + J^4 y^2} ξ q ∗ = J q ∫ k ( μ q − t Δ ( k ) ) 2 + J 4 y 2 μ q − t Δ ( k )
相境界はLarkin質量により決定される:
1 = J 2 ∫ k 1 ( μ c − t Δ ( k ) ) 2 1 = J^2\int_k \frac{1}{(\mu_c - t\Delta(k))^2} 1 = J 2 ∫ k ( μ c − t Δ ( k ) ) 2 1
成長率:
Σ q = − 1 2 ( ξ q ∗ ) 2 + q f ( ξ q ∗ + μ / J ) − q ∫ k log ( μ − t Δ ( k ) ) \Sigma_q = -\frac{1}{2}(\xi^*_q)^2 + qf(\xi^*_q + \mu/J) - q\int_k \log(\mu - t\Delta(k)) Σ q = − 2 1 ( ξ q ∗ ) 2 + q f ( ξ q ∗ + μ / J ) − q ∫ k log ( μ − t Δ ( k ))
大偏差率関数(ルジャンドル変換を通じて):
Φ ( e ) = max q [ q e − Σ q ] = 1 2 ( ξ q ∗ ) 2 \Phi(e) = \max_q[qe - \Sigma_q] = \frac{1}{2}(\xi^*_q)^2 Φ ( e ) = max q [ q e − Σ q ] = 2 1 ( ξ q ∗ ) 2
ここで e e e と q q q の関係:
e = ∂ q Σ q = f ( ξ q ∗ + μ / J ) − ∫ k log ( μ − t Δ ( k ) ) e = \partial_q\Sigma_q = f(\xi^*_q + \mu/J) - \int_k \log(\mu - t\Delta(k)) e = ∂ q Σ q = f ( ξ q ∗ + μ / J ) − ∫ k log ( μ − t Δ ( k ))
方法の等価性 :2つの方法(異なる極限順序)が同じ結果をもたらすことを証明:リッカチ/DBM法:lim N → ∞ lim L → ∞ \lim_{N\to\infty}\lim_{L\to\infty} lim N → ∞ lim L → ∞ 鞍点法:lim L → ∞ lim N → ∞ \lim_{L\to\infty}\lim_{N\to\infty} lim L → ∞ lim N → ∞ 行列値演算子の場合、局所統計は極限順序に敏感であるため、これは自明ではない ポテンシャル障壁貫通の新機構 :臨界点付近の新しい普遍性クラスを識別標準的なリフシッツテール:exp ( − c ∣ E ∣ 3 / 2 ) \exp(-c|E|^{3/2}) exp ( − c ∣ E ∣ 3/2 ) (N = 1 N=1 N = 1 ) 新しい指数:exp ( − N ( E ∗ − E ) 5 / 4 / J ~ 2 ) \exp(-N(E^*-E)^{5/4}/\tilde{J}^2) exp ( − N ( E ∗ − E ) 5/4 / J ~ 2 ) (大 N N N 、端付近) 臨界点での非従来的な密度消失行動に由来(∝ ( E ∗ − E ) 3 / 2 \propto (E^*-E)^{3/2} ∝ ( E ∗ − E ) 3/2 通常の E ∗ − E \sqrt{E^*-E} E ∗ − E ではなく) ゼロモード支配 :大 N N N 極限では、log det ( H − E ) \log\det(H-E) log det ( H − E ) の揺らぎは主にゼロモード ξ ˉ \bar{\xi} ξ ˉ により支配されることを証明(c > 0 c>0 c > 0 の場合のみ存在)本論文は純粋な理論研究であり、数値実験部分がない。すべての結果は解析計算と漸近解析を通じて得られている。
次元 :d = 0 , 1 d=0, 1 d = 0 , 1 および任意 d d d 行列サイズ :N N N (N → ∞ N\to\infty N → ∞ 極限に焦点)システムサイズ :L L L (または M = L d M=L^d M = L d )無秩序強度 :J J J (または J ~ \tilde{J} J ~ )相関パラメータ :c ≥ 0 c\geq 0 c ≥ 0 (c = 0 c=0 c = 0 または c = 1 c=1 c = 1 )著者は以下の方法で結果の一貫性を検証:
d = 1 d=1 d = 1 での2つの方法の結果比較既知の N = 1 N=1 N = 1 スカラー場合との対比 Allez-Dumaz 1 のDBMに関する結果との対比 特殊点(q = 0 , 1 q=0, 1 q = 0 , 1 など)の物理的意味の検査 これは単一のGOE行列に対角ノイズを加えたものに対応する。谱密度は半円:ρ K ( α ) = 1 2 π J 2 4 J 2 − α 2 \rho_K(\alpha) = \frac{1}{2\pi J^2}\sqrt{4J^2-\alpha^2} ρ K ( α ) = 2 π J 2 1 4 J 2 − α 2
単純相 (μ > ( 2 − q ) J \mu > (2-q)J μ > ( 2 − q ) J ):
ξ q ∗ = 2 J q μ + μ 2 + 4 J 2 ( q − 1 ) \xi^*_q = \frac{2Jq}{\mu + \sqrt{\mu^2 + 4J^2(q-1)}} ξ q ∗ = μ + μ 2 + 4 J 2 ( q − 1 ) 2 J q Σ q = − q ( μ 2 − μ μ 2 + 4 J 2 ( q − 1 ) ) 4 J 2 ( q − 1 ) + q log ( 1 + 1 + 4 J 2 ( q − 1 ) μ 2 ) − q 2 ( 1 + 2 log 2 ) \Sigma_q = -\frac{q(\mu^2 - \mu\sqrt{\mu^2+4J^2(q-1)})}{4J^2(q-1)} + q\log\left(1+\sqrt{1+\frac{4J^2(q-1)}{\mu^2}}\right) - \frac{q}{2}(1+2\log 2) Σ q = − 4 J 2 ( q − 1 ) q ( μ 2 − μ μ 2 + 4 J 2 ( q − 1 ) ) + q log ( 1 + 1 + μ 2 4 J 2 ( q − 1 ) ) − 2 q ( 1 + 2 log 2 )
複雑相 (0 ≤ μ ≤ ( 2 − q ) J 0 \leq \mu \leq (2-q)J 0 ≤ μ ≤ ( 2 − q ) J ):
μ q = μ 2 − q , ξ q ∗ = q 2 − q μ J , y 2 = 1 J 2 ( 1 − 1 ( 2 − q ) 2 μ 2 J 2 ) \mu_q = \frac{\mu}{2-q}, \quad \xi^*_q = \frac{q}{2-q}\frac{\mu}{J}, \quad y^2 = \frac{1}{J^2}\left(1 - \frac{1}{(2-q)^2}\frac{\mu^2}{J^2}\right) μ q = 2 − q μ , ξ q ∗ = 2 − q q J μ , y 2 = J 2 1 ( 1 − ( 2 − q ) 2 1 J 2 μ 2 ) Σ q = q 2 ( 2 − q ) μ 2 J 2 + q log J μ − q 2 \Sigma_q = \frac{q}{2(2-q)}\frac{\mu^2}{J^2} + q\log\frac{J}{\mu} - \frac{q}{2} Σ q = 2 ( 2 − q ) q J 2 μ 2 + q log μ J − 2 q
率関数 :μ < 2 J \mu < 2J μ < 2 J のとき、
Φ ( e ) = μ 2 J 2 φ ( J 2 μ 2 δ e ) , φ ( x ) = 1 + 2 x − 1 + 4 x \Phi(e) = \frac{\mu^2}{J^2}\varphi\left(\frac{J^2}{\mu^2}\delta e\right), \quad \varphi(x) = 1 + 2x - \sqrt{1+4x} Φ ( e ) = J 2 μ 2 φ ( μ 2 J 2 δe ) , φ ( x ) = 1 + 2 x − 1 + 4 x
Larkin質量:μ c = ( J / 2 ) 4 / 3 \mu_c = (J/2)^{4/3} μ c = ( J /2 ) 4/3
相境界:μ b = ( 3 − 2 q ) ( J / 2 ) 4 / 3 \mu_b = (3-2q)(J/2)^{4/3} μ b = ( 3 − 2 q ) ( J /2 ) 4/3
重要な観察 :q → 0 q\to 0 q → 0 のとき、μ b ( q = 0 ) = 3 ( J / 2 ) 4 / 3 = − E ∗ \mu_b(q=0) = 3(J/2)^{4/3} = -E^* μ b ( q = 0 ) = 3 ( J /2 ) 4/3 = − E ∗ 、これはまさにDBM相転移点である!
単純相 (μ > μ b \mu > \mu_b μ > μ b ):μ q \mu_q μ q は立方方程式を満たす
μ q − J 2 2 ( q − 1 ) 1 μ q 1 / 2 = μ \mu_q - \frac{J^2}{2}(q-1)\frac{1}{\mu_q^{1/2}} = \mu μ q − 2 J 2 ( q − 1 ) μ q 1/2 1 = μ
成長率:
Σ q = q ( μ q 1 / 2 − μ 1 / 2 − J 2 ( q − 1 ) 8 μ q ) \Sigma_q = q\left(\mu_q^{1/2} - \mu^{1/2} - \frac{J^2(q-1)}{8\mu_q}\right) Σ q = q ( μ q 1/2 − μ 1/2 − 8 μ q J 2 ( q − 1 ) )
相境界で:
Σ q ∣ μ = μ b = q ( 1 − 3 − 2 q − q − 1 2 ) ( J / 2 ) 2 / 3 \Sigma_q|_{\mu=\mu_b} = q\left(1 - \sqrt{3-2q} - \frac{q-1}{2}\right)(J/2)^{2/3} Σ q ∣ μ = μ b = q ( 1 − 3 − 2 q − 2 q − 1 ) ( J /2 ) 2/3
複雑相 (μ < μ b \mu < \mu_b μ < μ b ):方程式はより複雑で、積分 I 1 ( y ^ ) , I 2 ( y ^ ) I_1(\hat{y}), I_2(\hat{y}) I 1 ( y ^ ) , I 2 ( y ^ ) を含む:
I 1 ( y ^ ) = ∫ 0 ∞ d z 2 π z 1 ( 1 + z ) 2 + y ^ 2 , I 2 ( y ^ ) = ∫ 0 ∞ d z 2 π z 1 + z ( 1 + z ) 2 + y ^ 2 I_1(\hat{y}) = \int_0^\infty \frac{dz}{2\pi\sqrt{z}}\frac{1}{(1+z)^2+\hat{y}^2}, \quad I_2(\hat{y}) = \int_0^\infty \frac{dz}{2\pi\sqrt{z}}\frac{1+z}{(1+z)^2+\hat{y}^2} I 1 ( y ^ ) = ∫ 0 ∞ 2 π z d z ( 1 + z ) 2 + y ^ 2 1 , I 2 ( y ^ ) = ∫ 0 ∞ 2 π z d z ( 1 + z ) 2 + y ^ 2 1 + z
恒等式を利用:
I 2 ( y ^ ) ∓ i y ^ I 1 ( y ^ ) = 1 2 1 ± i y ^ I_2(\hat{y}) \mp i\hat{y}I_1(\hat{y}) = \frac{1}{2\sqrt{1\pm i\hat{y}}} I 2 ( y ^ ) ∓ i y ^ I 1 ( y ^ ) = 2 1 ± i y ^ 1
有効ポテンシャル障壁 (禁閉相、a > a ∗ a > a^* a > a ∗ ):
U = ∫ γ − z a d λ ( λ − z a ) ( γ − − λ ) ( γ + − λ ) U = \int_{\gamma_-}^{z_a} d\lambda (\lambda - z_a)\sqrt{(\gamma_- - \lambda)(\gamma_+ - \lambda)} U = ∫ γ − z a d λ ( λ − z a ) ( γ − − λ ) ( γ + − λ )
ここで z a z_a z a は二重根、γ ± = − z a ± 2 ( a − z a 2 ) \gamma_\pm = -z_a \pm \sqrt{2(a-z_a^2)} γ ± = − z a ± 2 ( a − z a 2 ) は他の根である。
明示的計算は以下を与える:
U = 2 3 a 6 z a 2 − 2 a − β sinh − 1 ( 4 ( − z a ) 3 / 2 β − 2 2 3 / 4 ) U = \frac{2}{3}a\sqrt{6z_a^2 - 2a} - \beta\sinh^{-1}\left(\sqrt{\frac{4(-z_a)^{3/2}}{\sqrt{\beta}}} - \frac{\sqrt{2}}{2^{3/4}}\right) U = 3 2 a 6 z a 2 − 2 a − β sinh − 1 ( β 4 ( − z a ) 3/2 − 2 3/4 2 )
臨界点付近の振る舞い (a → a + ∗ a \to a^*_+ a → a + ∗ ):
U ≃ 4 5 2 × 3 1 / 4 β 1 / 6 ( a − a ∗ ) 5 / 4 U \simeq \frac{4}{5}\sqrt{2\times 3^{1/4}}\beta^{1/6}(a-a^*)^{5/4} U ≃ 5 4 2 × 3 1/4 β 1/6 ( a − a ∗ ) 5/4
対応するエネルギー:E → E − ∗ E \to E^*_- E → E − ∗ U ≃ 4 5 2 × 3 1 / 4 β 1 / 6 ( E ∗ − E ) 5 / 4 U \simeq \frac{4}{5}\sqrt{2\times 3^{1/4}}\beta^{1/6}(E^*-E)^{5/4} U ≃ 5 4 2 × 3 1/4 β 1/6 ( E ∗ − E ) 5/4
大きな負のエネルギー (a = − E → + ∞ a = -E \to +\infty a = − E → + ∞ ):
U ≃ 2 3 a 4 a = 4 3 ∣ E ∣ 3 / 2 U \simeq \frac{2}{3}a\sqrt{4a} = \frac{4}{3}|E|^{3/2} U ≃ 3 2 a 4 a = 3 4 ∣ E ∣ 3/2
これは N = 1 N=1 N = 1 スカラー場合のリフシッツテール指数 3 / 2 3/2 3/2 を回復する。
状態密度推定 :
ρ K ( α ) ∼ τ crossing − 1 ∼ exp ( − N U J ~ 2 ) \rho_K(\alpha) \sim \tau_{\text{crossing}}^{-1} \sim \exp\left(-\frac{NU}{\tilde{J}^2}\right) ρ K ( α ) ∼ τ crossing − 1 ∼ exp ( − J ~ 2 N U )
ここで貫通時間 τ crossing ∼ exp ( N U / J ~ 2 ) \tau_{\text{crossing}} \sim \exp(NU/\tilde{J}^2) τ crossing ∼ exp ( N U / J ~ 2 ) (有効温度 T = J ~ 2 / N T = \tilde{J}^2/N T = J ~ 2 / N )。
新しい普遍性クラス :指数 5 / 4 5/4 5/4 は大 N N N 臨界点付近領域の新しい特徴であり、以下に由来する:臨界点での液滴密度の非従来的な消失:ρ ( λ ) ∝ ( λ − γ − ) 3 / 2 γ + − λ \rho(\lambda) \propto (\lambda-\gamma_-)^{3/2}\sqrt{\gamma_+-\lambda} ρ ( λ ) ∝ ( λ − γ − ) 3/2 γ + − λ 標準的な半円端の 1 / 2 1/2 1/2 指数と異なる 谱端対応 :厳密に証明した
E ∗ = − 3 ( J / 2 ) 4 / 3 = − μ c = α ∗ E^* = -3(J/2)^{4/3} = -\mu_c = \alpha^* E ∗ = − 3 ( J /2 ) 4/3 = − μ c = α ∗
DBM相転移点 = シュレーディンガー演算子谱端流と状態密度 :正確な対応を確立
ρ K ( α ) = 1 π d d E Im J ( E ) ∣ E = α \rho_K(\alpha) = \frac{1}{\pi}\frac{d}{dE}\text{Im}\mathcal{J}(E)\bigg|_{E=\alpha} ρ K ( α ) = π 1 d E d Im J ( E ) E = α N = 1 N=1 N = 1 の振動定理を大 N N N に推広自己平均検証 :2つの方法の一貫性は自己平均特性を独立検証(q = 1 q=1 q = 1 の場合は9,10 により証明済み)複雑度 :q = 1 q=1 q = 1 の場合、Σ 1 \Sigma_1 Σ 1 は無秩序弾性系の臨界点の総複雑度を与え、独立推導を提供古典的ランダム行列 :Hermitian/ユニタリ行列の特性多項式モーメントは広く研究されている25-35 Riemann zeta関数との関連性 Fisher-Hartwig特異点、Toeplitz行列式 ガウス自由場、凍結転移、ガウス乗法カオス 非Hermitian場合 :最近の拡張46-51 疎/帯状構造 :52-55 の研究本論文の貢献 :行列値シュレーディンガー演算子の谱行列式は以前あまり研究されておらず、本論文は系統的に理論を発展させた
一般化されたリャプノフ指数 Λ ( q ) \Lambda(q) Λ ( q ) :56-58 弾性弦(d = 1 d=1 d = 1 ):Fyodorov等11 による任意 q > − 1 q>-1 q > − 1 の研究 本論文は行列値場合(N > 1 N>1 N > 1 )に推広 ヘッシアン行列 :7,8,9,10 はエネルギー汎関数のヘッシアンとしての演算子を研究q = 1 q=1 q = 1 の複雑度 :
8 :N → ∞ N\to\infty N → ∞ 、任意 d d d 10 :厳密な証明本論文:独立推導、一般的な q q q に拡張 標準DBM :Dyson 63 、外部ポテンシャルなしキューブポテンシャル中のDBM :Allez-Dumaz 1 が相転移を発見
本論文は行列シュレーディンガー演算子との正確な対応を確立 有限 N N N のポテンシャル障壁貫通を計算 多チャネル無秩序導体:Dorokhov 2 、Mello-Pereyra-Kumar 3 数学的側面:5,6 インスタントン計算:13 本論文は新しい視点を提供 二重方法論の成功 :行列リッカチ/DBM法(d = 1 d=1 d = 1 、任意 N N N ) 鞍点法(任意 d d d 、大 N N N 極限) 両者は d = 1 d=1 d = 1 、N → ∞ N\to\infty N → ∞ で一貫した結果を与える 深い物理-数学的関連性 :行列シュレーディンガー演算子谱特性 ↔ キューブポテンシャル中のDBM 谱端 ↔ DBM相転移点 状態密度 ↔ 粒子流 リフシッツテール ↔ ポテンシャル障壁貫通 新しい普遍性 :臨界点付近領域の 5 / 4 5/4 5/4 指数は大 N N N 特有の現象完全な大偏差理論 :任意 q q q の成長率 Σ q \Sigma_q Σ q ルジャンドル変換を通じた率関数 Φ ( e ) \Phi(e) Φ ( e ) 相図(単純相vs複雑相) 物理的応用 :無秩序弾性系の複雑度(q = 1 q=1 q = 1 ) 準一次元導体の状態密度テール アンダーソン局在の行列推広 次元制限 :行列リッカチ法は d = 1 d=1 d = 1 のみに適用可能 高次元の場合は鞍点法のみ(大 N N N 極限) パラメータ範囲 :主要結果は q > − 1 q > -1 q > − 1 に対して q ≤ − 1 q \leq -1 q ≤ − 1 は追加の正則化が必要(未扱い)有限 N N N 効果 :ほとんどの明示的結果は N → ∞ N\to\infty N → ∞ 極限 ポテンシャル障壁貫通は有限 N N N 推定を与えるが、単一粒子近似に基づく 数値検証 :純粋な解析研究、数値シミュレーション検証がない 特にポテンシャル障壁貫通公式と 5 / 4 5/4 5/4 指数 複雑相の詳細 :d = 1 d=1 d = 1 複雑相の完全な解析解は複雑文中では方程式形式のみ 有限 N N N 補正 :1 / N 1/N 1/ N 展開の系統的研究単一粒子ポテンシャル障壁貫通近似を超える 集団励起の寄与 数値検証 :行列シュレーディンガー演算子の直接数値シミュレーション DBMシミュレーションによるポテンシャル障壁貫通公式検証 5 / 4 5/4 5/4 指数の普遍性検査推広 :他の対称クラス(GUE、GSE) 非Hermitian推広 より一般的なポテンシャル(立方を超える) 応用 :具体的な物理系の定量的予測 実験データとの比較(介観導体など) ガラス状態系の複雑度景観 数学的厳密性 :一般的な q q q に対する自己平均特性の厳密な証明 鞍点一意性の厳密な論証 大偏差原理の厳密な確立 境界効果 :方法の革新性 :二重アプローチ :2つの独立した方法が相互検証され、結果の信頼性を向上学際的橋梁 :一見無関係な分野(ランダム演算子、DBM、無秩序系)を接続技術的深さ :行列リッカチ法の応用は高度な技巧を示す物理的洞察 :相転移対応 :DBM相転移と谱端の正確な対応は重要な発見新しい普遍性クラス :5 / 4 5/4 5/4 指数は大 N N N 系の新しい物理を明らかにする統一的枠組み :複数の物理問題(アンダーソン局在、弾性系、導体)を統一理論に組み込む数学的厳密性 :推導は詳細で論理的に明確 既知の結果(N = 1 N=1 N = 1 、Allez-Dumaz)との一貫性検査 複数の独立した論証(谱端の3つの識別方法など) 完全性 :任意 d d d (鞍点法)と d = 1 d=1 d = 1 (リッカチ法)をカバー d = 0 d=0 d = 0 (単一粒子)から一般次元への系統的発展相図、大偏差、有限 N N N 効果の包括的分析 執筆品質 :構造が明確で、導入が充分 技術的詳細は付録に、主文は核心思想に焦点 物理的直感と数学的厳密性のバランスが良い 数値検証の欠如 :すべての結果は純粋に解析的で、数値シミュレーション支援がない 特に 5 / 4 5/4 5/4 指数のような新しい予言は数値検証が必要 ポテンシャル障壁貫通公式の単一粒子近似の有効性が未検証 有限 N N N 理論の不完全性 :主要結果は N → ∞ N\to\infty N → ∞ 極限 ポテンシャル障壁貫通は有限 N N N 推定を与えるが、ヒューリスティック論証に基づく 系統的な 1 / N 1/N 1/ N 展開が欠けている 技術的敷居の高さ :ランダム行列理論、確率過程、大偏差理論の深い背景が必要 いくつかのステップ(自己平均特性など)は他の文献の証明に依存 複雑相の完全な解析解は与えられていない 応用議論の不足 :複数の物理的応用が言及されているが、具体的な系の定量分析が欠けている 実験測定可能量との関連が不明確 実際の材料パラメータ下での推定が欠けている 推広の制限 :主にGOE対称性を考慮 c = 0 c=0 c = 0 と c = 1 c=1 c = 1 の2つの場合、中間値は未議論非ガウスノイズへの推広が不明確 理論的貢献 :開拓的 :行列値シュレーディンガー演算子の谱行列式大偏差の初めての系統的研究方法論 :行列リッカチ法は他の問題の研究を刺激する可能性概念的革新 :DBMとランダム演算子の関連性は新しい研究方向を開く学際的影響 :ランダム行列理論 :谱行列式研究を新しいカテゴリーに推広無秩序系 :複雑度計算に新しいツールを提供数学物理 :異なる数学構造(リッカチ方程式、Pastur方程式、DBM)を接続実用的価値 :介観物理への理論的予測を提供 ガラス状態系の複雑度理論 数値アルゴリズムの開発を刺激する可能性(DBM関連性を利用) 再現性 :高い :解析推導は詳細で段階的に検証可能主要な方程式と結果は明確に与えられている ただし数値実装には専門知識が必要 直接応用 :準一次元多チャネル無秩序導体の状態密度計算 無秩序弾性系(弦、膜)の平衡状態統計 アンダーソン局在の行列推広 方法の参考 :他の行列値ランダム演算子の研究 大偏差分析が必要なランダム系 リッカチ方程式を含む問題 概念的推広 :DBM思想を他の物理問題に応用 鞍点-ゼロモード構造を他の系の分析に利用 相転移対応関係の探索 理論発展 :他の対称クラスまたは次元への推広 有限 N N N 理論の発展 レプリカ対称性破れ理論との関連性 1 R. Allez and L. Dumaz. Random matrices in non-confining potentials. J. Stat. Phys. 160(3), 681–714 (2015) - DBMキューブポテンシャル相転移
7 Y.V. Fyodorov, P. Le Doussal. Manifolds pinned by a high-dimensional random landscape: Hessian at the global energy minimum. J. Stat. Phys 179(1), 176–215 (2020) - 大N極限状態密度
8 Y.V. Fyodorov and P. Le Doussal. Manifolds in a high-dimensional random landscape: Complexity of stationary points and depinning. Phys. Rev. E 101(2), 020101 (2020) - 複雑度計算
10 G. Ben Arous, P. Bourgade, and B. McKenna. Landscape complexity beyond invariance and the elastic manifold. Commun. Pure Applied Math. 77(2) 1302–1352 (2024) - q=1厳密証明
11 Y.V. Fyodorov et al. Exponential number of equilibria and depinning threshold for a directed polymer in a random potential. Annals of Physics 397, 1–64 (2018) - N=1場合
62 A. Ossipov. Gelfand-Yaglom formula for functional determinants in higher dimensions. J Phys A 51(49), 495201 (2018) - 行列リッカチ法
総合評価 :これは理論物理学/数学物理学の高品質な論文であり、行列値ランダム演算子の谱理論に重要な貢献をしている。ダイソンブラウン運動との深い関連性の確立を通じて、著者は具体的な技術的問題を解決するだけでなく、異なる数学構造間の内在的な関連性を明らかにした。新発見の 5 / 4 5/4 5/4 指数と相転移対応関係は重要な物理的意義を持つ。主な不足は数値検証と有限 N N N 理論の完全性である。本研究は当該分野の重要な参考文献となり、後続研究を刺激する可能性が高い。