Probing millisecond magnetar formation in binary neutron star mergers through X-ray follow-up of gravitational wave alerts

기술 요약: 중력파 경보에 따른 X-선 후속 관측을 통한 밀리초 마그네타 형성 탐사

문제 제기
이진 중성자별(BNS) 병합 이후 형성되는 잔해의 본질은 여전히 불확실하다. 가장 질량이 큰 병합은 즉각적으로 블랙홀(BH)로 붕괴할 것으로 예상되지만, 단기 감마선 폭발(sGRB) 잔광의 X-선 광도 곡선은 빠르게 회전하며 강한 자기장을 가진 중성자별(NS)—즉, 밀리초 마그네타—가 상당수의 트랜지언트(transient)에서 중심 엔진 역할을 할 수 있음을 시사한다. 안정적인 NS, 스핀다운 후 붕괴하는 초질량 NS(SMNS), 또는 즉각적인 BH 붕괴를 구분하는 것은 초고밀도 물질의 상태 방정식(EoS)과 감마선 폭발의 물리학을 이해하는 데 매우 중요하다. 그러나 중력파(GW) 검출만으로는 잔해의 본질을 항상 규명할 수 없고, 전자기파(EM) 대응 천체는 흔히 이젝트(ejecta)에 의해 가려지거나 오프-액시스(off-axis) 방향에서 관측되기 때문에 이러한 잔해를 식별하는 것은 도전적인 과제이다. 본 논문은 현재와 미래의 X-선 관측 성능을 평가함으로써, 중력파 경보에 대한 X-선 후속 관측을 통해 밀리초 마그네타 잔해를 탐지할 수 있는 역량을 평가하고, 특히 중력파 검출기 구성과 X-선 기기의 민감도가 탐지 전망에 미치는 영향을 분석하는 것을 목표로 한다.

방법론
저자들은 중력파 신호 생성과 상세한 X-선 광도 곡선 모델링을 결합한 종합적인 시뮬레이션 프레임워크를 개발하였다:

1. 합성 BNS 집단 및 중력파 검출 가능성:

   • LIGO-Virgo-KAGRA(LVK) O3 실행에서 유도된 질량 분포(특히 "Peak" prior)를 사용하여 합성 BNS 시스템 집단을 생성하였다.
   • 중력파 신호는 현재(LVK O4 및 O5) 및 차세대(Einstein Telescope [ET] 및 Cosmic Explorer [CE]) 간섭계에 대해 PyCBC 소프트웨어를 사용하여 시뮬레이션되었다.
   • X-선 후속 관측이 용이하도록 정밀하게 국소화된 이벤트($\Delta\Omega \le 50$ deg$^2$)를 필터링하여 하늘 위치 국소화 정확도를 계산하였다.
   • 다양한 EoS 모델(예: Shen, BSk21, DDME2)을 고려하여 최대 비회전 중성자별 질량($M_{TOV}$)과 그에 따른 잔해의 운명을 결정하였다.

2. X-선 광도 곡선 시뮬레이션:

   • 전용 수치 코드를 사용하여 자기 쌍극자 스핀다운, 프로펠러 효과 및 중력파 방출을 고려한 신생 밀리초 마그네타의 X-선 방출을 모델링하였다.
   • 이젝트 모델링: 이젝트를 동역학적 성분과 병합 후 성분의 중첩으로 취급하는 새로운 다성분 이젝트 프레임워크를 구현하였다. 이젝트의 질량과 불투명도는 최첨단 수치 상대론 시뮬레이션(Nedora et al. 2021a; Mösta et al. 2020a)의 다항식 피팅을 통해 도출되었으며, 이를 통해 이젝트 특성을 BNS 궤도 매개변수 및 EoS와 직접 연결하였다.
   • 기하학적 구조: 시뮬레이션은 두 가지 시선 방향을 구분하였다: "자유 구역"(제트 축 근처, 낮은 불투명도)과 "포획 구역"(적도 부근, 높은 불투명도). 모델은 이젝트에 의한 마그네타 복사의 재처리 과정을 고려하여 광학적 깊이와 흑체 복사를 계산하였다.
   • 붕괴 기준: 시뮬레이션은 스핀다운 진화를 추적하여, EoS 및 임계 회전 주기에 기반한 블랙홀로의 붕괴 시간을 결정하였다.

3. 관측 제약 조건:

   • 예측된 플럭스를 SVOM/MXT, Einstein Probe (EP)/FXT, 그리고 EP/WXT를 포함한 현재 및 미래 X-선 기기의 민감도 임계값과 비교하였다.
   • BNS 병합율 $R = 100$ Gpc$^{-3}$ yr$^{-1}$을 가정하여 검출률을 추정하였으며, 통계적 불확실성은 비매개변수 부트스트랩 방법을 통해 도출하였다.

주요 기여

• 통합 시뮬레이션 프레임워크: 본 논문은 단순한 매개변수화를 넘어, 현실적인 EoS 의존적 이젝트 모델링을 포함하여 중력파 검출 확률과 상세한 X-선 광도 곡선 진화를 연결하는 자기 일관적인 시뮬레이션을 제시한다.
• 시선 방향 편향 분석: 연구는 중력파 검출 가능성과 전자기파 가시성 사이의 트레이드오프를 정량화한다. "자유 구역" 시스템은 X-선에서 더 밝지만 먼 거리에서 중력파를 통한 국소화가 어려운 반면, "포획 구역" 시스템은 더 수가 많고 가까운 거리에서 국소화하기 쉽지만 높은 이젝트 불투명도의 영향을 받는다는 점을 입증한다.
• 최적 검출 창: 저자들은 최대 검출 가능성을 보이는 특정 시간 창을 식별하였으며, 우주론적 시간 지연으로 인해 차세대 검출기에서는 피크 검출 시점이 병합 후 약 2시간에서 3~4시간으로 이동함을 주목하였다.
• 기기 전략: 연구는 관측 전략을 제시하며, 중력파 검출기 세대에 따라 포획 구역 이벤트에는 협대역 민감 기기(예: MXT, FXT)를, 자유 구역 이벤트에는 광대역 모니터(예: WXT)를 사용하는 것을 권장한다.

결과

• 형성율: 시뮬레이션 결과, EoS(특히 $M_{TOV}$)에 따라 다르지만, BNS 병합의 2%에서 16% 사이에서 장수명 밀리초 마그네타가 형성된다.
• 검출률 (현재 세대):
  • 경직된(stiff) EoS(DDME2)를 가정할 때, LVK O5 실행에서 밀리초 마그네타의 예측 검출률은 $1.0^{+0.3}_{-0.3}$ yr$^{-1}$이다. 만약 EoS가 연질(soft, Shen)이라면 현재 기기로는 검출이 기대되지 않는다.
  • 최적의 검출 가능성은 병합 후 약 2시간 후에 나타난다.
  • 잘 국소화된 이벤트의 경우, 병합 1시간 후 검출 가능한 마그네타의 비율은 0%에서 6% 사이이다 (가장 경직된 EoS의 경우 최대 ~12%).
• 검출률 (차세대 세대):
  • ET와 CE를 사용할 경우, 검출률은 최대 3 자릿수 증가하여, ET+CE의 경직된 EoS 시나리오에서 연간 ~750개에 달할 수 있다.
  • 차세대 검출기의 경우, 피크 검출 가능성이 병합 후 3~4시간으로 이동한다.
• 자기장 의존성:
  • 자기장이 $\sim 7.5 \times 10^{14}$ G인 마그네타는 병합 1시간 후 가장 높은 검출 가능성(잘 국소화된 이벤트의 16%)을 보인다.
  • 더 높은 자기장($10^{16}$ G)은 빠른 스핀다운과 조기 붕괴(수 분 이내)를 초래하여 관측 창을 제한한다.
  • 더 낮은 자기장($10^{14}$ G)은 수명이 길지만 초기 광도가 낮아, 이젝트가 광학적으로 얇아지는 병합 수 시간 후에야 검출될 수 있다.
• EoS 제약: 결과는 밀리초 마그네타의 검출(또는 미검출)이 EoS에 대한 제약을 줄 수 있음을 시사한다. 구체적으로, 장기 지속 X-선 플래토(plateau)의 부재는 연질 EoS 또는 낮은 자기장을 지지할 것이다.

의의 및 주장
저자들은 밀리초 마그네타의 스핀다운 방출 검출이 EoS가 충분히 경직되어 있다면 현재의 중력파 네트워크(O5 실행 기간 포함)로도 "손에 닿을 수 있는 범위 내에 있음"을 주장한다. 이들은 단 한 번의 검출이라도 병합 물리, 감마선 폭발의 중심 엔진, 그리고 중성자별 물리학에 심오한 영향을 미칠 것이라고 강조한다.

본 논문은 병합율, EoS, 기기 응답의 불확실성으로 인해 검출이 어렵다는 점을 겸허히 인정하면서도, 제안된 시뮬레이션 기반 관측 전략(특히 병합 후 2~4시간의 타이밍 및 시선 방향에 따른 기기 선택)이 성공 확률을 크게 높인다고 결론짓는다. 이 연구는 BNS 잔해의 본질을 해결하기 위한 지속적인 이론적 및 관측적 노력의 필요성을 강조한다.