[논문 리뷰] On the Sound Speed in Neutron Stars

Sinan Altiparmak, Christian Ecker외 1명|arXiv (Cornell University)|2022. 03. 28.

Experimental and Theoretical Physics Studies인용 수 10

한 줄 요약

이 연구는 중성자별에 대한 1,000만 개 이상의 연속적인 음속 상태방정식(EOS)을 구축하며, 핵물리이론, 양자색역학(QCD)의 섭동 이론, 천문학적 제약 조건을 통합한다. 연구 결과, 음속 제곱(c²ₛ < 1/3)이 하위동형성인 경우는 통계적으로 매우 드물게 발생하며(0.03%), 별 내부에서 c²ₛ > 1/3인 모델이 강하게 지배함을 시사한다. 95% 신뢰도 추정치로는 R₁.₄ = 12.42⁺⁰.⁵²₋₀.⁹⁹ km, R₂.₀ = 12.12⁺¹.¹¹₋₁.²³ km이며, GW170817 사건에 대한 타이드 변형도 Ũ₁.₁₈₆ = 485⁺²²⁵₋₂₁₁로 도출되었다.

ABSTRACT

Determining the sound speed $c_s$ in compact stars is an important open question with numerous implications on the behaviour of matter at large densities and hence on gravitational-wave emission from neutron stars. To this scope, we construct more than $10^7$ equations of state (EOSs) with continuous sound speed and build more than $10^8$ nonrotating stellar models consistent not only with nuclear theory and perturbative QCD, but also with astronomical observations. In this way, we find that EOSs with sub-conformal sound speeds, i.e. with $c^2_s &lt; 1/3$ within the stars, are possible in principle but very unlikely in practice, being only $0.03\%$ of our sample. Hence, it is natural to expect that $c^2_s &gt; 1/3$ somewhere in the stellar interior. Using our large sample, we obtain estimates at $95\%$ credibility of neutron-star radii for representative stars with $1.4$ and $2.0$ solar masses, $R_{1.4}=12.42^{+0.52}_{-0.99}\,{ m km}$, $R_{2.0}=12.12^{+1.11}_{-1.23}\,{ m km}$, and for the binary tidal deformability of the GW170817 event, $ ildeΛ_{1.186}=485^{+225}_{-211}$. Interestingly, our lower-bounds on the radii are in very good agreement with the prediction derived from very different arguments, namely, the threshold mass. Finally, we provide simple analytic expressions to determine the minimum and maximum values of $ ildeΛ$ as a function of the chirp mass.

연구 동기 및 목표

• 실제 물리적 및 관측 제약 조건 하에서 중성자별 내에서 하위동형성 음속(c²ₛ < 1/3)이 발생할 가능성을 규명하는 것.
• 질량이 2 M☉ 이상인 중성자별이 존재하는 바탕으로, 별 내부에서 음속이 어떤 밀도에서라도 1/3을 초과해야 하는지 평가하는 것.
• 큰 규모의 물리적으로 일관된 상태방정식 샘플을 기반으로 관측에 기반한 신뢰할 수 있는 중성자별 반경 및 타이드 변형도 추정치를 제공하는 것.
• 조각별 다항식 모델에서 관찰되는 불연속성 문제를 피하기 위해 음속 매개변수화 기반의 통계적으로 안정된 상태방정식 특성 추정 프레임워크를 개발하는 것.

제안 방법

• Annala 등(2020)의 매개변수화 방법을 활용해 10⁷개 이상의 연속적인 음속을 가진 상태방정식(EOS)을 구축하며, 열역학적 및 원인성 안정성을 확보한다.
• 핵물리이론(BPS 모델, n < 0.5n₀), 효과적 장 이론(n ≲ n₀), 섭동 QCD(n ≫ n₀)를 통합하여 물리적 일관성을 확보한다.
• 최대 7개의 밀도 세그먼트를 사용한 패치워크 접근법을 통해 상태방정식를 모델링하고, 비물리적 해를 제거하고 수렴성을 향상시키기 위해 필터링을 적용한다.
• 모든 상태방정식 샘플에 대해 비회전 중성자별의 별 모형을 계산하여 질량-반경 관계 및 타이드 변형도 관계를 유도한다.
• 베이지안 신뢰도 구간을 적용하여 반경 및 타이드 변형도에 대한 95% 신뢰도 범위를 추정하며, 통계적 변동을 최소화하기 위해 대규모 샘플을 활용한다.
• 음속의 확률밀도함수(PDF)를 분석하여 다양한 시나리오(단조증가, 비단조증가, 하위동형성, 초동형성)의 통계적 가능성을 규명한다.

실험 결과

연구 질문

• RQ1 실제 물리적 및 관측 제약 조건 하에서 중성자별 내 음속 제곱(c²ₛ < 1/3)이 발생할 통계적 가능성이 얼마인가?
• RQ2 질량이 2 M☉ 이상인 중성자별이 존재하는 것으로 볼 때, 별 내부에서 어떤 밀도에서라도 음속이 1/3을 초과해야 하는가?
• RQ3 GW170817 및 중성자별 질량에서의 관측 제약 조건이 중성자별 반경 및 타이드 변형도의 신뢰도 범위에 어떤 영향을 미치는가?
• RQ4 고밀도 QCD 영역의 경계 조건에 따라 음속 프로파일의 특성(예: 단조성, 피크 위치)이 어느 정도 영향을 받는가?

주요 결과

• 하위동형성 음속(c²ₛ < 1/3)은 분석한 10⁷개 이상의 상태방정식 중 단지 0.03%에서만 발생하여 자연계에서 매우 불가능한 것으로 나타났다.
• 1.4 M☉ 중성자별의 95% 신뢰도 반경은 R₁.₄ = 12.42⁺⁰.⁵²₋₀.⁹⁹ km이며, 2.0 M☉ 별의 경우 R₂.₀ = 12.12⁺¹.¹¹₋₁.²³ km이다.
• GW170817 이벤트에 대해 타이드 변형도의 95% 신뢰도 추정치는 Ũ₁.₁₈₆ = 485⁺²²⁵₋₂₁₁로, 현재 관측 제약 조건과 일치한다.
• 이 분석에서 도출된 중성자별 반경 하한선은 임계 질량 추론에 의해 유도된 예측과 뛰어난 일치를 보이며, 독립된 이론적 접근법 간의 강건성을 시사한다.
• 음속이 밀도에 따라 증가하는 단조증가 프로파일(c²ₛ 증가)은 1.99 M☉를 초과하는 별을 지지하기에는 부족하여, 이러한 모델은 관측된 질량이 큰 중성자별과 호환되지 않음을 시사한다.
• 중성자별 밀도 근처에서 음속 PDF의 피크는 동형값 자체가 아니라, 원인성적으로 섭동 QCD 근사에 접근해야 한다는 요구 조건에 의해 주로 유도된다.