[论文解读] Mixing-length calibration from field stars: An investigation on statistical errors, systematic biases, and spurious metallicity trends

G. Valle, M. Dell’Omodarme+2 more|arXiv (Cornell University)|Jan 1, 2019

Stellar, planetary, and galactic studies参考文献 63被引用 17

一句话总结

本研究利用理论模拟数据集，调查了在实际恒星中混合长度参数（αml）校准的可靠性。研究发现，即使在理想条件下，αml估计值仍表现出较大的随机离散度和系统性偏差，且由于重元素混合物、氦-金属丰度比（∆Y/∆Z）以及[Fe/H]标度的不匹配，会引致虚假的金属丰度趋势，从而对αml与金属丰度之间存在本征依赖关系的主张提出质疑。

ABSTRACT

Aims. We critically analysed the theoretical foundation and statistical reliability of the mixing-length calibration by means of standard (Teff, [Fe/H]) and global asteroseismic observables (Δν, νmax) of field stars. We also discussed the soundness of inferring a possible metallicity dependence of the mixing-length parameter from field stars.
 
 Methods. We followed a theoretical approach based on mock datasets of artificial stars sampled from a grid of stellar models with a fixed mixing-length parameter αml. We then recovered the mixing-length parameter of the mock stars by means of SCEPtER maximum-likelihood algorithm. We finally analysed the differences between the true and recovered mixing-length values quantifying the random errors due to the observational uncertainties and the biases due to possible discrepancies in the chemical composition and input physics between artificial stars and the models adopted in the recovery.
 
 Results. We verified that the αml estimates are affected by a huge spread, even in the ideal configuration of perfect agreement between the mock data and the recovery grid of models. While the artificial stars were computed at fixed solar-calibrated αml = 2.10, the recovered values had a mean of 2.20 and a standard deviation of 0.52. Then we explored the case in which the solar heavy-element mixture used to compute the models is different from that adopted in the artificial stars. We found an estimated mixing-length mean of 2.24 ± 0.48 and, more interestingly, a metallicity relationship in which αml increases by 0.4 for an increase of 1 dex in [Fe/H]. Thus, a simple heavy-element mixture mismatch induced a spurious, but statistically robust, dependence of the estimated mixing-length on metallicity. The origin of this trend was further investigated considering the differences in the initial helium abundance Y – [Fe/H] – initial metallicity Z relation assumed in the models and data. We found that a discrepancy between the adopted helium-to-metal enrichment ratio ΔY/ΔZ caused the appearance of spurious trends in the estimated mixing-length values. An underestimation of its value from ΔY/ΔZ = 2.0 in the mock data to ΔY/ΔZ = 1.0 in the recovery grid resulted in an increasing trend, while the opposite behaviour occurred for an equivalent overestimation. A similar effect was caused by an offset in the [Fe/H] to global metallicity Z conversion. A systematic overestimation of [Fe/H] by 0.1 dex in the recovery grid of models forced an increasing trend of αml versus [Fe/H] of about 0.2 per dex. We also explored the impact of some possible discrepancies between the adopted input physics in the recovery grid of models and mock data. We observed an induced trend with metallicity of about Δαml = 0.3 per dex when the effect of the microscopic diffusion is neglected in the recovery grid, while no trends originated from a wrong assumption on the effective temperature scale by ±100 K. Finally, we proved that the impact of different assumptions on the outer boundary conditions was apparent only in the RGB phase.
 
 Conclusions. We showed that the mixing-length estimates of field stars are affected by a huge spread even in an ideal case in which the stellar models used to estimate αml are exactly the same models as used to build the mock dataset. Moreover, we proved that there are many assumptions adopted in the stellar models used in the calibration that can induce spurious trend of the estimated αml with [Fe/H]. Therefore, any attempt to calibrate the mixing-length parameter by means of Teff, [Fe/H], Δν, and νmax of field stars seems to be statistically poorly reliable. As such, any claim about the possible dependence of the mixing-length on the metallicity for field stars should be considered cautiously and critically.

研究动机与目标

• 本文研究了利用场星观测量校准混合长度参数（αml）时的统计误差与系统性误差。
• 研究旨在检验αml估计值中观测到的金属丰度趋势是否真实，还是模型与观测不匹配所致的人为效应。
• 本研究旨在量化化学成分、∆Y/∆Z以及输入物理假设的差异对αml恢复的影响。
• 评估使用标准（Teff, [Fe/H]）和全局星震学（Δν, νmax）观测量进行αml校准的稳健性。
• 目标是批判性评估关于场星中αml与金属丰度依赖关系的主张的有效性。

提出的方法

• 作者从一组固定αml = 2.10的恒星模型网格中生成了模拟的人工恒星数据集。
• 他们应用SCEPtER最大似然算法，从这些模拟观测中恢复αml值。
• 恢复过程使用与模拟数据在化学成分、∆Y/∆Z或输入物理方面可能存在差异的恒星模型网格。
• 他们系统性地改变模型假设，如重元素混合物、∆Y/∆Z、[Fe/H]标度偏移、微观扩散和边界条件，以隔离误差来源。
• 通过分析真实αml值与恢复αml值之间的差异，量化随机误差与系统性偏差。
• 研究重点在于这些不匹配如何在αml与[Fe/H]的关系中诱导虚假趋势，特别是在场星星震学校准的背景下。

实验结果

研究问题

• RQ1 观测数据中的统计不确定性在多大程度上影响了场星中αml估计的可靠性？
• RQ2 模型假设中的系统性偏差（如化学成分或∆Y/∆Z）是否会导致αml与[Fe/H]之间出现虚假趋势？
• RQ3 当模拟数据与恢复模型之间的[Fe/H]丰度标度存在不匹配时，会对αml校准产生何种影响？
• RQ4 假设的∆Y/∆Z比值在生成人工αml-金属丰度趋势中起到何种作用？
• RQ5 输入物理假设的差异（如微观扩散、有效温度标度）如何影响恢复的αml值及其与金属丰度的依赖关系？

主要发现

• 在理想情况下（模型与数据完全一致），恢复的αml平均值为2.20，标准差为0.52，尽管真实值固定为2.10。
• 重元素混合物不匹配（如GS98与AS09之间）导致αml出现每[Fe/H]每十年增加0.4的虚假趋势，平均估计值为2.24±0.48。
• 错误的∆Y/∆Z比值（如假设为1.0而非2.0）会在αml与[Fe/H]之间产生显著的上升趋势，而相反假设则导致下降趋势。
• [Fe/H]标度存在0.1 dex的系统性偏移，导致αml出现约每[Fe/H]每十年0.2的虚假趋势。
• 在恢复网格中忽略微观扩散，会引致约每[Fe/H]每十年0.3的虚假αml趋势。
• 有效温度标度存在±100 K的误差，会在αml中产生刚性偏移，但不会在[Fe/H]上产生系统性趋势。