The CO snow line favours strong clumping by the streaming instability in protoplanetary discs with porous grains

기술 요약: 다공성 입자를 포함한 원시 행성 원반에서 CO 스노우 라인이 스트리밍 불안정성에 의한 강한 응집을 촉진하는 현상

문제 제기
원시 행성 원반 내의 먼지 입자로부터 미행성체가 형성되는 과정은 파쇄(fragmentation), 바운싱(bouncing), 그리고 반경 방향 표류(radial-drift) 한계라는 중대한 장벽에 직면해 있다. 스트리밍 불안정성(streaming instability, SI)은 현재 직접적인 미행성체 형성을 위한 주요 메커니즘으로 간주되지만, 그 효능은 특정 조건, 즉 국지적 먼지-가스 부피 비($\epsilon$)가 임계값($\epsilon_{crit}$)을 초과하고, 먼지와 가스의 결합이 미미한 상태(스토크스 수 $St \lesssim 1$)를 충족해야 한다는 점에 의존한다. 기존 연구들은 $\epsilon_{crit}$가 외부 난류와 입자 크기 분포에 민감하다는 것을 입증했으나, 다공성 입자를 포함하여 진화하는 원시 행성 원반에서 이러한 조건들이 자연스럽게 충족되는지는 아직 검증되지 않았다. 특히, 입자의 성장을 가속화하고 항력 체계를 변화시키는 입자 다공성과 스노우 라인(특히 CO 스노우 라인)의 존재가 강한 응집을 유발하는 능력에 대해 어떻게 상호작용하는지는 충분히 평가되지 않았다.

방법론
저자들은 Michoulier 등(2024b, "Paper I")에서 제시된, 성장, 파쇄, 압축 및 바운싱 과정을 모델링하는 다공성 규산염 입자의 진화를 다룬 3D SPH(Smoothed Particle Hydrodynamics) 시뮬레이션의 하위 집합을 활용하였다. 시뮬레이션은 $1 M_\odot$의 별을 공전하는 $0.01 M_\odot$의 원반을 대상으로 하며, 범위는 10 au에서 400 au까지이고 난류 점성 파라미터는 $\alpha = 5 \times 10^{-3}$이다.

두 가지 특정 사례가 분석되었다:

1. CO 스노우 라인이 없는 경우: 입자들이 균일한 파쇄 임계값($v_{frag} = 20$ m s$^{-1}$)을 가진다.
2. CO 스노우 라인이 있는 경우: CO 스노우 라인 너머($T = 20$ K, $r \gtrsim 100$ au)에서는 CO 얼음 코팅으로 인해 파쇄 임계값이 $v_{frag} = 5$ m s$^{-1}$로 감소한다.

수직 적분된 표면 밀도(금속 함량 $Z$) 대신, 저자들은 수직 침강을 고려하기 위해 국지적 먼지-가스 부피 비 $\epsilon$에 집중하였다. 저자들은 시뮬레이션된 $\epsilon$ 및 $St$ 값을 외부 난류를 고려한 Lim 등(2024)의 임계값과 비교하였다:
$$\log \epsilon_{crit} = 0.42 (\log St)^2 + 0.72 \log St + 0.37$$
강한 응집은 시뮬레이션된 $\epsilon > \epsilon_{crit}$인 영역에서 식별되었다.

주요 결과
$t = 15 de 150$ kyr에서의 시뮬레이션 분석 결과, CO 스노우 라인의 포함 여부에 따라 뚜렷한 차이가 나타났다:

• CO 스노우 라인이 없는 경우: 강한 응집을 위한 조건은 내부 원반 가장자리(약 15 au 이내)의 제한된 영역에서만 충족된다. 이 영역은 초기에 ($\sim 20$ kyr) 형성되지만, 300 kyr까지 축소되어 사라진다. 두 번째 강한 응집 영역은 중간 거리($\sim 25-90$ au)에서 결국 출현하지만, 이는 약 200 kyr 이후에 일어난다.
• CO 스노우 라인이 있는 경우: CO 스노우 라인의 존재는 미행성체 형성이 가능한 영역을 크게 확장시킨다. 강한 응집은 초기에 ($\sim 70$ kyr) 발생하며, 내부 가장단에서 $\sim 80$ au까지 확장된다. 이 영역은 15 au와 40 au를 중심으로 하는 두 개의 높은 먼지 농축 존( $\epsilon$ 피크 $\sim 5$)을 포함한다. 외곽 영역은 중간면 위로 몇 au까지 확장된다.

이러한 차이를 유도하는 메커니즘은 CO 스노우 라인이 장벽 역할을 한다는 것이다. 이는 입자가 너무 빠르게 성장하는 것을 방지하고 스토크스 수를 $St \lesssim 0.1$로 제한한다. 결과적으로 입자들은 더 느리게 표류하며 원반의 더 먼 곳까지 축적되어, 비-스노우라인 사례에 비해 원반 진화 초기 단계부터 훨씬 더 넓은 영역에서 향상된 $\epsilon$을 생성한다.

의의 및 주장
본 논문은 다공성 입자의 진화가 CO 스노우 라인과 결합될 때, 이러한 요인들을 무시한 모델들보다 훨씬 더 이른 시기에, 그리고 더 넓은 영역에서 스트리밍 불안정성에 의한 강한 응집을 위한 유리한 조건을 조성한다고 주장한다.

저자들은 다공성 입자만으로도 응집을 유도할 수 있지만, 다공성 입자를 포함한 원반에서 미행성체 형성을 촉진하는 데 있어 CO 스노우 라인의 포함이 "핵심(key)"이라고 단언한다. 이는 스노우 라인이 더 큰 입자들을 외곽 원반에 유지하도록 도와, 입자들이 느리게 표류하게 하고 SI가 효과적으로 작동할 수 있는 영역에 축적되도록 하기 때문이다. 저자들은 이 메커니즘이 물 스노우 라인의 효과와 물리적으로는 구별되지만, 행성 형성을 촉발하는 데 있어 유사한 중요성을 공유한다고 언급한다.

한계 및 주의사항
저자들은 사용된 임계값 $\epsilon_{crit}$가 단분산 입자 분포(mSI)를 위해 도출된 것인 반면, 본인들의 시뮬레이션은 다분산 분포(PSI)를 포함하고 있다는 점을 명시적으로 인정하였다. 비록 큰 입자에 정점을 두는 분포가 mSI와 유사한 비선형 포화 상태와 선형 성장을 보일 수 있다는 문헌을 인용하였으나, 다분산에 의한 감쇠(dampening)가 성장률을 낮출 수 있음을 인정하였다. 또한, 다분산 시스템에서 SI를 안정화하거나 불안정화하는 자기장의 역할은 추가 조사가 필요한 미결 과제로 남아있다고 언급하였다. 본 연구는 진화하는 원반 내에서 직접적인 SI 시뮬레이션을 수행하는 것이 아니라, 도출된 먼지 특성을 바탕으로 불안정성의 *실행 가능성(viability)*을 평가하는 것이다.