Inferring Neutron Star Nuclear Properties from Gravitational-Wave and Gamma-Ray Burst Observations

기술 요약: 중력파 및 감마선 폭발 관측을 통한 중성자별 핵 성질의 추론

문제 정의
최근 킬로노바 방출(예: GRB 211211A, GRB 230307A)을 동반한 장기 지속 감마선 폭발(lGRB) 관측은 장기 지속 감마선 폭발이 오직 거대 질량 별의 붕괴로부터만 발생한다는 전통적인 패러다임에 도전하고 있다. 만약 이러한 사건들이 이중 중성자별(BNS) 또는 중성자별-블랙홀(NSBH) 병합에서 기원한다면, 폭발의 지속 시간을 결정하는 물리적 성질은 무엇인가라는 중요한 질문이 제기된다. 이론적 모델들은 병합의 총 질량이 병합 잔해의 수명을 결정한다고 제안한다. 가벼운 병합은 단기 지속 감마선 폭발(sGRB)을 일으키는 장수(long-lived) 중성자별 잔해를 형성할 수 있는 반면, 무거운 병합은 단기 수명(short-lived) 잔해를 생성하여 블랙홀(BH)로 붕괴하며, 이는 거대 강착 원반을 통해 lGRB를 일으킬 수 있다. 그러나 이러한 결과를 가르는 구체적인 질량 임계값과 톨만-오펜하이머-볼코프(TOV) 질량($M_{\text{TOV}}$) 사이의 관계, 즉 근본적인 핵 성질과의 관계는 아직 관측적으로 규명되지 않았다.

방법론
저자들은 LIGO-Virgo-Kaka(LVK) 협업체의 중력파(GW) 병합률을 관측된 sGRB 및 킬로노바 관련 lGRB의 발생률과 비교함으로써 특성 질량 척도를 추론하기 위한 베이지안 프레임워크를 개발하였다.

1. 특성 질량 정의: 본 연구는 TOV 질량($M_{\text{TOV}}$)에 대한 세 가지 임계 총 질량을 다음과 같이 매개변수화한다:

   • $M_{\text{VL}} = k_1 M_{\text{TOV}}$: 초장수(very long-lived) 잔해와 장수 잔해 사이의 경계.
   • $M_{\text{LS}} = k_2 M_{\text{TOV}}$: 장수 잔해와 단기 수명 잔해 사이의 경계. $M_{\text{VL}} < M_{\text{tot}} \le M_{\text{LS}}$인 병합은 sGRB를 생성하는 것으로 가설을 세운다.
   • $M_{\text{PC}} = k_3 M_{\text{TOV}}$: 단기 수명 잔해와 즉각적 붕괴(prompt collapse) 사이의 경계. $M_{\text{LS}} < M_{\text{tot}} \le M_{\text{PC}}$인 병합은 lGRB를 생성할 가능성이 있다. $M_{\text{tot}} > M_{\text{PC}}$인 병합은 즉각적 붕괴를 겪으나, 질량비($q$)가 충분한 원반 형성을 허용하는 경우($q$가 BNS의 경우 $1.2 \le q \le 3$, NSBH의 경우 $q \le 3$)에만 lGRB가 생성된다.

2. 데이터 소스:

   • GW 발생률: LVK GWTC-4 카탈로그(Abac et al. 2025a)와 FullPop-4.0 인구 모델(Abac et al. 2025b)은 구성 질량에 따른 국소 병합률을 제공한다.
   • GRB 발생률: sGRB의 고유 발생률($R_{\text{sGRB}}$)과 발생률 비 $\eta = R_{\text{lGRB}}/R_{\text{sGRB}}$는 상당한 불확실성을 가진 관측 제약 조건으로 취급되며, 정규 분포를 통해 모델링된다.

3. 추론: 저자들은 베이즈 정리를 사용하여 매개변수 집합 $\vec{\Theta} = (k_1, k_2, k_3, M_{\text{TOV}})$에 대한 사후 분포를 추론한다. 우도 함수(likelihood function)는 질량 임계값을 통해 유도된 GW 유도 GRB 발생률을 관측된 GRB 발생률과 비교한다. $M_{\text{TOV}}$와 TOV 질량에서의 중성자별 반지름($R_{\text{TOV}}$)에 대한 사전 확률(priors)은 펄서 및 GW 관측 결과에 기반한 Legred et al. (2021)의 데이터를 바탕으로 추출된다.

주요 결과

• 특성 질량비 ($k_2$): 분석 결과 $k_2 = 1.36^{+0.08}_{-0.09}$ (68% 신뢰 구간)의 중앙값을 얻었으며, $k_2 > 1.24$ (90% 신뢰 구간)이다. 이는 장수 잔해에서 단기 수명 잔해로의 전이가 약 $1.36 M_{\text{TOV}}$에서 발생함을 나타낸다.
• 잔해 수명의 함의: 높은 $k_2$ 값은 무거운 중성자별이 병합 후 즉시 붕괴하는 것이 아니라, 수백 밀리초 동안 연장된 기간 동안 생존할 수 있음을 시사한다.
• 상관관계: 연구는 $k_2$와 $k_3$ 사이의 강한 상관관계, 그리고 $k_2$와 $M_{\text{TOV}}$ 사이의 역상관관계를 식별하였다. 이 역상관관계 덕분에 $k_2$에 대한 제약 조건을 $M_{\text{TOV}}$의 상한치로 직접 매핑할 수 있다.
• 강건성: $k_2$의 추론은 가정된 sGRB 발생률($30$~$2000 \text{ Gpc}^{-3}\text{yr}^{-1}$) 및 발생률 비($\eta$가$0.01$~$3$인 경우)의 넓은 범위에서 강건하게 유지된다. 질량비 임계값($q_{\text{low}}$및$q_{\text{high}}$)의 변화와$M_{\text{VL}}$경계의 포함 여부는 추론된$k_2$에 미미한 영향을 미친다.

의의 및 주장
본 논문은 GW 병합 인구와 GRB 지속 시간을 연결함으로써 중성자별 상태 방정식(EoS)을 제약하는 새롭고 독립적인 방법을 확립한다. sGRB와 lGRB의 기원을 구분하는 특성 질량 척도($1.36 M_{\text{TOV}}$)를 식별함으로써, 저자들은 중성자별의 최대 질량에 대한 직접적인 관측적 제약을 제공한다.

저자들은 자신들의 결과가 $1.3$과$1.4 M_{\text{TOV}}$ 사이의 전이를 예측하는 수치 상대론 시뮬레이션(예: Perna et al. 2025)과 일치한다고 언급한다. 또한 현재의 결과가 관측 불확실성에 대해 강건하지만, 향후 GW 및 GRB 발생률 측정(특히 GRB와 GW 이벤트의 공동 검출)이 개선되면 특성 질량과 $M_{\text{TOV}}$ 사이의 상관관계가 더욱 긴밀해져 $M_{\text{TOV}}$에 대한 불확실성을 $\sim 0.1 M_{\odot}$ 수준까지 줄일 수 있을 것이라고 강조한다. 이 접근법은 기존 기술을 보완하며 중성자별 EoS에 대한 더 엄격한 한계를 설정하는 경로를 제시한다.