An Increase in the Galactic Planet Host Fraction Fails to Reproduce the Galactic Height Trend in Planet Occurrence

기술 요약: 은하계 행성 보유 분율의 증가가 은하 고도에 따른 행성 발생률의 경향을 재현하는 데 실패함

문제 제기
최근의 관측 결과, 구체적으로 Zink 등(2023)(이하 Z23)은 은하 평면으로부터의 최대 수직 진동 진폭($Z_{max}$)이 증가함에 따라 FGK 왜성 주위의 작은 행성(슈퍼-지구 및 서브-해 Neptune) 발생률이 유의미하게 감소하는 강력한 경향을 확인했다. 항성 금속함량이 행성 형성과 상관관계가 있다는 것은 알려져 있으나, 은하 금속함량 구배(gradient)만으로는 관찰된 급격한 감소를 설명하기에 불충분하다. 본 연구가 다루는 핵심 질문은 은하계 역사의 특정 시점에 발생한 시간 의존적인 행성 보유 분율($f$)의 증가가 관찰된 occurrence-$Z_{max}$ 경향을 재현할 수 있는가 하는 점이다. 저자들은 과거 은하계의 어느 시점에서 발생한 행성 형성 효율의 "부스트(boost)"가 시간에 따른 항성 궤도의 운동학적 가열(kinematic heating)과 결합될 때, 높은 은하 고도에서의 행성 결핍을 설명할 수 있는지 조사한다.

방법론
저자들은 합성 행성계 인구 생성 코드를 시뮬레이션하기 위해 psps를 개발하였다. 방법론은 다음과 같은 단계로 구성된다:

1. 항성 샘플: 연령 불확실성에 대한 결과의 견고성을 테스트하기 위해 두 가지 별개의 항성 샘플을 구축하였다:

   • 실제 샘플 (Real Sample): Hardegree-Ullman 등(2025) 카탈로그에서 추출되고 Berger 등(2020b, 2025)의 등온 연령(isochrone ages)과 교차 검증된 Kepler 및 K2 합동 샘플(22,184개 항성). $Z_{max}$ 값은 Gaia DR3 운동학과 Gala 패키지를 통한 은하 퍼텐셜 모델을 사용하여 계산되었다.
   • 합성 샘플 (Synthetic Sample): Kepler 및 K2 영역을 모사하도록 조정된 TRILEGAL 인구 생성 샘플(49,109개 항성). 이 샘플은 임의로 정밀한 항성 연령과 알려진 운동학적 가열 모델을 제공하여, 저자들이 연령 정밀도의 효과를 물리적 모델로부터 분리할 수 있게 한다.

2. 시간 의존적 발생 모델: 항성 연령(역행 시간, lookback time)의 함수로서 행성 보유 분율($f$)의 진화를 테스트하기 위해 두 종류의 모델을 시험하였다:

   • 계단 함수 (Step Functions): 특정 역행 시간 $t$에서 초기 분율 $f_1$에서 현재의 분율 $f_2$로 급격히 도약하는 형태.
   • 구간 함수 (Piecewise Functions): 시간 $t$에서 시작하여 $f_1$에서 $f_2$로 점진적으로 선형 증가하는 형태.
   • 모든 모델은 관측된 제약 조건과 일치하도록 현재의 행성 보유 분율을 약 30~35%로 제한하였다.

3. 인구 생성 및 관측: 저자들은 규정된 $f$를 바탕으로 psps를 사용하여 합성 행성계를 생성하였다. 시스템에는 기존의 인구 통계적 제약(예: Lam & Ballard 2024)에 따라 구조(조밀한 다행성계 vs 단일 행성)와 궤도 매개변수(이심률, 경사각)가 할당되었다. 생성된 인구는 기하학적 통과 확률과 Kepler 및 K2 미션의 특정 감도 함수를 고려한 역 탐지 효율 방법(IDEM)을 사용하여 "관측"되었다.

4. 분석: 시뮬레이션된 발생률을 $Z_{max}$ (100–1000 pc)별로 빈(bin)을 나누고, 기울기($\tau$)와 정규화 상수($\eta$)를 결정하기 위해 멱법칙($y \propto Z_{max}^\tau$)으로 피팅하였다. 이 결과들을 Z23의 관측된 경향($\tau \approx -0.28$)과 비교하였다.

주요 결과

• 보유 분율 증가의 불충분성: 본 연구는 행성 보유 분율($f$)을 시간에 따라 증가시키는 것(즉각적이든 점진적이든)이 일반적으로 관찰된 occurrence-$Z_{max}$ 경향의 가파른 정도를 재현하기에는 불충분하다는 것을 발견했다.
• 모델 제약:
  • 실제 Kepler-K2 샘플의 경우, 대다수의 모델이 관측된 $\tau \approx -0.28$보다 훨씬 완만한 기울기(0에 가까운 값)를 나타냈다. 관측된 경향에 가장 근접한 모델들은 극단적이고 물리적으로 타당하지 않은 매개변수(예: $f_1 \approx 1\%$ 또는 $f_2 \approx 100\%$)와 2.2~3.2 Gyr 사이의 좁은 창에 제한된 계단 시간($t$)을 요구했다.
  • 정밀한 연령을 가진 TRILEGAL 합성 샘플의 경우에도 결과는 유사했다. 비록 후기의 즉각적인 $f$ 상승이 Z23의 기울기를 맞출 수는 있었으나, 이 역시 이론적 극한치($f_1=1\%$ 또는 $f_2=100\%$)로 튜닝을 필요로 했다. 이는 완벽한 연령 정밀도를 갖추더라도 단순한 $f$의 증가만으로는 관찰된 경향을 주도할 수 없음을 시사한다.
• Kepler 전용 분석: 다른 연령 분포를 가지고 K2 데이터를 포함하지 않는 Kepler 영역만을 분석했을 때, 도출된 occurrence-$Z_{max}$ 기울기는 합동 샘플과 유사했다. $f$의 계단식 증가가 이 Kepler 전용 경향을 약간 더 잘 맞추었으나, 가파른 정도를 재현하는 근본적인 어려움은 여전히 남았다.
• 연령 정밀도의 역할: 실제 샘플과 합성 샘플 간의 비교는 모델이 경향을 맞추지 못하는 이유가 주로 항성 연령 측정의 불충분한 정밀도 때문이 아님을 나타낸다. 이상적인 연령을 사용하더라도 "점진적 증가" 모델은 배제되었으며, "계단" 모델은 비현실적인 매개변수를 요구했다.

의의 및 주장
본 논문은 소형 행성의 발생률과 은하 고도 사이의 강력한 음의 상관관계가 행성을 보유한 별의 분율이 시간에 따라 증가한다는 사실만으로는 설명될 수 없다고 결론짓는다. 저자들은 후기의 $f$ 부스트(특히 2~3 Gyr 전의 계단 함수)가 경향을 미약하게나마 재현할 수 있다고 언급하지만, 이는 모델의 물리적 한계(예: 오늘날 행성을 보유한 별이 거의 100%인 경우)까지 튜닝을 요구한다.

저자들은 이러한 불일치가 다른 요인들이 작용하고 있음을 시사한다고 주장한다. 그들은 향후 모델들이 다음을 통해 더 나은 적합도를 달erm 수 있을 것이라고 제안한다:

1. 슈퍼-지구와 서브-해 Neptune을 구분하는 것 (Z23이 이 경향이 주로 서브-해 Neptune에 의해 주도됨을 발견했듯이).
2. "온전한(intact)" 조밀한 다행성계의 분율을 일정하게 유지하는 대신 시간에 따라 진화하도록 허용하는 것.
3. 행성 형성을 강화하거나 오래된 시스템을 점진적으로 파괴하는 물리적 메커니즘(예: 항성 플라이바이, 은하 조석, 또는 성간 물질(ISM) 난류의 변화)을 조사하는 것.

이 연구는 국소적인 금속 함량 농축이 은하 맥락과 행성 인구 통계를 연결하는 유일한 고리라는 전통적인 관점이 Gaia 시대에는 불충분하다는 점을 강조한다. 결과는 행성 인구 통계학을 완전히 설명하기 위해 은하 규모의 역학적 역사(병합, 위성 유입 등)를 원시행성 원반 및 시스템 안정성에 영향을 미치는 국소적 과정과 연결할 필요가 있음을 보여준다.