Evolution of proto-neutron stars to pulsars, magnetars and central compact objects

기술 요약: 원시 중성자별에서 펄서, 마그네타 및 중심 컴팩트 천체(CCO)로의 진화

문제 제기
마그네타에서 관측되는 초강력 자기장($B \sim 10^{14}–10^{15}$ G)의 기원은 여전히 논쟁의 대상이다. "화석 자기장(fossil field)" 가설(전구체 별로부터의 자속 보존)은 일반적인 중성자별과 중심 컴팩트 천체(CCO)의 자기장은 설명할 수 있지만, 극단적으로 희귀한 전구체 조건을 가정하지 않고서는 마그네타의 자기장을 설명하는 데 어려움을 겪는다. 반대로, 마그네타를 위한 다이내모 가설은 원시 중성자별(PNS) 단계에서의 빠른 회전과 격렬한 대류를 필요로 하며 널리 받아들여지고 있으나, 마그네타의 추정 탄생율이 코어 붕괴 초신성 발생률과 맞먹는다는 점이 문제다. 이는 마그네타 형성이 극단적인 매개변수 요구 조건이 시사하는 것만큼 드문 일이 아님을 의미한다. 또한, 표준 자속 보존 계산에 따르면 철 핵의 붕괴 과정에서도 결과적인 자기장($B \sim 10^{10}$ G)은 일반적인 펄서 자기장($B \sim 10^{12}$ G)을 설명하기에 불충분하며, 이는 일반 펄서에게도 다이내모 메커니즘이 필요함을 시사한다. 핵심 문제는 다양한 종류의 중성자별의 높은 탄생율을 통일된 진화적 틀 안에서 어떻게 조화시킬 것인가 하는 점이다.

방법론
저자들은 백색 왜성 병합에 적용되었던 Wickramasinghe 등(2014)의 모델을 응용한, 축소된 1차원 전단 구동(shear-driven) $\alpha-\Omega$ 다이내모 모델을 사용한다. 이 모델은 원시 중성자별(PNS) 단계 동안 폴로이달(poloidal, $B_p$) 및 토로이달(toroidal, $B_\phi$) 자기장의 진화를 추적한다. 이 방법론의 주요 특징은 다음과 같다:

• 지배 방정식: 모델은 $B_p$와 $B_\phi$에 대한 결합된 미분 방정식을 푼다. 토로이달 자기장은 $\Omega$-효과(폴로이달 자기장에 작용하는 차등 회전)를 통해 생성되며, 폴로이달 자기장은 토로이달 자기장에 작용하는 $\alpha$-효과(난류 대류)를 통해 재생된다.
• PNS 진화: 정적인 모델과 달리, 이 프레임워크는 초기 반지름 $R_0 \sim 40$ km에서 최종 반지름 $R_\infty \approx 12$ km까지의 PNS의 역동적인 수축을 포함한다. 반지름 진화는 다양한 질량과 상태 방정식(EoS)에 대해 Camelio 등(2017)의 방식을 기반으로 매개변수화된다.
• 질량 및 각운동량: 모델은 중성미자 냉각에 따른 중력 질량의 감소와 각운동량 보존을 고려하며, 이는 PNS의 스핀업(spin-up)을 유도한다.
• 점성 및 전단 감쇠: 차등 회전의 감쇠(전단율 $q$)는 자기 회전 불안정성(MRI)에 의한 점성과 대류 점성의 조합을 통해 모델링된다. 전단율은 시간에 따라 감쇠하며, 결국 다이내모를 종료시킨다.
• 자속 손실 및 안정성: 모델은 자기 자속 손실을 고려하기 위한 재결합 항(초기값 20%로 설정)을 포함하며, 자기장 불안정성이 $B_p$와 $B_\phi$의 붕괴 시간 척도를 결정하는 안정성 기준(Braithwaite 2009)을 통합한다.
• 매개변수 공간 스캔: 저자들은 룬게-쿠타(Runge-Kutta) 방법을 사용하여 방정식을 수치적으로 적분하며, 초기 전단율($q_0$), 초기 질량($M_0$), 초기 반지름($R_0$), 초기 씨앗 자기장(seed fields), 그리고 최종 회전 주기($P_\infty$)를 포함한 초기 매개변수들을 체계적으로 변화시킨다.

주요 결과
시뮬레이션 결과, 최종 자기장 구성은 초기 회전율과 전단에 매우 민iente하게 반응하며, 이를 통해 단일 다이내모 메커니즘이 관측되는 전체 중성자별 유형의 스펙트럼을 생성할 수 있음을 보여준다:

1. 전단율 임계값: 토로이달 자기장의 포화 값($B_{\phi, \infty}$)은 초기 전단율 $q_0$에 강하게 의존한다. $q_0 \simeq 0.23$에서 임계값이 존재한다.

   • 임계값 미만 ($q_0 < 0.23$): $\alpha$-효과가 충분히 유발되지 않는다. 폴로이달 자기장은 자속 보존에 의해 결정된 값($B_p \approx 2.5 \times 10^{10}$ G) 근처에 머물며, 토로이달 자기장은 유의미하게 성장하지 않는다. 이 영역은 중심 컴팩트 천체(CCO)에 해당한다.
   • 임계값 초과: 다이내모가 활성화된다. 더 높은 전단율은 더 강한 토로이달 자기장을 생성한다.

2. 서로 다른 집단의 형성:

   • 마그네타: 가장 짧은 초기 주기(빠른 회전, 예: $P_\infty \sim 2.5$ ms)의 경우, 모델은 $B_p \simeq 10^{15}$ G 및 $B_\phi \simeq 10^{16}$ G의 마그네타급 강도를 생성한다.
   • 일반 펄서: 더 느린 회전자의 경우, 모델은 전형적인 펄서 자기장($B_p \simeq 10^{12}$ G, $B_\phi \simeq 10^{14}$ G)을 산출한다.
   • 저-B 마그네타/고-B 펄서: 중간 정도의 회전율은 저-B 마그네타 또는 고-B 펄서 범위의 자기장을 생성한다.

3. 자속 보존의 역할: 초기 폴로이달 자기장 강도는 주로 붕괴 중의 자속 보존에 의해 결정된다. 그러나 관측된 펄서 및 마그네타 수준으로 이 자기장들을 증폭시키기 위해서는 다이내모 과정(특히 $\alpha$-효과)이 필수적이다. 토로이달 자기장은 일반적으로 $\Omega$-효과에 의해 지배된다.

4. 중심 컴팩트 천체 (CCOs): 모델은 CCO가 느린 초기 회전($P_0 \sim 100$ ms)으로 인해 격렬한 다이내모 과정을 거치지 않는다고 제안한다. 그들의 폴로이달 자기장은 자속 보존에 의해 결정되지만, 오직 $\Omega$-효과에 의해서만 증폭된 토로이달 성분을 가질 수 있으며, 이는 작은 쌍극자 자기장을 가지면서도 잠재적으로 유의미한 내부 토로이달 성분을 갖는 이유를 설명한다.

의의 및 주장
본 논문은 마그네타의 높은 탄생율과 마그네타 형성에 흔히 요구되는 극단적인 조건 사이의 모순을 해결한다고 주장한다. 초기 회전율에 따라 결정되는 다이내모 과정이 중성자별의 다양한 자기장 강도를 자연스럽게 생성할 수 있음을 입증함으로써, 저자들은 다음과 같이 제안한다:

• 일반 펄서와 마그네타 모두 다이내모 단계를 거칠 가능성이 높으며, 차이점은 초기 회전에 의해 구동되는 이 과정의 효율성에 있다.
• "화석 자기장" 가설은 전형적인 펄서 자기장을 설명하기에 불충분이므로, 비마그네타 중성자에게도 다이내모가 필요하다.
• CCO는 다이내모가 억제되어 자속 보존에 의해 자기장이 결정되는 스펙트럼의 저회전 종단(low-rotation end)을 나타낸다.

저자들은 자신들의 1차원 모델이 단순화되었음에도 불구하고 회전, 전단, 부력에 대한 핵심적인 의존성을 포착하고 있으며, 더 복式한 3차원 시뮬레이션(예: Raynaud 등 2020)과 질적으로 일치한다고 강조한다. 결론적으로, $\alpha-\Omega$ 다이내모는 마그네타를 위해 희귀하고 극단적인 전구체 조건을 요구하지 않고도, 젊은 중성자별의 다양한 자기장을 이해하기 위한 통일된 프레임워크를 제공한다.