Control variates from Eulerian and Lagrangian perturbation theory: Application to the bispectrum

技術的サマリー：オイラー型およびラグランジュ型摂動論に基づく制御変量：ビスペクトルへの応用

問題定義
宇宙論的 N 体シミュレーションは、特にビスペクトルのような高次統計量に対して、宇宙のウェブの非ガウス構造をモデル化する上で不可欠である。しかし、これらのシミュレーションは初期条件のランダム性に起因する「宇宙論的分散（cosmic variance）」の影響を受ける。第 IV 段階の銀河サーベイに必要な高精度を達成するためには、通常、多数のシミュレーション実現の平均化が必要となるが、これは計算コスト的に実行不可能である。「ペアリングと固定（pairing and fixing）」のような分散低減技術は存在するが、パワースペクトルに比べビスペクトルに対しては効果が低い。制御変量法は、安価で相関のあるシミュレーション（代理モデル）を利用して観測量の標本分散を低減するという、有望な代替手段を提供する。重要な課題は、広範なスケールにわたって完全な N 体ビスペクトルと高い相関を持ち、かつ解析的に扱いやすい（高価なアンサンブルなしで平均を決定できる）代理モデルを特定することである。

手法
著者らは、$z=0.5$ における物質ビスペクトルを事例研究として、オイラー型摂動論（EPT）およびラグランジュ型摂動論（LPT）から導出された摂動的制御変量を調査した。本研究では、基準となる真実（ground truth）として Quijote N 体シミュレーション・スイート（$N=1000$ の実現）を利用した。

1. 玩具モデル解析：著者らはまず、「真の」場がゼリコビッチ近似（ZA）によって生成される簡略化されたモデルを解析した。彼らは、平均がゼロである「線形（ガウス）」制御変量が、非線形場に対して $\exp[-\frac{1}{2}\Sigma^2(k_1^2 + k_2^2 + k_3^2)]$ に比例する指数関数的な非相関を示すことを実証した。ここで $\Sigma$ は変位分散である。
2. オイラー型制御変量：著者らは密度の 5 次まで（$O(\delta^5)$）の EPT 場を用いて制御変量を構築した。彼らは、カウントパラメータ $\epsilon$ を用いて、完全な格子実現から樹形レベルの寄与（$O(\epsilon)$）と平均ゼロの高次ダイアグラム（$O(\epsilon^2)$）を分離した。解析的予測と格子測定を一致させる課題に対処するため、彼らはビン化効果を補正する「有効三角形（effective triangle）」および離散平均化スキームを開発した。
3. ラグランジュ型制御変量：著者らは ZA および第 2 次ラグランジュ型摂動論（2LPT）のビスペクトルを評価した。彼らは、2LPT が正しい樹形レベル構造と大規模なバルクフローを捉える一方で、その完全なビスペクトルは解析的に計算できず、経験的な平均推定を必要とすることを示した。一方、ZA ビスペクトルは解析的に計算可能であるが、EPT と比較して樹形レベル構造が誤っており、低-$k$ 領域で相関が低いという結果となった。
4. シフト制御変量：両アプローチの利点を組み合わせるため、著者らは「シフト制御変量」を導入した。これは「シフト演算子」基底を用いて構築され、場は線形（ゼリコビッチ）変位のみが指数化されたラグランジュ型バイアストレーサとして定義される。この構成により、場は正しい樹形レベルビスペクトル（EPT と一致）を保持しつつ、線形変位を再総和（resum）し、LPT の相関特性を維持することが保証される。
5. 解析的計算：ZA ビスペクトルについて、著者らは樹形レベル計算から導出された重みと、広幅ビンに対するモンテカルロサンプリングを用いて、ビン平均された解析的予測を計算する手法を開発し、格子測定との間で 1% 未満の一致を達成した。

主要な貢献と結果

• 非相関の限界：本研究は、EPT ベースの制御変量が変位分散より小さいスケール（$k \gtrsim \Sigma^{-1}$）において指数関数的な非相関に陥り、その有用性が大規模スケールに限定されることを確認した。
• ZA と 2LPT の比較：ゼリコビッチ近似は解析的に扱いやすいが、その樹形レベルビスペクトルが EPT の予測と異なるため、パワースペクトルに比べビスペクトルに対する制御変量として最適ではない。2LPT はより良い相関を提供するが、平均の解析的扱いやすさがない。
• シフト制御変量の性能：提案された「シフト制御変量」は、他のすべての摂動手法を上回る性能を示した。これは、EPT の正しい樹形レベル構造と、LPT からの線形変位の再総和を組み合わせる。
  • 分散低減：シフト制御変量を単一の $V = 1 (\text{Gpc}/h)^3$ N 体シミュレーションに適用することで、最も低-$k$ の三角形に対する物質ビスペクトルの分散を、$10^4$ 回のシミュレーションの平均化に相当する因子で低減した。
  • 精度：これにより、単一のシミュレーションを用いて、$k_{\text{min}} = 0.04$ から $k_{\text{max}} = 0.47 \, h\text{Mpc}^{-1}$ の範囲にあるすべての三角形構成に対して、2% 未満の精度で物質ビスペクトルを測定することが可能となった。
• 解析的扱いやすさ：著者らは、シフト制御変量は原理的に任意の次数まで解析的に計算可能であることを示したが、完全なビスペクトルに対する効率的な数値実装は今後の課題として残されている。現在、彼らはその性能を数値的に推定している。

意義と主張
本論文は、正確なビスペクトルエミュレータの開発が、数千回の N 体シミュレーションの平均化という莫大なコストなしに可能であると主張している。シフト制御変量を利用することで、シミュレーションベースのモデリングに内在する分散は、標本分散ではなく、小規模な三角形構成の測定によってのみ精度が制限されるレベルまで抑制可能である。

著者らは、この研究を制御変量を $n$ 点階層全体へ拡張するための理論的基盤を築くものとして位置づけている。彼らは明示的に、この手法はバイアストレーサ（銀河）および赤方偏移空間への拡張に十分であることを述べているが、これらの応用は今後の課題として残されている。本論文は、完全なシフトビスペクトルの解析的計算を解決したとは主張していないが、その理論的な妥当性と数値テストにおける卓越した性能を確立している。主な意義は、今後の大規模構造サーベイに適した、高精度なシミュレーションベースの宇宙論的モデリングを可能にする点にある。