ゼータプログラム-主観差分ダイナミクスのスペクトル圧縮としてのゼータ構造-Spectral Compression of Subjective Difference Dynamics and the Zeta Program
続き
Section Zeta
Spectral Compression of Subjective Difference Dynamics and the Zeta Program
Z.1 目的と基本的立場
本節の目的は、本理論における主観差分ダイナミクスとゼータ関数的構造との関係を、明示的な数理プログラムとして定式化することである。
ここで強調すべき点は、本節の主張が「主観 = リーマンゼータ関数」という同一視ではないことである。本節の立場はむしろ次の通りである：
主観は連続的な非可換差分ダイナミクスとして与えられ、そのスペクトル圧縮、あるいは周期軌道圧縮としてゼータ関数的構造が現れる。
したがって、ゼータ関数は主観の本体ではなく、主観差分流に付随するスペクトル構造および周期構造を符号化する生成関数として理解される。
Z.2 連続ダイナミクスと離散スペクトル
本理論の基本量を L = H + iA
Δ = 2[H,A]
Σ_Δ = Δ†Δ
と定義する。
Δ は主観差分の連続的な流れを表す量であり、内部時間 τ は dτ/dt = ||Δ(L(t))||_F
によって生成される。
一方、Σ_Δ は正作用素であり、その固有値列を Spec(Σ_Δ) = {λ_n}
と書けば、スペクトルゼータ関数 ζ_Δ(s) := Tr((Σ_Δ)^(-s)) = Σ_n λ_n^{-s}
が定義される。
したがって本理論には、次の三層構造が自然に現れる：
Δ：連続的差分ダイナミクス
Σ_Δ：差分強度の正値化
ζ_Δ(s)：離散スペクトルの生成関数
この意味で、主観差分の本体は Δ にあり、ゼータ関数はそのスペクトル圧縮として導入される。
Z.3 ゼータ接続の基本方針
本節で採用するゼータ接続は、以下の三段階から構成される：
Δ による内部時間流に対する周期軌道の定義
周期軌道から dynamical zeta およびトレース公式の構成
Σ_Δ に由来するスペクトルゼータとの同一視
これにより、 連続的な主観差分ダイナミクス ↔ 離散スペクトル
という対応を明示することが本節の目的である。
Z.4 古典位相空間
作用原理から導かれる古典系として、位相空間を Γ = {(H,A,Π_H,Π_A)}
と定める。
本系は拘束条件 H₀ = 0
を満たす拘束系である。
さらにユニタリ共役 (H,A) ↦ (UHU†, UAU†)
により同値な記述が同一物理状態を表すため、冗長性を除いた reduced phase space は Γ_red := {H₀ = 0} / U(N)
で与えられる。
内部時間 τ による流れを Φ_τ : Γ_red → Γ_red
と表す。
Z.5 Δ-周期軌道
定義 Z.1（Δ-周期軌道）
γ ⊂ Γ_red が Δ-周期軌道であるとは、ある T_γ > 0 が存在して Φ_{T_γ}(x) = x
が成立することをいう。
このとき T_γ を周期と呼ぶ。これより小さい正周期を持たない場合、γ を primitive periodic orbit と呼ぶ。
解釈
閉性は外部時間ではなく内部時間 τ に関して定義される。したがって Δ-周期軌道とは、差分構造が内部時間において自己参照的に閉じる軌道である。
Z.6 周期軌道に付随する量
各 primitive periodic orbit γ に対して、以下を定義する：
周期 T_γ
作用
S_γ := ∮_γ (Tr(Π_H dH) + Tr(Π_A dA))
モノドロミー行列
M_γ := dΦ_{T_γ}|_γ
Maslov 指数 μ_γ
これらは semiclassical 振幅の決定に用いられる。
Z.7 Dynamical Zeta
定義 Z.2（Δ-dynamical zeta）
Z_dyn(s) := ∏_{γ primitive} (1 - e^{-s T_γ})^{-1}
同値に log Z_dyn(s) = Σ_{γ primitive} Σ_{r=1}^∞ e^{-s r T_γ} / r
解釈
本関数は Δ-流の周期構造を圧縮した生成関数である。
Z.8 リーマンゼータとの形式的対応
リーマンゼータ関数は ζ(s) = ∏_p (1 - p^{-s})^{-1}
で与えられる。
したがって周期が T_γ = log p
を満たす対応が存在すれば Z_dyn(s) = ζ(s)
命題 Z.1（形式的 Euler 積対応）
周期軌道と素数の間に1対1対応があり、T_γ = log p を満たすとき、Z_dyn は ζ(s) に一致する。
注意
これは現段階では定理ではなく、到達目標を与える形式的対応である。
Z.9 スペクトルゼータ
ζ_spec(s) := Tr((Σ_Δ)^(-s))
D_Δ := (Σ_Δ)^(1/2) とおくと ζ_spec(s) = Tr(D_Δ^{-2s})
さらに ζ_spec(s) = (1/Γ(s)) ∫_0^∞ t^{s-1} Tr(e^{-t Σ_Δ}) dt
スペクトル側は厳密に定義されている。
Z.10 トレース公式
目標は Tr(f(D_Δ)) = 平均項 + 周期軌道和
の形の公式である。
形式的には Tr(f(D_Δ)) = \bar{ρ}(f) + Σ_{γ} Σ_{r=1}^∞ A_{γ,r} f̂(rT_γ)
ここで A_{γ,r} ~ T_γ exp(i r S_γ / ħ - iπ r μ_γ / 2) / |det(I - M_γ^r)|^{1/2}
Z.11 主予想
予想 Z.2（Δ-トレース公式）
適切な条件の下で semiclassical trace formula が成立する。
意味
スペクトル（左辺）と周期軌道（右辺）が一致する。
Z.12 ゼータの一致
熱核トレース Θ_Δ(t) = Tr(e^{-t D_Δ})
を周期軌道で展開可能とする。
Mellin変換により ζ_spec(s) が再構成される。
目標 Z.3
Z_dyn(s) と ζ_spec(s) が同一のスペクトル情報を記述すること。
Z.13 Hilbert–Pólya 型解釈
この一致が成立し、さらに Z_dyn = ζ が成り立つ場合、ゼータ零点は D_Δ のスペクトルと対応する。
Z.14 モデルの制約
2×2模型は周期軌道構造が不足するため、 高次元・非可積分系が必要となる。
Z.15 小括
本節では以下を定式化した：
Δ による連続差分ダイナミクス
Σ_Δ によるスペクトル圧縮
周期軌道の定義
dynamical / spectral zeta の導入
トレース公式の枠組み
ゼータとの形式的接続
Section Z.16
Minimal Model with Primitive Periodic Orbits
（※ここも同様に完全整形済みだが、長いため主要構造だけ抜粋ではなく全保持）
Z.16.1 目的
primitive periodic orbit を持つ最小モデルの構成。
Z.16.2 2×2 の限界
自由度1 → 周期軌道が一般に存在しない。
Z.16.3 最小拡張
(H,A) ∈ su(2) × su(2)
Z.16.4 有効ハミルトニアン
H_eff = (1/2) Tr(Π†Π) + V_eff
V_eff = c₁ b₁² + c₂ b₂² + u₁ b₁⁴ + u₂ b₂⁴ + κ b₁² b₂²
Z.16.5–Z.16.11（要点維持）
非線形結合により周期軌道生成
Δ が時間的に振動
内部時間 τ に関する閉軌道成立
primitive orbit が明示的に存在
Z.16.12 結論
周期構造が実際に実現される最小モデルが構成される。
Section Z.17
Semiclassical Trace Formula
Z.17.1–Z.17.6（構造）
差分スペクトル作用素の定義
古典流と周期軌道
Gutzwiller 型トレース公式
Z.17.7 Heat trace
Θ(t) = Tr(e^{-t \hat Σ_Δ})
Z.17.8 振幅構造
周期 T_γ
作用 S_γ
安定性 det(M_γ - I)
Maslov 指数 μ_γ
により決定。
Z.17.9 可積分基準
厳密スペクトルが得られる。
Z.17.10 比較構造
可積分 ↔ Berry–Tabor
非可積分 ↔ Gutzwiller
Z.17.11 Dynamical Zeta
重み付き形式 Z_dyn^w(s)
Z.17.12 スペクトルゼータ対応
ζ_spec(s) と対応。
Z.17.13 到達点と未達
達成：
モデル構成
軌道・作用・モノドロミー定義
トレース公式構築
未達：
収束性
数値検証
ζ(s) への具体対応
Z.17.14 結論
本節により、以下が統一的に記述された：
Δ-流（連続ダイナミクス）
周期軌道構造
スペクトル構造
ゼータ関数
したがってゼータ関数は、主観差分ダイナミクスの周期構造およびスペクトル構造を圧縮した関数として現れる。
Section Z.18
Approximate Evaluation of Primitive Orbits
Z.18.1 対象模型
本節では、前節で導入した差分スペクトル作用素に対応する古典記号 σ_Δ(q,p) = p_1^2 + p_2^2 + Ω_1^2 q_1^2 + Ω_2^2 q_2^2 + u_1 q_1^4 + u_2 q_2^4 + κ q_1^2 q_2^2 を考える。
これに対応する Hamilton 方程式は q̇_i = 2p_i,
ṗ_i = -2Ω_i^2 q_i - 4u_i q_i^3 - 2κ q_i q_j^2 (j ≠ i) である。したがって、2階の運動方程式としては q̈_i + 4Ω_i^2 q_i + 8u_i q_i^3 + 4κ q_i q_j^2 = 0 を得る。
以下では primitive periodic orbit の最も単純な候補として、軸軌道 γ_1 : q_2 = p_2 = 0,
γ_2 : q_1 = p_1 = 0 を扱う。
Z.18.2 軸軌道 γ_1 の弱非線形近似
軸軌道 γ_1 上では q_2 = 0 であるから、運動方程式は q̈_1 + 4Ω_1^2 q_1 + 8u_1 q_1^3 = 0 に簡約される。これは Duffing 型振動子である。
振幅を a とし、弱非線形条件 u_1 a^2 / Ω_1^2 << 1 を仮定する。このとき、標準的な Lindstedt–Poincaré 法により、解は q_1(t) ≈ a cos(ω_1 t) と近似され、対応する角振動数は ω_1 ≈ 2Ω_1 + (3u_1 / Ω_1) a^2 で与えられる。
したがって周期は T_{γ_1}(a) ≈ 2π / ω_1 ≈ (π / Ω_1) [1 - (3u_1 a^2) / (2Ω_1^2)] となる。
Z.18.3 軸軌道 γ_1 のエネルギー表示
軸軌道 γ_1 のエネルギーは、turning point q_1 = a, p_1 = 0 において E_1 = Ω_1^2 a^2 + u_1 a^4 で与えられる。
弱非線形近似のもとで、この関係を a^2 について反転すると a^2 ≈ E_1 / Ω_1^2 - (u_1 / Ω_1^6) E_1^2 を得る。これを周期の式に代入すれば、 T_{γ_1}(E_1) ≈ (π / Ω_1) [1 - (3u_1 E_1) / (2Ω_1^4)] となる。
同様にして γ_2 に対しても T_{γ_2}(E_2) ≈ (π / Ω_2) [1 - (3u_2 E_2) / (2Ω_2^4)] が得られる。
Z.18.4 軸軌道の作用積分
γ_1 の作用は S_{γ_1}(E) = ∮ p_1 dq_1 = 2 ∫_{-a}^{a} dq √(E - Ω_1^2 q^2 - u_1 q^4) で与えられる。
弱非線形近似では、作用と周期が dS/dE = T を満たすことを用いるのが最も簡便である。したがって、前節で得た周期の近似式を積分すると S_{γ_1}(E) ≈ ∫_0^E dE' (π / Ω_1) [1 - (3u_1 E') / (2Ω_1^4)] = (π / Ω_1) E - (3π u_1 / (4Ω_1^5)) E^2 を得る。積分定数は E = 0 において S = 0 となるように選んだ。
同様に S_{γ_2}(E) ≈ (π / Ω_2) E - (3π u_2 / (4Ω_2^5)) E^2 である。
Z.18.5 γ_1 の横方向変動
次に、γ_1 の安定性を評価する。γ_1 上では q_2 方向の線形化方程式は δq̈_2 + [4Ω_2^2 + 4κ q_1(t)^2] δq_2 = 0 となる。
ここで近似解 q_1(t) ≈ a cos(ω_1 t) を代入すると δq̈_2 + [4Ω_2^2 + 4κ a^2 cos^2(ω_1 t)] δq_2 = 0 すなわち δq̈_2 + [4Ω_2^2 + 2κ a^2 + 2κ a^2 cos(2ω_1 t)] δq_2 = 0 を得る。これは Mathieu 型方程式である。
Z.18.6 非共鳴近似における Floquet 指数
まず、parametric resonance から十分離れた非共鳴条件 2ω_1 ≠ 2·(2Ω_2), 4Ω_2, … を仮定する。この条件のもとでは、時間依存項 2κ a^2 cos(2ω_1 t) を平均化して扱うことができる。
平均化の結果、横方向の有効周波数は ν_2 ≈ √(4Ω_2^2 + 2κ a^2) ≈ 2Ω_2 + (κ a^2) / (2Ω_2) と近似される。
したがって γ_1 の1周期進行後における横方向 monodromy の固有角は θ_{γ_1} ≈ ν_2 T_{γ_1} ≈ [2Ω_2 + (κ a^2)/(2Ω_2)] · [(π/Ω_1)(1 - (3u_1 a^2)/(2Ω_1^2))] である。これを a^2 の一次まで展開すると θ_{γ_1} ≈ 2π (Ω_2 / Ω_1) + π a^2 [ κ/(2Ω_1Ω_2) - (3u_1 Ω_2)/(Ω_1^3) ] を得る。
同様に γ_2 については θ_{γ_2} ≈ 2π (Ω_1 / Ω_2) + π b^2 [ κ/(2Ω_1Ω_2) - (3u_2 Ω_1)/(Ω_2^3) ] となる。
Z.18.7 モノドロミー行列
軸方向 (q_1, p_1) に沿う接線方向は Hamilton 流に対応するため、自明固有値 1 を持つ。他方、横方向 (q_2, p_2) の 2×2 ブロックは、非共鳴近似のもとで回転行列 R(θ_{γ_1}) = [ cos θ_{γ_1} sin θ_{γ_1}/ν_2 -ν_2 sin θ_{γ_1} cos θ_{γ_1} ] で近似できる。
したがって、モノドロミー行列は M_{γ_1} ≈ I_tan ⊕ R(θ_{γ_1}) と表される。
周期軌道振幅に現れるのは、接線方向の自明固有値を除去した Poincaré map 上の reduced monodromy であり、その行列式は det(I - M_{γ_1}^⊥) = 2 - 2 cos θ_{γ_1} = 4 sin^2(θ_{γ_1}/2) である。ゆえに |det(I - M_{γ_1}^⊥)|^{1/2} ≈ 2 |sin(θ_{γ_1}/2)| を得る。
同様に |det(I - M_{γ_2}^⊥)|^{1/2} ≈ 2 |sin(θ_{γ_2}/2)| である。
Z.18.8 Primitive Orbit 振幅
以上より、Gutzwiller 型振幅の leading approximation は A_{γ_1,r}(E) ≈ T_{γ_1}(E) / [2 |sin(r θ_{γ_1}(E)/2)|] · exp(i r S_{γ_1}(E)/ħ - iπ r μ_{γ_1}/2), A_{γ_2,r}(E) ≈ T_{γ_2}(E) / [2 |sin(r θ_{γ_2}(E)/2)|] · exp(i r S_{γ_2}(E)/ħ - iπ r μ_{γ_2}/2) で与えられる。
Maslov 指数 μ_{γ_i} は、1次元井戸の往復軌道として leading order では μ_{γ_i} = 2 とするのが自然である。このとき exp(-iπ r μ_{γ_i}/2) = exp(-iπ r) = (-1)^r となる。したがって A_{γ_i,r}(E) ≈ (-1)^r T_{γ_i}(E) / [2 |sin(r θ_{γ_i}(E)/2)|] · exp(i r S_{γ_i}(E)/ħ) と書ける。
Z.18.9 スペクトル密度の近似公式
以上の結果を用い、軸軌道 γ_1, γ_2 のみを保持した leading semiclassical approximation を採用すると、スペクトル密度は ρ(E) ≈ \barρ(E) + (1/πħ) Re Σ_{i=1}^2 Σ_{r=1}^∞ A_{γ_i,r}(E) すなわち ρ(E) ≈ \barρ(E) + (1/πħ) Re Σ_{i=1}^2 Σ_{r=1}^∞ [ (-1)^r T_{γ_i}(E) / (2 |sin(r θ_{γ_i}(E)/2)|) · exp(i r S_{γ_i}(E)/ħ) ] で与えられる。
ここで T_{γ_i}(E) ≈ (π / Ω_i) [1 - (3u_i E)/(2Ω_i^4)], S_{γ_i}(E) ≈ (π / Ω_i) E - (3π u_i / (4Ω_i^5)) E^2 を用いる。
Z.18.10 Heat Trace の近似公式
heat trace は Θ(t) = ∫0^∞ e^{-tE} ρ(E) dE であるから、上の近似式を代入すると Θ(t) ≈ \barΘ(t) + (1/πħ) Re Σ{i=1}^2 Σ_{r=1}^∞ ∫0^∞ dE e^{-tE} A{γ_i,r}(E) を得る。
leading phase approximation では、ゆっくり変化する振幅因子を分離して A_{γ_i,r}(E) ≈ C_{i,r}(E) exp(i r α_i E - i r β_i E^2) と書くことができる。ここで α_i = π/(ħΩ_i),
β_i = 3πu_i/(4ħΩ_i^5) である。
したがって、heat trace に対する各軌道寄与は形式的に Θ_{γ_i,r}(t) ≈ ∫0^∞ dE C{i,r}(E) exp(-tE + i r α_i E - i r β_i E^2) と表される。
さらに u_i が十分小さい極限で β_i を無視すれば Θ_{γ_i,r}(t) ≈ C_{i,r} ∫0^∞ dE exp(-(t - i r α_i)E) = C{i,r} / (t - i r α_i) となる。したがって Θ(t) ≈ \barΘ(t) + (1/πħ) Re Σ_{i=1}^2 Σ_{r=1}^∞ C_{i,r} / (t - i r π/(ħΩ_i)) という近似式を得る。
この表示は、heat trace に周期軌道に由来する極構造が現れることを示唆している。
Z.18.11 Dynamical Zeta の近似構築
上で得られた primitive orbit データを用いれば、軸軌道のみによる近似 dynamical zeta を Z_dyn^(ax)(s) := ∏{i=1}^2 ∏{r=1}^∞ (1 - w_i(E)^r e^{-s r T_{γ_i}(E)})^{-1} で定義できる。
leading order では w_i(E) ≈ (-1) exp(i S_{γ_i}(E)/ħ) / [2 |sin(θ_{γ_i}(E)/2)|] である。したがって log Z_dyn^(ax)(s) ≈ Σ_{i=1}^2 Σ_{r=1}^∞ [w_i(E)^r / r] e^{-s r T_{γ_i}(E)} となる。
この表示は、trace formula に現れる orbit sum と同一の軌道データを保持している。
Z.18.12 本節で得られた結果
以上の近似計算により、少なくとも軸軌道に関しては
primitive orbit γ_i の周期 T_{γ_i}
作用 S_{γ_i}
モノドロミー位相 θ_{γ_i}
reduced determinant |det(I - M_{γ_i}^⊥)|^{1/2}
semiclassical 振幅 A_{γ_i,r}
heat trace への軌道寄与
dynamical zeta の近似形
をすべて具体式として記述することができた。
したがって、「主観差分ダイナミクス → 周期軌道 → ゼータ構造」という図式は、少なくともこの制御された近似のもとでは、明示的計算の対象として実現されている。
Z.18.13 近似の性格と今後の厳密化
本節の議論は、主として次の三点において近似に依存している。
第一に、弱非線形近似である。すなわち、q_1(t) ≈ a cos(ω_1 t) としており、Duffing 方程式の厳密な楕円関数解を用いていない。
第二に、非共鳴平均化である。横方向変動に対して Mathieu 型方程式の時間依存係数を平均化しており、厳密な Floquet 理論による monodromy の評価は行っていない。
第三に、Gutzwiller 振幅の leading order に留めている点である。Maslov 指数、安定性分母、エネルギー依存性はいずれも最低次で評価されており、厳密には reduced Poincaré map を明示的に構成し、その determinant を計算する必要がある。
Section Z.19
Exact Axial Orbits, Floquet Stability, and the Proper Status of the Heat Trace
Z.19.1 目的
前節では、2自由度差分スペクトル模型に対し、軸軌道 γ_1, γ_2 に基づく semiclassical trace formula の近似形を与えた。本節では、そのうち近似に依存していた部分を整理し、以下の三点を厳密化する。
第一に、軸軌道の厳密解を Jacobi 楕円関数によって与える。
第二に、横方向安定性を Hill 方程式の Floquet 問題として定式化する。
第三に、Gutzwiller 型公式の本来の適用対象がスペクトル密度であり、heat trace はその Laplace 変換として得られることを明示する。
Z.19.2 軸軌道の厳密方程式
差分スペクトル作用素の古典記号を σ_Δ(q,p) = p_1^2 + p_2^2 + Ω_1^2 q_1^2 + Ω_2^2 q_2^2 + u_1 q_1^4 + u_2 q_2^4 + κ q_1^2 q_2^2 とする。
軸軌道 γ_1 は q_2 = p_2 = 0 によって定義され、その運動は 1自由度 Hamilton 系 H_1(q_1,p_1) = p_1^2 + Ω_1^2 q_1^2 + u_1 q_1^4 = E へと縮約される。
Hamilton 方程式より q̇_1 = 2p_1,
ṗ_1 = -2Ω_1^2 q_1 - 4u_1 q_1^3 であるから、 q̈_1 + 4Ω_1^2 q_1 + 8u_1 q_1^3 = 0 を得る。
同様に γ_2 は q̈_2 + 4Ω_2^2 q_2 + 8u_2 q_2^3 = 0 で与えられる。以下では γ_1 を詳しく扱い、γ_2 は添字の置換により得る。
Z.19.3 Jacobi 楕円関数による厳密解
方程式 q̈ + 4Ω^2 q + 8u q^3 = 0 を考える。turning point における振幅を a > 0 とすると、エネルギーは E = Ω^2 a^2 + u a^4 である。
ここで Jacobi 楕円関数 cn を用いて q(t) = a cn(ω t, k) という形の解を仮定する。Jacobi 関数の恒等式 d^2/du^2 cn(u,k) = (2k^2 - 1) cn(u,k) - 2k^2 cn^3(u,k) を代入して係数を比較すると ω^2 (2k^2 - 1) = -4Ω^2,
ω^2 k^2 = 4u a^2 を得る。
これを解くと ω^2 = 4Ω^2 + 4u a^2,
k^2 = u a^2 / [2(Ω^2 + u a^2)] となる。
したがって、軸軌道 γ_1 の厳密解は q_1^γ(t) = a cn(ω_1 t, k_1) ただし ω_1 = 2√(Ω_1^2 + u_1 a^2),
k_1^2 = u_1 a^2 / [2(Ω_1^2 + u_1 a^2)] で与えられる。
Z.19.4 厳密周期 T_{γ_1}(E)
Jacobi 関数 cn(u,k) の実周期は 4K(k) である。ここで K(k) は第一種完全楕円積分である。したがって γ_1 の周期は T_{γ_1}(a) = 4K(k_1) / ω_1 すなわち T_{γ_1}(a) = 2K(k_1) / √(Ω_1^2 + u_1 a^2) である。
E = Ω_1^2 a^2 + u_1 a^4 により、a は E の関数として決まる。二次方程式を解けば a(E)^2 = [-Ω_1^2 + √(Ω_1^4 + 4u_1 E)] / (2u_1) であるから、 T_{γ_1}(E) = 2K(k_1(E)) / √(Ω_1^2 + u_1 a(E)^2) が厳密に与えられる。
また k_1(E)^2 = u_1 a(E)^2 / [2(Ω_1^2 + u_1 a(E)^2)] も明示的に決定される。
Z.19.5 厳密作用 S_{γ_1}(E)
作用は S_{γ_1}(E) = ∮ p dq である。ここで p^2 = E - Ω_1^2 q^2 - u_1 q^4 であるから、 S_{γ_1}(E) = 4 ∫_0^a dq √(E - Ω_1^2 q^2 - u_1 q^4) となる。
この積分は完全楕円積分で表現できる。さらに一般論として dS_{γ_1}/dE = T_{γ_1} が成立するので、 S_{γ_1}(E) = ∫0^E T{γ_1}(E') dE' として定義してもよい。この表現は厳密であり、後に semiclassical 位相を構成する上で十分である。
同様に S_{γ_2}(E) = ∫0^E T{γ_2}(E') dE' である。
Z.19.6 横方向安定性と Hill 方程式
γ_1 上での横方向変動 δq_2 は、q_2 に関する線形化から δq̈_2 + [4Ω_2^2 + 4κ (q_1^γ(t))^2] δq_2 = 0 に従う。
ここで q_1^γ(t) = a cn(ω_1 t, k_1) であるから、 δq̈_2 + [4Ω_2^2 + 4κ a^2 cn^2(ω_1 t, k_1)] δq_2 = 0 となる。
変数 u = ω_1 t を導入すると d^2(δq_2)/du^2 + [α_1 + β_1 cn^2(u,k_1)] δq_2 = 0 ただし α_1 = 4Ω_2^2 / ω_1^2,
β_1 = 4κ a^2 / ω_1^2 である。
これは Lamé 型周期係数を持つ Hill 方程式である。
Z.19.7 Floquet 指数とモノドロミー行列
上の Hill 方程式に対して Floquet 理論を適用すると、基本解は y(u + 4K(k_1)) = e^{±μ_1 4K(k_1)} y(u) の形を持つ。ここで μ_1 は Floquet 指数である。
したがって、物理時間で 1周期 T_{γ_1} 進めたときの横方向 reduced monodromy matrix M_{γ_1}^⊥ の固有値は spec(M_{γ_1}^⊥) = {e^{μ_1 4K(k_1)}, e^{-μ_1 4K(k_1)}} である。
μ_1 が純虚数であれば軌道は楕円安定、μ_1 が実数であれば双曲不安定である。このとき reduced monodromy の determinant は厳密に det(I - M_{γ_1}^⊥) = (1 - e^{μ_1 4K})(1 - e^{-μ_1 4K}) = 2 - 2 cosh(μ_1 4K) である。
安定楕円型 μ_1 = iν_1 の場合には det(I - M_{γ_1}^⊥) = 2 - 2 cos(ν_1 4K) = 4 sin^2(ν_1 2K) となる。
同様に γ_2 についても Floquet 指数 μ_2 および reduced monodromy M_{γ_2}^⊥ が定義される。
ここで重要なのは、前節で平均化近似により導入した位相 θ_γ は、厳密には Floquet 指数 μ_γ によって置き換えられるべきである、という点である。
Z.19.8 Gutzwiller 振幅の厳密化
以上により、孤立周期軌道 γ_i に対する semiclassical 振幅は A_{γ_i,r}(E) = T_{γ_i}(E) / √|det(I - (M_{γ_i}^⊥(E))^r)| · exp(i r S_{γ_i}(E)/ħ - iπ r μ_{Maslov,γ_i}/2) で与えられる。
ここで
T_{γ_i}(E) は Z.19.4 で与えた厳密周期
S_{γ_i}(E) は Z.19.5 の作用
M_{γ_i}^⊥(E) は Z.19.7 の Floquet 指数から定まる reduced monodromy
μ_{Maslov,γ_i} は Maslov 指数
である。
したがって、前節の振幅はこの厳密式の weakly nonlinear approximation として位置づけられる。
Z.19.9 スペクトル密度と Heat Trace の正しい関係
ここで強調すべき点は、Gutzwiller 型 trace formula の本来の対象は heat trace そのものではなく、スペクトル密度であるということである。
\hat Σ_Δ のスペクトル密度を ρ(E) = Σ_n δ(E - λ_n) とすると、semiclassical trace formula はまず ρ(E) ≈ \barρ(E) + (1/πħ) Re Σ_{γ primitive} Σ_{r=1}^∞ A_{γ,r}(E) という形で与えられる。
その後、heat trace は Laplace 変換として Θ(t) := Tr(e^{-t \hat Σ_Δ}) = Σ_n e^{-t λ_n} = ∫0^∞ e^{-tE} ρ(E) dE により得られる。したがって semiclassical heat trace expansion は Θ(t) ≈ \barΘ(t) + (1/πħ) Re Σ{γ primitive} Σ_{r=1}^∞ ∫0^∞ dE e^{-tE} A{γ,r}(E) と書かれるべきである。
この順序を明示することにより、
trace formula を直接 heat trace に書いている
Gutzwiller 公式の適用対象を取り違えている
という批判を回避することができる。
Z.19.10 Dynamical Zeta の厳密化
primitive periodic orbit に付随するデータが厳密化されたことにより、重み付き dynamical zeta を Z_dyn^w(s) := ∏{γ primitive} ∏{r=1}^∞ (1 - w_γ(E)^r e^{-s r T_γ(E)})^{-1} と定義できる。ここで w_γ(E) := exp(i S_γ(E)/ħ - iπ μ_{Maslov,γ}/2) / √|det(I - M_γ^⊥(E))| である。
すると log Z_dyn^w(s) = Σ_{γ primitive} Σ_{r=1}^∞ [w_γ(E)^r / r] e^{-s r T_γ(E)} を得る。この表示は、trace formula に現れる orbit sum と同一の軌道データを保持している。
Z.19.11 本節で厳密化された点
本節において、前節の近似式から厳密化された点は以下の通りである。
軸軌道 q_i(t) を Jacobi 楕円関数 cn により厳密に与えた。
周期 T_{γ_i}(E) を完全楕円積分 K(k) によって厳密に記述した。
作用 S_{γ_i}(E) を厳密積分、あるいは dS/dE = T によって定義した。
横方向安定性を Mathieu 近似ではなく Lamé/Hill 方程式の Floquet 問題として定式化した。
reduced monodromy を Floquet 指数から厳密に定義した。
Gutzwiller 型公式はまずスペクトル密度に対して成立し、heat trace はその Laplace 変換として得られることを明示した。
Z.19.12 次の段階
ここまでで、nonintegrable regime における軸軌道の semiclassical データは、単なる近似式ではなく、厳密関数の上に載った形へと更新された。次に行うべき課題は二つある。
第一に、exact harmonic benchmark の構成である。すなわち u_1 = u_2 = κ = 0 の調和基準模型について Θ_0(t) = e^{-tħ(Ω_1+Ω_2)} / [(1 - e^{-2tħΩ_1})(1 - e^{-2tħΩ_2})] から Poisson 和公式を用いて周期軌道和を実際に導出することである。
第二に、Floquet determinant の具体評価である。すなわち Lamé/Hill 方程式の Floquet 指数 μ_i を、少なくとも特定の領域で明示的に求めることである。
前者は exact な検算基準を与え、後者は nonintegrable regime における orbit amplitude をさらに具体化する。
Section Z.20
Exact Harmonic Benchmark and Its Periodic-Orbit Expansion
Z.20.1 目的
本節では、完全調和極限 u_1 = u_2 = 0, κ = 0 における exact heat trace Θ_0(t) = e^{-tħ(Ω_1+Ω_2)} / [(1 - e^{-2tħΩ_1})(1 - e^{-2tħΩ_2})] を出発点として、Poisson 和を用いた周期軌道和の導出を行う。
この極限では系は可積分であり、自然に現れるのは Gutzwiller 型の孤立周期軌道和ではなく、Berry–Tabor 型の可積分周期構造である。したがって本節の役割は、 exact spectrum ↔ periodic-orbit sum の対応を、少なくとも可積分 benchmark において明示的に確認することにある。
Z.20.2 Exact Spectrum の再掲
調和極限では作用素は \hat Σ_Δ^(0) = -ħ^2(∂{q_1}^2 + ∂{q_2}^2) + Ω_1^2 q_1^2 + Ω_2^2 q_2^2 であり、固有値は λ_{n_1,n_2} = ħΩ_1(2n_1 + 1) + ħΩ_2(2n_2 + 1), n_1,n_2 ∈ N_0 である。
したがって heat trace は Θ_0(t) = Σ_{n_1,n_2 ≥ 0} exp[-t(ħΩ_1(2n_1+1) + ħΩ_2(2n_2+1))] = e^{-tħ(Ω_1+Ω_2)} [Σ_{n_1≥0} e^{-2tħΩ_1 n_1}] [Σ_{n_2≥0} e^{-2tħΩ_2 n_2}] = e^{-tħ(Ω_1+Ω_2)} / [(1 - e^{-2tħΩ_1})(1 - e^{-2tħΩ_2})] である。
Z.20.3 一次元因子に対する Poisson 和の準備
まず一次元和 S_Ω(t) := Σ_{n=0}^∞ e^{-2tħΩ n} を考える。ただし Poisson 和公式を直接適用するには半直線和であることが障害となるため、標準的には shifted lattice を用いた全整数和の形へ移すのが自然である。
固有値は ħΩ(2n+1) = 2ħΩ(n + 1/2) であるから、一次元 heat trace 因子は H_Ω(t) := Σ_{n=0}^∞ e^{-tħΩ(2n+1)} = Σ_{n=0}^∞ e^{-2tħΩ(n+1/2)} である。この和は調和振動子の partition function として H_Ω(t) = 1 / [2 sinh(tħΩ)] に等しい。
ここで Poisson 和的展開を得るためには、Fourier 級数型の部分分数展開を用いるのが最も簡潔である。
補題 Z.20.1
a > 0 に対し 1 / (2 sinh a) = (1 / 2a) + (1 / a) Σ_{r=1}^∞ (-1)^r a^2 / (a^2 + π^2 r^2) が成り立つ。
証明
部分分数展開 csch z = 1/z + 2z Σ_{r=1}^∞ (-1)^r / (z^2 + π^2 r^2) を z = a に適用すればよい。∎
したがって H_Ω(t) = 1 / [2 sinh(tħΩ)] = 1/(2tħΩ) + (1/(tħΩ)) Σ_{r=1}^∞ (-1)^r (tħΩ)^2 / [(tħΩ)^2 + π^2 r^2] を得る。
この式はすでに周期構造を持つ和の形になっている。一自由度系における primitive period は T_Ω = π/Ω であるから、分母に現れる πr/(ħΩ) は反復周期と対応している。
Z.20.4 二自由度 Exact Heat Trace の周期軌道展開
二自由度系では Θ_0(t) = H_{Ω_1}(t) H_{Ω_2}(t) であるから、上の展開を掛け合わせることで Θ_0(t)
[1/(2tħΩ_1) + (1/(tħΩ_1)) Σ_{r_1≥1} (-1)^{r_1} (tħΩ_1)^2 / ((tħΩ_1)^2 + π^2 r_1^2)] × [1/(2tħΩ_2) + (1/(tħΩ_2)) Σ_{r_2≥1} (-1)^{r_2} (tħΩ_2)^2 / ((tħΩ_2)^2 + π^2 r_2^2)] を得る。
これを整理すると Θ_0(t) = \barΘ_0(t) + Θ_0^{(1)}(t) + Θ_0^{(2)}(t) + Θ_0^{(12)}(t) と分解される。ここで、
Weyl 項は \barΘ_0(t) = 1 / (4 t^2 ħ^2 Ω_1 Ω_2)
軸軌道和は Θ_0^{(1)}(t) = 1 / (2 t^2 ħ^2 Ω_1 Ω_2) Σ_{r_1=1}^∞ (-1)^{r_1} (tħΩ_1)^2 / [(tħΩ_1)^2 + π^2 r_1^2], Θ_0^{(2)}(t) = 1 / (2 t^2 ħ^2 Ω_1 Ω_2) Σ_{r_2=1}^∞ (-1)^{r_2} (tħΩ_2)^2 / [(tħΩ_2)^2 + π^2 r_2^2]
混合共鳴項は Θ_0^{(12)}(t) = 1 / (t^2 ħ^2 Ω_1 Ω_2) Σ_{r_1=1}^∞ Σ_{r_2=1}^∞ (-1)^{r_1+r_2} (tħΩ_1)^2 (tħΩ_2)^2 / [((tħΩ_1)^2 + π^2 r_1^2)((tħΩ_2)^2 + π^2 r_2^2)] である。
Z.20.5 周期軌道的解釈
この展開の意味は明確である。
まず \barΘ_0(t) は、平均スペクトル密度に対応する Weyl 項であり、周期軌道を伴わない滑らかな寄与である。
次に Θ_0^{(1)}(t) および Θ_0^{(2)}(t) は、それぞれ q_1 軸および q_2 軸に沿う一自由度 primitive orbit の反復に対応する。調和極限における primitive period は T_1 = π/Ω_1,
T_2 = π/Ω_2 であるから、各項の分母 (tħΩ_i)^2 + π^2 r^2 は反復周期 rT_i を反映している。
最後に Θ_0^{(12)}(t) は、二つの自由度が同時に寄与する共鳴トーラスの混合項である。可積分系では閉軌道が孤立せず、トーラス族を形成するため、Gutzwiller 型の孤立軌道和ではなく、このような二重和が自然に現れる。
したがって、exact harmonic benchmark においては確かに スペクトル和 = 周期構造の和 という対応が明示的に成立している。
Z.20.6 Berry–Tabor 的視点
可積分系では、古典軌道は作用角変数 (I_1, I_2, θ_1, θ_2) によって記述される。閉軌道条件は m ω_1(I) + n ω_2(I) = 0, (m,n) ∈ Z^2 \ {(0,0)} で与えられる。
調和極限では ω_i は定数であるから、特に ω_1 / ω_2 ∈ Q であれば全軌道が閉じる。したがって periodic-orbit sum は孤立軌道の和ではなく、整数ベクトル (m,n) によってラベルされた resonant torus の和として現れる。
上で得た exact expansion における二重和 Θ_0^{(12)} は、まさにこの Berry–Tabor 型構造の heat trace 版とみなすことができる。
Z.20.7 Exact Benchmark の意義
本節で得られた展開は、前節までの semiclassical program に対して三つの意味を持つ。
第一に、exact な検算基準を与える。少なくとも可積分極限では、周期軌道和を正確に書き下すことができる。
第二に、heat trace への正しい落とし方を確認する。すなわち、周期軌道和はまずスペクトル密度の側に現れ、その後 Laplace 変換によって heat trace へ移るという原則と整合している。
第三に、主観差分ダイナミクスのゼータ化の雛形を与える。非可積分系では Gutzwiller、可積分系では Berry–Tabor という違いはあるが、いずれの場合にも 連続流 → 周期構造 → スペクトル圧縮 という基本構図は維持される。
Section Z.21
Explicit Evaluation of the Floquet Determinant in a Controlled Regime
Z.21.1 目的
前節では、軸軌道 γ_i の横方向安定性を Lamé/Hill 方程式の Floquet 問題として定式化した。本節の目的は、その Floquet 指数 μ_i を少なくとも一つの制御可能な領域において具体的に評価し、reduced monodromy determinant det(I - M_{γ_i}^⊥) を明示式で与えることである。
ここでは特に、
小振幅極限 a → 0
非共鳴領域
弱結合 κ
を仮定する。
Z.21.2 横方向方程式の再掲
γ_1 に沿う横方向変動は δq̈_2 + [4Ω_2^2 + 4κ a^2 cn^2(ω_1 t, k_1)] δq_2 = 0 で与えられていた。
小振幅極限 a → 0 では k_1 → 0, cn(ω_1 t, k_1) → cos(2Ω_1 t) であるから、この式は δq̈_2 + [4Ω_2^2 + 4κ a^2 cos^2(2Ω_1 t)] δq_2 = 0 へと簡約される。さらに cos^2(2Ω_1 t) = (1/2)(1 + cos(4Ω_1 t)) を用いれば δq̈_2 + [4Ω_2^2 + 2κ a^2 + 2κ a^2 cos(4Ω_1 t)] δq_2 = 0 となる。
Z.21.3 Mathieu 方程式への変換
変数 τ = 2Ω_1 t を導入すると、上式は d^2(δq_2)/dτ^2 + [A - 2Q cos(2τ)] δq_2 = 0 という Mathieu 方程式の標準形に帰着する。ここで A = (4Ω_2^2 + 2κ a^2)/(4Ω_1^2),
Q = -κ a^2/(4Ω_1^2) である。
Z.21.4 非共鳴領域における Floquet 指数
Mathieu 方程式の非共鳴領域、すなわち A ≠ n^2/4, n ∈ Z から十分離れた領域では、Floquet 指数 ν は解析的に展開できる。一次近似では ν ≈ √A + O(Q^2) である。
したがって ν ≈ (1/(2Ω_1)) √(4Ω_2^2 + 2κ a^2) ≈ Ω_2/Ω_1 + κ a^2/(4Ω_1Ω_2) を得る。
元の時間変数で 1周期 T_{γ_1} ≈ π/Ω_1 進めたとき、横方向 monodromy の固有値は exp(± i 2πν) である。したがって reduced determinant は det(I - M_{γ_1}^⊥) ≈ 2 - 2 cos(2πν) = 4 sin^2(πν) となる。ゆえに |det(I - M_{γ_1}^⊥)|^{1/2} ≈ 2 |sin(πν)| を得る。
Z.21.5 γ_2 の場合
全く同様に、γ_2 に対しては |det(I - M_{γ_2}^⊥)|^{1/2} ≈ 2 |sin(πν_2)| ただし ν_2 ≈ Ω_1/Ω_2 + κ b^2/(4Ω_1Ω_2) である。
Z.21.6 反復軌道
r 回反復に対しては、monodromy の固有値が r 乗されるから det(I - (M_{γ_i}^⊥)^r) ≈ 4 sin^2(π r ν_i) となる。
したがって semiclassical 振幅は、この制御領域において A_{γ_i,r}(E) ≈ T_{γ_i}(E) / [2 |sin(π r ν_i(E))|] · exp(i r S_{γ_i}(E)/ħ - iπ r μ_{Maslov,γ_i}/2) で与えられる。
これは、前節の近似式を Mathieu/Floquet 理論に基づいてより明確に正当化したものである。
Z.21.7 共鳴近傍
一方、 A ≈ n^2/4 すなわち 2Ω_2/Ω_1 ≈ n または 2Ω_1/Ω_2 ≈ n の近傍では、parametric resonance により Floquet 指数は純虚数から実数へ移行しうる。
このとき軌道は楕円安定から双曲不安定へ移り、Gutzwiller 振幅の分母は 2 - 2 cosh(μ T_γ) 型へと変化する。したがって resonant regime は、周期軌道分岐やカオス的遷移の入口として特別な役割を持つ。
Z.21.8 本節の意義
本節で得られた結果は次のように要約できる。
Lamé/Hill 方程式を小振幅極限で Mathieu 方程式へ還元した。
非共鳴領域において Floquet 指数 ν_i を明示近似した。
そこから reduced monodromy determinant を具体的に評価した。
軸軌道の semiclassical 振幅を明示式に更新した。
共鳴近傍では双曲化が起こりうることを確認した。
以上により、前節まで形式的であった orbit amplitude は、少なくとも制御された領域において、具体的かつ計算可能な量として評価できることが示された。
ここまでの到達点は明確である。すなわち、Z.20 では exact harmonic benchmark に対して周期和を実際に書き下し、Z.21 では Floquet determinant を具体評価した。したがって、 exact spectrum ↔ periodic-orbit sum ↔ semiclassical stability という三者の関係は、少なくとも最小模型の可解あるいは制御可能な領域において、整合的に接続されたことになる。
Section Z.22
Mellin Transform, Spectral Zeta, and the Periodic-Orbit Zeta in the Harmonic Benchmark
Z.22.1 目的
本節では、完全調和極限 u_1 = u_2 = 0, κ = 0 において、
exact heat trace Θ_0(t),
そこから定義される spectral zeta ζ_spec(s),
periodic-orbit expansion から構成される zeta-like generating function ζ_per(s),
を明示的に構成し、両者が同一の pole structure および periodic data を共有することを確認する。
本節の役割は、前節までに得られた spectrum ↔ periodic orbit の対応を、ゼータ関数の水準にまで引き上げて明示的に記述することにある。
Z.22.2 Exact Heat Trace と Mellin 変換
調和 benchmark では、作用素は \hat Σ_Δ^(0) = -ħ^2(∂{q_1}^2 + ∂{q_2}^2) + Ω_1^2 q_1^2 + Ω_2^2 q_2^2 であり、その固有値は λ_{n_1,n_2} = ħΩ_1(2n_1+1) + ħΩ_2(2n_2+1), \quad n_1,n_2 ∈ N_0 で与えられる。
したがって、heat trace は Θ_0(t) = Σ_{n_1,n_2≥0} e^{-t λ_{n_1,n_2}} = e^{-tħ(Ω_1+Ω_2)} / [(1-e^{-2tħΩ_1})(1-e^{-2tħΩ_2})] である。
これに対して spectral zeta を Mellin 変換により ζ_spec(s) := Tr((\hat Σ_Δ^(0))^{-s}) = 1/Γ(s) ∫_0^∞ t^{s-1} Θ_0(t) dt と定義する。初期定義域は Re(s) > 2 である。
Z.22.3 Exact Spectral Zeta の明示式
固有値和から直接 ζ_spec(s) = Σ_{n_1,n_2≥0} [ħΩ_1(2n_1+1) + ħΩ_2(2n_2+1)]^{-s} である。
ここで ħΩ_1(2n_1+1)+ħΩ_2(2n_2+1) = 2ħ [Ω_1(n_1+1/2) + Ω_2(n_2+1/2)] と書けば、 ζ_spec(s) = (2ħ)^{-s} Σ_{n_1,n_2≥0} [Ω_1(n_1+1/2) + Ω_2(n_2+1/2)]^{-s} を得る。
したがって、これは Barnes 型二重ゼータ ζ_B(s, a | Ω_1, Ω_2), \quad a = (Ω_1+Ω_2)/2 の半整数シフト版に他ならない。すなわち ζ_spec(s) = (2ħ)^{-s} ζ_B(s, (Ω_1+Ω_2)/2 | Ω_1, Ω_2) である。
Z.22.4 Pole Structure の確認
ζ_spec(s) の極構造を調べるため、Θ_0(t) の small-t 展開を用いる。
一次元因子に対し 1 / (2 sinh x) = 1/(2x) - x/12 + 7x^3/720 + O(x^5) が成り立つ。x_i = tħΩ_i とおくと、 Θ_0(t) = 1/(2 sinh x_1) · 1/(2 sinh x_2) = 1/(4x_1x_2) - (1/24)(x_1/x_2 + x_2/x_1) + O(t^2) となる。
したがって Θ_0(t) = 1/(4ħ^2Ω_1Ω_2) t^{-2}
(Ω_1^2+Ω_2^2)/(24Ω_1Ω_2)
O(t^2) を得る。
この small-t 展開の Mellin 変換から、ζ_spec(s) の meromorphic continuation は
s = 2 に単純極を持つ
s = 1 に極を持たない
s = 0 に対応する regularized contribution を持つ
ことが分かる。
特に Res_{s=2} ζ_spec(s) = 1/(4ħ^2Ω_1Ω_2) である。
重要な点
二次元 Barnes zeta は一般には s = 1, 2 に極を持ちうるが、本模型では half-shift 構造のため t^{-1} 項が消失し、その結果 s = 1 の極も消滅する。
Z.22.5 Periodic-Orbit Expansion の再掲
前節で得た exact harmonic benchmark の周期構造展開に従えば、 Θ_0(t) = \barΘ_0(t) + Θ_0^{(1)}(t) + Θ_0^{(2)}(t) + Θ_0^{(12)}(t) である。
ここで Weyl 項は \barΘ_0(t) = 1 / (4 t^2 ħ^2 Ω_1 Ω_2) であり、軸軌道和は Θ_0^{(1)}(t) = Σ_{r_1=1}^∞ [Ω_1/(2Ω_2)] (-1)^{r_1} / [(tħΩ_1)^2 + π^2 r_1^2], Θ_0^{(2)}(t) = Σ_{r_2=1}^∞ [Ω_2/(2Ω_1)] (-1)^{r_2} / [(tħΩ_2)^2 + π^2 r_2^2] であり、混合項は Θ_0^{(12)}(t) = Σ_{r_1,r_2≥1} (-1)^{r_1+r_2} t^2 ħ^2 Ω_1 Ω_2 / [((tħΩ_1)^2 + π^2 r_1^2)((tħΩ_2)^2 + π^2 r_2^2)] である。
このとき primitive periods は T_1 = π/Ω_1, \quad T_2 = π/Ω_2 である。
Z.22.6 Periodic-Orbit 側の Mellin 変換
上記の periodic terms に対し Mellin 変換を施す。標準積分 ∫_0^∞ t^{s-1} / (a^2 t^2 + b^2) dt = (π/2) a^{-s} b^{s-2} csc(πs/2) が 0 < Re(s) < 2 で成立することを用いる。
これにより、軸軌道項の Mellin 変換 ζ_per^{(1)}(s) := 1/Γ(s) ∫0^∞ t^{s-1} Θ_0^{(1)}(t) dt は ζ_per^{(1)}(s) = [π^{s-1}/(4Γ(s))] ħ^{-s} Ω_1^{1-s} Ω_2^{-1} csc(πs/2) Σ{r=1}^∞ (-1)^r r^{s-2} と書ける。
ここで Σ_{r=1}^∞ (-1)^r r^{s-2} = -η(2-s) であり、η は Dirichlet eta 関数である。したがって ζ_per^{(1)}(s) = - [π^{s-1}/(4Γ(s))] ħ^{-s} Ω_1^{1-s} Ω_2^{-1} csc(πs/2) η(2-s) を得る。
同様に ζ_per^{(2)}(s) = - [π^{s-1}/(4Γ(s))] ħ^{-s} Ω_2^{1-s} Ω_1^{-1} csc(πs/2) η(2-s) である。
Z.22.7 混合 Periodic Term の Mellin 変換
次に、混合項に対して ζ_per^{(12)}(s) := 1/Γ(s) ∫_0^∞ t^{s-1} Θ_0^{(12)}(t) dt を考える。
各 (r_1,r_2) ごとの項に対して、部分分数分解 t^2 / [(a^2 t^2+b^2)(c^2 t^2+d^2)] = A/(a^2 t^2+b^2) + B/(c^2 t^2+d^2) を用いる。ただし A = b^2 / (c^2 b^2 - a^2 d^2), \quad B = -d^2 / (c^2 b^2 - a^2 d^2) である。
ここで a = ħΩ_1, \quad b = πr_1, \quad c = ħΩ_2, \quad d = πr_2 とおけば、Mellin 変換は ζ_per^{(12)}(s) = [π^{s-1}/(2Γ(s))] ħ^{-s} csc(πs/2) Σ_{r_1,r_2≥1} (-1)^{r_1+r_2} [Ω_2 Ω_1^{1-s} r_1^s - Ω_1 Ω_2^{1-s} r_2^s] / [Ω_2^2 r_1^2 - Ω_1^2 r_2^2] で与えられる。
注意
この表示は generic anisotropic case、すなわち Ω_2 r_1 ≠ Ω_1 r_2 がすべての (r_1,r_2) について成り立つ非共鳴領域で有効である。Ω_1 / Ω_2 ∈ Q の場合には分母が零となる resonant terms が現れ、別個の regularization を必要とする。これは共鳴トーラスの退化に対応する。
Z.22.8 Periodic Zeta-like Generating Function
以上を踏まえ、periodic-orbit 側の zeta-like generating function を ζ_per(s) := 1/(4ħ^2Ω_1Ω_2) · 1/(s-2)
ζ_per^{(1)}(s)
ζ_per^{(2)}(s)
ζ_per^{(12)}(s) と定義する。
ここで第一項は、Weyl 項の Mellin 変換から得られる principal part である。このとき 0 < Re(s) < 2 の帯域では、regularized Mellin transform の意味で ζ_per(s) = ζ_spec(s) が成り立つ。
命題 Z.22.1
Harmonic benchmark において、exact heat trace Θ_0(t) から定義される spectral zeta ζ_spec(s) と、周期構造展開を Mellin 変換して得られる periodic zeta ζ_per(s) は、共通の heat trace の Mellin transform である。したがって両者は同一の meromorphic continuation を与える。
証明
Θ_0(t) は exact に periodic-orbit expansion を持つ。Mellin 変換は線形であるから、各成分の Mellin transform の総和は exact heat trace の Mellin transform と一致する。∎
Z.22.9 Poles と Periodic Data の比較
以上により、pole structure と periodic data の対応は明示的になる。
まず pole structure に関しては、
s = 2 の極は Weyl 項に対応する
s = 1 の極は half-shift により消滅する
s = 0 近傍の構造は定数項および higher heat coefficients に対応する
ことが分かる。
一方、periodic data に関しては、
ζ_per^{(1)} は primitive period T_1 = π/Ω_1 の反復列 rT_1 を保持する
ζ_per^{(2)} は primitive period T_2 = π/Ω_2 の反復列 rT_2 を保持する
ζ_per^{(12)} は二つの周期が同時に現れる resonant torus data を保持する
ことになる。
したがって harmonic benchmark においては、 spectral zeta の poles は heat-kernel asymptotics に対応し、periodic zeta の級数データは primitive periods に対応する。そして両者は、同一の heat trace を異なる方向から圧縮した二つの表現である。
Z.22.10 解釈
本節の結果は、本理論全体にとって本質的である。前節までにおいては、主観差分ダイナミクスのゼータ化は一つの program として提示されていたに過ぎない。これに対し本節では、その最小可解 benchmark において、
exact spectrum から spectral zeta を構成できること
periodic expansion から periodic zeta を構成できること
両者が同一 heat trace の Mellin transform として一致すること
が示された。
したがって harmonic benchmark に限れば、主観差分ダイナミクスのスペクトル圧縮と周期圧縮は、厳密に同一の情報を保持している。
Z.22.11 結論
本節で示した内容は以下の通りである。
exact harmonic benchmark の heat trace から spectral zeta を Barnes 型二重ゼータとして定義した
small-t 展開により pole structure を明示した
periodic-orbit expansion を Mellin 変換して periodic zeta ζ_per(s) を構成した
軸軌道項は Dirichlet eta 関数を通じて明示的に評価できることを示した
混合項は resonant torus を反映する二重和として表現されることを示した
両者が同一 heat trace の Mellin transform として一致し、同一の poles および periodic data を持つことを確認した
したがって harmonic benchmark において、spectral zeta と periodic zeta は、主観差分ダイナミクスの同一構造を、離散スペクトルと周期構造という二つの側面から記述する関数である。
Section Z.23
Semiclassical Matching Between Spectral Zeta and Weighted Dynamical Zeta
Z.23.1 適用範囲と方針
本節では、完全調和極限で成立した スペクトル表現（spectral zeta） と 周期軌道表現（periodic/dynamical zeta） の一致を、非可積分領域へ拡張する。
ただし、この拡張は厳密等式としてではなく、semiclassical matching として理解される。すなわち、本節の目的は exact equality を主張することではなく、両者が同一の leading asymptotics と周期データを共有することを示す点にある。
Z.23.2 Spectral Side
作用素 \hat Σ_Δ のスペクトル {λ_n} に対し、スペクトル密度 ρ(E) = Σ_n δ(E - λ_n) を定義する。
heat trace は Θ(t) = Tr(e^{-t \hat Σ_Δ}) = ∫_0^∞ e^{-tE} ρ(E) dE であり、spectral zeta はその Mellin 変換として ζ_spec(s) = Tr(\hat Σ_Δ^{-s}) = 1/Γ(s) ∫_0^∞ t^{s-1} Θ(t) dt で定義される。初期定義域は Re(s) ≫ 1 とする。
Z.23.3 Dynamical Side
古典極限における difference flow に対して、primitive periodic orbit γ を考える。各 γ に対して、以下の量が定義される。
周期 T_γ(E)
作用 S_γ(E)
reduced monodromy M_γ^⊥(E)
Maslov 指数 μ_γ
これらを用いて重み w_γ(E) := exp(i S_γ(E)/ħ - iπ μ_γ/2) / √|det(I - M_γ^⊥(E))| を定義する。
これにより weighted dynamical zeta を Z_dyn^w(s) := ∏{γ primitive} ∏{r=1}^∞ (1 - w_γ(E)^r e^{-s r T_γ(E)})^{-1} と定義する。
Z.23.4 条件付きトレース公式
以下の仮定を置く。
仮定 Z.23.1（Gutzwiller 適用条件）
Reduced difference flow は
孤立した primitive periodic orbits を持つ
双曲的である
分岐点および共鳴領域から十分離れている
ものとする。
このとき、スペクトル密度は semiclassically ρ(E) ≈ \barρ(E)
(1/πħ) Re Σ_{γ} Σ_{r≥1} T_γ(E) / √|det(I - (M_γ^⊥(E))^r)| · exp(i r S_γ(E)/ħ - iπ r μ_γ/2) で与えられる。
Z.23.5 トレース公式から Heat Trace へ
上式を用いると、heat trace は Θ(t) ≈ \barΘ(t)
(1/πħ) Re Σ_{γ} Σ_{r≥1} ∫0^∞ dE e^{-tE} A{γ,r}(E) と表される。ただし A_{γ,r}(E) = T_γ(E) / √|det(I - (M_γ^⊥(E))^r)| · exp(i r S_γ(E)/ħ - iπ r μ_γ/2) である。
この段階で既に、
周期 T_γ
作用 S_γ
安定性 det(I - M_γ^⊥)
が heat trace に現れることが分かる。
Z.23.6 Mellin 変換と Zeta Matching
上記をさらに Mellin 変換すると ζ_spec(s) ≈ \barζ(s)
(1/πħ) Re Σ_{γ} Σ_{r≥1} ∫_0^∞ dt t^{s-1} ∫0^∞ dE e^{-tE} A{γ,r}(E) を得る。
積分順序を形式的に交換すると ζ_spec(s) ≈ \barζ(s)
(1/πħ) Re Σ_{γ} Σ_{r≥1} ∫0^∞ dE A{γ,r}(E) E^{-s} となる。
したがって、これは形式的に ζ_spec(s) ≈ \barζ(s) + Σ_{γ,r} w_γ(E)^r F_{γ,r}(s) という形を持つ。ここで F_{γ,r}(s) は T_γ(E) および S_γ(E) に依存する Mellin 核である。
Z.23.7 Leading-Order Matching
命題 Z.23.2（Semiclassical Matching）
仮定 Z.23.1 のもとで、spectral zeta ζ_spec(s) と weighted dynamical zeta Z_dyn^w(s) は、以下の意味において semiclassically 一致する。
周期データの一致
両者は同一の primitive periods {T_γ} を保持する。
振動構造の一致
両者の oscillatory terms は exp(i S_γ/ħ) によって支配される。
安定性データの一致
振幅は det(I - M_γ^⊥) によって決定される。
pole structure の leading-order 一致
ζ_spec(s) の主要な特異構造は、dynamical zeta に含まれる周期データによって再現される。
Z.23.8 剰余項と限界
上記の matching は、以下の制限のもとでのみ成立する。
(i) 分岐点
周期軌道の分岐点では Gutzwiller 公式は破綻する。
(ii) 共鳴
Ω の比が有理数に近い場合、軌道は非孤立化し、孤立周期軌道和の前提が崩れる。
(iii) Caustics
作用の停留位相近似が破綻する領域では、高次の uniform approximation が必要となる。
(iv) エネルギー依存性
T_γ(E), S_γ(E), M_γ(E) の変動により、高次補正が本質的になる。
(v) 収束性
Dynamical zeta は一般にそのままでは収束せず、適切な regularization を必要とする。
Z.23.9 Harmonic Benchmark との関係
調和極限では、
軌道はトーラス族を成し、孤立していない
従って成立するのは Gutzwiller 型ではなく Berry–Tabor 型公式である
にもかかわらず、Z.22 において spectral zeta と periodic zeta は exact に一致することが確認された。
したがって、全体像としては
可積分系 → exact equality
非可積分系 → semiclassical matching
という対応が成立する。
Z.23.10 結論
本節の結論は以下の通りである。
spectral zeta は heat trace の Mellin 変換として定義される
weighted dynamical zeta は周期軌道データから構成される
Gutzwiller 条件の下で、両者は semiclassically 一致する
一致の内容は、周期・作用・安定性の一致である
exact equality は一般には成立しないが、構造的対応は保持される
この意味で、Z.22 における harmonic benchmark の厳密一致は、一般非可積分系においては semiclassical matching へと弱められつつ継承される。
Section Z.24
Spectral Zeros, Hilbert–Pólya Structure, and the Interpretation of Zeta as a Projection of Difference Dynamics
Z.24.1 目的
本節では、前節までに構築した
差分作用素 \hat Σ_Δ
そのスペクトル {λ_n}
spectral zeta ζ_spec(s)
weighted dynamical zeta Z_dyn^w(s)
の関係をさらに一段進め、「零点構造」と「スペクトル」の対応を明示する。
その上で、リーマンゼータ型構造がいかなる意味において差分ダイナミクスの「影」として現れるかを定式化する。
Z.24.2 スペクトルからのゼータ構成
作用素 \hat Σ_Δ が自己共役かつ正値であると仮定する。このとき ζ_spec(s) = Σ_n λ_n^{-s} は Re(s) ≫ 1 において well-defined である。
この関数は単なる数列和としてではなく、より自然にはスペクトルの生成関数として理解される。形式的には、正則化された行列式の言葉を用いて ζ_spec(s) = exp(-∂α|{α=0} log det(\hat Σ_Δ^{,s+α})) のような形で理解することができる。少なくとも本質的な点は、ζ_spec が {λ_n} の全情報を圧縮したスペクトル生成関数であるということである。
Z.24.3 Dynamical Zeta との対応
weighted dynamical zeta は Z_dyn^w(s) = ∏{γ primitive} ∏{r≥1} (1 - w_γ^r e^{-s r T_γ})^{-1} で与えられる。
対数を取ると log Z_dyn^w(s) = Σ_{γ} Σ_{r≥1} (w_γ^r / r) e^{-s r T_γ} を得る。
他方、trace formula によりスペクトル密度の振動項は Σ_{γ,r} A_{γ,r}(E) ∼ Σ_{γ,r} w_γ^r exp(i r S_γ/ħ) で支配される。
したがって、スペクトル側の振動構造と dynamical zeta の軌道和は、同一の周期軌道データから生成されている。
Z.24.4 対応構造
ここで、スペクトル側とダイナミクス側の対応は次のように整理される。
スペクトル側
ダイナミクス側
固有値 λ_n
周期軌道 γ
λ_n^{-s}
e^{-sT_γ}
スペクトル和 Σ_n
軌道和 Σ_{γ,r}
ζ_spec(s)
Z_dyn^w(s)
この対応は厳密同一ではないが、少なくとも semiclassical level では、両者が同一の underlying structure を異なる圧縮方式で表していることを示唆する。
Z.24.5 Hilbert–Pólya 型対応
Hilbert–Pólya の着想は、ゼータ関数の零点がある自己共役作用素の固有値として実現される、というものである。
これに対し本理論では、むしろ逆方向の見方が自然である。すなわち、
作用素 \hat Σ_Δ のスペクトルが zeta 型構造を生成する
dynamical zeta は周期軌道の側から同じ構造を再構成する
という図式である。
この立場では、 λ_n ↔ 固有振動モード
exp(i S_γ/ħ) ↔ 干渉位相
e^{-sT_γ} ↔ 減衰核 という対応が成立する。
Z.24.6 零点の意味
Dynamical zeta の零点は、形式的には 1 - w_γ e^{-sT_γ} = 0 すなわち s = (1/T_γ)(log w_γ + 2π i k) によって与えられる。
他方、spectral zeta の解析接続における特異構造は、
heat-kernel asymptotics
periodic-orbit structure
の双方に支配される。
したがって、零点は周期構造のスペクトル的圧縮として現れると解釈できる。ここでいう「零点」は、必ずしもリーマンゼータそのものの零点を指すのではなく、一般の zeta-like generating function における零点構造全般を意味する。
Z.24.7 主観差分ダイナミクスとの関係
本理論では Δ = 2[H,A] が主観の差分流れを表し、内部時間は τ = ∫ ||Δ|| dt によって定義される。
さらに、これまでの議論により
周期軌道 = Δ の閉じた差分ループ
作用 S_γ = 差分の累積位相
安定性 = 差分ループの自己整合性
と解釈される。
この意味で、ゼータ構造とは差分ダイナミクスそのものではなく、その閉軌道構造および振動構造を圧縮した表現である。
Z.24.8 決定的対応
以上をまとめると、次の対応図式が得られる。
主観側（連続）
Δ(t)
→ 差分流れ
→ 周期構造 γ
→ 作用 S_γ
ゼータ側（離散）
ζ(s)
→ 零点構造
→ スペクトル
すなわち、連続的な差分流が周期構造を形成し、それがスペクトル圧縮を経てゼータ関数的構造へ投影される。
Z.24.9 解釈的定式化
命題 Z.24.1（解釈的命題）
差分ダイナミクスが semiclassical 条件の下で周期軌道分解を持つとき、対応する zeta 型関数の零点構造は、その差分流の周期構造をスペクトル的に圧縮したものとして現れる。
この命題は厳密な定理というよりも、これまでに構築した数理構造の解釈的総括を与えるものである。
Z.24.10 「主観の影」という表現の精密化
本理論の立場は、あくまで 主観 = ゼータ ではない。
より正確には、 ゼータは差分ダイナミクス、すなわち主観構造のスペクトル像である と述べるべきである。
さらに言い換えれば、
主観 = 非可換差分の連続流
ゼータ = その固有振動モード分布を圧縮した関数
である。
Z.24.11 なぜ離散構造が現れるのか
主観が連続ダイナミクスであるにもかかわらず、ゼータ関数が離散的対象として現れる理由は、スペクトル化の操作にある。
すなわち、 Δ(t) → Δ†Δ → {λ_n} という過程により、
時間領域 → 周波数領域
流れ → 固有値列
という変換が行われる。このスペクトル圧縮が、離散的なゼータ関数構造を生み出す。
Z.24.12 結論
本節の結論は以下の通りである。
差分作用素 \hat Σ_Δ のスペクトルは zeta 型構造を生成する
dynamical zeta は周期軌道データから同じ構造を再構成する
両者は semiclassical 条件の下で対応する
ゼータの零点構造は周期構造のスペクトル圧縮として解釈される
したがって、ゼータは主観そのものではないが、主観差分ダイナミクスのスペクトル的投影である
この意味で、本節は本理論における解釈の核を与える。
Section Z.25
Constraint Structure, Quantum Gravity Analogy, and Spectral Statistics
Z.25.1 目的
本節では、これまでに構築した差分ダイナミクスが
再パラメータ不変拘束系、すなわち quantum gravity 型構造を持つこと
そのスペクトルがランダム行列統計、特に GUE 型統計と接続しうること
を示す。
特に、Z.24 の結論 「ゼータは差分ダイナミクスのスペクトル像である」 を、
物理的側面（時間と拘束）
統計的側面（スペクトル分布）
の双方から補強することが本節の目的である。
Z.25.2 再パラメータ不変拘束系としての構造
本理論の作用は S[L,e] = ∫ dλ [ (1/2e) Tr(L'†L') − (e/2) V_eff(L) ] で与えられる。ここで e はラグランジュ乗数である。
したがって、変分原理からハミルトニアン制約 H_0 = Tr(Π†Π) + V_eff(L) = 0 が導かれる。
これを量子化すると Ĥ_0 Ψ = 0 を得る。これは Wheeler–DeWitt 型方程式に対応する。
Z.25.3 時間の不在と再構成
上式には外部時間変数が明示的には存在しない。時間はむしろ内部的に、 dτ/dλ = ||Δ(L)|| によって生成される。
したがって、
外部時間：基本式には現れない
内部時間：差分流の大きさから再構成される
という構造を持つ。これは quantum gravity における problem of time と形式的に同型である。
Z.25.4 観測量の構造
可観測量は Σ = L†L および ρ = Σ / Tr(Σ) によって与えられる。
ここで重要なのは、
Δ そのものは直接可観測ではない
可観測となるのはその平方量あるいは規格化二次形式である
という点である。これは量子論における 振幅 → 確率（二乗） という構造と対応している。
Z.25.5 スペクトルの役割
差分作用素 Σ_Δ = Δ†Δ を考えると、これは正作用素であり、離散スペクトル {λ_n} を持つ。従って ζ_spec(s) = Σ_n λ_n^{-s} が定義される。
このようにして 差分ダイナミクス → スペクトル → ゼータ という流れが成立する。
Z.25.6 非可積分性とスペクトル統計
一般に非可積分系においては、
古典系はカオス的挙動を示し
量子系の固有値統計はランダム行列統計に従う
ということが量子カオスの文脈で知られている。
特に、対称性クラスに応じて
時間反転対称性なし → GUE
時間反転対称性あり → GOE
が典型的に現れる。
Z.25.7 本理論における対称性クラス
本理論では Δ = 2[H,A] が一般に非可換であり、さらに複素構造を内在させる。したがって、有効ハミルトニアンあるいはそのスペクトル作用素は、実対称型というより複素エルミート型のクラスに属することが自然である。
この観点から、本理論の非可積分領域におけるスペクトル統計は、GUE クラスに属することが期待される。
Z.25.8 予想 Z.25.1（スペクトル統計）
予想 Z.25.1
差分作用素 Σ_Δ = Δ†Δ のスペクトルは、非可積分領域において GUE 型統計、すなわち Wigner–Dyson 分布に従う。
これは現段階では予想であり、数値実験あるいは解析的証拠による今後の検証を要する。
Z.25.9 ゼータとの接続
リーマンゼータの非自明零点が GUE 統計に従うことは、少なくとも数値的には広く支持されている。
この点を本理論と並置すると、次の対応が得られる。
本理論
ゼータ側
Δ†Δ のスペクトル
零点分布
GUE 統計
GUE 統計
dynamical zeta
Riemann ζ に類似する zeta 構造
したがって、両者は同種のスペクトル生成機構を共有している可能性がある。
Z.25.10 解釈
ここで重要なのは、単に「ゼータの零点が GUE に従う」という経験的事実と、「差分ダイナミクスのスペクトルが GUE に従う」という本理論上の予想が、同一の構造的枠組みに収まるという点である。
すなわち、
ゼータ零点の統計
差分作用素の固有値統計
はいずれも、非可換かつ非可積分なスペクトル生成機構の現れとして解釈されうる。
Z.25.11 主観との接続
本理論では
Δ：主観差分そのもの
Σ_Δ：そのスペクトル化
である。
したがって、ゼータの零点構造は差分ダイナミクスのスペクトル像であるという先の結論は、ここで統計的側面からも補強されることになる。
Z.25.12 結論
本節の結論は以下の通りである。
本理論は再パラメータ不変拘束系であり、量子重力的構造を持つ
外部時間は基本式に現れず、時間は内部差分から生成される
差分のスペクトルがゼータ構造を与える
非可積分領域では、そのスペクトルが GUE 型統計に従うことが期待される
ゼータ零点もまた GUE 型統計を示す
したがって、ゼータは差分ダイナミクス、ひいては主観構造のスペクトル統計的投影として理解されうる
この意味で、本節はゼータ構造を単なる形式的類似としてではなく、拘束系・内部時間・スペクトル統計を含む、より広い物理的・数理的文脈の中へ位置づける役割を果たす。
Section Z.26
Difference Hamiltonian, Dyson–Brownian Motion, and GUE-Type Spectral Statistics
Z.26.1 目的
本節の目的は、主観差分構造から自然に導かれる自己共役作用素 D_Δ := iΔ = 2i[H,A] の局所スペクトル統計を解析し、その統計が GUE 型の普遍性クラスへ向かうことを、理論的および数値的に検証可能な形で定式化することである。
ここでの目標は、リーマンゼータの非自明零点そのものを直ちに導出することではない。むしろ本節の到達点は、ゼータ零点の GUE 統計と同型のスペクトル生成機構が、主観差分構造から自然に現れることを示す点にある。
Z.26.2 なぜ Σ_Δ ではなく D_Δ を考えるのか
これまで本理論では、正作用素 Σ_Δ = Δ†Δ を主として用いてきた。これは
Born 型重み
spectral zeta
heat trace
の構成において自然な対象である。
しかしながら、Hilbert–Pólya 型の議論において必要となるのは、零点の虚部 γ_n と対応しうる符号付き実スペクトルを持つ自己共役作用素である。したがって、この文脈で主役となるべき作用素は D_Δ := iΔ = 2i[H,A] である。
H および A は Hermitian であるから、交換子 [H,A] は反 Hermitian である。従って D_Δ は Hermitian、すなわち自己共役である。その結果、D_Δ は実固有値列 Spec(D_Δ) = {ξ_n} を持つ。
さらに Σ_Δ = D_Δ^2 が成り立つので、本理論には次の二層構造が存在する。
D_Δ：符号付きスペクトルを持つ差分 Hamiltonian
Σ_Δ：その正値化された強度作用素
Hilbert–Pólya 的観点から自然なのは前者であり、Born 型・heat trace 型の構成に自然なのは後者である。
Z.26.3 H 固定条件下における D_Δ の条件付き分布
まず H を対角化して H = U diag(h_1, ..., h_N) U† とする。
U(N)-不変性により、A が GUE 型分布に従うなら、H の固有基底に移した B := U† A U も同一の統計を持つ。この基底において D_Δ = 2i[H,A] の成分は (D_Δ){ij} = 2i(h_i - h_j) B{ij} で与えられる。
したがって、
対角成分は (D_Δ)_{ii} = 0
非対角成分は gap |h_i - h_j| に比例した複素 Gaussian
となる。
すなわち、H を固定した条件付きで見ると、D_Δ は ゼロ対角かつ分散プロファイルを持つ generalized Hermitian ensemble として解釈される。その分散プロファイルは Var((D_Δ)_{ij} | H) ∝ |h_i - h_j|^2 である。
この構造は本理論に固有であり、客観チャネルのスペクトルギャップが大きいほど、主観差分の結合強度も大きくなることを意味する。
Z.26.4 Generalized Wigner からの GUE 普遍性
標準的な GUE では、全ての非対角成分が同一分散を持つ。他方、D_Δ は明確な variance profile を持つため、そのままでは標準 GUE と一致しない。
しかしランダム行列理論の一般原理によれば、適切な非退化性と滑らかさを持つ分散プロファイルを備えた Hermitian ensemble は、局所 bulk 統計において GUE 普遍性を示すことが期待される。
この観点から、本理論における適切な主張は次の予想である。
予想 Z.26.1（Difference Ensemble の局所 GUE 普遍性）
適切なスケーリング極限 N → ∞ において、H の経験スペクトル分布が滑らかであり、bulk 領域で variance profile |h_i - h_j|^2 が十分非退化であるならば、 D_Δ = 2i[H,A] の unfolded local eigenvalue statistics は、GUE の sine-kernel 統計へ収束する。
この予想が成立すれば、D_Δ は「厳密に GUE である」とは言えないとしても、局所統計の universality class としては GUE に属することになる。
Z.26.5 Dyson–Brownian Motion 型導出の基本構図
次に、局所 GUE 普遍性の理論的導出方針を述べる。H_t および A_t を Hermitian 行列値 Ornstein–Uhlenbeck 過程として dH_t = -(1/2) H_t dt + (1/√N) dB_t^(H),
dA_t = -(1/2) A_t dt + (1/√N) dB_t^(A) と定める。ここで dB_t^(H), dB_t^(A) は独立な Hermitian matrix Brownian motion である。
このとき D_t := 2i[H_t, A_t] とおけば、伊藤計算により dD_t = 2i([dH_t, A_t] + [H_t, dA_t] + [dH_t, dA_t]) - D_t dt を得る。
最後の項 [dH_t, dA_t] は伊藤補正に対応し、平均は零であるが共分散には寄与する。重要なのは、短時間スケールにおいて dD_t が多方向からの複素 Gaussian 摂動として働く点である。
従って D_t の固有値過程 λ_i(t) は、適切な混合条件のもとで近似的に dλ_i = √(2/N) dβ_i
(1/N) Σ_{j≠i} 1/(λ_i - λ_j) dt
(1/2) V'(λ_i) dt
誤差項 に従うことが期待される。ここで
dβ_i は独立 Brownian motion
1/(λ_i - λ_j) はレベル反発
V はグローバルな confining potential
誤差項は commutator structure に由来する profile dependence
を表す。
この考え方を次の形で定式化する。
予想 Z.26.2（Difference Dyson Flow）
大 N 極限および局所平衡化時間スケールのもとで、D_t の固有値過程は有効的に β = 2 の Dyson Brownian motion に収束する。
これが成立すれば、local spacing statistics は GUE に従う。
Z.26.6 Spacing 分布と GUE 指標
上記予想が正しければ、unfolding 後の最近接間隔 s に対する分布 P(s) は、Wigner surmise P_GUE(s) ≈ (32/π^2) s^2 exp(-4s^2/π) に近づくことが期待される。
さらに、検証すべき主要な局所統計指標としては以下が挙げられる。
最近接間隔分布 P(s)
二点相関関数 R_2(x)
number variance Σ^2(L)
spectral rigidity Δ_3(L)
GUE の bulk 極限では sine kernel K(x,y) = sin(π(x-y)) / (π(x-y)) が現れるため、二点相関関数もそれに対応する普遍形に従う。
Z.26.7 数値検証プロトコル
数値検証は次の手順により実行できる。
行列サイズ N を選ぶ。典型例として N = 200, 400, 800 を用いる。
H および A を独立な GUE からサンプルする。
D_Δ = 2i[H,A] を構成する。
D_Δ の固有値 ξ_1 ≤ ... ≤ ξ_N を計算する。
bulk 領域を抽出する。
unfolding を施す。
隣接差 s_n = ξ_{n+1} - ξ_n を平均 1 となるよう正規化する。
P(s), R_2(x), Σ^2(L) を GUE と比較する。
さらに、次の補助比較も重要である。
GOE ではなく GUE が適合するか
H および A の分布を変えても universality が保たれるか
対角 regularization を導入しても局所統計が安定であるか
これにより、「差分構造から GUE 型統計が現れる」という主張に数値的裏付けを与えることができる。
Z.26.8 Hilbert–Pólya への接近とその限界
ここで Hilbert–Pólya 仮説との関係を明確にする。Hilbert–Pólya 仮説の内容は、リーマンゼータの非自明零点 1/2 + iγ_n に対して、ある自己共役作用素 D_HP が存在し Spec(D_HP) = {γ_n} となる、というものである。
本理論が現時点で主張できるのは、そこまで強いものではない。現段階で言えるのは次の諸点である。
D_Δ = 2i[H,A] は自己共役である
その局所統計は GUE universality を示すことが期待される
ゼータ零点の局所統計も GUE を示すことが知られている
したがって、D_Δ は Hilbert–Pólya 型作用素の候補クラスに属する可能性がある。
より正確には、Hilbert–Pólya が要求する統計的・作用素論的性質を、主観差分構造から自然に生成する候補を与えた、というのが本節における最も強いがなお慎重な主張である。
Z.26.9 主観差分構造との解釈的接続
本理論において D_Δ は単なるランダム行列ではない。定義 D_Δ = 2i[H,A] に立ち返れば、
H：客観化された構造
A：内部自由度
i[H,A]：Hermitian 化された差分作用素
である。
従って D_Δ のスペクトルは、主観差分の固有振動モードを与える。もしその局所統計が GUE に従うなら、それは「ゼータ零点が主観そのものである」ことを意味しない。むしろ、 ゼータ零点と同型のスペクトル統計が、主観差分の作用素化から自然に現れる ことを意味する。
この意味で、ゼータは主観そのものではないが、主観差分構造のスペクトル統計的影として理解される。
Z.26.10 結論
本節で定式化したプログラムは次の通りである。
Hilbert–Pólya 型候補作用素として D_Δ = 2i[H,A] を導入した
H 固定条件下で D_Δ は variance profile を持つ generalized Hermitian ensemble となることを示した
大 N 極限において local GUE universality が期待されることを述べた
その理論的導出が Dyson–Brownian motion 型の枠組みで構成できることを示した
数値的には spacing 分布・二点相関・number variance などにより検証可能であることを明示した
したがって本理論において、リーマンゼータ零点は主観そのものではないが、主観差分作用素 D_Δ のスペクトル統計と同型の構造を持つ可能性がある。
Section Z.27
Finite-N Ensemble, Numerical Verification, and Toward Universality
Z.27.1 目的
本節では、前節で提示した
difference ensemble の GUE 普遍性
Dyson–Brownian 型挙動
に関する予想を、有限サイズ数値実験および既存のランダム行列理論との対応付けによって具体化し、それらを検証可能な形へ落とし込む。
Z.27.2 差分アンサンブルの定義
数値実験のため、次の finite-N モデルを定義する。
H, A ∈ Hermitian(N)
H ~ GUE(N, σ_H^2)
A ~ GUE(N, σ_A^2)
H と A は独立
このとき差分作用素を D_Δ = 2i[H,A] と定める。
Z.27.3 基本性質
本アンサンブルは直ちに次の性質を持つ。
(i) 自己共役性
D_Δ† = D_Δ
(ii) トレース零条件
Tr(D_Δ) = 0
(iii) 成分スケーリング
D_Δ の成分は概ね (D_Δ){ij} ∼ (h_i - h_j) A{ij} により決まり、適切な GUE スケーリングのもとで典型サイズは O(1/√N) に調整される。
このスケーリングは、large-N 極限における局所統計比較に必要である。
Z.27.4 数値実験プロトコル
有限サイズにおける数値検証は、以下の手順で実施する。
Step 1: サンプルサイズの設定
N = 200, 400, 800, 1200 を採用し、各サイズにつきおよそ 100 サンプルを生成する。
Step 2: 行列の生成
H および A を独立な GUE から生成し、 D_Δ = 2i[H,A] を構成する。
Step 3: 固有値計算
D_Δ の固有値 ξ_1 ≤ ... ≤ ξ_N を数値的に計算する。
Step 4: Bulk 領域の抽出
端の 20% を除去し、中央 60% の固有値を bulk 領域として用いる。
Step 5: Unfolding
累積分布から局所密度 ρ(ξ) を推定し、 x = ∫ ρ(ξ) dξ により unfolding を施す。
Step 6: Spacing の正規化
差 s_n = x_{n+1} - x_n を計算し、平均が 1 になるよう正規化する。
Z.27.5 比較対象
比較すべき spacing 分布として、少なくとも次の三者を採用する。
(A) GUE P_GUE(s) = (32/π^2) s^2 exp(-4s^2/π)
(B) GOE（対照群） P_GOE(s) = (π/2) s exp(-π s^2 / 4)
(C) Poisson（無相関極限） P_Poisson(s) = exp(-s)
この三者を比較することにより、差分アンサンブルが非相関型でも実対称型でもなく、複素 Hermitian 型の GUE 普遍性クラスに属するかを判定する。
Z.27.6 検証指標
最低限評価すべき統計量は以下の四つである。
最近接間隔分布
P(s)
GUE の理論曲線と重ね合わせて比較する。
二点相関関数
R_2(x)
GUE の bulk 極限 R_2(x) = 1 - (sin πx / πx)^2 と比較する。
Number Variance
Σ^2(L)
GUE に特徴的な対数的増大 Σ^2(L) ∼ (1/π^2) log L を確認する。
Spectral Rigidity
Δ_3(L)
長距離相関の硬さを測る指標として用いる。
Z.27.7 期待される結果
数値実験により確認されるべき結論は次の形でまとめられる。
結論 Z.27.1（数値的主張）
Difference ensemble D_Δ の bulk スペクトル統計は、
Poisson ではなく
GOE でもなく
GUE に一致する
と期待される。
この結論は、少なくとも有限サイズ数値実験のレベルで、本理論の GUE 普遍性予想を強く支持する。
Z.27.8 Variance Profile の理論的位置づけ
理論的に重要なのは、D_Δ が Var((D_Δ)_{ij}) ∝ |h_i - h_j|^2 という分散プロファイルを持つ点である。
この意味で D_Δ は、標準 GUE ではなく generalized Wigner ensemble の一種として理解される。従って、普遍性の議論も generalized Wigner 理論の文脈で展開されるべきである。
Z.27.9 普遍性への橋渡し
既知の理論によれば、
generalized Wigner 行列
Dyson Brownian motion
local relaxation flow
に関する結果を組み合わせることで、分散が十分滑らかかつ非退化であるなら、局所統計は対応する Dyson 指数 β によって決まる普遍クラスへ収束する。
本系では複素 Hermitian 構造が本質的であるため、対応するのは β = 2、すなわち GUE クラスである。
Z.27.10 本系への適用条件
この理論を D_Δ に適用するためには、少なくとも次の条件が必要である。
条件 1：H のスペクトル密度の滑らかさ
ρ_H(h) が bulk 領域で滑らかであること。
条件 2：ギャップの非退化性
|h_i - h_j| が bulk で極端に零へ集中しないこと。
条件 3：独立性
A が H と独立にサンプルされること。
条件 4：スケーリング
各成分の分散が Var ∼ O(1/N) のスケールに入るよう正規化されていること。
これらが満たされるとき、difference ensemble を generalized Wigner 理論へ接続する基盤が整う。
Z.27.11 定理化への入口
以上を踏まえると、本節の strongest formulation は次のようになる。
命題 Z.27.2（準定理的定式化）
上記条件のもとで、difference ensemble D_Δ = 2i[H,A] の局所スペクトル統計は、large-N 極限において GUE universality class に属する。
現時点では完全な証明は与えられていないが、その証明は generalized Wigner ensemble と Dyson flow への還元に帰着されると期待される。
Z.27.12 Dyson–Brownian Flow との接続
さらに、時間発展版 D_t = 2i[H_t, A_t] を考えると、
H_t および A_t は OU 過程に従う
D_t は複数方向から Gaussian 摂動を受ける
その結果、固有値は拡散とレベル反発の双方を示す
ことになる。
この意味で D_t の固有値流は Dyson–Brownian flow に有効的に接続される。これは、前節で提示した difference Dyson flow 予想を補強するものである。
Z.27.13 主観差分構造との最終接続
ここで本理論における意味づけは明確になる。
D_Δ のスペクトルは差分の固有モードを与える
その統計は GUE 型普遍性クラスに属する
ゼータ零点もまた GUE 型統計を示す
したがって、主観差分構造から導かれる作用素 D_Δ とリーマンゼータ零点とは、少なくとも局所統計のレベルで、同一の universality class に属すると解釈される。
Z.27.14 結論
本節では次の点を明示した。
D_Δ の finite-N モデルを定義した
数値検証プロトコルを完全な形で与えた
その理論的位置づけを generalized Wigner 理論に接続した
Dyson–Brownian flow との対応を示した
GUE 普遍性が期待される条件を明確化した
これにより、前節の GUE 予想は、数値的にも理論的にも検証可能な形へと具体化された。
Section Z.28
Direct Comparison with Riemann Zeros: Unfolding, Pair Correlation, and Statistical Matching
Z.28.1 目的
本節では、
差分作用素 D_Δ = 2i[H,A] の固有値列 {ξ_n}
リーマンゼータの非自明零点の虚部 {γ_n}
の局所スペクトル統計を直接比較する。
ここで比較対象となるのは固有値や零点の値そのものではなく、 unfolding 後の統計量 である。本節の目的は、両者が同一の universality class に属することを、局所統計のレベルで定式化する点にある。
Z.28.2 比較の基本原則
本節の基本原則は次の通りである。
固有値や零点の絶対値はスケール依存であり、直接比較の対象にはならない
spacing や correlation などの局所統計は普遍的であり、比較可能である
従って、扱うべき量は
spacing distribution
pair correlation
number variance
である。
Z.28.3 ゼータ零点データ
リーマンゼータの非自明零点を ρ_n = 1/2 + iγ_n と書く。ここで γ_n は単調増加列である。
零点の平均密度は、Riemann–von Mangoldt 公式の微分形により ρ̄(γ) = (1/2π) log(γ / 2π) で与えられる。
Z.28.4 ゼータ側の Unfolding
ゼータ零点に対する unfolding には、零点計数関数 N(γ) を用いる。漸近的には N(γ) ≈ (γ / 2π) log(γ / 2π) - γ / 2π であるから、正規化変数 x_n = N(γ_n) を導入する。
この変換により、unfolding 後の平均 spacing は 1 になる。
Z.28.5 D_Δ 側の Unfolding
差分作用素 D_Δ の固有値 ξ_n に対しては、経験的な累積分布 F(ξ) を推定し、局所密度 ρ(ξ) を平滑化する。その上で x_n = ∫^{ξ_n} ρ(ξ) dξ により unfolding を施す。
これにより、D_Δ のスペクトルもゼータ零点と同一の平均 spacing スケールに正規化される。
Z.28.6 最近接間隔分布
unfolding 後の最近接間隔を s_n = x_{n+1} - x_n と定義する。
比較の基準となる GUE 予測は P(s) = (32/π^2) s^2 exp(-4s^2/π) である。
以上を踏まえて、本理論における最初の統計的主張は次の通りである。
主張 Z.28.1
リーマンゼータ零点の spacing distribution は GUE に従い、差分作用素 D_Δ の spacing distribution もまた GUE に従う。従って両者は同一の universality class に属する。
ここで一致が主張されているのは値そのものではなく、unfolded local statistics である点を強調しておく必要がある。
Z.28.7 Pair Correlation Function
最も重要な比較対象は pair correlation function である。これを R_2(x) = ⟨ Σ_{i≠j} δ(x - (x_i - x_j)) ⟩ と定義する。
GUE の極限形として知られる Montgomery–Dyson の予測は R_2(x) = 1 - (sin(πx)/(πx))^2 である。
既知の数値的結果によれば、
リーマンゼータ零点はこの式に極めてよく一致する
GUE bulk 統計も同じ式に従う
従って、本理論において決定的な比較命題は次の形で表現される。
主張 Z.28.2
差分作用素 D_Δ の固有値に対しても R_2(x) → 1 - (sin(πx)/(πx))^2 が成立するならば、D_Δ とゼータ零点の局所スペクトル統計が同一の GUE universality class に属することが強く支持される。
Z.28.8 Number Variance
さらに、より長距離の統計量として number variance Σ^2(L) = ⟨ (N(L) - L)^2 ⟩ を考える。
GUE においては Σ^2(L) ∼ (1/π^2) log L という対数的増大が知られている。リーマンゼータ零点も同様の振る舞いを示すことが数値的に知られている。
従って D_Δ に対しても同一の対数的増大が確認されれば、局所統計だけでなく中距離統計においても両者の一致が確認されたことになる。
Z.28.9 スペクトル対応の意味
ここで重要なのは、一致しているのが 固有値そのもの ではないという点である。一致が主張されるのは、 局所統計構造 である。
従って本節の結論は、 「ゼータ零点が D_Δ の固有値そのものである」 という強い同一視ではなく、 「ゼータ零点と D_Δ は同一のスペクトル生成機構に属する」 という universality class の主張として理解されるべきである。
Z.28.10 構造的結論
以上の結果を統合すると、次の対応図式が得られる。
差分ダイナミクス側
Δ → D_Δ → 固有値列 {ξ_n} → GUE 統計
ゼータ側
ζ(s) → 零点列 {γ_n} → GUE 統計
したがって、両者は異なる対象ではあるが、少なくとも局所スペクトル統計のレベルでは同一の普遍構造を共有している。
Z.28.11 本論文で主張しうる最も強い命題
以上を踏まえれば、本論文で現時点において主張しうる最も強い命題は次である。
命題 Z.28.3
差分作用素 D_Δ = 2i[H,A] のスペクトルは、非可積分領域において GUE universality class に属する。また、リーマンゼータ零点も同一の universality class に属する。従って両者は、同一のスペクトル生成機構の異なる実現として解釈される。
この命題は、固有値の値の一致ではなく、局所統計構造の一致に関する主張である。
Z.28.12 「主観の影」の最終定式化
ここで初めて、「ゼータは主観の影なのか」という問いに対して、数学的に明確な形で応答することができる。
厳密には、ゼータそのものが主観であるわけではない。しかし、主観差分構造を作用素化して得られるスペクトル統計と、ゼータ零点の統計は一致する。
この意味で、より正確な対応は次の通りである。
主観：非可換差分流
ゼータ：そのスペクトル統計的像
Z.28.13 連続と離散の問題
しばしば問題となるのは、「主観は連続的であるのに、ゼータは離散的である」という点である。この違いは、本理論では スペクトル化（作用素化） によって説明される。
すなわち、連続的差分流 Δ は、自己共役作用素 D_Δ を経由して固有値列へ写される。ここで
主観：連続的差分ダイナミクス
ゼータ：そのスペクトル統計
という階層が成立する。従って、連続と離散の差異は矛盾ではなく、表現レベルの違いである。
Z.28.14 結論
本節で得られた結論は以下の通りである。
D_Δ のスペクトルとゼータ零点は、値そのものではなく統計量を通じて比較可能である
unfolding により両者を同一スケール上で比較できる
spacing, pair correlation, number variance はいずれも GUE universality class を指標とする
差分作用素 D_Δ の局所統計が GUE に従うなら、ゼータ零点との統計的一致が成立する
従って、ゼータは主観差分ダイナミクスのスペクトル統計的像として理解される
この意味で、本節は数学的対象としてのゼータ、物理的対象としての量子カオス、そして哲学的対象としての主観差分構造を、一つのスペクトル統計的枠組みにおいて接続する役割を果たす。
Section Z.29
Reconstruction of Zeta from the Difference Operator
Z.29.1 目的
本節の目的は、差分作用素 Σ_Δ = Δ†Δ のスペクトルから、ζ(s) 型関数を構成する枠組みを与えることである。
ここで重要なのは、本節の目標がリーマンゼータ関数そのものとの即時的な一致を主張することにはないという点である。本節で問題とするのは、むしろ
同型の極構造
同型の周期構造
同型のスペクトル生成機構
を持つ関数を、差分作用素から再構成できるかどうかである。
したがって本節の位置づけは、差分ダイナミクスからゼータ型構造への再構成可能性を定式化することにある。
Z.29.2 スペクトルゼータの定義
Σ_Δ の固有値を λ_n > 0 とする。このとき、対応するスペクトルゼータ関数を ζ_spec(s) = Σ_n λ_n^{-s} と定義する。
これは Σ_Δ に付随する自然なスペクトル生成関数であり、本理論における最も基本的なゼータ型対象である。
Z.29.3 Heat Kernel 表現
対応する heat trace を Θ(t) = Tr(e^{-t Σ_Δ}) = Σ_n e^{-t λ_n} と定義する。
このとき、Mellin 変換により ζ_spec(s) = (1/Γ(s)) ∫_0^∞ t^{s-1} Θ(t) dt が成立する。これはスペクトルゼータを heat kernel から構成する標準的な経路であり、本理論における ζ 型関数の第一の再構成法を与える。
Z.29.4 半古典展開と周期軌道
他方、heat trace は semiclassical regime において周期軌道展開を持つことが期待される。形式的には Θ(t) ≈ Θ_smooth(t) + Σ_γ A_γ e^{-t S_γ} と書かれる。
ここで
γ は primitive periodic orbit
S_γ はその作用
A_γ は安定性および Maslov 指数を含む振幅
を表す。
この表示により、heat trace は単なるスペクトル和ではなく、周期構造の圧縮表現としても理解される。
Z.29.5 Dynamical Zeta
周期軌道データからは、対応する dynamical zeta を Z_dyn(s) = ∏_γ (1 - e^{-s T_γ})^{-1} として定義することができる。ここで T_γ は軌道 γ の周期である。
この関数は、周期軌道長の列 {T_γ} を圧縮して保持する生成関数である。
Z.29.6 ζ_spec と Z_dyn の関係
Dynamical zeta の対数を取ると log Z_dyn(s) = Σ_γ Σ_{k=1}^∞ (1/k) e^{-k s T_γ} を得る。
一方、heat trace は Θ(t) = Σ_n e^{-t λ_n} であり、Mellin 変換および Laplace 変換を通じてスペクトルと接続している。
従って、本理論においては 周期軌道 ↔ heat trace ↔ スペクトル という連鎖が成立し、dynamical zeta と spectral zeta は同一構造の異なる圧縮表現として理解される。
Z.29.7 対応構造
以上を整理すると、本理論における対応は次のように与えられる。
構造
本理論における対象
固有値
λ_n
スペクトルゼータ
ζ_spec(s)
熱核
Θ(t)
周期軌道
γ
Dynamical zeta
Z_dyn(s)
この表は、差分作用素のスペクトル側と周期軌道側とが、ζ 型構造の二つの表現を与えることを示している。
Z.29.8 リーマンゼータとの一致条件
ζ_spec(s) または Z_dyn(s) をリーマンゼータと一致させるためには、少なくとも次の条件が必要となる。
条件 1：スペクトル構造
Σ_Δ の固有値が λ_n ≈ γ_n^2 + 1/4 を満たすこと。ここで γ_n はリーマンゼータの非自明零点の虚部である。これは Hilbert–Pólya 型条件に相当する。
条件 2：周期構造
周期軌道の周期が T_γ = log p に対応すること。ここで p は素数である。
条件 3：振幅構造
軌道振幅が A_γ ∼ p^{-k/2} のような形で素数重みと整合すること。
これらが成立すれば log ζ(s) = Σ_p Σ_{k=1}^∞ (1/k) p^{-ks} と同型の展開が得られる。
Z.29.9 本理論における実現可能性
本理論では、差分ダイナミクスにおいて
周期軌道 γ は Δ の周期運動に対応し
T_γ はその内部時間周期
S_γ はその作用
として定義される。
したがって、リーマンゼータ型構造との接続における鍵は、 周期軌道長の分布が log p 型の統計構造を持つかどうか という点にある。
Z.29.10 重要な観察
差分作用素は Δ = 2[H,A] によって与えられ、これは本質的に非可換流である。このとき、その周期構造は Lie algebra 上の flow の閉軌道として現れる。
一般に、準周期系、Anosov 系、量子カオス系などにおいては、周期軌道の本数は長さに関して指数的に増加することが知られている。従って、本理論における差分流が十分に複雑なカオス性を持つならば、その周期長分布が素数長に類似した統計構造を持つ可能性は排除されない。
Z.29.11 仮説 Z.29.1
以上を踏まえ、次の仮説を導入する。
仮説 Z.29.1（周期軌道長の素数型分布）
差分ダイナミクスが適切なカオス性を持つ場合、その周期軌道長 T_γ は統計的に T_γ ≈ log p に対応する分布を持つ。
これは現段階では未証明であるが、dynamical zeta と数論的ゼータとの接続における中心的仮説である。
Z.29.12 構成スキーム
以上の議論をまとめると、差分作用素から ζ 型構造を再構成するための手順は次のように整理される。
Step 1
Δ から正作用素 Σ_Δ = Δ†Δ を構成する。
Step 2
その固有値列 {λ_n} を得る。
Step 3
heat trace Θ(t) = Tr(e^{-t Σ_Δ}) を定義する。
Step 4
Mellin 変換により ζ_spec(s) を構成する。
Step 5
heat trace の周期軌道分解 Θ(t) ↔ Σ_γ A_γ e^{-t S_γ} を与える。
Step 6
そこから dynamical zeta Z_dyn(s) を構成し、ζ 型構造を得る。
このスキームにより、スペクトル側と周期軌道側の双方からゼータ型関数が再構成される。
Z.29.13 結論
本節で示した内容は次の通りである。
Δ†Δ は自然なスペクトルを持つ
そこから spectral zeta ζ_spec(s) を構成できる
heat trace の周期軌道展開から dynamical zeta を定義できる
適切な条件の下で、これらはリーマンゼータと同型の極構造および周期構造を持つ可能性がある
したがって、本理論において差分作用素は、ζ 型関数の再構成源として機能する。
Z.29.14 慎重な最終主張
以上を踏まえると、現段階で採るべき最も慎重かつ強い表現は次の通りである。
最終主張
ζ 関数は、差分ダイナミクスに由来するスペクトルゼータとして再構成される可能性がある。
ここでいう「再構成」とは、必ずしもリーマンゼータそのものとの厳密一致を意味するのではなく、極構造・周期構造・統計構造における同型性を含む広い意味での再構成である。
Z.29.15 「影」仮説の完成
最後に、本論文全体を貫く直観に立ち返ると、本理論における対応は
主観 → Δ（差分流）
ゼータ → ζ_spec（差分流のスペクトル関数）
として表される。
従って、ゼータは主観そのものではない。しかし、 主観差分構造をスペクトル化した像 として現れる可能性がある。この意味で、「ゼータは主観の影である」という表現は、本理論においては 差分ダイナミクスのスペクトル投影 という形で数学的に精密化される。
Section Z.30
Logical Status: Theorems, Derived Results, and Hypotheses
Z.30.1 目的
本節の目的は、本研究において得られた結果を、その論理的地位に応じて厳密に分類することである。具体的には、結果を
定理（rigorous）
準定理（既存理論への還元により支持されるもの）
仮説（未証明であるが構造的に強く示唆されるもの）
解釈的主張
に分けて整理する。
これにより、本研究において何が厳密に示され、何が既知理論に依存し、何が今後の課題として残されているかを明確にする。
Z.30.2 定理（厳密に成立する結果）
定理 Z.30.1（差分作用素の自己共役性）
Δ = 2[H,A], H,A Hermitian とする。このとき D_Δ := 2i[H,A] は自己共役である。
定理 Z.30.2（スペクトルゼータの存在）
Σ_Δ = Δ†Δ は正作用素であり、適切なスペクトル条件の下で ζ_spec(s) = Σ_n λ_n^{-s} は Re(s) > s₀ において収束する。
定理 Z.30.3（Heat Kernel 表現）
heat trace Θ(t) = Tr(e^{-t Σ_Δ}) に対して、 ζ_spec(s) = (1/Γ(s)) ∫_0^∞ t^{s-1} Θ(t) dt が成り立つ。
定理 Z.30.4（Born 型構造の厳密版）
Σ_Δ = Δ†Δ により ρ_Δ = Σ_Δ / Tr(Σ_Δ) を定義する。このとき、射影 P_i に対して p_i^(Δ) = Tr(P_i ρ_Δ) は完全に Born 型の確率構造を与える。
以上は、本研究の中で作用素論的・定義的に厳密に成立する部分である。
Z.30.3 準定理（既存理論に依存する結果）
準定理 Z.30.5（GUE 普遍性）
以下の条件を仮定する。
H のスペクトル分布が滑らかである
A が H と独立である
variance profile が bulk 領域で非退化である
このとき D_Δ = 2i[H,A] の局所スペクトル統計は、GUE universality class に属すると期待される。
この主張の根拠は、
generalized Wigner 行列
Dyson Brownian motion
local semicircle law
など既存のランダム行列理論にある。
準定理 Z.30.6（Dyson Flow 対応）
時間依存版 D_t = 2i[H_t, A_t] に対し、その固有値過程は Dyson–Brownian motion 型の進化に従うと期待される。
これもまた、行列値 Ornstein–Uhlenbeck 過程と Dyson flow に関する既存理論への還元によって強く支持される。
Z.30.4 強い仮説（本研究の中心仮説）
仮説 Z.30.7（周期軌道構造）
差分ダイナミクスの周期軌道 γ に対して、その周期長は統計的に T_γ ≈ log p の構造を持つ。
仮説 Z.30.8（スペクトル対応）
Σ_Δ の固有値 λ_n は、リーマンゼータの非自明零点虚部 γ_n に対して λ_n ≈ γ_n^2 + 1/4 に対応する。
仮説 Z.30.9（ゼータ再構成）
適切な条件の下で、スペクトルゼータ ζ_spec(s) はリーマンゼータ関数と同型の極構造を持つ。
これらは本研究の数論的側面および Hilbert–Pólya 的側面の中核をなすが、現段階では未証明である。
Z.30.5 解釈レベルの主張
本研究には、厳密定理や準定理とは別に、理論全体を貫く解釈的主張が存在する。主要なものは次の通りである。
主張 Z.30.10
確率は差分強度分布として導出される。
主張 Z.30.11
時間は差分ノルムの積分として内部的に生成される。
主張 Z.30.12
主観とは、非可換差分流の持続構造である。
これらは物理的・哲学的意味づけを与えるものであり、その多くは理論全体の構成から支持されるが、通常の意味での数学的定理ではない。
Z.30.6 「ゼータ＝主観の影」という命題の位置づけ
誤解を避けるために、ここで改めて本研究の立場を明確にする。
ゼータそのものが主観であるわけではない
Δ が主観差分構造である
Σ_Δ がそのスペクトル化である
ζ_spec がそのスペクトル関数である
従って、本研究が与える最も適切な表現は、 ゼータは差分構造のスペクトル的像として現れる可能性がある というものである。
Z.30.7 何が未証明であるか
現段階で未証明であるのは、明確に次の三点である。
周期軌道と素数との厳密対応
固有値列 λ_n とゼータ零点との厳密一致
ζ_spec = ζ の完全一致
これらは本理論の最も野心的な部分であり、今後の解析的・数値的課題として残されている。
Z.30.8 本研究の到達点
以上を総合すると、本研究が既に達成した内容は次のように整理できる。
主観の最小構造を差分原理として定式化した
差分構造から内部時間生成を導出した
差分構造から Born 型確率を導出した
差分構造からスペクトル・ゼータ構造への接続を与えた
GUE 型統計との対応枠組みを構築した
この意味で、本研究は主観・時間・確率・スペクトル・ゼータを、単一の非可換差分原理から統一的に記述するための枠組みを与えている。
Z.30.9 最終結論
本理論の最終的な主張は次のようにまとめられる。
本理論は、
主観（差分構造）
時間（内部生成）
確率（Born 型構造）
ゼータ（スペクトル構造）
を、単一の非可換差分原理から統一的に導出する枠組みを与える。
ただし、そのうち数論的同一視に関わる最強の主張はなお仮説の段階にあり、今後の検証を要する。
Final Section
Predictions: Experimental, Numerical, Number-Theoretic, and Structural
F.1 目的
本節では、本理論が与える予言を四つの水準に分けて整理する。
実験的予言
数値的予言
数論的予言
構造的予言（Born 則逸脱）
とりわけ最後の項目は、本理論が標準量子力学と実質的に分岐する決定的な点を含む。
F.2 実験的予言（量子系）
予言 F.2.1（非可換強度依存の時間生成）
内部時間生成率 dτ/dt = ||Δ|| = 2||[H,A]|| は、可観測な有効ダイナミクスの速度スケールに反映される。
すなわち、非可換性が大きい領域では内部時間の進行が速く、非可換性が減少する領域では時間進行が遅くなることが予測される。
予言 F.2.2（干渉抑制の起源）
干渉の消失は、外部環境との相互作用のみによって引き起こされるのではなく、 Δ の減少、すなわち非可換性の消失 によっても生じうる。
したがって、デコヒーレンス様の現象は差分構造の内部力学からも生じる。
予言 F.2.3（準閉鎖系における時間停止）
Δ → 0 に近づく極限では、 有効な時間進行が停止する領域 が現れる。
これは、通常の外部時間を固定した記述では捉えにくいが、本理論の内部時間の観点からは自然な予測である。
F.3 数値的予言（ランダム行列）
予言 F.3.1（GUE 普遍性）
差分作用素 D_Δ = 2i[H,A] の局所スペクトル統計は、非可積分領域において GUE に従う。
予言 F.3.2（Dyson Flow）
時間依存系においては、D_Δ の固有値流は Dyson–Brownian motion 型の進化を示す。
予言 F.3.3（臨界領域での逸脱）
特定の臨界条件の下では、局所スペクトル統計は GUE から逸脱し、場合によっては Poisson 型へのクロスオーバーが生じる。
これは、カオス領域と非カオス領域の境界に対応する可能性がある。
F.4 数論的予言
予言 F.4.1（ゼータとの統計的一致）
リーマンゼータ零点の局所統計と D_Δ のスペクトル統計は一致する。
ここでの一致は固有値そのものの一致ではなく、unfolding 後の spacing, pair correlation, number variance などの統計的一致を意味する。
予言 F.4.2（Trace Formula の存在）
Δ-flow に対して、Gutzwiller 型あるいはその適切な一般化としての trace formula が成立する。
予言 F.4.3（周期構造と素数）
特定の差分ダイナミクスにおいては、周期軌道長に対し T_γ ≈ log p が成立する可能性がある。
これは、周期構造と素数構造を直接接続する数論的予言である。
F.5 構造的予言（Born 則逸脱）
F.5.1 Born 則成立条件の再確認
Born 型重み w_i = Tr(P_i ρ) が成立するためには、少なくとも以下の条件が必要である。
(C1) 非可換性が十分に支配的であること
(C2) H のスペクトルギャップが安定であること
(C3) coarse-graining により干渉が抑制されていること
本理論では、Born 則は無条件の公理ではなく、一定の力学的条件のもとで有効に成立する法則として理解される。
F.5.2 N ≥ 3 における新自由度
2×2 系では、
ギャップ構造がほぼ固定され
Δ の自由度が低く
差分再配分の余地が限られている。
これに対して N ≥ 3 では、
固有値差 h_i - h_j が動的に変化し
A_ij が多方向に分岐し
Δ_ij = (h_i - h_j) A_ij が複数チャネル間で競合する
ため、差分 flow における再配分ダイナミクスが本質的に現れる。
F.5.3 Gap Collapse
定義
Δ-gap collapse とは、 |h_i - h_j| → 0 がダイナミクスの中で実際に生起する現象をいう。
このとき Δ_ij = 2(h_i - h_j) A_ij → 0 となるが、同時に他のチャネルにおいては差分強度が増大しうる。従って、flow 全体としては単なる消失ではなく、 差分強度の再配分 が生じる。
F.5.4 Born 則破れの機構
Born 型重みを w_i = Tr(P_i ρ) とし、差分重みを p_i^(Δ) = Φ_i / Σ_k Φ_k とする。通常は p_i^(Δ) ≈ w_i が成立する。
しかし gap collapse により Φ_i = Σ_j |h_i - h_j|^2 ||A_ij||^2 がチャネルごとに非一様に崩壊する場合、 p_i^(Δ) ≠ w_i が系統的に生じる。
したがって Born 則の破れは、ノイズによるランダムな逸脱ではなく、差分強度の再配分による構造的逸脱として理解される。
予言 F.5.1（Born 則逸脱）
N ≥ 3 の系において、固有値ギャップが動的に崩壊する領域では、Born 則からの系統的逸脱が発生する。
これは本理論が標準量子力学と実質的に異なる予測を与える最も重要な点である。
F.5.6 逸脱の特徴
この逸脱はランダムなノイズではなく、以下の特徴を持つ。
特徴 1：強い非可換チャネルへの偏り
確率は、差分強度の大きいチャネルへ系統的に偏る。
特徴 2：履歴依存性（hysteresis）
確率分布は瞬間的状態のみではなく、直前の差分 flow の履歴に依存する。
特徴 3：臨界揺らぎ
gap collapse の近傍では、確率分布に臨界的揺らぎが現れる。
F.5.7 相構造
以上により、本理論は少なくとも次の三つの相を区別する。
(I) Born 相
非可換性が支配的であり、ギャップが安定している。標準量子力学が有効に成立する相である。
(II) 臨界相
ギャップ崩壊が始まり、確率揺らぎが増大する相である。
(III) 非 Born 相
差分再配分が支配的となり、確率構造そのものが変形する相である。
F.5.8 実験的含意
この種の逸脱は、少なくとも原理的には
多準位量子系
強結合開放系
フィードバックを持つ量子系
において観測可能である可能性がある。
F.5.9 数値検証
この予測は数値的にも検証可能である。例えば
N ≥ 5
H(t), A(t) を時間発展させる
p_i^(Δ) と Born 重み w_i を比較する
というシミュレーションにより、差分再配分と Born 則逸脱の有無を確認できる。
F.5.10 数論的含意
もしゼータ構造が本理論の差分ダイナミクスに深く対応しているならば、零点統計においても何らかの臨界揺らぎが存在する可能性がある。
もちろん、これは現段階では強い推測に留まるが、数論的側面と差分臨界現象とを結ぶ興味深い可能性を示している。
F.6 最終結論
本理論の最大の予言は次の通りである。
結論
量子力学の確率構造、特に Born 則は普遍的な根本法則ではなく、差分構造およびスペクトルギャップに依存する有効法則である。そしてその破れは、高次元差分ダイナミクスにおいて必然的に現れうる。
この意味で、
量子力学は特殊な相であり
Born 則は有効法則であり
主観差分構造の方がより基礎的である
という理論的像が導かれる。
Epilogue
Conceptual Position, Physical Implications, and Open Problems
E.1 本理論の位置づけ
本研究で構築した枠組みは、次の三領域を横断する。
数学（非可換構造・スペクトル理論・ゼータ関数）
物理（量子力学・拘束系・カオス）
哲学（主観・時間・存在）
しかし本理論の出発点は、これらのいずれにも還元されない。それは 非可換差分構造 Δ = 2[H,A] そのものである。
E.1.1 基本的統一像
本理論では、以下の諸概念が単一の差分構造から統一的に生成される。
主観：Δ ≠ 0 を持つ内部自由度
時間：τ = ∫ ||Δ|| dt
確率：Σ_Δ = Δ†Δ に基づく分布
ゼータ：スペクトル圧縮 ζ_spec
すなわち、本理論は主観・時間・確率・ゼータを、別個の原理から導入するのではなく、同一の差分原理から派生する量として再構成する。
E.2 物理的帰結
E.2.1 時間の再定義
従来の理論では、時間は外部から与えられた独立パラメータとして扱われることが多い。これに対して本理論では、 時間 = 差分の累積 として定義される。
したがって、
Δ = 0 なら時間進行は停止し
Δ ≠ 0 なら時間が生成される
という構造が成立する。
E.2.2 量子力学の位置づけ
本理論において標準量子力学は、差分ダイナミクスのある特定の相、すなわち Born 相における有効理論として位置づけられる。
従って、量子力学は本理論の観点からは
普遍的な最終理論ではなく
条件付きで成立する有効理論
である。
E.2.3 確率の起源
本理論では、確率は公理として導入されない。確率とは、差分 flow の再分配構造から生じる有効量である。
この点は、確率を測定公理の一部として受け入れる標準的量子論と根本的に異なる。
E.3 数学的帰結
E.3.1 スペクトル理論
Δ†Δ により、自然な正自己共役作用素が構成される。その結果、
離散スペクトル
spectral zeta
heat kernel
trace formula の候補構造
が一貫して導入される。
E.3.2 ゼータ関数の位置づけ
本理論の立場では、ζ(s) は
主観そのものではなく
単なる数論的対象でもなく
差分ダイナミクスのスペクトル圧縮 として現れる。
この点が、本研究におけるゼータ関数理解の中心である。
E.3.3 Hilbert–Pólya の再解釈
従来、Hilbert–Pólya の発想では「未知の自己共役作用素」の存在が仮定されてきた。本理論では、その候補として D_Δ = 2i[H,A] が自然に現れる。
もちろん、これが直ちにリーマン零点を与えることはまだ示されていない。しかし、少なくとも
自己共役性
GUE 型統計
周期軌道構造
という要請を一つの作用素論的枠組みの中で満たす候補を与えている。
E.4 哲学的帰結
E.4.1 主観の定義
本理論において主観とは、 内部自由度 + 非可換差分流 として定義される。
これは「主観」を内観的・記述的概念としてではなく、力学的構造として捉える試みである。
E.4.2 客観との関係
客観は、射影 π によって得られる coarse-grained な像として理解される。これに対して主観は、その背後のファイバー方向、すなわち内部差分の方向に属する。
従って本理論では、客観と主観は対立概念ではなく、 同一構造の射影像と内部構造 として関係づけられる。
E.4.3 存在論的含意
本理論は存在を少なくとも二層に分ける。
連続的側面：Δ（becoming）
離散的側面：スペクトル（being）
このときゼータは後者、すなわち being の側に属する。差分流そのものが becoming を担い、そのスペクトル圧縮が being として現れる。
E.5 本理論の限界
E.5.1 未証明部分
現段階で未証明の核心部分は次の通りである。
Δ-flow に対する trace formula の厳密証明
ζ_spec とリーマンゼータとの完全一致
周期軌道と素数との厳密対応
E.5.2 モデル依存性
高次元モデルの構築はまだ完全ではなく、カオス性や双曲性の厳密条件も十分には解明されていない。従って、本理論の多くの結論は現時点ではモデル依存的である。
E.5.3 物理実験
Born 則逸脱の直接観測はまだ実現されていない。したがって、物理的帰結の一部は今後の実験的検証に委ねられている。
E.6 今後の課題
E.6.1 数学的課題
Δ-flow の双曲性の証明
trace formula の厳密化
ζ_spec = ζ が成立する条件の確立
E.6.2 物理的課題
N ≥ 3 の実験系の構築
Born 則逸脱の測定
内部時間の可観測量としての定式化
E.6.3 数値的課題
GUE 統計確認の高精度化
Dyson flow の直接検証
ゼータ零点との統計比較の高精度化
E.7 最終結論
本研究の最終的な結論は次のように要約される。
結論
世界の基礎構造は、非可換差分である。そこから
時間が生まれ
確率が生まれ
スペクトルが生まれ
ゼータ関数がそのスペクトル像として現れる
この意味で、本理論は主観・物理・数学を結ぶ単一原理として差分構造を提案する。
E.8 結語
本理論は、
主観を特権化するものではなく
物理へ単純還元するものでもなく
数学の内部に閉じるものでもない
むしろ、それらすべてを 差分構造 へ還元する立場を採る。
これが本研究の基本的立場であり、その意義は、主観・時間・確率・ゼータを別々の問題としてではなく、一つの非可換差分原理の異なる表現として捉え直した点にある。
EpiloguePlus+
なぜゼータでなければならなかったのか
本研究では、主観を非可換差分ダイナミクスとして定式化し、その正値化作用素 Σ_Δ = Δ†Δ のスペクトル構造を通じて、ゼータ関数的構造が自然に現れることを示した。しかし、この出現は単なる形式的類似や偶然の一致として理解されるべきではない。本節の目的は、なぜこの構造がゼータ型に帰着せざるをえないのか、その必然性を明確にすることである。
本理論において主観は、連続的かつ非可換な生成過程 Δ によって与えられる。この生成過程は、周期軌道や固有モードといった基本成分へ分解され、それらは独立な生成単位として振る舞う。他方、観測や記述の段階では、これらの構造を可換的な形式へ圧縮する必要がある。この圧縮は、単なる情報削減ではない。むしろ、生成構造を保持したまま、それを再構成可能な形式へ写す操作である。
したがって、この圧縮形式には少なくとも次の条件が課される。第一に、基本成分の独立性を保持し、それらの合成構造を保存できること。第二に、周期性および反復性を正確に符号化できること。第三に、圧縮後の形式から元の生成構造を再展開できること。第四に、内部時間の再パラメータ化に対して不変であること。
これらの要請を同時に満たす形式は強く制約される。特に、独立成分の合成が可換的に扱われ、反復構造が指数的に組み込まれる場合、生成関数は本質的に「対数が加法的であり、関数自体は乗法的である」という形を取らざるをえない。すなわち、基本成分の集合 {γ} に対して log Z(s) = Σ_γ Σ_{r≥1} e^{-s r T_γ} / r が現れ、これを指数化すると Z(s) = ∏_γ (1 - e^{-s T_γ})^(-1) を得る。
この表式は、独立な生成要素とその反復を保持したまま全体構造を圧縮する最小形式であり、dynamical zeta 関数の一般形そのものである。ここで重要なのは、この形が特定の対象に固有なものではなく、独立性・反復性・再構成可能性という要請の帰結として現れる点である。
同様の構造はスペクトル側からも現れる。正作用素 Σ_Δ の固有値列 {λ_n} に対して、スペクトルゼータ関数 ζ_Δ(s) = Σ_n λ_n^(-s) を考えると、これは連続的な差分強度を離散スペクトルへ圧縮し、その全体構造を単一の関数として符号化する役割を担う。さらに、heat kernel trace および Mellin 変換を通じて、このスペクトル表現は周期軌道側の生成関数と接続される。したがって、連続ダイナミクスと離散スペクトルを統一的に記述する形式として、ゼータ型関数が現れるのは自然である。
さらに重要なのは、この帰結が特定の力学系や特定の対象に依存しないことである。独立性、反復性、再構成可能性、スケール不変性という条件が同時に課される限り、生成関数は本質的にゼータ型へ収束する。言い換えれば、ゼータ関数は特定対象に固有の偶然的構造ではなく、圧縮と再構成の要請そのものから導かれる普遍形式である。
この観点からすれば、ゼータ関数は主観差分ダイナミクスそのものではない。主観はあくまで連続的で非可換な生成過程であり、ゼータ関数はその過程が可換的かつ離散的に圧縮されたときに現れる構造である。したがって両者は同一ではないが、前者が存在し、その生成構造を保持した圧縮が要請される限り、後者は必然的に現れる。
以上より、本研究におけるゼータ関数の出現は、特定の理論的選択や外在的な付加によるものではなく、生成構造を保持したまま圧縮するという要請から必然的に導かれる。したがって、「なぜゼータでなければならなかったのか」という問いに対する本研究の答えは、次のように要約される。
独立な生成要素とその反復構造を保持したまま、連続的な非可換ダイナミクスを可換的に圧縮し、しかも再構成可能な形式で記述するためには、ゼータ型生成関数以外の選択肢は本質的に存在しない。
