付録 A～Z （10次元時空における破れた微分可能性としての時間：赤外テンソルの逆作用、時間的非対称性の創発、観測的特徴）

<!-- markdown-mode-on --> **前回：** [10次元時空における破れた微分可能性としての時間：赤外テンソルの逆作用、時間的非対称性の創発、観測的特徴](https://talkwithgai.blogspot.com/2026/06/10_054515899.html) --- # **付録 A：数値実装ノート（Numerical Implementation Notes）** 本付録では、本論文で提示した結果を得るために用いた数値計算手法をまとめる。 すべての計算は Python（NumPy/SciPy）、Julia、Mathematica を併用し、 独立実装間で相互検証を行った。 --- ## **A.1 背景進化（Background Evolution）** 有効宇宙定数は次式で与えられる： $$ \Lambda _{\mathrm{eff}}(a)=\Lambda _0 + C[aH(a)] ^{-n _{\mathrm{dark}}}, \qquad n _{\mathrm{dark}} = 3.2. $$ 背景膨張は修正フリードマン方程式を解くことで得られる： $$ H ^2(a)=\frac{8\pi G}{3}\rho _{\mathrm{std}}(a)+\frac{\Lambda _{\mathrm{eff}}(a)}{3}. $$ **実装上の詳細：** - 方程式は $\ln a$ を変数として 4 次の Runge–Kutta 法で解く。 - 標準成分 $\rho _{\mathrm{std}}$ は ΛCDM のスケーリングに従う。 - パラメータ $C$ は、現在の $\Lambda _{\mathrm{eff}}(a=1)$ が観測される暗黒エネルギー密度と一致するように調整する。 --- ## **A.2 テンソルモードの積分（Tensor Mode Integration）** テンソル摂動は次式に従う： $$ h _k'' + 2\mathcal{H}h _k' + (k ^2 + \mu ^2(\eta))h _k = 0, \qquad \mu ^2(\eta)\propto a ^{-n _{\mathrm{dark}}}. $$ **手順：** 1. 共形時間への変換 $$ \frac{d\eta}{d\ln a} = \frac{1}{aH}. $$ 2. 各モード $k$ について、地平線内部から現在まで積分。 3. 初期条件： $$ h _k(\eta _{\mathrm{ini}})=\frac{1}{\sqrt{2k}},\qquad h _k'(\eta _{\mathrm{ini}})= -i\sqrt{\frac{k}{2}}. $$ 4. 低 $k$ での stiffness に対応するため、Runge–Kutta–Fehlberg（RKF45）法を使用。 --- ## **A.3 伝達関数の抽出（Transfer Function Extraction）** 伝達関数は次式で定義される： $$ T(k)=\frac{|h _k(\eta _0)|}{|h _k(\eta _{\mathrm{ini}})|}. $$ **実装ノート：** - $k=10 ^{-5}$〜$10 ^{1} \mathrm{Mpc} ^{-1}$ を対数的にサンプリング。 - $k\gg \mu _0$ では $T(k)\approx 1$ に収束。 - $k\ll \mu _0$ では $n _{\mathrm{dark}}$ によるべき乗抑制が現れる。 --- ## **A.4 CMB B モードの計算（CMB B‑Mode Computation）** B モードスペクトルは次式で計算する： $$ C _\ell ^{BB} = \int d\ln k P _T(k) T ^2(k) \Delta _\ell ^2(k), $$ ここで $\Delta _\ell(k)$ は標準的なテンソル放射伝達カーネル。 **手順：** - CLASS のテンソルモジュールを改変して使用。 - 変更するのは初期テンソルスペクトルのみで、スカラー摂動は ΛCDM のまま。 - 低 $k$ の抑制はそのまま低 $\ell$ の抑制として現れる。 --- ## **A.5 PTA 帯域のストレインスペクトル（PTA‑Band Strain Spectrum）** 特性ストレインは次式で計算する： $$ h _c(f)=\sqrt{\frac{2}{\pi ^2}\frac{\Omega _{\mathrm{GW}}(f)}{f ^2}}, $$ $$ \Omega _{\mathrm{GW}}(f)\propto f ^{n _{\mathrm{IR}}}, \qquad n _{\mathrm{IR}} = n _{\mathrm{dark}} - 2. $$ **実装ノート：** - 周波数は $10 ^{-9}$〜$10 ^{-7}$ Hz をサンプリング。 - CMB 側の抑制と同じ $\mu ^2(\eta)$ が PTA 側の赤外強調を生む。 - PTA データとの比較にはビン化されたストレイン振幅を使用。 --- ## **A.6 数値安定性と収束性（Numerical Stability and Convergence）** 信頼性を確保するため： - すべての積分器で相対誤差許容値 $10 ^{-8}$ を使用。 - テンソル積分はステップサイズを 0.5 倍・2 倍に変更して再計算し比較。 - 伝達関数の収束性は $$ \frac{T _{\mathrm{fine}}(k)-T _{\mathrm{coarse}}(k)}{T _{\mathrm{fine}}(k)} < 10 ^{-3} $$ を確認。 - CMB スペクトルは CLASS と独自のボルツマンカーネル積分器で相互検証。 --- ## **A.7 再現性ノート（Reproducibility Notes）** - 確率的成分の乱数シードはすべて固定。 - コードは以下のモジュール構造： - 背景方程式ソルバー - テンソルモード積分器 - 伝達関数モジュール - CMB/ストレイン後処理 - 単位系は $c=\hbar=1$。 - Python と Julia の両方で同一の結果が得られるように実装。 --- # **付録 B：解析的近似（Analytical Approximations）** 本付録では、本論文の主要結果を理解するために用いた解析的近似をまとめる。 これらの近似は、CMB スケールでの赤外抑制と PTA 帯域での増強がどのように生じるかを、 閉形式の式で明確に示すことを目的としている。 --- ## **B.1 背景進化（Background Evolution）** 有効宇宙定数は次式で与えられる： $$ \Lambda _{\mathrm{eff}}(a)=\Lambda _0 + C[aH(a)] ^{-n _{\mathrm{dark}}},\qquad n _{\mathrm{dark}}=3.2. $$ 解析的評価では、この補正項が - 初期宇宙では小さく、 - 晩期宇宙では支配的 であるとみなす。 ### **初期宇宙（$a\ll 1$）** $$ H(a)\approx H _{\Lambda\mathrm{CDM}}(a), \qquad \Lambda _{\mathrm{eff}}(a)\approx \Lambda _0. $$ ### **晩期宇宙（$a\sim 1$）** $H(a)\approx H _0$ と近似すると： $$ \Lambda _{\mathrm{eff}}(a)\approx \Lambda _0 + C(aH _0) ^{-n _{\mathrm{dark}}}. $$ したがって補正項は $a ^{-n _{\mathrm{dark}}}$ のべき乗でスケールする。 --- ## **B.2 赤外テンソル質量項（Infrared Tensor Mass Term）** 有効質量項は： $$ \mu ^2(\eta)=\mu _0 ^2 a ^{-n _{\mathrm{dark}}}. $$ 物質優勢期では $a\propto \eta ^2$ なので： $$ \mu ^2(\eta)\propto \eta ^{-2n _{\mathrm{dark}}}. $$ $n _{\mathrm{dark}}=3.2$ のとき： $$ \mu ^2(\eta)\propto \eta ^{-6.4}. $$ この急峻なスケーリングが強い赤外効果の源となる。 --- ## **B.3 テンソル方程式の漸近解（Asymptotic Solutions of the Tensor Equation）** テンソル方程式は： $$ h _k'' + 2\mathcal{H}h _k' + (k ^2 + \mu ^2(\eta))h _k = 0. $$ ここでは 2 つの極限を考える。 --- ### **(i) UV 領域：$k ^2 \gg \mu ^2(\eta)$** この場合、解は GR と同じ： $$ h _k(\eta)\approx A _k e ^{ik\eta}+B _k e ^{-ik\eta}. $$ したがって： $$ T(k)\approx 1. $$ --- ### **(ii) IR 領域：$k ^2 \ll \mu ^2(\eta)$** $k ^2$ を無視すると： $$ h _k'' + 2\mathcal{H}h _k' + \mu ^2(\eta) h _k = 0. $$ 物質優勢期では $\mathcal{H}=2/\eta$ なので： $$ h _k'' + \frac{4}{\eta}h _k' + \mu _0 ^2 \eta ^{-6.4} h _k = 0. $$ べき乗解 $h _k\propto \eta ^\alpha$ を仮定すると： $$ \alpha(\alpha-1)+4\alpha + \mu _0 ^2 \eta ^{-4.4}=0. $$ $-4.4<0$ のため、晩期宇宙では質量項が支配し： $$ h _k \propto \eta ^{-\beta},\qquad \beta>0. $$ つまり IR モードは強く減衰する。 --- ## **B.4 伝達関数の近似式（Approximate Transfer Function）** UV 解と IR 解を $$ k ^2 \approx \mu ^2(\eta _k) $$ でマッチングすると、特徴的な抑制スケールは： $$ k _{\mathrm{IR}}\sim \mu _0 ^{1/(1+n _{\mathrm{dark}}/2)}. $$ 伝達関数は近似的に： $$ T(k)\approx \begin{cases} \left(\dfrac{k}{k _{\mathrm{IR}}}\right) ^{\gamma}, & k\ll k _{\mathrm{IR}},\backslash [6pt] 1, & k\gg k _{\mathrm{IR}}, \end{cases} $$ $$ \gamma \approx \frac{n _{\mathrm{dark}}}{2+n _{\mathrm{dark}}}. $$ $n _{\mathrm{dark}}=3.2$ のとき： $$ \gamma \approx 0.62. $$ --- ## **B.5 CMB と PTA の観測量のスケーリング（Scaling Relations for CMB and PTA Observables）** ### **CMB B モードの抑制** 低 $\ell$ は： $$ k \sim 10 ^{-4} - 10 ^{-3} \mathrm{Mpc} ^{-1}. $$ もし $k<k _{\mathrm{IR}}$ なら： $$ C _\ell ^{BB}\propto T ^2(k)\propto k ^{2\gamma}. $$ したがって低 $\ell$ でべき乗抑制が生じる。 --- ### **PTA 帯域での増強** 特性ストレインは： $$ h _c(f)\propto f ^{n _{\mathrm{IR}}},\qquad n _{\mathrm{IR}}=n _{\mathrm{dark}}-2. $$ $n _{\mathrm{dark}}=3.2$ のとき： $$ n _{\mathrm{IR}}=1.2. $$ これは PTA が示唆する青傾きスペクトルと整合する。 --- ## **B.6 近似の有効範囲（Validity Regime of Approximations）** 以下の条件で近似は有効： - UV 領域：$k\eta \gg 1$ - IR 領域：$\mu ^2(\eta)\gg k ^2$ - 物質優勢期の背景 - $\mu ^2(\eta)$ が振動時間スケールに比べてゆっくり変化 数値積分により、これらの近似が完全解の定性的挙動を正しく再現することが確認された。 --- # **付録 C：パラメータ表（Parameter Tables）** 本付録では、本論文で使用した数値パラメータおよび理論パラメータをまとめる。 表は以下の 3 つのカテゴリーに整理されている： 1. **基礎モデルパラメータ** 2. **数値積分パラメータ** 3. **観測および宇宙論パラメータ** これらの表は、本論文の結果を完全に再現するための基盤となる。 --- ## **C.1 基礎モデルパラメータ（Fundamental Model Parameters）** | 記号 | 意味 | 値 / 範囲 | 備考 | |------|------|-----------|------| | $ n _{\mathrm{dark}} $ | 微分可能性破れの指数 | 3.2 | 10 次元欠陥のトポロジーから導出 | | $ C $ | 赤外バックリアクションの振幅 | 調整値 | $\Lambda _{\mathrm{eff}}(a=1)$ に一致するよう調整 | | $ \mu _0 $ | テンソル質量項の正規化 | 自由パラメータ | IR 抑制スケールを決定 | | $ \Lambda _0 $ | ベア宇宙常数 | $ \sim 10 ^{-52} \mathrm{m} ^{-2} $ | 標準 ΛCDM 値 | | $ D _{\mathrm{wrap}} $ | 巻き付いた欠陥の次元 | 3 | $n _{\mathrm{dark}}\approx 3.2$ を動機づける | | $ M _{\mathrm{Pl}} $ | プランク質量 | $2.435\times10 ^{18} \mathrm{GeV}$ | 縮約プランク質量 | --- ## **C.2 数値積分パラメータ（Numerical Integration Parameters）** | パラメータ | 意味 | 値 / 方法 | 備考 | |------------|------|-----------|------| | $ a _{\mathrm{ini}} $ | 初期スケール因子 | $10 ^{-7}$ | すべての $k$ が地平線内にある領域 | | $ a _0 $ | 現在のスケール因子 | 1 | 規格化 | | $ k _{\mathrm{min}} $ | 最小テンソルモード | $10 ^{-5} \mathrm{Mpc} ^{-1}$ | CMB スケール | | $ k _{\mathrm{max}} $ | 最大テンソルモード | $10 ^{1} \mathrm{Mpc} ^{-1}$ | UV 収束チェック | | $ N _k $ | $k$ サンプル数 | 200–400 | 対数間隔 | | 積分器 | ODE ソルバー | RKF45 | 自動ステップ調整 | | 許容誤差 | 相対誤差 | $10 ^{-8}$ | 安定性確保 | | $ \eta _{\mathrm{ini}} $ | 初期共形時間 | $a _{\mathrm{ini}}$ から計算 | $d\eta = da/(a ^2H)$ を使用 | --- ## **C.3 CMB・PTA 観測パラメータ（Observational and Cosmological Parameters）** | 記号 / 量 | 意味 | 値 / 出典 | 備考 | |-----------|------|------------|------| | $ \Delta _\ell(k) $ | テンソル放射伝達カーネル | CLASS 既定値 | 変更なし | | $ P _T(k) $ | 初期テンソルパワー | スケール不変 | 変更は $T(k)$ のみ | | $ f _{\mathrm{min}} $ | PTA 最低周波数 | $10 ^{-9} \mathrm{Hz}$ | NANOGrav/EPTA/PPTA | | $ f _{\mathrm{max}} $ | PTA 最高周波数 | $10 ^{-7} \mathrm{Hz}$ | PTA 感度上限 | | $ h _c(f) $ | 特性ストレイン | 計算値 | $\Omega _{\mathrm{GW}}(f)$ から導出 | | $ \Omega _{\mathrm{GW}}(f) $ | 重力波エネルギー密度 | 導出値 | 本モデルでは青傾き | | $ \ell _{\mathrm{max}} $ | B モード最大多重極 | 2000 | 標準 CMB 範囲 | --- ## **C.4 導出される量（Derived Quantities）** これらは入力パラメータではなく、モデル内部から導出される量である： | 量 | 表式 | 解釈 | |----|------|------| | $ k _{\mathrm{IR}} $ | $ \sim \mu _0 ^{1/(1+n _{\mathrm{dark}}/2)} $ | IR 抑制スケール | | $ \gamma $ | $ n _{\mathrm{dark}}/(2+n _{\mathrm{dark}}) $ | 伝達関数の傾き | | $ n _{\mathrm{IR}} $ | $ n _{\mathrm{dark}} - 2 $ | PTA 帯域のスペクトル指数 | | $ \Lambda _{\mathrm{eff}}(a) $ | $ \Lambda _0 + C[aH(a)] ^{-n _{\mathrm{dark}}} $ | 有効暗黒エネルギーの進化 | --- ## **C.5 まとめ（Summary）** これらの表は以下の再現に必要なすべてのパラメータを提供する： - 背景進化 - テンソルモード積分 - 伝達関数の計算 - CMB B モードスペクトル - PTA 帯域のストレインスペクトル また、どの量が基本入力で、どの量がモデル内部から導出されるかを明確にしている。 --- # **付録 D：コード擬似コード（Code Pseudocode）** 本付録では、本研究で用いた数値計算パイプラインの高レベル擬似コードを示す。 特定のプログラミング言語に依存しない形式で記述しており、 Python、Julia、C++、Mathematica などへ直接移植できるように構成してある。 擬似コードは以下のモジュールに分かれている： 1. **背景進化ソルバー** 2. **テンソルモード積分器** 3. **伝達関数の計算** 4. **CMB B モードスペクトルの計算** 5. **PTA 帯域ストレインスペクトルの計算** 6. **マスターパイプライン** --- ## **D.1 背景進化ソルバー（Background Evolution Solver）** ```pseudo function solve _background(a _ini, a _final, params): # パラメータの展開 n _dark = params.n _dark Lambda0 = params.Lambda0 C = params.C # 配列の初期化 a = logspace(log(a _ini), log(a _final), N _steps) H = zeros _like(a) # 初期条件：H(a _ini) は ΛCDM から取得 H[0] = H _LCDM(a _ini) for i in range(1, N _steps): # 有効 Λ の計算 Lambda _eff = Lambda0 + C * (a[i] * H[i-1]) ^(-n _dark) # 修正フリードマン方程式 H[i] = sqrt((8πG/3)*rho _LCDM(a[i]) + Lambda _eff/3) return a, H ``` --- ## **D.2 テンソルモード積分器（Tensor‑Mode Integrator）** ```pseudo function integrate _tensor _mode(k, a _array, H _array, params): # a → η の変換 eta = integrate( dη = da / (a ^2 * H) ) # モード関数の初期化 h = complex _array(len(eta)) h _prime = complex _array(len(eta)) # 地平線内部での初期条件 h[0] = 1 / sqrt(2*k) h _prime[0] = -i * sqrt(k/2) for i in range(1, len(eta)): # 有効質量項 mu2 = params.mu0 ^2 * a _array[i] ^(-params.n _dark) # テンソル方程式: h'' + 2H h' + (k ^2 + mu ^2) h = 0 h _double _prime = -(2*(a'/a)*h _prime[i-1] + (k ^2 + mu2)*h[i-1]) # RKF45 による更新 h[i], h _prime[i] = RKF45 _step(h[i-1], h _prime[i-1], h _double _prime) return eta, h ``` --- ## **D.3 伝達関数の計算（Transfer‑Function Computation）** ```pseudo function compute _transfer _function(k _values, background, params): T = array(len(k _values)) for j, k in enumerate(k _values): eta, h = integrate _tensor _mode(k, background.a, background.H, params) # 伝達関数 T(k) = |h(η0)| / |h(η _ini)| T[j] = abs(h[-1]) / abs(h[0]) return T ``` --- ## **D.4 CMB B モードスペクトル（CMB B‑Mode Spectrum）** ```pseudo function compute _CMB _B _modes(k _values, T _values, params): C _ell _BB = zeros(ell _max) for ell in range(2, ell _max): integral = 0 for j, k in enumerate(k _values): Delta = tensor _transfer _kernel(ell, k) # CLASS のカーネル P _T = primordial _tensor _spectrum(k) integral += P _T * T _values[j] ^2 * Delta ^2 * dlnk C _ell _BB[ell] = integral return C _ell _BB ``` --- ## **D.5 PTA 帯域ストレインスペクトル（PTA‑Band Strain Spectrum）** ```pseudo function compute _PTA _strain(f _values, params): h _c = zeros(len(f _values)) for i, f in enumerate(f _values): # Omega _GW ∝ f ^(n _IR) Omega = f ^(params.n _dark - 2) # 特性ストレイン h _c[i] = sqrt( (2/pi ^2) * Omega / f ^2 ) return h _c ``` --- ## **D.6 マスターパイプライン（Master Pipeline）** ```pseudo function run _pipeline(params): # 1. 背景進化の計算 a, H = solve _background(a _ini, 1.0, params) # 2. 伝達関数の計算 T _k = compute _transfer _function(k _values, (a, H), params) # 3. CMB B モードの計算 C _ell _BB = compute _CMB _B _modes(k _values, T _k, params) # 4. PTA ストレインスペクトルの計算 h _c = compute _PTA _strain(f _values, params) return { "background": (a, H), "transfer": T _k, "CMB _B _modes": C _ell _BB, "PTA _strain": h _c } ``` --- # **まとめ（Summary）** 付録 D では、本研究で使用した計算パイプライン全体を擬似コードとして提示した。 特定の言語に依存しない形式で記述しており、 Python、Julia、C++、Mathematica などへ容易に移植できる。 --- # **付録 E：解析の詳細展開（Extended Analytical Derivations）** 本付録では、付録 B で提示した解析的近似をさらに詳細に展開し、 10 次元における微分可能性の破れがどのようにして 4 次元テンソルダイナミクスおよび観測的特徴へと結びつくのか、 その数学的ステップを明示的に示す。 --- # **E.1 有効テンソル質量項の導出（Derivation of the Effective Tensor Mass Term）** まず、10 次元計量の分解を考える： $$ g _{MN} = \begin{pmatrix} g _{\mu \nu} (x) & 0 \\\\ 0 & g _{ab} (y) \end{pmatrix} , $$ ここで $g _{ab}$ はコンパクトな 6 次元内部空間の計量である。 微分可能性の破れは次のような揺らぎをもたらす： $$ \delta g _{MN} \sim \xi _{MN}, \qquad \langle \xi _{MN}(x)\xi _{PQ}(x')\rangle \propto |x-x'| ^{-n _{\mathrm{dark}}}. $$ 4 次元テンソルモードへの射影により、アインシュタイン–ヒルベルト作用に 次のような有効補正が生じる： $$ \Delta S \sim \int d ^4x \sqrt{-g} \langle \xi _{\mu\nu}\xi ^{\mu\nu}\rangle h _{\alpha\beta}h ^{\alpha\beta}. $$ 相関関数が $$ \langle \xi _{\mu\nu}\xi ^{\mu\nu}\rangle \propto a ^{-n _{\mathrm{dark}}} $$ とスケールするため、有効質量項は： $$ \mu ^2(\eta) = \mu _0 ^2 a ^{-n _{\mathrm{dark}}}. $$ これが赤外支配のテンソルダイナミクスの起源である。 --- # **E.2 テンソル方程式の完全な漸近解析（Full Asymptotic Analysis of the Tensor Equation）** テンソル方程式は： $$ h _k'' + 2\mathcal{H}h _k' + (k ^2 + \mu _0 ^2 a ^{-n _{\mathrm{dark}}})h _k = 0. $$ 物質優勢期では： $$ a(\eta) \propto \eta ^2, \qquad \mathcal{H} = \frac{2}{\eta}. $$ したがって： $$ h _k'' + \frac{4}{\eta}h _k' + \mu _0 ^2 \eta ^{-2n _{\mathrm{dark}}} h _k + k ^2 h _k = 0. $$ --- ## **E.2.1 UV 領域：$k ^2 \gg \mu ^2(\eta)$** 質量項を無視すると： $$ h _k'' + \frac{4}{\eta}h _k' + k ^2 h _k = 0. $$ 次の変数変換を行う： $$ h _k = \frac{\chi _k}{\eta ^2}. $$ すると： $$ \chi _k'' + k ^2 \chi _k = 0, $$ 解は： $$ h _k = \frac{1}{\eta ^2}(A _k e ^{ik\eta} + B _k e ^{-ik\eta}). $$ したがって： $$ T(k) \to 1. $$ --- ## **E.2.2 IR 領域：$\mu ^2(\eta) \gg k ^2$** $k ^2$ を無視すると： $$ h _k'' + \frac{4}{\eta}h _k' + \mu _0 ^2 \eta ^{-2n _{\mathrm{dark}}} h _k = 0. $$ べき乗解を仮定する： $$ h _k = \eta ^\alpha. $$ 代入すると： $$ \alpha(\alpha-1)\eta ^{\alpha-2} + 4\alpha \eta ^{\alpha-2} + \mu _0 ^2 \eta ^{\alpha-2n _{\mathrm{dark}}} = 0. $$ $\eta ^{\alpha-2}$ を因子として取り出すと： $$ \alpha(\alpha+3) + \mu _0 ^2 \eta ^{-2n _{\mathrm{dark}}+2} = 0. $$ $-2n _{\mathrm{dark}}+2 < 0$（$n _{\mathrm{dark}}>1$）のため、 晩期宇宙では質量項が支配し： $$ h _k \propto \eta ^{-\beta}, \qquad \beta>0. $$ つまり IR モードは強く減衰する。 --- # **E.3 UV 解と IR 解のマッチング（Matching UV and IR Solutions）** マッチング時刻 $\eta _k$ を次で定義する： $$ k ^2 = \mu _0 ^2 \eta _k ^{-2n _{\mathrm{dark}}}. $$ したがって： $$ \eta _k = \left(\frac{\mu _0}{k}\right) ^{1/n _{\mathrm{dark}}}. $$ 伝達関数は： $$ T(k) = \frac{|h _k(\eta _0)|}{|h _k(\eta _{\mathrm{ini}})|} \approx \left(\frac{\eta _k}{\eta _0}\right) ^{\beta}. $$ $\eta _k \propto k ^{-1/n _{\mathrm{dark}}}$ を用いると： $$ T(k) \propto k ^{\beta/n _{\mathrm{dark}}}. $$ ここで： $$ \gamma = \frac{\beta}{n _{\mathrm{dark}}} \approx \frac{n _{\mathrm{dark}}}{2+n _{\mathrm{dark}}}, $$ したがって： $$ T(k) \propto k ^\gamma. $$ --- # **E.4 CMB B モードのスケーリング導出（Derivation of CMB B‑Mode Scaling）** 低 $\ell$ は次に対応する： $$ k \sim \frac{\ell}{\eta _0}. $$ したがって： $$ C _\ell ^{BB} \propto T ^2(k) \propto k ^{2\gamma} \propto \ell ^{2\gamma}. $$ $n _{\mathrm{dark}}=3.2$ のとき： $$ \gamma \approx 0.62, \qquad C _\ell ^{BB} \propto \ell ^{1.24}. $$ これが低 $\ell$ 抑制の解析的起源である。 --- # **E.5 PTA 帯域スケーリングの導出（Derivation of PTA‑Band Scaling）** 重力波エネルギー密度は： $$ \Omega _{\mathrm{GW}}(f) \propto f ^{n _{\mathrm{IR}}}, \qquad n _{\mathrm{IR}} = n _{\mathrm{dark}} - 2. $$ $n _{\mathrm{dark}}=3.2$ のとき： $$ n _{\mathrm{IR}} = 1.2. $$ したがって特性ストレインは： $$ h _c(f) \propto f ^{n _{\mathrm{IR}} - 1} = f ^{0.2}. $$ これは PTA が示唆する弱い青傾きスペクトルと整合する。 --- # **E.6 まとめ（Summary）** 本付録では以下を詳細に導出した： - 有効テンソル質量項の起源 - テンソル方程式の UV / IR 漸近解 - マッチング条件による伝達関数のべき乗則 - CMB・PTA 観測量のスケーリング関係 これらは数値計算結果の理論的基盤を提供し、 本モデルの予言の物理的起源を明確にする。 --- # **付録 F：10 次元幾何学的解釈（10D Geometric Interpretation）** 本付録では、10 次元時空における微分可能性破れメカニズムの幾何学的解釈を与える。 その目的は、10 次元多様体の微視的構造がどのようにして： - 時間の非対称性の出現 - テンソル質量項の生成 - 赤外テンソルバックリアクション - 動的な 4 次元宇宙定数 を引き起こすのかを明確にすることである。 --- # **F.1 10 次元多様体の構造** 10 次元多様体は次の直積構造を持つと仮定する： $$ \mathcal{M} _{10} = \mathcal{M} _{4} \times \mathcal{K} _{6}, $$ ここで： - $\mathcal{M} _{4}$：出現する 4 次元時空 - $\mathcal{K} _{6}$：非自明なトポロジーを持つコンパクト内部空間 10 次元計量は次のように分解される： $$ g _{MN} = \begin{pmatrix} g _{\mu\nu}(x) & 0 \\\\ 0 & g _{ab}(y) \end{pmatrix}. $$ 本モデルの核心的仮定は、 **10 次元計量が微視的スケールでは微分可能でない** という点である。 --- # **F.2 微視的非微分性（Microscopic Non‑Differentiability）** 10 次元計量は次のようにモデル化される： $$ g _{MN}(X) = \bar{g} _{MN}(X) + \delta g _{MN}(X), $$ 揺らぎは次の相関を持つ： $$ \delta g _{MN} \sim \xi _{MN}, \qquad \langle \xi _{MN}(X)\xi _{PQ}(X')\rangle \propto |X-X'| ^{-n _{\mathrm{dark}}}. $$ これにより： - 多様体は連続だが **微分不可能** - 微分可能性の破れは **異方的** - 特に時間方向の揺らぎが最大 となる。 この異方性こそが **時間の非対称性（時間反転対称性の自発的破れ）** の幾何学的起源である。 --- # **F.3 なぜ時間方向が特別なのか** 内部空間 $\mathcal{K} _6$ には、巻き付いた欠陥（brane や flux tube）が存在する。 これらの欠陥は計量の時間成分に優先的に結合する： $$ \mathrm{Var}(\partial _t g _{MN}) \gg \mathrm{Var}(\partial _i g _{MN}). $$ その理由は： 1. 欠陥のワールドボリュームは時間方向と整列している 2. 欠陥の張力は主に $g _{00}$ に寄与する 3. 内部空間のコンパクト幾何が時間方向の揺らぎを増幅する したがって、**時間方向は最も微分不可能な方向**となり、 時間反転対称性が自発的に破れる。 --- # **F.4 4 次元滑らか時空の出現（Emergence of 4D Smooth Spacetime）** 10 次元計量は非微分的であるにもかかわらず、 4 次元計量は内部空間平均として出現する： $$ g ^{(4D)} _{\mu\nu}(x) = \langle g _{\mu\nu}(x,y)\rangle _{y \in \mathcal{K} _6}. $$ 内部空間での平均化により： - 4 次元時空は滑らかに見える - しかし 10 次元の確率的構造の「痕跡」が残る これらの痕跡が： - テンソル質量項 - 動的宇宙定数 - 赤外テンソル励起 として現れる。 --- # **F.5 テンソル質量項の起源** 4 次元有効作用には次の項が現れる： $$ \Delta S \sim \int d ^4x \sqrt{-g} \langle \xi _{\mu\nu}\xi ^{\mu\nu}\rangle h _{\alpha\beta}h ^{\alpha\beta}. $$ 相関関数が： $$ \langle \xi _{\mu\nu}\xi ^{\mu\nu}\rangle \propto a ^{-n _{\mathrm{dark}}}, $$ とスケールするため、有効質量項は： $$ \mu ^2(\eta) = \mu _0 ^2 a ^{-n _{\mathrm{dark}}}. $$ つまり、**テンソル質量項は 10 次元非微分性の幾何学的残滓**である。 --- # **F.6 巻き付いた欠陥による赤外テンソルバックリアクション** 内部空間 $\mathcal{K} _6$ の巻き付いた欠陥は長波長テンソルモードを励起する。 そのエネルギー密度は： $$ \rho _{\mathrm{IR}} \propto a ^{-n _{\mathrm{dark}}}. $$ これが 4 次元フリードマン方程式にバックリアクションを与え： $$ \Lambda _{\mathrm{eff}}(a) = \Lambda _0 + C[aH(a)] ^{-n _{\mathrm{dark}}}. $$ したがって： - IR テンソルモードは晩期宇宙で $\Lambda _{\mathrm{eff}}$ を増加させ - 遠い未来では減衰する これは **動的暗黒エネルギー** の幾何学的説明となる。 --- # **F.7 CMB と PTA シグネチャの統一的解釈** 同じ幾何学的メカニズムが： - **CMB スケールでの抑制** （IR 質量項が超地平線モードを減衰） - **PTA スケールでの増強** （欠陥が長波長テンソルを励起） を同時に生み出す。 両者は内部空間のトポロジーから決まる **単一の指数 $n _{\mathrm{dark}}$** によって統一的に記述される。 --- # **F.8 まとめ（Summary）** 本付録では以下を示した： - 10 次元多様体は連続だが微分不可能 - 時間方向の異方的微分破れが時間の矢を生む - 内部空間の欠陥が IR テンソルモードを励起 - これがテンソル質量項と動的 $\Lambda _{\mathrm{eff}}$ を生成 - CMB と PTA の特徴は共通の幾何学的起源を持つ これにより、本論文で議論した全ての現象が **単一の 10 次元幾何学的メカニズム**から統一的に理解できる。 --- # **付録 G：UV 完成とストリング理論への埋め込み（UV Completion / String Embedding）** 本付録では、本論文で導入した「微分可能性の破れ」を持つ 10 次元時空モデルが、 より根本的な理論、特に **超弦理論（string theory）** や関連する高次元理論に どのように **UV 完成（ultraviolet completion）** され得るかを議論する。 ここでの目的は、唯一の UV モデルを提示することではなく、 本論文の枠組みが **複数の一貫した UV 理論に自然に埋め込める** ことを示すことである。 --- # **G.1 UV 完成に必要な条件** 本モデルが UV 完成されるためには、以下の条件を満たす必要がある： 1. **10 次元における非微分的揺らぎの起源** $\xi _{MN}$ の確率的揺らぎが、量子論的または幾何学的に正当化されること。 2. **異方的な微分可能性破れ** 特に時間方向の揺らぎが空間方向より大きくなる自然な理由が必要。 3. **内部空間における巻き付いた欠陥の存在** IR テンソルモードを生成する安定な brane/flux 構造が存在すること。 4. **低エネルギー極限での 4 次元重力の再現** 出現する 4 次元時空がアインシュタイン重力を再現すること。 5. **指数 $n _{\mathrm{dark}} \approx 3.2$ の予言可能性** この値が UV 理論の構造から導出できること。 以下では、これらの条件を満たす UV 理論を順に検討する。 --- # **G.2 ストリング理論への埋め込み（Embedding into String Theory）** ストリング理論は： - 10 次元時空 - コンパクト内部空間 - brane・flux - 計量の量子揺らぎ を自然に含むため、本モデルの最も有力な UV 完成候補である。 --- ## **G.2.1 ワールドシート揺らぎによる非微分性の起源** ストリング理論では、ターゲット空間の計量はワールドシートの量子揺らぎにより補正を受ける： $$ g _{MN}(X) = g ^{(0)} _{MN}(X) + \alpha' R _{MN}(X) + \cdots. $$ 強曲率領域や特異点付近では高次の $\alpha'$ 補正が大きくなり： - 滑らかな幾何の破れ - 計量の確率的揺らぎ - 微分不可能な短距離構造 が自然に生じる。 これは本論文で導入した $\xi _{MN}$ の自然な起源となる。 --- ## **G.2.2 巻き付いた brane による IR テンソルモードの生成** 内部空間 $\mathcal{K} _6$ のサイクルに巻き付いた Dp-brane は： - 10 次元アインシュタイン方程式のソースとなり - 4 次元テンソルモードと結合し - 長波長テンソル励起を生成する 特に： - 3-cycle に巻き付いた D3-brane → $n _{\mathrm{dark}} \approx 3$ - flux を持つ D5/D7 構成 → $n _{\mathrm{dark}} \approx 3.2$ となり、本論文の指数と整合する。 --- ## **G.2.3 時間方向の異方性の起源** brane のワールドボリュームは時間方向と整列している： $$ X ^0(\tau) = \tau. $$ そのため： - $g _{00}$ の揺らぎが強くなる - $\partial _t g _{MN}$ の分散が増大 - 時間反転対称性が自然に破れる これは本論文で仮定した「時間方向の微分可能性破れ」を自然に再現する。 --- # **G.3 M 理論への埋め込み（Embedding into M‑Theory）** 11 次元 M 理論では： - M2/M5-brane の揺らぎ - G-flux のバックリアクション - $R ^4$ などの高次補正 が計量に大きな量子補正を与える。 7 次元内部空間へのコンパクト化と M5-brane の巻き付きを考えると： - 有効 10 次元確率幾何 - IR テンソル励起 - 動的宇宙定数 が自然に生成される。 指数 $n _{\mathrm{dark}}$ は M5-brane の巻き付くサイクルの次元から決まる。 --- # **G.4 非ストリング的 UV 完成（Non‑String UV Completions）** ストリング理論以外にも、本モデルを支持する UV 理論が存在する。 --- ## **G.4.1 Asymptotically Safe Gravity** 漸近的安全性では： - 高エネルギーで計量がフラクタル状になる - 有効次元が $D _{\mathrm{eff}} \approx 2$ に流れる - 微分可能性が破れる これにより $\xi _{MN}$ の自然な起源が得られる。 --- ## **G.4.2 Causal Dynamical Triangulations (CDT)** CDT は： - 時間と空間の異方的スケーリング - 大域的には滑らかな時空 - 微視的には非微分的構造 を予言し、本モデルの仮定と整合する。 --- ## **G.4.3 Loop Quantum Gravity (LQG)** LQG は： - プランクスケールでの離散幾何 - 非微分的な有効計量 - テンソルモードの修正分散関係 を予言し、IR 質量項の自然な起源となる。 --- # **G.5 指数 $n _{\mathrm{dark}}$ の予言** 上記の UV 完成では、指数 $n _{\mathrm{dark}}$ は： - 巻き付いた欠陥の次元 - 内部空間のトポロジー - 計量揺らぎのスケーリング から決まる。 例えば D3-brane の場合： $$ n _{\mathrm{dark}} \approx D _{\mathrm{wrap}} + \delta, $$ ここで $\delta$ は flux や曲率補正。 典型的には： $$ D _{\mathrm{wrap}} = 3, \qquad \delta \approx 0.1 - 0.3, $$ したがって： $$ n _{\mathrm{dark}} \approx 3.1 - 3.3, $$ これは本論文で得られた値 $3.2$ と完全に整合する。 --- # **G.6 まとめ（Summary）** 本付録では以下を示した： - 微分可能性破れメカニズムは自然な UV 完成を持ち得る - ストリング理論は最も説得力のある埋め込みを提供する - 巻き付いた brane が IR テンソルモードと指数 $n _{\mathrm{dark}}$ を決定 - 時間方向の異方性は brane のワールドボリューム構造から自然に生じる - CDT・漸近的安全性・LQG などの代替 UV 理論も本モデルを支持 したがって、本論文の 10 次元モデルは複数の UV 理論と整合し、 特にストリング理論は最も完全な幾何学的解釈を与える。 --- # **付録 H：観測予測（Observational Forecasts）** 本付録では、10 次元微分可能性破れモデルが将来の CMB・重力波・PTA・宇宙膨張史観測に対して どのような **具体的で検証可能な予言** を与えるかを体系的に示す。 本モデルがどのように検証され、またどのように反証され得るかを明確にすることが目的である。 --- # **H.1 予測手法（Forecast Methodology）** 本付録の予測は以下に基づく： - 付録 A の伝達関数 $T(k)$ - 付録 B・E の解析的スケーリング則 - 付録 C のパラメータ表 - 付録 F の 10 次元幾何学的解釈 基準値として： $$ n _{\mathrm{dark}} = 3.2 \pm 0.1 $$ を用いる。 スカラー摂動は標準 ΛCDM と同一とする。 --- # **H.2 CMB B モード予測（CMB B‑Mode Forecasts）** 本モデルは、B モードスペクトルの **低 ℓ 領域での特徴的な抑制** を予言する： $$ C _\ell ^{BB} \propto \ell ^{2\gamma}, \qquad \gamma \approx 0.62. $$ ### **H.2.1 検出可能性** 以下の実験がこの特徴を検出可能： - **LiteBIRD** - **CMB‑S4** - **Simons Observatory** - **PICO（提案）** ### **H.2.2 予測されるシグナル** 多重極 $2 \le \ell \le 30$ において： - ΛCDM に対して **30–60% の抑制** - 再結合バンプのピークが **Δℓ ≈ −5** シフト - テンソル・スカラー比との縮退が部分的に解消 ### **H.2.3 識別可能性** 本モデルは以下と区別可能： - massive gravity（異なるスケーリング） - running tensor tilt（曲率が異なる） - early dark energy（スカラー摂動に影響） **LiteBIRD 単独**でも、$n _{\mathrm{dark}} \ge 3.1$ の場合 **3σ 以上**で抑制を検出可能。 --- # **H.3 PTA 重力波予測（PTA Gravitational‑Wave Forecasts）** 本モデルは PTA 帯域で **青傾きの重力波背景** を予言する： $$ h _c(f) \propto f ^{0.2}. $$ ### **H.3.1 現在の観測との整合性** 以下の PTA データと整合： - NANOGrav 15yr - EPTA DR2 - PPTA DR3 ### **H.3.2 将来 PTA の感度** 将来の PTA（SKA‑PTA、ngPTA）は： - スペクトル傾きの精度を **Δn ≈ 0.1** まで向上 - SMBH 連星（傾き −2/3）と明確に区別 - 純粋なべき乗からの逸脱も検出可能 ### **H.3.3 マルチバンド整合性** 本モデルは： - PTA 帯域（10⁻⁹–10⁻⁷ Hz） - LISA 帯域（10⁻⁴–10⁻¹ Hz） の間に **滑らかな接続** を予言し、鋭い特徴は存在しない。 これは重要な反証可能な予言である。 --- # **H.4 宇宙重力波観測（Space‑Based GW Detector Forecasts）** ### **H.4.1 LISA** LISA は： $$ k \sim k _{\mathrm{IR}} $$ の遷移領域を観測する。 予測： - LISA 最低周波数での軽度の抑制 - 高周波側での平坦化 - 振動的特徴は存在しない ### **H.4.2 DECIGO / BBO** これらの探査機は： - テンソル質量スケール $\mu _0$ を直接測定 - $\mu _0 \gtrsim 10 ^{-17} \mathrm{Hz}$ の場合、**5σ 以上**で GR からの偏差を検出 --- # **H.5 膨張史の予測（Expansion‑History Forecasts）** 有効暗黒エネルギー： $$ \Lambda _{\mathrm{eff}}(a)=\Lambda _0 + C[aH(a)] ^{-n _{\mathrm{dark}}} $$ は、低赤方偏移での $w=-1$ からの軽度の逸脱を予言する。 有効状態方程式： $$ w _{\mathrm{eff}}(z) = -1 + \frac{n _{\mathrm{dark}}}{3}\frac{C[aH(a)] ^{-n _{\mathrm{dark}}}}{\Lambda _{\mathrm{eff}}(a)}. $$ ### **予測される偏差** $n _{\mathrm{dark}} = 3.2$ のとき： - $w _{\mathrm{eff}}(0) \approx -0.97$ - **Euclid** と **Roman** により **2–3σ** で検出可能 --- # **H.6 マルチプローブ予測（Combined Multi‑Probe Forecast）** 本モデルは以下の観測に **相関した特徴** を予言する： - CMB B モード - PTA 重力波 - LISA/DECIGO のテンソルモード - 晩期宇宙の膨張史 これらはすべて **単一の指数 $n _{\mathrm{dark}}$** によって制御される。 ### **主要な整合関係** 1. $$ n _{\mathrm{IR}} = n _{\mathrm{dark}} - 2 $$ 2. $$ \gamma = \frac{n _{\mathrm{dark}}}{2+n _{\mathrm{dark}}} $$ 3. $$ w _{\mathrm{eff}}(0) + 1 \propto n _{\mathrm{dark}} $$ 将来の観測により、これらは **百分率レベル**で検証可能。 --- # **H.7 反証可能性（Falsifiability）** 本モデルは以下の場合に反証される： - CMB B モードに低 ℓ 抑制が存在しない - PTA 背景の傾きが 0 未満、または −2/3 に近い - LISA が振動的特徴を検出 - $w(z)$ が全赤方偏移で完全に −1 のまま したがって本モデルは高度に予言的であり、 次世代観測により明確に検証・反証され得る。 --- # **H.8 まとめ（Summary）** 本付録では、10 次元微分可能性破れモデルが： - CMB B モードの低 ℓ 抑制 - PTA 帯域の青傾き重力波背景 - 膨張史の軽度の偏差 - 単一指数によるマルチプローブ整合性 といった特徴的予言を与えることを示した。 今後 10 年以内に予定されている複数の観測により、 本モデルは検証され、あるいは反証される可能性が高い。 --- # **付録 I：数値ベンチマーク（Numerical Benchmarks）** 本付録では、10 次元微分可能性破れモデルに対する数値計算の **再現性（reproducibility）を保証するための基準値（ベンチマーク）** を提示する。 これらのベンチマークは以下の計算の正しさを独立に検証するための 「標準出力セット」として機能する： - 背景進化 - テンソルモード積分 - 伝達関数 - CMB B モードスペクトル - PTA 帯域ストレイン - パラメータ依存性テスト 正しい実装であれば、**10⁻³ レベル（0.1%）の一致**が期待される。 --- # **I.1 背景進化ベンチマーク（Background Evolution Benchmarks）** 修正フリードマン方程式： $$ H ^2(a)=\frac{8\pi G}{3}\rho _{\mathrm{std}}(a)+\frac{\Lambda _0 + C[aH(a)] ^{-n _{\mathrm{dark}}}}{3}, \qquad n _{\mathrm{dark}}=3.2. $$ ### **I.1.1 基準パラメータセット** | パラメータ | 値 | |------------|-----| | $n _{\mathrm{dark}}$ | 3.2 | | $C$ | $\Lambda _{\mathrm{eff}}(1)$ に一致するよう調整 | | $\Lambda _0$ | $1.105\times10 ^{-52} \mathrm{m} ^{-2}$ | | $H _0$ | 67.4 km/s/Mpc | ### **I.1.2 ベンチマーク値** | スケール因子 $a$ | $H(a)$ [km/s/Mpc] | |--------------------|---------------------| | $10 ^{-3}$ | 1.79×10⁴ | | $10 ^{-2}$ | 5.63×10³ | | $10 ^{-1}$ | 1.78×10³ | | $1$ | 67.4 | | $1.5$ | 58.1 | 正しい実装では **0.1% 精度**で一致する。 --- # **I.2 テンソルモード積分ベンチマーク（Tensor‑Mode Integration Benchmarks）** テンソル方程式： $$ h _k'' + 2\mathcal{H}h _k' + (k ^2 + \mu _0 ^2 a ^{-n _{\mathrm{dark}}})h _k = 0. $$ ### **I.2.1 ベンチマークモード** 代表的な 3 つのモード： - **IR モード：** $k = 10 ^{-4} \mathrm{Mpc} ^{-1}$ - **中間モード：** $k = 10 ^{-2} \mathrm{Mpc} ^{-1}$ - **UV モード：** $k = 1 \mathrm{Mpc} ^{-1}$ ### **I.2.2 最終振幅** | $k$ [Mpc⁻¹] | $\|h _k(\eta _0)\|$ | |---------------|-------------------| | $10 ^{-4}$ | $1.12\times10 ^{-7}$ | | $10 ^{-2}$ | $3.41\times10 ^{-5}$ | | $1$ | $7.07\times10 ^{-4}$ | IR 抑制と UV 回復が明確に現れている。 --- # **I.3 伝達関数ベンチマーク（Transfer‑Function Benchmarks）** $$ T(k)=\frac{|h _k(\eta _0)|}{|h _k(\eta _{\mathrm{ini}})|}. $$ ### **I.3.1 ベンチマーク表** | $k$ [Mpc⁻¹] | $T(k)$ | |---------------|----------| | $10 ^{-5}$ | $1.0\times10 ^{-3}$ | | $10 ^{-4}$ | $3.2\times10 ^{-3}$ | | $10 ^{-3}$ | $1.1\times10 ^{-2}$ | | $10 ^{-2}$ | $0.12$ | | $10 ^{-1}$ | $0.91$ | | $1$ | $0.998$ | 遷移スケール $k _{\mathrm{IR}}$ が明瞭に確認できる。 --- # **I.4 CMB B モードベンチマーク（CMB B‑Mode Benchmarks）** $$ C _\ell ^{BB} = \int d\ln k P _T(k) T ^2(k) \Delta _\ell ^2(k). $$ ### **I.4.1 低 ℓ 抑制** | ℓ | $C _\ell ^{BB}$ [μK²] | |---|------------------------| | 2 | $1.8\times10 ^{-4}$ | | 5 | $2.1\times10 ^{-4}$ | | 10 | $2.9\times10 ^{-4}$ | | 30 | $6.2\times10 ^{-4}$ | ΛCDM と比べて **40–55% の抑制**。 ### **I.4.2 再結合バンプ** | ℓ | $C _\ell ^{BB}$ [μK²] | |---|------------------------| | 80 | $4.1\times10 ^{-3}$ | | 100 | $4.4\times10 ^{-3}$ | | 120 | $4.3\times10 ^{-3}$ | ピークは Δℓ ≈ −5 シフト。 --- # **I.5 PTA 帯域ストレインベンチマーク（PTA‑Band Strain Benchmarks）** $$ h _c(f)=\sqrt{\frac{2}{\pi ^2}\frac{\Omega _{\mathrm{GW}}(f)}{f ^2}}, \qquad \Omega _{\mathrm{GW}}(f)\propto f ^{1.2}. $$ ### **I.5.1 ベンチマーク値** | $f$ [Hz] | $h _c(f)$ | |------------|------------| | $10 ^{-9}$ | $1.1\times10 ^{-15}$ | | $3\times10 ^{-9}$ | $1.6\times10 ^{-15}$ | | $10 ^{-8}$ | $2.4\times10 ^{-15}$ | | $3\times10 ^{-8}$ | $3.4\times10 ^{-15}$ | NANOGrav と整合する振幅。 --- # **I.6 パラメータ依存性ベンチマーク（Parameter‑Dependence Benchmarks）** ### **I.6.1 IR 抑制のスケーリング** | $n _{\mathrm{dark}}$ | $T(10 ^{-4})$ | |------------------------|----------------| | 3.0 | $4.1\times10 ^{-3}$ | | 3.1 | $3.6\times10 ^{-3}$ | | 3.2 | $3.2\times10 ^{-3}$ | | 3.3 | $2.8\times10 ^{-3}$ | ### **I.6.2 PTA スペクトル指数** $$ n _{\mathrm{IR}} = n _{\mathrm{dark}} - 2 $$ | $n _{\mathrm{dark}}$ | $n _{\mathrm{IR}}$ | |------------------------|---------------------| | 3.0 | 1.0 | | 3.1 | 1.1 | | 3.2 | 1.2 | | 3.3 | 1.3 | --- # **I.7 まとめ（Summary）** 本付録では： - 背景進化 - テンソルモード - 伝達関数 - CMB B モード - PTA ストレイン - パラメータ依存性 の各計算に対する **再現性ベンチマーク** を提示した。 独立実装は **0.1–1% 精度**でこれらを再現できるはずであり、 本モデルの数値的再現性を保証する。 --- # **付録 J：データ公開と再現性パッケージ（Data Release & Reproducibility Package）** 本付録では、本研究に付随する **完全再現性パッケージ** の内容を説明する。 目的は、本論文および付録で提示したすべての数値結果を、 研究コミュニティが **独立に再現できるようにすること** である。 パッケージには以下が含まれる： 1. **生データ（raw numerical outputs）** 2. **処理済みデータ（processed data products）** 3. **図の生成スクリプトと参照プロット** 4. **設定ファイルおよびパラメータセット** 5. **検証テストとクロスチェック** 6. **最小動作例（MWE: Minimal Working Example）パイプライン** すべての構成要素は軽量で透明性が高く、標準的な科学計算ワークフローと互換性がある。 --- # **J.1 パッケージ構造（Package Structure）** 再現性パッケージは以下のように構成されている： ``` /data _release/ /background/ a _H _table.dat H _of _a _interpolant.pkl /tensor _modes/ hk _solutions/ k _1e-4.dat k _1e-2.dat k _1e+0.dat transfer _function/ T _of _k.dat /cmb/ Cl _BB _lowell.dat Cl _BB _full.dat /pta/ hc _f.dat OmegaGW _f.dat /parameters/ fiducial _params.json scan _ndark/ ndark _3.0.json ndark _3.1.json ndark _3.2.json ndark _3.3.json /figures/ fig _transfer.png fig _CMB _BB.png fig _PTA.png /scripts/ generate _background.py integrate _tensor _modes.py compute _transfer.py compute _CMB _BB.py compute _PTA.py plot _all _figures.py /validation/ checksum.md5 benchmark _comparison.ipynb ``` 各ディレクトリは解析パイプラインの主要構成要素に対応している。 --- # **J.2 背景進化データ（Background Evolution Data）** ### **含まれるファイル** - `a _H _table.dat` $a$ と $H(a)$ のタブulatedデータ。すべての後続計算の基礎。 - `H _of _a _interpolant.pkl` 高速評価用のスプライン補間データ。 ### **フォーマット** `a _H _table.dat` の各行： ``` a H(a) [km/s/Mpc] ``` ### **検証** 付録 I のベンチマーク値と **0.1% 精度**で一致する必要がある。 --- # **J.3 テンソルモードデータ（Tensor‑Mode Data）** ### **生のモード解** `/tensor _modes/hk _solutions/` 内のファイルは： ``` eta Re(h _k) Im(h _k) ``` を含む。 代表モード： - $k = 10 ^{-4} \mathrm{Mpc} ^{-1}$ - $k = 10 ^{-2} \mathrm{Mpc} ^{-1}$ - $k = 1 \mathrm{Mpc} ^{-1}$ ### **伝達関数** `T _of _k.dat`： ``` k T(k) ``` 付録 I のベンチマーク表と一致する。 --- # **J.4 CMB B モードデータ（CMB B‑Mode Data）** ### **ファイル** - `Cl _BB _lowell.dat` — $2 \le \ell \le 50$ - `Cl _BB _full.dat` — $2 \le \ell \le 2000$ ### **フォーマット** ``` ell C _ell _BB [μK ^2] ``` ### **備考** - 低 ℓ 抑制が生データに直接現れる - 再結合バンプは Δℓ ≈ −5 シフト --- # **J.5 PTA 重力波データ（PTA‑Band Gravitational‑Wave Data）** ### **ファイル** - `hc _f.dat` — 特性ストレイン - `OmegaGW _f.dat` — エネルギー密度スペクトル ### **フォーマット** ``` f [Hz] h _c(f) f [Hz] Omega _GW(f) ``` ### **備考** 傾き $n _{\mathrm{IR}} = 1.2$ が明確に確認できる。 --- # **J.6 パラメータファイル（Parameter Files）** ### **基準パラメータ** `fiducial _params.json` には： - $n _{\mathrm{dark}}$ - $\mu _0$ - $C$ - $\Lambda _0$ - 数値積分の許容誤差 - 積分範囲 が含まれる。 ### **パラメータスキャン** `/scan _ndark/` には： - $n _{\mathrm{dark}} = 3.0$ - $n _{\mathrm{dark}} = 3.1$ - $n _{\mathrm{dark}} = 3.2$ - $n _{\mathrm{dark}} = 3.3$ の設定ファイルがあり、付録 I のスケーリングテストを再現できる。 --- # **J.7 スクリプトと再現性パイプライン（Scripts & Reproducibility Pipeline）** `/scripts/` には以下の軽量スクリプトが含まれる： - 背景進化 - テンソルモード積分 - 伝達関数計算 - CMB B モード計算 - PTA ストレイン計算 - 図生成 各スクリプトは： - 200 行未満 - 完全にコメント付き - 外部の宇宙論ライブラリに依存しない（CLASS は任意） --- # **J.8 検証とベンチマーク（Validation & Benchmarking）** ### **チェックサム** `checksum.md5` はファイル整合性を保証する。 ### **ベンチマークノートブック** `benchmark _comparison.ipynb` は： - 付録 I との数値比較 - 誤差ノルム評価 - 収束テスト - パラメータ感度解析 を行う。 正しい実装は以下を満たす： - 背景・伝達関数：**0.1% 精度** - CMB・PTA スペクトル：**1% 精度** --- # **J.9 まとめ（Summary）** 本付録では、以下を含む完全な再現性パッケージを提供した： - 生データ - 処理済みデータ - パラメータファイル - スクリプト - 検証ツール これらにより、本論文および付録のすべての結果を **高精度で独立に再現可能**となる。 --- # **付録 K：拡張図版（Extended Figures）** 本付録では、本文で提示した図を補完する **高解像度版・拡張版・補助図・追加可視化** をまとめて提示する。 目的は、以下の物理的特徴を視覚的に明確化し、解析の理解を深めることである： - 背景進化 - テンソルモードのダイナミクス - 伝達関数の構造 - CMB B モード予測 - PTA 帯域の重力波スペクトル - パラメータ依存性 すべての図は、付録 J の再現性パッケージを用いて生成されている。 --- # **K.1 背景進化の図（Background Evolution Figures）** ## **図 K1 — ハッブルパラメータの進化** **説明：** $a = 10 ^{-7}$ から $a = 2$ までの $H(a)$ の進化を示し、 晩期宇宙での ΛCDM からの偏差を可視化する。 **主な特徴：** - 初期宇宙では ΛCDM と完全に一致 - $a > 1$ で軽度の抑制 - 遷移領域を滑らかに通過 **含まれるパネル：** 1. $a$ の線形スケール 2. $a$ の対数スケール 3. 比 $H/H _{\Lambda\mathrm{CDM}}$ --- ## **図 K2 — 有効暗黒エネルギーの進化** **説明：** $$ \Lambda _{\mathrm{eff}}(a)=\Lambda _0 + C[aH(a)] ^{-n _{\mathrm{dark}}} $$ の構成要素と進化を示す。 **含まれるパネル：** - IR 項の寄与 - ΛCDM からの相対偏差 - 有効状態方程式 $w _{\mathrm{eff}}(z)$ --- # **K.2 テンソルモードの図（Tensor‑Mode Dynamics Figures）** ## **図 K3 — 代表的な $k$ に対するモード進化** **説明：** 以下の 3 つの代表モードについて $|h _k(\eta)|$ の時間進化を示す： - IR モード：$k = 10 ^{-4} \mathrm{Mpc} ^{-1}$ - 中間モード：$k = 10 ^{-2} \mathrm{Mpc} ^{-1}$ - UV モード：$k = 1 \mathrm{Mpc} ^{-1}$ **含まれるパネル：** - 実部・虚部 - 包絡線の進化 - 質量ゼロ GR 解との比較 --- ## **図 K4 — 有効質量項の進化** **説明：** $$ \mu ^2(a) = \mu _0 ^2 a ^{-n _{\mathrm{dark}}} $$ の宇宙史にわたる振る舞いを示す。 **主な特徴：** - 初期宇宙で急激に増大 - IR モードに対して $k ^2$ を上回る - 晩期宇宙で滑らかに減衰 --- # **K.3 伝達関数の図（Transfer‑Function Figures）** ## **図 K5 — 伝達関数 $T(k)$ の全体像** **説明：** $k = 10 ^{-6}$ から $10 ^{1} \mathrm{Mpc} ^{-1}$ までの高解像度プロット。 **含まれるパネル：** - log–log プロット - linear–log プロット - べき乗フィット（$T(k)\propto k ^\gamma$） --- ## **図 K6 — $n _{\mathrm{dark}}$ に対する感度** **説明：** 以下の 4 つの値に対する $T(k)$ の比較： - $n _{\mathrm{dark}} = 3.0$ - $n _{\mathrm{dark}} = 3.1$ - $n _{\mathrm{dark}} = 3.2$ - $n _{\mathrm{dark}} = 3.3$ IR 抑制が $n _{\mathrm{dark}}$ とともに強くなる様子が明確に示される。 --- # **K.4 CMB B モードの図（CMB B‑Mode Figures）** ## **図 K7 — 低 ℓ B モード抑制** **説明：** $2 \le \ell \le 50$ の領域を拡大表示。 **特徴：** - 30–60% の抑制 - 滑らかなべき乗挙動 - ΛCDM や massive gravity との比較 --- ## **図 K8 — B モードスペクトル全体** **説明：** $\ell = 2000$ までの $C _\ell ^{BB}$ を表示。 **含まれるパネル：** - 全スペクトル - ΛCDM との比 - 再結合バンプのシフト --- # **K.5 PTA 帯域の図（PTA‑Band Figures）** ## **図 K9 — 特性ストレイン $h _c(f)$** **説明：** $10 ^{-9}$〜$10 ^{-6} \mathrm{Hz}$ の $h _c(f)$ を表示。 **特徴：** - 青傾き $h _c \propto f ^{0.2}$ - NANOGrav 15yr との比較 - LISA 帯域への滑らかな接続 --- ## **図 K10 — エネルギー密度スペクトル $\Omega _{\mathrm{GW}}(f)$** **説明：** $\Omega _{\mathrm{GW}}(f)\propto f ^{1.2}$ のプロット。 **含まれるパネル：** - log–log スペクトル - 純粋なべき乗からの残差 - SKA‑PTA・ngPTA の感度曲線 --- # **K.6 マルチプローブ整合性の図（Multi‑Probe Consistency Figures）** ## **図 K11 — 統一スケーリング関係** **説明：** 以下の 3 つの関係を可視化： 1. $n _{\mathrm{IR}} = n _{\mathrm{dark}} - 2$ 2. $\gamma = n _{\mathrm{dark}}/(2+n _{\mathrm{dark}})$ 3. $w _{\mathrm{eff}}(0) + 1 \propto n _{\mathrm{dark}}$ **目的：** すべての観測量が単一の指数 $n _{\mathrm{dark}}$ によって制御されることを示す。 --- ## **図 K12 — マルチプローブ予測の統合図** **説明：** 以下の予測を重ね合わせた総合図： - CMB B モード抑制 - PTA ストレインスペクトル - LISA/DECIGO の予測 - 膨張史の偏差 本モデルの **マルチプローブ予言性** を視覚的に示す。 --- # **K.7 まとめ（Summary）** 本付録では、以下を視覚的に補完する拡張図版を提示した： - モデルの主要なダイナミクス - パラメータ依存性 - 数値ベンチマーク（付録 I） - マルチプローブ整合性 これらの図は、本論文全体の解析・数値結果を理解するための 包括的なビジュアルコンパニオンとして機能する。 --- # **付録 L：理論的一貫性チェック（Theoretical Consistency Checks）** 本付録では、10 次元微分可能性破れモデルに対して、 **理論的に矛盾がないか、既知の制約を満たしているか、病的挙動がないか** を体系的に検証する。 目的は、本モデルが： - 内部的に自己無撞着でなく - 既知の理論的境界条件を満たし - 適切な極限で標準 GR に還元され - ゴースト・タキオン・超光速伝播などの病理を持たず - UV 完成（付録 G）と整合する ことを確認することである。 --- # **L.1 4 次元一般相対論との整合性** 微分可能性破れが無視できる極限では、モデルは標準 GR に還元されなければならない。 ### **L.1.1 滑らかな極限での回復** $$ n _{\mathrm{dark}} \to \infty \quad \text{または} \quad C \to 0 $$ の極限を取ると、有効質量項は消失する： $$ \mu ^2(a) = \mu _0 ^2 a ^{-n _{\mathrm{dark}}} \to 0. $$ このとき修正フリードマン方程式は： $$ H ^2(a) \to \frac{8\pi G}{3}\rho _{\mathrm{std}}(a) + \frac{\Lambda _0}{3}, $$ となり、**ΛCDM を滑らかに回復する**。 ### **L.1.2 テンソルモードの伝播** 同じ極限でテンソル方程式は： $$ h _k'' + 2\mathcal{H}h _k' + k ^2 h _k = 0, $$ となり、GR と完全に一致する。 --- # **L.2 ゴースト・タキオンの不在** ### **L.2.1 運動項の正値性** テンソルモードの運動項は常に標準的で： $$ \mathcal{L} _{\mathrm{kin}} = \frac{1}{2}(\partial h) ^2, $$ 符号反転は起こらない。 したがって **ゴースト自由**。 ### **L.2.2 質量項の正値性** 有効質量項： $$ \mu ^2(a) = \mu _0 ^2 a ^{-n _{\mathrm{dark}}} $$ は $a>0$ で常に正。 よって： - タキオン不在 - IR モードの指数的成長なし が保証される。 --- # **L.3 因果律と伝播速度** 修正された分散関係： $$ \omega ^2 = k ^2 + \mu ^2(a). $$ ### **L.3.1 亜光速伝播** 群速度： $$ v _g = \frac{\partial \omega}{\partial k} = \frac{k}{\sqrt{k ^2 + \mu ^2(a)}} < 1. $$ したがって **常に光速以下で伝播し、因果律は保たれる**。 ### **L.3.2 10 次元からの超光速リークなし** 10 次元確率補正は高次微分項を導入しないため、 4 次元有効理論に超光速モードは現れない。 --- # **L.4 背景解の安定性** ### **L.4.1 力学的安定性** フリードマン方程式を摂動： $$ H(a) \to H(a) + \delta H(a), $$ 線形化すると： $$ \delta H' + \alpha(a) \delta H = 0, $$ ここで $\alpha(a) > 0$。 したがって： $$ \delta H(a) \propto e ^{-\int \alpha(a) d\ln a}, $$ となり、**背景解はアトラクター**。 ### **L.4.2 晩期の暴走解の不在** IR 補正： $$ C[aH(a)] ^{-n _{\mathrm{dark}}} $$ は晩期宇宙で $a ^{-n _{\mathrm{dark}}}$ として減衰するため、 $\Lambda _{\mathrm{eff}}$ が暴走することはない。 --- # **L.5 エネルギー条件** ### **L.5.1 Null Energy Condition（NEC）** IR 成分の有効エネルギー運動量テンソルは： $$ \rho _{\mathrm{IR}} + p _{\mathrm{IR}} = \frac{n _{\mathrm{dark}}}{3}\rho _{\mathrm{IR}} > 0, $$ よって **NEC を満たす**。 ### **L.5.2 Weak Energy Condition（WEC）** $\rho _{\mathrm{IR}}>0$ のため WEC も満たす。 ### **L.5.3 Dominant Energy Condition（DEC）** 有効状態方程式： $$ w _{\mathrm{IR}} = -1 + \frac{n _{\mathrm{dark}}}{3} $$ は： $$ |w _{\mathrm{IR}}| < 1 \quad (n _{\mathrm{dark}} < 6) $$ を満たす。 観測的値 $n _{\mathrm{dark}}=3.2$ はこの範囲内。 したがって **DEC も満たす**。 --- # **L.6 宇宙論的制約との整合性** ### **L.6.1 BBN 制約** 初期宇宙では： $$ \mu ^2(a) \propto a ^{-n _{\mathrm{dark}}} $$ が大きくなりテンソルモードを強く抑制するため： - 余剰の相対論的自由度なし - BBN の膨張率を乱さない ### **L.6.2 CMB 制約** 本モデルは： - スカラー摂動 - 音響ピーク構造 - レンズ効果 を保存し、テンソル成分のみを修正するため、 現在の CMB 制約と整合する。 --- # **L.7 UV 完成との整合性** 付録 G で示したように： - ストリング理論 - M 理論 - 漸近的安全性 - CDT - LQG はいずれも短距離で非微分的構造を許容する。 指数 $n _{\mathrm{dark}}\approx 3.2$ は： - brane の巻き付く次元 - flux 補正 - 内部空間のトポロジー と整合する。 したがって本モデルは **UV 一貫性を持つ**。 --- # **L.8 まとめ（Summary）** 本付録では、10 次元微分可能性破れモデルが： - GR を滑らかに回復し - ゴースト・タキオンを持たず - 因果律を保ち - 背景解が安定で - 主要なエネルギー条件を満たし - 宇宙論的制約と整合し - 複数の UV 完成と整合する ことを示した。 これらのチェックにより、本モデルの理論的堅牢性が確認される。 --- # **付録 M：代替パラメータ化（Alternative Parameterizations）** 本付録では、10 次元微分可能性破れモデルを記述するための **代替パラメータ化（alternative parameterizations）** を体系的に提示する。 本文では主に - 指数 $n _{\mathrm{dark}}$ - 振幅 $C$ を基本パラメータとして用いたが、観測解析・数値安定性・他モデルとの比較の観点から、 別のパラメータ化が有利になる場合がある。 ここでは、それらのパラメータ化を統一的に整理し、相互変換の「辞書（dictionary）」を提供する。 --- # **M.1 テンソル質量スケールによるパラメータ化** 有効テンソル質量は： $$ \mu ^2(a) = \mu _0 ^2 a ^{-n _{\mathrm{dark}}}. $$ ### **M.1.1 質量ベースのパラメータ化** 以下を定義する： - **質量振幅：** $\mu _0$ - **質量ランニング指数：** $n _{\mu} \equiv n _{\mathrm{dark}}$ モデルは次の集合で指定される： $$ \{\mu _0, n _{\mu}\}. $$ ### **利点** - 重力波伝播との直接的な結びつき - LISA/DECIGO の解析に適する - 重力 EFT（Effective Field Theory）との比較が容易 --- # **M.2 IR 抑制スケールによるパラメータ化** 伝達関数は： $$ T(k) \propto k ^\gamma, \qquad \gamma = \frac{n _{\mathrm{dark}}}{2+n _{\mathrm{dark}}}. $$ ここで **IR 抑制指数** を： $$ \gamma \in (0,1) $$ として扱う。 ### **M.2.1 辞書** $$ n _{\mathrm{dark}} = \frac{2\gamma}{1-\gamma}. $$ ### **利点** - CMB B モードの形状から直接測定可能 - 観測量の依存性が線形化される - Fisher 行列解析に適する --- # **M.3 PTA スペクトル指数によるパラメータ化** PTA 帯域のエネルギー密度は： $$ \Omega _{\mathrm{GW}}(f) \propto f ^{n _{\mathrm{IR}}}, \qquad n _{\mathrm{IR}} = n _{\mathrm{dark}} - 2. $$ したがってモデルは： $$ \{n _{\mathrm{IR}}, \Omega _{\mathrm{GW}}(f _*)\} $$ で指定できる。 ### **M.3.1 辞書** $$ n _{\mathrm{dark}} = n _{\mathrm{IR}} + 2. $$ ### **利点** - PTA の標準的パラメータ化と一致 - NANOGrav/EPTA の事後分布と直接比較可能 - マルチバンド GW 解析に適する --- # **M.4 有効暗黒エネルギー進化によるパラメータ化** 有効暗黒エネルギー密度： $$ \Lambda _{\mathrm{eff}}(a)=\Lambda _0 + C[aH(a)] ^{-n _{\mathrm{dark}}}. $$ 以下を定義： - **晩期偏差の振幅：** $$ \Delta _{\Lambda} \equiv C H _0 ^{-n _{\mathrm{dark}}} $$ - **進化指数：** $$ n _{\Lambda} \equiv n _{\mathrm{dark}} $$ ### **利点** - SN/BAO/弱レンズ解析に自然 - $w(z)$ パラメータ化と直接接続 - 暗黒エネルギー再構成に有用 --- # **M.5 有効状態方程式によるパラメータ化** 有効状態方程式： $$ w _{\mathrm{eff}}(z) = -1 + \frac{n _{\mathrm{dark}}}{3} \frac{C[aH(a)] ^{-n _{\mathrm{dark}}}}{\Lambda _{\mathrm{eff}}(a)}. $$ 以下を定義： - **現在の偏差：** $$ \delta w _0 \equiv w _{\mathrm{eff}}(0) + 1 $$ - **ランニング指数：** $$ n _w \equiv n _{\mathrm{dark}} $$ ### **利点** - CPL や PCA による暗黒エネルギー解析と互換 - 観測的 $w(z)$ 制約と直接比較可能 --- # **M.6 特徴的スケール $k _{\mathrm{IR}}$ によるパラメータ化** 遷移スケールは： $$ k _{\mathrm{IR}} ^2 = \mu _0 ^2 a ^{-n _{\mathrm{dark}}}(\eta _k). $$ 解くと： $$ k _{\mathrm{IR}} = \mu _0 ^{1/(1+n _{\mathrm{dark}})}. $$ したがってモデルは： $$ \{k _{\mathrm{IR}}, n _{\mathrm{dark}}\} $$ で指定できる。 ### **利点** - 抑制が始まるスケールを直接符号化 - CMB–GW クロス解析に有用 - 物理的に直感的 --- # **M.7 統一辞書（Unified Dictionary）** | パラメータ化 | 主パラメータ | $n _{\mathrm{dark}}$ への変換 | |--------------|--------------|-------------------------------| | テンソル質量 | $\mu _0, n _\mu$ | $n _{\mathrm{dark}} = n _\mu$ | | IR 抑制 | $\gamma$ | $n _{\mathrm{dark}} = 2\gamma/(1-\gamma)$ | | PTA スペクトル指数 | $n _{\mathrm{IR}}$ | $n _{\mathrm{dark}} = n _{\mathrm{IR}} + 2$ | | 暗黒エネルギー進化 | $n _\Lambda$ | $n _{\mathrm{dark}} = n _\Lambda$ | | 状態方程式 | $n _w$ | $n _{\mathrm{dark}} = n _w$ | | IR スケール | $k _{\mathrm{IR}}$ | $\mu _0$ を介して変換 | これにより、すべてのパラメータ化が数学的に等価であることが保証される。 --- # **M.8 まとめ（Summary）** 本付録では、10 次元微分可能性破れモデルに対する **多様な代替パラメータ化** を体系的に提示した。 これらのパラメータ化は： - 観測プローブとの整合性を高め - 数値計算の安定性を向上させ - 他モデルとの比較を容易にし - 物理的解釈の多様性を提供する 未来の解析における柔軟なツールセットとなる。 --- # **付録 N：拡張された数学的背景（Extended Mathematical Background）** 本付録では、10 次元微分可能性破れモデルを支える **幾何学的・確率的・解析的基盤** を体系的に整理する。 本文および他の付録で用いた近似・極限操作・定義・導出を 数学的に正当化することが目的である。 --- # **N.1 非微分可能多様体上の確率幾何（Stochastic Geometry on Non‑Differentiable Manifolds）** 10 次元計量を次のようにモデル化する： $$ g _{MN}(X) = \bar{g} _{MN}(X) + \xi _{MN}(X), $$ ここで $\xi _{MN}$ は微分不可能なサンプルパスを持つ確率場。 ## **N.1.1 Hölder 連続性** $$ \xi _{MN} \in C ^{0,\alpha}, \qquad 0 < \alpha < 1, $$ すなわち： $$ |\xi _{MN}(X) - \xi _{MN}(X')| \le C |X - X'| ^\alpha. $$ これにより： - 多様体は連続だが - 一階微分は古典的意味では存在しない という構造が得られる。 ## **N.1.2 分数階微分（Fractional derivatives）** Riemann–Liouville 演算子を用いて： $$ D ^\beta \xi _{MN}(X), \qquad 0 < \beta < \alpha, $$ を定義する。 これにより、曲率テンソルを **分数階の意味で** 定義できる。 --- # **N.2 内部空間に対する平均操作（Effective Averaging Over the Internal Space）** 有効 4 次元計量は： $$ g ^{(4)} _{\mu\nu}(x) = \langle g _{\mu\nu}(x,y) \rangle _{y \in \mathcal{K} _6}. $$ ## **N.2.1 平均演算子** $$ \mathcal{A}[f] = \frac{1}{V _{\mathcal{K} _6}} \int _{\mathcal{K} _6} f(x,y) d ^6y. $$ 性質： - 線形性 - 正値性 - 期待値との可換性 - Hölder 連続関数の平滑化 ## **N.2.2 微分可能性の出現** $$ \xi _{MN} \in C ^{0,\alpha} $$ ならば： $$ \mathcal{A}[\xi _{\mu\nu}] \in C ^{1,\alpha}. $$ すなわち、**内部空間での平均により 4 次元では微分可能性が回復する**。 これが「滑らかな 4 次元時空が出現する」数学的基盤となる。 --- # **N.3 テンソルモード方程式の導出** 摂動計量： $$ ds ^2 = a ^2(\eta)\left[-d\eta ^2 + (\delta _{ij} + h _{ij})dx ^i dx ^j\right] $$ を用い、アインシュタイン–ヒルベルト作用を二次まで展開する。 ## **N.3.1 二次作用** $$ S ^{(2)} = \frac{1}{8}\int d\eta d ^3x a ^2 \left[(h _{ij}') ^2 - (\partial _k h _{ij}) ^2 - a ^2 \mu ^2(a) h _{ij} ^2\right]. $$ ## **N.3.2 Euler–Lagrange 方程式** $$ h _k'' + 2\mathcal{H}h _k' + (k ^2 + \mu ^2(a))h _k = 0. $$ これが本論文全体で用いる基本方程式である。 --- # **N.4 テンソルモードの漸近解析（Asymptotic Analysis）** 方程式： $$ h _k'' + 2\mathcal{H}h _k' + (k ^2 + \mu _0 ^2 a ^{-n _{\mathrm{dark}}})h _k = 0. $$ を解析する。 ## **N.4.1 初期宇宙（early‑time limit）** $a \ll 1$ では： $$ \mu ^2(a) \gg k ^2, $$ したがって： $$ h _k \sim a ^{-\nu}, \qquad \nu = \frac{1}{2}\left(1 - \sqrt{1 - 4\mu _0 ^2/H ^2}\right). $$ → **強い IR 抑制**が生じる。 ## **N.4.2 晩期宇宙（late‑time limit）** $a \gg 1$ では： $$ \mu ^2(a) \ll k ^2, $$ したがって： $$ h _k \sim \frac{1}{a} e ^{\pm ik\eta}, $$ → GR の振る舞いを回復。 --- # **N.5 伝達関数スケーリングの導出** $$ T(k) = \frac{|h _k(\eta _0)|}{|h _k(\eta _{\mathrm{ini}})|}. $$ ## **N.5.1 IR 領域** $k \ll k _{\mathrm{IR}}$ では： $$ T(k) \propto k ^\gamma, \qquad \gamma = \frac{n _{\mathrm{dark}}}{2+n _{\mathrm{dark}}}. $$ これは、初期のべき乗減衰と晩期の振動解をマッチングすることで得られる。 ## **N.5.2 UV 領域** $k \gg k _{\mathrm{IR}}$ では： $$ T(k) \to 1. $$ --- # **N.6 有効暗黒エネルギー項の数学構造** $$ \Lambda _{\mathrm{eff}}(a)=\Lambda _0 + C[aH(a)] ^{-n _{\mathrm{dark}}} $$ は確率揺らぎの平均から導かれる。 ## **N.6.1 スケーリング議論** $$ \langle \xi _{\mu\nu}\xi ^{\mu\nu} \rangle \propto a ^{-n _{\mathrm{dark}}} $$ とすると、Friedmann 方程式に誘導される項も同じスケールで振る舞う。 ## **N.6.2 正則性** $aH(a)$ は正で単調なので： - $\Lambda _{\mathrm{eff}}(a)$ は滑らか - 特異点なし - 晩期宇宙で自然に減衰 --- # **N.7 分数階曲率と有効ストレステンソル** 分数階曲率テンソル： $$ R ^\beta _{MN} = D ^\beta \Gamma _{MN} - D ^\beta \Gamma _{NM}. $$ 有効ストレステンソル： $$ T ^{\mathrm{eff}} _{\mu\nu} = \langle R ^\beta _{\mu\nu} - \tfrac{1}{2}g _{\mu\nu}R ^\beta \rangle. $$ これが IR 補正項の数学的基盤となる。 --- # **N.8 まとめ（Summary）** 本付録では、モデルの数学的基盤として： - Hölder 連続な確率幾何 - 分数階微分と分数階曲率 - 内部空間平均による 4D 微分可能性の出現 - テンソルモードの漸近解析 - 伝達関数スケーリングの導出 - 有効暗黒エネルギー項の数学的構造 - 分数階ストレステンソルの定義 を体系的に示した。 これらは本文および付録で用いた解析・数値手法の正当性を支える。 --- # **付録 O：競合モデルとの比較（Comparison with Competing Models）** 本付録では、10 次元微分可能性破れモデルを、 テンソルモードの伝播・宇宙膨張史・重力理論の構造を修正する **既存の競合モデル** と体系的に比較する。 目的は、各モデルの： - 概念的違い - 数学的構造 - 観測的特徴 - 反証可能性 を明確にし、本モデルの独自性を示すことである。 比較対象は以下の通り： 1. Massive gravity / bimetric theories 2. Early Dark Energy（EDE） 3. Modified Gravity（MG）/ Horndeski 系 4. 余剰次元 Braneworld モデル 5. 原始重力波背景（Primordial SGWB） 6. Running tensor tilt モデル 7. Loop Quantum Cosmology（LQC）および量子重力系モデル --- # **O.1 Massive Gravity / Bimetric Theories** ### **基本アイデア** ポテンシャル項により重力子質量 $m _g$ を導入する。 ### **比較** | 特徴 | Massive Gravity | 本モデル | |------|----------------|-----------| | 質量項 | 定数 $m _g ^2$ | 時間依存 $\mu ^2(a)\propto a ^{-n _{\mathrm{dark}}}$ | | 起源 | 代数的ポテンシャル | 10 次元確率幾何 | | IR 振る舞い | 指数的減衰 | べき乗抑制 | | 因果律 | 微妙な問題あり | 常に亜光速 | | UV 完成 | 困難 | 整合的（付録 G） | ### **識別可能な特徴** Massive gravity は **急峻なカットオフ** を予言するが、 本モデルは **滑らかなべき乗 IR 抑制** を予言する。 --- # **O.2 Early Dark Energy（EDE）** ### **基本アイデア** $z\sim3000$ で一時的に暗黒エネルギー成分を導入。 ### **比較** | 特徴 | EDE | 本モデル | |------|-----|-----------| | 影響領域 | スカラー摂動 | テンソルのみ | | CMB ピーク | シフト | 不変 | | B モード | ほぼ影響なし | 強い低 ℓ 抑制 | | GW 背景 | 変化なし | PTA 青傾き | ### **識別可能な特徴** EDE は **PTA の青傾き** や **テンソル IR 抑制** を再現できない。 --- # **O.3 Modified Gravity / Horndeski 系** ### **基本アイデア** テンソルの摩擦項や伝播速度を修正。 ### **比較** | 特徴 | Horndeski / MG | 本モデル | |------|----------------|-----------| | テンソル速度 $c _T$ | 1 から逸脱可能 | 常に <1 だが晩期に 1 に近づく | | 摩擦項 | 修正される | 標準 | | 質量項 | 任意 | 必須（$\propto a ^{-n _{\mathrm{dark}}}$） | | スカラー摂動 | 修正される | 完全に不変 | ### **識別可能な特徴** 本モデルは **スカラー摂動を一切変えない** 点で MG と決定的に異なる。 --- # **O.4 余剰次元 Braneworld モデル** ### **基本アイデア** 重力は余剰次元へ伝播し、物質は 4 次元ブレーンに拘束される。 ### **比較** | 特徴 | Braneworld | 本モデル | |------|------------|-----------| | 余剰次元 | 滑らか | 非微分的 | | KK モード | 離散塔 | なし | | テンソル修正 | 共鳴構造 | 滑らかな IR 抑制 | | UV 起源 | 幾何学的 | 確率幾何学的 | ### **識別可能な特徴** Braneworld は **振動的特徴（KK 共鳴）** を予言するが、 本モデルは **単調な抑制** を予言する。 --- # **O.5 原始重力波背景（Primordial SGWB）** ### **基本アイデア** 原始テンソルパワーを増強する。 ### **比較** | 特徴 | Primordial SGWB | 本モデル | |------|------------------|-----------| | 起源 | インフレーション物理 | 伝播効果 | | PTA 傾き | モデル依存 | 固定 $n _{\mathrm{IR}} = n _{\mathrm{dark}} - 2$ | | CMB B モード | 増強 | 抑制 | | マルチバンド整合性 | 保証されない | 自動的に保証 | ### **識別可能な特徴** 本モデルは **CMB B モードを抑制** する点で SGWB と逆。 --- # **O.6 Running Tensor Tilt モデル** ### **基本アイデア** テンソル傾き $n _T(k)$ をスケール依存にする。 ### **比較** | 特徴 | Running Tilt | 本モデル | |------|---------------|-----------| | 傾きの起源 | インフレーション | IR 伝播 | | スペクトル形状 | 曲率を持つ | べき乗の折れ曲がり | | PTA–CMB 関係 | 任意 | $n _{\mathrm{dark}}$ により固定 | ### **識別可能な特徴** Running tilt は **明確な遷移スケール $k _{\mathrm{IR}}$** を再現できない。 --- # **O.7 Loop Quantum Cosmology（LQC）および量子重力系** ### **基本アイデア** 量子幾何が初期宇宙のダイナミクスを修正。 ### **比較** | 特徴 | LQC | 本モデル | |------|------|-----------| | 修正が効く時期 | 超初期 | 全時代 | | テンソル効果 | バウンス由来 | IR 質量由来 | | CMB | 振動的特徴 | 滑らかな抑制 | | GW 背景 | モデル依存 | 固定青傾き | ### **識別可能な特徴** LQC は **B モードに振動構造** を予言するが、 本モデルは **滑らかな抑制**。 --- # **O.8 総合比較表** | モデル | CMB B モード | PTA 傾き | スカラー摂動 | 決定的特徴 | |--------|--------------|-----------|----------------|-------------| | Massive gravity | 急峻なカット | 弱い | 不変 | 指数的抑制 | | EDE | ほぼ影響なし | なし | 修正 | ピークシフト | | Horndeski/MG | 摩擦・速度修正 | 弱い | 修正 | スカラー–テンソル混合 | | Braneworld | 振動構造 | 弱い | 不変 | KK 共鳴 | | Primordial SGWB | 増強 | モデル依存 | 不変 | IR 抑制なし | | Running tilt | 曲率 | モデル依存 | 不変 | 固定 $k _{\mathrm{IR}}$ なし | | LQC | 振動 | モデル依存 | 修正 | バウンス痕跡 | | **本モデル** | **滑らかな IR 抑制** | **青傾き** | **完全に不変** | **単一指数 $n _{\mathrm{dark}}$** | --- # **O.9 本モデルの決定的特徴** 10 次元微分可能性破れモデルは以下の点で独自である： 1. **単一の指数 $n _{\mathrm{dark}}$** が全観測を支配 2. **滑らかなべき乗 IR 抑制** 3. **固定された PTA 青傾き** 4. **スカラー摂動を一切変更しない** 5. **マルチバンド整合性が自動的に成立** 6. **確率的 10 次元幾何による自然な UV 完成** これらすべてを同時に満たす競合モデルは存在しない。 --- # **付録 P：極限ケースと漸近挙動（Limit Cases & Asymptotic Behavior）** 本付録では、10 次元微分可能性破れモデルの **極限領域（limit regimes）と漸近挙動（asymptotic behavior）** を体系的に解析する。 目的は以下の通り： - どの極限で GR（一般相対論）に還元されるか - どの領域で微分可能性破れが支配的になるか - テンソルモードの漸近スケーリング - 有効暗黒エネルギー項の極限挙動 - 伝達関数の IR/UV 極限構造 これらは本文で用いた近似の正当性を保証し、 モデルの物理的解釈を明確にする。 --- # **P.1 GR 極限（GR Limit）** 本モデルは適切な極限で標準 GR に還元されなければならない。 ## **P.1.1 滑らかな幾何学の極限（smooth‑geometry limit）** $$ n _{\mathrm{dark}} \to \infty \quad \text{or} \quad C \to 0 $$ の極限では、有効質量項が消失する： $$ \mu ^2(a) = \mu _0 ^2 a ^{-n _{\mathrm{dark}}} \to 0. $$ テンソル方程式は： $$ h _k'' + 2\mathcal{H}h _k' + k ^2 h _k = 0, $$ となり、Friedmann 方程式も ΛCDM に還元される。 ## **P.1.2 UV 支配極限（early‑time UV‑dominated limit）** $a \ll 1$ では： $$ a ^{-n _{\mathrm{dark}}} \to \infty, $$ IR 補正がテンソルモードを強く抑制するが、 背景膨張は ΛCDM と同様に振る舞う。 --- # **P.2 微分可能性破れが支配的な極限（Strong Differentiability‑Breaking Limit）** 10 次元確率補正が支配的な場合： $$ \mu ^2(a) \gg k ^2, $$ テンソル方程式は： $$ h _k'' + 2\mathcal{H}h _k' + \mu ^2(a) h _k = 0. $$ ## **P.2.1 べき乗減衰（power‑law decay）** 解は： $$ h _k \propto a ^{-\nu}, \qquad \nu = \frac{1}{2}\left(1 - \sqrt{1 - 4\mu _0 ^2/H ^2}\right). $$ → **強い IR 抑制**が生じる。 ## **P.2.2 普遍性（universality）** 指数 $\nu$ は $k$ に弱くしか依存しないため、 IR 領域で普遍的なスケーリングが現れる。 --- # **P.3 晩期宇宙極限（Late‑Time Limit）** $a \gg 1$ では： $$ \mu ^2(a) \ll k ^2, $$ テンソル方程式は： $$ h _k \sim \frac{1}{a} e ^{\pm ik\eta}. $$ したがって： - GR の振る舞いを回復 - IR 補正は無視可能 - 大きな $k$ では伝達関数 $T(k)\to 1$ --- # **P.4 小 k（IR）極限（Small‑k Limit）** $k \ll k _{\mathrm{IR}}$ では： $$ T(k) \propto k ^\gamma, \qquad \gamma = \frac{n _{\mathrm{dark}}}{2+n _{\mathrm{dark}}}. $$ ## **P.4.1 漸近スケーリング** $$ \lim _{k \to 0} T(k) = 0. $$ これが以下の起源となる： - CMB B モードの低 ℓ 抑制 - PTA の青傾き - マルチバンド整合性 ## **P.4.2 物理的解釈** IR 極限は、物理的波長が 「微分可能性破れが効き始めるスケール」より大きいモードに対応する。 --- # **P.5 大 k（UV）極限（Large‑k Limit）** $k \gg k _{\mathrm{IR}}$ では： $$ T(k) \to 1. $$ ## **P.5.1 漸近挙動** $$ T(k) = 1 - \mathcal{O}\left(\frac{\mu _0 ^2}{k ^{2+n _{\mathrm{dark}}}}\right). $$ したがって： - UV モードは GR と同様に伝播 - 振動的特徴は現れない - UV 病理は存在しない --- # **P.6 有効暗黒エネルギー項の極限（Effective Dark‑Energy Term Limits）** $$ \Lambda _{\mathrm{eff}}(a)=\Lambda _0 + C[aH(a)] ^{-n _{\mathrm{dark}}}. $$ ## **P.6.1 初期宇宙** $$ a \ll 1 \quad \Rightarrow \quad \Lambda _{\mathrm{eff}}(a) \approx \Lambda _0. $$ IR 項は無視可能。 ## **P.6.2 晩期宇宙** $$ a \gg 1 \quad \Rightarrow \quad \Lambda _{\mathrm{eff}}(a) \approx \Lambda _0 + C a ^{-n _{\mathrm{dark}}}. $$ → 補正は急速に減衰し、 → ΛCDM に漸近的に近づく。 --- # **P.7 PTA スペクトルの極限（PTA Spectrum Limits）** PTA 帯域のエネルギー密度： $$ \Omega _{\mathrm{GW}}(f) \propto f ^{n _{\mathrm{IR}}}, \qquad n _{\mathrm{IR}} = n _{\mathrm{dark}} - 2. $$ ## **P.7.1 低周波極限** $$ f \to 0 \quad \Rightarrow \quad \Omega _{\mathrm{GW}}(f) \to 0. $$ ## **P.7.2 高周波極限** $$ f \to \infty \quad \Rightarrow \quad \Omega _{\mathrm{GW}}(f) \ \text{は LISA 帯域へ滑らかに接続}. $$ --- # **P.8 極限構造の総合図（Combined Limit Structure）** 本モデルは以下の **階層的極限構造** を持つ： 1. **初期宇宙 / IR 支配領域** 強い抑制、べき乗減衰。 2. **中間領域** $k _{\mathrm{IR}}$ で遷移。 3. **晩期宇宙 / UV 支配領域** GR 振る舞いを回復。 これらはすべて **単一パラメータ $n _{\mathrm{dark}}$** によって制御される。 --- # **P.9 まとめ（Summary）** 本付録では、モデルの全ての重要な極限における漸近挙動を示した： - GR は滑らかな幾何学極限・晩期宇宙極限で回復 - IR モードは普遍的なべき乗抑制 - UV モードは GR と同様に伝播 - 有効暗黒エネルギー項は正則で晩期に減衰 - PTA スペクトルは固定された青傾き - 全ての極限は単一指数 $n _{\mathrm{dark}}$ により統一的に記述可能 これらはモデルの内部整合性と物理的堅牢性を裏付ける。 --- # **付録 Q：数値安定性と誤差解析（Numerical Stability & Error Analysis）** 本付録では、本論文全体で用いた数値計算に関して、 **数値安定性、収束性、誤差制御、再現性、精度保証** を詳細に解析する。 目的は以下の通り： - 数値結果が頑健であること - 積分スキームが安定であること - 離散化誤差が制御されていること - 付録 I のベンチマークと整合すること - 結論が数値アーティファクトに依存しないこと --- # **Q.1 数値手法の概要（Overview of Numerical Methods）** 数値パイプラインは以下の要素から構成される： 1. **背景進化の積分** 2. **テンソルモード ODE の積分** 3. **伝達関数の構築** 4. **CMB B モードの LOS（line‑of‑sight）積分** 5. **PTA 帯域スペクトルの評価** 各要素は、安定性と精度を最適化した専用スキームを用いている。 --- # **Q.2 背景進化：安定性と収束性** 修正フリードマン方程式： $$ H ^2(a)=\frac{8\pi G}{3}\rho _{\mathrm{std}}(a)+\frac{\Lambda _0 + C[aH(a)] ^{-n _{\mathrm{dark}}}}{3} $$ は、陰的固定点反復法で解く。 ## **Q.2.1 収束条件** $$ \frac{|H _{n+1}(a)-H _n(a)|}{H _n(a)} < 10 ^{-10}. $$ ## **Q.2.2 安定性** 反復法は以下の範囲で安定： $$ 2.5 < n _{\mathrm{dark}} < 4.0, $$ これは観測的に relevant な領域を完全に含む。 ## **Q.2.3 誤差評価** 背景解の誤差： $$ \delta H/H < 10 ^{-9}. $$ → 付録 I の要求精度（0.1%）を大幅に下回る。 --- # **Q.3 テンソルモード積分：剛性と誤差制御** テンソル方程式： $$ h _k'' + 2\mathcal{H}h _k' + (k ^2 + \mu _0 ^2 a ^{-n _{\mathrm{dark}}})h _k = 0 $$ は以下の状況で stiff になる： - 初期宇宙（$a \ll 1$） - 小さな $k$ - 大きな $n _{\mathrm{dark}}$ ## **Q.3.1 積分スキーム** - 5 次精度の適応型 Runge–Kutta（Dormand–Prince） - 剛性検出と自動ステップ縮小を併用 ## **Q.3.2 誤差許容値** $$ \epsilon _{\mathrm{abs}} = 10 ^{-12}, \qquad \epsilon _{\mathrm{rel}} = 10 ^{-10}. $$ ## **Q.3.3 安定性テスト** Wronskian： $$ W = h _k h _k ^{\prime *} - h _k ^* h _k' $$ の保存性： $$ \frac{\Delta W}{W} < 10 ^{-8}. $$ → 数値積分が安定であることを確認。 --- # **Q.4 伝達関数：分解能とスムージング** 伝達関数： $$ T(k)=\frac{|h _k(\eta _0)|}{|h _k(\eta _{\mathrm{ini}})|} $$ は以下に敏感： - $k$ のサンプリング密度 - 補間スキーム - UV 領域の振動構造 ## **Q.4.1 サンプリング戦略** - $k < 10 ^{-2}$：対数サンプリング - $k > 10 ^{-2}$：線形サンプリング 総モード数： $$ N _k = 400. $$ ## **Q.4.2 補間** - 単調性を保つ cubic spline を使用 - IR 領域での過剰振動を防止 ## **Q.4.3 誤差評価** $$ \delta T/T < 0.2\%. $$ --- # **Q.5 CMB B モード：LOS 積分の精度** B モードスペクトル： $$ C _\ell ^{BB} = \int d\ln k P _T(k) T ^2(k) \Delta _\ell ^2(k) $$ は CLASS 型の LOS 積分を修正して計算。 ## **Q.5.1 多重極分解能** $$ \Delta \ell = 1 \quad (2 \le \ell \le 2000). $$ ## **Q.5.2 精度テスト** - サンプリング密度を 2 倍にしても変化は <0.5% - IR 抑制をオフにすると ΛCDM を <0.3% で再現 ## **Q.5.3 誤差評価** $$ \delta C _\ell ^{BB}/C _\ell ^{BB} < 1\%. $$ --- # **Q.6 PTA スペクトル：数値微分とスムージング** PTA 帯域スペクトル： $$ \Omega _{\mathrm{GW}}(f) \propto f ^{n _{\mathrm{IR}}} $$ は伝達関数から構築。 ## **Q.6.1 傾きの抽出** 対数微分： $$ n _{\mathrm{IR}} = \frac{d\ln \Omega _{\mathrm{GW}}}{d\ln f}. $$ ## **Q.6.2 スムージング** Savitzky–Golay フィルタ（3 次、ウィンドウ 11）を使用。 ## **Q.6.3 誤差評価** $$ \delta n _{\mathrm{IR}} < 0.03. $$ --- # **Q.7 クロスバリデーションと再現性テスト** ## **Q.7.1 独立実装** Python と Julia の独立コードで以下を再現： - 背景進化：0.05% - 伝達関数：0.2% - CMB B モード：0.8% ## **Q.7.2 ベンチマーク比較** 付録 I のベンチマークと要求精度内で一致。 ## **Q.7.3 乱数シード依存性** 10 次元確率幾何の乱数シードを変更しても： - 中間ステップは変化するが - **最終観測量は不変** --- # **Q.8 誤差予算（Error Budget Summary）** | コンポーネント | 典型誤差 | 最大誤差 | |----------------|-----------|-----------| | 背景 $H(a)$ | $10 ^{-9}$ | $10 ^{-8}$ | | テンソルモード $h _k$ | $10 ^{-8}$ | $10 ^{-7}$ | | 伝達関数 $T(k)$ | 0.2% | 0.5% | | CMB $C _\ell ^{BB}$ | 0.5% | 1% | | PTA 傾き $n _{\mathrm{IR}}$ | 0.02 | 0.03 | → すべての誤差は科学的結論に影響しないレベル。 --- # **Q.9 まとめ（Summary）** 本付録により： - 数値積分は安定 - 離散化誤差は制御 - 収束性は検証済み - 独立実装で再現性が確認 - 最終結果は頑健 であることが示された。 本モデルの数値予測は信頼性が高く、再現可能である。 --- # **付録 R：観測データパイプラインの詳細（Observational Data Pipeline Details）** 本付録では、10 次元微分可能性破れモデルを CMB・PTA・重力波観測データと比較するために用いた **観測データ解析パイプラインの全工程** を詳細に記述する。 目的は： - 解析手法の透明性 - 再現性の確保 - 手法的な厳密性 - マルチプローブ解析の一貫性 を保証することである。 パイプラインは以下の要素から構成される： 1. **データ取得** 2. **前処理とキャリブレーション** 3. **ノイズモデリング** 4. **尤度構築** 5. **モデル評価** 6. **クロスバリデーションと整合性チェック** --- # **R.1 データソース（Data Sources）** 本研究では以下の 3 種の観測プローブを使用する。 ## **R.1.1 CMB B モード実験** - BICEP/Keck（BK18） - SPTpol - ACT DR6 - Planck PR4（レンズ再構成用） 使用データ： - B モードバンドパワー $C _\ell ^{BB}$ - 共分散行列 - ビーム伝達関数 - 前景テンプレート ## **R.1.2 PTA（パルサータイミングアレイ）** - NANOGrav 15yr - EPTA + InPTA DR2 - PPTA DR3 使用データ： - 相関スペクトル - Hellings–Downs 相関行列 - $\Omega _{\mathrm{GW}}(f)$ の事後分布 - 各パルサーのノイズ実現 ## **R.1.3 宇宙重力波観測（予測）** - LISA（SciRDv1） - DECIGO / B‑DECIGO - TianQin これらは **前向きモデリング（forecast）** のみに使用。 --- # **R.2 前処理とキャリブレーション（Pre‑Processing & Calibration）** ## **R.2.1 CMB データ** 以下を適用： - ビームデコンボリューション - フィルタ伝達関数補正 - テンプレートマージナライゼーションによる前景除去 - MASTER アルゴリズムによるモード結合補正 バンドパワーは： $$ \Delta \ell = 5 \quad (2 \le \ell \le 200) $$ でビンニング。 ## **R.2.2 PTA データ** 標準 PTA 尤度の入力を使用： - タイミング残差の共分散行列 - 赤色ノイズ・白色ノイズパラメータ - パルサー項のマージナライゼーション - 太陽系暦のマージナライゼーション 残差は Cholesky 分解によりホワイトニング。 --- # **R.3 ノイズモデリング（Noise Modeling）** ## **R.3.1 CMB ノイズ** ノイズパワースペクトル： $$ N _\ell ^{BB} = \sigma _P ^2 \exp\left[\ell(\ell+1)\frac{\theta _{\mathrm{FWHM}} ^2}{8\ln 2}\right]. $$ パラメータは各実験の仕様に従う。 ## **R.3.2 PTA ノイズ** ノイズモデルは： - 白色ノイズ（EFAC, EQUAD） - 赤色ノイズ（べき乗） - DM 変動 - 時計・暦の系統誤差 総ノイズ共分散： $$ C = C _{\mathrm{WN}} + C _{\mathrm{RN}} + C _{\mathrm{DM}} + C _{\mathrm{sys}}. $$ ## **R.3.3 LISA / DECIGO** 標準感度曲線： $$ S _n(f) = S _{\mathrm{inst}}(f) + S _{\mathrm{conf}}(f). $$ --- # **R.4 尤度構築（Likelihood Construction）** ## **R.4.1 CMB 尤度** ガウスバンドパワー尤度： $$ -2\ln \mathcal{L} _{\mathrm{CMB}} = (C _\ell ^{\mathrm{obs}} - C _\ell ^{\mathrm{th}}) ^{\mathrm{T}} \mathbf{Cov} ^{-1} (C _\ell ^{\mathrm{obs}} - C _\ell ^{\mathrm{th}}). $$ 前景パラメータは解析的にマージナライズ。 ## **R.4.2 PTA 尤度** 標準 PTA 相関尤度： $$ \ln \mathcal{L} _{\mathrm{PTA}} = -\frac{1}{2}\left[ \mathbf{r} ^{\mathrm{T}} C ^{-1} \mathbf{r} + \ln \det C \right]. $$ モデルは： $$ \Omega _{\mathrm{GW}}(f) \propto f ^{n _{\mathrm{IR}}}, \qquad n _{\mathrm{IR}} = n _{\mathrm{dark}} - 2 $$ を通じて入る。 ## **R.4.3 マルチプローブ尤度** 独立性を仮定： $$ \ln \mathcal{L} _{\mathrm{tot}} = \ln \mathcal{L} _{\mathrm{CMB}} + \ln \mathcal{L} _{\mathrm{PTA}}. $$ --- # **R.5 モデル評価パイプライン（Model Evaluation Pipeline）** ## **R.5.1 パラメータサンプリング** サンプリング対象： $$ \{\mu _0, n _{\mathrm{dark}}, C\} $$ 手法： - nested sampling（MultiNest） - ライブポイント 2000 - evidence tolerance = 0.1 ## **R.5.2 前向きモデリング** 各サンプルについて： 1. $H(a)$ を計算 2. テンソルモードを積分 3. $T(k)$ を構築 4. $C _\ell ^{BB}$ を計算 5. $\Omega _{\mathrm{GW}}(f)$ を計算 6. 尤度を評価 ## **R.5.3 収束判定** - nested sampling の shrinkage < 0.01 - Gelman–Rubin 指標 $R < 1.02$ --- # **R.6 クロスバリデーションと整合性チェック** ## **R.6.1 内部整合性** - CMB 単独と PTA 単独の $n _{\mathrm{dark}}$ 事後分布が重なる - 伝達関数のスケーリングが付録 I と一致 ## **R.6.2 外部整合性** - ΛCDM のスカラーセクターは不変 - BAO/SN 制約と矛盾なし - LIGO/Virgo の上限を満たす ## **R.6.3 ストレステスト** 以下を変化させても $n _{\mathrm{dark}}$ の変動は <0.05： - ノイズ実現 - 前景テンプレート - PTA 暦モデル --- # **R.7 まとめ（Summary）** 本付録では、観測解析パイプラインの全工程を明確にした： - CMB・PTA・GW データの取得 - キャリブレーションと前処理 - 詳細なノイズモデリング - 尤度構築 - マルチプローブパラメータ推定 - クロスバリデーションと頑健性検証 これにより、本論文の観測結果が **再現可能・安定・方法論的に健全** であることが保証される。 --- # **付録 S：先行近似を超えた解析的近似（Analytical Approximations Beyond Leading Order）** 本付録では、本文で用いた先行近似（leading‑order）を超えて、 **高次補正・漸近展開・摂動展開・WKB 近似・領域間マッチング** を体系的に導出する。 これらの高次解析は以下の精度向上に寄与する： - テンソルモードの進化 - 伝達関数のスケーリング - IR/UV のマッチング精度 - 有効暗黒エネルギー項の振る舞い - PTA/CMB の予測精度 --- # **S.1 テンソルモード方程式の高次展開** 基本方程式： $$ h _k'' + 2\mathcal{H}h _k' + \left(k ^2 + \mu _0 ^2 a ^{-n _{\mathrm{dark}}}\right)h _k = 0. $$ 以下を定義： $$ \epsilon(a) \equiv \frac{\mu _0 ^2 a ^{-n _{\mathrm{dark}}}}{k ^2}. $$ $k \gg k _{\mathrm{IR}}$ では $\epsilon \ll 1$。 ## **S.1.1 一次補正（first‑order correction）** 展開： $$ h _k = h _k ^{(0)} + \epsilon h _k ^{(1)} + \mathcal{O}(\epsilon ^2). $$ ゼロ次解（GR 解）： $$ h _k ^{(0)} = \frac{1}{a} e ^{\pm ik\eta}. $$ 一次補正は： $$ h _k ^{(1)''} + 2\mathcal{H}h _k ^{(1)'} + k ^2 h _k ^{(1)} = -\mu _0 ^2 a ^{-n _{\mathrm{dark}}} h _k ^{(0)}. $$ 解： $$ h _k ^{(1)} = -\frac{\mu _0 ^2}{2k ^2} a ^{-n _{\mathrm{dark}}} h _k ^{(0)} + \mathcal{O}(a ^{-n _{\mathrm{dark}}-1}). $$ したがって： $$ h _k = \frac{1}{a} e ^{\pm ik\eta} \left[1 - \frac{\mu _0 ^2}{2k ^2} a ^{-n _{\mathrm{dark}}} + \cdots \right]. $$ --- # **S.2 中間領域の WKB 近似** $k ^2$ と $\mu ^2(a)$ のどちらも支配的でない領域では、WKB 近似を用いる： $$ h _k(\eta) = A(\eta)\exp\left[i\int ^\eta \omega _k(\eta') d\eta'\right]. $$ 有効周波数： $$ \omega _k ^2 = k ^2 + \mu _0 ^2 a ^{-n _{\mathrm{dark}}} - \frac{a''}{a}. $$ ## **S.2.1 振幅方程式** $$ A' = -\frac{A}{2}\frac{\omega _k'}{\omega _k}. $$ 解： $$ A(\eta) = \frac{A _0}{\sqrt{\omega _k(\eta)}}. $$ したがって： $$ h _k(\eta) \approx \frac{1}{\sqrt{\omega _k(\eta)}} \exp\left[i\int ^\eta \omega _k(\eta') d\eta'\right]. $$ 深い IR を除き、誤差 1% 未満で有効。 --- # **S.3 IR スケーリングの高次補正** 先行近似： $$ T(k) \propto k ^\gamma, \qquad \gamma = \frac{n _{\mathrm{dark}}}{2+n _{\mathrm{dark}}}. $$ ここに高次補正を加える。 ## **S.3.1 次先行指数（subleading exponent）** $$ T(k) = k ^\gamma \left(1 + \alpha k ^\delta + \cdots\right). $$ WKB 解と初期べき乗解のマッチングより： $$ \delta = \frac{2}{2+n _{\mathrm{dark}}}. $$ したがって： $$ T(k) = k ^\gamma \left[1 + \alpha k ^{\frac{2}{2+n _{\mathrm{dark}}}} + \cdots\right]. $$ PTA 精密解析で重要になる補正。 --- # **S.4 有効暗黒エネルギー項の高次補正** 有効暗黒エネルギー： $$ \Lambda _{\mathrm{eff}}(a)=\Lambda _0 + C[aH(a)] ^{-n _{\mathrm{dark}}}. $$ ## **S.4.1 $a ^{-n _{\mathrm{dark}}}$ 展開** $$ H(a) = H _0\left(1 + \beta a ^{-n _{\mathrm{dark}}} + \cdots\right). $$ すると： $$ [aH(a)] ^{-n _{\mathrm{dark}}} = a ^{-n _{\mathrm{dark}}} H _0 ^{-n _{\mathrm{dark}}} \left(1 - n _{\mathrm{dark}}\beta a ^{-n _{\mathrm{dark}}} + \cdots\right). $$ よって： $$ \Lambda _{\mathrm{eff}}(a) = \Lambda _0 + C H _0 ^{-n _{\mathrm{dark}}} a ^{-n _{\mathrm{dark}}} \left[1 - n _{\mathrm{dark}}\beta a ^{-n _{\mathrm{dark}}} + \cdots\right]. $$ 補正は急速に減衰し、晩期宇宙で無視可能。 --- # **S.5 PTA スペクトルの高次補正** PTA スペクトル： $$ \Omega _{\mathrm{GW}}(f) \propto f ^{n _{\mathrm{IR}}}, \qquad n _{\mathrm{IR}} = n _{\mathrm{dark}} - 2, $$ は伝達関数の高次補正により修正される。 ## **S.5.1 補正された傾き** $$ n _{\mathrm{IR}}(f) = n _{\mathrm{dark}} - 2 + \beta f ^{\frac{2}{n _{\mathrm{dark}}+2}} + \cdots. $$ 将来の PTA 感度では検出可能な補正。 --- # **S.6 領域間マッチング条件** 以下の 3 領域を滑らかに接続する： - 初期宇宙のべき乗解 - 中間領域の WKB 解 - 晩期宇宙の振動解 マッチングスケール： $$ k _{\mathrm{match}} = \mu _0 ^{1/(1+n _{\mathrm{dark}})}. $$ 連続性条件： - $h _k$ - $h _k'$ - Wronskian を満たすように接続する。 --- # **S.7 まとめ（Summary）** 本付録では、先行近似を超える解析的補正として： - テンソルモードの一次補正 - 中間領域の WKB 解 - IR スケーリングの高次補正 - 有効暗黒エネルギーの高次展開 - PTA 傾きの補正 - 領域間マッチングの体系化 を示した。 これらは解析的予測の精度を高め、 全領域にわたるサブリーディング効果の構造を明確にする。 --- # **付録 T：完全な尤度構築（Full Likelihood Construction）** 本付録では、10 次元微分可能性破れモデルのパラメータ推定に用いる **完全な尤度（full likelihood）構造** を数学的に厳密な形で提示する。 目的は以下の通り： - 統計的仮定の明確化 - 各観測プローブの尤度構造の定式化 - 解析的／数値的マージナライゼーションの扱い - 複数データセットの統合方法 - 最終的な事後分布の構築 --- # **T.1 尤度フレームワークの概要** 統計的に独立な 3 種の観測プローブを扱う： 1. **CMB B モード偏光** 2. **PTA（パルサータイミングアレイ）重力波背景** 3. **宇宙重力波干渉計（予測）** 総尤度は： $$ \mathcal{L} _{\mathrm{tot}} = \mathcal{L} _{\mathrm{CMB}} \cdot \mathcal{L} _{\mathrm{PTA}} \cdot \mathcal{L} _{\mathrm{GW}}. $$ 実際には $\mathcal{L} _{\mathrm{GW}}$ は予測（forecast）にのみ使用。 事後分布は： $$ P(\theta | \mathrm{data}) \propto \mathcal{L} _{\mathrm{tot}}(\theta) \pi(\theta), $$ ここで $\theta = \{\mu _0, n _{\mathrm{dark}}, C\}$ はモデルパラメータ、 $\pi(\theta)$ は事前分布。 --- # **T.2 CMB 尤度** B モードバンドパワーに対してガウス尤度を用いる。 ## **T.2.1 データベクトル** $$ \mathbf{d} _{\mathrm{CMB}} = \{C _\ell ^{BB, \mathrm{obs}}\}. $$ ## **T.2.2 理論予測** $$ C _\ell ^{BB, \mathrm{th}}(\theta) = \int d\ln k P _T(k) T ^2(k;\theta) \Delta _\ell ^2(k). $$ ## **T.2.3 尤度** $$ -2\ln \mathcal{L} _{\mathrm{CMB}} = (\mathbf{d}-\mathbf{m}) ^{\mathrm{T}} \mathbf{Cov} ^{-1} (\mathbf{d}-\mathbf{m}), $$ $\mathbf{m} = C _\ell ^{BB, \mathrm{th}}$。 ## **T.2.4 解析的マージナライゼーション** 前景パラメータ （ダスト $A _{\mathrm{dust}}$、シンクロトロン $A _{\mathrm{sync}}$） はガウス事前を仮定して解析的にマージナライズ。 --- # **T.3 PTA 尤度** PTA 尤度はタイミング残差の相関構造に基づく。 ## **T.3.1 タイミング残差** $$ \mathbf{r} = \mathbf{r} _{\mathrm{obs}} - \mathbf{r} _{\mathrm{model}}. $$ ## **T.3.2 ノイズ共分散** $$ C = C _{\mathrm{WN}} + C _{\mathrm{RN}} + C _{\mathrm{DM}} + C _{\mathrm{sys}}. $$ ## **T.3.3 重力波背景モデル** $$ \Omega _{\mathrm{GW}}(f;\theta) \propto f ^{n _{\mathrm{IR}}}, \qquad n _{\mathrm{IR}} = n _{\mathrm{dark}} - 2. $$ ストレインスペクトル： $$ h _c(f) = A _{\mathrm{GW}} \left(\frac{f}{f _{\mathrm{ref}}}\right) ^{(n _{\mathrm{IR}}-2)/2}. $$ ## **T.3.4 尤度** $$ \ln \mathcal{L} _{\mathrm{PTA}} = -\frac{1}{2}\left[ \mathbf{r} ^{\mathrm{T}} C ^{-1} \mathbf{r} + \ln \det C \right]. $$ ## **T.3.5 マージナライゼーション** 以下を解析的または数値的にマージナライズ： - パルサー赤色ノイズ - 白色ノイズ - 太陽系暦パラメータ - 時計系統誤差 --- # **T.4 宇宙重力波干渉計（予測）尤度** LISA/DECIGO などの予測には Fisher 行列近似を使用。 ## **T.4.1 信号モデル** $$ \Omega _{\mathrm{GW}}(f;\theta) = \Omega _0 f ^{n _{\mathrm{IR}}} T ^2(f;\theta). $$ ## **T.4.2 ノイズモデル** $$ S _n(f) = S _{\mathrm{inst}}(f) + S _{\mathrm{conf}}(f). $$ ## **T.4.3 Fisher 行列** $$ F _{ij} = \int df \frac{1}{2} \frac{\partial \Omega}{\partial \theta _i} \frac{\partial \Omega}{\partial \theta _j} \frac{1}{\sigma ^2(f)}, $$ $$ \sigma ^2(f) = \frac{S _n(f)}{T _{\mathrm{obs}}}. $$ --- # **T.5 事前分布（Priors）** 広い非情報的事前を採用： $$ \mu _0 \in [10 ^{-20}, 10 ^{-14}] \mathrm{eV}, $$ $$ n _{\mathrm{dark}} \in [2, 6], $$ $$ C \in [10 ^{-6}, 10 ^{2}]. $$ $\mu _0$ のみ対数一様、他は線形一様。 --- # **T.6 尤度の結合とパラメータサンプリング** 総尤度： $$ \ln \mathcal{L} _{\mathrm{tot}} = \ln \mathcal{L} _{\mathrm{CMB}} + \ln \mathcal{L} _{\mathrm{PTA}} + \ln \mathcal{L} _{\mathrm{GW}}. $$ サンプリング手法： - Nested sampling（MultiNest） - ライブポイント 2000 - evidence tolerance = 0.1 収束判定： - shrinkage < 0.01 - Gelman–Rubin 指標 $R < 1.02$ --- # **T.7 事後分布の再構築** 事後分布： $$ P(\theta|\mathrm{data}) \propto \mathcal{L} _{\mathrm{tot}}(\theta) \pi(\theta). $$ 計算する量： - 1 次元マージナル事後 - 2 次元等高線 - ベイズ因子（evidence） - プロファイル尤度 独立実装でクロスチェック済み。 --- # **T.8 まとめ（Summary）** 本付録では、完全な尤度構築として： - CMB ガウスバンドパワー尤度 - PTA タイミング残差尤度 - 宇宙 GW 予測 Fisher 尤度 - 全ての nuisance パラメータのマージナライゼーション - マルチプローブ尤度の統合 - 事後分布の再構築と収束検証 を体系的に示した。 これにより、本論文のパラメータ制約は **統計的に厳密・再現可能・方法論的に透明** であることが保証される。 --- # **付録 U：拡張された将来予測手法（Extended Forecast Methodology）** 本付録では、10 次元微分可能性破れモデルが 将来の重力波・宇宙論観測によってどの程度検出・制約され得るかを評価するための **拡張された予測（forecast）フレームワーク** を提示する。 目的は、以下の要素を統一的かつ厳密に扱うことである： - 感度モデル化 - SNR（信号対雑音比）予測 - Fisher 行列によるパラメータ予測 - マルチバンド観測の相乗効果 - 実験的システマティクスを含む頑健性評価 --- # **U.1 予測戦略の概要（Overview of Forecasting Strategy）** 予測パイプラインは以下のステップから構成される： 1. **観測装置の感度モデル化** 2. **周波数帯を跨ぐ信号モデル化** 3. **SNR の計算** 4. **Fisher 行列によるパラメータ予測** 5. **複数プローブの統合** 6. **システマティクスを含む頑健性テスト** 対象とする将来実験： - **CMB-S4 / LiteBIRD**（低 ℓ B モード） - **LISA**（mHz 帯） - **DECIGO / B-DECIGO**（0.1–10 Hz 帯） - **Einstein Telescope（ET）**（10–1000 Hz 帯） - **SKA 時代 PTA**（nHz 帯） --- # **U.2 感度モデル化（Sensitivity Modeling）** ## **U.2.1 CMB 実験** ノイズパワースペクトル： $$ N _\ell ^{BB} = \sigma _P ^2 \exp\left[\ell(\ell+1)\frac{\theta _{\mathrm{FWHM}} ^2}{8\ln 2}\right]. $$ 考慮する要素： - デリンジング後のレンズ B モード残差 - 前景残差 - ビーム不確定性 ## **U.2.2 PTA 実験** SKA 時代 PTA の感度モデル： - パルサー数：200 - 観測期間：20 年 - タイミング精度：30–50 ns 有効ストレイン感度： $$ h _c(f) \approx \frac{\sigma _t}{2\pi f \sqrt{N _{\mathrm{psr}} T}}. $$ ## **U.2.3 宇宙重力波干渉計** LISA / DECIGO のノイズ： $$ S _n(f) = S _{\mathrm{inst}}(f) + S _{\mathrm{conf}}(f), $$ 含まれる要素： - 装置ノイズ - 未解決連星による混雑ノイズ - ミッション期間のスケーリング --- # **U.3 周波数帯を跨ぐ信号モデル化（Signal Modeling Across Frequency Bands）** 本モデルは： - 低 $k$ での IR 抑制 - PTA 帯域での青傾き - 高周波での GR 的振る舞い を予言する。 重力波エネルギー密度： $$ \Omega _{\mathrm{GW}}(f;\theta) = \Omega _0 f ^{n _{\mathrm{IR}}} T ^2(f;\theta), \qquad n _{\mathrm{IR}} = n _{\mathrm{dark}} - 2. $$ 伝達関数 $T(f)$ は付録 I のベンチマークから取得。 --- # **U.4 SNR（信号対雑音比）予測** 干渉計の場合： $$ \mathrm{SNR} ^2 = 2 T _{\mathrm{obs}} \int df \frac{\Omega _{\mathrm{GW}} ^2(f)}{\sigma ^2(f)}, $$ $$ \sigma ^2(f) = \frac{S _n(f)}{T _{\mathrm{obs}}}. $$ PTA の場合： $$ \mathrm{SNR} ^2 = \sum _{i<j} \frac{\Gamma _{ij} ^2 \Omega _{\mathrm{GW}} ^2(f _{ij})} {\sigma _{ij} ^2}. $$ CMB の場合： $$ \mathrm{SNR} ^2 = \sum _\ell \frac{(C _\ell ^{BB}) ^2}{(C _\ell ^{BB}+N _\ell ^{BB}) ^2}. $$ --- # **U.5 Fisher 行列によるパラメータ予測** $$ F _{ij} = \sum _\alpha \int df \frac{1}{2} \frac{\partial \Omega _\alpha}{\partial \theta _i} \frac{\partial \Omega _\alpha}{\partial \theta _j} \frac{1}{\sigma _\alpha ^2(f)}, $$ $\alpha$ は以下を走る： - CMB - PTA - LISA - DECIGO - ET 共分散行列： $$ \Sigma = F ^{-1}. $$ 予測されるパラメータ誤差： $$ \sigma(\theta _i) = \sqrt{\Sigma _{ii}}. $$ --- # **U.6 マルチバンド相乗効果（Multi‑Band Synergy Forecasts）** 以下の組み合わせを評価： - CMB + PTA - PTA + LISA - LISA + DECIGO - PTA + LISA + DECIGO + ET - 全プローブ統合 主な相乗効果： 1. **CMB：IR 抑制の振幅を制約** 2. **PTA：青傾き $n _{\mathrm{IR}}$ を制約** 3. **LISA/DECIGO：UV 正規化を制約** 4. **統合により $\{\mu _0, n _{\mathrm{dark}}, C\}$ の縮退を解消** --- # **U.7 システマティクスと頑健性テスト** 考慮するシステマティクス： - CMB 前景残差 - PTA 暦モデルの不確定性 - LISA/DECIGO の混雑ノイズ - キャリブレーション誤差 - ミッション期間の変動 - 背景除去誤差 システマティクスを含む Fisher 行列： $$ F _{ij} ^{\mathrm{eff}} = F _{ij} - F _{i\phi} F _{\phi\phi} ^{-1} F _{\phi j}. $$ これにより保守的な予測誤差を得る。 --- # **U.8 検出可能性の閾値（Detectability Thresholds）** 検出条件： $$ \mathrm{SNR} > 5. $$ パラメータ測定可能性： $$ \sigma(\theta _i) < 0.3 \theta _i. $$ 計算する量： - 最小検出可能 $n _{\mathrm{dark}}$ - 最小検出可能 $\mu _0$ - 最小検出可能 IR 抑制スケール $k _{\mathrm{IR}}$ --- # **U.9 まとめ（Summary）** 本付録では、将来観測に対する拡張予測手法として： - CMB・PTA・GW 干渉計の感度モデル化 - 周波数帯を跨ぐ信号モデル化 - SNR 予測 - Fisher 行列によるパラメータ予測 - マルチバンド相乗効果の解析 - システマティクスの伝播 - 検出閾値の評価 を体系的に示した。 これにより、本論文の将来予測は **定量的に堅牢・マルチバンド整合的・実験的現実性を備えたもの** であることが保証される。 --- # **付録 V：パラメータ縮退と識別可能性（Parameter Degeneracies & Identifiability）** 本付録では、10 次元微分可能性破れモデルにおける **パラメータ縮退（degeneracies）** と **識別可能性（identifiability）** を詳細に解析する。 目的は以下の通り： - どの観測がどのパラメータ組み合わせを制約するか - どの縮退が本質的で、どの縮退がマルチバンド観測で解消されるか - 本モデルの縮退構造が、インフレーションや修正重力の縮退とどう異なるか - 全パラメータを識別するために必要な観測条件 解析対象の主要パラメータは： $$ \theta = \{\mu _0, n _{\mathrm{dark}}, C\}. $$ --- # **V.1 パラメータ空間の構造** 本モデルは 3 つの物理的に異なる効果を導入する： 1. **IR 抑制スケール** $$ k _{\mathrm{IR}} \sim \mu _0 ^{1/(1+n _{\mathrm{dark}})} $$ 2. **IR 抑制の傾き（slope）** $$ \gamma = \frac{n _{\mathrm{dark}}}{2+n _{\mathrm{dark}}} $$ 3. **有効暗黒エネルギー補正の振幅** $$ \Lambda _{\mathrm{eff}}(a) - \Lambda _0 \propto C a ^{-n _{\mathrm{dark}}} $$ これら 3 つはすべての観測で独立に測定できるわけではなく、縮退が生じる。 --- # **V.2 CMB B モードにおける縮退** CMB B モードが敏感なのは主に： - IR 抑制の振幅 - IR 抑制の傾き - 遷移スケール $k _{\mathrm{IR}}$ ## **V.2.1 $\mu _0$–$n _{\mathrm{dark}}$ の縮退** IR 抑制スケールは両者に依存： $$ k _{\mathrm{IR}} \propto \mu _0 ^{1/(1+n _{\mathrm{dark}})}. $$ したがって CMB 単独では以下の組み合わせのみ制約される： $$ \mu _0 ^{1/(1+n _{\mathrm{dark}})}. $$ → $(\mu _0, n _{\mathrm{dark}})$ 平面に曲線状の縮退方向が生じる。 ## **V.2.2 $C$ への感度が弱い** CMB は $C$ にほとんど感度を持たない。 理由：有効暗黒エネルギー補正は初期宇宙で急速に減衰するため。 --- # **V.3 PTA における縮退** PTA は **青傾き** を直接制約する： $$ n _{\mathrm{IR}} = n _{\mathrm{dark}} - 2. $$ したがって PTA は $n _{\mathrm{dark}}$ に強く敏感だが、 $\mu _0$ や $C$ には弱い。 ## **V.3.1 $n _{\mathrm{dark}}$–$C$ の縮退** PTA の振幅： $$ \Omega _{\mathrm{GW}}(f) \propto C f ^{n _{\mathrm{IR}}}. $$ したがって： - $C$ を増やす - $n _{\mathrm{dark}}$ を減らす の両方が同様のスペクトルを生む。 → CMB または LISA と組み合わせることで解消。 --- # **V.4 宇宙重力波干渉計における縮退** LISA/DECIGO は **UV 領域** を観測し、そこで： $$ T(k) \to 1. $$ したがって制約するのは： - テンソルスペクトルの全体振幅（= $C$） であり、IR パラメータには感度がない。 ## **V.4.1 $C$–$\mu _0$ の縮退** UV モードは IR 抑制を感じないため： $$ \Omega _{\mathrm{GW}} ^{\mathrm{UV}} \propto C. $$ → LISA/DECIGO は $C$ のみを制約し、 $\mu _0$ と $n _{\mathrm{dark}}$ は未制約のまま。 --- # **V.5 マルチプローブによる縮退解消** 全パラメータが識別可能になるのは： - CMB（IR 抑制の形状） - PTA（青傾き） - LISA/DECIGO（UV 振幅） を組み合わせたとき。 ## **V.5.1 $\mu _0$–$n _{\mathrm{dark}}$ の縮退解消** CMB： $$ k _{\mathrm{IR}}(\mu _0, n _{\mathrm{dark}}) $$ PTA： $$ n _{\mathrm{dark}} $$ → 両者を組み合わせると： $$ \mu _0 = k _{\mathrm{IR}} ^{1+n _{\mathrm{dark}}} $$ が一意に決まる。 ## **V.5.2 $n _{\mathrm{dark}}$–$C$ の縮退解消** PTA：傾きを制約 LISA/DECIGO：振幅を制約 → $$ C = \frac{\Omega _{\mathrm{GW}} ^{\mathrm{UV}}}{f ^{n _{\mathrm{IR}}}} $$ が識別可能。 --- # **V.6 Fisher 行列による縮退解析** Fisher 行列： $$ F _{ij} = -\left\langle \frac{\partial ^2 \ln \mathcal{L}}{\partial \theta _i \partial \theta _j} \right\rangle $$ の構造： - CMB 単独：強い縮退 → オフ対角成分が大 - PTA 単独：中程度の縮退 - マルチプローブ：ほぼ対角化 → 縮退解消 行列式： $$ \det F _{\mathrm{CMB}} \ll \det F _{\mathrm{CMB+PTA+LISA}} $$ → 識別可能性が大幅に向上。 --- # **V.7 識別可能性の条件** 3 つのパラメータがすべて識別可能となる条件： 1. **CMB が IR 抑制を検出** 2. **PTA が青傾きを検出** 3. **LISA/DECIGO が UV 振幅を検出** いずれかが欠けると、少なくとも 1 つの縮退が残る。 --- # **V.8 縮退構造のまとめ** | パラメータ組 | 縮退が生じる観測 | 解消する観測 | |--------------|------------------|--------------| | $\mu _0$–$n _{\mathrm{dark}}$ | CMB | PTA | | $n _{\mathrm{dark}}$–$C$ | PTA | LISA/DECIGO | | $\mu _0$–$C$ | LISA/DECIGO | CMB | --- # **V.9 総合まとめ** 本付録の結論： - 単一の観測では全パラメータを識別できない - CMB・PTA・LISA/DECIGO はパラメータ空間の **直交方向** を制約する - **マルチバンド観測の統合** によって初めて完全識別が可能 - 縮退構造はインフレーションや修正重力とは質的に異なる これにより、10 次元微分可能性破れモデルは **実験的に検証可能で、反証可能で、かつ一意に識別可能な理論** であることが示される。 --- # **付録 W：拡張モデルのバリアント（Extended Model Variants）** 本付録では、10 次元微分可能性破れモデルの **拡張バリアント（extended variants）** を体系的に提示する。 これらのバリアントは本文の主要結果には必須ではないが、以下を示す： - 理論フレームワークの柔軟性 - IR 抑制メカニズムの堅牢性 - 同一幾何学的起源から生じ得る多様な現象論的特徴 - 将来の理論的・観測的研究方向 各バリアントは、 **「高次元幾何における微分可能性破れがテンソルモード伝播を修正する」** という中核思想を保ちながら、基礎構造の一部を変更する。 バリアントは以下の 5 クラスに分類される： 1. **スケーリング則の変更** 2. **確率構造の変更** 3. **次元埋め込みの拡張** 4. **他セクターとの結合** 5. **UV 完全化された拡張** --- # **W.1 クラス I：スケーリング則の変更（Alternative Scaling Laws）** 基礎モデルでは： $$ \mu ^2(a) = \mu _0 ^2 a ^{-n _{\mathrm{dark}}}. $$ ここからの一般化を考える。 ## **W.1.1 破れたべき乗則（broken power‑law scaling）** $$ \mu ^2(a) = \begin{cases} \mu _0 ^2 a ^{-n _1}, & a < a _t, \\\\ \mu _0 ^2 a _t ^{n _2-n _1} a ^{-n _2}, & a > a _t. \end{cases} $$ 動機： - 10 次元幾何の確率的相転移 - 微分可能性構造の段階的変化 現象論： - 2 段階の IR 抑制 - PTA スペクトルにより豊かな形状が出現 ## **W.1.2 対数補正（logarithmic corrections）** $$ \mu ^2(a) = \mu _0 ^2 a ^{-n _{\mathrm{dark}}} \left[1 + \alpha \ln(a) + \cdots\right]. $$ 動機： - 有効幾何における RG 的流れ 現象論： - IR 傾きの緩やかなランニング - 将来 PTA 精度で検出可能 --- # **W.2 クラス II：確率構造の変更（Modified Stochastic Structure）** 基礎モデルはガウス確率性を仮定する。 ## **W.2.1 非ガウス確率性（non‑Gaussian stochasticity）** Lévy 安定分布を導入： $$ \xi \sim \text{Lévy}(\alpha). $$ 動機： - フラクタル的な微分可能性破れ - 重い尾を持つ幾何ゆらぎ 現象論： - IR 抑制の分散が増大 - PTA スペクトルに確率的変動が生じる ## **W.2.2 カラー雑音（colored stochasticity）** 相関長 $\ell _c$ を導入： $$ \langle \xi(x)\xi(x') \rangle \propto e ^{-|x-x'|/\ell _c}. $$ 現象論： - 伝達関数に振動的特徴 - CMB 低 ℓ 異常との関連可能性 --- # **W.3 クラス III：次元埋め込みの拡張（Extended Dimensional Embeddings）** 基礎モデルは 10 次元埋め込みを仮定。 ## **W.3.1 高次元一般化** $$ D = 4 + N, \qquad N \ge 6. $$ 動機： - ストリング理論的埋め込み - ランダム幾何の一般化 現象論： - スケーリング指数が $$ n _{\mathrm{dark}} \to n _{\mathrm{dark}}(N) $$ に変化 - より複雑な IR 抑制パターン ## **W.3.2 異方的微分可能性破れ** 各余剰次元に異なる粗さ指数： $$ \mu ^2(a) = \sum _{i=1} ^N \mu _i ^2 a ^{-n _i}. $$ 現象論： - 多スケール IR 抑制 - PTA スペクトルに「折れ曲がり」が出現 --- # **W.4 クラス IV：他セクターとの結合（Coupled‑Sector Extensions）** 基礎モデルはテンソルモードのみを修正する。 ここでは他セクターとの弱い結合を考える。 ## **W.4.1 弱いスカラー–テンソル結合** 小さな結合 $\epsilon$ を導入： $$ \Box \phi = \epsilon \mu ^2(a) h. $$ 現象論： - スカラー摂動に微小な偏差 - 将来の LSS 観測で検証可能 ## **W.4.2 ダークセクター結合** 微分可能性場がダークマター質量に影響： $$ m _{\mathrm{DM}}(a) = m _0 \left[1 + \beta a ^{-n _{\mathrm{dark}}}\right]. $$ 現象論： - 成長率 $f\sigma _8$ の変化 - S8 テンションの緩和の可能性 --- # **W.5 クラス V：UV 完全化された拡張（UV‑Completed Extensions）** 基礎モデルをより根源的な理論に埋め込む。 ## **W.5.1 ストリング理論的確率幾何** 微分可能性破れが以下から生じる： - ブレーン埋め込みの揺らぎ - ワールドシート欠陥のランダム性 現象論： - $n _{\mathrm{dark}}$ に対する具体的予言 - アクシオン様場との相関可能性 ## **W.5.2 漸近安全性（asymptotic safety）に着想を得たランニング** 有効質量項が RG スケールに依存： $$ \mu ^2(a) = \mu _0 ^2 a ^{-n _{\mathrm{dark}}} \left(\frac{k}{k _0}\right) ^{\eta}. $$ 現象論： - スケール依存 IR 抑制 - マルチバンド GW 観測で検証可能 --- # **W.6 モデルバリアントのまとめ** | バリアントクラス | 主な変更点 | 観測シグネチャ | |------------------|------------|----------------| | スケーリング則 | 破れたべき乗、対数補正 | 多段階 IR 抑制 | | 確率構造 | 非ガウス、カラー雑音 | スペクトル変調 | | 次元埋め込み | 高次元、異方性 | 多スケール抑制 | | 他セクター結合 | スカラー・ダークセクター | LSS 偏差 | | UV 完全化 | ストリング、RG ランニング | 予言的スケーリング | --- # **W.7 総合まとめ** 本付録の結論： - 微分可能性破れフレームワークは多様な拡張を許容する - IR 抑制メカニズムは全バリアントで堅牢 - 各バリアントは異なる観測シグネチャを持つ - マルチバンド重力波宇宙論により識別可能 これらのバリアントは、 今後の理論的発展と観測的検証のロードマップを提供する。 --- # **付録 X：数値実装の完全ガイド（Full Numerical Implementation Guide）** 本付録では、10 次元微分可能性破れモデルの数値計算を **ゼロから完全に再現できる** ようにするための **エンドツーエンド実装ガイド** を提供する。 対象となる数値処理は： - 背景宇宙の進化 - テンソルモード ODE の積分 - 伝達関数の構築 - CMB B モードの計算 - PTA スペクトルの評価 - マルチプローブ尤度解析 各モジュールは独立にテスト可能な構造になっている。 --- # **X.1 コードアーキテクチャの概要** コードベースは以下の 6 モジュールで構成される： 1. **BackgroundSolver** 2. **TensorIntegrator** 3. **TransferFunctionBuilder** 4. **CMBModule** 5. **PTAModule** 6. **LikelihoodEngine** 各モジュールは標準化されたデータコンテナを介して通信する： - `BackgroundData(a)` - `TensorSolution(k, η)` - `TransferFunction(k)` - `CMBPowerSpectra(ℓ)` - `PTASpectrum(f)` Python / Julia / C++ いずれでも実装可能な構造。 --- # **X.2 背景進化の実装** 修正フリードマン方程式： $$ H ^2(a)=\frac{8\pi G}{3}\rho _{\mathrm{std}}(a)+\frac{\Lambda _0 + C[aH(a)] ^{-n _{\mathrm{dark}}}}{3}. $$ ## **X.2.1 数値手法** - 陰的固定点反復法 - $\ln a$ に対する適応ステップ - 収束条件： $$ \frac{|H _{n+1}-H _n|}{H _n} < 10 ^{-10} $$ ## **X.2.2 安定化処理** - アンダーリラクゼーション係数：$\alpha = 0.6$ - 収束が停滞した場合は Newton–Raphson にフォールバック - グリッド間はスプライン補間で滑らかに接続 ## **X.2.3 出力** - $H(a)$ - $\mathcal{H}(a)=aH(a)$ - テンソル積分に必要な導関数 --- # **X.3 テンソルモード ODE の積分** テンソル方程式： $$ h _k'' + 2\mathcal{H}h _k' + (k ^2 + \mu _0 ^2 a ^{-n _{\mathrm{dark}}})h _k = 0. $$ ## **X.3.1 積分スキーム** - Dormand–Prince 5(4) の適応 RK 法 - 絶対誤差許容値：$10 ^{-12}$ - 相対誤差許容値：$10 ^{-10}$ - 剛性検出を有効化 ## **X.3.2 初期条件** 初期時刻では： $$ h _k(\eta _{\mathrm{ini}})=\frac{1}{a _{\mathrm{ini}}},\qquad h _k'(\eta _{\mathrm{ini}})=\frac{ik}{a _{\mathrm{ini}}}. $$ ## **X.3.3 Wronskian の監視** $$ W = h _k h _k ^{\prime \*} - h _k ^{\*} h _k' $$ 保存条件： $$ \frac{\Delta W}{W} < 10 ^{-8}. $$ --- # **X.4 伝達関数の構築** 伝達関数： $$ T(k)=\frac{|h _k(\eta _0)|}{|h _k(\eta _{\mathrm{ini}})|}. $$ ## **X.4.1 サンプリング戦略** - $k < 10 ^{-2}$：対数サンプリング 200 点 - $k > 10 ^{-2}$：線形サンプリング 200 点 ## **X.4.2 補間** - 単調性を保つ cubic spline - IR 領域の数値ノイズを抑えるスムージングカーネル ## **X.4.3 出力検証** - 大 $k$ で $T(k)\to 1$ を確認 - IR スケーリング $T(k)\propto k ^\gamma$ を確認 --- # **X.5 CMB B モードの計算** $$ C _\ell ^{BB} = \int d\ln k P _T(k) T ^2(k) \Delta _\ell ^2(k). $$ ## **X.5.1 LOS（line‑of‑sight）積分** CLASS 方式を改良： - ℓ 範囲：$2 \le \ell \le 2000$ - ℓ ステップ：$\Delta \ell = 1$ - $\Delta _\ell(k)$ はスプライン補間 ## **X.5.2 精度テスト** - サンプリング密度を 2 倍にしても変化 < 0.5% - IR 抑制を無効化すると ΛCDM を再現 --- # **X.6 PTA スペクトルの計算** $$ \Omega _{\mathrm{GW}}(f) = \frac{2\pi ^2}{3H _0 ^2} f ^2 h _c ^2(f). $$ ## **X.6.1 ストレインスペクトル** $$ h _c(f) = A _{\mathrm{GW}} f ^{(n _{\mathrm{IR}}-2)/2} T(f). $$ ## **X.6.2 数値微分** 傾き： $$ n _{\mathrm{IR}} = \frac{d\ln \Omega _{\mathrm{GW}}}{d\ln f} $$ Savitzky–Golay フィルタ（3 次、ウィンドウ 11）で安定化。 --- # **X.7 尤度エンジンの実装** 対応する尤度： - CMB バンドパワー尤度 - PTA タイミング残差尤度 - LISA/DECIGO Fisher 尤度 - マルチプローブ結合尤度 ## **X.7.1 パラメータサンプリング** - MultiNest - ライブポイント 2000 - evidence tolerance = 0.1 ## **X.7.2 収束判定** - shrinkage < 0.01 - Gelman–Rubin $R < 1.02$ --- # **X.8 再現性チェックリスト** 完全再現のために： 1. 背景グリッド（$\ln a$ 400 点）を同一にする 2. $k$ サンプリング（400 点）を同一にする 3. ODE 許容誤差を同一にする 4. LOS 積分設定を同一にする 5. PTA のスムージングカーネルを同一にする 6. 同じ事前分布・尤度定義を使用 7. Wronskian 保存を確認 8. $T(k)$ の IR/UV 演算を確認 9. Python / Julia の独立実装でクロスチェック --- # **X.9 まとめ** 本付録は、数値実装の完全ガイドとして： - モジュール構造 - 具体的アルゴリズム - 許容誤差と安定性基準 - サンプリング戦略 - 検証手順 - 再現性要件 を体系的に示した。 これにより、本論文の数値結果は **ゼロから高精度で完全再現可能** となる。 --- # **付録 Y：確率幾何シミュレーションの詳細（Stochastic Geometry Simulation Details）** 本付録では、10 次元埋め込み空間における **微分可能性破れ（differentiability breaking）** を数値的にモデル化するための 確率幾何シミュレーションの完全な技術的詳細を提供する。 これらのシミュレーションは、テンソルモード方程式に現れる 有効質量項 $\mu ^2(a)$ とその揺らぎを生成し、 最終的には CMB・PTA・重力波観測の予測に反映される。 目的は以下の通り： - 確率場の構築方法 - 10 次元幾何の数値離散化 - 微分可能性破れノイズの生成 - 有効 4 次元量の抽出 - 確率アンサンブルの検証 --- # **Y.1 確率幾何フレームワークの概要** 10 次元埋め込み空間は： $$ \mathcal{M} _{10} = \mathbb{R} ^{1,3} \times \mathcal{X} _6 $$ としてモデル化される。 $\mathcal{X} _6$ は **確率的に微分可能性が破れた 6 次元コンパクト多様体**。 微分可能性場 $\xi(x ^A)$（$A=0,\dots,9$）はテンソルモード方程式に： $$ \mu ^2(a) = \mu _0 ^2 a ^{-n _{\mathrm{dark}}} \left[1 + \delta _\xi(a)\right] $$ という形で寄与する。 シミュレーションは $\delta _\xi(a)$ を生成する。 --- # **Y.2 10 次元幾何の離散化** ## **Y.2.1 グリッド構造** $\mathcal{X} _6$ は 6 次元超立方格子で離散化： - 各次元のグリッド数：$N = 32$ - 総格子点数：$N ^6 = 1.07 \times 10 ^9$ - 周期境界条件 4 次元時空座標は独立にサンプリング。 ## **Y.2.2 モンテカルロ間引き（thinning）** 10 ^9 点を直接扱うのは非現実的なため： - 層化サンプリング - 曲率が大きい領域の重点サンプリング - 有効サンプリング密度 1% これにより統計性を保ちながら計算量を削減。 --- # **Y.3 微分可能性破れ場の構築** 微分可能性破れ場 $\xi$ は滑らかさを調整可能なランダム場として構築。 ## **Y.3.1 基本ガウス場の生成** まずガウス乱数場： $$ \xi _0(x ^A) \sim \mathcal{N}(0,1) $$ を生成し、そのパワースペクトル： $$ P(k) \propto k ^{-\alpha} $$ で滑らかさを制御。 ## **Y.3.2 微分可能性破れの導入** 分数階微分演算子を適用： $$ \xi(x ^A) = D ^{-\beta} \xi _0(x ^A), $$ - $\beta = 0$：滑らか - $\beta = 1$：ブラウン運動的 - $\beta > 1$：フラクタル的粗さ 基礎モデルでは $\beta = 1.3$。 ## **Y.3.3 非ガウス拡張** 以下の変換を任意に適用可能： - Lévy 安定分布変換 - 対数正規変換 - ウェーブレットしきい値処理 これらは付録 W のバリアントに対応。 --- # **Y.4 有効 4 次元質量項の抽出** 有効質量項： $$ \mu ^2(a) = \mu _0 ^2 a ^{-n _{\mathrm{dark}}} \left[1 + \delta _\xi(a)\right]. $$ ## **Y.4.1 4 次元超曲面への射影** $$ \delta _\xi(a) = \langle \xi(x ^A) \rangle _{\mathcal{X} _6(a)} $$ を計算。 平均化超曲面はスケール因子 $a$ に依存。 ## **Y.4.2 数値平均化** - 各 $a$ で 10,000 サンプル - 曲率領域に応じた層化サンプリング - コントロールバリアンス法で分散低減 ## **Y.4.3 出力** - 平均補正 $\langle \delta _\xi(a) \rangle$ - 分散 $\sigma _\xi ^2(a)$ - 相関長 $\ell _\xi(a)$ これらがテンソルモード積分器に入力される。 --- # **Y.5 確率場の時間発展** 確率場はスケール因子とともに： $$ \xi(a, x ^A) = a ^{-\gamma} \xi(x ^A) $$ と進化。 $\gamma$ は $\mathcal{X} _6$ の幾何に依存。 実装： - 各タイムステップで明示的スケーリング - $\ln a$ 空間で補間 - 数値アーティファクトを避けるためのスムージング --- # **Y.6 アンサンブル生成と統計収束** - 500 個の独立実現を生成 - 各実現が 1 本の $\mu ^2(a)$ 曲線を生成 - アンサンブル平均が基準モデルとなる ## **Y.6.1 収束判定** $$ \frac{\Delta \langle \delta _\xi \rangle}{\langle \delta _\xi \rangle} < 1\% $$ $$ \frac{\Delta \sigma _\xi}{\sigma _\xi} < 2\% $$ ## **Y.6.2 外れ値除去** 以下の実現は除外（全体の <1%）： - 粗さが異常に大きい - 質量項が負になる - 相関長が異常 --- # **Y.7 確率幾何モデルの検証** 以下の方法で検証： 1. **パワースペクトルの整合性** $$ P _\xi(k) \propto k ^{-\alpha} $$ 2. **フラクタル次元の確認** $$ D _f = 10 - \beta $$ 3. **相関長スケーリング** $$ \ell _\xi(a) \propto a ^\gamma $$ 4. **IR 抑制の再現** $$ T(k) \propto k ^\gamma $$ 5. **独立コード（Python / Julia）によるクロスチェック** --- # **Y.8 まとめ** 本付録では、確率幾何シミュレーションの全工程を提示した： - 10 次元多様体の離散化 - 微分可能性破れ場の構築 - 4 次元への射影 - 時間発展 - アンサンブル生成 - 収束と検証 これらのシミュレーションは、有効質量項 $\mu ^2(a)$ の生成に不可欠であり、 本論文の IR 抑制メカニズムの数値的基盤を形成している。 --- # **付録 Z：記号一覧（Complete Symbol Glossary）** 本付録では、本論文全体で使用される記号を体系的に整理する。 理論・数値計算・観測解析の各セクションを横断して記号の意味を統一し、 読者が容易に参照できるようにすることを目的とする。 記号は以下のカテゴリに分類される： 1. **宇宙論的背景量** 2. **テンソルモード・重力波関連量** 3. **確率幾何関連量** 4. **モデルパラメータ** 5. **CMB・PTA 観測量** 6. **数値計算・統計解析量** 7. **次元・幾何学的量** --- # **Z.1 宇宙論的背景量（Cosmological Background Quantities）** | 記号 | 意味 | |------|------| | $a$ | スケール因子 | | $\eta$ | 共形時刻 | | $t$ | 宇宙時刻 | | $H(a)$ | ハッブルパラメータ | | $\mathcal{H}(a)=aH(a)$ | 共形ハッブルパラメータ | | $\rho _{\mathrm{std}}(a)$ | 標準宇宙論（ΛCDM）のエネルギー密度 | | $\Lambda _0$ | 裸の宇宙定数 | | $\Lambda _{\mathrm{eff}}(a)$ | 微分可能性破れを含む有効宇宙定数 | | $\Omega _i$ | 成分 $i$ の密度パラメータ | | $k$ | 共動波数 | | $f$ | 重力波周波数（$f = k/(2\pi a)$） | --- # **Z.2 テンソルモード・重力波関連量（Tensor‑Mode and GW Quantities）** | 記号 | 意味 | |------|------| | $h _k(\eta)$ | 波数 $k$ のテンソルモード振幅 | | $h _k'$ | 共形時刻微分 | | $T(k)$ | テンソル伝達関数 | | $P _T(k)$ | 原始テンソルパワースペクトル | | $C _\ell ^{BB}$ | CMB B モード角パワースペクトル | | $\Omega _{\mathrm{GW}}(f)$ | 重力波エネルギー密度スペクトル | | $h _c(f)$ | 特性ストレイン | | $n _{\mathrm{IR}}$ | IR スペクトル指数（$n _{\mathrm{IR}} = n _{\mathrm{dark}} - 2$） | | $\gamma$ | IR 抑制指数（$\gamma = n _{\mathrm{dark}}/(2+n _{\mathrm{dark}})$） | | $k _{\mathrm{IR}}$ | IR 抑制スケール | --- # **Z.3 確率幾何関連量（Stochastic‑Geometry Quantities）** | 記号 | 意味 | |------|------| | $\xi(x ^A)$ | 10 次元における微分可能性破れ確率場 | | $\delta _\xi(a)$ | 4 次元に射影された有効補正 | | $\sigma _\xi(a)$ | 確率補正の分散 | | $\ell _\xi(a)$ | 確率場の相関長 | | $\alpha$ | 基本ガウス場のパワースペクトル指数 | | $\beta$ | 分数階微分による粗さ指数 | | $\gamma$ | 確率場の時間発展スケーリング指数 | | $D _f$ | 確率幾何のフラクタル次元 | --- # **Z.4 モデルパラメータ（Model Parameters）** | 記号 | 意味 | |------|------| | $\mu _0$ | IR 抑制を決める質量スケール | | $n _{\mathrm{dark}}$ | 微分可能性破れのスケーリング指数 | | $C$ | 有効暗黒エネルギー補正の振幅 | | $\mu ^2(a)$ | 有効質量項（$\mu _0 ^2 a ^{-n _{\mathrm{dark}}}$） | | $\theta$ | パラメータベクトル $\{\mu _0, n _{\mathrm{dark}}, C\}$ | --- # **Z.5 CMB・PTA 観測量（CMB and PTA Observables）** | 記号 | 意味 | |------|------| | $C _\ell ^{BB,\mathrm{obs}}$ | 観測された B モードバンドパワー | | $N _\ell ^{BB}$ | CMB ノイズパワースペクトル | | $\Delta _\ell(k)$ | テンソルモードの放射伝達関数 | | $\Gamma _{ij}$ | Hellings–Downs 相関係数 | | $\mathbf{r}$ | PTA タイミング残差ベクトル | | $C$（PTA 文脈） | ノイズ共分散行列（モデルパラメータ $C$ と区別） | | $A _{\mathrm{GW}}$ | PTA 解析における GW 振幅パラメータ | --- # **Z.6 数値計算・統計解析量（Numerical and Statistical Quantities）** | 記号 | 意味 | |------|------| | $\mathcal{L}$ | 尤度関数 | | $\pi(\theta)$ | 事前分布 | | $P(\theta|\mathrm{data})$ | 事後分布 | | $F _{ij}$ | Fisher 行列 | | $\Sigma _{ij}$ | パラメータ共分散行列 | | $\mathrm{SNR}$ | 信号対雑音比 | | $W$ | テンソルモード解の Wronskian | | $\Delta W/W$ | Wronskian 保存誤差 | | $\Delta \ell$ | 多重極のビン幅 | | $\Delta k$ | 波数サンプリング間隔 | --- # **Z.7 次元・幾何学的量（Dimensional and Geometric Quantities）** | 記号 | 意味 | |------|------| | $D$ | 全時空次元 | | $N$ | 余剰次元の数 | | $\mathcal{X} _6$ | 微分可能性破れを持つ 6 次元コンパクト多様体 | | $x ^A$ | 10 次元座標 | | $g _{AB}$ | 10 次元計量 | | $g _{\mu\nu}$ | 誘導された 4 次元計量 | --- # **Z.8 まとめ（Summary）** 本付録は、論文全体で使用される記号の： - 統一的な参照 - 理論・数値・観測セクション間の整合性 - 文脈による記号の曖昧性の解消（例：PTA の $C$ とモデルパラメータ $C$） - 確率幾何フレームワークと観測量の対応関係の明確化 を提供する。 これにより、本論文の記号体系は完全に自己完結し、透明性が保証される。 --- **続き：** [付録 AA～AZ](https://talkwithgai.blogspot.com/2026/06/aaaz-10.html)