[论文解读] FUNDPAR: A Program for deriving Fundamental Parameters from Equivalent Widths

C. Saffe|arXiv (Cornell University)|Sep 14, 2010

Stellar, planetary, and galactic studies参考文献 9被引用 7

一句话总结

FUNDPAR 是一个 FORTRAN 程序，通过强制执行电离平衡和激发能平衡，利用铁线等效宽度推导类太阳恒星的基本恒星参数（Teff, log g, [Fe/H], ξ）。它通过 χ² 最小化过程迭代调整参数，结合 ATLAS 模型大气和 MOOG 元素丰度合成，成功确定了 58 颢系系外行星候选恒星的基本参数，结果与文献值一致。

ABSTRACT

We implemented a fortran code that determine fundamental parameters of solar type stars from a list of Fe line equivalent widths. The solution should verify 3 conditions in the standard method: ionization equilibrium, excitation equilibrium and independence between metallicity and equivalent widths. We added the condition that the input metallicity of the model atmosphere should be similar to the output metallicity derived with equivalent widths. Solar-scaled Kurucz model atmospheres with NEWODF opacities are calculated with an independent program. Parameter files control different details, such as the mixing-length parameter, the overshooting, the damping of the lines and the weight factors in the definition of the chi2 function. FUNDPAR derive the uncertainties following 2 methods: the criteria of Gonzalez & Vanture (1998) and the dispersion using the chi2 function. The code use the 2009 version of the MOOG program. The results derived with FUNDPAR are in agreement with previous determinations in the literature. In particular we obtained the fundamental parameters of 58 exoplanet host stars. The program is freely available from the web (http://icate-conicet.gob.ar/saffe/fundpar/).

研究动机与目标

• 开发一种稳健、自动化的工具，用于从类太阳恒星的铁线等效宽度推导基本恒星参数（Teff, log g, [Fe/H], ξ）。
• 强制执行三项标准平衡条件：电离平衡（Fe I/Fe II）、激发能独立性，以及等效宽度的金属丰度独立性。
• 确保输入模型大气的金属丰度与从等效宽度推导出的金属丰度保持一致。
• 使用 Gonzalez & Vanture (1998) 方法和 χ² 方差分析两种方法提供不确定度估计。
• 将该方法应用于 58 颗主序系外行星候选恒星样本，确认其为富金属恒星。

提出的方法

• 该程序在 Teff、log g、[Fe/H] 和微湍流速度（ξ）四个维度上进行迭代最小化，初始值默认为 5000 K、4.00 dex、-0.10 dex、1.00 km/s。
• 对每个参数组合，通过 atlas.launcher 接口调用 ATLAS9 生成 Kurucz 模型大气，从 Kurucz 网格中选取最接近的模型并进行优化。
• 通过子程序 'kurucz2moog' 将模型大气转换为 MOOG 可读格式，以确保与 MOOG 元素丰度分析代码的兼容性。
• 调用 MOOG 的 'abfind' 驱动程序，从观测到的等效宽度计算 Fe I 和 Fe II 丰度，返回推导出的 [Fe/H] 值。
• χ² 函数计算为三项残差的加权和：电离平衡（c₁²）、激发能平衡（c₂²）和金属丰度独立性（c₃²），权重由用户定义（w₁, w₂, w₃）。
• 算法通过下降法迭代调整参数，直至 χ² 值收敛，确保满足全部三项平衡条件。

实验结果

研究问题

• RQ1 是否能够通过完全自动化的 FORTRAN 代码，从等效宽度可靠地推导恒星基本参数，同时强制执行电离平衡和激发能平衡？
• RQ2 与已有研究中已知的系外行星候选恒星相比，推导出的参数在准确性上如何？
• RQ3 模型大气中输入与输出金属丰度之间的一致性在多大程度上提升了参数的可靠性？
• RQ4 推导出的参数的不确定度是多少？与 Gonzalez & Vanture (1998) 方法和 χ² 方差估计相比有何差异？
• RQ5 该方法是否能成功识别出如文献中先前报道的系外行星候选恒星的富金属趋势？

主要发现

• FUNDPAR 成功推导出 58 颗主序系外行星候选恒星的基本参数，确认其为富金属恒星，与先前文献结果一致。
• 通过强制执行电离平衡、激发能平衡和金属丰度独立性，程序在全部 58 颗恒星上实现收敛，χ² 函数作为主要收敛判据。
• 通过两种方法估计不确定度：Gonzalez & Vanture (1998) 方法和 χ² 函数的方差分析，两种方法均得到一致且可靠的误差估计。
• 对 58 颗恒星推导出的参数（Teff, log g, [Fe/H], ξ）与文献中先前的测定结果高度一致，验证了代码的准确性。
• 该代码在处理广泛多样的恒星参数方面表现出稳健性，包括金属丰度范围从 [Fe/H] ≈ -0.54 dex 到 +0.43 dex 的恒星。
• FUNDPAR 可免费获取：http://icate-conicet.gob.ar/saffe/fundpar/，并提供完整文档、安装指南以及输入/输出格式说明。