Black Hole Spin-down in Collapsars in 3D Neutrino Transport GRMHD Simulations

1. 연구 배경 및 문제 제기 (Problem)

• 배경: 장기간 감마선 폭발 (LGRB) 은 거대 항성의 핵이 블랙홀 (BH) 로 붕괴하는 '콜랩사 (Collapsar)' 모델에 의해 구동되는 것으로 알려져 있습니다. 이 과정에서 블랙홀의 스핀 ($a$) 과 자기장 세기가 제트의 에너지를 결정하는 핵심 요소입니다.
• 기존 연구의 한계: 이전의 비방사성 (non-radiative) 단열 두꺼운 원반 (MAD, Magnetically Arrested Disk) 시뮬레이션 연구들은 블랙홀이 매우 낮은 스핀 ($a_{eq} \approx 0.035 \sim 0.07$) 으로 빠르게 감속된다고 결론지었습니다. 이러한 낮은 스핀은 제트의 출력이 매우 약해져 항성 외피를 뚫고 나가는 데 실패할 수 있음을 시사했습니다.
• 문제점: 실제 콜랩사 원반은 고밀도 및 고온 환경으로 인해 중성미자 (neutrino) 방출을 통해 냉각됩니다. 중성미자 냉각은 원반의 두께를 줄이고 각운동량 공급을 증가시키며, 이는 블랙홀 스핀 진화에 중대한 영향을 미칠 수 있습니다. 그러나 기존 연구에서는 중성미자 수송을 고려하지 않아 실제 물리 과정을 정확히 반영하지 못했습니다.
• 목표: 중성미자 냉각이 포함된 3 차원 일반상대론적 자기유체역학 (GRMHD) 시뮬레이션을 통해, 중성미자 냉각된 원반에서 블랙홀의 스핀 진화와 최종 평형 스핀 ($a_{eq}$) 을 규명하는 것.

2. 방법론 (Methodology)

• 시뮬레이션 코드: GPU 가속화된 중성미자 -GRMHD 코드인 $\nu$H-AMR을 사용했습니다.
• 물리 모델:
  • 중성미자 수송: 에너지 적분형 ('gray') M1 근사법을 사용하여 전자 중성미자 ($\nu_e$), 반전자 중성미자 ($\bar{\nu}_e$), 그리고 무거운 렙톤 중성미자 ($\nu_x$) 의 3 종을 유체처럼 처리했습니다.
  • 상태 방정식: 헬름홀츠 (Helmholtz) 상태 방정식을 사용했습니다.
  • 초기 조건: 두 가지 다른 항성 밀도 프로파일을 가진 콜랩사 전구체 모델을 설정했습니다.
    • p# 모델: 전형적인 콜랩사 프로파일 ($\rho \propto r^{-1.5}$), 질량 흡수율 $\dot{m} \sim 0.1-1 M_\odot/s$.
    • c# 모델: 중심부에 일정한 밀도 코어를 가진 모델 ($\alpha_p=0$), 초기 높은 질량 흡수율 ($\dot{m} \gtrsim 10 M_\odot/s$).
  • 스핀 범위: 다양한 초기 블랙홀 스핀 ($a = 0.0 \sim 0.8$) 을 가진 시뮬레이션을 수행했습니다. 또한 중성미자 수송을 끄는 대조군 (c1N0, p1N0) 시뮬레이션도 수행하여 중성미자의 영향을 분리 분석했습니다.
• 분석: 블랙홀 사건의 지평선에서 질량 흡수율, 자기 플럭스, 에너지 및 각운동량 플럭스를 측정하여 스핀 업 파라미터 ($s$) 와 평형 스핀을 계산했습니다.

3. 주요 기여 및 결과 (Key Contributions & Results)

A. 중성미자 냉각의 영향 및 새로운 평형 스핀 발견

• 중성미자 냉각의 효과: 중성미자 냉각은 원반을 더 얇게 만들고 ($h/r$ 감소), 각운동량 수송 효율을 높입니다.
• 새로운 평형 스핀 ($a_{eq}$): 중성미자 냉각이 포함된 MAD 상태에서는 블랙홀이 $a_{eq} \approx 0.13$의 평형 스핀에 도달하는 것으로 확인되었습니다.
  • 이는 비방사성 (비냉각) MAD 모델의 $0.035 \sim 0.07$보다 2~4 배 높습니다.
  • 방사성 냉각이 포함된 얇은 원반 모델의 예측치 ($a_{eq} \approx 0.3$) 와 비방사성 모델 사이에서 중간 값을 가지며, 원반 두께가 얇아질수록 평형 스핀이 증가한다는 기존 경향을 지지합니다.

B. 제트 출력 및 LGRB 성공 가능성

• 제트 파워 증가: $a_{eq} \approx 0.13$은 비냉각 모델보다 4~16 배 더 강력한 LGRB 제트를 생성합니다.
• 항성 탈출: 이 정도의 출력은 제트가 항성 외피를 뚫고 나오는 (breakout) 데 훨씬 유리하며, 관측된 LGRB 의 에너지 분포와 더 잘 부합합니다.
• 스핀 감속 시간: 일반적인 LGRB 지속 시간 ($t \gtrsim 30s$) 동안 블랙홀은 충분한 질량을 흡수하여 $a_{eq} \approx 0.13$에 도달할 수 있습니다. 다만, 더 짧은 지속 시간이나 낮은 흡수율의 경우 더 빠르게 회전하는 블랙홀이 남을 수 있습니다.

C. 중성미자의 직접적 vs 간접적 역할

• 직접적 토크: 블랙홀 스핀 진화에 대한 중성미자 플럭스의 직접적인 토크 기여도는 무시할 수 있을 정도로 작았습니다.
• 간접적 영향: 중성미자는 원반의 열역학적 구조 (두께) 를 변화시켜 전자기적 (EM) 및 유체역학적 (HD) 토크의 세기를 조절함으로써 스핀 진화에 결정적인 영향을 미쳤습니다.

D. 불안정성과 변동성

• 저스핀 모델의 불안정성: 초기 스핀이 낮은 ($a \lesssim 0.2$) 모델, 특히 얕은 밀도 프로파일을 가진 p# 모델에서는 제트가 **글로벌 킥 불안정성 (kink instability)**을 겪으며 붕괴하거나 뒤틀리는 현상이 관찰되었습니다.
• 결과: 이로 인해 원반과 블랙홀 스핀 축이 정렬되지 않고 (misalignment), MAD 상태가 일시적으로 깨지거나 제트 효율이 변동하는 현상이 발생했습니다. 이는 저에너지 또는 실패한 제트 (subluminous/failed GRBs) 를 설명할 수 있는 메커니즘을 제공합니다.

4. 의의 및 결론 (Significance)

1. LGRB 엔진의 물리적 제약: 이 연구는 LGRB 를 구동하는 블랙홀 엔진이 초기에 매우 빠르게 회전하더라도, 중성미자 냉각된 MAD 상태에서는 빠르게 감속되어 $a \approx 0.13$ 근처의 낮은 스핀으로 수렴함을 보였습니다.
2. 중력파 관측과의 일치: LIGO/Virgo/KAGRA 가 관측한 블랙홀 쌍성계 (BBH) 병합의 스핀 분포는 대체로 낮은 값 ($a \le 0.3$) 에 집중되어 있습니다. 본 연구에서 도출된 $a_{eq} \approx 0.13$은 이러한 중력파 관측 데이터와 매우 잘 일치하며, 콜랩사가 BBH 의 전구체일 가능성을 강력히 지지합니다.
3. 관측적 함의:
   • 최대 광도 제한: 블랙홀 스핀이 감속되는 과정이 제트 파워를 제한하여, 관측된 LGRB 의 최대 광도 상한선을 설명할 수 있습니다.
   • 변동성: 저스핀 상태에서의 제트 불안정성은 LGRB 의 광도 변동성이나 실패한 폭발 (failed jets) 을 설명하는 데 기여할 수 있습니다.
4. 모델링의 발전: 중성미자 수송을 포함한 3D GRMHD 시뮬레이션을 통해 콜랩사 물리 모델의 정합성을 크게 높였으며, 향후 블랙홀 스핀 진화와 제트 역학 연구의 새로운 기준을 제시했습니다.

요약: 본 논문은 중성미자 냉각이 포함된 정교한 시뮬레이션을 통해, 콜랩사 블랙홀이 비냉각 모델보다 높은 평형 스핀 ($a \approx 0.13$) 에 도달함을 증명했습니다. 이는 LGRB 가 항성 외피를 뚫고 나올 수 있을 만큼 강력한 에너지를 가지면서도, 중력파 관측에서 보이는 낮은 스핀 분포와 모순되지 않는 물리적 메커니즘을 제공합니다.