벨라형 펄서 글리치의 각운동량 전달 실현 모델

초록

펄서 Vela와 유사한 글리치를 설명하기 위해, 구형 중성자별을 n=1 폴리트로프 밀도 프로파일과 밀도 의존 핀ning 힘을 적용한 물리적으로 타당한 초유체 소용돌이 모델을 구축하였다. 최대 핀ning 힘 ≈ 10¹⁵ dyn cm⁻¹ 로 저장된 초유체 각운동량 ≈ 10⁴⁰ erg·s 가 약 3 년 주기로 방출되며, 전체 소용돌이의 10 % 이하만이 글리치 순간에 크러스트와 결합한다는 가정 하에 관측된 ∆Ω와 ∆Ω̇ 값을 재현한다.

상세 분석

이 논문은 Vela형 펄서 글리치를 초유체 소용돌이(vortex) 메커니즘으로 직접 정량화한 최초의 연구로, 기존의 원통형·균일밀도 가정에서 벗어나 구형 중성자별을 n=1 폴리트로프(압력 ∝ 밀도²) 모델로 설정하였다. 별의 중심밀도 λ와 차원화 반경 ξ = r/Rₛ 를 이용해 ρ(ξ)=λ·sin(πξ)/(πξ) 로 표현하고, 핵융합이 시작되는 핵심 반경 R_c 를 ρ_c = 0.6 ρ₀ 로 정의하였다. 이 구조는 내부 크러스트(ξ > x_c)와 핵심(ξ < x_c) 영역을 명확히 구분한다.

핀이 힘 f_pin(ρ)는 실험·이론적 계산에 기반해 최대값 f_m≈10¹⁵ dyn cm⁻¹ 를 갖고, 핵심 경계와 중성자 드립(ρ_d≈0.0015 ρ₀)에서 0 으로 감소한다. 핀ning 힘의 밀도 의존성은 그림 1(우)과 유사한 형태로, 초유체 짝짓기 갭이 최대가 되는 ρ≈0.2 ρ₀ 부근에서 피크를 보인다.

소용돌이 라인에 작용하는 마그누스 힘과 핀ning 힘을 축대칭으로 적분하면, 임계 라그 ω_cr(x) = (∫dz f_pin)/(κ x ∫dz ρ) 로 얻어진다. 여기서 κ=πħ/m_N 은 양자화된 순환량이다. 결과적으로 ω_cr(x)는 크러스트 영역에서 급격히 상승해 x_m≈1−u_m(≈0.8) 부근에 피크 ω_max를 형성하고, 핵심 영역에서는 ω_min≈10⁻² ω_max 수준의 거의 일정한 값을 가진다.

핀이 강하게 작용하는 크러스트에서는 소용돌이가 연속적으로 고정돼 라그가 ω_max에 도달하면 대량 탈핀ning이 일어나며, 이는 글리치의 급격한 각속도 상승 ∆Ω_gl을 초래한다. 반면 핵심에서는 소용돌이가 양쪽 끝만 핀ning되고 중간은 자유이므로, 실질적으로 약한 드래그(weak‑drag) 한계가 적용된다. 따라서 핵심 소용돌이는 짧은 시간(τ_c≈10⁰–10¹ s) 내에 지속적으로 탈핀·재핀을 반복하며, 평균 라그는 |Ω̇| τ_c 수준으로 유지된다.

글리치 직전, 별이 감속하면서 라그 ω(t)=|Ω̇_∞| t 가 증가하고, 이에 대응하는 임계 반경 x(t) 가 외부로 이동한다. 라그가 ω_max에 도달하면, x(t)=x_m 에서 축적된 소용돌이 시트가 동시에 방출돼 전체 각운동량 ∆L_gl=I_c ∆Ω_gl≈10⁴⁰ erg·s 를 전달한다. 이 값은 기존 모델이 요구한 ∆L≈10⁴¹ erg·s 보다 한 차수 낮지만, 전체 소용돌이 중 <10 %만이 글리치 순간에 크러스트와 결합한다는 가정 하에 관측된 ∆Ω_gl≈10⁻⁶ Ω와 ∆Ω̇_gl≈10⁻² Ω̇ 를 충분히 재현한다.

핵심‑크러스트 사이에 정상 물질 층이 없다는 최신 미시적 계산(Zhou et al. 2004)을 반영함으로써, 소용돌이 라인이 연속적으로 존재한다는 물리적 일관성을 유지한다. 또한, 초전도성 양성자에 대한 약한 마찰 가정을 채택했지만, 강한 마찰(강형 초전도) 경우에도 ∆L_gl와 관측값 사이의 차이는 비슷한 수준으로 유지된다. 최종적으로, 이 모델은 Vela 글리치의 주기(≈3 년), 각운동량 규모, 그리고 포스트 글리치 회복에 필요한 다중 지수 감쇠를 모두 일관된 물리적 메커니즘으로 설명한다.