COSMOS: A numerical relativity code specialized for PBH formation

COSMOSは、特殊な非デカルト座標系によるスケールアップと、非線形重力ダイナミクスを扱うための固定メッシュリファインメントメントを用いて、3+1次元におけるアインシュタイン方程式を解くことにより原始ブラックホールの形成をシミュレートするために設計された、スタンドアロンのOpenMP並列化されたC++数値相対論コードである。

要約

ビッグピクチャー：ゼロからのブラックホール構築

初期宇宙を、穏やかで滑らかな海ではなく、巨大な波がうねる荒れた海だと想像してみてください。時として、これらの波はあまりに高く重くなりすぎ、自らの重みで崩壊し、**原始ブラックホール（PBH）**を形成します。私たちが普段耳にするブラックホール（死んだ星の残骸）とは異なり、PBHはビッグバンの直後、時空の「しわ」から瞬時に誕生しました。

科学者たちは、これらのブラックホールがどのように形成されるかを理解したいと考えていますが、その数学は非常に困難です。それは、水たまりに水滴が落ちたときに、具体的にどのように水しぶきが上がるかを予測しようとするようなものです。ただし、その水たまりは膨張しており、しかも水自体が純粋な重力でできているのです。

解決策：COSMOS（デジタル実験室）

この論文では、これらのブラックホールの誕生をシミュレートするために特別に設計されたコンピュータプログラム（C++で記述）であるCOSMOSを紹介しています。COSMOSを、重力の**「ハイテクなデジタル風洞」**だと考えてください。エンジニアが車のスケールモデルを作って空気の流れをテストするように、物理学者はCOSMOSを使って初期宇宙の「スケールモデル」を作り、物事が混沌とした状況にあるときに重力がどのように振る舞うかを確認します。

仕組み：「スマートズーム」のトリック

このシミュレーションにおける最大の課題の一つは、2つの非常に異なるサイズが同時に存在していることです。

1. 極小の領域： ブラックホールが崩壊している特定の領域（非常に小さい）。
2. 大きな全体像： その周囲で膨張している宇宙全体（非常に大きい）。

もし顕微鏡で宇宙全体を見ようとすれば、全体像を見失います。逆に、広角レンズで宇宙全体を見ようとすれば、崩壊の微細なディテールを見ることができません。

COSMOSは、「スマートズーム」機能によってこれを解決します。
崩壊する星の映画を見ていると想像してください。画面の大部分は広大な空虚な宇宙を示しています。しかし、星が縮み始めると、カメラは自動的にその一点へと超至近距離までズームインし、必要な場所に正確に「ピクセル（解像度）」を追加していきます。これにより、コンピュータは、周囲の空っぽな空間を計算するために都市ひとつ分のようなスーパーコンピュータを必要とすることなく、小さく激しい崩壊を扱うことができるのです。

材料：シミュレーションの中身は？

シミュレーションを実行するために、COSMOSは、パン職人が特定の種類の生地を作るように、2つの主要な材料を混ぜ合わせます。

1. 流体： 初期宇宙の物質を表す、完璧で滑らかな流体。これは単純なルール（風船が膨らんだり萎んだりするようなルール）に従います。
2. スカラー場： 空間を波のように伝わる、幽霊のようなエネルギー場。

プログラムは、これらの材料がどのように相互作用するかを見るために、アインシュタイン方程式（重力のルールブック）を解きます。そしてこう問いかけます。「もし宇宙に特定の『コブ』があった場合、それは滑らかに消えるのか、それとも自らの重みで押しつぶされてブラックホールになるのか？」

「準備不要」のメリット

通常、物理シミュレーションのセットアップは、オーブンをつける前に複雑な数学方程式を解いて「生地をあらかじめ混ぜておく」必要があるケーキ作りによく似ています。

COSMOSは異なります。すでに流体で満たされた宇宙をシミュレートするため（空っぽの空間ではなく）、「最初から生地が混ざった状態」でスタートできます。初期条件が完璧に設定されているため、コンピュータは開始時に難しい「楕円型」の方程式を解くために時間を費やす必要がありません。ただ「実行」ボタンを押し、物語が展開するのを眺めるだけです。これにより、コードはより軽量で高速になり、他の科学者がインストールして使用することも容易になります。

論文が示すもの（例題）

論文には、コードがどのように機能するかを示す3つの「試運転」が含まれています。

1. 単純な波： 既知の単純な数学と一致することを示すために、空間の小さな穏やかな波紋をシミュレートします。
2. 完全な球体（断熱的）： 完璧に丸い、空間の泡の崩壊をシミュレートします。論文では、「ラプス関数」（時間の流れの尺度）の図と、ブラックホールが形成されるにつれてコンピュータが中心部へどのようにズームしていくかが示されています。
3. ゴースト・ウェーブ（等曲率）： 先ほど述べた「幽霊のような」エネルギー場によって引き起こされる崩壊をシミュレートします。

これらのケースすべてにおいて、コードはブラックホールが誕生した瞬間（技術的には「見かけの地平線」と呼ばれます）を正常に検出し、その形状をマッピングすることに成功しています。

なぜこれが重要なのか

著者たちは単なるおもちゃを作っているのではありません。彼らは特定の仕事のために設計された「専用のレンチ」を作っているのです。一般的なブラックホールの衝突（合体する星など）のためのツールは他に存在しますが、COSMOSは、初期宇宙のユニークで混沌とした、かつ膨張する環境のために特別に設計された**「専用のレンチ」**なのです。

このコードをオープンで使いやすいものにすることで、著者たちは他の科学者がこれを拡張し、「ダークマターは何でできているのか」や「なぜ今日、重力波が観測されるのか」といった、宇宙の隠された歴史に関する新たな問いを投げかけられるようになることを期待しています。

---

要約すると： この論文は、初期宇宙のための「スマートズーム」カメラとして機能する、特殊なコンピュータプログラムであるCOSMOSを提示しています。これにより、科学者は時空の小さな「しわ」がどのように崩壊して最初のブラックホールを形成するかを、計算能力を必要な場所に集中させることで、効率的にシミュレートできるようになります。

技術概要

COSMOSの技術概要：原始ブラックホール形成のための数値相対論コード

問題提起
原始ブラックホール（PBH）は、初期宇宙において、非線形に大きな初期振幅を持つ超地平線的な原始揺らぎから、恒星進化を経ることなく形成されると仮定されている。その形成をシミュレートするには、崩壊領域の微視的なスケールと、宇宙膨張のマクロなスケールの両方が共存する条件下で、アインシュタイン方程式を解く必要がある。標準的な数値的手法では、これらの極端に異なるスケールを効率的に解像するための計算コストが課題となることが多い。さらに、漸近的に平坦な真空領域ではなく、流体とスカラー場で満たされた時空における初期データの取得は、拘束方程式に関して特有の困難を伴う。

手法
本論文では、PBH形成のために特別に設計された、アインシュタイン方程式を3+1次元で解くためのC++パッケージであるCOSMOSを紹介する。このコードは、もともと異なる目的のために開発されたSACRAコードのC++への移植および発展版である。

• 物理モデル: COSMOSは、物質を線形状態方程式（$P = w\rho$）を持つ完全流体と無質量スカラー場としてモデル化する。エネルギー・運動量テンソルは、流体（$T^{FL}_{\mu\nu}$）とスカラー場（$T^{SC}_{\mu\nu}$）の寄与の和である。系は、アインシュタイン方程式（$G_{\mu\nu} = 8\pi G/c^4 T_{\mu\nu}$）、スカラー場の方程式（$\nabla_\mu\nabla^\mu\phi = 0$）、および流体の保存方程式（$\nabla_\mu T^{FL}_{\mu\nu} = 0$）を解く。流体の進化は、Kurganov & Tadmor (2000) および Shibata & Font (2005) によって確立されたスキームに従う。
• 数値技術: マルチスケール問題に対処するため、COSMOSは以下を実装している：
  • 非デカルト座標によるスケールアップ: 崩壊領域を効率的に解像するため。
  • 固定メッシュ細分化 (FMR): 必要に応じて解像度を動的に調整するため。
  • 並列化: 実用的な計算速度を達成するためにOpenMPを利用。
  • 依存関係: インストールや拡張を容易にするため、外部依存関係を持たない設計となっている。
• 初期データ構築: 漸近的に平坦な真空時空のための標準的な手法（楕円型微分方程式の求解を必要とする）とは異なり、COSMOSは長波長成長モード解（膨張パラメータ $\epsilon = k/(aH) \ll 1$ における次項までの近似）を採用することで初期データを生成する。
  • 断熱的揺らぎの場合、初期データは曲率揺らぎ $\zeta(\vec{x})$ によって特徴付けられる。
  • 等曲率揺らぎの場合、無質量スカラー場に関連する $\Upsilon(\vec{x})$ によって特徴付けられる。
  • 流体分布は、楕円型ソルバーを必要とせずに、機械精度内でアインシュタインの拘束方程式を満たすように生成される。

主な貢献

1. 特化したソフトウェアの公開: 本論文は、一般的な数値相対論コードとは一線を画す、PBH形成のシミュレーションに特化したツールであるCOSMOSの公開を発表するものである。
2. 効率的な初期データ生成: 楕円型ソルバーによる宇宙論的流体背景における初期データの計算オーバーヘッドを回避するため、長波長近似を通じて解析的に拘束方程式を満たす手法を実装している。
3. ドキュメントと例題: パッケージには、教育目的の3つのデモンストレーション例が含まれている：
   • 単一モードの正弦波摂動の進化（線形摂動論との比較検証）。
   • 断熱的な球対称初期揺らぎ（見かけの地平線の形成の提示）。
   • 無質量スカラー場によって駆動される球対称な等曲率揺らぎ。
   • これらの例題は、高精度な結果を示すためではなく、機能性を実証するために低解像度を使用している。

結果
論文では、含まれる例題からの視覚的および定性的な結果を提供している：

• 線形検証: 小さな正弦波揺らぎに対する外積の跡（$trK$）の時間進化は、線形摂動方程式の解と一致することが示されている。
• 地平線形成: 断熱的および等曲率摂動のシミュレーションは、見かけの地平線の形成を正常に示している。コードはデータファイル（例：ini_all.dat）を出力し、ユーザーが見かけの地平線が検出された瞬間の幾何学構造と物質分布を再構成することを可能にする。
• メッシュ細分化: 固定メッシュ細分化レイヤーの動作を確認する可視化により、崩壊領域の周囲に明確な解像度領域（第0、第1、第2レイヤー）が存在することが示されている。

意義と主張
著者らは、特に重力波観測の文脈におけるPBHへの関心の高まりを支えるための特化したツールとしてCOSMOSを位置づけている。本論文は、以下の点を通じて、初期宇宙における非線形重力ダイナミクスのシミュレーションに対する合理的なアプローチを提供すると主張している：

• ソースコードを理解すれば、様々な科学的設定に対して容易にインストールおよび拡張可能な、自己完結型のC++環境を提供すること。
• 楕円型ソルバーを用いずに、宇宙論的流体に関する特定の初期値問題を扱う堅牢な手法を提供すること。
• 著者らおよび共同研究者による数多くの出版物（例：Yoo et al., 2013; Okawa et al., 2014; Escrivà & Yoo, 2025）に示される通り、過去および将来の研究の基礎として機能すること。

論文では、含まれる例題はデモンストレーション用であり、解像度は意図的に低く抑えられていることが謙虚に記されている。また、著者らの過去の出版物で使用された一部の高度な機能は、公開版には含まれていない可能性があることを明記しており、過去の結果をすべて再現するには、非公開の追加実装が必要になる可能性があることを示唆している。