Kick Velocities and Mass Function of Free-Floating Planets from Dynamical Ejection in Hierarchical Three-Body Systems

기술 요약: 계층적 삼체계의 역학적 방출에 의한 자유 부유 행성의 킥 속도 및 질량 함수

문제 정의
자유 부유 행성(FFP) 또는 떠돌이 행성은 모항성에 묶여 있지 않은 아항성 천체이다. 이들의 속도 및 질량 분포는 이들의 형성 및 방출 이력을 보여주는 운동학적 화석 역할을 한다. 중력 미세중력렌즈 현상은 지구 질량에서 목성 질량에 이르는 FFP 집단을 발견해 왔으나, 방출 속도와 그에 따른 질량 함수를 지배하는 역학적 메커니즘은 여전히 불완전하게 규명되어 있다. 구체적으로, 계층적 삼체계(모항성, 거대 질량 교란자, 그리고 더 가벼운 내부 행성)의 방출 결과를 관측 가능한 은하계 FFP 집단의 특성, 즉 점근적 속도($v_\infty$) 및 방출 시간($t_{ejec}$)과 연결하는 작업이 필요하다.

방법론
저자들은 고정밀 근접 조우 시의 정확도를 보장하기 위해 rebound 적분 패키지의 IAS15 적분기(15차 적응형 타임스텝)를 사용한 직접 $N$-체 시뮬레이션 세트를 제시한다. 본 연구는 내부 행성을 태양 질량 모항성($M_\star = 1 M_\odot$)과 거대 질량 교란자 사이의 테스트 입자로 취급하는 제한된 삼체 문제를 다룬다.

시뮬레이션 그리드는 통제되고 체계적이다:

• 고정 파라미터: 모항성 질량($1 M_\odot$), 거대 행성 장반경($a_J = 5.2$ AU, 목성 유사체), 초기 궤도 경사각($0^\circ$–$3^\circ$).
• 변화 파로미터:
  • 방출되는 행성 질량 ($m_p$): $10^{-5}$ ~ $10^{-1} M_{Jupiter}$.
  • 거대 행성 질량 ($M_J$): $1$~$75 M_{Jupiter}$.
  • 방출되는 행성의 장반경 ($a_p$): $0.7$~$1.5 a_J$.
  • 두 행성의 이심률 ($e_p$ 및 $e_J$): $0.0$~$0.9$.
• 실현(Realizations): 각 구성은 파라미터 효과를 격리하기 위해 평균 근점 이각(mean anomalies)에 공유된 무작위 시드를 사용하여 $N=10,000$개의 독립적인 실현으로 시뮬레이션된다.
• 지표: 연구는 점근적 방출 속도($v_\infty$), 방출 확정 시점의 유한 반경 킥 속도($v_{kick}$), 그리고 방출 시간($t_{ejec}$)을 측정한다. 방출은 특정 에너지 양의 양수 상태, 시스템의 최대 장반경을 초과하는 거리, 그리고 10 궤도 주기 동안 지속되는 외향 반경 방향 속도에 의해 정의된다.

분석 프레임워크
수치적 결과는 다음을 결합한 반분석적 프레임워크를 통해 해석된다:

1. 에너지 예산: 방출에 필요한 에너지 증분은 묶인 궤도 에너지와 점근적 운동 에너지 사이의 차이로부터 도출된다.
2. 중력 슬링샷: 최대 킥 속도는 교환 속도 $u$와 거대 행성의 궤도 속도에 의존하는 거대 행성의 정지 좌표계에서의 슬링샷 메커니즘에 의해 제한된다.
3. 힐-규모 산란(Hill-Scale Scattering): 중간값 방출 속도는 거대 행성의 힐 구(Hill sphere) 내에서의 전형적인 조우에 의해 결정되며, $v_H \propto M_J^{1/3}$로 스케일링된다.
4. 티세랑 파라미터(Tisserand Parameter): 궤도 요소와 교환 속도를 연관시키기 위한 보존 진단 도구로 사용된다.
5. 세큘러 진화(Secular Evolution): 불안정성 시간은 세큘러 펌핑 및 평균 운동 공명(MMR) 중첩과 연결된다.

주요 결과

1. 테스트 입자 한계: 방출되는 행성의 질량($m_p$)은 네 자릿수 범위($10^{-5}$ ~ $10^{-1} M_{Jupiter}$)에 걸쳐 방출 통계에 미치는 영향이 미미하다. 방출 분율은 2.3% 미만으로 변하며, 평균 속도는 안정적으로 유지된다. 이는 $m_p/M_J \lesssim 10^{-2}$인 경우 역학이 거대 교란자에 의해 지배되며, 방출되는 물체가 테스트 입자처럼 행동함을 확인시켜 준다.
2. 거대 행성 질량 스케일링: 거대 행성의 질량($M_J$)은 방출 규모와 시간을 결정한다.
   • 평균 방출 속도는 $\langle v_\infty \rangle \propto M_J^{0.42}$로 스케일링된다.
   • 방출 시간은 $\langle t_{ejec} \rangle \propto M_J^{-1.09}$로 스케일링되며, 이는 세큘러 예측인 $t_{sec} \propto M_J^{-1}$와 일치한다.
   • 방출 분율은 비단조적(non-monotonic)이며, 매우 낮은 질량과 매우 높은 질량 영역에서 각각 세큘러 시간 제한과 힐 구 포획으로 인해 정점을 찍고 감소한다(정점은 약 $10 M_{Jupiter}$ 근처).
3. 장반경 의존성: 공명 중첩 구역 내의 구성($|a_p - a_J|/a_J \lesssim 0.40$, $10 M_{Jupiter}$ 교란자 기준)에서 방출 시간은 세큘러 진화보다는 빠른 카오틱 중첩 시간 timescale에 의해 지배되어 평탄하게 나타난다 ($\sim 1.9 \times 10^4$ 년). 평균 속도는 $a_p = 0.9 a_J$ 근처에서 정점에 도달한다.
4. 이심률 효과 (이중 속도 응답):
   • 내부 행성 이심률 ($e_p$): 주로 고속 꼬리(high-velocity tail)에 영향을 미친다. $e_p$가 0에서 0.9로 증가하면 최대 킥 속도는 약 3배(약 7.8에서 23 km/s로) 증가하지만, 평균 속도는 약 1.7배 증가하는 데 그친다.
   • 거대 행성 이심률 ($e_J$): 고속 꼬리에 훨씬 더 강력한 영향을 미친다. $e_J$가 0에서 0.9로 증가하면 최대 킥 속도는 약 10배(약 7.8에서 80 km/s로) 증가한다. 이는 거대 행성의 근점 속도($v_{peri} \propto \sqrt{(1+e_J)/(1-e_J)}$)의 비선형적 성장에 의해 구동되며, 이를 통해 드물지만 고에너지 슬링샷 조우가 가능해진다.
5. 분석적 상한 포락선(Analytic Upper Envelope): 저자들은 슬링샷 메커니즘과 티세ран 파라미터를 기반으로 최대 방출 속도에 대한 분석 공식을 도출하였다. 이 포락선은 특히 $e_p \gtrsim 0.5$인 경우 시뮬레이션된 속도 분포의 천장을 정확하게 예측한다.

의의 및 시사점
본 논문은 이러한 결과들이 행성계의 구조와 관측 가능한 은하계 FFP 집단의 특성 사이에 직접적인 연결 고리를 제공한다고 주장한다:

• 질량 함수: 방출 분율이 행성 질량에 독립적이므로, 현재의 FFP 질량 함수는 근본적인 행성 발생률($dN_{FFP}/dm_p \approx \bar{f}_{eject} \cdot dN_{planet}/dm_p$)을 직접 추적한다. 이는 만약 행성 발생률이 낮은 질량 쪽으로 높아진다면(Kepler 데이터가 시사하는 바와 같이), 하위 지구 질량의 FFP 집단이 수적으로 우세할 것임을 의미한다.
• 미세중력렌즈 신호: 변환된 미세중력렌즈 시간 분포($dN/d\log t_E$)는 행성 발생률이 평탄하면 $\log t_E$에 대해 평탄할 것이며, 발생률이 낮은 질량 쪽으로 높아지면 짧은 시간대에서 강화될 것이라고 예측된다. 이는 Roman Wide Field Survey와 Euclid를 위한 테스트 가능한 예측을 제공한다.
• 운동학: 방출된 행성 대부분은 점근적 속도가 불과 몇 km/s에 불 불과하여, 모항성 집단과 운동학적으로 구별하기 어렵다(3D 속도 분산을 약 1% 변화시킴). 그러나 매우 이심률이 높은 거대 교란자를 가진 시스템($e_J \gtrsim 0.5$)은 고속 꼬리($v_\infty \sim 80$ km/s)를 생성할 수 있으며, 이는 Gaia나 Roman을 통해 탐지 가능한 고유 운동 이상치(proper-motion outliers)를 만들 수 있다.
• 방출 시간: 질량이 크거나 이심률이 높은 시스템으로부터의 대부분의 방출은 매우 빠르게($10^3$–$10^4$ 년) 발생하며, 종종 전주계열 단계 동안 발생한다. 결과적으로 일반적인 FFP 집단의 공간 분포는 일반적인 항성계와 구별하기 어려우며, 오직 최근에 방출된 젊은 천체들만이 탄생 환경의 기억을 유지한다.

한계점
저자들은 정확히 세 개의 천체로 제한한 것이 주요 한계점이라고 언급한다. 실제 시스템은 연쇄적인 불안정성을 유발하는 여러 행성을 포함할 수 있다. 또한, 본 시뮬레이션에는 원시행성계 원반, 항성 질량 손실, 또는 다양한 경사각(특히 코자이-리도프 메커니즘)이 포함되지 않았으며, 이는 방출 통계에 변화를 줄 수 있다. 향-후 연구에서는 이러한 요소들을 완전한 인구 합성 모델(population synthesis model)에 통합하는 작업이 필요하다고 명시하였다.