극비대칭 중성자별 물질의 상태방정식 액체방울 모델과 메존 폴트로프

초록

본 논문은 자유 전자·중성자 가스의 열역학과 액체방울 모델을 결합해 비대칭 중성자별 물질의 상태방정식을 구축하고, 고밀도에서는 양자 하드로다이내믹스 기반 메존 평균장(polytrope)으로 짧은 거리 반발을 추가한다. 얻어진 EOS를 토울먼‑오펜하이머‑볼크 방정식에 적용해 질량‑반경, 압축도, 표면 적색편이 등을 계산했으며, 관측된 X‑선 타이밍 및 중력파 데이터와 일치함을 보였다.

상세 분석

이 연구는 중성자별 내부 물질을 네 개의 밀도 구간으로 나누어 각각의 물리적 메커니즘을 명시적으로 모델링한다. 저밀도 구간(ρ≲10⁸ g cm⁻³)에서는 전자가 비상대론적 가스에서 초상대론적 퇴화 가스로 전이하며, 전자와 양성자·중성자 사이의 β‑평형 µ_p+µ_e=µ_n을 적용한다. 액체방울 모델(LDM)을 이용해 핵의 부피·표면·쿨롱·대칭 에너지를 파라미터화하고, 최소화 조건을 통해 핵질량수 A와 양성자 비율 Y_p를 결정한다. 이 과정에서 전자 드립(ρ≈4×10¹¹ g cm⁻³) 이후 자유 중성자 비율 Y_n이 도입되어, 핵 내부와 외부의 자유 중성자 기여를 구분한다.

고밀도 구간(ρ≳5×10¹³ g cm⁻³)에서는 핵 클러스터가 겹쳐 ‘핵 파스타’ 상태가 되며, 전통적인 표면·쿨롱 항을 무시한다. 여기서는 양성자 비율 Y_p와 자유 중성자 비율 Y_n을 β‑평형과 에너지 최소화식으로 연결하고, 핵 포화 밀도 ρ₀ 근처에서는 Y_n=1−ρ₀/ρ라는 간단한 보간식을 도입해 압력에 기여하지 않는 물질을 제외한다.

이후 자유 페르미온 압력만으로는 중성자별의 기계적 안정성을 보장할 수 없으므로, 양자 하드로다이내믹스(QHD) 평균장 모델을 차용해 σ, ω, ρ 메존 교환에 의한 짧은 거리 반발을 폴트로프 형태(P=K (Y_n ρ/ρ₀)^γ)로 표현한다. 메존 평균장은 스칼라 장 σ가 유효 중성자 질량 m*_n=m_n−g_σ σ를 감소시키고, 벡터 장 ω와 ρ가 압력에 양의 기여를 하도록 설계된다. 비선형 스칼라 자기상호작용을 최소화하고, 대신 밀도 의존적인 유효 질량을 도입해 파라미터 수를 줄이면서도 핵 포화 특성을 재현한다.

구축된 EOS는 연속성과 인과성을 보장하도록 각 구간 사이에 매끄러운 연결을 수행한다. 이를 토울먼‑오펜하이머‑볼크(TOV) 방정식에 입력하면 질량‑반경 곡선이 관측된 2 M⊙ 이상의 중성자별과 10–14 km 범위의 반경을 동시에 만족한다. 또한 표면 적색편이 z= (1−2GM/Rc²)^−½−1와 압축도 C=GM/Rc²도 최신 X‑ray 타이밍 및 중력파 파라미터와 일치한다. 모델의 단순성에도 불구하고, 메존 폴트로프의 강성도(γ)와 상수 K를 조절함으로써 EOS의 경도를 손쉽게 탐색할 수 있어, 보다 복잡한 핵 물리(예: 하이퍼온, 탈콘드)나 초밀도 상전이를 포함한 확장 연구의 베이스라인으로 활용 가능하다.