[論文レビュー] Accelerated Lambda Iteration in Rapidly Expanding Envelopes

P. Hoêflich|arXiv (Cornell University)|Jul 4, 2002

Gamma-ray bursts and supernovae被引用数 3

ひとこと要約

本稿では、超新星の核崩壊や熱核爆発、新星、WR星における急速に拡張する低密度のエンベロープをモデル化するため、Hydra放射流体力学コード内に高速化されたラムダ反復（ALI）法を提示する。主な遷移を優先するための「リーディングエレメント」、ネットレート、レベルロック、予測補正スキームを用いることで、放射性崩壊と時間依存効果が支配的な非局所平衡（non-LTE）状態においても安定的かつ高速に収束する。この手法により、テストケースとしてのSN 99byにおける光曲線、フラックススペクトル、偏光スペクトルの正確な合成が可能になる。

ABSTRACT

We discuss the current implementation of the ALI method into our HYDrodynamical RAdiation(HYDRA) code for rapidly expanding, low density envelopes commonly found in core collapse and thermonuclear supernovae, novae and WR stars. Due to the low densities, non-thermal excitation by high energy photons (e.g. by radioactive decays) and the time dependence of the problem, large departures from LTE are common throughout the envelope even at large optical depths. ALI is instrumental for both the coupling of the statistical equations and the hydrodynamical equations with the radiation transport (RT). We employ several concepts to improve the stability, and convergence rate/control including the concept of leading elements, the use of net rates, level locking, reconstruction of global photon redistribution functions, equivalent-2-level approach, and predictive corrector methods. For appropriate conditions, the solution of the time-dependent rate equations can be reduced to the time-independent problem plus an analytic solution of an ODE For the 3-D problem, we solve the radiation transport via the moment equations. To construct the Eddington tensor elements, we use a Monte Carlo scheme to determine the deviation of the solution of the RT equation from the diffusion approximation (ALI of second kind). At the example of a thermonuclear supernova (SN99by),we show an analysis of light light curves, flux and polarization spectra and discuss the limitations of our approach.

研究の動機と目的

• 急速に拡張する低密度の天体的エンベロープにおける非局所平衡（non-LTE）放射輸送のラムダ反復の安定性と収束速度を向上させること。
• 特にSNe Iaおよび低光度イベントにおいて、放射性崩壊によって駆動される時間依存的非局所平衡効果を正確にモデル化できること。
• 高度なALI技術を用いて、放射輸送、原子準位分布、流体力学を反復的に結合し、数値的効率を向上させること。
• 原子準位構造を簡略化し、周波数分解能を制限した条件下でも、3次元モデルへの放射流体力学シミュレーションの適用範囲を拡大すること。
• SN 99byの詳細なスペクトルおよび光曲線解析を通じて、手法を検証し、遅い時期の挙動を捉える際の限界を評価すること。

提案手法

• 時間依存的かつ急速に拡張するエンベロープにおける非局所平衡統計平衡方程式を解くために、加速ラムダ反復（ALI）法を実装する。
• 源関数の更新において主要な遷移を優先するための「リーディングエレメント」を用い、収束性の向上と数値的不安定性の低減を図る。
• ネットレートとレベルロックを用いて、原子準位分布と放射輸送の反復的解法を安定化する。
• 全光子再分配関数の再構築と、等価2準位法を用いることで、レート方程式の解法を単純化する。
• 予測補正法を適用し、1階線形常微分方程式（ODE）を解析的に解くことで、時間依存レート方程式を扱い、完全な時間ステッピングに依存するのを減らす。
• 3次元放射輸送における拡散近似からの逸脱を補正するため、モンテカルロ法を用いてエディントンテンソル要素を計算する。

実験結果

リサーチクエスチョン

• RQ1 非局所平衡（non-LTE）、時間依存的放射流体力学シミュレーションにおける急速に拡張する超新星エンベロープにおいて、ラムダ反復をどのように高速化・安定化できるか。
• RQ2 時間依存的非局所平衡問題が、時間に依存しない解に加えて、準位分布の時間発展を記述する解析的ODE解にどれほど還元可能か。
• RQ3 低密度・高速度流におけるALI法の主な数値的不安定要因は何か。また、リーディングエレメントとレート補正を用いることで、それらをどのように緩和できるか。
• RQ4 本手法は、SN 99byのような低光度熱核超新星において、観測された光曲線、フラックススペクトル、偏光をどれほど正確に再現できるか。
• RQ5 特に禁制線寄与の影響により、遅い時期の進化において、この手法にどのような限界があるか。

主な発見

• リーディングエレメント、ネットレート、予測補正を組み合わせたALI法は、急速に拡張するエンベロープにおける非局所平衡計算において、安定的かつ高速に収束することが確認された。
• 束縛自由遷移が支配的である場合には、時間依存レート方程式を1階線形ODEとして解析的に解くことができ、完全な時間ステッピングに比べて計算コストを削減できる。
• SN 99byの例では、観測された光曲線、フラックススペクトル、偏光がモデルにより正確に再現され、手法の予測能力が妥当であることが検証された。
• 低密度と放射性崩壊による高エネルギー光子の吸収により、高光学厚さであっても強力な非局所平衡効果が再現された。
• 遅い時期の光曲線において、モデルに反映されていない禁制線寄与の影響により限界が生じており、禁制遷移の取り扱いの改善が求められる。
• 本手法により、原子モデルを簡略化し、周波数分解能を制限した条件下でも3次元放射流体力学が可能となり、複雑な天体的対象のシミュレーションが現実可能になった。