[논문 리뷰] Hydrodynamics and survivability during post-main-sequence planetary engulfment

Ricardo Yarza, Naela Razo Lopez외 7명|arXiv (Cornell University)|2022. 03. 20.

Astro and Planetary Science인용 수 4

한 줄 요약

이 연구는 3D 유체역학 시뮬레이션을 사용하여 별성장의 진화 단계에서 별성체(행성 및 브라운 왜성)가 별의 대기 중으로 삽입되는 과정을 모델링하며, 동반체를 강체이고 반사성인 물체로 간주하여 유동 역학을 정확히 반영한다. 연구 결과, 행성의 삽입은 별의 밝기를 수십만 배에서 수천만 배까지 증가시킬 수 있으며, 약 10 $M_{\text{Jup}}$ 정도의 작은 동반체라도 적색거성 가지의 정점에서 별의 대기층을 탈출시킬 수 있으며, 이는 마찰력과 궤도 감쇠가 측정 가능한 마찰 계수에 의해 결정됨을 밝혀냈다.

ABSTRACT

Data to reproduce the results of Yarza et al. 2022 (submitted to ApJ). For the main reproducibility repository, see https://gitlab.com/ryarza/planetary-engulfment-hydro After downloading all the files, run unpack.py to extract them. The script will also put all grid2 files into the correct folder structure. As a result of size limitations, this repository contains only the hydrodynamical simulations used in the figures in the paper. Simulations at other times or for other parameters are available upon request.

연구 동기 및 목표

• 거대 별의 후주계열에서 별성체(SBs)의 유체역학적 거동과 생존 가능성에 대해 연구하기 위해.
• 유성체가 별의 대기층 내에 잠겨 있을 때 작용하는 마찰력(압력 마찰력과 중력 마찰력)을 정량화하기 위해.
• 삽입 과정 중 궤도 감쇠와 에너지 전달을 모델링하여 밝기 증가와 대기층 탈출 가능성 평가하기 위해.
• 3D 유체역학에서 강체이고 반사성인 SB 주위의 유동 형태를 모의하기 위해 FLASH를 사용한 수치 프레임워크 개발하기 위해.
• 별의 반지름에 따라 대기층을 탈출시키기 위한 최소 SB 질량을 결정하는 경험적 스케일링 법칙 유도하기 위해.

제안 방법

• 적응 메쉬 해상도를 사용한 FLASH 4.6.2 코드를 활용하여 $1M_{\bigodot}$ 거대 별 내에 강체이고 반사성인 별성체(SB)의 3D 유체역학을 시뮬레이션한다.
• 유동 형태의 특징을 정확히 반영하기 위해 SB를 응축하지 않고 통과할 수 없는 물체로 모델링하며, 반사 경계 조건을 적용한다.
• 다양한 해상도에서 시간 평균 마찰 계수($C_p$: 압력 마찰력, $C_g$: 중력 마찰력)를 유체역학 시뮬레이션으로 측정한다.
• 시간에 따라 변화하는 마찰력과 궤도 내부에 포함된 질량의 변화를 고려하여 SB의 운동 방정식을 통합한다.
• 측정된 마찰 계수와 궤도 에너지 변화를 이용하여 궤도 감쇠 timescale과 대기층 탈출 임계 조건을 계산한다.
• 시뮬레이션 결과에 기반하여 대기층 탈출을 위한 최소 SB 질량과 복사권 아래 생존을 위한 최소 질량을 별의 반지름 함수로 경험적 공식에 맞춘다.

실험 결과

연구 질문

• RQ1 압력 마찰력과 중력 마찰력이 별성체가 삽입되는 동안 궤도 감쇠에 어떤 영향을 미치는가?
• RQ2 행성의 삽입이 $1M_{\bigodot}$ 거대 별의 밝기를 얼마나 증가시키며, 원래 상태로 회복하는 데 얼마나 오래 걸리는가?
• RQ3 다른 진화 단계에서 $1M_{\bigodot}$ 별의 대기층을 탈출시키기 위해 필요한 최소 질량의 별성체는 얼마인가?
• RQ4 별성체는 삽입 과정을 생존할 수 있으며, 복사권 아래 생존을 결정짓는 조건은 무엇인가?
• RQ5 강체이고 반사성인 물체로 간주할 때와 자유롭게 물질을 응축하는 물체일 때, 삽입된 SB 주위의 유체역학적 유동 구조는 어떻게 다를까?

주요 결과

• 궤도 감쇠에 의한 에너지 소산으로 인해 $1M_{\bigodot}$ 별의 밝기는 몇 개의 주기 수준에서 수천만 배까지 증가할 수 있다.
• 별이 약 ~10 $R_{\bigodot}$에 도달할 경우 원래 밝기로 복귀하는 데 최대 약 3,000년이 소요되며, 적색거성 가지 정점에서는 단지 약 30년이면 충분하다.
• 궤도 에너지가 유일한 탈출 원천이라면 별이 약 ~10 $R_{\bigodot}$에 도달하기 이전에는 어떤 SB도 대기층을 탈출시킬 수 없다; 탈출 가능성은 이후 단계에서만 나타난다.
• 약 10 $M_{\text{Jup}}$의 작은 별성체도 적색거성 가지 정점에서 별의 대기층을 탈출시킬 수 있으며, 이때 별은 약 ~190 $R_{\bigodot}$로 팽창해 있다.
• 대기층 탈출을 위한 최소 SB 질량은 별의 반지름에 비례하며, $11.45 \leq R_{\star}/R_{\odot} \leq 193.7$ 범위에서 시뮬레이션 데이터와 1.1%의 일치도를 보이는 피팅 공식을 유도하였다.
• 복사권 아래 생존을 위한 최소 SB 질량은 로그 공식으로 예측되며, $2 \leq R_{\star}/R_{\odot} \leq 193.7$ 범위에서 시뮬레이션 데이터와 7.9%의 일치도를 보였다.