NICER와 GW170817 관측이 제시하는 중성자별 압축 시나리오

초록

NICER와 GW170817의 질량·반지름·조석 변형률 데이터를 베이지안 분석에 적용해 속도-음속 파라미터화된 세 종류의 무모델 방정식(EoS) 군집을 제약하였다. 관측은 저밀도 영역을 강하게 제한하고, 중성자별이 중간 밀도에서는 부드럽고 고밀도에서는 급격히 강직해지는 구조를 선호한다는 결과를 보여준다. 이는 최대 질량을 내는 EoS보다 최대 압축성을 내는 EoS와 유사한 형태이며, 현재 데이터가 “압축된” 별을 더 선호한다는 해석을 가능하게 한다.

상세 분석

본 연구는 중성자별 내부 물질의 상태 방정식(EoS)을 모델에 의존하지 않고, 속도-음속(c_s^2)과 화학 퍼텐셜(µ)이라는 두 물리량을 조각별 선형 보간으로 파라미터화하는 방법을 채택하였다. 저밀도 영역(0.5–1.1 n₀)은 기존의 BPS·폴리트로프와 χEFT 기반의 압력 범위에 맞추어 고정하고, 그 위의 고밀도 구간은 네 개의 자유 구간으로 나누어 각각 c_s^2와 µ를 무작위로 샘플링한다. 이렇게 구성된 10개의 자유 파라미터는 베이지안 프레임워크에서 UltraNest 중첩 샘플링을 통해 약 5 × 10³개의 사후 샘플을 얻었다.

EoS는 크게 세 가지 형태—단조 증가(monotonic), 비단조(non‑monotonic), 불연속(discontinuous)—로 구분되었으며, 이는 각각 속도‑음속이 중심에서 감소, 외부에서 피크, 혹은 급격한 점프를 보이는지를 기준으로 한다. 각 형태에 대해 최대 질량(M_TOV)과 최대 압축성(최대 M/R) 구성을 계산하고, 관측 제약을 단계별로 적용했다. NICER 3개 펄서(PSR J0030+0451, J0740+6620, J0437‑4715)와 최신 PSR J0614‑3329, 그리고 GW170817의 조석 변형률 데이터를 독립적인 likelihood로 구성해 결합하였다.

결과적으로 모든 EoS 군집에서 가장 낮은 밀도 구간의 속도‑음속(c_s^2₁)은 약 0.04±0.01로 강하게 수축되었으며, 이는 χEFT와의 일치를 의미한다. 중간 밀도 구간(c_s^2₂, c_s^2₃)은 관측에 따라 크게 달라졌다. NICER와 GW170817를 모두 포함하면 c_s^2₂는 0.4 ~ 0.5 수준으로 낮아져 “부드러운” 구간을 형성하고, 고밀도 구간(c_s^2₃)은 0.6 ~ 0.7까지 급격히 상승한다. 이는 압축성을 극대화하는 EoS와 구조적으로 동일하며, 최대 질량을 내는 EoS가 보이는 초기 강직화와는 대조적이다.

불연속형 EoS에서는 밀도 점프(Δn)가 0 ~ 0.8 fm⁻³ 범위 내에서 작은 값으로 제한되었으며, 큰 점프는 관측과 충돌한다는 점이 확인되었다. 또한, 최대 질량을 달성하는 EoS는 중간 밀도에서 비교적 높은 c_s^2를 유지해 압력 지지를 제공하지만, 관측이 선호하는 “최대 압축성” EoS는 중간 밀도에서 부드럽고 고밀도에서 급격히 강직해지는 형태를 보인다.

이러한 결과는 현재 관측이 중성자별의 질량 상한보다는 반지름·조석 변형률을 통해 별의 압축성을 더 강하게 제약하고 있음을 시사한다. 즉, 중성자별이 중력에 의해 “압축된” 구조를 갖는 것이 데이터에 의해 자연스럽게 선택되고 있다. 이는 고밀도 QCD 상태(예: 강한 상전이)와 직접 연결되지는 않지만, 중성자별 내부 물질이 고밀도에서 급격히 강직해지는 보편적 패턴을 암시한다.