[논문 리뷰] Constraining the properties of the potential embedded planets in the disk around HD 100546

Max Ackermann Pyerin, Timmy N. Delage외 4명|arXiv (Cornell University)|2021. 10. 07.

Astrophysics and Star Formation Studies참고 문헌 53인용 수 13

한 줄 요약

이 연구는 자가교정된 ALMA 0.9 mm 먼지 연속분포 관측을 통해 HD 100546의 태양계 형성 중인 원반에서 두 개의 내재된 행성의 질량과 궤도 반경을 제약한다. 유체역학적 시뮬레이션(FARGO3D)과 합성 관측(RADMC-3D)을 조합하여 최적의 맞춤 모델은 내부 행성(13 au, 8 MJup)과 외부 행성(143 au, 3 MJup)을 규명하며, 외부 고리의 낮은 빛의 세기를 내부 고리와 비교해 낮게 유지하기 위해 자가유사적이고 지수적으로 퇴색하는 기체 표면 밀도 프로파일이 필요하다.

ABSTRACT

Context. The protoplanetary disk around the star HD 100546 displays prominent substructures in the form of two concentric rings. Recent observations with the Atacama Large Millimeter/sub-millimeter Array (ALMA) have revealed these features with high angular resolution and have resolved the faint outer ring well. This allows us to study the nature of the system further. Aims. Our aim is to constrain some of the properties of potential planets embedded in the disk, assuming that they induce the observed rings and gaps. Methods. We present the self-calibrated $0.9\,$mm ALMA observations of the dust continuum emission from the circumstellar disk around HD 100546. These observations reveal substructures in the disk that are consistent with two rings, the outer ring being much fainter than the inner one. We reproduced this appearance closely with a numerical model that assumes two embedded planets. We varied planet and disk parameters in the framework of the planet-disk interaction code FARGO3D and used the outputs for the gas and dust distribution to generate synthetic observations with the code RADMC-3D. Results. From this comparison, we find that an inner planet located at $r_1 = 13\,$au with a mass $M_1 = 8 M_{ m{Jup}}$ and an outer planet located at $r_2 = 143\,$au with a mass $M_2 = 3 M_{ m{Jup}}$ leads to the best agreement between synthetic and ALMA observations (deviation less than $3\sigma$ for the normalized radial profiles). To match the very low brightness of the outer structure relative to the inner ring, the initial disk gas surface density profile needs to follow an exponentially tapered power law (self-similar solution), rather than a simple power-law profile.

연구 동기 및 목표

• HD 100546 원반에 존재할 잠재적 내재된 행성의 질량과 궤도 간격을 제약하는 것.
• 관측된 하위구조, 특히 어두운 외부 고리를 재현하기 위해 하나 또는 두 개의 행성이 필요한지 확인하는 것.
• 내부 고리와 외부 고리 간의 관측된 밝기 비율을 일치시키기 위해 필요한 원반의 초기 기체 표면 밀도 프로파일을 규명하는 것.
• 현재 ALMA 및 직접 영상 장치로 이러한 행성을 탐지할 수 있는지 평가하는 것.

제안 방법

• 자기교정된 0.9 mm ALMA 관측을 통해 HD 100546의 먼지 연속분포를 추출하고 하위구조를 식별하였다.
• FARGO3D 코드를 사용하여 기체와 먼지의 유체역학적 시뮬레이션을 수행하였으며, 다양한 행성 질량과 궤도 반경을 가진 행성-원반 상호작용을 모델링하였다.
• RADMC-3D 코드를 사용하여 시뮬레이션 출력물에서 합성 관측을 생성하여 ALMA 유사 관측 조건을 시뮬레이션하였다.
• 이미지 형태와 반경 방향 밝기 프로파일을 분석하여 모델 출력물을 실제 ALMA 데이터와 비교하였다.
• 자기유사적이고 지수적으로 퇴색하는 거듭제곱법 초기 기체 표면 밀도 프로파일이 외부 고리의 낮은 밝기를 재현하기 위해 필요하다는 것이 밝혀졌다.
• 디스크 점도(α), 행성 질량, 궤도 위치를 변화시켜 반복적으로 모델 피팅을 수행하여 관측된 반경 방향 프로파일과의 편차를 최소화하였다.

실험 결과

연구 질문

• RQ1 관측된 하위구조, 특히 어두운 외부 고리를 가장 잘 재현하는 행성 질량과 궤도 반경은 무엇인가?
• RQ2 두 개의 행성 모델이 내부 고리와 외부 고리를 모두 설명할 수 있는가, 아니면 하나의 행성으로도 충분한가?
• RQ3 내부 고리와 외부 고리 간의 관측된 밝기 비율을 일치시키기 위해 필요한 초기 원반 기체 표면 밀도 프로파일은 무엇인가?
• RQ4 디스크 점도(α)는 추정된 행성 질량과 모델 적합도에 어떤 영향을 미치는가?
• RQ5 예측된 행성 질량은 현재의 ALMA 또는 직접 영상 장치로 탐지 가능한가?

주요 결과

• 합성 관측과 ALMA 관측 간의 최적 일치는 내부 행성(13 au, 8 MJup)과 외부 행성(143 au, 3 MJup)을 가진 모델에서 달성되며, 정규화된 반경 방향 프로파일에서 3σ 이내의 편차를 보였다.
• 외부 고리의 낮은 밝기를 내부 고리와 비교해 낮게 유지하기 위해 초기 기체 표면 밀도 프로파일이 자가유사적이고 지수적으로 퇴색하는 거듭제곱법을 따라야 한다.
• 디스크 점도 α = 10⁻³일 경우 외부 행성 질량은 3 MJup로 추정되었으며, 더 낮은 점도(α = 10⁻⁴)일 경우 질량은 1.7 MJup로 추정되었다.
• 모델은 외부 고리의 애자일 비대칭성이 없음을 성공적으로 재현하였으며, 이는 행성 유도 압력 함정 메커니즘을 지지한다.
• 추정된 행성 질량은 현재 ALMA 탐지 한계(4 Myr된 시스템 기준 33 MJup) 이하이므로 기존 장비로 직접 탐지하기는 어려울 것이다.
• 결과는 산란된 빛 영상(약 100 au에서 나선형 패턴)과 분자선 관측(약 15 AU에서의 동반 천체) 결과와 일치하며, 두 개의 거대 행성이 존재한다는 것을 뒷받침한다.