[論文レビュー] A calibration of the mixing-length for solar-type stars based on hydrodynamical simulations - I. Methodical aspects and results for solar metallicity

H.‐G. Ludwig, B. Freytag+1 more|arXiv (Cornell University)|Nov 11, 1998

Stellar, planetary, and galactic studies参考文献 1被引用数 35

ひとこと要約

この論文は、太陽型星の対流表面層を対象とした2次元放射移動流体力学（RHD）シミュレーションを用いて、混合長パラメータ（α_MLT）を補正する。α_MLTの値はF型矮星で約1.3、K準巨星で1.75であり、太陽近辺で平坦なピークを示す。太陽振動学的データとの整合性を確認し、内部精度±0.05を達成し、標準的な星形成モデルへの応用を想定したフィッティング式を提供する。

ABSTRACT

Based on detailed 2D numerical radiation hydrodynamics (RHD) calculations of time-dependent compressible convection, we have studied the dynamics and thermal structure of the convective surface layers of solar-type stars. The RHD models provide information about the convective efficiency in the superadiabatic region at the top of convective envelopes and predict the asymptotic value of the entropy of the deep, adiabatically stratified layers. This information is translated into an effective mixing-length parameter alpha suitable to construct standard stellar structure models. We validate the approach by a detailed comparison to helioseismic data. The grid of RHD models for solar metallicity comprises 58 simulation runs with a helium abundance of Y=0.28 in the range of effective temperatures 4300K < Teff < 7100K and gravities 2.54 < logg < 4.74. We find a moderate, nevertheless significant variation of alpha between about 1.3 for F-dwarfs and 1.75 for K-subgiants with a dominant dependence on Teff. In the close neighbourhood of the Sun we find a plateau where alpha remains almost constant. The internal accuracy of the calibration of alpha is estimated to be +/- 0.05 with a possible systematic bias towards lower values. An analogous calibration of the convection theory of Canuto & Mazzitelli gives a different temperature dependence but a similar variation of the free parameter. For the first time, values for the gravity-darkening exponent beta are derived independently of mixing-length theory: beta=0.07...0.10. We show that our findings are consistent with constraints from stellar stability considerations and provide compact fitting formulae for the calibrations.

研究の動機と目的

• 直接的な流体力学的シミュレーションを用いて、ヒューリスティックな仮定に依存しない太陽型星の混合長パラメータ（α_MLT）を補正すること。
• 特に太陽近辺を含むHertzsprung-Russell図全体におけるα_MLTの変動を定量化し、星形成構造モデルの精度を向上させること。
• シミュレーションに基づく補正が太陽振動学的データと整合することを検証し、観測された太陽構造と一貫性を持つことを確認すること。
• 混合長理論に依存せずに重力減光指数βを独立して導出することにより、新たな物理的制約を提供すること。
• 標準的な星形成モデルに適用可能な、α_MLTおよびα_CMT（Canuto-Mazzitelli理論）の簡潔で実用的なフィッティング式を提供すること。

提案手法

• 太陽金属量（Y=0.28）を仮定した、58回の時間依存的かつ圧縮性の2次元放射移動流体力学（RHD）シミュレーションを実施。対象は4300 K ≤ T_eff ≤ 7100 K、2.54 ≤ log g ≤ 4.74の太陽型星をカバーする。
• RHDシミュレーションから、断熱的層化が成立する深部のエントロピー（s_env）を抽出し、状態方程式を用いて熱力学的変数を関連付ける。
• シミュレートされた対流効率と標準的な混合長理論との一致をもとに、s_envを有効な混合長パラメータα_MLTに変換する。
• 同様の手順をCanuto & Mazzitelli（1991, 1992）の対流理論（CMT）に適用し、比較のためのα_CMTを導出する。
• 数値的に使いやすい簡潔なフィッティング関数を導出するため、正規化された温度（T̃ = (T_eff - 5770)/1000）および重力（g̃ = log(g/27500)）を用いる。
• 太陽構造の予測値（s_env、α_MLT）を太陽振動学的制約と比較することで、補正の妥当性を検証し、小さな不確実性の範囲内で一貫性があることを確認する。

実験結果

リサーチクエスチョン

• RQ1 太陽金属量を有する太陽型星において、Hertzsprung-Russell図全体にわたる混合長パラメータα_MLTの変動はどのように変化するか？
• RQ2 RHDシミュレーションは、太陽振動学的制約に基づく太陽構造を正確に再現できるか？ これにより、シミュレーション手法の妥当性が検証されるか？
• RQ3 α_MLTは効用温度および表面重力にどのように依存するか？ また、標準的な混合長理論の仮定と比較するとどうなるか？
• RQ4 RHDシミュレーションから、混合長理論に依存せずに重力減光指数βを独立して導出できるか？
• RQ5 RHDシミュレーションの結果を、標準的な星形成モデルに容易に統合可能な形で簡潔にパラメータ化できるか？

主な発見

• 混合長パラメータα_MLTは、F型矮星で約1.3、K準巨星で1.75にまで変動し、主に効用温度に依存する。
• 太陽近辺（T_eff ≈ 5770 K）では、α_MLTはほぼ一定の平坦なピークを示し、約1.587の値を示す。これは太陽振動学的制約と整合的である。
• α_MLT補正の内部精度は±0.05と推定され、やや低い値への系統的バイアスの可能性がある。
• 重力減光指数βは独立してβ = 0.07〜0.10として導出され、星の安定性の観点からも整合的である。
• Canuto & Mazzitelli（CMT）対流理論の補正では、温度依存性は標準混合長理論とは異なるが、自由パラメータの変動は類似している。
• 正規化された温度および重力変数を用いた簡潔なフィッティング式がα_MLTおよびα_CMTについて提供され、星形成モデルへの直接統合が可能である。