DSWIM:Efficient and Stable Deterministic Computation of Warm Inflation Perturbations

技術要約：DSWIM – ウォーム・インフレーション摂動の効率的かつ安定的な決定論的計算

問題提起
ウォーム・インフレーション（WI）の摂動は、量子ゆらぎと熱ゆらぎの両方によって駆動されるため、結合された確率微分方程式を解く必要がある。公開されているコードであるSWIMは、堅牢な確率論的アプローチを実装しているが、統計的な収束を得るために多くの実現例に対してアンサンブル平均を行う必要があり、パラメータ推論や尤度評価において計算コストが高くなる。

フォッカー・プランク法に基づく決定論的な代替手法は、摂動変数の相関行列を常微分方程式（ODE）を通じて進化させることで、計算コストの大幅な削減を実現できる。しかし、著者らは、直接的な決定論的進化が、多くのWIモデルにおいて深刻な数値的不安定性（ill-conditioning）に陥ることを指摘している。この不安定性は、摂動変数が極めて異なるスケールで進化するために発生し、相関行列の要素が多くの桁数（テストされたモデルでは13から90桁）にわたって広がることに起因する。この階層構造は、浮動小数点演算の丸め誤差、補正係数 $G(Q)$ の人工的な歪み、そして数値進化の破綻（特に強い散逸領域において）を引き起こす。さらに、従来の決定論的実装（例：WI2easy）は、相関のある熱雑音の不適切な取り扱いにより、中程度の散逸領域において確率論的結果との系統的な不一致を示していた。

手法
これらの課題に対処するため、著者らはSWIMフレームワーク内の新しいモジュールとして、決定論的進化方程式にスケーリング行列フレームワークを実装したDSWIM（Deterministic SWIM）を導入する。

1. スケーリング変換： 著者らは、スケーリングされた状態ベクトル $\tilde{\Phi} = S\Phi$ を定義する。ここで、$S$ はハッブルパラメータ $H$ の累乗からなる対角スケーリング行列である。このスケーリングは、インフレーション中の摂動変数の有効なスケーリング挙動に基づき、物理的に動機付けられている：

   • メトリック摂動 $\psi$（無次元）はスケーリングされない。
   • インフラトンゆらぎ $\delta\phi$ およびその導関数は $H^{-1}$ でスケーリングされる。
   • 放射エネルギー密度摂動 $\delta\rho_r$ は $H^{-2}$ でスケーリングされる。
   • 放射運動量摂動 $\delta q_r$ は、摂動方程式における出現に基づき $H^{-1}$ でスケーリングされる。
     この変換は変数を正規化し、進化する相関行列の要素が数値的に比較可能な状態に保たれるようにすることで、ダイナミックレンジを維持し、精度の喪失を防ぐ。

2. 物理的観測量の保存： 著者らは、このスケーリング変換が共動曲率摂動 $R$、ひいては原始スカラーパワースペクトル $P_R$ を正確に不変に保つことを証明している。変換されたドリフト行列 $\tilde{A}$ および拡散行列 $\tilde{D}$ は、システムの物理的一貫性を維持するように導出されている。

3. ノイズ相関の修正： 決定論的形式においては、拡散行列 $D$ における非対角成分の相関を明示的に組み込むよう更新されている。著者らは、インフラトンと放射の熱雑音項は、エネルギー・運動量保存則により独立ではなく、互いに相関していることを示している。これらの非ゼロの非対角項を組み込むことで、中程度の散逸領域で以前観察されていた確率論的コードと決定論的コードの間の系統的な不一致が解消される。

主な結果

• 数値的安定性： スケーリングされた決定論的フレームワークは、スケーリングされていない実装で見られる数値的アーティファクトを効果的に抑制する。ランナウェイ指数関数ポテンシャルモデルにおいて、このスケーリング変換は共分散行列の階層を約40〜50桁減少させ、決定論的ソルバーの安定進化範囲を、スケーリングなしの手法と比較して散逸比 $Q$ に関してほぼ1桁拡大させた。
• 確率論的手法との一致： DSWIMの結果は、弱い、中間、および強い散逸領域において、確率論的なSWIM実装と極めて良好に一致しており、スケーリングされた決定論的アプローチの物理的な正確性を検証している。
• 計算効率： 決定論的アプローチは大幅な高速化をもたらす。代表的なモデルにおいて、DSWIMは補正係数 $G(Q)$ の計算時間を、確率論的なSWIMと比較して約2桁短縮し、既存の決定論的コードであるWI2easyを大幅に上回る性能を示した。例えば、4次ポテンシャルモデルの実行時間は、約116秒（確率論的）から約1.56秒（決定論的）へと減少した。
• 不一致の解消： 拡散行列における相関のある熱雑音構造を正しく実装することにより、DSWIMは中程度の散逸領域（$0.01 \lesssim Q \lesssim 10$）における確率論的結果と決定論的結果の間の系統的な偏差を排除した。

意義と主張
本論文は、DSWIMを、ウォーム・インフレーションのスカラ摂動を計算するための、数値的に堅牢かつ計算効率の高いフレームワークとして確立するものである。著者らは、DSWIMが、Cobaya尤度解析フレームワークと直接インターフェースを持つ、数値的に安定化された決定論的形式を提供する、現在唯一の公開されているWI摂動ソルバーであると主張している。この能力により、確率論的な平均化のコストや、スケーリングされていない決定論的手法の数値的不安定性のために、以前は計算的に困難であったウォーム・インフレーション・モデルの効率的なパラメータ推論研究が可能となる。本研究は新しい物理モデルを提案するものではなく、既存のウォーム・インフレーション・シナリオのパラメータ空間を正確かつ効率的に探索するために必要な計算ツールを提供するものである。