헬리컬 자기장 MHD 충격이 만든 제트 회전

초록

이 논문은 헬리컬 자기장을 가진 제트에 MHD 충격이 발생하면, 압축된 자기장이 토러스 방향의 로렌츠 힘을 생성해 제트 물질을 회전시킨다는 가설을 수치 시뮬레이션으로 검증한다. 1.5 D·2.5 D 축대칭 PLUTO 모델에서 초기 회전이 없는 제트를 다양한 마그네틱·동역학 파라미터로 주입했으며, 충격 후 토러스 속도가 전체 흐름 속도의 0.1–1 % 수준으로 증가함을 확인했다. 이는 원시성별 제트에서 관측된 회전 속도와 일치한다.

상세 분석

본 연구는 기존의 ‘초기 회전 전달’ 모델이 관측된 제트 회전과 일치하지 않을 가능성을 지적하고, 충격 전후의 자기 토션 전이를 통한 각운동량 변환 메커니즘을 제시한다. 헬리컬 자기장은 폴로이드 성분 Bₚ와 토러스 성분 B_φ가 결합된 구조로, 충격이 발생하면 Bₚ가 압축되어 B_φ가 강화된다. 이때 로렌츠 힘 F_φ = (J × B)_φ 은 양(또는 음) 방향의 토러스 가속을 제공한다. 수식적으로는 MHD 점프 조건을 적용해 전단면에서 Δ(B_φ)·Δ(v_φ) ≈ (Bₚ·Δρ)/ρ 와 같은 관계가 도출되며, 이는 자기 각운동량 L_m = r B_φ Bₚ / (4πρ vₚ) 가 물질 각운동량 L_k = r v_φ 로 전환되는 과정을 의미한다.

시뮬레이션 설정은 PLUTO 코드의 고전적 TVD 스킴을 이용해 축대칭 1.5 D(반구면 좌표)와 2.5 D(R–Z 평면)에서 수행되었다. 초기 조건은 내부 알벤 마흔수 M_A, int 와 빠른 음속 마흔수 M_f 를 조절해 다양한 충격 강도를 구현했으며, 제트와 주변 매질의 밀도 대비 η = ρ_jet/ρ_amb 도 변형하였다. 모든 경우에서 제트는 v_φ = 0 로 시작했으며, 충격 전후에 토러스 속도 v_φ 가 급격히 상승하는 것을 확인했다. 특히 M_A, int ≈ 2–5, M_f ≈ 5–10 인 경우에 v_φ / v_jet ≈ 10⁻³–10⁻² 로, 관측된 원시성별 제트(예: DG Tau, HH 212)의 회전 속도와 일치한다.

또한 파라미터 스터디를 통해 η 가 낮을수록(제트가 가벼울수록) 전이된 각운동량이 크게 나타나며, 이는 외부 매질과의 상호작용이 강할수록 충격에 의한 자기 토션이 효율적으로 전달된다는 물리적 해석을 뒷받침한다. 반면, 알벤 마흔수가 매우 높으면 자기장 자체가 강하게 꼬여 충격 전후의 B_φ 변동이 제한되어 회전 효율이 감소한다.

이 메커니즘은 헬리컬 자기장이 존재하는 모든 제트—원시성별, 마이크로퀘이사, AGN 등—에 적용 가능하다. 특히 외부 충격(제트와 주변 가스의 충돌)과 내부 충격(제트 내부의 knot) 모두에서 동일한 각운동량 전환이 일어나므로, 관측된 회전이 반드시 초기 디스크 회전에서 직접 유래한 것이 아니라 충격에 의한 2차적인 현상일 수 있음을 시사한다.