[논문 리뷰] The maximum mass of neutron stars may be higher than expected: an inference from binary systems

L. S. Rocha, R. R. A. Bachega외 2명|arXiv (Cornell University)|2021. 07. 15.

Seismology and Earthquake Studies참고 문헌 4인용 수 7

한 줄 요약

이 연구는 96개의 중성자별 이중성계 샘플에 대한 베이지안 추론을 통해 중성자별의 최대 질량을 재평가하며, 마르코프 체인 몬테카를로(MCMC) 및 사후 예측 검증을 적용하여 $ m_{\text{max}} \sim 2.5-2.6\,M_\odot $를 추론한다. 이 결과는 GW170817에서 유도된 낮은 추정치를 도전하며, GW190814에서 관측된 약 $ 2.6\,M_\odot $의 성분이 중성자별일 가능성을 지지하며, 이는 이전에 생각했던 것보다 더 높은 본질적 최대 질량을 시사한다.

ABSTRACT

We have analyzed in this work the updated sample of neutron star masses derived from the study of a variety of 96 binary systems containing at least one neutron star using Bayesian methods. After updating the multimodality of the distributions found in previous works, we determined the maximum mass implied by the sample using a robust truncation technique, with the result $m_{max} \sim 2.5-2.6 \, M_{\odot}$. We have checked that this mass is actually consistent by generating synthetic data and employing a Posterior Predictive Check. A comparison with seven published $m_{max}$ values inferred from the remnant of the NS-NS merger GW170817 was performed and the tension between the latter and the obtained $m_{max}$ value quantified. Finally, we performed a Local Outlier Factor test and verified that the result for $m_{max}$ encompasses the highest individual mass determinations with the possible exception of PSR J1748-2021B. The conclusion is that the whole distribution already points toward a high value of $m_{max}$, while several lower values derived from the NS-NS merger event are disfavored and incompatible with the higher binary system masses. A large $m_{max}$ naturally accommodates the lower mass component of the event GW190814 as a neutron star.

연구 동기 및 목표

• 최소한 하나의 중성자별을 포함하는 96개의 이중성계에서의 최신 관측 데이터를 활용하여 중성자별의 최대 질량을 재평가한다.
• 이중성계에서 유도된 중성자별 질량 추정치와 GW170817과 같은 중력파 사건에서 유추된 추정치 사이의 갈등을 해결한다.
• 로컬 아웃라이어 요인(Locally Outlier Factor, LOF)과 같은 이상치 탐지 기법을 사용하여 고질량 중성자별 후보의 통계적 탄력성을 시험한다.
• 관측된 질량 분포가 이전에 보고된 바보다 더 높은 $ m_{\text{max}} $를 지지하는지 평가한다. 특히 GW190814 사건을 고려할 때 그러한 경향이 있는지 살펴본다.

제안 방법

• 96개의 이중성계에서 중성자별 질량의 사후 분포를 모델링하기 위해 마르코프 체인 몬테카를로(MCMC) 샘플링을 활용한 베이지안 추론을 적용한다.
• 질량 분포의 상한을 추정하기 위해 강력한 절단 기법을 사용하며, $ m_{\text{max}} $는 추론된 최대 질량으로 정의한다.
• 관측된 데이터와의 일관성을 검증하기 위해 5,000개의 시뮬레이션 데이터셋을 생성하여 사후 예측 검증(Posterior Predictive Check, PPC)을 수행한다.
• 질량 분포에서 잠재적 이상치를 탐지하기 위해 로컬 아웃라이어 요인(Locally Outlier Factor, LOF) 알고리즘을 활용하며, 극단적 질량 측정치가 주 분포에서 얼마나 떨어져 있는지 평가한다.
• 모델 적합도를 정량화하기 위해 임계값을 초과하는 시뮬레이션 질량 수를 테스트 통계량 $ T $로 정의하고, 베이지안 p-값을 계산한다.
• 마테매티카를 사용하여 표 2의 주변 사후 분포에서 새로운 데이터를 시뮬레이션함으로써 추론된 $ m_{\text{max}} $의 통계적 검증을 가능하게 한다.

실험 결과

연구 질문

• RQ1 최소한 하나의 중성자별을 포함하는 96개의 이중성계를 종합적으로 분석한 결과, 중성자별의 추론된 최대 질량은 얼마인가?
• RQ2 이중성계에서 추론된 $ m_{\text{max}} $는 GW170817 중성자별-중성자별 병합 사건에서 유도된 값과 어떻게 비교되는가?
• RQ3 관측된 고질량 중성자별 후보 PSR J1748-2021B는 주 분포와 통계적으로 일치하는가, 아니면 이상치인가?
• RQ4 추론된 $ m_{\text{max}} \sim 2.5-2.6\,M_\odot $는 GW190814 사건에서 관측된 약 $ 2.6\,M_\odot $ 성분이 중성자별일 수 있도록 수용 가능한가?
• RQ5 사후 예측 검증을 통해 검증된 바에 따르면, 더 높은 $ m_{\text{max}} $ 값은 통계적 불일치에 얼마나 강인한가?

주요 결과

• 분석 결과, 중성자별 최대 질량은 $ m_{\text{max}} \sim 2.5-2.6\,M_\odot $로 추론되었으며, 이는 GW170817에서 유도된 일부 이전 추정치보다 뚜렷이 높다.
• 사후 예측 검증은 추론된 $ m_{\text{max}} $가 관측된 데이터와 통계적으로 일관됨을 확인하여 그 신뢰성을 지지한다.
• $ m_{\text{max}} \sim 2.5-2.6\,M_\odot $ 값은 GW170817에서 유도된 여러 낮은 $ m_{\text{max}} $ 추정치와 불일치하며, 두 접근법 사이의 통계적 갈등을 시사한다.
• 로컬 아웃라이어 요인(Locally Outlier Factor, LOF) 테스트 결과, 가장 높은 개별 질량 측정치인 PSR J1748-2021B는 이상치가 아니며 추론된 분포와 일치함을 보였다.
• 더 높은 $ m_{\text{max}} $ 값은 GW190814 사건에서 관측된 약 $ 2.6\,M_\odot $ 성분이 중성자별일 수 있도록 자연스럽게 수용하며, 질량 간극에 대한 오랫동안 지속된 애매함을 해소한다.
• 이 연구는 전체 관측 질량 분포가 높은 $ m_{\text{max}} $를 향해 지시하고 있으며, 이는 이중성계에서 더 높은 질량을 가지는 것이 관측된 바에 비해 GW170817 기반의 낮은 값들이 통계적으로 불리하다고 결론내린다.