목성 질량 원시행성의 무거운 원소 함량과 궤도 위치 관계

초록

이 연구는 원반 불안정에 의해 형성된 목성 질량 원시행성이 별로부터의 거리, 원반 질량, 표면밀도 지수(α)에 따라 얼마나 많은 무거운 원소를 흡수할 수 있는지를 계산한다. 0.01, 0.05, 0.1 M☉ 원반과 α = 0.5, 1, 1.5를 적용해 1 × 10⁵ 년 동안의 초기 팽창 단계에서의 고체 흡수를 모델링하였다. 결과는 행성당 1 ~ 110 M⊕의 무거운 원소를 획득할 수 있음을 보여주며, 원반 표면밀도가 높을수록 풍부한 금속 함량이 나타난다. 이는 관측된 외계 가스행성들의 금속 함량 다양성과 별의 금속성 상관관계를 설명한다. 또한 목성도 현재보다 더 바깥에서 형성돼도 내부 모델이 요구하는 무거운 원소량을 충분히 흡수할 수 있음을 시사한다.

상세 분석

본 논문은 디스크 불안정(DI) 모델 하에서 목성 질량(J ≈ 1 M_J) 원시행성이 형성 후 초기 10⁵ 년 동안 고체 물질을 얼마나 효율적으로 포획할 수 있는지를 정량적으로 평가한다. 핵심 가정은 원시행성이 아직 팽창된 상태이며, 이때는 반경이 수십 R_J에 달해 고체 입자와의 충돌 단면적이 크게 증가한다는 점이다. 원반의 고체 표면밀도 σ_s는 σ_s = σ₀ a^(−α) 형태로 기술되며, 여기서 a는 별로부터의 거리, α는 0.5, 1, 1.5의 세 값을 탐색한다. α가 2보다 작을 경우 σ_s는 거리와 반비례하지만, 전체 고체 질량은 원반 면적이 a²에 비례해 증가하므로, 원거리일수록 가용 고체 총량이 더 많다. 그러나 고체 입자들의 궤도 주기와 상대 속도가 감소하면서 포획 효율은 급격히 떨어진다. 저자들은 이 두 효과를 결합해, 고체 입자들의 체적 밀도와 상대 속도, 그리고 원시행성의 부피 팽창률을 고려한 충돌 확률 함수를 도입하였다.

시뮬레이션은 세 가지 원반 질량(0.01, 0.05, 0.1 M☉)과 세 가지 α값을 조합해 9가지 경우를 만든다. 각 경우마다 a = 5, 10, 20, 30 AU에서의 포획 가능한 무게를 계산하고, 포획이 효과적인 기간을 10⁵ 년으로 제한한다. 이 기간 이후에는 원시행성이 수축하고, 원반에 기공(gap)이 형성되어 물질 공급이 급격히 감소한다는 가정이다. 결과는 다음과 같다. 낮은 원반 질량(0.01 M☉)에서는 α = 0.5일 때도 최대 5 M⊕ 정도만 흡수되지만, 높은 원반 질량(0.1 M☉)과 α = 1.5 조합에서는 30 AU에서 80 M⊕ 이상을 포획한다. 전반적으로 α가 클수록(표면밀도가 급격히 감소) 내부에서의 고체 농도가 높아져 근거리에서의 포획 효율이 상승하지만, 원거리에서는 전체 고체량이 충분히 커서 α가 작아도 상당량을 흡수한다.

이러한 결과는 관측된 외계 가스행성들의 금속 함량이 10 M⊕에서 100 M⊕까지 넓은 범위에 걸쳐 있음을 자연스럽게 설명한다. 별의 금속성이 높은 시스템은 일반적으로 원반 표면밀도가 높아지므로, DI 모델에서도 높은 무거운 원소 함량을 기대할 수 있다. 또한 목성 자체에 적용하면, 현재 궤도인 5.2 AU보다 외곽(≈10–20 AU)에서 형성되었더라도, 초기 팽창 단계에서 충분히 20–30 M⊕ 정도의 무거운 원소를 흡수해 내부 구조 모델이 요구하는 금속량을 만족시킬 수 있다.

결론적으로, DI 모델에서 최종 행성의 금속 함량은 원반의 질량·표면밀도 분포·형성 거리라는 세 가지 환경 변수에 크게 좌우된다. 따라서 금속 함량만으로는 행성 형성 메커니즘(핵 획득 vs. 디스크 불안정)을 구분하기 어렵다.