[논문 리뷰] Catching drifting pebbles II. A stochastic equation of motions for pebbles

Chris W. Ormel, Beibei Liu|arXiv (Cornell University)|2018. 03. 16.

Astrophysics and Star Formation Studies참고 문헌 88인용 수 31

한 줄 요약

이 논문은 원형행성계 원반 내 페블 크기의 입자에 대한 확률적 운동 방정식(SEOM)을 개발하며, 난류 기체 저항력과 확산을 포함하여 수직 방향 이동과 축적 효율성을 모델링한다. SEOM은 난류가 주로 수직 확산을 통해 축적을 억제하며, 높은 상대 속도에서는 정착 메커니즘이 실패함으로써 영향을 미친다는 것을 드러내며, 특히 외부 원반 영역과 낮은 질량을 가진 별 주변에서 더 강한 영향을 미친다.

ABSTRACT

Turbulence plays a key role in the transport of pebble-sized particles. It also affects the ability of pebbles to be accreted by protoplanets, because it stirs pebbles out of the disk midplane. In addition, turbulence can suppress pebble accretion once the relative velocities become too large for the settling mechanism to be viable. Following Paper I, we aim to quantify these effects by calculating the pebble accretion efficiency $\varepsilon$ using three-body simulations. To model the effect of turbulence on the pebbles, we derive a stochastic equation of motion (SEOM) applicable to stratified disk configurations. In the strong coupling limit (ignoring particle inertia) the limiting form of this equation agrees with previous works. We conduct a parameter study and calculate $\varepsilon$ in 3D, varying pebble and gas (turbulence) properties and accounting for the planet inclination. We find that strong turbulence suppresses pebble accretion through turbulent diffusion, agreeing within factors of order unity with previous works. Another reduction of $\varepsilon$ occurs when the turbulent rms motions are large and the settling mechanism fails. Efficiency-wise, the outer disk regions are more affected by turbulence than the inner regions. At the location of the H$_2$O iceline, planets around low-mass stars achieve much higher efficiencies. Including the results from Paper I, we present a framework to obtain $\varepsilon$ under general circumstances.

연구 동기 및 목표

• 기체와의 에어로다이내믹 결합과 난류 확산을 고려하여 원형행성계 원반 내 페블의 확률적 운동을 모델링하기.
• 난류가 중간 평면 농도를 해체하고 상대 속도를 증가시켜 페블 축적 효율성(ε)을 감소시키는 방식을 정량화하기.
• 입자 정지 시간, 기체 확산도, 상관 시간, 중력력 등을 포함하는 일반적인 ε 프레임워크를 유도하기.
• 원반의 분층 구조, 난류 특성(예: αz) 및 행성 기울기가 축적 결과에 미치는 영향을 조사하기.
• 논문 I의 결과를 확장하여 다양한 원반 및 행성 조건에 적용 가능한 3차원, 난류 해상도가 확보된 ε에 대한 규정을 제공하기.

제안 방법

• 난류 기체 결합으로 인한 확률적 힘을 포함한 뉴턴 역학을 확장하여 페블에 대한 확률적 운동 방정식(SEOM)을 유도한다.
• 두 가지 핵심 난류 매개변수인 기체 확산도(D_gas)와 상관 시간(t_corr)을 도입하여 입자 운동의 통계적 성격을 제어한다.
• 수직으로 변화하는 난류 확산도(αz(z))를 가진 분층 구조의 원반 구성을 적용하여 중간 평면 집중과 표면 확산을 허용한다.
• 세 체계 시뮬레이션을 통해 SEOM을 수치적으로 해결하여 다양한 난류 조건, 입자 크기(τs), 행성 질량에서의 페블 축적 효율성(ε)을 계산한다.
• 강한 결합 근사(limit t_stop → 0)에서 SEOM을 검증하여 이전 연구에서 백색 잡음(상관 없음) 난류를 가정한 결과와 일치함을 보였다.
• 논문 I의 2차원 표현과 결과를 통합하여 페블, 행성, 원반 매개변수에 따라 변화하는 ε에 대한 통합적이고 일반적인 규정을 구축한다.

실험 결과

연구 질문

• RQ1 난류 확산은 원형행성계 원반 내 페블의 수직 분포와 중간 평면 집중에 어떤 영향을 미치는가?
• RQ2 정착 메커니즘이 실패할 경우 난류 상대 속도가 페블 축적을 얼마나 억제하는가?
• RQ3 난류 특성의 변화(예: αz, t_corr)가 다양한 원반 위치와 행성 질량에서 축적 효율성 ε에 어떤 영향을 미치는가?
• RQ4 난류 확산과 난류 속도 분산 중 어느 것이 ε 감소에 더 중요한가?
• RQ5 행성 질량, 기울기, 별 질량이 난류 존재 하에서 축적 효율성 ε에 어떻게 영향을 미치는가?

주요 결과

• 페블 축적 효율성(ε)을 억제하는 주요 메커니즘은 확산으로, 이는 페블을 중간 평면에서 떨어지게 하여 국소 밀도를 감소시킨다.
• τs ≈ 0.1–1 인 입자에 대해 이상적인 MRI 실행에서는 암페어디온 확산 실행보다 ε가 낮은 편이지만, 이는 더 강한 난류 속도 효과(σ ≳ v_*)로 인해 정착 메커니즘이 약간 실패하기 때문이다.
• SEOM 프레임워크는 t_stop → 0 근사에서 강한 결합 근사(Scale)를 회복하여 이전 모델이 백색 잡음 난류를 가정한 결과와 일관성을 보였다.
• 중력 집중 효과가 강하고 페블 밀도가 높기 때문에 내부 원반에서 축적이 더 효율적이며, 외부 영역은 난류로 인해 더 심하게 억제된다.
• 낮은 질량의 별 주변에서는 물의 얼음선(H2O iceline)에서 더 큰 페블 포착 반경과 낮은 난류 속도 분산으로 인해 더 높은 ε를 달성한다.
• 제시된 프레임워크는 난류 효과(확산 및 속도)와 중력 집중을 통합하는 일반적인 3차원 ε 규정을 가능하게 하며, 다양한 원반 및 행성 조건에 적용 가능하다.