Formation of Isotopically Heterogeneous Molecular Cloud Cores in Filamentary Molecular Clouds

기술 요약: 필라멘트 구조 분자운 내 동위원소 불균질적 분자운 핵의 형성

문제 제기
운석 분석 결과, 초기 태양계 고체 물질(예: 칼슘-알루미늄 풍부 함유물, CAI)은 동위원소 불균질성(예: $^{54}\text{Cr}/^{52}\text{Cr}$ 및 $^{50}\text{Ti}/^{47}\text{Ti}$ 비율의 변화)을 나타낸다. 열적 처리 과정이 원시 행성계 원반에서 발생했다는 가설과 불균질한 모 분자운 핵으로부터 물려받았다는 가설이 모두 제시되었으나, 필라멘트형 분자운이 중력 붕괴하여 밀집된 핵을 형성하는 과정에서 초기 동위원소 불균질성이 어느 정도까지 유지되는지는 명확하지 않다. 구체적으로, 필라멘트형 구름에 존재하는 초기 동위원소 변이가 중력 파쇄 및 핵 형성 과정 동안 얼마나 보존되거나 평활화(smoothing out)되는지, 그리고 이러한 불균질성이 결과적인 성주위 원반(circumstellar disk)으로 어떻게 전달되는지를 정량화할 필요가 있다.

방법론
저자들은 필라멘트형 분자운의 중력 파쇄를 모델링하기 위해 Smoothed Particle Hydrodynamics(SPH) 시뮬레이션, 구체적으로 Godunov smoothed particle magnetohydrodynamic(GSPM) 방법을 사용하였다.

• 시뮬레이션 설정: 시뮬레이션은 관측값과 일치하는 임계 선질량($M_{\text{line,crit}} \approx 18 M_\odot \text{pc}^{-1}$)과 폭($\sim 0.1$ pc)을 가진 필라멘트를 모델링하였다. 초기 조건에는 정수압 평형 밀도 프로파일, Kolmogorov 난류($\sigma = 2c_s$), 그리고 필라멘트 축에 수직인 $10 \mu\text{G}$의 초기 자기장이 포함되었다. 비교를 위해 유체역학 전용(hydrodynamic-only) 케이스도 실행되었다.
• 동위원소 추적: 동위원소 종은 수동 추적자(passive tracers)로 취급되었다. 저자들은 서로 다른 축($z$, $x, y$ 축 및 이들의 조합)을 따라 초기 동위원소 비율 구배($C_0$)를 부여하기 위해 6개의 별도 모델(표 1)을 구축하였다. 동위원소 비율은 계산 영역 전체에 대해 0과 1 사이로 정규화되었다.
• 분석: SPH 입자를 추적함으로써 동위원소 비율의 진화를 추적하였다. 저자들은 형성되는 핵의 밀도 정점 주변에 구형 껍질(spherical shells)을 정의하여 평균 동위원소 비율($\bar{C}$)과 중심으로부터의 차이($|\Delta C|$)를 계산하였다. 또한, 난류 속도장과 밀도 파워 스펙트럼을 기반으로 껍질의 초기 질량 중심 변위($\delta z_{\text{ini}}$)를 추정하기 위해 반해석적(semi-analytical) 접근법을 채택하였다.

주요 결과

1. 불균질성의 생존: 시뮬레이션은 필라멘트 파쇄 과정에서 동위원소 불균질성이 완전히 평활화되지 않음을 보여준다. 핵은 초기 필라멘트 길이 1 pc에 걸쳐 존재했던 초기 동위원소 불균질성의 약 **0.1%에서 1%**를 유지한다.
2. 기하학적 의존성: 잔류하는 불균질성의 크기는 초기 동위원소 구배의 필라멘트 기하학적 구조에 대한 방향에 따라 달라진다.
   • 종방향 축($z$-방향, 모델 $M_z$)을 따르는 구배는 껍질 간 가장 큰 동위원소 차이($|\Delta C| \sim 0.01$)를 생성한다.
   • 단축(또는 부축, $x$ 또는 $y$)을 따르는 구배는 현저히 작은 차이를 보인다. 이는 필라멘트의 기하학적 구조가 단축 방향으로의 유체 요소 변위를 제한하여, 서로 다른 동위원소 비율을 가진 물질 간의 혼합을 제한하기 때문이다.
3. 불균질성의 메커니즘: 껍질 간 잔류 동위원소 차이의 주요 동인은 난류 속도장에 의한 껍질의 초기 질량 중심 차이이다. 초기 유체 요소의 변위를 고려한 반해석적 추정치는 시뮬레이션 결과와 2배 이내의 오차로 일치하며, 이는 초기 난류 구성이 최종 동위원소 프로파일을 결정함을 확인시켜 준다.
4. 자기장의 영향: 자기장($10 \mu\text{G}$)의 존재는 순수 유체역학적 케이스와 비교했을 때 동위원소 비율의 내부 프로파일을 크게 변화시키지 않는다.
5. 성주위 원반에 대한 함의: 유입되는 물질의 원심 반지름을 추정함으로써, 저자들은 동위원소 불균질성이 성주위 원반에도 지속됨을 발견하였다. 원반은 핵의 불균질성을 물려받지만, 각운동량의 비등방적 분포로 인해 가스는 다양한 구형 껍질로부터 유입된다. 결과적으로 초기 필라멘트 불균질성의 약 **1%**가 원반에 남아 있음을 시사한다.

의의 및 주장
본 논문은 운석에서 관찰되는 동위원소 이상(anomalies)과 모 분자운의 천체물리학적 조건 사이의 정량적 연결 고리를 제공한다고 주장한다. 저자들은 자신들의 연구 결과와 Jacquet 등(2019)이 보고한 감소 계수(핵에서 원반으로 갈 때 약 3배 감소)를 결합하여, 최종 성주위 원반에 남은 불균질성은 초기 필라멘트 변이의 약 0.03%에서 0.3% 수준임을 추론한다.

저자들은 이 정도의 잔류 불균질성이 극단적인 열적 처리나 붕괴 후의 이례적인 핵합성 주입 시나리오 없이도 태양계 물질(예: $^{54}\text{Cr}/^{52}\text{Cr}$ 및 $^{50}\text{Ti}/^{47}\text{Ti}$ 변이인 1~10 parts in 10,000)에서 관찰되는 동위원소 이상을 설명하기에 충분하다고 단언한다. 본 연구는 적절한 수준의 초기 불균질성을 가진 모 구름이 전제된다면, 태양계 구성 성분의 동위원소 불균질성이 필라멘트 파쇄에서 원반 형성까지의 전 과정을 자연스럽게 통과하여 생존할 수 있음을 결론짓는다. 저자들은 향och적인 유입 및 유출(infall and outflows)의 비등방적 특성을 완전히 해결하기 위해 고해 resolution 시뮬레이션이 향후 작업에 필요하다고 언급하였다.