블랙홀‑토러스 시스템의 파팔루지‑프링클리 불안정이 만든 중력파

초록

3차원 일반상대론적 시뮬레이션을 통해 블랙홀 주변 토러스에서 m=1 파팔루지‑프링클리 불안정(PPI)이 넓은 매개변수 범위에서 성장하고, 포화 후에도 장시간 유지되며 강한 준주기 중력파를 방출한다는 것을 확인하였다. 질량 10 M⊙(거리 100 Mpc)와 10⁶ M⊙(거리 10 Gpc) 규모의 시스템에 대해 각각 LIGO‑VIRGO와 LISA 탐지 가능성을 제시한다.

상세 분석

본 연구는 블랙홀‑토러스 시스템을 초기조건으로 삼아, 고정 메쉬 정밀화(Fixed Mesh Refinement, FMR)와 SACRA 코드 기반의 3차원 일반상대론 수치 시뮬레이션을 수행하였다. 블랙홀은 무자이 회전(스핀 0)으로 설정하고, 토러스는 자기 중력이 중요한 질량비 R = M_tor/M_BH = 0.06–0.10 범위에서 조사하였다. 각 모델은 두 종류의 각운동량 분포를 사용했는데, (i) j = const(‘C’ 모델)와 (ii) j ∝ Ω^(-1/6) 형태의 비정상적인 분포(‘NC’ 모델)이다. 방정식 상태는 Γ = 4/3 폴리트로프(EOS)로, 방사선 혹은 고에너지 전자 가스 상황을 근사한다.

시뮬레이션 결과, m = 1 비축대칭 모드가 대부분의 모델에서 가장 빠르게 성장함을 확인하였다. 성장률 Im(ω₁) ≈ 5–25 % Ω_c 로, 토러스 질량이 클수록, 그리고 j = const인 경우에 더 크게 나타났다. 반면, NC06 모델처럼 각운동량이 외부로 급격히 감소하는 경우에는 PPI가 억제되어 불안정이 발생하지 않았다. 성장 단계에서는 토러스 중심밀도가 거의 일정하게 유지되지만, 비선형 포화 시점에 최대 밀도가 급증하고 이후 새로운 준평형 상태로 수렴한다. 포화 후에도 m = 1 구조가 지속적으로 유지되며, 이는 블랙홀-토러스 전체의 질량중심이 블랙홀 중심에 고정되지 못하고 비대칭적인 질량분포를 형성하기 때문이다.

이 비대칭 구조는 중력파 방출에 직접적인 원천이 된다. 파동의 주요 모드는 (l,m) = (2,2)이며, 이는 토러스 내부 물질이 블랙홀 주위를 비대칭적으로 회전하면서 발생한다. 파동 신호는 포화 직후 급격한 ‘버스트’와 그 후 수십 회에 걸친 준주기 진동으로 구성된다. GW의 유효 진폭 h_eff는 h_eff =