[論文レビュー] 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を有する星のモデルグリッドから人工星のモックデータセットを生成する。
• 彼らは、これらのモック観測からαml値を回復するためにSCEPtERの最尤法を適用する。
• 回復プロセスでは、モックデータと異なる化学組成、∆Y/∆Z、または入力物理を有する星のモデルグリッドを用いる。
• 重元素混合物、∆Y/∆Z、[Fe/H]スケールオフセット、微視的拡散、境界条件といったモデル仮定を体系的に変化させ、誤差要因を特定する。
• 真の値と回復されたαml値の差を分析し、ランダム誤差と系統的バイアスを定量化する。
• 本研究は、これらの不一致がαml vs. [Fe/H] における誤った傾向を引き起こすメカニズム、特に場の星の太陽系外惑星キャリブレーションの文脈において焦点を当てる。

実験結果

リサーチクエスチョン

• RQ1 観測データにおける統計的不確実性が、場の星からのαml推定値の信頼性にどの程度影響を及えるか？
• RQ2 化学組成や∆Y/∆Zといったモデル仮定の系統的バイアスが、[Fe/H] におけるαmlの誤った傾向を生じるか？
• RQ3 モックデータと回復モデル間の[Fe/H]濃度スケールにおける不一致が、αmlキャリブレーションにどのような影響を及えるか？
• RQ4 仮定された∆Y/∆Z比が人工的なαml-金属量傾向を生成する役割を果たすか？
• RQ5 入力物理（例：微視的拡散、効果的温度スケール）の不一致が、回復されたαml値およびその金属量依存性に与える影響は何か？

主な発見

• 完全なモデル-データ一致の理想状態においても、回復されたαmlの平均は2.20、標準偏差は0.52であったが、真の値は固定された2.10であった。
• 重元素混合物の不一致（例：GS98対AS09）により、[Fe/H] 1デキサドごとにαmlが0.4上昇する誤った傾向が生じ、平均推定値は2.24±0.48となった。
• 誤った∆Y/∆Z比（例：2.0ではなく1.0を仮定）は、[Fe/H] とともに増加するαml傾向を生じさせたが、逆に仮定した場合、減少傾向が生じた。
• [Fe/H]スケールに0.1デキサドの系統的オフセットがあると、[Fe/H] 1デキサドあたり約0.2の誤ったαml傾向が生じた。
• 回復グリッドで微視的拡散を無視すると、[Fe/H] 1デキサドあたり約0.3の誤ったαml傾向が生じた。
• 効果的温度スケールに±100 Kの誤差があると、αmlに剛体オフセットが生じたが、[Fe/H] との系統的傾向は生じなかった。