Benchmarking wall velocities in cosmological phase transitions: Fluid Ansatz and WallGo

技術要約：宇宙論的相転移における壁速度のベンチマーク解析

問題提起
宇宙論的な一次相転移における膨張する真空バブルの終端速度（$v_w$）を正確に決定することは、重力波（GW）やバリオン非対称性などの観測的シグネチャーを予測する上で極めて重要である。この速度は、真空による駆動圧と、プラズマによる流体力学的摩擦のバランスによって決定される。この摩擦を計算するには、バブル壁近傍における非平衡粒子分布関数（$\delta f$）に関する半古典的なボルツマン方程式を解く必要がある。

しかし、完全な積分微分形式のボルツマン方程式を解くことは、モデル構築者にとって計算コストが非常に高い。そのため、コミュニティでは、この問題を扱いやすい方程式系へと簡略化する近似手法に依存している。現在、主に2つのアプローチが存在する：

1. 流体仮説（Fluid Ansatz）： 非平衡分布が平衡分布（ボース＝アインシュタイン／フェルミ＝ディラック分布）の関数形式を保持していると仮定し、変動するパラメータ（化学ポテンシャル、温度、速度）を用いる手法。これにより、ゆらぎをマクロな流体力学的特性へと写像する。
2. スペクトル法（WallGo）： 非平衡分布関数を直交チェビシェフ多項式の基底で展開する手法であり、公開コードであるWallGoに実装されている。

これらの手法が広く用いられているにもかかわらず、両者の整合性を検証したり、乖離が生じる領域を特定したりするための系統的な比較はこれまで行われてこなかった。

手法
著者らは、以下の2つの理論的枠組みを用いて、流体仮説とWallGoスペクトル法の厳密なベンチマークを実施した：

• SMEFTモデル： 標準模型有効場理論（SMEFT）に次元6の演算子（$\phi^6$）を拡張したモデル。これにより、物理的なヒッグス質量（$m_h = 125$ GeV）を維持しつつ、強い一次相転移を可能にしている。
• ライト・ヒッグス・モデル： ヒッグス質量を人工的に低減させた（$m_h \lesssim 80$ GeV）標準模型に近いポテンシャルを用い、より弱い一次相転移を生成するモデル。

本研究は以下の3段階で進行する：

1. 線形比較： 両手法を、トップクォークの対消滅（$t\bar{t} \to gg$）のみを考慮する場合と、その後に散乱過程（$tq \to tq$, $tg \to tg$）を含める場合の両方で適用し、$v_w$を計算する。流体仮説は、運動量展開の次数を3次（10個のゆらぎ）とするグリーン法を実装したカスタムPythonコードを用いて解かれ、一方WallGoは$N=11$および$M_{grid}=30$のチェビシェフ格子を使用する。
2. 摂動論的解析： 無次元のゆらぎの総和として定義される摂動パラメータ$R$を導入する。これにより、流体仮説に内在する線形近似（$\delta f \ll 1$）の破綻を監視する。
3. 非線形拡張： 強固な相転移における不一致の要因を調査するため、分布関数のテイラー展開において2次の項（$O(\delta^2)$）を含むよう流体仮説の形式を拡張する。これには、二次的なゆらぎが壁速度に与えるバックリアクション（背後からの影響）を定量化するために、非線形境界値問題（BVP）を解くプロセスが含まれる。

主な結果

• 弱い領域での一致： 相転移の強さパラメータ $\alpha \lesssim 0.01$ の場合、両手法は実質的に同一の終端速度を与える。この一致は、トップクォークの対消滅チャネルのみを考慮した場合に特に強固である。
• より強い領域での乖離： 相転移の強さが増大するにつれ（$\alpha \gtrsim 0.1$）、不一致が生じる。散乱過程を含めると、流体仮説とWallGoの結果は$O(40\% - 60\%)$の乖離を示す。
• 線形化の破綻： この乖離は、摂動パラメータ$R$が1に近づく（$R \gtrsim 0.1$）ことと相関しており、これは流体仮説における分布関数の線形近似が破綻していることを示している。同時に、壁の厚さパラメータ $L_w T$ も1に近づいており、ボルツマン方程式の基礎となるWKB近似自体の妥当性に疑問を投げかけている。
• 非線形項の影響： 流体仮説に$O(\delta^2)$項を含めると、強固な相転移における予測終端速度が大きく変化し、一般に$v_w$が減少（摩擦が増大）する。極めて重要な点として、二次的な寄与が圧力全体に与える直接的な影響は定量的に小さい（全圧力の約4%）ものの、一次のゆらぎに対する非線形なバックリアクションが、マクロな壁の動力学に大きなシフトを引き起こすことが判明した。このシフトは、結果としてWallGoとの一致を改善するのではなく、むしろ悪化させる。
• 収束性： 本研究は、WallGoが十分な格子解像度とともに収束すること、また流体仮説が弱い領域において運動量展開の次数 $\kappa=3$（10パラメータ）で収束することを確認している。

意義と主張
本論文は、流体仮説とWallGoはともに弱から中程度の相転移においては計算効率が高く信頼できるが、線形化された流体仮説は強い相転移領域（$\alpha \gtrsim 0.01$）において根本的に破綻するということを主張している。

著者らは、この不一致は単なる数値的なアーティファクトではなく、マクロなゆらぎが大きくなった際に線形展開が機能しなくなることに起因すると強調している。さらに、これらの不一致が最も顕著になる強い相転移領域において、壁の厚さが縮まり（$L_w T \lesssim 3$）、ボルツマン方程式の基礎である半古典的なWKB近似がもはや有効ではない限界に近づいているという、より深い理論的限界についても指摘している。

したがって、本論文は、強固な相転移（まさに将来の重力波観測所がターゲットとする領域）におけるバブル壁速度の高精度な予測を実現するためには、現在の半古典的なアプローチを洗練させるだけでは不十分であると結論付けている。これは、現在の手法（流体仮説およびWallGoの両方）が、最も物理的に興味深いシナリオにおいて、その有効領域の外側で動作している可能性があることを示唆しており、WKB近似を超えた取り扱いが必要であることを示している。