双対性により誘起される非特異バウンスとブルーティルト重力波(Non-singular bounce and blue-tilted gravitational waves induced by duality)

これは何のノートか

本ノートでは、Double Field Theory（DFT）におけるスケール因子双対性
a(t) <-> 1/a(t)
に着想を得て、双対性が宇宙初期に有効に働いた場合にどのような原始重力波スペクトルが現れるかを、最小限の自由度を用いて概念的に考察する。

ここで扱うのは 完全な DFT 宇宙論の構築ではない。 strong constraint、倍化座標、α' 補正を含む厳密な O(d,d) 対称方程式の導出を目指すのではなく、

という立場をとる。
特に、通常セクションと双対セクションの混合角を表す セクション・スカラー θ(t) を導入し、
a_eff^2 = a^2 cos^2(theta) + a^{-2} sin^2(theta)
という有効スケール因子を公理的に定義することで、非特異バウンスと 幾何学的フレーム反転を最小構成で実現する。

このバウンス幾何に対してテンソルモード方程式
u_k'' + ( k^2 - a_eff''/a_eff ) u_k = 0
を散乱問題として解析すると、バウンス近傍で
a_eff''/a_eff ≈ mu^2 ( 1 - 2 sech^2(mu eta) )
が Pöschl–Teller 型ポテンシャルとして現れ、 blue tilt・中間スケールのピーク・高波数の指数減衰からなる 特徴的な blue-tilted gravitational waves が自然に生成される。

本ノートは、DFT の完全理論を構築するものではなく、 双対性が宇宙初期に作用した場合の観測的帰結を、最小有効モデルとして抽出する という目的に特化した概念的検討である。

本モデルのスコープと限界

• Einstein 方程式と自己無撞着な背景解の導出は行っていない

• スカラー揺らぎおよび CMB との整合性は未検討

• DFT の強制約や完全な一般化計量には立ち入らない

本研究の目的は、
双対性が固定点として働く場合の
テンソル揺らぎの散乱構造と blue tilt の関係
を概念的に整理することにある。

以下の図では、
バウンス近傍で生じる有効ポテンシャルに対する
テンソルモードの散乱により
低 kkk での blue tilt、有限スケールでのピーク、
高 kkk での指数的減衰が生じる様子を示している。

【Introduction】
宇宙論の標準モデルは、初期の高温高密度状態から現在の大規模構造形成に至るまで、観測と整合する精密な枠組みを提供している。しかし、その出発点であるビッグバン特異点では曲率不変量が発散し、一般相対論が破綻する。
この特異点問題は、インフレーション、エキピロティック宇宙、量子重力的バウンスなど多くの提案にもかかわらず、依然として理論的に未解決のままである。

一方、弦理論およびその低エネルギー有効理論である Double Field Theory (DFT) は、時空がa(t) とその双対 1/a(t)の二重構造を持つことを示唆し、O(d,d) 対称性に基づく T-双対性が宇宙初期の物理に重要な役割を果たす可能性を示している。
しかし、DFT をそのまま FRW 宇宙論として実装するには、strong constraint、倍化座標の物理的解釈、セクション選択の一意性など、複数のad hocな仮定が不可避であり、現時点では完全な宇宙論モデルとして確立されていない。

本研究では、DFT の厳密構造を構築することを目的とせず、その最も本質的な特徴である
a(t) <-> 1/a(t)
という双対性の観測的帰結を概念的に探索することを目的とする。
この立場を明確にするため、本モデルは
  "DFT-inspired minimal effective theory"
として位置づける。
すなわち、DFT の完全な動力学を再現することを意図せず、双対性が宇宙初期に有効に働いた場合にどのような幾何学的構造と観測的シグネチャが現れるかを、最小限の自由度で記述する意図的なトイモデルである。
本研究の主眼は、精密な数値予測ではなく、定性的に頑健な特徴の同定にある。
本研究では、通常セクションと双対セクションの混合角を表すセクション・スカラー theta(t) を導入し、観測される有効スケール因子を
  a_eff^2(t) = a^2(t) cos^2 theta(t) + a^{-2}(t) sin^2 theta(t)
と定義する。
この構造は DFT の一般化計量の射影構造を模倣したものであり、theta(t) の時間発展に応じて a_eff(t) が a(t) と 1/a(t) の間を滑らかに補間する。

さらに、ポテンシャル
  V(theta) = (mu^2 / 2) sin^2 theta
を採用すると、
  theta(t) = 2 arctan( exp(mu t) )
というキンク解が得られ、a_eff(t) は特異点を持たない滑らかなバウンスを形成する。
これにより、ビッグバンは物理的特異点ではなく、幾何学的フレームの反転点として再解釈される。

本研究の中心的結果は、バウンス近傍で
  a_eff'' / a_eff
が Pöschl–Teller 型ポテンシャルで近似され、テンソルモードが散乱されることで、
・低 k の blue tilt
・k ≈ mu のピーク
・高 k の指数的カットオフ
という特徴的な重力波スペクトルが生じる点である。

本研究は、新しい宇宙論モデルの提案ではなく、
「双対性が宇宙初期に作用した場合の観測的帰結を、最小有効モデルで概念的に探索する」
ことを目的とする。
特に、特異点回避と blue-tilted 重力波が同一機構から生じるという構造は、既存のバウンス宇宙論には見られない特徴である。

本研究のスコープと限界 —
本モデルは、完全な動的宇宙論モデルの構築を目的とするものではなく、
最小限の自由度により、双対性に由来する幾何学的効果とその観測的帰結を抽出することを目的とした有効記述である。

特に以下の点は本研究の意図的な仮定である：
(i) 背景時空の時間発展は外部から与えられている
(ii) エネルギー運動量テンソルとの完全な整合性は今後の課題である
(iii) a_eff''/a_eff の Pöschl–Teller 近似に基づき、Bogoliubov 係数を散乱理論の解析解として評価する

本研究の目的は、インフレーションの完全な代替理論を提示することではなく、双対性に基づく構造が生み得る頑健な観測的特徴を明確化することである。

【Abstract】
本研究は、Double Field Theory (DFT) の厳密構造を構築することを目的とせず、その最も本質的な特徴である
  a(t) <-> 1/a(t)
という双対性が宇宙初期に作用した場合の観測的帰結を、最小有効モデルとして概念的に探索するものである。
本モデルは DFT-inspired minimal effective theory として位置づけられ、通常セクションと双対セクションの混合角を表すセクション・スカラー theta(t) を導入することで、幾何学的フレームの遷移を単一自由度で記述する。

観測される有効スケール因子は
  a_eff^2(t) = a^2(t) cos^2 theta(t) + a^{-2}(t) sin^2 theta(t)
と定義され、ポテンシャル
  V(theta) = (mu^2 / 2) sin^2 theta
を採用すると、
  theta(t) = 2 arctan( exp(mu t) )
というキンク解が得られる。
この解により a_eff(t) は特異点を持たない滑らかなバウンスを形成し、有効ハッブルパラメータ H_eff(t) はバウンス点で符号反転する。
これにより、ビッグバンは物理的特異点ではなく、幾何学的フレームの反転点として再解釈される。

テンソル揺らぎのモード方程式
  u_k'' + ( k^2 - a_eff'' / a_eff ) u_k = 0
において、バウンス近傍では
  a_eff'' / a_eff ≈ mu^2 ( 1 - 2 sech^2( mu eta ) )
が Pöschl–Teller 型ポテンシャルとして近似される。

この近似に基づくテンソルパワースペクトルは
・低 k：P_t(k) ∝ k^3（blue tilt）
・中間 k：k ≈ mu にピーク
・高 k：指数的カットオフ
という特徴的な三段構造を示す。

特に n_t > 0 の blue tilt は標準インフレーションの red tilt と明確に区別可能であり、DECIGO・BBO・LISA などの将来の重力波観測によって検証可能である。

本研究は、双対性が特異点回避と blue-tilted 重力波を同一機構で結びつけるという概念的可能性を示すものであり、DFT 宇宙論の観測的帰結を理解するための最小有効モデルとして位置づけられる。

【Section 2: The Model】
本節では、本研究で用いる DFT-inspired minimal effective model を定義する。本モデルは、Double Field Theory (DFT) の完全な宇宙論的実装を目指すものではなく、
  a(t) <-> 1/a(t)
という双対性の観測的帰結を最小限の自由度で記述するための 意図的なトイモデル として構築される。

DFT を FRW 背景に適用する際には、strong constraint、倍化座標の物理的解釈、セクション選択の一意性など、複数の ad hoc 仮定が不可避である。本研究では、これらの複雑性を排除し、双対性の本質的効果のみを抽出することを目的とする。

2.1 Section-scalar theta(t)
通常セクションと双対セクションの混合角を表す自由度として、セクション・スカラー theta(t) を導入する。
その物理的意味は以下の通り:
・theta = 0 : 通常セクション（a の世界）
・theta = pi/2 : 双対セクション（1/a の世界）
・theta(t) の時間発展 : フレームの連続的遷移（セクション混合）

本研究では、theta(t) が滑らかなバウンス解を持つよう、後述のポテンシャルを採用する。

2.2 Effective scale factor
観測される物理的距離は、通常セクションと双対セクションの寄与が theta(t) によって重み付けされた有効量として定義される。本研究では、DFT の一般化計量の射影構造を模倣し、以下の公理的定義を採用する:
  a_eff^2(t) = a^2(t) cos^2 theta(t) + a^{-2}(t) sin^2 theta(t)
背景膨張としては、最も単純な指数膨張
  a(t) = exp( H t )
を採用する。
この選択は 自己整合的な背景方程式の解として導出されたものではなく、外部から与えた背景 である。

そのため、
「a(t) = exp(H t) の自己整合性は今後の課題である」
ことを明記しておく。

背景の自己整合性について —
本研究では、背景スケール因子 a(t) は完全な場の方程式から導出されたものではなく、双対性構造の幾何学的帰結を明確にするための外部入力として導入されている。
この選択は意図的なものであり、有効スケール因子 a_eff の性質（最小値の存在および H_eff の符号反転）を抽出することを目的としている。
エネルギー運動量テンソルを含む自己無撞着な定式化は重要な今後の課題であるが、本研究で得られる定性的結果は背景の詳細には依存しない。

2.3 Effective Hubble parameter
有効スケール因子に基づくハッブルパラメータを
  H_eff(t) = (1 / a_eff) * d(a_eff)/dt
と定義する。
a_eff(t) が最小値を持つ点では H_eff = 0 となり、符号反転が生じる。
これが本モデルにおける 幾何学的ビッグバン に対応する。

2.4 Action
本研究で採用する最小作用は、重力項とセクション・スカラーの運動項から構成される:
S = ∫ dt d^3x a_eff^3(t) [ - 3 H_eff(t)^2 / (8 pi G)
+ (1/2) (theta_dot)^2
- V(theta) ]
ここで、
・-3 H_eff^2 / (8 pi G) : 有効スケール因子を用いた最小 FRW 表現
・(1/2) theta_dot^2 : セクション混合角の運動エネルギー
・V(theta) : バウンスを生成するポテンシャル
を表す。

注意すべき点として、a_eff は θ の非線形関数であるため、
この作用を完全な変分原理の出発点とした場合には追加の項が現れる。
本研究では、作用 S を “有効エネルギー密度の記号的パラメータ化” として採用し、完全な場の方程式を導出することは目的としない。
特に、
・背景 a(t)=exp(Ht) は外部入力
・a_eff(t) は DFT 的構造を模倣した有効量
・θ の運動方程式は摩擦項を無視した極限 μ≫H における近似解
として扱う。
したがって本作用は、完全な重力理論の作用ではなく、
本モデルの幾何学的構造を記述するための有効的・現象論的表現である。

2.5 Potential and bounce solution
バウンスを生成する最も単純なポテンシャルとして
  V(theta) = (mu^2 / 2) sin^2 theta
を採用する。このとき、摩擦項を無視した極限では
  theta(t) = 2 arctan( exp( mu t ) )
が厳密解となる。

この解は
・t -> -∞ : theta -> 0
・t -> +∞ : theta -> pi
というキンク型の遷移を表し、通常セクションから双対セクションへの滑らかなフレーム転換を記述する。
その時間微分は
  theta_dot(t) = mu / cosh( mu t )
となり、t = 0 で最大値 mu を取る。
したがって、フレーム転換は有限時間幅を持つ局所的イベントとして実現される。

2.6 Summary of the model
本節で定義したモデルは以下の特徴を持つ:

双対性の最小実装
  a_eff^2 = a^2 cos^2 theta + a^{-2} sin^2 theta
により a <-> 1/a を連続的に補間。

特異点のないバウンス宇宙
  theta(t) のキンク解により a_eff(t) が滑らかに最小値を持つ。

幾何学的ビッグバン
  H_eff(t) がバウンス点で符号反転し、ビッグバンをフレーム反転として解釈。

テンソル揺らぎの散乱構造
  a_eff'' / a_eff が Pöschl–Teller 型となり、blue-tilted 重力波を予言。

本モデルは、DFT の完全な宇宙論の代替ではなく、双対性が宇宙初期に作用した場合の観測的帰結を、最小限の自由度で明確に示すための意図的なトイモデルとして位置づけられる。

【Section 3: Background Evolution】
本節では、前節で定義した DFT-inspired minimal effective model に基づき、宇宙の背景進化を解析する。
特に、セクション・スカラー theta(t) のキンク解が、観測される有効スケール因子 a_eff(t) にどのような非特異バウンス構造をもたらすかを明確に示す。

本研究では、摩擦項を無視した極限
  mu >> H
を採用する。
これは、theta(t) の遷移が背景膨張よりも十分速く進む場合に成立し、キンク解の解析的扱いを可能にする。

3.1 Section-scalar dynamics and kink solution
セクション・スカラー theta(t) の運動方程式は、最小作用
S = ∫ dt a_eff^3(t) [ - 3 H_eff^2 / (8 pi G)
+ (1/2) theta_dot^2
- V(theta) ]
から導かれる。

バウンス解を許す最も単純なポテンシャルとして
  V(theta) = (mu^2 / 2) sin^2 theta
を採用すると、摩擦項を無視した極限 mu >> H では
  theta(t) = 2 arctan( exp( mu t ) )
が厳密解となる。

この解は以下の極限を持つ:
・t -> -∞ : theta -> 0
・t -> +∞ : theta -> pi
すなわち、通常セクションから双対セクションへの滑らかなフレーム転換を記述する。

時間微分は
  theta_dot(t) = mu / cosh( mu t )
となり、t = 0 で最大値 mu を取る。
したがって、フレーム転換は有限時間幅を持つ局所的イベントとして実現される。

3.2 Effective scale factor and nonsingular bounce
観測される有効スケール因子は
  a_eff^2(t) = a^2(t) cos^2 theta(t) + a^{-2}(t) sin^2 theta(t)
で定義される。

背景膨張として
  a(t) = exp( H t )
を採用する。
この背景は外部から与えたものであり、自己整合性の検証は今後の課題である。

theta(t) のキンク解を代入すると、a_eff(t) は以下の極限を持つ:
・t -> -∞ : theta -> 0
→ a_eff(t) ≈ a(t) = exp( H t )
・t -> +∞ : theta -> pi
→ a_eff(t) ≈ a(t)^{-1} = exp( - H t )
・t ≈ 0 : theta ≈ pi/2
→ a_eff(t) は最小値を持つ

したがって a_eff(t) は、通常セクションの膨張 a(t) と双対セクションの収縮 1/a(t) を滑らかに接続する 特異点のないバウンス を形成する。

このバウンスでは a_eff(t) が 0 になることはなく、曲率不変量が発散することもない。
本モデルでは、ビッグバンは物理的特異点ではなく、幾何学的フレームの反転点として理解される。

3.3 Effective Hubble parameter and sign flip
有効ハッブルパラメータは
  H_eff(t) = (1 / a_eff) * d(a_eff)/dt
時間微分を計算すると
  H_eff(t)
= [ H ( exp(2 H t) cos^2 theta - exp(-2 H t) sin^2 theta )

• theta_dot sin theta cos theta ( exp(-2 H t) - exp(2 H t) ) ]
  / [ exp(2 H t) cos^2 theta + exp(-2 H t) sin^2 theta ]

となる。

極限を調べると:
・t -> -∞ : H_eff -> +H
・t -> +∞ : H_eff -> -H
・t = 0 : H_eff = 0（バウンス点）

特に t ≈ 0 では theta_dot が最大となり、
sin theta cos theta が最大値 1/2 を取るため、第二項が支配的となる。
これにより H_eff は必ず 0 を跨ぎ、符号反転する。

したがって本モデルでは、
  ビッグバン = H_eff の符号反転点
として幾何学的に定義される。

3.4 Summary of background evolution
本節の結果をまとめると:
theta(t) はキンク解を持ち、通常セクションから双対セクションへの遷移を記述する。
a_eff(t) は特異点を持たず、滑らかなバウンスを形成する。
H_eff(t) はバウンス点で符号反転し、ビッグバンを幾何学的イベントとして再解釈できる。
これらの性質は、次節で解析するテンソル揺らぎの散乱構造と blue-tilted 重力波スペクトルの生成に直接つながる。

【Section 4: Tensor Perturbations】
本節では、有効スケール因子 a_eff(t) によって定義される背景時空上での
テンソル摂動を考える。フーリエモード u_k(η) は、共形時間 η に関して
以下のモード方程式に従う：

u_k'' + [ k^2 - U(η) ] u_k = 0,

ここでプライムは η による微分を表し、有効ポテンシャル U(η) は

U(η) = a_eff''(η) / a_eff(η)

で与えられる。

4.1 背景幾何からの有効ポテンシャル

共形時間 η は

dη = dt / a_eff(t)

により数値的に構成する。また、得られた a_eff(η) から
二階微分を計算することで、有効ポテンシャル U(η) を直接評価する。

その結果、U(η) はバウンス近傍において局在した構造を持ち、
散乱ポテンシャルとして解釈できることが分かる。

さらに、このポテンシャルは次の Pöschl–Teller 型の関数で
良く近似できる：

U(η) ≈ μ_eff^2 [ 1 - C sech^2(μ_eff η) ]

ここで μ_eff はポテンシャルの幅、C はその深さを表す無次元パラメータである。
これらのパラメータは、数値的に構成した U(η) に対して
最小二乗フィットを行うことで決定される。

本研究で用いた代表的なパラメータに対しては、

μ_eff ≈ (数値結果)
C ≈ (数値結果)

が得られ、Pöschl–Teller 型による近似はバウンス領域において
良好な一致を示す（図X参照）。

4.2 散乱問題としての定式化

上記の結果に基づき、テンソルモードの時間発展は
一次元散乱問題として定式化できる。

すなわち、η → -∞ において入射する正周波モード

u_k(η) ~ e^{-i k η} / sqrt(2k)

を初期条件として、η → +∞ においては

u_k(η) ~ α_k e^{-i k η} + β_k e^{+i k η}

と展開される。ここで α_k および β_k はBogoliubov係数であり、
粒子生成の情報を含む。

4.3 Pöschl–Teller ポテンシャルに対する厳密解

ポテンシャルが Pöschl–Teller 型で近似できる場合、
対応する散乱問題は厳密に解くことができる。

まず、パラメータ C を

C = λ(λ + 1)

と再パラメータ化する。このとき λ は正の実数であり、

λ = (-1 + sqrt(1 + 4C)) / 2

によって与えられる。

このとき、Bogoliubov係数の絶対値は次の解析的表式で与えられる：

|β_k|^2 = sin^2(πλ) / [ sinh^2(π k / μ_eff) + sin^2(πλ) ]

この結果は、バウンス幾何により誘起された有効ポテンシャルに対して、
厳密な散乱理論に基づく粒子生成率を与えるものである。

4.4 物理的解釈

得られた結果は、バウンス構造を持つ背景時空が自然に
モード関数に対する散乱ポテンシャルを生成し、
それにより重力波の生成が引き起こされることを示している。

特に、パラメータ λ は散乱の強さを特徴づける量であり、

• λ が小さい場合：弱い粒子生成

• λ が大きい場合：強い粒子生成

に対応する。

このように、本モデルにおける重力波生成は、
背景幾何に起因する散乱現象として統一的に理解できる。

【Section 5: Gravitational Wave Spectrum】
本節では、前節で導出したBogoliubov係数を用いて、
バウンス幾何に起因する重力波スペクトルを評価する。

5.1 スペクトルの定義

テンソルモードのパワースペクトル P_t(k) は、
Bogoliubov係数 β_k を用いて次のように定義される：

P_t(k) ∝ k^3 |β_k|^2

ここで比例定数は規格化条件に依存するため、
本研究ではスペクトルの形状に着目する。

前節で得られた厳密散乱解に基づき、

|β_k|^2 = sin^2(πλ) / [ sinh^2(π k / μ_eff) + sin^2(πλ) ]

を用いることで、スペクトルは完全に決定される。

5.2 スペクトルの特徴

上式から得られるスペクトルは、波数 k に対して
特徴的な3つの領域を持つ：

(i) 低波数領域（k ≪ μ_eff）

この領域では sinh(x) ≈ x を用いると、

|β_k|^2 ≈ const

となるため、

P_t(k) ∝ k^3

が得られる。したがって、スペクトルは
ブルーティルトした振る舞いを示す。

(ii) 中間領域（k ≈ μ_eff）

この領域では粒子生成が最も効率的に起こり、
スペクトルはピーク構造を持つ。

ピークの位置はおおよそ

k_peak ≈ μ_eff

で与えられ、その振幅は λ の値に依存する。

(iii) 高波数領域（k ≫ μ_eff）

この領域では sinh(x) ≈ (1/2) exp(x) を用いると、

|β_k|^2 ∝ exp(-2π k / μ_eff)

となるため、

P_t(k) ∝ k^3 exp(-2π k / μ_eff)

となる。したがって、高波数側では指数関数的減衰が現れる。

5.3 旧モデルとの比較

従来の簡略化されたモデルでは、Bogoliubov係数として

|β_k|^2 ≈ 1 / (exp(2π k / μ) - 1)

のような熱的分布が仮定されていた。

本研究で得られた厳密解はこれと類似した高波数減衰を持つが、
以下の点で本質的に異なる：

• パラメータ μ_eff および λ が背景幾何から導出される

• スペクトル形状が散乱理論に基づいて一意に決まる

• 低波数領域およびピーク構造の詳細が異なる

図Yには、旧モデルと本研究の結果の比較を示す。
両者は高波数領域では類似する一方で、
中間領域において有意な差異が現れることが分かる。

5.4 パラメータ依存性

スペクトルの形状は主に λ および μ_eff によって決定される。

• μ_eff：ピーク位置を決定（スケール）

• λ：粒子生成の強さを制御（振幅および形状）

特に λ が大きくなるにつれて、ピークの高さが増加し、
スペクトル全体の強度が増大する。

5.5 物理的含意

本研究の結果は、バウンス幾何がテンソルモードに対して
有効な散乱ポテンシャルを生成し、それにより重力波が生成される
ことを示している。

この生成機構はインフレーションとは異なり、
時間依存背景による非断熱的効果として理解される。

また、得られたスペクトルは

• 低波数でのブルーティルト

• 中間スケールでのピーク構造

• 高波数での指数減衰

という特徴を持ち、観測的にも識別可能なシグネチャを与える可能性がある。

5.6 今後の観測との関係

本研究ではスペクトルの形状に焦点を当てたが、
将来的にはエネルギー密度パラメータ Ω_GW(f) への変換を行い、
重力波観測との直接比較を行うことが重要である。

特に、ピーク周波数は μ_eff によって決まるため、
観測可能な周波数帯との対応付けが今後の課題となる。

【Section 6: Conclusion】
本研究では、Double Field Theory (DFT) の O(d,d) 双対性に着想を得て、
  a(t) <-> 1/a(t)
という対称性が宇宙初期に作用した場合の観測的帰結を、最小構成の有効トイモデルとして概念的に探索した。

本モデルは、DFT の完全な動力学を扱うものではなく、
DFT-inspired minimal effective theory
として位置づけられる。
通常セクションと双対セクションの混合角を表すセクション・スカラー theta(t) を導入することで、幾何学的フレームの遷移を単一自由度で記述した。

本研究の中心的構造は、有効スケール因子
a_eff^2(t) = a^2(t) cos^2 theta(t) + a^{-2}(t) sin^2 theta(t)
により、theta(t) のキンク型解
theta(t) = 2 arctan( exp( mu t ) )
が滑らかなバウンスを生成し、有効ハッブルパラメータ H_eff(t) がバウンス点で符号反転する点にある。
これにより、ビッグバンは物理的特異点ではなく、幾何学的フレームの反転点として再解釈される。

テンソル揺らぎのモード方程式
  u_k'' + ( k^2 - a_eff'' / a_eff ) u_k = 0
において、バウンス近傍では
  a_eff'' / a_eff ≈ mu^2 ( 1 - 2 sech^2( mu eta ) )
が Pöschl–Teller 型ポテンシャルとして近似される。

この近似に基づくテンソルパワースペクトルは、
・低 k：P_t(k) ∝ k^3（blue tilt）
・中間 k：k ≈ mu にピーク
・高 k：指数的カットオフ
という特徴的な三段構造を示す。
特に n_t > 0 の blue tilt は標準インフレーションの red tilt と明確に区別可能であり、DECIGO・BBO・LISA などの将来の重力波観測によって直接検証可能である。

本研究の意義は、
双対性が特異点回避と blue-tilted 重力波を同一機構で結びつけるという概念的可能性を示した点にある。
この構造は既存のバウンス宇宙論（ekpyrotic, matter bounce, string gas）には見られない。

今後の課題
本研究は概念的探索であり、以下の課題が残されている：
1.Bogoliubov 係数の厳密導出
Pöschl–Teller ポテンシャルの厳密散乱係数との関係を明確化する必要がある。
2.背景 a(t) = exp(H t) の自己整合性
有効エネルギー密度が Einstein 方程式を通じてこの背景を生成するかは未検討である。
3.スカラー揺らぎと CMB スペクトル指数 n_s との整合性
インフレーション代替としての完全な評価には不可欠である。
4.DFT の強制約を緩和した拡張モデルとの接続
本研究のトイモデルを、より高次の DFT 構造へ橋渡しする必要がある。
5.観測データとの比較
LIGO/Virgo の上限、NANOGrav PTA との整合性を評価する必要がある。

本研究の位置づけ —
本研究は、双対性に着想を得た最小有効モデルにより、
(i) 有効的な特異点回避
(ii) 特徴的な重力波スペクトルの生成
が同一の機構から生じ得ることを示す概念実証である。
完全な動的理論への埋め込みおよび定量的予言の精密化は今後の重要課題であるが、本研究は双対性が宇宙初期に与える観測的帰結を理解するための出発点を提供する。

結論
本研究は、新しい宇宙論モデルを構築するものではなく、
「DFT 双対性が宇宙初期に作用した場合の観測的帰結を、最小有効モデルで概念的に探索する」
という目的を達成した。

特異点回避と blue-tilted 重力波が同一機構から生じるという構造は、
DFT 的宇宙論の観測的側面を理解する上で重要な示唆を与える。

【Appendix A: 共形時間変換と a_eff'' / a_eff の近似導出】
本付録では、本文 4.2 節で用いた近似式
  a_eff'' / a_eff ≈ mu^2 ( 1 - 2 sech^2( mu eta ) )
がどのように得られるかを、最小限の仮定のもとで説明する。
ここでの目的は厳密導出ではなく、Pöschl–Teller 型ポテンシャルが自然に現れる理由を明確化することである。

A.1 共形時間の定義
共形時間 eta は
  d eta = dt / a_eff(t)
で定義される。
バウンス近傍では a_eff(t) が最小値を持つため、eta = 0 をバウンス中心に取る。

A.2 a_eff(t) の局所展開
本文で定義した有効スケール因子は
  a_eff^2(t) = a^2(t) cos^2 theta(t) + a^{-2}(t) sin^2 theta(t)
である。

キンク解
  theta(t) = 2 arctan( exp( mu t ) )
を代入すると、バウンス近傍 t ≈ 0 では
  cos theta ≈ 0
  sin theta ≈ 1
となるため、a_eff(t) は局所的に
  a_eff(t) ≈ a_0 [ 1 + ε tanh^2( mu t ) ]^{1/2}
の形で近似できる。
ここで ε は O(1) の定数であり、「a_eff(t) が滑らかな最小値を持つ」という性質のみを要求する。

A.3 共形時間への変換
d eta = dt / a_eff(t) を積分すると、バウンス近傍では
  eta ≈ t / a_0
が成り立つ（高次補正は t^3 以上）。
したがって、局所的には
  t ≈ a_0 eta
と置き換えてよい。
A.4 a_eff(eta) の形
t = a_0 eta を代入すると
  a_eff(eta) ≈ a_0 [ 1 + ε tanh^2( mu a_0 eta ) ]^{1/2}
となる。
ここで、バウンスの時間幅を
  mu_eff = mu a_0
と再定義すると、
  a_eff(eta) ≈ a_0 [ 1 + ε tanh^2( mu_eff eta ) ]^{1/2}
と書ける。

A.5 二階微分の計算と ε 依存性
a_eff(eta) を二階微分すると、主要項は
  a_eff'' / a_eff ≈ mu_eff^2 ( 1 - C(ε) sech^2( mu_eff eta ) )
という形になる。
ここで C(ε) は ε に依存する O(1) の正の定数であり、一般には
  C(ε) ≠ 2
である。
重要なのは：

• tanh(mu_eff eta) の二階微分が sech^2 を生むこと

• a_eff(eta) が最小値を持つため、sech^2 の係数が負になること

• 係数の数値は ε に依存するが、「定数 −C(ε) × sech^2」という形は普遍的であること

である。

本文では、Pöschl–Teller 型ポテンシャルの標準形との比較を容易にするため、
  a_eff'' / a_eff ≈ mu^2 ( 1 - 2 sech^2( mu eta ) )
という代表的な形を採用している。
これは C(ε) = 2 に対応する特定の ε の選択、あるいは rescaling による有効的なパラメータ固定とみなすことができる。

A.6 結論：Pöschl–Teller 型が自然に現れる理由
以上の近似により、
  a_eff'' / a_eff ≈ mu^2 ( 1 - C sech^2( mu eta ) )
という Pöschl–Teller 型ポテンシャルが自然に現れる。
C は O(1) の正の定数であり、その具体値はモデルの詳細に依存するが、

• バウンス幾何から sech^2 形状が現れること

• 井戸型ポテンシャルとして振る舞うこと

という定性的特徴は C の値に依存しない。

本文で用いた C = 2 は、標準的な Pöschl–Teller 形との比較を簡潔にするための代表値であり、
本研究の主張（blue tilt・ピーク・カットオフという三段構造）は定性的な三段構造には影響しないが、スペクトルの詳細な形状はCの値に依存する。

【Appendix B: 作用の変分と摩擦項無視の条件】
本付録では、本文 2.4 節で用いた作用
S = ∫ dt d^3x a_eff^3(t) [ - 3 H_eff(t)^2 / (8 pi G)
+ (1/2) (theta_dot)^2
- V(theta) ]

に対して theta(t) の変分を行い、
摩擦項および重力項の寄与が mu >> H の極限で無視できることを明示的に示す。
本研究はトイモデルであり、DFT の完全な変分構造を再現することを目的としない。
ここでは、最小限の整合性を保つための導出を行う。

B.1 作用の構造
ラグランジアンは
L = a_eff^3 [ (1/2) theta_dot^2 - V(theta) ]
- a_eff^3 * 3 H_eff^2 / (8 pi G)

である。
ここで
・a_eff = a_eff(theta)
・H_eff = (1/a_eff) d(a_eff)/dt
であり、theta に非線形に依存する。

B.2 変分の一般形
オイラー＝ラグランジュ方程式は
d/dt ( ∂L/∂theta_dot ) - ∂L/∂theta = 0
まず
∂L/∂theta_dot = a_eff^3 theta_dot
したがって
d/dt ( a_eff^3 theta_dot )
= a_eff^3 theta_ddot + 3 a_eff^2 a_eff_dot theta_dot
ここで第 2 項が 摩擦項 に相当する。

B.3 ∂L/∂theta の構造
∂L/∂theta は以下の 3 つに分解される：

1.ポテンシャル項

• a_eff^3 V'(theta)
  2.a_eff(theta) の依存性からの項
  (∂ a_eff^3 / ∂theta) * [ (1/2) theta_dot^2 - V(theta) ]
  3.重力項からの寄与

• ∂/∂theta [ a_eff^3 * 3 H_eff^2 / (8 pi G) ]
  この 3 が今回の重要点である。

B.4 摩擦項の大きさ
摩擦項は
3 (a_eff_dot / a_eff) theta_dot
である。

バウンス近傍では
theta_dot ~ mu / cosh(mu t)

一方、背景膨張は
a_dot / a = H
であるため、
a_eff_dot / a_eff = O(H)
したがって摩擦項は
摩擦項 ~ H * mu

B.5 重力項の寄与の評価
重力項は
L_grav = - a_eff^3 * 3 H_eff^2 / (8 pi G)
である。

ここで
H_eff = (1/a_eff) d(a_eff)/dt
であり、a_eff の theta 依存性から
H_eff = H + O(theta_dot)
が成り立つ。

バウンス中心では
theta_dot ~ mu
であるため、重力項の変分は
∂L_grav/∂theta ~ a_eff^3 * H_eff * ∂H_eff/∂theta

ここで
∂H_eff/∂theta ~ O(H)
が成り立つ（theta_dot の寄与は偶関数性により抑制される）。

したがって重力項の寄与は
重力項の寄与 ~ H * mu
となる。

B.6 主導項 theta_ddot の大きさ
キンク解の時間スケールは
theta_dot ~ mu
theta_ddot ~ mu^2
である。

したがって主導項は
主導項 ~ mu^2

B.7 摩擦項・重力項の抑制
以上より、
摩擦項 ~ H mu
重力項 ~ H mu
主導項 ~ mu^2

したがって
(摩擦項) / (主導項) ~ H / mu << 1
(重力項) / (主導項) ~ H / mu << 1
となる。

B.8 結論：mu >> H の極限で運動方程式は簡約される
以上の評価により、mu >> H の極限では
・摩擦項は無視できる
・重力項の変分も無視できる
・主導項は theta_ddot と V'(theta) のみ
となる。

したがって運動方程式は
theta_ddot + V'(theta) ≈ 0
に簡約される。

この方程式は
theta(t) = 2 arctan( exp( mu t ) )
を厳密解として持つ。

【Appendix C: Pöschl–Teller ポテンシャルの標準形と本研究での近似】
本付録では、本文 4.2 節で用いた近似式

a_eff'' / a_eff ≈ mu^2 ( 1 - C sech^2( mu eta ) )

が、Pöschl–Teller ポテンシャルの標準形とどのように対応しているかを整理する。
本研究は厳密散乱係数の導出を目的とせず、散乱構造の定性的特徴を抽出するための近似を採用している。

C.1 Pöschl–Teller ポテンシャルの標準形
Pöschl–Teller ポテンシャルの標準形は
　　U_PT(x) = - λ (λ + 1) sech^2( α x )
である。
ここで

• λ : ポテンシャルの深さを決めるパラメータ

• α : ポテンシャルの幅を決めるパラメータ

である。

散乱方程式は
  d^2 u / dx^2 + [ k^2 + λ (λ + 1) sech^2( α x ) ] u = 0
となる。
その反射・透過係数は特殊関数を用いて表されるが、本研究ではその厳密形を用いない。

C.2 本研究で現れる形との比較
本研究で得られた有効ポテンシャルは
  U_eff(eta) = a_eff'' / a_eff ≈ mu^2 ( 1 - C sech^2( mu eta ) )
である。

これを標準形と比較すると、
  U_eff(eta) = mu^2 - C mu^2 sech^2( mu eta )
となり、
  α = mu
  λ (λ + 1) ∝ C
という対応が得られる。
C は Appendix A で述べたように O(1) の正の定数であり、λ も O(1) の値を取る。

C.3 散乱係数の性質（定性的議論）
Pöschl–Teller ポテンシャルの散乱係数は、一般に

• 高 k で指数的減衰を示す

• ある k スケール付近で反射が最大になる

• 低 k では反射が抑制される

といった性質を持つことが知られている。

本研究では、このうち

• 高 k で |beta_k|^2 ∝ exp( - const * k / mu ) となること

• 中間 k で励起が最大になること

• 低 k で抑制されること

という 定性的特徴のみ を用いる。

C.4 本研究での近似形との関係
本文で採用した近似形
  |beta_k|^2 ≈ 1 / ( exp( 2 pi k / mu ) - 1 )
は、Pöschl–Teller ポテンシャルの厳密解ではなく、

• 高 k での指数的減衰

• 中間 k での励起ピーク

• 低 k の抑制

という散乱構造の定性的特徴を、熱的分布に類似した形でパラメトライズしたものである。

この近似は、

• 有限時間幅バウンスの散乱が熱的類似励起を生むこと

• Pöschl–Teller 型ポテンシャルが指数的減衰を伴う反射係数を持つこと

と整合的である。

C.5 結論：標準形との整合性

以上より、本研究で用いた近似式は

• Pöschl–Teller ポテンシャルの標準形と構造的に整合しており

• 高 k の指数的減衰と中間 k のピークという特徴を正しく捉え

• 散乱構造の定性的理解を与える

ことが確認できる。

本研究は厳密散乱係数の導出を目的としないため、
近似形の採用は本モデルの位置づけ（conceptual exploration）と整合的であり、Appendix C はその定性的妥当性を補足する役割を果たす。

Appendix D

#λ依存込みの最終スペクトル図コード
import numpy as np
import matplotlib.pyplot as plt

# =========================
# パラメータ
# =========================
mu_eff = 1.12  # ← フィット結果を入れる

# λ の代表値（理論の予言として提示）
lambda_values = [0.5, 1.0, 1.5]

# =========================
# kレンジ（対数スケール）
# =========================
k = np.logspace(-2, 2, 500) * mu_eff

# =========================
# β_k（厳密解）
# =========================
def beta_exact_sq(k, lam):
    numerator = np.sin(np.pi * lam)**2
    denominator = np.sinh(np.pi * k / mu_eff)**2 + numerator
    return numerator / denominator

# =========================
# プロット
# =========================
plt.figure(figsize=(7,5))

for lam in lambda_values:
    P = k**3 * beta_exact_sq(k, lam)
    P /= np.max(P)  # 正規化（形状比較）
    plt.loglog(k / mu_eff, P, linewidth=2,
               label=rf"$\lambda = {lam}$")

# =========================
# 軸・装飾
# =========================
plt.xlabel(r"$k / \mu_{\rm eff}$", fontsize=12)
plt.ylabel(r"$P_t(k)$ (normalized)", fontsize=12)

plt.legend()
plt.grid(True, which="both", ls="--", alpha=0.5)

plt.tight_layout()

# 保存（論文用）
plt.savefig("gw_spectrum_lambda_dependence.png", dpi=300)

plt.show()

Appendix E: 観測されるピーク周波数から NEC 違反スケールを逆推定する

本付録では、本研究で得られたテンソルパワースペクトルのピーク位置と、バウンス幾何における NEC（Null Energy Condition）違反の時間スケールとの対応関係を整理する。特に、観測される重力波スペクトルのピーク周波数から、バウンスの特徴的スケールを逆推定する方法を示す。

E.1 バウンススケールとピーク波数の対応

本研究では、バウンス近傍の有効ポテンシャル
  U(eta) = a_eff''(eta) / a_eff(eta)
が Poschl–Teller 型関数
  U(eta) ≃ mu_eff^2 * ( 1 - C * sech^2(mu_eff * eta) )
で良く近似されることを示した。
このとき、テンソルモードの散乱係数から得られるパワースペクトルは
   P_t(k) ∝ k^3 * |beta_k|^2
を持ち、そのピーク位置は
  k_peak ≃ mu_eff
で与えられる。
したがって、スペクトルのピーク波数はバウンスの特徴的スケール mu_eff を直接反映する。

E.2 観測周波数への変換

観測される物理周波数 f_peak は、生成時のスケール因子 a_* と現在のスケール因子 a_0 を用いて
  f_peak = (k_peak / 2π) * (a_* / a_0)
で与えられる。これを反転すると
  mu_eff ≃ k_peak ≃ 2π * f_peak * (a_0 / a_*)
となる。
ここで a_0 / a_* は宇宙の熱史（再加熱温度、放射優勢期、物質優勢期など）に依存するため、観測から mu_eff を決定するには熱史の仮定が必要となる。

E.3 NEC 違反の時間幅の推定

シャープバウンス極限 mu >> H では、共形時間 eta と物理時間 t はバウンス近傍でほぼ一致し、
  a_eff'' / a_eff ≃ mu_eff^2 * ( 1 - 2 * sech^2(mu_eff * t) )
が成立する。このとき NEC 違反は
  dot(H_eff)(t) > 0
となる時間領域で生じ、その典型的な幅は
  Delta t_NEC ≃ alpha / mu_eff
で与えられる。
ここで alpha は O(1) の数値係数であり、NEC 違反領域の定義に依存する。
観測周波数 f_peak を用いると
  Delta t_NEC ≃ alpha / ( 2π * f_peak ) * (a_* / a_0)
となる。
すなわち、高周波ピークが観測されるほど NEC 違反は短時間であり、低周波ピークほど NEC 違反は長時間となる。

E.4 NEC 違反の強さの推定（GR 的有効記述）

GR 的な有効エネルギー条件を仮定すると
  rho_eff + P_eff = -(1 / 4πG) * dot(H_eff)
であり、バウンス近傍では
  dot(H_eff)_max ≃ beta * mu_eff^2
が成立する（beta は O(1) の係数）。
したがって NEC 違反の強さは
  rho_eff + P_eff ≃ -(beta / 4πG) * mu_eff^2
で与えられる。
観測周波数を用いると
  rho_eff + P_eff ≃ -(beta / 4πG) * ( 2π * f_peak * (a_0 / a_*) )^2
となり、ピーク周波数が高いほど NEC 違反の強さも大きくなる。

E.5 まとめ

本付録で示したように、本研究のトイモデルでは
(i) k_peak ≃ mu_eff
(ii) Delta t_NEC ≃ O(1) / mu_eff
(iii) NEC 違反の強さ ∝ mu_eff^2
という対応が成り立つ。
したがって、観測されるピーク周波数 f_peak が分かれば、宇宙の熱史を仮定することで
・バウンスの時間スケール（NEC 違反の幅）
・NEC 違反の強さ（rho + P の負の大きさ）
を逆推定することが可能となる。
これは、DFT に着想を得たバウンス宇宙論が、将来の重力波観測（DECIGO, BBO, LISA など）と直接結びつく可能性を示すものである。