균일 회전 중성자별·쿼크별의 케플러 주파수와 경험식 검증

초록

본 연구는 다중 영역 스펙트럴 방법을 이용한 2차원 정밀 계산으로, 실질적인 핵물질 방정식(EOS)과 크러스트를 포함한 쿼크 별을 대상으로 케플러(질량 탈락) 회전 한계 구성을 구한다. Lattimer‑Prakash가 제시한 경험식 f_K = C (M/M_⊙)^½ (R/10 km)^‑3⁄2의 유효성을 검증했으며, 중성자별에서는 C≈1.08 kHz, 쿼크별(최대 크러스트)에서는 C≈1.15 kHz가 얻어졌다. 또한, 케플러 구성의 적도 반지름 R_K는 비회전 별의 반지름 R에 대해 R_K≈1.44 R이라는 일정한 비율을 보이며, 이는 고전적인 로셰 모델의 a=1.5보다 작다.

상세 분석

이 논문은 강체 회전하는 중성자별과 쿼크별의 질량 탈락(케플러) 한계 주파수를 정밀하게 계산하고, 기존에 제안된 경험식의 정확성을 체계적으로 검증한다. 저자들은 다중 영역 스펙트럴 방법을 기반으로 하는 2‑D 코드(LORENE 계열)를 사용해, 중성자별에 대해 10여 개의 현실적인 핵물질 방정식(EOS)과, 쿼크별에 대해 MIT bag model, NJL model 등 다양한 EOS를 적용하였다. 특히, 쿼크별 표면에 얇은 핵 크러스트(최대 질량 ≈10⁻⁴ M_⊙)를 포함시켜 실제 천체와의 비교 가능성을 높였다.

계산은 각 EOS마다 비회전(정지) 구성을 먼저 구하고, 동일 질량 M에 대해 회전 속도를 점진적으로 증가시켜 원심력에 의해 적도에서 물질이 탈락하는 순간을 찾는다. 이때 얻어지는 케플러 주파수 f_K와 반지름 R_K를, 비회전 별의 반지름 R(M)와 질량 M을 이용한 경험식
f_K(M)=C (M/M_⊙)^½ (R/10 km)^‑3⁄2
에 대입해 C 값을 역산한다. 결과는 0.5 M_⊙ < M < 0.9 M_max,stat 구간에서 2~3 % 이내의 오차로 경험식이 성립함을 보여준다. 중성자별의 경우 C_NS≈1.08 kHz, 쿼크별은 크러스트 유무에 크게 좌우되지 않는 C_QS≈1.15 kHz가 도출되었으며, 이는 로셰 모델이 예측한 C_Roche=1.00 kHz보다 현저히 큰 값이다.

또한, R_K와 R(M) 사이의 비례 관계 R_K≈a R(M)도 확인되었다. 여기서 a_NS≈a_QS≈1.44이며, 이는 로셰 모델의 a_Roche=1.5보다 작다. 이 비율은 회전으로 인한 원심 팽창이 EOS에 크게 의존하지 않으며, 질량 구간이 위와 같이 제한될 때 일정하게 유지된다.

이러한 결과는 고밀도 물질의 강체 회전 동역학을 이해하는 데 중요한 기준을 제공한다. 특히, 관측 가능한 맥스웰 주파수(예: 716 Hz 펄서)와 비교했을 때, 대부분의 실제 별은 케플러 한계보다 훨씬 낮은 회전 속도를 가지고 있음을 확인한다. 또한, 쿼크별에 대한 C 값이 중성자별보다 크게 나오므로, 만약 초고속 펄서가 관측된다면 EOS 구분에 활용될 가능성이 있다.