[論文レビュー] 2.5D global-disk oscillation models of the Be shell star ζ Tauri I. Spectroscopic and polarimetric analysis

C. Escolano, A. C. Carciofi+4 more|arXiv (Cornell University)|Mar 2, 2015

Astrophysics and Star Formation Studies参考文献 21被引用数 4

ひとこと要約

本研究では、Be星ζ Tauriの分光的および偏光的変動をよりよく再現するために、2.5次元のグローバルディスク振動モデルを導入した。水平方向の構造を組み込みつつも、等張的ダイナミクスを保持することで、HαおよびBr 15 V/R変動の一致が向上したが、特に高い粘性係数（α = 0.8）を用いた場合に顕著であった。一方で、シェル線プロファイルの変動や偏光振幅の再現には依然課題が残る。

ABSTRACT

Context. Various Be stars exhibit intensity variations of the violet and red emission peaks in their HI lines observed in emission. This so-called $V/R$ phenomenon is usually explained by the precession of a one-armed spiral density perturbation in the circumstellar disk. That global-disk oscillation scenario was confirmed, both observationally and theoretically, in the previous series of two papers analyzing the Be shell star ζ Tauri. The vertically averaged (2D) global-disk oscillation model used at the time was able to reproduce the $V/R$ variations observed in Hα, as well as the spatially resolved interferometric data from AMBER/VLTI. Unfortunately, that model failed to reproduce the $V/R$ phase of Br15 and the amplitude of the polarization variation, suggesting that the inner disk structure predicted by the model was incorrect. Aims. The first aim of the present paper is to quantify the temporal variations of the shell-line characteristics of ζ Tauri. The second aim is to better understand the physics underlying the $V/R$ phenomenon by modeling the shell-line variations together with the $V/R$ and polarimetric variations. The third aim is to test a new 2.5D disk oscillation model, which solves the set of equations that describe the 3D perturbed disk structure but keeps only the equatorial (i.e., 2D) component of the solution. This approximation was adopted to allow comparisons with the previous 2D model, and as a first step toward a future 3D model. Results. The new 2.5D formalism improves the agreement with the observed $V/R$ variations of Hα and Br15, under the proviso that a large value of the viscosity parameter, α = 0.8, be adopted. Nonetheless, it remains challenging for the models to reproduce consistently the amplitude and the average level of the polarization data, whatever formalism is adopted.

研究の動機と目的

• ζ Tauriのシェル線特性（ピーク間隔、吸収深さ、ドップラー速度シフト）の時間的変動を定量化すること。
• 分光的・偏光的・干渉計的データを同時にモデル化することで、V/R現象の物理的起源を解明すること。
• 水平方向の構造を捉えつつも等張的ダイナミクスを保持する新しい2.5次元形式をテストし、今後の3次元モデル化への一歩とすること。
• 従来の2次元モデルの限界、特にBr 15 V/R位相と偏光振幅の再現性の不足を評価すること。
• 非軸対称ディスク構造の探査に向けたシェル線変動の診断的潜在能力を同定すること。

提案手法

• Hα、Hβ、Hγ、Hδ、Br 15線の広範な分光的解析を実施し、発光ピーク強度、シェル吸収深さ、非対称性、および径方向ドップラー速度を測定した。
• HDUST放射線輸送コードに2.5次元グローバルディスク振動モデルを実装し、3次元の摂動ディスク方程式を解いたが、解のうち等張的（2次元）成分のみを保持した。
• 2.5次元モデルを用いて線プロファイルおよび偏光をシミュレートし、AMBER/VLTIおよび他の機器による複数エポック分の分光的・干渉計的観測データと比較した。
• 粘性係数（α）、ディスクスケール高さ、基盤密度、および傾斜角といった主要パラメータを変化させ、観測量に与える影響を評価した。
• 2.5次元モデルの予測を2次元モデルと観測結果と比較し、V/Rサイクルおよび偏光の再現性向上を評価した。
• HαおよびHβプロファイルにおける三重ピーク位相を分析し、連星による駆動の摂動とその関連を検討した。

実験結果

リサーチクエスチョン

• RQ1 シェル線特性（吸収深さ、ドップラー速度、非対称性）はV/Rサイクルと位相を同じくして変動するか？また、それらはグローバルディスク振動の診断として有効か？
• RQ2 2.5次元形式は、従来の2次元モデルに比べてBr 15のV/R位相をよりよく再現できるか？
• RQ3 なぜモデルは観測された偏光変動の振幅を再現できないのか？また、ディスクスケール高さはその要因としてどのような役割を果たすか？
• RQ4 HαおよびHβにおける三重ピーク位相の原因は何か？また、それらの出現は連星伴と関連しているか？
• RQ5 粘性係数αの選択が、らせん状密度波の閉じ込め方および観測結果との一致度にどのように影響するか？

主な発見

• HαおよびBr 15のピーク間隔は、V/Rサイクルと同期したサイン波に類似した変動を示し、グローバルディスク振動の診断としての有効性が裏付けられた。
• シェル吸収深さにはわずかな変動が認められ、V/R最大付近でやや浅くなる傾向が見られたが、データ不足のため明確な結論は得られなかった。
• シェル線の非対称性は時間的にほとんど一定であり、モデル予測の大きな変動とは矛盾しており、これはモデル化の主要な限界を示している。
• 2.5次元モデルは、高い粘性係数（α = 0.8）を用いることで、Br 15 V/R位相および連続スペクトル偏光の観測とより一致するよう改善された。これは2次元モデルの欠陥を是正した。
• 2.5次元モデルおよび2次元モデルの両方とも、ピーク間隔の変動振幅を過小評価しており、モデル内での摂動の閉じ込めが不十分であることが示唆された。
• 2.5次元形式は2次元モデルに比べてより閉じたらせん状密度波を生成するが、観測と一致させるにはより高いα値を必要とし、他のBe星に関する最近の知見とも整合的である。