블랙홀 충돌 계산기
쌍성 블랙홀의 병합 시간, 중력파 에너지, 최종 블랙홀 성질을 계산합니다.
두 블랙홀의 질량, 초기 궤도 분리, 이심률, 경사각을 입력하면 Peters의 일반상대론 공식을 사용해 나선형 나선접근과 병합 매개변수를 계산할 수 있습니다.
블랙홀 충돌 계산기
쌍성 블랙홀의 병합 시간, 중력파 에너지, 최종 블랙홀 성질을 계산합니다.
블랙홀 충돌 계산기 소개
두 블랙홀이 중력적으로 묶인 쌍성계를 이루면, 중력파로 에너지와 각운동량을 방출하면서 점차 서로 나선형으로 가까워집니다. 이 과정은 나선접근(inspiral)이라 불리며, 1964년 Peter Peters가 아인슈타인의 일반상대론을 바탕으로 처음 유도한 수학적 틀을 따릅니다. 블랙홀 충돌 계산기는 Peters의 공식을 구현해 쌍성이 병합하는 데 걸리는 시간, 방출하는 에너지, 그리고 최종 병합체의 모습을 추정합니다.
중력파 천문학에서 가장 중요한 유도량은 처프 질량이며, M_chirp = (m₁ m₂)^(3/5) / (m₁+m₂)^(1/5)로 정의됩니다. 처프 질량은 중력파 주파수가 얼마나 빠르게 증가하는지, 즉 신호가 낮은 주파수에서 높은 주파수로 휩쓸리듯 올라가는 특징을 결정합니다. 대칭 질량비 η = m₁m₂/(m₁+m₂)²와 함께, 처프 질량은 1차 나선접근 동역학을 계산하는 데 필요한 모든 정보를 담고 있습니다.
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 태양질량인 두 블랙홀이 13억 광년 거리에서 병합한 사건이었습니다. 이 사건은 1초도 채 안 되는 시간에 약 3 태양질량에 해당하는 에너지를 중력파로 방출해, 중력파 광도 기준으로 잠시 관측 가능한 우주 전체를 능가했습니다. 이후 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에서 잘 작동합니다. 이심률이 매우 높거나 분리가 매우 작은 궤도는 정확한 추정을 위해 수치상대론이 필요합니다.
왜 이심률이 높을수록 병합이 빨라지나요?
이심 궤도의 가장 가까운 접근점(근일점)에서는 천체가 가장 빠르고 가장 가까워지며, 그 순간 방출되는 중력파 파워가 크게 증가합니다. 궤도 한 바퀴당 방출되는 에너지가 더 많아져 궤도가 더 빨리 줄어들고, 같은 평균 분리를 가진 원궤도보다 병합 시간이 짧아집니다.
ISCO는 무엇이며 왜 최대 중력파 주파수를 정의하나요?
가장 안쪽 안정 원궤도(ISCO)는 슈바르츠실트(비회전) 블랙홀 주위에 더 이상 안정적인 원궤도가 존재하지 않는 경계입니다. 나선접근이 이 지점에 도달하면 작은 천체는 급격히 안쪽으로 떨어집니다. ISCO의 궤도 주파수를 중력파 주파수로 두 배 변환한 값이 나선접근 신호의 최고 주파수이며, 병합 링다운의 시작을 나타냅니다.
중력파로 얼마나 많은 에너지가 방출되나요?
질량이 비슷한 블랙홀 병합에서는 수치상대론 시뮬레이션이 총질량의 약 4~8%가 중력파로 방출된다고 보여줍니다. 이 계산기는 감소 질량 에너지의 약 5%를 대략적인 값으로 사용합니다. GW150914에서는 약 3 태양질량(총질량의 약 5%)이 1초도 안 되어 중력파 에너지로 바뀌었습니다.
이 계산기를 중성자별 병합에도 사용할 수 있나요?
이 나선접근 공식은 중성자별–중성자별(BNS)과 중성자별–블랙홀(NSBH) 쌍성에도 마찬가지로 적용됩니다. 다만 BNS 사건에서는 조석 파괴와 중성자별 상태방정식이 추가 보정을 도입하며, 여기에는 포함되어 있지 않습니다. 대략적인 추정에는 사용할 수 있지만, 정확한 BNS 결과는 포스트 뉴턴 또는 수치상대론 파형 모델을 사용해야 합니다.