하이브리드 별의 스핀다운 압축에 의한 냉각 메커니즘

초록

스핀다운으로 중심 밀도가 증가하면서 변형된 우라카, 핵·쿼크 직접 우라카 등 다양한 중성미자 방출 과정이 순차적으로 활성화된다. 저자들은 각 과정의 방출률을 계산하고 회전하는 하이브리드 별의 냉각 곡선을 시뮬레이션한다. 결과는 자기장에 따라 직접 우라카가 발생 시 냉각이 크게 달라지며, 회전 모델이 정적 전이 모델보다 온도가 높고 복잡한 행동을 보인다는 것을 보여준다. 일부 회전 하이브리드 별의 예측 온도는 관측된 펄사 데이터와 일치한다.

상세 분석

본 논문은 하이브리드 별(핵 물질과 자유 쿼크 코어를 동시에 갖는 중성자 별)의 열진화에 스핀다운에 의한 구조 변화를 결합한 새로운 모델을 제시한다. 스핀다운은 주로 자기쌍극자 복사에 의해 발생하며, 각속도 Ω의 감소율은 (\dot Ω = -K B^{2} Ω^{3}) (K는 상수, B는 표면 자기장) 로 기술된다. Ω가 감소함에 따라 원심력 감소가 일어나 별의 중심 밀도가 점진적으로 상승한다. 이 밀도 상승은 핵 물질과 쿼크 물질의 화학적 포텐셜 차이를 변화시켜, 특정 임계 밀도에서 수정된 우라카(modified Urca) → 핵 직접 우라카(direct Urca) → 쿼크 직접 우라카(quark direct Urca) 순으로 중성미자 방출 메커니즘을 전환시킨다.

각 반응의 중성미자 방출률(ε)은 온도 T에 대한 지수적 의존성을 갖는다. 수정된 우라카는 (\epsilon_{MU} \propto T^{8}), 핵 직접 우라카는 (\epsilon_{DU}^{N} \propto T^{6}), 쿼크 직접 우라카는 (\epsilon_{DU}^{Q} \propto T^{6}) 로 표현된다. 따라서 직접 우라카가 활성화되면 냉각 속도가 급격히 가속된다. 저자들은 이러한 ε 값을 최신 핵·쿼크 상태 방정식(EOS)과 전이 임계 조건을 이용해 계산하고, 열전달 방정식
(C(T) \frac{dT}{dt} = -L_{\nu}(T,ρ(t)) - L_{\gamma}(T) + H_{comp}(t))
에 시간 의존적인 밀도 ρ(t)를 삽입해 수치 적분하였다. 여기서 (H_{comp})는 압축에 의한 내부 가열(압축 가열) 항이며, 스핀다운 속도가 빠를수록 이 항이 크게 기여한다.

시뮬레이션 결과는 두 가지 모델을 비교한다. 첫 번째는 ‘전이 정적 모델(transitional static model)’로, 스핀다운에 따라 별이 새로운 정적 평형 상태로 즉시 전이한다고 가정한다. 두 번째는 ‘회전 동적 모델(rotational dynamic model)’로, 구조 변화와 온도 변화를 동시에 추적한다. 회전 모델에서는 직접 우라카가 스핀다운 과정 중에 급격히 등장할 경우, 압축 가열과 결합해 초기 온도가 상승하고, 이후에도 높은 온도 유지가 가능함을 보인다. 특히 자기장이 (B \gtrsim 10^{13}) G인 경우, 스핀다운 속도가 빨라 직접 우라카가 일찍 발현하므로 냉각 곡선이 크게 변형된다.

이러한 결과는 관측된 젊은 펄사(예: Vela, Geminga)의 표면 온도와 비교했을 때, 전이 정적 모델이 예측하는 온도보다 회전 모델이 더 일치함을 보여준다. 즉, 스핀다운에 의한 압축 효과와 직접 우라카의 시기적 발현을 고려해야 하이브리드 별의 실제 냉각 경로를 정확히 설명할 수 있다.