중성자별 토로이드 자기장 비축대칭 불안정 연구

초록

이 논문은 강한 토로이드 자기장을 가진 등온 중성자별이 비축대칭 섭동에 대해 얼마나 불안정한지를 3차원 일반 상대론적 MHD 시뮬레이션으로 조사한다. 비회전 및 강체 회전 모델 모두에서 파커·타일러·교환 불안정이 나타나며, 회전이 파커 불안정을 완전히 억제하지 못한다는 점을 확인한다. 불안정이 발달하면 난류가 지속적으로 유지되고, 자기장은 축대칭 경우와 달리 새로운 평형에 도달하지 못한다.

상세 분석

본 연구는 토로이드 자기장이 지배적인 중성자별이 실제 천체에서 가질 수 있는 안정성 한계를 정량적으로 제시한다. 저자들은 먼저 등온 폴리트로프 방정식을 이용해 토로이드 자기장이 내부에만 존재하도록 설계된 평형 모델을 구축하였다. 자기장 프로파일은 Bφ∝ρhα²γφφ^k·γ^(-½)φφ 형태이며, k=1인 경우만을 선택해 축대칭 교환 불안정이 억제된 상태를 가정하였다. 그러나 비축대칭 섭동을 허용한 3D GRMHD 시뮬레이션 결과, 파커 불안정이 가장 먼저 성장함을 확인했다. 파커 불안정은 자기장 선이 수직으로 상승하는 영역에서 부압이 감소하면서 발생하는 부력형 불안정으로, 알펜 시간 척도(τ_A≈R/v_A) 정도의 성장률을 보인다. 특히, 회전이 존재하더라도 회전축을 중심으로 각운동량 기울기가 음수인 영역이 형성되면 파커 불안정이 여전히 발달한다. 이는 기존 선형 이론이 급격한 회전이 모든 비축대칭 모드를 억제한다는 가정과 달리, 실제 별 내부에서 발생하는 펄스와 각운동량 재분포가 불안정을 촉진한다는 점을 시사한다.

시뮬레이션에서는 회전 모델에서도 타일러 불안정과 교환 불안정이 부분적으로 나타났으며, 특히 회전이 빠를수록 교환 불안정이 억제되는 대신 파커 불안정이 더욱 두드러졌다. 불안정이 전개된 뒤에는 자기장과 유체가 복잡한 난류 상태에 머물며, 자기 에너지는 지속적으로 감소하지만 새로운 정적 평형에 도달하지 못한다. 이는 축대칭 2D 시뮬레이션에서 관찰된 ‘자기장 재구성 후 안정화’와는 근본적으로 다른 거동이다.

또한, 저자들은 자기장 강도를 10¹⁶–10¹⁷ G 수준으로 설정했지만, 알펜 시간 척도에 대한 스케일링 관계를 이용해 실제 10¹⁴–10¹⁵ G 수준에서도 동일한 물리적 메커니즘이 작동할 것이라고 주장한다. 즉, 자기장 강도가 절반이면 불안정 성장 시간은 두 배가 되지만, 불안정의 종류와 최종 난류 상태는 변하지 않는다.

이러한 결과는 마그네터 형성 과정에서 토로이드 자기장이 장기적으로 유지되기 어렵고, 비축대칭 불안정에 의해 빠르게 난류화될 가능성을 제시한다. 따라서 관측적으로 마그네터가 보이는 거대한 자기장은 내부에서 지속적인 재생 메커니즘(예: MRI 혹은 대류)과 결합되어야 함을 암시한다.