[论文解读] Allesfitter: Flexible Star and Exoplanet Inference From Photometry and Radial Velocity

Maximilian N. Günther, Tansu Daylan|arXiv (Cornell University)|Mar 31, 2020

Stellar, planetary, and galactic studies参考文献 130被引用 99

一句话总结

allesfitter 是一个灵活的开源 Python 软件包，用于从光度和径向速度数据中联合推断恒星与系外行星，支持掩食时间变体、相位曲线、恒星变异性、黑子、耀斑以及通过高斯过程建模的系统性噪声。它支持使用马尔可夫链蒙特卡洛（MCMC）和嵌套采样进行参数估计与贝叶斯模型比较，可生成出版级输出，并提供用户友好的图形界面或编程接口，适用于编程用户。

ABSTRACT

We present allesfitter, a public and open-source python software for flexible and robust inference of stars and exoplanets given photometric and radial velocity data. Allesfitter offers a rich selection of orbital and transit/eclipse models, accommodating multiple exoplanets, multi-star systems, transit-timing variations, phase curves, stellar variability, star spots, stellar flares, and various systematic noise models including Gaussian Processes. It features both parameter estimation and Bayesian model selection, allowing to easily run either a Markov Chain Monte Carlo (MCMC) or Nested Sampling fit. For novice users, a graphical user interface allows to specify all input and perform analyses; for python users, all modules can be readily imported in any existing scripts. Allesfitter also produces publication-ready tables, LaTeX commands, and figures. The software is publicly available (https://github.com/MNGuenther/allesfitter), pip-installable (pip install allesfitter) and well documented (www.allesfitter.com). Finally, we demonstrate the software's capabilities on several examples and provide updates to the literature where possible for Pi Mensae, TOI-216, WASP-18, KOI-1003 and GJ 1243.

研究动机与目标

• 为应对从多任务光度和径向速度数据中对系外行星与恒星系统建模的可扩展性、可重复性和自动化推断日益增长的需求。
• 克服现有软件通常仅关注特定问题（如仅掩食或径向速度）且缺乏稳健贝叶斯模型比较能力的局限性。
• 将多种前向模型（如掩食、食变、相位曲线、耀斑和噪声模型）统一到一个可扩展的框架中。
• 为新手用户提供图形界面，为高级用户提供可编程 API，确保广泛可及性并可无缝集成到现有工作流中。
• 提供出版级输出，包括 LaTeX 表格、图表和命令，提升可重复性和传播性。

提出的方法

• 通过贝叶斯定理使用贝叶斯推断计算给定数据和先验分布下模型参数的后验概率。
• 采用假设噪声为高斯分布的似然函数，对光度和径向速度数据进行建模，对数似然定义为卡方项之和。
• 支持使用 MCMC（通过 emcee）和嵌套采样（通过 dynesty）进行后验抽样，后者可实现稳健的贝叶斯证据估计。
• 集成来自成熟软件包的模块化前向模型：使用 ellc 进行掩食/食变和径向速度建模，使用 aflare 进行恒星耀斑建模，使用 celerite 进行高斯过程建模，以及其他模型。
• 通过多项式、样条或高斯过程对系统性噪声进行建模，支持多历元和多源数据的灵活配置。
• 支持注入-恢复测试，并可处理多行星、多恒星及包含掩食时间变体、相位曲线和黑子建模的复杂系统。

实验结果

研究问题

• RQ1 一个单一、可扩展的软件框架能否稳健地处理从光度和径向速度数据中联合推断多种系外行星与恒星现象？
• RQ2 allesfitter 在参数估计和贝叶斯模型比较方面，对包含复杂噪声和变异性现象的多样化系统，其性能如何？
• RQ3 allesfitter 的输出在多大程度上与独立分析在真实系外行星系统（如 TOI-216、WASP-18、KOI-1003 和 GJ 1243）中发表的结果一致？
• RQ4 allesfitter 在建模恒星耀斑、黑子和相位曲线时，其统计严谨性和计算效率表现如何？
• RQ5 allesfitter 的图形界面和 API 能否同时满足新手用户和专家开发者的需要，实现统一、可重复的工作流？

主要发现

• allesfitter 成功复现了 TOI-216 的已发表结果，其轨道参数和掩食时间变体（TTVs）的后验估计与 Dawson 等人（2019）的结果高度一致。
• 在 WASP-18 系统中，allesfitter 恢复的光学相位曲线振幅（A1 = 0.0221 ± 0.0045，B1 = 0.1443 ± 0.0056，B2 = 0.1922 ± 0.0059）与 Shporer 等人（2019）的结果高度一致。
• 在 KOI-1003 系统中，allesfitter 对双星系统的后验估计（如 RA/RB = 0.1634 ± 0.0021，P = 8.35992 ± 0.00094 天）与 Roettenbacher 等人（2016）的结果一致，微小偏差归因于黑子穿越效应。
• 在 GJ 1243 系统中，allesfitter 在 1.8 天的 TESS 观测段中成功建模了两个黑子和三个耀斑，恢复的耀斑峰值时间（如 t_peak,flare1 = 2456519.028417 ± 0.000038 BJD）和振幅（A_flare1 = 0.0973 ± 0.0034）与预期一致。
• 该软件可生成出版级输出，包括 LaTeX 表格和图表，且用户操作极少，所有测试案例均得到验证。
• allesfitter 的嵌套采样实现提供了低方差的贝叶斯证据估计，支持统计上稳健的模型比较，而仅使用标准 MCMC 无法实现这一点。