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

技術要約：銀河系における惑星主星比率の増加は、惑星出現率における銀河系高度の傾向を再現できない

問題提起
近年の観測結果、具体的には Zink et al. (2023) (以下 Z23) は、FGK型主系列星における小型惑星（スーパー・アースおよびサブ・ネプチューン）の出現率が、銀河面からの最大垂直振動振幅 ($Z_{max}$) の増加に伴って著しく減少するという強い傾向を特定した。恒星の金属量は惑星形成と相関することが知られているが、銀河系の金属量勾配だけでは、観察された急激な減少を説明するには不十分である。本研究が取り組む中心的な問いは、天の川銀河の歴史における時間依存的な惑星主星比率 ($f$) の増加が、観察された出現率-$Z_{max}$ の傾向を再現できるかどうかである。著者らは、過去のある時点における惑星形成効率の「ブースト」と、時間の経過に伴う恒星軌道の動的加熱（キネマティック・ヒーティング）を組み合わせることで、銀河系高度における惑星の不足を説明できるかどうかを調査している。

手法
著者らは、FGK型主系列星のサンプルに対して合成惑星集団をシミュレートするために、新しい惑星系集団合成コード psps を開発した。手法には以下のステップが含まれる：

1. 恒星サンプル: 年齢の不確実性に対する結果の堅牢性を検証するため、2つの異なる恒星サンプルを作成した：

   • 実サンプル: Hardegree-Ullman et al. (2025) のカタログから抽出され、Berger et al. (2020b, 2025) のアイソクロネ年齢とクロスマッチングされた、Kepler および K2 の結合サンプル（22,184個の星）。$Z_{max}$ 値は、Gaia DR3 の運動学と Gala パッケージを用いた銀河ポテンシャルモデルを用いて計算された。
   • 合成サンプル: Kepler および K2 の視野に似せるよう調整された TRILEGAL 集団合成サンプル（49,109個の星）。これは、年齢精度が任意に高く、既知の動的加熱モデルを持つ「楽観的」なシナリオを提供し、年齢精度による影響を物理モデルから分離することを可能にした。

2. 時間依存的な出現モデル: 恒星の年齢（ルックバックタイム）の関数としての惑星主星比率 ($f$) の進化に関する2つのモデル群をテストした：

   • ステップ関数: 特定のルックバックタイム $t$ において、初期の比率 $f_1$ から現在の比率 $f_2$ へと急激に変化するモデル。
   • 区分関数: 時刻 $t$ から $f_1$ から $f_2$ へと緩やかに線形に増加するモデル。
   • すべてのモデルは、現在の惑星主星比率が観測的制約と一致する $\sim 30-35\%$ となるように制約されている。

3. 集団合成と観測: psps を用いて、規定された $f$ に基づく惑星系を生成した。惑星系には、既存の人口統計学的制約（例：Lam & Ballard 2024）に基づいたアーキテクチャ（コンパクトな多惑星系 vs 単一惑星系）および軌道パラメータ（離心率、傾斜角）を割り当てた。生成された合成集団は、幾何学的なトランジット確率と Kepler および K2 ミッションの特定の感度関数を考慮した逆検出効率法（IDEM）を用いて「観測」された。

4. 解析: シミュレーションされた出現率は $Z_{max}$ (100–1000 pc) ごとにビン分割され、べき乗則 ($y \propto Z_{max}^\tau$) でフィッティングされ、傾き ($\tau$) と正規化定数 ($\eta$) が決定された。これらの結果は、Z23 による観測された傾向 ($\tau \approx -0.28$) と比較された。

主な結果

• 主星比率の増加の不十分さ: 本研究では、惑星主星比率 ($f$) を時間の経過とともに増加させても（瞬間的であれ漸進的であれ）、観察された出現率-$Z_{max}$ の傾向の急峻さを再現するには一般的に不十分であることが示された。
• モデルの制約:
  • 実際の Kepler-K2 サンプルの場合: ほとんどのモデルは、観測された $\tau \approx -0.28$ と比較して、傾き ($\tau$) が緩やかすぎる（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 視野のみを解析した場合、導出された出現率-$Z_{max}$ の傾きは結合サンプルと同様であった。$f$ のステップ増加は、この Kepler 単独の傾向にはわずかに良く一致したが、急峻さを再現するという根本的な困難は残った。
• 年齢精度の役割: 実サンプルと合成サンプルの比較は、モデルが傾向に一致しない主な原因が、恒星年齢測定の精度不足によるものではないことを示している。理想的な年齢を用いた場合でも、「漸進的な増加」モデルは否定され、「ステップ」モデルは非現実的なパラメータを必要とした。

意義と主張
本論文は、小惑星の出現率と銀河系高度の間の強い負の相関は、惑星を保有する恒星の比率の時間依存的な増加のみでは説明できないと結論付けている。著者らは、最近の $f$ のブースト（具体的には 2–3 Gyr 前のステップ関数）は、モデルの物理的限界（例：現在、恒星のほぼ 100% が惑星を保有している状態）へのチューニングを行えば、傾向をわずかに再現できる可能性があると述べている。

著者らは、この不一致は他の要因が作用していることを示唆していると主張している。彼らは、将来のモデルがより良い適合を得るためには、以下の点が必要であると考えている：

1. スーパー・アースとサブ・ネプチューンの区別を行うこと（Z23 が、この傾向は主にサブ・ネプチューンによって駆動されていることを発見した通り）。
2. 「完全な」コンパクトな多惑星系の割合を、一定に保つのではなく、時間の経過とともに進化させること。
3. 惑星形成を強化するか、あるいは古い惑星系を段階的に破壊する物理的メカニズム（例：恒星のフライバイ、銀河潮汐、または ISM 乱流の変化）を調査すること。

本研究は、局所的な金属濃縮を、銀河論的文脈と惑星デモグラフィクスを結びつける唯一のリンクとする標準的な図式は、ガイア（Gaia）時代の観測においては不十分であることを強調している。結果は、惑星のデモグラフィクスを完全に説明するためには、銀河規模の動力学的履歴（合体、衛星銀河の流入）と、原始惑星系円盤や系の安定性に影響を与える局所的なプロセスとを結びつける必要があることを浮き彫りにしている。