Search for Long-Transient Gravitational Waves from Supernova SN2023ixf using GFH-v2 Pipeline

기술 요약: GFH-v2를 이용한 SN 2023ixf의 장기 과도 중력파 탐색

문제 및 동기
핵붕괴 초신성(CCSNe)에서 형성되는 빠르게 회전하는 신생 마그네타는 장기 과도 연속 중력파(tCWs)를 방출할 것으로 예측된다. 붕괴 자체와 관련된 단기 지속 버스트(burst)와 달리, 이러한 신호는 새롭게 형성된 비축퇴 중성자별 또는 마그네타의 스핀다운(spin-down)으로부터 발생한다. 이 방출은 회전 에너지의 손실에 의해 구동되며, 멱법칙(power law)에 따라 감소하는 주파수와 시간에 따라 감쇠하는 변형 진폭(strain amplitude)을 생성한다. 이러한 신호를 검출하는 것은 중성자별의 탄생, 내부 구조 및 스핀 진화에 대한 직접적인 통찰을 제공할 것이다.

인근 은하 M101(거리 $\approx$ 6.85 Mpc)에서 발견된 제2형 핵붕괴 초신성인 SN 2023ixf는 이러한 탐색에 있어 독특한 기회를 제공한다. 이 사건의 근접성과 전자기적 검출의 정밀한 타이밍(약 2시간 이내의 창)은 지향성 탐색(directed search)을 가능하게 한다. 결정적으로, 이 사건의 시점은 LIGO의 엔지니어링 런 15(ER15)와 겹치며, 이는 네 번째 관측 실행(O4) 직전의 핸포드(H1) 및 리빙스턴(L1) 검출기로부터 얻은 일치하는 데이터셋을 제공한다.

방법론
본 분석은 멱법칙 주파수 진화에 특화된 일반화된 주파수 허프 변환(Generalized Frequency Hough transform)에 기반한 계층적 준결맞음(semi-coherent) 탐색 방법인 GFH-v2 파이프라인을 사용한다.

• 신호 모델: 탐색은 주로 중력파 방출을 통해 스핀다운되는 강체 비전정(non-precessing) 중성자별을 대상으로 한다. 이는 브레이킹 인덱스(braking index) $n=5$에 해당한다. 중력파 주파수 $f_{gw}$는 $\dot{f}_{gw} = -k f_{gw}^n$로 진화하며, 변형 진폭 $h_0(t)$는 $t^{-1/2}$로 감쇠한다. 탐색은 표준적인 중성자별 파라미터($M=1.4 M_\odot$, $R=12$ km)를 가정하며, 전자기적 국지화에 기반한 고정된 하늘 위치를 사용한다.
• 데이터: 분석에는 ER15 기간 동안 H1과 L1에서 얻은 일치하는 과학 세그먼트를 사용하며, 보수적인 온소스 윈도우(on-source window) 내에서 총 약 48,619초(0.56일)의 데이터를 포함한다. 데이터 품질은 600–2000 Hz 범위에서 O4a 감도와 유사함을 확인하였다.
• 탐색 구성: 파라미터 공간은 초기 주파수($f_0$) 600에서 2000 Hz(단계 10 Hz) 및 타원율($\epsilon$) $3 \times 10^{-4}$에서 $3 \times 10^{-3}$ (5개의 하위 구간으로 분할)를 포괄한다. 이는 총 700개의 독립적인 탐색 구성을 생성한다.
• 파이프라인 단계:
  1. 데이터 준비: 단시간 푸리에 변환 데이터베이스(SFDB) 데이터에 밴드패스 필터를 적용하고 복소 해석적 시계열로 변환한다.
  2. 피크맵(Peakmap) 생성: 등가 전력 스펙트럼(equalized power spectrum)에서의 국소 극대값을 선택하여 시간-주파수 피크맵을 형성한다.
  3. 허프 변환(Hough Transform): 피크맵을 $(t, f)$ 좌표에서 변환된 좌표 $(t, x)$로 매핑한다 (여기서 $x = 1/f^{n-1}$). 이 공간에서 멱법칙 진화는 선형이 된다. 허프 변환은 이를 $(x_0, k)$ 파라미터 공간(초기 주파수 및 스핀다운 파라미터)으로 매핑한다.
  4. 후보 선택: 후보는 국소 배경 대비 과잉 계수를 가진 빈(bin)으로 식별된다. 비정상적 노이즈를 처리하기 위해 국소적으로 추정된 노이즈 통계량(10 Hz 간격 내의 중앙값 및 분산)을 사용하여 임계 비율(CR)을 계산한다.
  5. 일치(Coincidence): H1과 L1의 후보들을 $(x_0, k)$ 파라미터 공간에서 비교한다. 거리 메트릭 $d_{coin} \le 2$를 만족하면 한 쌍의 후보를 일치하는 것으로 간見한다.
  6. 후속 처리(Follow-up): $CR \ge 3.5$인 일치 후보들은 헤테로다인(heterodyne) 후속 처리를 거친다. 예상되는 위상 진화를 데이터에서 제거한 후, 신호의 유의성이 증가하는지 확인하기 위해 원래보다 더 긴 결맞음 시간($3\times$ 및 $6\times$)을 사용하여 신호를 재분석한다. 이는 실제 천문학적 신호의 존재를 나타낸다.

주요 결과

• 후보 선택: 일치 요구 사항과 CR 임계값($CR \ge 3.5$)을 적용한 결과, 131개의 일치 후보가 남았다.
• 후속 처리 결과: 131개의 후보 중 어느 것도 헤테로다인 후속 처리 단계에서 예상되는 임계 비율의 증가를 보이지 않았다. 결과적으로, 어떤 후보도 두 번째 후속 처리 단계로 승격되지 않았다.
• 검출: 본 분석은 ER15 데이터에서 SN 2023ixf와 연관된 장기 과도 중력파 방출의 증거를 발견하지 못했다.
• 상한선(Upper Limits): 최대 검출 가능 거리에 대한 상한선이 초기 주파수 및 타원율의 함수로서 설정되었다. 가장 높은 타원율 구간($\epsilon \in [2.5 \times 10^{-3}, 3.0 \times 10^{-3}]$)에 대해, 90% 효율 상한선은 분석된 대부분의 주파수 대역에서 약 1 ~ 2.5 Mpc 범위를 나타냈다. 낮은 타원율의 경우, 제한된 데이터 기간 및 낮은 신호 대 잡음비로 인해 상한선을 결정할 수 없거나 더 작은 거리로 제한되었다.

의의 및 주장
저자들은 도출된 상한선(1–2.5 Mpc)이 실제 M101까지의 거리($\approx$ 6.85 Mpc)보다 낮지만, 본 연구가 가까운 핵붕괴 초신성을 위한 지향성 탐색에 GFH-v2 파이프라인을 적용한 첫 사례라고 밝히고 있다. 본 연구는 실제 커미셔닝 시대의 검출기 데이터(ER15)에 대한 파이프라인의 성능을 성공적으로 규명하였다. 이는 GFH-v2 방법이 마그네타 형성 이벤트로부터 발생하는 tCW를 탐색하는 데 효과적으로 배치될 수 있음을 입증하며, 향로의 관측 실행(O5) 및 차세대 검출기를 사용하는 미래 탐색을 위한 프레임워크를 구축한다. 결론적으로, 본 탐색은 ER15 데이터에 대한 장기 과도 탐색의 직접적인 민감도 추정치를 제공하며 향후 노력을 위한 벤치마크 역할을 한다.