ブラックホール衝突計算機
連星ブラックホールの合体時間、重力波エネルギー、最終ブラックホールの性質を計算します。
2つのブラックホールの質量、初期軌道間隔、離心率、傾斜角を入力すると、Peters の一般相対論式を用いてインスパイラルと合体のパラメータを計算できます。
ブラックホール衝突計算機
連星ブラックホールの合体時間、重力波エネルギー、最終ブラックホールの性質を計算します。
ブラックホール衝突計算機について
2つのブラックホールが重力的に束縛された連星系を作ると、重力波としてエネルギーと角運動量を放出しながら、ゆっくり互いに螺旋状に近づいていきます。この過程はインスパイラルと呼ばれ、1964 年に Peter Peters がアインシュタインの一般相対論を用いて初めて導いた数式の枠組みに従います。ブラックホール衝突計算機は Peters の式を実装し、連星が合体するまでの時間、放射するエネルギー、そして最終的な合体天体の姿を見積もります。
重力波天文学で最も重要な導出量はチャープ質量で、M_chirp = (m₁ m₂)^(3/5) / (m₁+m₂)^(1/5) と定義されます。チャープ質量は重力波周波数がどれだけ速く上昇するか、つまり低周波から高周波へ掃き上がる特徴的な変化を決めます。対称質量比 η = m₁m₂/(m₁+m₂)² と合わせることで、一次近似のインスパイラル動力学を計算するのに必要な情報をすべて表します。
Peters（1964）が示した円軌道の合体時間式は T = (5/256) × c⁵ a⁴ / (G³ m₁ m₂ M_total) で、a は初期長半径、G は重力定数、c は光速です。離心軌道では合体時間が概ね (1−e²)^(7/2) の因子で短くなります。つまり、離心率が高い連星ほど、同じ初期距離の円軌道より効率よくエネルギーを放射し、より速く合体します。この近似は e ≲ 0.6 で特に有効です。
インスパイラル中に放射される重力波エネルギーは、換算質量エネルギー（μc²）の約 5% と見積もられます。これは同程度の質量比を持つ連星の数値相対論シミュレーションと整合的です。残りの質量が最終的なブラックホールになり、そのシュワルツシルト半径は r_s = 2 G M_final / c² です。最内安定円軌道（ISCO）でのピーク重力波周波数は f_peak = c³ / (π × 6√6 × G × M_total) で、インスパイラルから急降下、そしてリングダウンへの移行点を示します。これは合体で最も大きく、最も検出しやすい瞬間です。
重力波の初の直接検出 GW150914 は、2015 年 9 月 14 日に LIGO により観測されました。これはおよそ 36 太陽質量と 29 太陽質量の2つのブラックホールが、13 億光年先で合体したものです。この事象では、約 3 太陽質量分のエネルギーが 1 秒にも満たない短時間で重力波として放射され、重力波光度では観測可能宇宙全体を一時的に上回りました。その後、LIGO–Virgo–KAGRA 共同観測は 90 件以上の連星合体事象を検出し、重力波天文学を精密科学へと変えました。
ブラックホール衝突の例
下の表は、代表的な連星ブラックホール系と主要な合体パラメータを示します。
| 系のパラメータ | 主要結果 | 事象 / 背景 |
|---|---|---|
| m₁=36 M☉, m₂=29 M☉, a=10,000,000 km, e=0 | T_merge ≈ 94.4 Myr, M_chirp ≈ 28.1 M☉, f_peak ≈ 67.6 Hz | GW150914 に類似（LIGO, 2015） |
| m₁=1000 M☉, m₂=800 M☉, a=100,000,000 km, e=0.3 | T_merge ≈ 32.0 Myr, M_chirp ≈ 778 M☉, f_peak ≈ 2.44 Hz | 中間質量ブラックホール連星 |
| m₁=20 M☉, m₂=20 M☉, a=5,000,000 km, e=0 | T_merge ≈ 25.0 Myr, M_chirp ≈ 17.4 M☉, f_peak ≈ 110 Hz | 等質量の恒星質量連星 |
ブラックホール衝突計算機の使い方
- 各ブラックホールの質量を太陽質量 (M☉) で入力します。恒星質量ブラックホールはおよそ 3〜100 M☉、超大質量ブラックホールは 10⁹ M☉ を超えることがあります。
- 初期軌道間隔を km で入力します。これは連星軌道の初期長半径です。
- 軌道離心率を 0（円軌道）から 0.99（ほぼ放射状）まで設定します。LIGO が検出する多くの事象は、検出帯域に入る頃にはほぼ円軌道になっています。
- 傾斜角を度で入力します（0° = 正面向き、90° = 斜め向き）。これは地球で観測される重力波振幅に影響しますが、合体時間には影響しません。
- 計算をクリックすると、チャープ質量、推定合体時間、重力波エネルギー、最終ブラックホール質量、シュワルツシルト半径、ピーク重力波周波数が表示されます。
よくある質問
チャープ質量とは何ですか。なぜ重要なのですか。
チャープ質量 M_chirp = (m₁m₂)^(3/5)/(m₁+m₂)^(1/5) は、重力波検出で最も重要な単一パラメータです。重力波周波数が上昇する速さ（チャープ率）を決めるため、個々の質量が分からなくても波形から非常に高精度に推定できます。
合体時間の推定精度はどの程度ですか。
ここで使う Peters の式は、分離がシュワルツシルト半径より十分大きい初期インスパイラルに対して有効です。離心率補正 (1−e²)^(7/2) は近似で、e ≲ 0.6 ではよく機能します。離心率が高い軌道や非常に近い分離では、正確な推定に数値相対論が必要です。
なぜ離心率が高いほど合体が速くなるのですか。
離心軌道の最接近点（近点）では、天体が最も速く、最も近くなり、その瞬間の重力波放射パワーが大きく増加します。1 周あたりに失うエネルギーが増えるため軌道はより速く縮み、同じ平均分離の円軌道より合体時間が短くなります。
ISCO とは何ですか。なぜピーク周波数を決めるのですか。
最内安定円軌道（ISCO）は、シュワルツシルト（非自転）ブラックホールの周囲で安定な円軌道が存在しなくなる境界です。インスパイラルがこの点に達すると、小さい天体は急速に落下します。ISCO の軌道周波数を重力波周波数に換算したものが、インスパイラル信号の最高周波数であり、合体リングダウンの始まりを示します。
重力波としてどれくらいのエネルギーが放出されますか。
同程度の質量のブラックホール合体では、数値相対論シミュレーションにより総質量の約 4〜8% が重力波として放射されます。この計算機では、換算質量エネルギーの約 5% を概算値として用います。GW150914 では、約 3 太陽質量（総質量の約 5%）が 1 秒未満で重力波エネルギーに変換されました。
この計算機は中性子星合体にも使えますか。
このインスパイラル式は、中性子星–中性子星（BNS）や中性子星–ブラックホール（NSBH）連星にも同様に使えます。ただし BNS では、潮汐破壊や中性子星の状態方程式による補正がここには含まれていません。概算には使えますが、正確な BNS の結果にはポスト・ニュートンまたは数値相対論の波形モデルを使ってください。