Wide Jets or Low Rates: Reconciling Short GRB and Gravitational-Wave Neutron Star Merger Rates

논문 개요

이 연구는 중력파 (GW) 관측과 짧은 감마선 폭발 (sGRB) 관측 데이터를 결합하여, 중성자별 - 중성자별 (BNS) 병합과 중성자별 - 블랙홀 (NSBH) 병합이 sGRB 의 주요 기원인지, 그리고 현재 관측된 sGRB 발생률과 LIGO/Virgo/KAGRA (LVK) 가 추정한 중력파 병합률이 서로 모순되는지 여부를 규명하는 것을 목표로 합니다. 특히 4 차 관측 주기 (O4) 에서 BNS 병합이 추가로 발견되지 않아 추정이 낮아진 LVK 의 BNS 병합률과, 문헌에서 보고된 다양한 sGRB 국소 발생률 사이의 일관성을 검증합니다.

1. 연구 배경 및 문제 제기

• 모순의 발생: GW170817 사건은 BNS 병합이 sGRB 를 생성할 수 있음을 입증했으나, 최근 LVK 의 4 차 관측 주기 (O4) 동안 고신뢰도 BNS 병합 신호가 추가로 발견되지 않았습니다. 이로 인해 LVK 가 추정한 국소 BNS 병합률 ($R_{BNS}$) 은 기존 추정치보다 현저히 낮아졌습니다 (약 $39^{+70}_{-29} \text{ Gpc}^{-3} \text{ yr}^{-1}$).
• sGRB 발생률의 불확실성: 반면, 전자기파 (EM) 관측을 기반으로 한 sGRB 국소 발생률 ($R_{sGRB}$) 은 문헌에 따라 약 $1 \sim 13 \text{ Gpc}^{-3} \text{ yr}^{-1}$까지 광범위하게 보고됩니다. (Coward et al. 2012 의 경우 중앙값은 $8 \text{ Gpc}^{-3} \text{ yr}^{-1}$이나, 상한 $1\sigma$ 구간은 약 $13 \text{ Gpc}^{-3} \text{ yr}^{-1}$에 달함).
• 핵심 질문: 낮은 중력파 병합률로 높은 sGRB 발생률을 설명할 수 있는가? 이를 위해서는 제트 (Jet) 의 각도 구조 (넓은 제트) 나 제트 발사 효율이 어떻게 되어야 하는가? NSBH 병합이 이 모순을 해결할 수 있는가?

2. 방법론 (Methodology)

저자들은 베이지안 추론 (Bayesian Inference) 프레임워크를 사용하여 sGRB 발생률과 중력파 병합률을 연결했습니다.

• 데이터 소스:
  • 중력파: LVK 의 4 차 관측 주기 (O4) 데이터 기반의 FullPop-4.0 모델을 사용. O4 기간 중 BNS 미발견 사실을 반영하여 BNS 병합률을 보정 (스케일링) 했습니다.
  • sGRB: 문헌에서 보고된 다양한 국소 sGRB 발생률 ($\lambda_{sGRB}$) 을 시나리오로 설정:
    1. 높은 경우: $8 \text{ Gpc}^{-3} \text{ yr}^{-1}$ (Coward et al. 2012 기준, GRB 080905A 포함).
    2. 중간/보수적 경우: $3 \text{ Gpc}^{-3} \text{ yr}^{-1}$ (GRB 080905A 제외 시).
    3. 낮은 경우: $1 \text{ Gpc}^{-3} \text{ yr}^{-1}$ (Ghirlanda et al. 2016, Salafia et al. 2023 등 최신 연구 기준).
• 모델링 시나리오:
  1. BNS-only 시나리오: 모든 sGRB 가 BNS 병합에서 기원한다고 가정.
  2. BNS + NSBH 시나리오: sGRB 가 BNS 와 NSBH 병합 모두에서 기원한다고 가정.
• 파라미터:
  • 제트 반개방각 ($\theta_j$): 제트의 기하학적 빔링 (beaming) 결정.
  • 제트 발사 비율 ($f_s$): 병합 사건 중 실제로 상대론적 제트를 성공적으로 발사하는 비율.
  • 중력파 기반 시뮬레이션 (GWTC-4) 을 사용하여 NSBH 병합 중 제트 발사 가능 비율에 대한 물리적 제약을 설정.
• 방정식: 관측된 sGRB 발생률은 $R_{sGRB} = f_s \cdot \eta \cdot R_{merger}$로 표현되며, 여기서 $\eta = 1 - \cos\theta_j$는 기하학적 빔링 인자입니다.

3. 주요 결과 (Key Results)

가. BNS-only 시나리오

• 높은 sGRB 발생률 ($\sim 8 \text{ Gpc}^{-3} \text{ yr}^{-1}$) 의 경우:
  • 현재 LVK 의 낮은 BNS 병합률과 일치하려면 매우 넓은 제트 ($\theta_{BNS} \gtrsim 25^\circ$) 와 높은 제트 발사 효율 ($f_{s,BNS} \approx 0.73$) 이 필요합니다.
  • 이는 후광 (afterglow) 관측을 통해 일반적으로 추정되는 좁은 제트 ($\sim 6^\circ$) 와는 상충되며, 물리적으로 불합리한 시나리오로 간주됩니다.
• 낮은 sGRB 발생률 ($\sim 1 - 3 \text{ Gpc}^{-3} \text{ yr}^{-1}$) 의 경우:
  • sGRB 발생률이 $1 \sim 3 \text{ Gpc}^{-3} \text{ yr}^{-1}$로 낮아지면, 좁은 제트 가정이 가능해집니다. 이 regime 에서 추정된 제트 개방각의 중앙값은 약 $12^\circ$에서 $19^\circ$ 사이입니다.
  • 특히 Table 1 의 $1\sigma$ 신뢰구간 하한을 고려할 때, 좁은 제트 regime 에서 $\theta_{BNS} \approx 3^\circ$까지의 값이 가능합니다.
  • 이러한 조건 하에서 LVK 의 BNS 병합률과 완벽하게 일치합니다. 특히 GRB 080905A 와 같이 불확실한 적색편이를 가진 사건을 제외하거나 최신 연구의 낮은 발생률을 채택할 때, BNS 병합이 sGRB 의 유일한 기원이라는 가정이 중력파 데이터와 모순되지 않음을 보여줍니다.

나. NSBH 병합의 기여도

• NSBH 병합이 sGRB 에 기여할 수는 있으나, 그 기여도는 BNS 에 비해 부차적입니다.
• 시뮬레이션 결과, NSBH 병합 중 제트를 발사할 수 있는 비율 ($f_{s,NSBH}$) 은 약 $23\%$ 이하로 제한됩니다.
• 결과적으로 NSBH 채널은 관측된 sGRB 의 약 $2\% \sim 16\%$ 만 기여하는 것으로 추정되며, 높은 sGRB 발생률 ($> 7 \text{ Gpc}^{-3} \text{ yr}^{-1}$) 을 설명하기 위해 NSBH 가 결정적인 역할을 하지는 못합니다.

다. 시스템적 불확실성 및 교정

• 제트 각도 편향 (Selection Bias): 기존 제트 각도 추정은 '제트 브레이크 (jet break)'가 관측된 좁은 제트 사건에 편향되어 있을 가능성이 큽니다. 넓은 제트는 브레이크가 늦게 나타나거나 감지되지 않아 샘플에서 누락될 수 있습니다. 만약 실제 sGRB 군집에 넓은 제트가 상당 부분 포함된다면, 현재 GW 제한 내에서 sGRB 발생률 $7 \text{ Gpc}^{-3} \text{ yr}^{-1}$까지도 설명 가능합니다.
• 적색편이 오차 및 발생률 과대평가: sGRB 발생률 추정치와 중력파 데이터 간의 불일치는 일부 문헌에서 국소 sGRB 발생률이 과대평가되었기 때문일 가능성이 높습니다. 이는 주로 부적절한 적색편이 할당과 발생률의 적색편이 진화 (redshift evolution) 를 올바르게 고려하지 않았기 때문입니다. GRB 080905A 는 이러한 요인 중 하나일 뿐, 더 큰 표본을 사용하는 연구들에서 불일치를 주도하는 주된 요인은 아닙니다.
• 적색편이 진화: sGRB 발생률은 우주적 거리 (적색편이) 에 따라 진화할 수 있습니다. LVK 가 관측하는 국소 우주 ($z \approx 0$) 와 sGRB 가 주로 관측되는 중적색편이 ($z \approx 0.5$) 영역 사이의 발생률 차이를 고려하면 모순이 완화될 수 있습니다.

4. 결론 및 의의

• 주요 결론:

  1. 낮은 sGRB 발생률 가정이 타당함: 현재 LVK 의 BNS 병합률 제한과 sGRB 관측을 조화시키기 위해서는, sGRB 발생률이 $1 \sim 3 \text{ Gpc}^{-3} \text{ yr}^{-1}$ 수준으로 낮아야 합니다. 이는 최신 문헌의 여러 추정치와 일치합니다.
  2. 넓은 제트 필요성: 만약 sGRB 발생률이 $8 \text{ Gpc}^{-3} \text{ yr}^{-1}$로 높다면, 물리적으로 비현실적인 넓은 제트 ($\theta > 25^\circ$) 를 가정해야 하므로, 이는 sGRB 발생률 자체가 과대평가되었거나 제트 구조에 대한 편향이 존재함을 시사합니다.
  3. BNS 의 지배적 역할: NSBH 병합은 sGRB 의 일부만 설명할 뿐, BNS 병합이 sGRB 의 주된 기원이라는 결론은 변하지 않습니다.
  4. 일관성: 적절한 제트 각도 분포 (넓은 제트 포함) 와 낮은 sGRB 발생률을 가정하면, 중력파와 전자기파 관측은 BNS 병합이 sGRB 의 주요 기원이라는 가설 하에 일관된 그림을 제공합니다.

• 의의:
  이 연구는 중력파 천문학의 정밀도가 높아짐에 따라 (BNS 병합률 감소), sGRB 발생률 추정치와 제트 물리 (각도, 효율) 에 대한 재평가가 필요함을 강조합니다. 또한, sGRB 발생률 추정의 불확실성 (특히 적색편이와 제트 브레이크 관측 편향) 이 GW-EM 일관성 논의에서 얼마나 중요한지를 정량화했습니다. NSBH 병합만으로는 sGRB 발생률이 높다고 가정하든 낮다고 가정하든 간에, 중력파와 sGRB 관측 사이의 긴장 (tension) 을 해결할 수 없습니다. 향후 더 많은 중력파 및 전자기파 동시 관측을 통해 제트 구조와 기원 천체의 비율에 대한 제약이 더욱 강화될 것으로 기대됩니다.