[논문 리뷰] Multi-messenger constraints on the neutron-star equation of state and the Hubble constant

Tim Dietrich, M. W. Coughlin외 6명|arXiv (Cornell University)|2020. 02. 26.

Gamma-ray bursts and supernovae인용 수 6

한 줄 요약

이 연구는 중성자별 상태방정식과 허블 상수를 제약하기 위해 GW170817, GW190425, AT2017gfo, GRB170817A 및 펄서 관측 데이터를 양성자 효과 이론 핵물리 계산과 융합한다. 중력파, 킬로노바 빛의 곡선, 질량 측정을 포함한 베이지안 순차 분석을 통해, 1.4 M⊙ 중성자별 반지름은 11.75+0.86−0.81 km이고, 허블 상수는 1σ 신뢰수준에서 66.2+4.4−4.2 km Mpc⁻¹ s⁻¹로 측정된다.

ABSTRACT

Observations of neutron-star mergers based on distinct messengers, including gravitational waves and electromagnetic signals, can be used to study the behavior of matter denser than an atomic nucleus, and to measure the expansion rate of the Universe described by the Hubble constant. We perform a joint analysis of the gravitational-wave signal GW170817 with its electromagnetic counterparts AT2017gfo and GRB170817A, and the gravitational-wave signal GW190425, both originating from neutron-star mergers. We combine these with previous measurements of pulsars using X-ray and radio observations, and nuclear-theory computations using chiral effective field theory to constrain the neutron-star equation of state. We find that the radius of a $1.4$ solar mass neutron star is $11.75^{+0.86}_{-0.81}\ m km$ at $90\%$ confidence and the Hubble constant is $66.2^{+4.4}_{-4.2}\ m km \,Mpc^{-1}\, s^{-1}$ at $1\sigma$ uncertainty.

연구 동기 및 목표

• 다양한 메시지 관측과 핵물리학 이론을 융합하여 중성자별 상태방정식(EOS)에 대한 제약을 향상시키기.
• 이중 중성자별 병합에서 발생한 중력파 및 전자기파 데이터를 이용해 허블 상수를 측정하기.
• 중력파, 킬로노바, 감마선 폭발, 펄서 질량-반지름 측정을 포함한 천체물리학 관측을 통합된 베이지안 프레임워크에 통합하기.
• 웨이브폼 모델, 분석 순서, 킬로노바 모델 가정에 대한 결과의 강인성 평가하기.
• 저밀도에서의 치르알 효과 이론 예측과 고밀도에서의 모델에 의존하지 않는 파arametric 확장 방식을 융합하여 인과성과 이론적 일관성을 확보하기.

제안 방법

• 핵 포화 밀도(nsat)의 1.5배 이하까지 치르알 효과 이론(EFT)을 사용해 5,000개의 중성자별 상태방정식(EOS) 모델을 구축한다.
• 1.5 nsat를 초과하는 영역에서는 인과성을 유지하기 위해 속도-소리 평면에서 모델에 의존하지 않는 파라메트릭 전개 방식을 사용해 EOS 모델을 연장한다.
• 순차적으로 제약 조건을 통합하기 위해 베이지안 추론을 적용한다: 펄서에서의 최대 중성자별 질량과 GW170817/AT2017gfo, NICER가 측정한 PSR J0030+0451의 질량-반지름 측정, IMRPhenomPv2_NRTidalv2를 사용한 GW170817 중력파 웨이브폼, AT2017gfo 킬로노바 빛의 곡선, 그리고 전자기 방출 비검출 여부를 고려한 GW190425의 중력파 데이터.
• GW170817과 GW190425에 대해 NRTidalv2 웨이브폼 모델을 사용하고, 결과를 SEOBNRv4T와 비교하여 웨이브폼 민감도를 테스트한다.
• 분석 단계의 순서를 변화시키고, AT2017gfo에 대해 다양한 킬로노바 모델(I–III)을 테스트함으로써 체계적 불확실성 전파를 수행한다.
• 후행 분포를 순차적 분석 단계에 걸쳐 전파하기 위해 커널 밀도 추정을 적용하고, 일관된 사전 분포 구축을 보장한다.

실험 결과

연구 질문

• RQ1 다양한 메시지 관측과 핵물리학 이론을 융합하여 1.4 태양질량 중성자별 반지름에 대해 가장 정밀한 제약 조건을 얻을 수 있는가?
• RQ2 AT2017gfo 킬로노바 빛의 곡선과 GW190425의 전자기 병합 부재 정보를 포함함으로써 중성자별 상태방정식에 대한 제약 조건은 어떻게 변화하는가?
• RQ3 다른 중력파 웨이브폼 모델(예: NRTidalv2 대비 SEOBNRv4T)이 유도된 중성자별 반지름과 허블 상수에 어떤 영향을 미치는가?
• RQ4 베이지안 분석에서 천체물리학적 및 이론적 제약 조건을 적용하는 순서에 따라 최종 제약 조건이 얼마나 강인한가?
• RQ5 GW170817, AT2017gfo, 펄서 질량 데이터와 핵물리학 EFT 제약 조건을 융합했을 때 허블 상수는 어떻게 측정되는가?

주요 결과

• 모든 다메시지어 데이터와 치르알 EFT 제약 조건을 통합한 결과, 1.4 M⊙ 중성자별 반지름은 90% 신뢰구간에서 11.75+0.86−0.81 km로 제약된다.
• 허블 상수는 GW170817, AT2017gfo, 펄서 질량 측정 데이터의 공동 분석을 통해 1σ 신뢰수준에서 66.2+4.4−4.2 km Mpc⁻¹ s⁻¹로 측정된다.
• GW170817 및 AT2017gfo 정보를 GW190425 분석에 통합함으로써 질량 가중 평균 휘어짐도 ˜Λ는 140.80+144.22−64.73에서 117.90+114.60−49.24로 감소하여 더 작은 반지름을 지지한다.
• 웨이브폼 모델 선택에 대해 결과가 강인하다: NRTidalv2 대신 SEOBNRv4T를 사용해도 반지름에 미치는 영향은 매우 작고, 불확도 범위 내에 국한된다.
• 분석 단계의 순서에 미치는 영향은 미미하며, NICER 데이터를 먼저 적용하는 대신 나중에 적용해도 후행 분포의 반지름 제약 조건은 안정하다.
• AT2017gfo에 대해 다양한 킬로노바 모델(I–III)을 사용했을 때 거리-기울기 측정치는 통계적 불확실성 내에서 일관되며, 이는 결과가 주로 통계적 효과에 의해 결정됨을 시사한다.