MRI 난류가 만든 비표준 α 디스크의 안정성

초록

본 논문은 MRI‑구동 난류의 효율이 미세 점성·전도도 비인 magnetic Prandtl number(Pm)에 따라 달라진다는 최신 시뮬레이션 결과를 도입해, α‑파라미터를 α∝Pm^δ 형태의 비표준 형태로 가정한 얇은 원형 디스크 모델을 구축한다. Spitzer식 미세 확산계수를 사용해 Pm이 반경에 따라 감소함을 보이고, δ값이 클수록 외부 디스크에서 α와 난류 가열이 크게 억제된다. 결과적으로 δ>2/3 일 때 방사압이 거의 없더라도 중간 반경에서 열·세속 불안정이 발생함을 확인한다.

상세 분석

이 연구는 MRI‑구동 난류가 전통적인 α‑디스크 모델의 핵심 가정인 “α는 일정하다”는 전제를 깨고, 미세 점성(ν)과 전기 저항성(η)의 비인 magnetic Prandtl number(Pm) 의 함수임을 실험·수치 결과에 기반해 제시한다. 기존 로컬 시어링 박스 시뮬레이션(Lesur & Longaretti 2007; Simon & Hawley 2009)은 α∝Pm^δ (δ>0) 관계를 발견했으며, 이는 점성 증가가 MRI 성장률을 억제함에도 불구하고 비선형 난류에서는 전반적인 전단 응력이 강화된다는 역설적인 현상을 반영한다. 논문은 Spitzer 전도 및 점성 공식 ν∝ρ⁻¹T^{5/2}, η∝T^{-3/2} 를 채택해 Pm∝T^4 n^{-1} 형태를 얻고, 온도가 감소함에 따라 Pm이 급격히 감소한다는 물리적 직관을 제공한다.

δ값을 매개변수화함으로써 디스크 전반에 걸친 α의 반경 의존성을 구체화한다. δ=0이면 전통적인 Shakura‑Sunyaev 모델을 회복하고, δ가 0.25, 0.5, 1.0 등으로 증가하면 외부 반경에서 α가 급격히 낮아져 점성 가열이 약화되고, 그 결과 표면 밀도 Σ가 증가하며 방사압이 지배적인 내부 영역과는 달리 가스 압력이 우세한 중·외부 영역이 확대된다. 특히, Pm이 1을 초과하는 임계 반경 r_c≈100 r_s 근처에서 α가 급격히 상승하는데, 이는 MRI가 효율적으로 작동하는 “Pm>1” 구역이 디스크 내부에 한정된다는 점을 시사한다.

안정성 분석에서는 중력 불안정(Q≈1), 열 불안정(∂Q⁺/∂T>∂Q⁻/∂T), 세속 불안정(∂(ν_tΣ)/∂Σ<0) 세 가지 고전적 기준을 적용한다. δ가 커질수록 α∝Pm^δ 로 인한 점성 감소가 외부에서 Σ를 높이고, 따라서 토러스 질량이 집중돼 중력 불안정이 쉽게 발생한다. 더 흥미로운 점은 δ>2/3 일 때, 방사압이 무시될 정도로 낮은 경우에도 열·세속 불안정이 중간 반경에서 발생한다는 것이다. 이는 전통적인 α‑모델에서는 방사압이 우세할 때만 불안정이 나타나는 것과 대조된다. 따라서 MRI‑구동 난류의 Pm 의존성이 디스크의 장기 진화와 변광 현상(예: dwarf nova, X‑ray transient)의 핵심 메커니즘이 될 가능성을 제시한다.

마지막으로, 실제 천체 디스크(AGN, X‑ray binary)에서는 Pm이 수십에서 수백까지 변동하므로, δ가 0에 가까운 약한 의존성을 유지해야만 관측된 안정적인 질량이송을 설명할 수 있다. 강한 Pm 의존성(δ≈1)은 디스크가 지속적인 평형 상태를 유지하기 어렵게 만들며, 주기적인 폭발이나 퇴화 현상을 야기할 수 있다.