목성 중간 자기권‑전리층 결합: 자기디스크 자기장 모델을 포함한 자기‑전류 상호작용 연구

초록

본 논문은 목성의 중간 자기권에서 플라즈마 디스크가 생성하는 자기장 구조를 자기‑전리층(M‑I) 결합 전류와 자기 일관적으로 결합시켜 계산한다. 전리층 페데르센 전도도와 이오니아 플라즈마 질량 유출률이 전류 강도와 필드‑정렬 전류 위치에 미치는 영향을 분석하고, 관측된 오벌 오로라 위치 변동을 설명하기 위해 플라즈마 질량 밀도의 큰 변화가 필요함을 제시한다.

상세 분석

이 연구는 기존의 고정된 자기장 모델(단극 혹은 경험적 전류시트) 대신, Caudal(1986)의 자기디스크 모델을 자기‑전리층 결합 방정식에 자체적으로 연동시켰다. 핵심은 Hill‑Pontius 이론을 이용해 플라즈마 각속도 ω(ρ) 를 구하고, 이를 자기디스크의 Euler 포텐셜 α(r,θ)와 연계함으로써 자기장 B와 전류 j_φ 를 동시에 수치해석한다. 전리층 페데르센 전도도 Σ*_P 와 플라즈마 질량 유출률 Ṁ 은 Hill 거리 ρ_H 를 결정하는 주요 파라미터이며, ρ_H ∝ (Σ*_P/Ṁ)^{1/4} 로 정의된다. Σ*_P 가 증가하거나 Ṁ 가 감소하면 ρ_H 가 커져 플라즈마가 더 멀리까지 강하게 회전하게 되고, 이에 따라 원심력에 의한 원심 전류가 외부 자기권에서 우세해진다. 반면, 내부 영역(ρ≲15 R_J)에서는 고온 플라즈마 압력에 의한 전류가 지배적이며, 이는 기존 연구와 일치한다.

전류 구조는 두 부분으로 나뉜다. (1) 하부 전류는 전리층에서 하강하는 페데르센 전류 i_P = 2Σ*P B_J Ω_J ρ_i(1−ω/Ω_J) 로 표현되며, 이는 전자기적 연속성에 의해 방사형 전류 i_ρ 로 전환된다. (2) 방사형 전류는 I_ρ = 8π Σ*P Ω_J F_e (1−ω/Ω_J) 로 계산되고, 그 발산에 의해 필드‑정렬 전류 j∥(θ) = (B_J/2πρ_e|B_z|) dI_ρ/dρ_e 가 도출된다. 자기디스크 모델을 적용하면 |B_z| 가 외부에서 플라즈마 디스크의 프린징 필드에 의해 반전되며, 이로 인해 외부 영역에서 j∥ 가 부호를 바꾸어 하강 전류가 상승 전류로 전환된다.

플라즈마 압력 프로파일은 두 구성요소로 나뉜다. 고온(~30 keV) 플라즈마는 압력 구배가 원심력을 압도해 전류를 주도하고, 저온(~100 eV) 플라즈마는 원심력에 의해 압축되어 원심 전류를 생성한다. 저온 플라즈마의 밀도 N_c(ρ) 와 온도 T_c(ρ)는 Voyager 관측값을 기반으로 구간별 선형 보간을 적용하였다.

모델 결과는 전리층 페데르센 전도도와 질량 유출률을 변화시켰을 때 전류 강도는 크게 변하지만, 필드‑정렬 전류의 피크가 발생하는 적도 거리 ρ_peak 은 비교적 일정하게 유지된다는 점을 보여준다. 이는 오벌 오로라의 위도 변동이 전도도·Ṁ 변화만으로는 설명되지 않으며, 플라즈마 질량 밀도(특히 저온 플라즈마)의 10배 이상 변동이 필요함을 시사한다. 또한, 외부 프린징 필드가 전류를 역전시키는 현상은 관측된 외부 오로라(극광) 구조와 일치한다.

결론적으로, 자기디스크와 플라즈마 각속도를 자기‑전리층 결합에 자기 일관적으로 포함함으로써, 전류 강도와 위치, 그리고 오로라 매핑 사이의 복잡한 상호작용을 보다 정밀하게 재현할 수 있었다. 이는 향후 목성 외곽 자기권의 동역학 및 외계 행성의 M‑I 결합 모델링에 중요한 기반을 제공한다.