장시간 감마선 폭발: 붕괴성 원반의 유체 불안정성

초록

이 논문은 실제 별 붕괴 모델을 시작점으로 3차원 SPH 시뮬레이션을 수행해, 장시간 감마선 폭발(LGRB)의 핵심인 블랙홀·원반 시스템을 재현한다. 고온·고밀도 원반은 뉴트리노 냉각을 받으며 m=1,2 비축축 모드와 클럼핑 현상을 보이고, 이러한 나선 구조가 동역학적 시간척도 내에 각운동량을 외부로 효율적으로 전달한다. 결과적으로 0.1–1 M☉ s⁻¹ 수준의 높은 질량이입률을 유지해 제트 형성에 필요한 에너지를 공급한다.

상세 분석

본 연구는 기존 2차원 축대칭 모델이 놓치기 쉬운 비축축 불안정성을 3차원 Smoothed Particle Hydrodynamics(SPH) 코드인 Gadget‑2를 기반으로 정밀하게 탐구한다. 초기 조건은 최신 핵융합 전 단계 별 모델을 사용해, 핵심 붕괴 직후 형성되는 블랙홀(≈3 M☉)과 그 주변에 형성되는 원반을 자연스럽게 재현한다. 코드 수정 사항으로는 (1) 뉴트리노 쿨링을 구현하기 위해 적절한 방출률(URCA, 전자‑양성자 캡처 등)을 적용하고, (2) 블랙홀의 강한 중력장을 근사하기 위해 Paczyński‑Wiita 의사‑포텐셜을 도입해 일반상대론적 효과를 반영하였다. 이러한 물리적 처리 덕분에 원반 내부 온도는 10¹⁰ K, 밀도는 10¹⁰–10¹² g cm⁻³에 이르는 초고온·초고밀도 상태를 유지한다.

시뮬레이션 결과, 원반은 즉시 수소와 헬륨의 핵융합이 억제된 상태에서 뉴트리노 냉각에 의해 압력이 감소하고, 그에 따라 토러스 형태의 고밀도 영역이 형성된다. 이때 원반은 토러스형 압력 구배와 회전 속도 구배가 강해지는 영역에서 Rayleigh‑Taylor 및 Kelvin‑Helmholtz와 유사한 유체 불안정성을 보이며, 특히 m=1(일차 나선)과 m=2(이중 나선) 모드가 급격히 성장한다. 이러한 비축축 모드는 원반 전체에 걸쳐 1–2 회전 주기(≈10⁻³ s) 이내에 비선형 단계에 도달하고, 고밀도 클럼프(덩어리)와 강한 나선 파동을 만든다.

핵심적인 물리적 효과는 각운동량의 외부 전달이다. 나선 구조는 전파되는 토르크를 통해 내부 물질의 각운동량을 외부로 빼앗아, 내부 물질이 급격히 붕괴해 블랙홀로 흡수되게 만든다. 측정된 질량이입률은 0.1–1 M☉ s⁻¹ 수준으로, LGRB 모델이 요구하는 최소값을 충분히 초과한다. 또한, 클럼프가 형성된 지역에서는 뉴트리노 방출이 집중되어 냉각 효율이 증가하고, 이는 다시 불안정성 성장에 피드백을 제공한다.

코드 검증을 위해 동일 초기 조건에서 2차원 축대칭 시뮬레이션과 비교했을 때, 비축축 모드가 없을 경우 질량이입률이 30 % 이하로 감소함을 확인했다. 이는 기존 2D 모델이 과소평가한 핵심 메커니즘을 강조한다. 마지막으로, 시뮬레이션은 블랙홀 주변에 형성된 저밀도 채널을 통해 고에너지 전자·양성자 쌍이 빠져나가며, 이들이 전자기적 과정(예: 블렌더-코어 전도)과 결합해 초고속 제트를 촉발할 가능성을 시사한다.