Lagrangian Bias as a Gaussian Random Field

この論文「Lagrangian Bias as a Gaussian Random Field（ガウス確率場としてのラグランジュ的バイアス）」は、宇宙論におけるハローバイアス（銀河団や暗黒物質ハローの物質分布に対する偏り）の理解に対して、従来の摂動論的な枠組みとは根本的に異なる新しい視点を提供するものです。

以下に、論文の技術的な要約を問題提起、手法、主要な貢献、結果、そして意義の観点から詳述します。

1. 問題提起 (Problem)

従来の宇宙論における「トレーサーバイアス（tracer bias）」の枠組みでは、バイアスは観測可能なトレーサー（ハローや銀河）の過密度を、密度場の摂動展開の係数として定義するアプローチが主流でした（ピーク・バッキングスプリット、Separate Universe 法、有効場理論など）。

• 従来の見方: バイアスは、ハローが形成される「崩壊過程（collapse process）」によって生成される、トレーサーに固有の性質である。
• 課題: このアプローチでは、バイアスがなぜ特定の値をとるのか、また固定された質量を持つハロー群の中でなぜバイアスの分布（アセンブリーバイアス）が存在するのかを、第一原理から完全に説明することが困難でした。

2. 手法と理論的枠組み (Methodology)

著者は、初期密度場が**ガウス確率場（Gaussian Random Field; GRF）**であるという性質に立ち返り、バイアスを「係数」ではなく「場（field）」として再定義しました。

• ラグランジュ的バイアス場の導出:

  • ガウス場における任意の点 $q$ における過密度 $\delta(q)$ を条件としたとき、その点に存在するすべての独立したフーリエモード（波数 $k$ のシェル）の平均値が、そのモードの分散 $\sigma_k^2$ に対して一定の比率でシフトすることを示しました。
  • 具体的には、条件付き平均 $\langle \delta_k(q) | \delta(q) \rangle$ は $\sigma_k^2 \frac{\delta(q)}{\sigma^2}$ となり、これを分散で規格化した「バイアス」$b(q) = \frac{\delta(q)}{\sigma^2}$ は波数 $k$ に依存しない（スケール非依存性）ことが導かれます。
  • この結果、ラグランジュ空間のすべての点に定義された、スケール非依存のバイアス場 $b(q)$ が存在することが証明されました。これは初期条件にエンコードされた統計的構造そのものです。

• ハローバイアスの幾何学的解釈:

  • ハローは、この事前に存在するバイアス場 $b(q)$ 内から、崩壊基準（collapse criterion）によって「パッチ（領域）」が幾何学的に選択された結果であると解釈されます。
  • ハローのバイアスは、そのハローを構成するラグランジュ的パッチ内の全粒子のバイアス値の平均値 $\langle b \rangle$ として定義されます。

• 二次バイアス場と「アクティブ」な演算子:

  • ラグランジュ密度場 $\delta$ に対する線形演算子（例：ラプラシアン $\nabla^2 \delta$ や潮汐テンソル）を適用すると、新たなガウス場（二次バイアス場）が生成されます。
  • 崩壊基準が $\delta$ だけでなく、これらの二次場（例：$\nabla^2 \delta$）とも相関している場合、それらの場は「アクティブ（active）」となり、ハローの集積特性に独立した情報を加えます。

3. 主要な貢献と結果 (Key Contributions & Results)

A. バイアス場のスケール非依存性と UV 安全性

• 数値シミュレーション（Quijote シミュレーションの初期条件）を用いた検証により、任意の点 $q$ において、異なる波数 $k$ のシェルから推定されるバイアス値が、理論予測 $\delta(q)/\sigma^2$ に一致し、$k$ に依存しないことを実証しました。
• これにより、バイアスは高波数（UV）領域の寄与によって不安定になるのではなく、ガウス構造によって厳密に定義され、UV に対して安全（UV safe）であることが示されました。

B. ハローバイアスの予測精度

• 質量 $M$ のハローのバイアス $b(M)$ は、そのハローのラグランジュ的パッチに含まれる粒子の平均バイアス値 $\langle b_{\text{halo}} \rangle$ を線形成長因子 $D(a)$ でスケーリングし、$+1$ を加えることで（$b_E = \langle b_{\text{halo}} \rangle/D + 1$）、パラメータなしで正確に予測できることを示しました。
• 図 2 に示されるように、シミュレーションで測定された $z=0$ でのハローバイアスと、この理論予測は完全に一致します。

C. アセンブリーバイアスと二次バイアス場の説明

• アセンブリーバイアスの起源: 固定された質量を持つハロー群において、バイアス値 $b_\delta$ がガウス分布に従うこと、およびその分布の幅が生じる理由は、崩壊基準が $\delta$ だけでなく $\nabla^2 \delta$（曲率）などの二次場とも相関しているためであると説明しました。
• アクティブな場の役割: 崩壊基準が $\nabla^2 \delta$ に依存する場合、ハロー群における $\delta$ と $\nabla^2 \delta$ の共分散構造が背景場とは異なり（$\langle \delta_s \nabla^2 \delta_s \rangle_{\text{halo}} \neq \langle \delta_s \nabla^2 \delta_s \rangle_V$）、これが「アクティブ」なバイアス場を生み出します。
• 実証: 図 3 と図 4 に示すように、$\delta_s$ を固定して $\nabla^2 \delta_s$ でハローを分割すると、大規模構造におけるバイアスに明確な差（バイアス・スプリット）が生じることが確認されました。これは、$\nabla^2 \delta$ が崩壊プロセスに能動的に関与していることを直接証明するものです。

D. 統一的理解

• この枠組みは、ピーク・バッキングスプリット、有効場理論（EFT）のバイアス展開、アセンブリーバイアス、および非ガウス性（local non-Gaussianity）によるスケール依存バイアスなど、一見異なる現象がすべて「ガウス演算子の共分散構造の修正」という単一のメカニズムから導かれることを示しました。

4. 意義 (Significance)

1. パラメータフリーの第一原理モデル:
   従来の摂動論的アプローチでは、バイアス係数をシミュレーションからフィッティングする必要がありましたが、この枠組みでは初期条件（ラグランジュ場）から直接、パラメータなしでハローの集積特性を予測する「ab initio モデル」への道筋を示しました。

2. 物理的直観の革新:
   バイアスを「崩壊によって生成される係数」から、「初期条件に内在する場から幾何学的に選択される量」へと概念を転換させました。これにより、ハローのクラスタリングが、重力による進化以前に決定されている側面を持つことが明確になりました。

3. 将来の応用:

   • このバイアス場はラグランジュ的変数であるため、非線形進化においても物質流に沿って輸送（advection）され、保存されます。これにより、非線形スケールにおける情報の混入や消失を新しい言語で記述する可能性があります。
   • 中性子星質量（massive neutrinos）によるスケール依存バイアスや、他の宇宙論的パラメータの影響を、このバイアス場フレームワーク内で統一的に扱うことが可能になります。
   • 最終的には、崩壊モデルとこのバイアス場を組み合わせた、完全なハロークラスタリングの予測モデルの構築が期待されます。

結論として、この論文はハローバイアスの本質を「統計的な場」の選択問題として再定義し、観測される複雑なバイアス現象（質量依存性、アセンブリーバイアス、スケール依存性）を、初期ガウス場の統計的構造と幾何学的選択の観点から統一的に説明する強力な理論的基盤を提供しています。