펄스 역방향 폭발 모델 디터전션 점화 조건 연구

초록

이 논문은 백색왜성의 비효율적 디플라그레이션이 펄스를 일으키고, 그 후 형성되는 물질 흡수 충격파가 중심 코어를 압축·가열하여 역방향 디터전션을 유도하는 ‘펄스 역방향 폭발(Pulsating Reverse Detonation, PRD)’ 메커니즘을 3차원 수치 시뮬레이션으로 검증한다. 디플라그레이션 동안 소모된 질량 M_defl 을 변수로 삼아 다양한 경우를 탐색했으며, M_defl 이 백색왜성 결합에너지에 근접(≈0.30 M☉)할 때는 충격파가 충분히 핵연료를 가열·구속하지 못해 디터전션이 불안정해진다는 결론을 얻었다.

상세 분석

본 연구는 PRD 시나리오의 핵심인 ‘역방향 디터전션’이 실제로 발생할 수 있는 물리적 조건을 정량적으로 규명하고자 3차원 스무딩 입자 수소역학(SPH) 코드를 사용하였다. 초기 조건은 Chandrasekhar 질량(≈1.38 M☉)의 탄소‑산소 백색왜성을 가정하고, 디플라그레이션 단계에서는 핵반응을 비활성화한 채로 유동역학만을 추적하였다. 이는 디플라그레이션이 방출한 에너지(≈E_defl)와 그에 따른 질량 손실 M_defl 을 정확히 제어하기 위함이다. 여러 M_defl 값(0.15–0.35 M☉)에 대해 시뮬레이션을 수행했으며, 각 경우에 대해 충격파가 형성되는 시점, 충격파 전파 속도, 핵심 코어의 압축 비율, 온도 상승을 측정하였다. 핵심 결과는 다음과 같다. 첫째, M_defl 이 작을수록(≈0.15 M☉) 디플라그레이션이 남긴 잔류 에너지가 충분히 작아 백색왜성 외피가 강하게 수축하고, 이때 형성되는 물질 흡수 충격파는 코어 표면에 도달해 압축·가열을 유도한다. 충격파가 코어 내부에 도달했을 때 온도는 ≈2×10⁹ K, 밀도는 ≈2×10⁷ g cm⁻³에 이르러 디터전션 임계조건을 만족한다. 둘째, M_defl 이 결합에너지와 거의 동등한 수준(≈0.30 M☉)에 이를 경우, 디플라그레이션이 방출한 에너지로 인해 외피가 과도하게 팽창하고, 재수축 과정에서 형성되는 충격파는 약해진다. 이때 충격파가 코어를 충분히 가열하지 못해 온도가 ≈1×10⁹ K 수준에 머물며, 디터전션이 자발적으로 발생하지 않는다. 셋째, M_defl 이 0.35 M☉를 초과하면 백색왜성 자체가 해체될 위험이 있어 PRD 메커니즘이 전제하는 ‘핵심 코어 보존’ 조건이 깨진다. 이러한 결과는 PRD가 ‘중간 정도’의 디플라그레이션(0.15–0.25 M☉)에서 가장 안정적으로 작동함을 시사한다. 또한, 충격파의 형성 및 전파 과정에서 비선형 유동, 회전, 그리고 작은 규모의 난류가 디터전션 임계조건에 미치는 영향을 정량화했으며, 이는 향후 고해상도 시뮬레이션에서 추가 검증이 필요함을 강조한다.