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

技术摘要：银河系行星宿主比例的增加无法重现行星发生率的银河系高度趋势

问题陈述
最近的观测结果，特别是 Zink 等人 (2023)（以下简称 Z23），确定了一个强烈的趋势：在围绕 FGK 型矮星的小型行星（超级地球和亚海王星）发生率随最大垂直振荡振幅（$Z_{max}$）增加而显著下降。虽然恒星金属度已知与行星形成相关，但仅靠银河系金属度梯度不足以解释观察到的这种剧烈下降。本研究旨在解决的核心问题是：在银河系历史上，一个随时间变化的行星宿主比例（$f$）的增加，是否能够重现观察到的 occurrence-$Z_{max}$ 趋势。作者研究了在银河系过去的某个时刻，行星形成效率的“提升”，结合随时间变化的恒星轨道动力学加热效应，是否可以解释高银河系高度处行星数量的缺失。

方法论
作者开发了一种新型行星系统群体合成代码 psps，用于模拟围绕 FGK 矮星的合成行星群体。其方法步骤如下：

1. 恒星样本： 构建了两个不同的恒星样本，以测试结果对年龄不确定性的稳健性：

   • 真实样本： 一个联合 Kepler 和 K2 样本（22,184 颗恒星），提取自 Hardegree-Ullman 等人 (2025) 的目录，并与来自 Berger 等人 (2020b, 2025) 的等龄线年龄进行交叉匹配。使用 Gaia DR3 运动学和通过 Gala 包构建的银河系势场模型计算了 $Z_{max}$ 值。
   • 合成样本： 一个经过调优以模拟 Kepler 和 K2 场特征的 TRILEGAL 群体合成样本（49,109 颗恒星）。该样本提供了一个“乐观”的情景，具有任意精确的恒星年龄和已知的动力学加热模型，使作者能够分离年龄精度与物理模型本身的影响。

2. 随时间变化的发生率模型： 研究测试了两类关于行星宿主比例（$f$）随恒星年龄（回溯时间）演化的模型：

   • 阶跃函数： 在特定的回溯时间 $t$，从初始比例 $f_1$ 突然跳跃到现今比例 $f_2$。
   • 分段函数： 从时间 $t$ 开始，从 $f_1$ 到 $f_2$ 进行逐渐的线性增加。
   • 所有模型都受到约束，以产生与观测约束一致的约 30-35% 的现今行星宿主比例。

3. 群体合成与观测： 使用 psps，作者根据预设的 $f$ 生成了行星系统。系统被分配了架构（紧凑的多行星系统 vs 单行星系统）以及轨道参数（离心率、倾角），这些参数基于现有的统计特征（例如 Lam & Ballard 2024）。随后，利用逆向检测效率方法 (IDEM) 对合成群体进行“观测”，该方法考虑了几何凌日概率以及 Kepler 和 K2 任务的具体灵敏度函数。

4. 分析： 将模拟的发生率按 $Z_{max}$（100–1000 pc）进行分箱，并使用幂律函数（$y \propto Z_{max}^\tau$）进行拟合，以确定斜率（$\tau$）和归一化系数（\eta}）。将这些结果与 Z23 的观测趋势（$\tau \approx -0.28$）进行对比。

主要结果

• 宿主比例增加的局限性： 研究发现，随着时间的推移增加行星宿主比例（$f$），无论是瞬时增加还是逐渐增加，通常不足以重现观察到的 occurrence-$Z_{max}$ 趋势的陡峭程度。
• 模型约束：
  • 对于真实的 Kepler-K2 样本，绝大多数模型产生的斜率（$\tau$）都过于平缓（接近 0），与观察到的 $\tau \approx -0.28$ 不符。最接近匹配该趋势的模型需要极端的、在物理上不切实际的参数（例如，$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 的分析： 当仅分析 Kepler 场（其年龄分布不同且缺乏 K2 数据）时，推导出的 occurrence-$Z_{max}$ 斜率与联合样本相似。一个 $f$ 的阶跃增加比联合样本稍好地匹配了 Kepler 仅有的趋势，但重现其陡峭度的根本困难依然存在。
• 年龄精度的作用： 真实样本与合成样本的比较表明，模型无法匹配该趋势的主要原因并非由于恒星年龄测量的精度不足。即使在理想化的年龄条件下，“逐渐增加”的模型也被排除在外，而“阶跃”模型则需要不切实际的参数。

意义与结论
论文得出结论，观察到的小行星发生率与银河系高度之间的强负相关关系，不能仅由随时间变化的行星宿主比例来解释。作者指出，虽然一个晚期的 $f$ 提升（特别是发生在 2–3 Gyr 前的阶跃函数）可以勉强重现该趋势，但这需要对模型的物理极限进行调优（例如，现今几乎 100% 的恒星都拥有行星）。

作者认为，这种差异意味着必然有其他因素在起作用。他们建议，未来的模型可以通过以下方式实现更好的拟合：

1. 区分超级地球和亚海王星（正如 Z23 所发现的，该趋势主要由亚海王星驱动）。
2. 允许“完整”的紧凑多行星系统的比例随时间演化，而非保持不变。
3. 调查可能增强近期行星形成或逐渐破坏较老系统的物理机制（例如，通过恒星飞掠、银河系潮汐或星际介质湍流的变化）。

该研究强调，在 Gaia 时代，将局部金属富集作为连接银河系背景与行星统计特征的唯一纽带的传统观点已不再充分。研究结果凸显了将银河系尺度的动力学历史（合并、卫星入侵）与影响原行星盘和系统稳定性的局部过程联系起来的必要性，以全面解释系外行星的人口统计特征。