초기 원시 원반 진화와 초기 구름 핵의 연결 고리

초록

본 논문은 원시 원반(프로토디스크) 단계에서 별이 형성된 직후, 원반 자체 중력이 지배적인 시기를 자기유사해(self‑similar) 모델로 분석한다. 베르누이‑에버트 구름 핵의 붕괴에 의해 공급되는 물질 유입률과 중력 토크, 원반 풍선에 의한 각운동량 손실을 주요 파라미터로 삼아, 토루멘 불안정성, 표면밀도·온도·특이각운동량 등의 시간·반경 프로파일을 도출한다. 핵의 불안정도가 클수록 원반은 토루멘 불안정에 빠지기 쉬우며, 특이각운동량은 10¹⁹–10²⁰ cm² s⁻¹ 범위에 머문다. 모델은 약 2 × 10³ yr 동안 신뢰할 수 있으며, 약 5–10 kyr까지는 약한 자기제동·강한 물질 유입 조건에서 연장 가능하다.

상세 분석

이 연구는 원시 원반 단계, 즉 원시 별이 형성된 직후 디스크 질량이 별 질량보다 크게 지배적인 시기를 대상으로 한다. 저자들은 기존의 이소열 붕괴 모델이나 단순 점성 파라미터화와 달리, 원반 자체 중력에 의해 유도되는 비선형 토러스 토크를 ‘중력 점성(ν_GI)’으로 명시하고, 이를 자기유사해 접근법에 포함시켰다. 핵심 방정식은 질량·운동량 보존식에 풍선(디스크 풍) 손실 항을 추가한 형태이며, 2 π r σ v_r와 2 · σ̇_w가 각각 질량 흐름과 풍선에 의한 질량 손실을 나타낸다.

자기유사 변수 x = K⁻¹/2 G^{γ‑1 / 2} r t^{γ‑2} 로 변환함으로써 편미분 방정식을 상미분 방정식(ODE) 체계로 축소한다. 여기서 K와 γ는 폴리트로픽 상수·지수이며, 저자는 γ = 1.1(≈등온)으로 설정해 디스크가 차가운 상태에서 중력 불안정을 유지하도록 했다.

중력 점성 ν_GI = η G σ² r⁶ j⁻³ 형태는 Lin & Pringle(1987)의 토러스 토크 모델을 차용했으며, η는 효율 파라미터(본 연구에서는 η′ = 10⁻²)로 설정돼 토러스에 의한 각운동량 전달 속도를 조절한다. 풍선 파라미터 Λ = Λ₀ x^s (s = 0.7, Λ₀ = 0.1)와 l = 1.5는 풍선이 각운동량을 얼마나 효율적으로 빼앗는지를 나타낸다.

경계 조건은 중심에서 v_r = 0, 외곽에서 물질 유입률 Ṁ_infall을 고정(베르누이‑에버트 구름 핵 붕괴에 기반)한다. 구름 핵의 안정도 파라미터 α₀ = 0.5(중간 불안정) 경우, Ṁ_infall ≈ 2.8 × 10⁻⁵ M_⊙ yr⁻¹이며, 이는 핵이 더 불안정(α₀ = 0.3)일 때는 6.1 × 10⁻⁵ M_⊙ yr⁻¹까지 증가한다. 이렇게 물질 공급이 강할수록 원반 표면밀도 σ가 급격히 상승하고, 토루멘 Q = c_s κ/(π G σ) 값이 1 이하가 되는 영역이 확대돼 파편화 가능성이 높아진다.

시뮬레이션 결과는 0.5, 1, 2, 3 kyr 시점에서 σ, v_r, v_φ, Q, M_disk/M_env 등의 프로파일을 제시한다. 초기에는 r ≈ 5 au 내부에서 σ가 급증하고, 외부(≈ 70 au)까지 디스크가 점차 팽창한다. 회전 속도 v_φ는 케플러(v ∝ r⁻¹/²)보다 얕은 v_φ ∝ r⁻⁰·³ 정도를 보이며, 이는 중력 토러스가 각운동량을 외부로 전달하기 때문이다.

또한, 특이각운동량 j(r) ≈ 10¹⁹–10²⁰ cm² s⁻¹ 범위에 머무는 것이 관측된 Class 0 원반과 일치한다는 점을 강조한다. 모델은 2 × 10³ yr 동안은 별 질량이 무시될 정도로 M_disk ≫ M_* 가 유지된다고 가정하지만, 약한 자기제동과 지속적인 고속 물질 유입이 유지될 경우 5–10 kyr까지 확장 가능하다고 제시한다.

결론적으로, 초기 구름 핵의 불안정도(α₀)와 물질 유입률이 원시 원반의 토루멘 안정성, 크기, 회전 프로파일을 결정짓는 핵심 요인임을 밝히며, 자기유사해 모델이 복잡한 중력·풍선 상호작용을 정량적으로 해석하는 유용한 도구임을 입증한다.