중성자별 급격 붕괴와 디스크 잔류: 전일반상대론 3차원 시뮬레이션

초록

이 연구는 주변에 약 7 % 질량을 가진 얇은 원반을 가진 중성자별이 압력 감소에 의해 붕괴해 회전 블랙홀을 형성하는 과정을 전일반상대론 3차원 시뮬레이션으로 최초로 조사한다. 원반이 없는 경우와 비교했을 때 중력파 파형은 거의 동일하지만, 원반이 남아 블랙홀 주변에 빠르게 흡수되면서 평균 ~10⁻² M☉ s⁻¹의 높은 질량 흡수율을 보인다. 이는 짧은 감마선 폭발(GRB)의 엔진 후보로서 전자기적 신호를 제공할 가능성을 시사한다.

상세 분석

본 논문은 중성자별(NS)이 주변 물질을 흡수해 질량이 최대 허용치를 초과하면서 블랙홀(BH)로 붕괴하는 ‘Accretion‑Induced Collapse(AIC)’ 현상을, 기존에 진공에서의 붕괴만을 다루던 연구와 달리 원반(torus)을 포함한 3차원 전일반상대론(General Relativistic) 시뮬레이션으로 최초 구현하였다. 초기 모델은 RNS 코드를 이용해 다중 회전 NS를 구축하고, TORERO 코드를 통해 질량이 NS 질량의 약 7 %에 해당하는 원반을 삽입하였다. 원반의 자가 중력은 무시했으며, 원반이 회전 BH의 스페이스타임에 존재한다는 가정 하에 Hamiltonian 제약 위반은 해상도 상승에 따라 최소화되었다. 붕괴는 인위적으로 압력을 0.1 % 감소시켜 유도했으며, 유체는 Γ=2 다항식 EOS를 사용한 이상유체로 모델링하였다. 자기장은 전혀 포함되지 않았으며, 이는 향후 연구 과제로 남겨졌다.

수치 해석은 Whisky 코드(플럭스 보존형 GRMHD)와 Ccatie 코드(형식 무시 트레이스프리 자유도 해석)를 결합해 수행했으며, Carpet 드라이버를 통한 고정 메쉬 재세분화(Fixed Mesh Refinement)로 7단계의 격자 정밀도를 적용하였다. 대칭 조건으로는 z=0 평면에 반사 대칭, x=0 평면에 π-대칭을 사용해 계산량을 절감했지만, 대칭을 해제한 실험에서도 축대칭성이 유지됨을 확인하였다.

시뮬레이션 결과, 원반이 없는 경우(‘vacuum collapse’)는 BH 형성 직후 모든 물질이 사건지평선 안으로 빨려 들어가며, 남은 질량이 거의 0에 수렴한다. 반면 원반이 있는 경우는 BH 형성 후에도 원반이 남아 지속적으로 흡수되며, 초기 흡수율은 약 10⁻² M☉ s⁻¹에 달한다. 이는 전형적인 NS‑NS 합병 후 형성되는 토러스와 비슷한 수준이며, 알파‑점성(α≈0.1)을 가정하면 실제 흡수율은 1 M☉ s⁻¹에 이를 수도 있다.

중력파 측면에서는 붕괴가 거의 축대칭적으로 진행되므로 l=2, m=0 모드가 지배적이며, 원반이 있든 없든 파형 차이는 미미하다. 원반이 있는 경우에만 붕괴 직전 NS 진동에 의해 작은 진동이 추가되지만, 전체적인 파형은 급격한 지수 증가와 BH 링다운(ring‑down)으로 구성된다. 현재 2세대 LIGO/Virgo는 은하계 내(≈10 kpc)에서만 검출 가능하지만, 차세대 Einstein Telescope은 ≈1 Mpc 거리까지 탐지할 수 있다.

전자기적 관점에서 원반의 지속적인 흡수는 높은 방출 효율(ε≈0.1)을 가정하면 ≈10⁵¹ erg s⁻¹ 수준의 γ‑레이 광도를 제공한다. 이는 짧은 감마선 폭발(short GRB)의 전형적인 피크 광도와 일치하므로, AIC가 짧은 GRB의 잠재적 엔진임을 뒷받침한다. 다만 자기장이 포함되지 않았기 때문에, 실제 제트 형성 메커니즘이나 플라즈마 방출 과정은 추가 연구가 필요하다.

결론적으로, 원반이 있는 AIC는 중력파 신호에서는 기존 진공 붕괴와 구별이 어려우나, 전자기적 신호—특히 짧은 GRB와 연관된 고에너지 방출—를 통해 두 현상을 구분할 수 있다. 향후 자기장, 복잡한 EOS, 더 큰 원반 질량 등을 포함한 시뮬레이션이 필요하며, 관측 측면에서는 다중 파장(멀티‑밴드) 관측 전략이 중요하다.