[논문 리뷰] Characterization of exoplanets from their formation I. Models of combined planet formation and evolution

C. Mordasini, Y. Alibert외 2명|arXiv (Cornell University)|2012. 01. 01.

Stellar, planetary, and galactic studies참고 문헌 96인용 수 150

한 줄 요약

이 논문은 행성 형성에서 수십억 년 된 시스템에 이르기까지 태양계 간행성의 질량, 반지름, 구성, 반지름, 복사원수를 동시에 계산하는 자기일관성 있고 결합된 행성 형성 및 진화 모델을 제시한다. 모델은 단순화에도 불구하고 이전 모델과 강한 정량적 일치를 보이며, 특히 목성의 형성에서 복사원수 피크와 핵질량 변화와 같은 핵심적 특징을 성공적으로 재현한다. 이는 다중 관측 제약 조건 통합을 위한 강력한 인구합성에 기여한다.

ABSTRACT

A first characterization of many exoplanets has recently been achieved by the observational determination of their radius. For some planets, a measurement of the luminosity has also been possible, with many more directly imaged planets expected in the future. The statistical characterization of exoplanets through their mass-radius and mass-luminosity diagram is thus becoming possible. This is for planet formation and evolution theory of similar importance as the mass-distance diagram. Our aim in this and a companion paper is to extend our formation model into a coupled formation and evolution model. We want to calculate in a self-consistent way all basic characteristics (M,a,R,L) of a planet and use the model for population synthesis calculations. Here we show how we solve the structure equations describing the gaseous envelope not only during the early formation phase, but also during gas runaway accretion, and during the evolutionary phase at constant mass on Gyr timescales. We then study the in situ formation and evolution of Jupiter, the mass-radius relationship of giants, the influence of the core mass on the radius and the luminosity both in the "hot start" and the "cold start" scenario. We put special emphasis on the comparison with other models. We find that our results agree very well with those of more complex models, despite a number of simplifications. The upgraded model yields the most important characteristics of a planet from its beginning as a seed embryo to a Gyr old planet. This is the case for all planets in a synthetic planetary population. Therefore, we can now use self-consistently the statistical constraints coming from all major observational techniques. This is important in a time where different techniques yield constraints on very diverse sub-populations of planets, and where its is challenging to put all these constraints together in one coherent picture.

연구 동기 및 목표

• 행성 형성에서 진화에 이르기까지 질량, 반지름, 복사원수, 구성, 궤도 거리와 같은 모든 기본 행성 성질을 동시에 계산하는 통합 모델을 개발하는 것.
• 특히 Pollack 등 (1996)의 획기적인 연구를 기준으로 모델의 물리적 일관성과 정확도를 검증하는 것.
• 계산 효율적이면서도 물리적으로 탄탄한 프레임워크를 제공하여 합성 행성 인구를 생성하는 데에 인구합성 연구를 가능하게 하는 것.
• 디스크 제한된 기체 응집, 핵의 분화, 방사성 가열과 같은 핵심 물리 과정을 통합하여 현실성 향상에 기여하는 것.
• 이식, 중력파, 직접 영상 촬영 등의 다중 기술적 관측 제약 조건을 통합하여 행성계의 이론적 그림을 통합적으로 제공하는 것.

제안 방법

• 초기 형성, 급속 기체 응집, 질량이 일정한 장기적 열진화 단계의 세 단계에 걸쳐 원형행성 기체 대류층의 구조 방정식을 해결한다.
• 디스크 절단과 점성 진화를 고려한 새로운 디스크 제한 기체 응집률 규정을 도입한다.
• 변동하는 밀도와 구성을 가진 분화된 내부 구조 모델을 사용하며, 표준 목성 모델에서 핵 밀도는 14.31 g/cm³이다.
• 핵 내부의 방사성 붕괴를 내부 열원으로 통합하여 열진화 및 복사원수 예측을 향상시킨다.
• 관측된 복사원수 및 반지름 진화와 더 잘 일치하기 위해 흩어진 입자 투과도를 조정한 처리 방법을 적용한다.
• 1000 km의 태양계 원시행성에서 46억 년까지의 자기일관성 있는 진화를 수행하며, 질량, 반지름, 복사원수, 핵 질량을 시간에 따라 추적한다.

실험 결과

연구 질문

• RQ1 단순화되었지만 물리적으로 일관된 모델이 목성 형성 중 복사원수 진화와 핵질량 증가를 얼마나 정확하게 재현할 수 있는가?
• RQ2 핵의 분화와 방사성 가열이 거대행성의 예측 반지름과 복사원수에 어떤 영향을 미치는가?
• RQ3 응집 시간 상수와 초기 조건의 변화가 급속 기체 응집의 시기와 그로 인한 행성 질량 및 반지름에 어떤 영향을 미치는가?
• RQ4 이러한 복잡한 시뮬레이션에서 관측된 복사원수 피크, 특히 '뜨거운 시작'과 '차가운 시작' 시나리오를 얼마나 잘 재현하는가?
• RQ5 형성 역사를 관측 가능한 성질인 질량, 반지름, 복사원수와 일관되게 연결할 수 있는가를 고려할 때, 이 모델이 신뢰성 있게 인구합성 연구에 사용될 수 있는가?

주요 결과

• 모델은 목성 형성 중 첫 번째 및 두 번째 복사원수 피크를 높은 정확도로 재현한다. 표준 모델에서 두 번째 피크는 log(L/L☉) = -2.77에 발생하며, Lissauer 등 (2009)이 발견한 -2.4에서 -2.3의 범위와 유사하다.
• 급속 응집의 시작 시점(교차 시간)은 표준 모델에서 0.817 Myr에 예측되며, 핵 질량은 16.27 M⊕이다. 이는 초기 조건의 차이를 고려할 때 Pollack 등 (1996)의 3.32 Myr와 12.24 M⊕과 유사하다.
• 표준 목성 모델의 최종 질량은 316.6 M⊕이며, 반지름은 0.99 R_J, 복사원수는 1.13 L_J이다. 이는 목성의 관측값(318 M⊕, 1.0 R_J, 1.73 L_J)과 뛰어난 일치를 보인다.
• 표준 모델에서 핵 밀도는 14.31 g/cm³로 예측되며, J1a 모델의 3.2 g/cm³보다 훨씬 높다. 이는 최종 구조를 결정하는 데 핵의 압축의 중요성을 보여준다.
• 고체 응집을 인위적으로 중단하는 효과도 성공적으로 재현한다. J1a 모델에서는 3.38 Myr에 급속 응집이 시작되며, J1 모델의 7.49 Myr보다 이르게 발생한다. 이는 Pollack 등 (1996)의 결과와 정량적으로 일치한다.
• 현존하는 중심 온도는 표준 모델에서 1.76 × 10⁴ K로 예측되며, J1 모델의 2.16 × 10⁴ K와 일치하며, 목성 내부의 기대 범위와도 일치한다.