Modelling dust coagulation, dynamical drag and turbulent mixing during star and disc formation

기술 요약: 별 및 원반 형성 과정 중의 먼지 응집, 역학적 항력 및 난류 혼합 모델링

문제 제기
먼지는 별과 행성 형성 모두에서 중요한 역할을 하며, 가스 냉각, 이온화 수준, 그리고 자기장과 가스의 결합에 영향을 미칩니다. 먼지가 더 큰 알갱이(자갈)로 응집되는 것은 행성 미행성 및 핵 성장에 필수적이지만, 기존의 별 형성 3차원 유체역학 시뮬레이션은 가스에 대한 먼지의 역학적 진화를 대체로 무시해 왔습니다. Bate (2022)의 초기 방법들은 먼지를 각 유체 요소(SPH 입자) 내에서 독립적으로 진화하는 것으로 취급하며, 먼지가 가스에 완벽하게 결합되어 있다고 가정했습니다. 이러한 근사는 초기 저밀도 붕괴 단계에서는 유효했으나, 밀도가 더 높고 입자 성장이 빠른 원시 행성 원반 내의 장기적인 먼지 이동, 침강 및 혼합을 설명하는 데는 실패했습니다. 또한, 기존 모델들은 종종 입자 성장에 필요한 상대 속도와 먼지를 수직으로 섞어주는 난류 확산을 유발하는 메커니즘을 어떻게 자가 일관적으로 처리할 것인가에 대한 처리가 부족했습니다.

방법론
저자들은 세 가지 뚜렷한 물리적 과정을 결합한 새로운 SPHNG 코드 내 구현 방식을 제시합니다:

1. 먼지 응집 (Dust Coagulation): 이산화된 평균장 Smoluchowski 방정식 방법(Hirashita & Yan 2009; Bate 2022)을 활용하여, 코드는 모든 SPH 가스 입자에 대해 $N$개의 로그 간격 크기 빈(bin)에 걸쳐 먼지 입자 집단의 진화를 추적합니다. 이 방법은 특정 응집 임계 속도 미만에서의 충돌 확률을 1로 가정하여 이진 충돌을 고려합니다.
2. 역학적 항력 (MULTIGRAIN): 저자들은 터미널 속도 근사법을 사용하여 가스와 다중 먼지 종의 혼합에 대한 연속체 유체 방정식을 푸는 MULTIGRAIN 알고리즘(Hutchison et al. 2018)을 통합합니다. 이를 통해 압력 구배와 중력으로 인한 먼지 크기별 벌크 차동 속도를 계산할 수 있습니다.
3. 암시적 적분 및 난류 확산: 핵심적인 혁신은 먼지 항력 방정식을 풀기 위해 최근 개발된 암시적 적분법(Elsender & Bate 2024)을 사용하는 것입니다. 이는 정지 시간이 긴 입자들에 대한 엄격한 타임스텝 제약을 피하고 비물리적인 음수 먼지 분율 발생을 방지합니다. 또한, 저자들은 난류 확산을 모델링하기 위한 새로운 암시적 방법을 도입했습니다. 그들은 먼지의 난류 교반을 확산 과정($D_i = \alpha_{SS} c_s H / Sc$)으로 취급하며, 여기서 슈미트 수($Sc$)는 서로 다른 입자 크기의 결합 효율을 나타냅니다. 이 확산 항은 Gauss-Seidel 연산자 분할(operator-splitting) 방식을 사용하여 해결되며, 이를 통해 먼지가 유체 요소 간에 섞일 수 있게 합니다.

코드는 분자 구름 핵의 붕괴부터 제1차 정수압 핵 단계를 거쳐 원시 행성 원반의 형성 및 진화까지 먼지의 진화를 추적하며, 별의 핵 형성 이후의 시스템을 추적하기 위해 싱크 입자(sink particles)를 사용합니다.

주요 기여

• 최초의 3D 결합 시뮬레이션: 본 연구는 별과 원반 형성의 초기 단계 동안 먼지 응집, 역학적 가스 항력, 그리고 난류 혼합을 동시에 모델링하는 최초의 3차원 유체역학 시뮬레이션입니다.
• 암시적 확산 스킴: 난류 먼지 확산에 대한 견고한 암시적 솔버의 개발은 명시적 확산 스킴과 관련된 심각한 타임스텝 제한 없이 안정적인 먼지 혼합 시뮬레이션을 가능하게 합니다.
• 통합된 처리: 코드는 동일한 난류 파라미터($\alpha_{SS}$)를 사용하여 상대 속도(성장을 유도)와 난류 확산(혼합을 유도)의 처리를 하나의 프레임워크 내에서 통합합니다.

결과
저자들은 테스트 계산(수직 기둥에서의 먼지 침강)과 전체 별 형성 시뮬레이션 결과를 모두 제시합니다:

• 테스트 계산: 코드는 먼지 침강과 난류 확산 사이의 평형에 대한 해석적 해를 성공적으로 재현합니다. 테스트 결과, 난류 확산을 포함하는 것이 비현실적으로 얇은 먼지 층 형성을 방지하고 유한한 먼지 척도 높이를 유지함을 확인했습니다.
• 별 형성 시뮬레이션:
  • 제1차 정수압 핵: Bate (2022)의 선행 연구와 일치하게, 제1차 정수압 핵 내부에서 급격한 먼지 성장이 일어나며, 핵 깊은 곳에서 크기가 100 $\mu$m를 초과하는 국지적 단분산(mono-disperse) 집단이 생성됩니다. 항력과 확산의 포함은 짧은 진화 시간으로 인해 이 초기 단계에 미미한 영향을 미쳤습니다.
  • 원시 행성 원반 진화: 원반이 형성되면 난류 확산의 포함이 먼지 진화를 크게 변화시킵니다.
    • 성장 촉진: 난류 확산은 먼지 입자의 성장률을 실질적으로 향상시킵니다. 확산은 작고 중간 크기의 입자들을 (큰 입자들이 성장하고 있는) 고밀도 영역으로 지속적으로 공급함으로써 응집을 위한 국지적 결핍을 방지합니다.
    • 난류 의존성: 난류 수준($\alpha_{SS}$)이 매우 중요합니다. 높은 난류($\alpha_{SS} = 10^{-3}$)는 원반 전체에 걸쳐 $\sim30-100$ $\mu$m 크기의 단분산 집단으로의 빠른 성장을 이끕니다. 낮은 난류($\alpha_{SS} = 10^{-4}$)는 더 느린 성장과 더 수직적으로 층이 나뉜 먼지 집단을 결과로 가져오며, 더 큰 입자들은 중간면(mid-plane) 근처에 갇히게 됩니다.
    • 먼지 침강: 약한 난류가 있는 시뮬레이션에서는 먼지 침강이 중간면의 먼지-가스 비율을 크게 증가시키는 반면(50% 이상), 작은 입자들은 중간면 근처에서 결핍되고 표면층에 갇히게 됩니다.
  • 항력 vs. 확산: 특정 과정의 포함 여부에 따른 시뮬레이션 비교 결과, 입자들이 잘 결합되어 있기 때문에 먼지-가스 역학적 항력은 이러한 시간 척도 내에서 먼지 분포에 거의 영향을 미치지 않았습니다. 반면, 난류 확산의 포함은 먼지 분포를 매끄럽게 하고 응집을 가속화하는 데 엄청난 효과를 보였습니다.

의의 및 주장
본 논문은 이 새로운 방법론이 별과 그 원반의 형성 과정 중 먼지 성장, 혼합 및 이동을 최초로 자가 일관적으로 모델링할 수 있게 한다고 주장합니다. 저자들은 항력과 확산의 포함이 제1차 정수압 핵에서의 먼지 성장에 관한 근본적인 발견을 바꾸지는 않지만, 후속 원시 행성 원반의 진화에는 결정적이라는 점을 강조합니다. 특히, 난류 먼지 확산이 큰 입자가 형성되는 지역에 작은 입자를 공급하는 원천 역할을 함으로써 먼지 입자의 성장률을 실질적으로 향상시킨다고 주장합니다.

저자들은 다음과 같은 점을 언급하며 연구 범위에 대해 겸허한 태도를 유지합니다:

• 결과는 채택된 난류 모델(예: $\alpha$-디스크 vs. MHD 난류)에 따라 달라집니다.
• 현재 시뮬레이션은 응집은 고려하지만 파쇄(fragmentation)는 고려하지 않습니다. 다른 응집 임계값을 사용하면 파쇄와 작은 입자의 재공급이 발생할 수 있습니다.
• 복사 전달 계산은 여전히 진화하는 국지적 입자 크기 분포 대신 표준 성간 불투명도를 가정하고 있으며, 이는 향후 반드시 수행해야 할 과제로 식별되었습니다.

궁극적으로, 이 연구는 엣지온(edge-on) 원반에서의 먼지 침강 및 난류 교반과 관련하여 ALMA, VLT-SPHERE, JWST와 같은 시설의 풍부한 새로운 관측 데이터와 3차원 유체역학 시뮬레이션을 비교할 수 있는 도구를 제공합니다.