Light Travel Time Effects in Kilonova Models

問題の背景
GW170817/AT2017gfoに代表されるキロノバは、極めて急速なスペクトル進化を示し、その特徴は数時間という時間スケールで変化する。この急速な変化は、より緩やかに進化する超新星などの天体とは対照的である。特定の観測上のパズルとして、AT2017gfoにおけるストロンチウム（Sr II）のP Cygniプロファイルが挙げられる。すなわち、その放出成分は、付随する吸収ディップよりも後に現れ、かつ、より長く持続するというものである。時間依存的な放射伝達コードは存在するものの、光の伝播時間の効果（リバーブレーション／残響）がこれらの特徴を形成する上で果たす具体的な役割については、これまでの研究では分離・特定されてこなかった。さらに、一部の研究では、有限の光速による遅延を無視し、光の到着を瞬時と仮定する時間独立的なコードを使用している（キロノバの膨張速度は最大で $\sim 0.3c - 0.7c$ に達する）。著者らは、光の伝播時間の遅延が、キロノバのスペクトル特徴の形状、進化、および持続性にどのように寄与しているかを定量化することを目的とし、これらの遅延のみによって、AT2017gfoで観測された長寿命の放出と赤方偏移へのドリフトを説明できるかどうかを検証する。

手法
著者らは、光子パケットの光球から遠方の観測者への旅路を追跡するために、簡略化されたモンテカルロ放射伝達（MCRT）モデルを用いる。

• 幾何学構造: 線形成領域は、二つの同心球（連続的に膨張する $v = r/t_{ph}$）の間の体積としてモデル化されている。内側の球はブラックボディ光球を表し、外側の球は噴出物の境界を定義する。
• 時間シーケンシング: 二つの時間変数、$t_{ph}$（パケットが光球で生成された合体後の時間）と $t_{det}$（パケットが観測者に到着した時間）が利用される。関係式 $t_{det} = t_{ph} + d/c$ は、経路長 $d$ が、放出位置（近傍側か遠方側か）や散乱の履歴に基づいて変化することを考慮している。
• 物理: モデルは、ソボレフ近似を用いた、Sr IIマルチプレット（静止波長 10,500 Å）を近似する単一の散乱線を対象としている。純粋な散乱（蛍光なし）を想定し、束縛・束縛相互作用を追跡する。
• パラメータ: シミュレーションは、AT2017gfoのX-shooterによる初期観測に一致するように、合体後 $t_{det} = 1.43$ から $4.40$ 日間の範囲をカバーする。光球温度と速度は、観測データ（Gillanders et al. 2022; Smartt et al. 2017）の経験的なフィットから導出される。温度は 10,000 K から 3,200 K へと進化し、速度は $0.30c$ から $0.12c$ へと減少する。
• 比較: 本研究では、完全な時間依存性（光の伝播遅延を考慮）を用いて生成されたスペクトルを、「時間独立」なアプローチ（光速を無限大と仮定）と比較し、さらに、時間独立の結果に単純な平均時間遅延を適用した第三のアプローチと比較する。

主な結果

1. 顕著な時間遅延: 光球から同時に放出された光子は、観測者に最大 1.25 日間の広がりを持って到着する。初期段階（例：$t_{ph} = 1.43$ 日）では、散乱光子の飛行時間は合体後の時間の 10% 以上を占めることがある。
2. 急速なスペクトル進化: 同じ $t_{det}$ ウィンドウ内に到着するものの、異なる $t_{ph}$ 値に由来する光子から合成されたスペクトルは、明確に異なる特性を示す。わずか数時間の間であっても、光球温度の変化とドップラーシフトにより、連続成分のピークや特徴の形状は大きく変化する。
3. 放出の持続性: モデルは、P Cygниプロファイルの放出成分が吸収ディップよりも長く持続するという観測結果を再現することに成功した。散乱光子はより長い経路を辿るため、より遅れて到着し、赤方偏移するため、これが「リバーブレーション（残響）」効果として作用し、放出特徴を維持させる。
4. 赤方偏移へのドリフト: 散乱特徴のピークは、時間の経過とともに長波長側へとシフトする。これは、連続膨張（後退側から来る光子が、より青方偏移していないこと）と、より初期の高温な時期からの、より長い経路を辿った光子の支配力の増大という組み合わせによって駆動される。
5. リバーブレーションによる説明の限界: リバーブレーション効果は、放出特徴の持続性と赤方偏移へのドリフトを説明するが、モデルによれば、散乱のみでは AT2017gfo における Sr II 特徴の強度と寿命を完全に説明することはできない。モデルにおける散乱フラックスは、観測と比較して減衰が早すぎるためであり、これは他の物理プロセスが特徴の持続に寄与していることを示唆している。
6. 時間独立モデルへの影響: 時間独立モデルとの比較により、特に急速な温度進化が行われている期間（例：1.43 日時）において、重大な相違があることが明らかになった。時間独立モデルは、異なる熱的エポックからの光子の混合を考慮していないため、フラックスを過大評価し、より青い温度を予測してしまう。単純な平均時間遅延による補正は一致度を改善するが、急速な進化フェーズにおいては依然として不十分である。

意義と主張
本論文は、光の伝播時間の効果は、特に急速な膨張と進化する熱条件を伴うイベントにおいて、キロノバのモデリングにおいて無視できないものであると結論付けている。

• 時間依存性の必要性: 急速に進化するイベントにおいて、時間独立なアプローチは、視覚的に異なるだけでなく、物理的に不正確なスペクトル（特に、推定される光球温度や連続成分のフラックスレベルに関して）を生み出す可能性がある。
• リバーブレーションの役割: リバーブレーション効果（散乱光子の到着の遅れ）は、放出特徴の長寿命化と赤方偏移へのドリフトに寄与する要因であることが確認された。しかし、著者らは、この効果単独では AT2017gfo の全観測データを説明するには不十分であると控えめに主張しており、それは複雑な要因の相互作用の一要素であるとしている。
• 方法論的ガイダンス: 本研究は、時間独立なコードも引き続き使用可能であるが、特にスペクトルの急速な進化段階を分析する際には、その限界を認識しておく必要があると助言している。最も正確な解釈のためには、計算資源が許す場合には、完全な時間依存コードの使用が推奨される。本論文は、キロノバの全物理を解決すると主張するものではなく、観測データの誤解を避けるために、有限の光速を考慮することの決定的な重要性を強調するものである。