자기장 중성자별에 대한 본드 호일 흡수 시뮬레이션

초록

본 연구는 초음속으로 이동하는 자기쌍극자 중성자별이 주변 별간 매질을 흡수하는 과정을 2차원 축대칭 MHD 시뮬레이션으로 조사한다. 자기반경이 본드-호일 반경보다 작을 때, 자기장은 흡수 흐름을 억제하고, 레위드(뒷면) 극으로 물질이 더 많이 흐른다. 약한 자기장에서는 자기권이 앞뒤로 진동하며 흡수율이 3배 정도 변동한다.

상세 분석

본 논문은 마하수 M = 3인 초음속 흐름 속에서, 자기쌍극자 형태의 자기장을 가진 고정된 중성자별이 별간 매질(ISM)을 흡수하는 과정을 축대칭 저항성 MHD 방정식으로 직접 계산하였다. 시뮬레이션은 별의 자기반경 Rₘ이 본드-호일 반경 R_BH보다 작아 중력 초점이 지배적인 상황을 가정했으며, 이는 전통적인 본드-호일 이론(˙M_BH = πR_BH² ρ v)과 비교할 수 있는 기준을 제공한다.

핵심 물리적 파라미터는 다음과 같다. 별 질량 M ≈ 1.4 M⊙, 사운드 속도 c_s, 별 속도 v = 100 km s⁻¹(정규화 v₁₀₀ = 1), 그리고 표면 자기장 B_* = 10¹⁰–10¹² G에 해당하는 자기쌍극자 모멘트 μ = 10³⁰ G cm³(무차원 μ ≈ 1.6 × 10⁻⁸). 이때 자기압력과 동압의 균형을 이루는 반경 Rₘ ≈ 0.02 R_B (시뮬레이션 단위)이며, 이는 수치별 반경 R_* = 0.0025 R_B보다 약 8배 크다.

시뮬레이션 결과는 크게 세 가지 현상을 보여준다. 첫째, 자기장이 존재하면 전체 흡수율 ˙M 이 비자기 경우보다 약 30 %~ 50 % 감소한다. 이는 자기권이 흐름을 차단하고, 물질이 별 표면에 도달하기 전까지 억제되기 때문이다. 둘째, 물질은 별의 양극(레위드) 쪽으로 집중적으로 흐르며, 이는 중력 초점이 뒤쪽에 형성된 고밀도 ‘뒤쪽 꼬리’와 연계된다. 셋째, 약한 자기장(μ < 2 × 10⁻⁸)에서는 자기권이 앞뒤로 진동하며, 이때 압력 균형이 주기적으로 깨지면서 흡수율이 ~ 3배까지 변동한다. 진동 주기는 별이 통과하는 시간 t₀ 에 비례하며, 전자기 확산계수 ηₘ = 10⁻⁵ (무차원)로 설정된 경우 약 5–10 t₀ 정도의 주기를 보인다.

또한, 시뮬레이션은 전단층(Shock) 구조와 온도·압력 프로파일을 상세히 제시한다. 앞쪽 충격파는 z ≈ ‑0.04 R_B에서 형성되고, 물질 압력 p와 자기압력 p_mag 의 비율 β 은 충격 전후로 급격히 변한다. 충격 뒤쪽에서는 압력이 상승하고, 물질이 극 방향으로 가속되어 초음속 흐름을 회복한다. 이러한 구조는 비자기 경우와 유사하지만, 자기장에 의해 ‘귀’ 모양의 변형이 발생하고, 이는 레위드 쪽으로 물질이 집중되는 메커니즘을 강화한다.

마지막으로, 저자들은 시뮬레이션 결과를 실제 고령 중성자별(IONS)의 X‑ray 방출 가능성에 연결한다. 비자기 경우보다 낮은 흡수율에도 불구하고, B ≈ 10¹¹ G 수준의 자기장을 가진 오래된 별은 여전히 ˙M ≈ 10⁸–10⁹ g s⁻¹ 수준의 물질을 흡수할 수 있으며, 이는 L_X ≈ 10³¹–10³² erg s⁻¹ 정도의 X‑ray 광도를 생성한다. 별의 회전축, 자기축, 이동축이 서로 다른 경우, 이러한 광도는 회전 주기에 따라 변조될 것으로 예상된다.