Similar Fermi-GBM sGRBs to GW/sGRB 170817A in MeV-GeV energies

기술 요약: MeV-GeV 에너지 대역에서 GW/sGRB 170817A 와 유사한 페르미-GBM sGRB

문제 제기
GW/sGRB 170817A 의 검출은 짧은 감마선 폭발 (sGRB) 이 쌍성 중성자별 (NS-NS) 병합에서 기원할 수 있음을 확인시켰습니다. 그러나 이 사건은 독특한 특성을 보였습니다: 낮은 광도 ($L_{iso} \approx 1.6 \times 10^{47}$ erg s$^{-1}$) 와 큰 축외 관측각 ($\sim 20^\circ$) 에서 관측된 독특한 이중 구성 요소 방출 구조 (비열적 피크 다음에 열적 꼬리가 따름) 입니다. 이론적 모델들은 많은 NS-NS 병합이 유사한 축외 sGRB 를 생성할 수 있음을 시사하지만, 대부분의 사건에 대한 거리 측정이 부재하기 때문에 역사적 데이터셋에서 이를 식별하는 것은 어렵습니다. 저자들은 전자기 (EM) 대응체의 GW 사건에 대한 발생률을 더 잘 추정하기 위해 페르미-GBM 카탈로그에서 sGRB 170817A 의 특정 스펙트럼 및 시간적 특징을 공유하는 sGRB 를 식별하는 것을 목표로 합니다.

방법론
이 연구는 페르미-GBM 데이터셋 (2008 년 7 월 17 일 – 2017 년 8 월 17 일) 에 적용된 다단계 필터링 및 클러스터링 접근법을 사용합니다:

1. 초기 선택: 저자들은 $T_{90} < 2.05$초인 635 개의 sGRB 를 선택했습니다.
2. 스펙트럼 필터링 (Zone-A): sGRB 170817A 의 첫 번째 피크는 특정 매개변수 ($E_{peak} \approx 215$ keV, $\alpha \approx 0.85$) 를 가진 콤프턴화 (Comp) 모델에 적합하므로, 샘플은 $E_{peak}$의 $1\sigma$ 상한선 (317 keV) 내에 있는 80 개의 사건으로 축소되었습니다.
3. 경도비 (HR) 계산: 광곡선은 두 시간 구역으로 나뉩니다: Zone-A ($T_{start}^{90}$부터 $T_{end}^{50}$까지) 와 Zone-B ($T_{end}^{50}$부터 $T_{end}^{90}$까지). 경도비 ($HR1$및$HR2$) 는 두 구역 모두에 대해 계산되었으며, 이는 50–300 keV 대역의 광자 수와 10–50 keV 대역의 광자 수 비율을 나타냅니다.
4. 이상치 제거: 마할라노비스 거리 지표를 사용하여 80 개 사건 샘플에서 두 개의 통계적 이상치 (sGRB 081101 및 sGRB 120814) 가 제거되어 78 개의 사건이 남았습니다.
5. 클러스터링: $HR1$대$HR2$매개변수 공간에 K-평균 클러스터링 알고리즘이 적용되었습니다. 최적의 클러스터 수 ($K=6$) 는 Kneedle 알고리즘 (관성 곡선에서 최대 곡률 감지) 을 사용하여 결정되었습니다.
6. 최종 선택: sGRB 170817A 와 다른 16 개의 사건을 포함한 클러스터 1 이 추가 분석되었습니다. 저자들은 두 가지 최종 기준을 적용했습니다:
   • 광곡선에서 뚜렷한 이중 피크의 존재.
   • Zone-B 를 온도 $kT < 12$ keV 인 흑체 (BB) 모델에 스펙트럼 적합.
7. 발생률 추정: 8 개의 선택된 사건 (이 중 6 개는 적색편이가 측정되지 않음) 에 대해, 저자들은 Amati 관계 ($E_{peak}$ 대 $E_{iso}$) 를 사용하여 의사 적색편이 ($z_p$) 를 추정했습니다. 이러한 거리는 사건 발생률을 계산하고 LIGO 관측 기간 ($O1$부터 $O5$까지) 에 대한 GW+sGRB 검출 수를 예측하는 데 사용되었습니다.

주요 결과

• 유사 사건의 식별: 이 연구는 페르미-GBM 카탈로그에서 이중 피크 광곡선 구조와 스펙트럼 진화 (비열적 방출 다음에 열적 방출) 측면에서 sGRB 170817A 와 매우 유사한 8 개의 sGRB를 식별했습니다. 식별된 사건은 다음과 같습니다: sGRB 150805, sGRB 150101, sGRB 131128, sGRB 131004, sGRB 130808, sGRB 120524, sGRB 090108, 그리고 sGRB 081122.
• 스펙트럼 특성: 분석은 이러한 사건의 첫 번째 구역 (Zone-A) 이 콤프턴화 모델을 따르는 반면, 두 번째 구역 (Zone-B) 은 sGRB 170817A 에서 관측된 열적 꼬리와 일치하는 $\sim 7.6$에서 $12.0$ keV 범위의 온도를 가진 흑체 모델로 잘 적합됨을 확인했습니다.
• 사건 발생률 예측: 유도된 광도 거리 (의사 적색편이 포함) 를 사용하여 저자들은 LIGO 관측 기간에 대한 예상 GW+sGRB 사건 수를 계산했습니다:
  • $O1$ 및 $O2$: 각각 약 0.58 개와 1.36 개의 사건이 예상됩니다. 이는 $O2$ 기간 중 sGRB 170817A 관측과 일치합니다.
  • $O3$ 및 $O4$: 합계 예상치는 약 6 개에서 8 개 사건입니다.
  • $O5$: 다가오는 $O5$ 기간 (예상 기간 $\sim 884$일, 거리 범위 240–330 Mpc) 으로 외삽할 때, 저자들은 약 7.49 $\pm$ 3.71개의 GW+sGRB 사건을 예측합니다.

의의 및 주장
이 논문은 이러한 8 개의 "쌍둥이"sGRB 를 식별함으로써 관측 가능한 EM 대응체를 생성하는 NS-NS 병합의 발생률을 정교하게 할 수 있는 통계적으로 견고한 샘플을 제공한다고 주장합니다. 저자들은 이러한 사건의 존재가 많은 NS-NS 병합이 높은 경사로 발생하여 관측 가능하지만 축내 폭발과 구별되는 축외 sGRB 를 생성한다는 가설을 지지한다고 논합니다.

이 연구는 예측된 사건 발생률에서의 편차 (예: 향후 기간에서 현저히 적거나 많은 사건이 검출되는 경우) 가 거리와 함께 이러한 사건의 진화나 제트의 기하학에 대한 현재의 이해를 도전할 수 있음을 강조합니다. 이 작업은 광곡선 형태뿐만 아니라 비열적에서 열적 구성 요소로의 스펙트럼 전이에 특히 초점을 맞춤으로써 페르미-GBM 협력팀과 Swift 팀의 이전 검색들에 대한 통계적 보완 역할을 합니다. 저자들은 $O4$ 기간이 끝날 무렵 거의 5 개의 GW+sGRB 사건이 있을 수 있으며, 이 수치에서의 편차는 거리와 함께 이러한 사건의 진화에 대한 우려를 제기할 것이라고 결론지었습니다.