Excursion-set for Primordial Black Holes I: white noise and moving barrier

この論文は、初期ブラックホールの質量分布を導出する際のエクスカーションセット形式が、ホワイトノイズと移動バリアを用いることで色ノイズやクラウド・イン・クラウド効果に関する最近の批判を克服し、現実的なシナリオにおいて堅牢かつ不可欠であることを示しています。

要約

🌌 物語の舞台：宇宙の「嵐」と「ブラックホール」

昔々、宇宙が生まれた直後、空間には「密度のむら（揺らぎ）」という嵐が吹き荒れていました。
この嵐が非常に激しい場所では、重力が暴走して**「原始ブラックホール」**という小さなブラックホールが生まれました。

科学者たちは、「どのくらいの大きさのブラックホールが、どれくらいたくさんできるのか？」を知りたがっています。これを調べるために使われてきたのが**「エクスカージョン・セット（Excursion Set）」**という計算方法です。

🎮 従来のルールと、そこに潜む「バグ」

この計算方法は、**「ランダムウォーク（酔っ払いの歩き方）」**というゲームに例えられます。

1. ゲームの仕組み：

   • 酔っ払い（密度のむら）が、ランダムに足踏みしながら歩きます。
   • 地面に**「壁（バリア）」**が立っています。
   • 酔っ払いがその壁を初めて越えた瞬間、そこで**「ブラックホールが誕生した！」**と判定されます。
   • 壁を越えるまでの歩数（距離）によって、ブラックホールの大きさが決まります。

2. これまでの問題点（バグ）：

   • 問題①：「色付きのノイズ」の誤解
     以前、このゲームを「ハッブル半径（宇宙の見える範囲）」という基準で計算すると、酔っ払いの動きが「過去の影響を受けすぎていて、予測不能（色付きノイズ）」だと指摘されました。これだと計算が複雑すぎて、正解が出せません。
     • この論文の解決策： 「ハッブル半径」ではなく、**「同じ時刻（同期した瞬間）」で観察すれば、酔っ払いの動きは単純な「白いノイズ（完全なランダム）」になります。つまり、「見るタイミングを変えるだけで、計算が劇的にシンプルになる」**という発見です。
   • 問題②：「雲の中の雲」の無視
     大きなブラックホールができると、その中に小さなブラックホールが含まれてしまうことがあります（雲の中に雲がある状態）。
     以前の計算では、「大きなものができたら、中の小さなものはカウントしない」というルール（クラウド・イン・クラウド）を無視する人がいましたが、「パワーが広い（バラエティに富んだ）嵐」の場合、このルールを無視すると計算が破綻し、物理的にありえない「負の数（マイナスのブラックホール）」が出てきてしまうことがわかりました。

🔧 新しい計算方法：「動く壁」と「確率の迷路」

この論文では、上記の問題を解決するために、以下の新しいアプローチを提案しています。

1. 「動く壁」を相手にする

従来の計算では、壁（ブラックホールになるための条件）は固定されていると仮定していました。しかし、正しい計算方法（同期した瞬間で見る）を使うと、壁の高さが「距離（スケール）」によって動いてしまうことがわかりました。

• 例え： 酔っ払いが歩いている間、ゴールの壁が近づいたり遠ざかったりしている状態です。
• 解決策： 著者たちは、この「動く壁」を相手にするための**「ボルテラ積分方程式」**という、効率的な数学のツールを開発しました。これを使えば、コンピューターで高速かつ正確にシミュレーションできます。

2. 「雲の中の雲」の重要性を再確認

「大きなブラックホールが小さなものを飲み込む」という現象（クラウド・イン・クラウド）は、「嵐の広がり（パワースペクトル）が広い場合」には無視できないと証明しました。

• 例え： 小さな雨粒（小さなブラックホール）が、大きな雨粒（大きなブラックホール）に飲み込まれて消えてしまう現象です。
• 発見： 嵐が広範囲にわたって激しい場合、この「飲み込み」が頻繁に起き、結果として**「小さなブラックホールの数が減り、大きなものが相対的に増える」**という、直感とは異なる結果が生まれます。これを無視すると、現実とかけ離れた予測をしてしまいます。

📊 具体的な成果：どんなブラックホールが生まれる？

この新しい計算方法を使って、さまざまな「嵐のパターン（パワースペクトル）」をシミュレーションしました。

• 狭いピークの場合： ほぼ同じ大きさのブラックホールが大量に生まれます（単一のサイズ）。
• 広いピークの場合： 非常に小さなものから巨大なものまで、多様なサイズのブラックホールが混在します。
• 二つのピークの場合： 2 つの異なるサイズにピークがある「二峰性」の分布になります。

特に驚くべきは、従来の方法（プレス・シェクター近似）では「負の数」が出てしまうようなケースでも、この新しい方法を使えば**「常に正の値（現実的な数）」**が得られることです。

🎯 まとめ：なぜこの論文は重要なのか？

1. 計算の誤りを正した： 「色付きノイズ」は計算のタイミングの間違いだったことを証明し、正しい「白いノイズ」の計算方法を確立しました。
2. 新しい計算ツールを提供： 「動く壁」を扱うための効率的な数値計算フレームワークを公開しました。
3. 「飲み込み」現象の重要性： 広い範囲でブラックホールが生まれるシナリオでは、「大きなものが小さなものを飲み込む」現象を無視できないことを示しました。

一言で言うと：
「宇宙のブラックホールの誕生シミュレーションにおいて、これまで見逃されていた『計算のタイミング』と『飲み込み効果』を正しく組み込むことで、より現実的で正確な予測ができるようになった」という画期的な研究です。

これにより、重力波観測やダークマターの研究など、現代宇宙論の重要な分野において、より信頼性の高いデータが得られるようになるでしょう。

技術概要

論文「Excursion-set for Primordial Black Holes I: white noise and moving barrier」の技術的サマリー

本論文は、原始ブラックホール（PBH）の形成と質量分布を記述するための「遊走集合（Excursion-set）」定式化に関する最近の批判に答えるとともに、その手法の堅牢性と必要性を再確認するものである。著者らは、PBH 形成シナリオにおいて、ランダムウォークのノイズ特性とバリアの時間依存性を正しく扱うことが不可欠であることを示し、移動バリアを有する初回通過時間（First-Passage Time）問題に対する効率的な数値枠組みを提案している。

以下に、問題提起、手法、主要な貢献、結果、そして意義について詳述する。

1. 背景と問題提起

PBH の質量分布を予測する際、遊走集合アプローチは密度揺らぎが閾値を越える「初回通過時間」をランダムウォークとして扱うことで、階層的構造形成（雲中雲問題：雲の中にさらに小さな雲が含まれる現象）を自動的に考慮できる利点を持つ。しかし、最近の研究 [29] において、このアプローチに対して以下の 2 点の批判がなされていた。

1. 有色ノイズとドリフトの問題:
   • 従来の手法では、ハッブル半径を横断する面（Hubble-crossing surface）上でサンプリングを行うため、密度揺らぎのランダムウォークが時間相関を持つ「有色ノイズ」にさらされ、ドリフト項が存在すると主張されていた。これにより、初回通過時間問題が複雑化し、物理的にあり得ない「負の質量分布」が得られる可能性が指摘されていた。
2. 雲中雲（Cloud-in-Cloud）の無視:
   • 別の研究 [31, 32] では、PBH の形成スケールが十分に離れている場合、雲中雲効果は無視できると主張されていた。これにより、遊走集合の複雑な計算（境界の進化を考慮する必要）が不要であるという見解が示唆されていた。

2. 手法と理論的枠組み

著者らは、上記の問題がサンプリング面の選び方とバリアの扱いに起因することを明らかにし、以下の理論的修正と数値手法を提案した。

A. サンプリング面の再定義：同期面（Synchronous Surface）

• ハッブル横断面の問題点: ハッブル半径を横断する面（$R = H^{-1}$）でサンプリングすると、モード関数の時間進化が「ドリフト」として現れ、かつスケール $R$ と分散 $S$ の関係が単調でなくなる場合がある（これが負の質量分布の原因）。また、このアプローチではノイズが有色に見えるが、これはドリフトとノイズの分離が不適切であるためである。
• 同期面の採用: 著者らは、一定の時間 $t=t_0$ を持つ「同期面（synchronous hypersurface）」上でサンプリングを行うことを提案する。
  • この場合、ドリフト項が完全に消滅し、ランダムウォークは純粋な「白色ノイズ（white noise）」に従うランジュバン方程式となる。
  • 分散 $S$ とスケール $R$ の関係が常に単調減少となるため、時間変数の再定義が常に可能であり、負の質量分布の問題が解消される。
• 代償としての移動バリア: 同期面を採用する代償として、PBH 形成の閾値（バリア）$\delta_c$ が、粗視化スケール $R$（または分散 $S$）に依存する「移動バリア」となる。

B. 移動バリアを持つ初回通過時間問題の数値解法

移動バリアを持つブラウン運動の初回通過時間分布を解析的に解くことは困難であるため、著者らは以下の数値手法を採用した。

• ボルテラ積分方程式（Volterra Integral Equations）: 吸収境界条件を持つブラウン運動の初回通過時間確率分布 $P_{\text{FPT}}$ を記述する、第 2 種のボルテラ積分方程式を導出した。
• 核関数の最適化: 数値的不安定性を除去するため、特異点を除去する特定の核関数 $K(S) = \delta_c'(S)$ を選択した。
• 行列方程式への離散化: 積分方程式を離散化し、下三角行列を用いた線形方程式系 $P_{\text{FPT}} = (I - J\Delta s)^{-1} P_{\text{PS}}$ として解く効率的なアルゴリズムを構築した。
  • この手法は、標準的なモンテカルロシミュレーションに比べて計算コストが低く、収束が速いことが確認された。

3. 主要な貢献と結果

A. 有色ノイズと負の質量分布の解決

• ハッブル横断面でのサンプリングが「ドリフトとノイズの誤った分離」および「$R$ と $S$ の非単調性」を引き起こすことを示し、同期面への切り替えですべての病理が解消されることを証明した。
• 移動バリアを考慮しない従来のプレス・シュヒター（Press-Schechter）近似では、バリアが急峻な場合に負の質量分布が予測されるが、遊走集合アプローチ（移動バリアを正しく扱う）では常に正の確率分布が得られることを示した。

B. 雲中雲効果の再評価

• 広範なパワースペクトルにおける重要性: 雲中雲効果は、ブラックホールの形成スケールが十分に離れている場合（狭いピークを持つパワースペクトル）には無視できるが、**広範なパワースペクトル（Broad power spectra）**や連続的なモードが強化されている場合には無視できないことを示した。
• 質量分布への影響: 雲中雲効果を考慮すると、低質量側の PBH 数が抑制され、高質量側の PBH 数が増加する（大きなブラックホールが小さなものを飲み込むため）。
• 非単調な振る舞い: 2 つのピークを持つパワースペクトル（ダブル・ログノーマル）を用いたシミュレーションでは、大規模ピークの振幅を増大させた際、低質量 PBH の存在割合が最初は増加し、その後減少するという非単調な振る舞いを示すことがわかった。これは、初期条件の拡散と、大規模ピークでのバリア越えによる「雲中雲」効果の競合によるものである。

C. 具体的なパワースペクトルモデルへの適用

• トプハット型、ログノーマル型、ダブル・ログノーマル型のパワースペクトルに対して、提案された枠組みを適用し、PBH の質量分布を計算した。
• 狭いピークでは質量分布が単色的になるが、広いピークでは数桁にわたる質量分布が得られることを確認した。
• 従来の定数バリア近似やプレス・シュヒター近似では得られない、移動バリアの効果を定量的に評価した。

4. 意義と結論

本論文の成果は、PBH 形成研究における遊走集合アプローチの信頼性を確立するものである。

1. 手法の正当化: 移動バリアを伴う初回通過時間問題を効率的に解く数値枠組みを提供し、有色ノイズや負の質量分布といった以前の問題が、サンプリング面の不適切な選択に起因していたことを解明した。
2. 物理的洞察: 広範なパワースペクトルを持つ現実的なシナリオにおいて、雲中雲効果が PBH の質量分布の形状（特に低質量側の抑制）に決定的な役割を果たすことを示した。
3. 将来展望: プレス・シュヒター近似の限界を明らかにし、より正確な PBH 質量分布の予測には、移動バリアを考慮した遊走集合アプローチが不可欠であることを示唆した。

今後は、非ガウス性や臨界スケーリング、圧縮関数（compaction function）の導入など、より現実的な PBH 形成条件への拡張が期待される。本論文は、そのための堅固な基礎を提供している。