Hydrodynamic simulations of the recurrent nova T Coronae Borealis: Nucleosynthesis predictions

기술 요약: 재발성 노바 T Coronae Borealis의 유체역학적 시뮬레이션: 핵합성 예측

문제 및 배경
재발성 노바(Recurrent novae, RNe)는 매우 짧은 재발 주기(통상 1년에서 100년 사이)를 갖는 백색왜성(white dwarf, WD) 표면에서 발생하는 열핵 폭주(thermonuclear runaways, TNRs)를 특징으로 하는 쌍성계이다. 이러한 짧은 주기는 매우 높은 질량 축적률, 샤드라세카르 한계에 근접하는 백색왜성의 질량, 그리고 높은 초기 백색왜성 광도를 필요로 한다. T Coronae Borealis(T CrB)는 1866년과 1946년에 폭발했으며, 약 80년의 재발 주기가 예측되는 저명한 은하계 내 재발성 노바이다. 다음 폭발이 임박함에 따라, 시스템을 구동하는 조건을 이해하기 위해 물리적 매개변수(구체적으로 백색왜성의 질량, 초기 광도, 축적률 및 금속 함량)를 제약해야 할 결정적인 필요성이 존재한다. 또한, 재발성 노바는 제Ia형 초신성의 전구체일 가능성이 있으므로, 이들의 핵합성과 질량 축적 과정을 규명하는 것은 매우 중요하다.

방법론
저자들은 T CrB의 폭발을 시뮬레이션하기 위해 유체역학적 라그랑주 유한차분 코드인 SHIVA를 사용하였다. 시뮬레이션에는 퇴화 전자 가스, 다성분 이온 플라즈마, 복사를 포함하는 일반 상태 방정식과 OPAL 로슬랜드 평균 불투명도가 활용되었다. 핵 에너지 생성은 STARLIB 데이터베이스의 최신 비율을 반영하여 120개 종($^1$H부터 $^{48}$Ti까지)과 630개 공정으로 연결된 반응 네트워크를 통해 계산되었다.

본 연구는 다양한 조건 하에서 80년의 재발 주기($\tau_{rec}$)를 재현하도록 설계된 11개의 새로운 유체역학 모델을 제시한다:

• 백색왜성 질량 ($M_{WD}$): 1.30, 1.35, 1.38 $M_{\odot}$.
• 초기 백색왜성 광도 ($L_{WD}$): 0.01, 0.1, 1.0 $L_{\odot}$.
• 축적된 금속 함량 ($Z_{acc}$): 태양 질량 ($1.0 Z_{\odot}$), 태양 미달 ($0.1 Z_{\odot}$), 태양 초과 ($10 Z_{\odot}$).
• 혼합(Mixing): 축적된 외포층과 하부 백색왜성 층(탄소-산소 및 산소-네온 조성) 사이의 혼합(10% 및 25%) 영향을 평가하기 위해 1.35 및 1.38 $M_{\odot}$ 백색왜성에 대해 8개의 추가 모델을 계산하였다.
• 저질량 백색왜성: 더 낮은 질량 추정치를 가진 재발성 노바를 고려하기 위해 1.20 $M_{\odot}$ 백색왜성에 대한 3개의 모델을 계산하였다.

각 구성에 대하여, 정확히 80년의 축적 후 폭발을 일으키도록 질량 축적률($\dot{M}_{acc}$)을 조정하였다.

주요 결과

• 축적률 및 물리적 조건: 태양 조성 물질로 80년의 재발 주기를 달성하기 위해서는, 1.30 - 1.38 $M_{\odot}$ 질량 및 0.01 - 1.0 $L_{\odot}$ 광도를 가진 백색왜성에 대해 $\dot{M}_{acc} \sim 10^{-8} - 10^{-7} M_{\odot} \text{yr}^{-1}$의 질량 축적률이 필요하다. 낮은 초기 백색왜스 광도, 낮은 백색왜스 질량, 또는 축적 물질의 낮은 금속 함량은 80년이라는 시간 범위 내에서 임계 점화 압력에 도달하기 위해 더 높은 축적률을 필요로 한다.
• 폭발 역학: 낮은 $L_{WD}$ 또는 낮은 $Z_{acc}$를 가진 모델은 더 많은 축적된 외포층 질량과 더 높은 점화 압력을 초래한다. 결과적으로, 이 모델들은 더 높은 최고 온도($T_{max}$), 더 큰 방출 질량($\Delta M_{eje}$), 그리고 더 큰 운동 에너지를 특징으로 하는 더 격렬한 폭발을 생성한다. 예를 들어, $L_{WD}=0.01 L_{\odot}$이고 태양 금속 함량을 가진 1.38 $M_{\odot}$ 백색왜성은 $T_{max} \approx 2.92 \times 10^8$ K에 도달하는 반면, $L_{WD}=1.0 L_{\odot}$인 1.30 $M_{\odot}$ 백색왜성은 $2.01 \times 10^8$ K에만 도달한다.
• 핵합성 및 원소 함량:
  • 동위원소 비율: 방출물은 CNO 순환 동위원소($^{12,13}$C, $^{14,15}$N)에 의해 지배된다. 서로 다른 모델 간에 동위원소 비율(예: $^{12}$C/$^{13}$C, $^{16}$O/$^{17}$O, $^{28}$Si/$^{29}$Si)의 상당한 변화가 관찰된다.
  • 백색왜스 질량 의존성: 동일한 초기 광도를 가졌으나 백색왜스 질량이 다른 모델들은 Ne, Na, Mg, Al, Si, P, S, Ar, K, Ca, Sc의 함량에서 상당한 차이를 보인다. 구체적으로, 1.38 $M_{\odot}$ 모델은 태양 함량에 비해 S, Cl, K, Ca, Sc에서 상당한 농축을 보이는 반면, 1.20 $M_{\odot}$ 및 1.30 $M_{\odot}$ 모델은 (모델 130C를 제외하고) 일반적으로 이 원소들에 대해 태양 수준 또는 약간 낮은 수준의 함량을 보인다. 실리콘(Si)은 1.20 및 1.30 $M_{\odot}$ 모델에서 농축되지만, 1.38 $M_{\odot}$ 모델에서는 결핍되거나 태양 수준에 가깝다.
  • 금속 함량 효과: 1.38 $M_{\odot}$ 모델에서 금속 함량의 변화는 조성을 급격히 변화시킨다. 낮은 금속 함량($0.1 Z_{\odot}$)은 상당한 $^{40}$Ca 생성을 유도하는 반면, 높은 금속 함량($10 Z_{\odot}$)은 N, C, Ne, Si로의 우세를 이동시킨다.
• 혼합 시나리오: 외포층과 백색왜스 기질(CO 또는 ONe) 사이의 혼합을 포함하는 모델은 비혼합 모델에 비해 축적 단계가 단축되고 총 방출 질량이 감소함을 보여준다. 결정적으로, 혼합은 감마선 방출체($^{7}$Be, $^{18}$F, $^{22}$Na, $^{26}$Al)의 생성을 크게 향상시킨다. 비혼합 모델은 무시할 만한 수준의 감마선 방출을 예측하는 반면, 혼합 모델은 검출 가능한 신호를 생성할 수 있어, CO와 ONe 백색왜스를 구분할 수 있는 잠재적 지표가 될 수 있다.
• 감마선 방출: 표준 비혼합 모델의 경우, 방사성 동위원소의 낮은 수율로 인해 감마선 선(511 keV 및 1275 keV)과 그와 관련된 연속 스펙트럼은 현재의 관측 기기들의 검출 한계보다 수십 배 이상 낮을 것으로 저자들은 예측한다.

의의 및 주장
본 논문은 방출물의 조성 차이, 특히 Ne, Na, Mg, Al, Si, P, S, Ar, K, Ca, Sc와 같은 종의 원소 함량이 잠재적인 진단 도구 역할을 한다고 주장한다. 이러한 변화는 매개변수 공간을 제약하고 본 연구에서 제안된 다양한 T CrB 모델들을 구별할 수 있다. 구체적으로, 저자들은 다가올 T CrB 폭발에 대한 초기 분광 측정값이 백색왜스의 질량과 혼합 정도를 결정함으로써, 시스템의 폭발을 지배하는 물리적 조건을 명확히 할 수 있다고 주장한다. 본 연구는 감마선 검출이 혼합 없이 이루어지기는 어렵지만, 원소 함량 패턴이 밑바닥에 있는 백색왜스의 본질과 작용하는 혼합 과정을 식별하는 강력한 방법이 될 것임을 결론짓는다.