ENGRAVE follow-up of a type IIb supernova spatially coincident with the sub-threshold gravitational wave trigger S250818k

기술 요약: 하위 역치 중력파 유발 S250818k 와 공간적으로 일치하는 IIb 형 초신성에 대한 ENGRAVE 후속 연구

문제 제기
LIGO-Virgo-KAGRA(LVK) 네트워크의 네 번째 관측 기간 (O4) 동안, 이항성 중성자별 (BNS) 병합 분류 확률이 29% 이고 국소화 영역이 매우 넓음 (90% 신뢰도 기준 949 제곱도) 으로 확인된 저신호 중력파 (GW) 후보 S250818k 가 식별되었습니다. Zwicky 임시 천체 시설 (ZTF) 은 이 영역 내에서 유발 후 약 2.7 시간 뒤에 광학 임시 천체 SN 2025ulz 를 발견했습니다. 이 임시 천체는 초기 급격한 감쇠와 특징이 없는 스펙트럼을 보였는데, 이는 초기에 중성자별 질량이 태양 질량 미만인 경우를 포함하는 킬로노바 (KN) 나 "슈퍼킬로노바" 시나리오와 유사한 특성이었습니다. 주요 과학적 과제는 GW 사건의 진정한 전자기 대응체인지 우연의 일치인지, 특히 KN 신호를 모방할 수 있는 젊은 초신성 군집과 구별하여 SN 2025ulz 를 명확하게 분류하는 것이었습니다.

방법론
ENGRAVE(초대형 망원경에서의 중력파 원천의 전자기 대응체) 협력단은 ESO 초대형 망원경 (VLT) 및 기타 시설을 광범위한 다중 파장 관측 캠페인에 활용했습니다.

• 영상 촬영: VLT/FORS2(R, I, z 대역), VLT/HAWK-I(K 대역), NOT/ALFOSC(r 대역) 를 사용하여 광학 및 근적외선 (NIR) 영상 촬영을 수행했습니다. 이는 은하계 배경 차단을 위한 심층 참조 템플릿 구축을 위해 허블 우주 망원경 (HST) 관측 (F336W, F606W, F110W, F160W) 으로 보완되었습니다. 차등 영상 촬영은 ECSNOoPY 및 HOTPANTS 파이프라인을 사용하여 임시 천체 플럭스를 분리하도록 실행되었습니다.
• 분광학: VLT/X-shooter(UVB, VIS, NIR), VLT/FORS2(300I 격자), VLT/MUSE(적분 필드 단위) 를 사용하여 분광 후속 관측을 획득했습니다. MUSE 데이터는 임시 천체를 은하계 배경으로부터 공간적으로 분리하는 데 결정적이었으며, 다중 구멍 및 고리 방법을 사용하여 은하계 배경을 모델링하고 차감함으로써 깨끗한 스펙트럼을 추출할 수 있었습니다.
• 전파 관측: VLA, MeerKAT, uGMRT, e-MERLIN 을 사용하여 전파 파장에서 이 천체를 모니터링하여 잔광 신호를 탐색하거나 제트 경사도에 대한 제약을 확인했습니다.
• 분석: 광곡선을 구성하고 색상 변화를 분석했습니다. SNID-SAGE 를 사용하여 스펙트럼 분류를 수행했습니다. 광곡선 모델링은 두 가지 시나리오를 비교했습니다. 즉, cocoon 냉각을 포함하는 2 성분 킬로노바 모델과 충격 냉각 plus 방사성 붕괴를 포함하는 IIb 형 초신성 모델입니다. 우연한 정렬 확률을 평가하기 위해 체적 속도를 계산했습니다.

주요 결과

• 분광학적 분류: 초기 스펙트럼 (8 월 20 일) 은 대부분 특징이 없었습니다. 그러나 8 월 24 일과 28 일에 획득한 MUSE 스펙트럼은 뚜렷한 넓은 Hα P-Cygni 프로파일과 Hβ 특징을 드러냈습니다. SNID-SAGE 분석은 SN 2025ulz 를 IIb 형 초신성 (특히 SN 2011dh 및 SN 2016gkg 의 초기 단계와 일치) 으로 높은 신뢰도로 분류했습니다. 은하계 적색편이는 $z = 0.0849 \pm 0.0004$로 정제되었습니다.
• 광곡선 진화: 임시 천체는 초기 급격한 감쇠 후 약 15~20 일째에 두 번째 피크로 재밝아지는 모습을 보였습니다. 이 "충격 냉각 꼬리" 거동은 확장된 외피를 가진 IIb 형 초신성의 특징이지만, 일반적으로 단조로운 감쇠를 보이는 표준 킬로노바 모델에서는 예측되지 않습니다.
• 전파 제약: 전파 관측은 은하계의 항성 형성 (SFR $\sim 1.3 M_\odot$ yr$^{-1}$) 과 일치하는 일정한 플럭스 밀도를 검출했으며, 변동성이나 임시 잔광에 대한 증거는 없었습니다. 상한선은 관측 각도 $\theta \lesssim 12.5^\circ$ 내에서 GW170817 과 유사한 잔광을 배제했습니다.
• 모델링:
  • 킬로노바 적합: 2 성분 KN 모델은 처음 4 일간의 데이터만 적합시켰으며, 비물리적인 매개변수 (예: BNS 병합 기대치와 일치하지 않는 분출물 질량과 불투명도) 를 요구했고 이후의 재밝아짐을 설명하지 못했습니다.
  • 초신성 적합: IIb 형 초신성 모델 (충격 냉각 + 방사성 붕괴) 은 60 일까지의 전체 광곡선을 성공적으로 설명했습니다. 도출된 매개변수에는 분출물 질량 $2.2 M_\odot$, $^{56}$Ni 질량 $0.14 M_\odot$, 그리고 GW 유발 시간보다 약 1.6 일 이전의 폭발 시간이 포함되었습니다.
• 은하계: 은하계는 중간 정도의 태양 질량 미만 금속성 ($12 + \log(\text{O/H}) \approx 8.51$) 은하로, 항성 형성률은 $1.6 \pm 0.5 M_\odot$ yr$^{-1}$입니다.

결론 및 의의
이 논문은 SN 2025ulz 가 저신호 GW 후보 S250818k 와 우연한 시간적 및 공간적 일치를 보인 IIb 형 초신성이라고 결론지었습니다. 이 분류는 GW 유발과 양립할 수 없는 분광학적 특징, 광곡선 형태, 그리고 폭발 시간 제약을 기반으로 강력합니다. "슈퍼킬로노바" 시나리오는 데이터와 오컴의 면도날에 의해 강력히 기각됩니다.

저자들은 향후 GW 대응체 탐색에 중요한 함의를 강조합니다: 젊은 IIb 형 초신성은 킬로노바 탐색에서 주요 오염원입니다. 약 30~50% 의 IIb 형 초신성이 킬로노바의 급격한 초기 감쇠와 특징 없는 스펙트럼을 모방하는 충격 냉각 꼬리를 보입니다. 저자들은 전형적인 GW 국소화 영역 (S250818k 와 같은) 에서 우연히 이러한 간섭체가 나타날 확률을 약 30% 로 추정합니다. 이는 Rubin/LSST 와 같은 고빈도 관측 시대에는 이러한 우연한 정렬이 빈번해질 것이므로, 오염원과 진정한 킬로노바를 구별하기 위한 신속한 분광 후속 관측의 필요성을 강조합니다. 이 연구는 충격 냉각 꼬리가 킬로노바 탐색에서 가장 두드러진 오염원 중 하나임을 재확인합니다.