붕괴성 별 낙하와 감마선 폭발 X선 빛 곡선

초록

축대칭 2차원 FLASH 시뮬레이션을 이용해 회전하는 장거리 감마선 폭발(롱 GRB) 전구 별인 콜랩서의 장기 낙하 과정을 연구하였다. 시뮬레이션은 점성 효과를 포함하고, 중심 흡수율이 초기 수십 초 동안 급격히 감소한 뒤 t⁻²·⁸ 형태의 전력법칙을 보이며 X선 빛 곡선의 급락 단계와 일치한다. 이후 저각운동량 물질이 축 방향으로 역류해 얇은 원반을 형성하고, 흡수율이 약 5×10⁻⁵ M☉ s⁻¹ 수준의 ‘플래토’ 단계로 전환된다. 저온(≈10¹⁰ K)에서 중성미자 냉각이 사라지는 것이 급락을 촉발한다는 분석 모델도 제시한다.

상세 분석

본 논문은 콜랩서 모델을 기반으로, 회전하는 대질량 별이 중심 블랙홀(또는 중성자별)로 낙하하면서 발생하는 복잡한 유체역학을 고해상도 적응격자(AMR) 코드 FLASH로 2차원 축대칭 시뮬레이션하였다. 점성 항을 명시적으로 포함함으로써 원반 내부의 각운동량 전달과 열 발생을 현실적으로 재현했으며, 이는 기존의 무점성 시뮬레이션보다 물리적 타당성을 크게 높인다. 초기 단계에서는 별의 외피가 자유낙하하면서 높은 각운동량을 보유한 물질이 원반 형태로 원형화(circularization)되고, 이때 발생하는 충격파가 별 내부를 관통한다. 충격파가 외부로 전파되면, 충격 가열된 물질은 압력지지된 대기와 얇은 원반을 형성한다. 특히, 충격 전후의 온도 변화를 추적한 결과, 온도가 ≈10¹⁰ K 이하로 떨어지는 순간 중성미자 냉각 효율이 급격히 감소한다는 점을 확인하였다. 이는 에너지 손실이 급감하면서 내부 압력이 상승하고, 충격파가 다시 팽창하게 되는 메커니즘을 제공한다.

시뮬레이션에서 관측된 중심 흡수율(Ṁ)의 시간 의존성은 세 단계로 구분된다. 첫 번째 ‘프롬프트’ 단계는 수십 초 동안 Ṁ이 급격히 감소하며, 이때의 감소율은 t⁻²·⁸에 근접한다. 이는 관측된 X선 빛 곡선의 L_X∝t⁻³ 급락과 거의 일치한다. 두 번째 단계에서는 충격파가 별 외부까지 도달한 뒤, 저각운동량 물질이 축 방향 퍼널을 통해 역류하면서 원반 표면에서 바람(아웃플로우)이 발생한다. 이때 Ṁ는 급격히 완만해져 약 5×10⁻⁵ M☉ s⁻¹ 수준의 ‘플래토’에 머무른다. 세 번째 단계는 시뮬레이션 시간 한계에 도달하기 전까지 유지되며, 실제 GRB에서는 수천 초에 걸친 장기 플래토와 일치한다는 점에서 의미가 크다.

또한, 저자들은 중성미자 냉각이 사라지는 임계 온도(≈10¹⁰ K)를 기준으로 한 분석 모델을 제시하였다. 이 모델에 따르면, 원반이 형성되는 반경이 충분히 커지면(≈10⁹ cm 이상) 물질의 비열이 급격히 감소하고, 따라서 충격파가 외부로 팽창하게 된다. 이 과정이 X선 급락의 물리적 트리거로 작용한다는 가설은, 기존에 제시된 ‘고에너지 전자기 방출’ 모델보다 중심 엔진의 변동성을 직접 연결시켜 설명한다는 점에서 혁신적이다.

마지막으로, 시뮬레이션 해상도에 따라 ‘프롬프트’ 단계의 지속 시간이 변한다는 점을 지적하면서, 실제 물리적 해상도(≈10⁶ cm 이하)에서는 약 20 s 정도가 될 것으로 예측한다. 이는 관측된 장거리 GRB의 프롬프트 감마선 지속시간과도 일치한다. 전반적으로, 본 연구는 콜랩서 낙하 과정에서 발생하는 다중 단계 흡수율 변동을 정량적으로 재현하고, 이를 X선 빛 곡선의 주요 특징과 연결함으로써 중앙 엔진 모델의 타당성을 크게 강화한다.