초전도 중성자 별의 비방사 진동 모드와 온도 의존성

초록

본 논문은 비회전 초전도 중성자 별의 비방사(non‑radial) 진동을 유한 온도에서 자체적으로 계산한다. APR 방정식과 밀도 의존적인 중성자·양성자 임계 온도 모델을 적용해, 중심과 외부 코어는 비초전도, 중간 구역만 초전도인 3‑층 구조를 포함한다. 초전도 모드의 고유진동수는 온도 변화에 매우 민감하며, 별이 냉각함에 따라 진동 스펙트럼이 급격히 변한다는 결과를 제시한다.

상세 분석

이 연구는 기존에 0 K 가정에 머물렀던 초전도 중성자 별의 진동 이론을 확장하여, 실제 별 내부 온도가 10⁸–10⁹ K 수준일 때의 동역학을 정밀히 다룬다. 핵심은 일반 상대성 하에서 초전도 흐름을 기술하는 4‑벡터 wᵘ(k)와 엔트레인먼트 행렬 Yᵢₖ를 이용해, 정상 물질과 초전도 물질 사이의 결합 파라미터 s가 |s|≲0.05임을 확인하고, s≈0인 제한에서 두 계를 완전히 분리한다는 전제 하에 초전도 전용 진동 방정식(식 5)을 도출한다. 이 방정식은 온도 의존적인 파라미터 y(T/T_cⁿ, T/T_cᵖ)를 통해 초전도 영역의 크기와 초전도 사운드 속도 v_sf를 연결한다.

핵심 물리 입력으로는 APR 방정식과, 중성자 T_cⁿ(r)와 양성자 T_cᵖ(r)의 밀도 의존 프로파일을 채택하였다. 특히 T_cⁿ는 핵심에서 약 5×10⁸ K의 최대값을 갖고, 별이 냉각될수록 중심에서부터 외부까지 초전도 영역이 확장되는 3‑층 구조(중심 비초전도‑중간 초전도‑외부 비초전도)를 만든다. 온도 T_∞가 일정하게 유지된다고 가정하면, 적색편이된 임계 온도 T_c^∞(r)로 초전도 영역을 정의할 수 있다.

수치 해석에서는 다중극성 l=0–3에 대해 방사적 노드 수 n=0–4인 고유모드를 계산하였다. 결과는 다음과 같다. (1) n이 클수록 고유진동수 ω_ln이 증가한다. (2) 온도가 T_cⁿ보다 높을 때는 초전도 영역이 없으므로 초전도 모드가 사라진다. (3) T_∞가 T_cⁿ의 최대값에 접근하면 초전도 영역이 급격히 확대하면서, ω_ln은 두 상반된 효과—영역 확대에 의한 진동수 감소와 v_sf 상승에 의한 진동수 증가—가 상쇄되어 거의 온도 독립적인 값을 보인다. 특히 n=0 모드(가장 낮은 방사적 노드)는 T_∞가 T_cⁿ에 가까워질수록 ω_l0이 선형적으로 0에 수렴한다. 이는 식 ω_l0∝l(l+1) v_sf / r_maxⁿ · (1−T_∞/T_cⁿ) 형태로 해석된다.

또한, 고온(≈4×10⁸ K)에서는 초전도 영역이 얇은 구형 쉘 형태이므로, 방정식 5의 주요 항은 h⁻¹ 항만 남아 레전드르 다항식 P_n에 비례하는 고유함수를 갖는다. 이 경우 ω_ln∝n(n+1) 형태의 간단한 스펙트럼이 나타난다.

시간적 진화 측면에서는 별이 냉각함에 따라 초전도 영역이 중심에서부터 전체 코어까지 확장되는 과정을 시뮬레이션하였다. 결과는 10⁴ 년 이내에 초전도 모드의 진동수가 수백 Hz 수준에서 수 kHz 수준으로 급변함을 보여준다. 이는 관측 가능한 X‑ray 진동수 변동이나 중성자 별의 그라비톤 방출과 연결될 가능성을 시사한다.

전반적으로 이 논문은 (i) 초전도 물질의 엔트레인먼트와 온도 의존성을 정확히 포함한 일반 상대성 진동 방정식을 제시하고, (ii) 현실적인 핵 물성 모델을 적용해 온도에 민감한 새로운 초전도 비방사 모드 스펙트럼을 예측함으로써, 중성자 별 내부 초전도 상태를 천체물리학적 관측으로 탐색할 수 있는 이론적 기반을 제공한다는 점에서 큰 의미를 가진다.