밀도 의존 초유체 에너지 갭을 고려한 회전화학 가열

초록

회전 속도가 감소하면서 원심력이 약해져 별이 수축하고, 압력 상승으로 베타 평형이 깨진다. 이때 초유체 에너지 갭이 밀도에 따라 변하는 모델들을 적용해 회전화학 가열을 계산하고, 예측된 표면 온도를 근접 밀리초 펄서 J0437‑4715의 관측값과 비교한다.

상세 분석

본 논문은 회전하는 중성자 별이 각운동량을 잃을 때 발생하는 회전화학 가열(rotochemical heating) 현상을, 핵물질이 초유체(superfluid) 상태에 있을 경우에 한정하여 재검토한다. 기존 연구(Fernández & Reisenegger, 2005)는 비초유체 물질을, Petrovich & Reisenegger(2010)는 초유체 물질을 가정했으나, 두 경우 모두 초유체 에너지 갭을 일정한 값으로 고정하였다. 그러나 핵물리학 이론에 따르면 중성자 별 내부의 초유체 에너지 갭은 밀도에 따라 크게 변동한다. 즉, 중심부에서는 큰 갭이, 외부 연성층에서는 작은 갭이 존재한다는 것이 최신 핵 상호작용 모델(Nambu–Jona-Lasinio, BCS 등)에서 제시된다. 이러한 공간적 변화를 무시하면 반응률, 즉 뮤온 및 전자 캡처와 베타 붕괴 과정에서 방출되는 열량을 과소·과대 평가하게 된다.

논문은 먼저 여러 제안된 갭 모델(예: “A”, “B”, “C” 유형) 중 3~4가지를 선택하고, 각 모델에 대해 밀도 프로파일을 별 구조 방정식(TOV)으로부터 얻는다. 그런 다음, 회전감소에 따른 압력 변화 ΔP와 그에 따른 베타 불균형 η(=μ_n−μ_p−μ_e) 를 계산한다. η는 초유체 갭 Δ(ρ) 에 의해 억제되며, 반응률 λ∝exp(−Δ/kT) 형태로 감소한다. 따라서 갭이 큰 영역에서는 η가 크게 축적되고, 열 생산이 지연된다. 반면, 갭이 작은 외부에서는 반응이 빠르게 진행되어 η를 완화시키고, 열이 더 빨리 방출된다.

시간 진화 방정식 dη/dt = −(∂η/∂t)_comp − λ(η,T) 를 수치적으로 적분해, 별이 수백 Myr에서 수 Gyr까지 도달하는 quasi‑steady state 를 찾는다. 여기서 quasi‑steady state는 압축에 의해 생성된 η와 초유체 억제 하의 반응에 의해 소멸되는 η가 균형을 이루는 상태이며, 이때 표면 온도 T_s는 내부 온도 T와 열 전도도 κ에 의해 결정된다.

핵심 결과는 다음과 같다. (1) 밀도 의존 갭을 도입하면 quasi‑steady state에 도달하는 시간과 최종 η 값이 크게 변한다. (2) 갭이 큰 중심부 모델에서는 내부 온도가 상대적으로 높게 유지되어, 표면 온도 T_s ≈ (0.5–1)×10^6 K 수준을 보인다. (3) 반면, 갭이 작은 모델은 η가 빠르게 소멸해 T_s가 ≈ (0.2–0.4)×10^6 K 로 낮아진다. (4) 관측된 가장 가까운 밀리초 펄서 J0437‑4715의 표면 온도(T_s ≈ 1.5×10^6 K)와 비교했을 때, 중심부에 큰 갭을 갖는 “A형” 모델이 가장 일치한다는 결론을 내렸다.

이러한 분석은 초유체 갭의 밀도 의존성을 고려하지 않은 이전 연구와 비교해, 회전화학 가열이 실제 중성자 별에서 어떻게 작동하는지를 보다 정밀하게 예측한다는 점에서 의미가 크다. 또한, 관측 가능한 X‑ray 스펙트럼을 통해 초유체 갭 모델을 검증할 수 있는 새로운 방법론을 제시한다는 점에서도 학문적 가치를 지닌다.