[論文レビュー] 2D non-LTE Modeling for Axi-symmetric Winds. II. A Short Characteristic Solution for Radiative Transfer in Rotating Winds

J. Zsargó, D. J. Hillier+1 more|ArXiv.org|Oct 22, 2005

Astrophysics and Star Formation Studies参考文献 37被引用数 9

ひとこと要約

本論文は、非単調な速度場と非局所熱平衡（non-LTE）状態を正確に扱える、回転対称な星風に対する新しい短特徴線放射移動コードを提示する。インパクトパラメータベクトルと特徴線に沿った適応的解像度を用いることで、高精度かつ分散処理可能な放射移動計算が可能となり、最小限の誤差で2次元回転モデルにおける期待される線幅特徴を正確に再現した。

ABSTRACT

We present a new radiative transfer code for axi-symmetric stellar atmospheres and compare test results against 1D and 2D models with and without velocity fields. The code uses the short characteristic method with modifications to handle axi-symmetric and non-monotonic 3D wind velocities, and allows for distributed calculations. The formal solution along a characteristic is evaluated with a resolution that is proportional to the velocity gradient along the characteristic. This allows us to accurately map the variation of the opacities and emissivities as a function of frequency and spatial coordinates, but avoids unnecessary work in low velocity regions. We represent a characteristic with an impact-parameter vector p (a vector that is normal to the plane containing the characteristic and the origin) rather than the traditional unit vector in the direction of the ray. The code calculates the incoming intensities for the characteristics by a single latitudinal interpolation without any further interpolation in the radiation angles. Using this representation also provides a venue for distributed calculations since the radiative transfer can be done independently for each p.

研究の動機と目的

• 複雑で非単調な風の流れを有する回転対称星の大気における、強固な放射移動手法の開発。
• 1次元および球対称モデルの限界を克服し、回転する大質量星、Be星、および連星系のシミュレーションを可能にする。
• 非局所熱平衡（non-LTE）モデルと分散処理をサポートする、柔軟で高精度な短特徴線手法の実装。
• 周波数依存の吸収係数と発光度の適切なマッピングにより、回転星風の正確なスペクトル合成を可能にする。
• 1次元コードと既知の2次元風の特徴との比較を通じて、今後の自己自己一致型シミュレーションにおける信頼性を保証する。

提案手法

• インパクトパラメータベクトルを方向ベクトルの代わりに用いる独自の表現を採用した短特徴線アプローチを採用。これにより、各レイトの独立計算が可能となる。
• 入射放射度の計算に1回の緯度方向補間のみを実行し、角度補間の必要性を排除することで、計算効率が向上する。
• 各特徴線に沿った解像度を、速度勾配に基づいて適応的に調整することで、必要な箇所での周波数マッピングの精度を確保する。
• 速度変動に応じて空間サンプリングを動的に調整することで、非単調な風の流れに対しても対応可能となる。
• インパクトパラメータと軸傾き角のペア、または周波数範囲に沿って作業を分割することで、分散処理に最適化された設計となっている。
• ASTAROTHコードにおける従来の連続体移動とSobolev近似を置き換え、2次元モデルにおける完全な非局所熱平衡（non-LTE）線移動を実現している。

実験結果

リサーチクエスチョン

• RQ1 非単調な速度場を有する回転対称星風において、どのように放射移動を高精度にモデル化できるか？
• RQ2 非球面的かつ非局所熱平衡状態の風構造を扱うにあたり、角度補間を回避できるように短特徴線手法をどのように適応可能にするか？
• RQ3 適応的解像度スキームは、高速度勾配領域における精度と計算効率をどのように向上させるか？
• RQ4 新しい手法は、2次元回転モデルにおける既知のスペクトル特徴をどの程度正確に再現できるか？
• RQ5 大規模な星の大気シミュレーションにおいて、分散処理がどのように効率的に並列化可能か？

主な発見

• 純散乱テストにおいて、CMFGENの基準1次元結果と比較して数パーセントの誤差で正確に再現した。
• 2次元回転モデルにおいて、青方シフトおよび歪んだ線幅特徴が正確に捉えられ、特にβ = π/4およびπ/2のような中間視点角で顕著に現れた。
• 光球面の角度サイズと端部の暗黙の効果を適切に扱ったため、He I 4713.17 Å線における観測された線中心の青方シフトとブルーサイドの傾きが正しく再現された。
• 光球面における前後非対称性が、後方側で広がったプロファイルと前方側で狭くなったプロファイルを生じさせ、結果として観測されたブルーサイドの傾きが誘発された。
• 散乱における周波数再分配を含んでいなくても、合成スペクトルは既存の報告結果と定性的に一致し、本手法の信頼性を確認した。
• 本実装により、今後のASTAROTHにおける自己自己一致型2次元非局所熱平衡（non-LTE）シミュレーションが可能となり、回転および非対称星風のモデリングが著しく向上した。