Multi-Epoch X-Ray Detection of SLSN-I 2018bsz: Constraints on the Powering Mechanism and Ejecta Structure

기술 요약: SLSN-I 2018bsz의 다중 에포크 X선 검출

문제 및 동기
초광도 초신성(SLSNe)은 일반적인 코어 붕괴 초신성보다 약 100배 더 큰 에너지를 가진 희귀한 항성 폭발이다. 이러한 사건을 구동하는 메커니즘은 여전히 논쟁 중이며, 주요 후보로는 이젝타(ejecta)와 주변 성간 물질(CSM) 사이의 강력한 상호작용, 또는 빠르게 회전하는 고자기장 중성자별(마그네타, magnetar)로부터의 에너지 주입이 있다. 광학 모델링은 많은 SLSN-I의 광도 곡선을 마그네타 모델로 성공적으로 재현해 왔으나, 서로 다른 메커니즘이 유사한 광학적 징후를 생성할 수 있기 때문에 동력을 구분하는 것은 어렵다. X선 관측은 이러한 퇴행성을 깨뜨릴 수 있는 결정적인 탐사 도구를 제공한다. 그러나 이젝타 흡수가 감소하고 마그네타 방출이 가시화될 수 있는 후기 에포크에서 X선 검출은 매우 드물다. 본 연구 이전에는 단 하나의 SLSN-I(SCP 06F6)만이 X선이 확정적으로 검출된 바 있다.

방법론
저자들은 가장 가까운 것으로 알려진 스트립드(stripped) SLSN-I(거리 ~111 Mpc)인 SN 2018bsz에 대한 다중 에포크 X선 분광 연구를 제시한다. 데이터셋은 폭발 후 87일에서 304일 사이에 촬영된 4개의 Chandra 관측(ACIS-S/S3)과 폭발 후 1253일째에 촬영된 1개의 XMM-Newton 관측(EPIC-pn/MOS)으로 구성된다.

데이터 감소(Data reduction)는 Chandra의 경우 CIAO(v4.13), XMM의 경우 SAS(v22.1.0)를 사용하여 표준 절차를 따랐다. 낮은 카운트율로 인해 분광 분석에는 C-통계량을 사용하였으며 XSPEC(v12.15.0c)을 활용하였다.

• 소스 식별: Chandra 이미지는 SN 소스와 XMM 점확산 함수(PSF) 내에 있는 세 개의 더 희미한 인근 X선 소스(S1–S3)를 밝혔다.
• 분광 모델링: 스펙트럼은 은하 흡수($N_H = 2.02 \times 10^{21}$ cm$^{-2}$)로 피팅되었다. 모델에는 파워 로(power law, 마그네타 스핀다운 또는 비열적 상호작용을 나타냄)와 열적 플라즈마 모델(mekal, CSM 충격 상호작용을 나타냄)이 포함되었다.
• 오염 처리: 후기 XMM 관측의 경우, 저자들은 Chandra 데이터로부터 유도된 "오염 환형 영역(contamination annulus)"을 사용하여 인근 소스(S1–S3)의 기여도를 모델링하였다. XMM 스펙트럼은 tbabs(powCont.A + pow) 또는 tbabs(powCont.A + mekal) 형태의 합성 모델로 피팅되었으며, 이때 오염 성분은 Chandra 결과에 따라 고정되었다.
• 동력원 제약 조건: 저자들은 두 가지 시나리오를 테스트하였다:
  1. 마그네타: 쌍극자 스핀다운 모델을 가정하여 초기 회전 주기($P_0$)와 자기장($B$)을 제약하기 위해 마르코프 연쇄 몬테카를로(MCMC) 피팅을 수행하였다. 이 결과들을 Gomez et al. (2024)의 광학 광도 곡선 피팅 결과와 비교하였으며, 이젝타 흡수 및 이온화 돌파(ionization breakout)의 영향을 고려하였다.
  2. CSM 상호작용: 질량 손실률($\dot{M}$)은 Chevalier & Fransson (2017)의 광도-질량 손실 관계를 사용하여 추정되었다.

주요 결과

1. 검출: SN 2018bsz는 4개의 모든 Chandra 에포크에서 >3$\sigma$ 유의성으로 검출되었다. 1253일째의 XMM 관측에서도 잠정적인 검출이 이루어졌으나, 이는 인근 소스의 오염 수준에 민감하다. 이로써 SN 2018bsz는 두 번째로 X선이 검출된 SLSN-I이자, 세 번째로 X선이 검출된 SLSN 전체가 되었다.
2. 분광 진화:
   • 초기 에포크 (Chandra): 스펙트럼은 광자 지수 $\Gamma \approx 0.9–1.9$인 단단한 파워 로(hard power laws)로 잘 설명된다 (동시 피팅 시 최적값 $\Gamma = 1.32$).
   • 후기 에포크 (XMM): 스펙트럼이 크게 부드러워진다(soften). 1253일째의 데이터는 부드러운 파워 로($\Gamma < 4.0$) 또는 $kT \approx 0.20$ keV인 열적 플라즈마 모델과 일치한다.
3. 광도 진화: X선 광도는 시간 경과에 따라 상대적으로 평탄하게 유지되며 ($L_{0.5-8 \text{ keV}} \approx 1.1 \times 10^{40}$ erg s$^{-1}$), 87일에서 304일 사이에 거의 감소하지 않았다.
4. 마그네타 제약 조건:
   • (이젝타 흡수를 무시한다고 가정한) X선 데이터에 대한 MCMC 피팅은 약한 자기장 ($\log B \approx 13.1$)과 긴 초기 회전 주기 ($P_0 \approx 16$ ms)를 산출하였다.
   • 이 매개변수들은 Gomez et al. (2024)이 광학 광도 곡선으로부터 도출한 "최적합" 마그네타 매개변수 ($P_0 \approx 0.9$ ms, $\log B \approx 14.4$)와 일치하지 않는다.
   • 광학 유도 매개변수($B_{opt}, P_{opt}$)를 사용할 경우, 마그네타 모델은 가장 초기 단계를 제외한 모든 에포크에서 관측된 X선 광도를 과소 예측하며 평탄한 광도 곡선을 재현하는 데 실패한다.
   • 이젝타 흡수를 고려하면 불일치는 더욱 심화된다. 이온화 돌파 모델은 유리한 매개변수를 고려하더라도 마그네타 방출이 폭발 후 ~2–72년이 지나야 이젝타를 탈출할 수 있음을 시사하며, 이는 관측 창을 훨씬 넘어서는 시점이다.
5. CSM 상호작용 제약 조건:
   • 스펙트럼의 연화(softening)와 평탄한 광도 진화는 시간이 지남에 따라 약화되는 초기 강한 이젝타-CSM 상호작용과 일치한다.
   • 유도된 질량 손실률은 $\sim 5.4 \times 10^{-3}$에서 $1.4 \times 10^{-2}$ $M_{\odot}$ yr$^{-1}$이다. 이 비율은 일반적인 스트립드 초신성 및 대부분의 SLSN-I를 초과하며, SLSN-II 또는 분출적 질량 손실을 통해 상호작용하는 시스템으로 진화하는 스트립드 초신성과 더 가깝다.
   • 추정된 전구체 질량 손실 에피소드의 시기는 폭발 전 5–18년이며, 이는 이전의 광학 연구에서 제안된 껍질형(shell-like) 또는 원반형(disk-like) CSM 구조와 일치한다.

의의 및 결론
본 논문은 SN 2018bsz의 X선 방출이 표준적인 밀리초 마그네타 중앙 엔진에 의해 구동될 가능성이 낮다고 결론짓는다. 왜냐하면 X선 데이터를 맞추기 위해 요구되는 매개변수가 광학 모델링에서 도출된 것과 모순되며, 예측된 이온화 돌파 시간도 너무 길기 때문이다.

대신, 관측 결과는 X선 방출이 밀도가 높고 비대칭적인(원반 형태의) CSM과 이젝타 사이의 상호작용에 의해 지배된다는 시나리오를 지지한다. 이 해석은 다음 사항들에 의해 뒷받받된다:

• X선 광도의 평탄한 시간적 진화.
• 단단한 스펙트럼에서 부드러운 스펙트럼으로의 진화.
• 분출적 전구체 활동과 일치하는 높은 유도 질량 손실률.
• 수소 방출 및 원반형 CSM과의 상호작용을 나타내는 이전의 광학 결과.

저자들은 SN 2018bsz가 CSM 상호작용이 강력한 방출에 핵심적인 역할을 하는, 그리고 만약 마그네타가 존재하더라도 관측된 고에너지 출력을 지배하지 않는 독특한 SLSN-I 하위 부류를 나타낸다고 주장한다. 이 연구는 극단적인 초신성을 구분하기 위해 다파장 관측이 필수적임을 강조한다.