[論文レビュー] A New Atmospheric Model for HD 189733 b

Jonathan Langton, Gregory Laughlin|ArXiv.org|Aug 22, 2008

Stellar, planetary, and galactic studies参考文献 16被引用数 3

ひとこと要約

本稿では、HD 189733 bの大気のための新しい2次元、静水的平衡の流体力学モデルを提示する。このモデルは、超音速で乱流的な風を再現し、混沌とした温度変動を引き起こす。モデルは、観測された8-μmの光曲線を成功裏に再現し、予期しない30°東方への熱帯のオフセットと寒冷帯の西寄りオフセットを、定常的流れではなく大規模な大気乱流に起因すると説明している。

ABSTRACT

We have developed a new two-dimensional hydrostatically-balanced isobaric hydrodynamic model for use in simulation of exoplanetary atmospheres. We apply this model to the infrared photosphere of the hot Jupiter HD 189733 b, for which an excellent 8-micron light curve has been obtained. For reasonable parameter choices, the results of our model are consistent with these observations. In our simulations, strongly turbulent supersonic flow develops, with wind speeds of approximately 5 km/s. This flow geometry causes chaotic variation of the temperature distribution, leading to observable variations in the light curve from one orbit to the next.

研究の動機と目的

• より物理的に整合性のある2次元の熱木星の大気モデルを開発すること。このモデルは静水的平衡を維持し、従来の浅水域モデルや圧縮性モデルの限界を克服すること。
• HD 189733 bの観測された8-μm位相曲線を説明すること。特に、予想に反して最も高温の領域が亜点から30°東方にあること、および最も寒冷な領域が反亜点から30°西寄りにあるという異常なオフセットを説明すること。
• 大規模な大気乱流がフラックス極値の軌道変動を説明できるかどうかを検証すること。観測された光曲線の形状に動的な説明を提供すること。
• 検証可能な仮説を提示すること：乱流が実際に存在する場合、フラックスの最小値と最大値は複数の軌道間で数時間ずれることになる。これは将来的な観測で検証可能である。

提案手法

• モデルは2次元、等圧、静水的平衡の流体力学フレームワークを採用しており、垂直方向の静水的平衡を保ちながら、高解像度で水平方向の力学を解明している。
• 完全圧縮性流体力学の手法を用い、静水的平衡を強制する修正された方程式系を採用することで、従来の圧縮性モデルが示す物理的に不自然な特徴を回避している。
• 調整可能なパラメータを用いて放射力（放射強制）を組み込み、星からの放射照射を模擬しているが、完全な多波長放射移動は含んでいない。
• 初期条件として強い温度勾配を設定し、強い放射照射によって駆動されることで、超音速で乱流的な風のパターンが発生する。
• モデルは時間に依存する温度場および風場を出力し、それらを用いて可視円盤全体のフラックスを統合することで、8 μmでの合成光曲線を計算している。
• 合成光曲線は、Knutsonら（2007年）の観測された8-μmフラックス曲線と直接比較され、フィットの質と軌道変動の程度が評価されている。

実験結果

リサーチクエスチョン

• RQ1 2次元、静水的平衡の流体力学モデルは、非定常的フラックス変動を含め、HD 189733 bの観測された8-μm位相曲線を再現できるか？
• RQ2 なぜHD 189733 bの最も高温の領域が亜点から30°東方にあるのか、また最も寒冷な領域が反亜点から30°西寄りにあるのか。これは、対称的でゼロナルフローの予想とは反する。
• RQ3 定常的風のパターンではなく、大規模な大気乱流が、観測されたフラックス極値の軌道変動を説明できるのか？
• RQ4 モデルが予測する時間に依存するフラックス最小値と最大値（複数の軌道間で数時間ずれる）は、将来的な観測で検証可能か？
• RQ5 8 μmと24 μmでの観測されたフラックス曲線の特徴の違いは、異なる気圧レベルでの大気力学的要因によるものか、それとも乱流変動によるものか、その程度はどの程度か？

主な発見

• モデルは約5 km s⁻¹の超音速で乱流的な風速を生成し、この風が惑星の昼側および夜側に混沌とした非定常的温度分布を引き起こしている。
• シミュレートされた光曲線は、フラックス極値における軌道変動を示しており、最小値と最大値のタイミングが複数の軌道間で数時間ずれることを示している。これは観測の不確実性（±1時間）と整合的である。
• モデルは、亜点から30°東方へのフラックス最大値のオフセットと、反亜点から30°西寄りへのフラックス最小値のオフセットを、観測結果とよく一致させている。
• 現在の観測ではモデルが除外されない。放射強制パラメータの制限はあるが、8-μm光曲線へのフィットは定性的には妥当である。
• モデルは、HD 189733 bの8-μm位相曲線における異常な寒冷帯オフセットを、対称的ゼロナルフローではなく乱流輸送に起因すると説明する最初の説明を提供している。
• モデルは、将来的な観測で複数の軌道にわたってフラックスタイミングをモニタリングすることで、乱流に起因するシフトを検証すべきであると示唆しており、仮説の直接的な観測的検証を提供する。