다중상태 ISM의 열불안정과 폴리바리 압력 효과

초록

본 연구는 2차원 MHD 시뮬레이션을 통해 약하게 이온화된 ISM에서 열불안정(TI)을 조사한다. 배경 자기장과 압력 이방성을 포함한 폴리바리(다중압) 모델을 적용했으며, 압력의 평행·수직 비율이 위상(따뜻한 중성 매질, 차가운 중성 매질, 불안정 중성 매질)의 질량·부피 분포에 미치는 영향을 분석하였다. 자기장 세기는 전체 진화에 미미하지만, 압력 이방성 파라미터는 위상 비율을 크게 변화시켜 장기적인 TI 진화에 중요한 역할을 함을 확인하였다.

상세 분석

이 논문은 전통적인 CGL 이중단열 모델 대신, 관측에 근거한 폴리바리 압력 모델을 도입함으로써 ISM의 압력 이방성을 보다 현실적으로 묘사한다. 압력 텐서는 (p_{\perp})와 (p_{\parallel}) 두 성분으로 분리되며, 이들의 비율을 나타내는 압력 이방성 파라미터 (a_p = p_{\parallel}/p_{\perp} - 1)를 일정값으로 가정한다. 이는 밀도와 자기장 세기가 관측 범위(0.1–100 cm⁻³) 내에서 크게 변하지 않는다는 전제와 일치한다.

수치 해법으로는 다중유체 플럭스‑보정 전송(mFCT) 코드를 사용했으며, 이는 고차 정확도와 인공 확산 최소화를 동시에 달성한다. 방정식은 일반화된 연속식 형태로 변환되어 FCT 스킴에 적합하도록 구성되었다. 시간 적분은 2차 정확도의 스플릿‑스텝 방식을 적용했으며, 모든 물리량은 평형값으로 정규화하여 수치적 안정성을 확보하였다.

열손실 함수 (\mathcal{L}(\rho,T))는 Field(1965)의 형태를 따르며, 냉각·가열 항을 각각 (\Lambda(T))와 (\Gamma)로 구분한다. 전도 계수는 자기장 방향에 따라 이방성인 (K_{\parallel}=K_{\perp}=K_c T^{1/2}) 로 설정하였다. 초신성 피드백은 무작위 강제항 (\mathbf{f}_{SN})으로 모델링했으며, 이는 솔레노이달·압축성 성분 비율을 조절할 수 있다.

선형 단계에서는 1차원 축소 모델을 이용해 성장률과 임계 파장을 검증했으며, 이는 기존 이론과 일치함을 확인하였다. 비선형 2차원 시뮬레이션 결과는 초기 균일 상태에서 작은 온도·밀도 요동이 성장해 따뜻한 중성 매질(WNM), 차가운 중성 매질(CNM), 그리고 불안정 중성 매질(UNM)으로 구분되는 다중상태 구조를 형성한다. 특히, (a_p) 값을 0.0(등방성)에서 2.0(강한 평행 압력)까지 변화시켰을 때, UNM의 부피 비율이 20 %에서 45 %로 크게 증가하고, CNM의 질량 비율이 감소하는 현상이 관찰되었다. 이는 압력 이방성이 열전달과 냉각 효율을 변화시켜 임계 온도 구간을 이동시키기 때문이다. 반면, 자기장 세기 (B_0)를 0 µG에서 5 µG까지 변화시켜도 위상 비율 변화는 5 % 이하에 머물렀으며, 주된 역할은 흐름을 자기장 방향으로 정렬시키는 정도에 국한되었다.

결과적으로, 폴리바리 압력 모델은 ISM의 다중상태 형성에 중요한 매개변수이며, 특히 압력 이방성은 별 형성률에 영향을 미칠 수 있는 UNM 비율을 조절한다는 점에서 기존 등방성 MHD 시뮬레이션보다 물리적 타당성이 높다.