[論文レビュー] Neutrino transport in type II supernovae: Boltzmann solver vs. Monte Carlo method

Shoichi Yamada, Han-Thomas Janka+1 more|arXiv (Cornell University)|Sep 2, 1998

Neutrino Physics Research被引用数 7

ひとこと要約

この論文は、静的後発コアモデルを用いたタイプII超新星におけるニュートリノ輸送について、有限差分ボルツマン解法とモンテカルロ法を比較している。全体的には両手法の一致が得られるが、ボルツマン解法は、限られた角度分解能のため、極めて非等方的な領域でフラックスおよびエディントン因子を低く見積もっており、ホットバブル領域におけるニュートリノ加熱の過大評価のリスクがある。これは、精度を向上させるために、適応的角グリッドの導入が不可欠であることを示している。

ABSTRACT

We have coded a Boltzmann solver based on a finite difference scheme (S_N method) aiming at calculations of neutrino transport in type II supernovae. Close comparison between the Boltzmann solver and a Monte Carlo transport code has been made for realistic atmospheres of post bounce core models under the assumption of a static background. We have also investigated in detail the dependence of the results on the numbers of radial, angular, and energy grid points and the way to discretize the spatial advection term which is used in the Boltzmann solver. A general relativistic calculation has been done for one of the models. We find overall good agreement between the two methods. However, because of a relatively small number of angular grid points (which is inevitable due to limitations of the computation time) the Boltzmann solver tends to underestimate the flux factor and the Eddington factor outside the (mean) ``neutrinosphere'' where the angular distribution of the neutrinos becomes highly anisotropic. This fact suggests that one has to be cautious in applying the Boltzmann solver to a calculation of the neutrino heating in the hot-bubble region because it might tend to overestimate the local energy deposition rate. A comparison shows that this trend is opposite to the results obtained with a multi-group flux-limited diffusion approximation of neutrino transport. The accuracy of the Boltzmann solver can be considerably improved by using a variable angular mesh to increase the angular resolution in the semi-transparent regime.

研究の動機と目的

• タイプII超新星におけるニュートリノ輸送のための有限差分ボルツマン解法の開発と検証。
• 現実的な超新星条件下におけるボルツマン解法の精度をモンテカルロ輸送コードと比較。
• グリッド分解能（径方向、角度方向、エネルギー）および空間的移流離散化法がボルツマン解法の性能に与える影響の評価。
• ボルツマン解法の限界がホットバブル領域におけるニュートリノ加熱に与える影響の評価。
• 可変角メッシュを用いた改善策の検討（半透明領域における応用を想定）。

提案手法

• ボルツマン解法は、位相空間（位置、エネルギー、角度）におけるニュートリノ輸送方程式を有限差分S_N法で解く。
• ボルツマン解法は、固定および可変グリッド間隔を用いて角度およびエネルギー次元を離散化し、半透明領域では適応的角分解能を導入。
• 空間的移流は、異なる離散化スキームで処理され、その解の精度への影響をテストした。
• 解法は、超新星シミュレーションから得られた静的後発コアモデルに適用され、固定背景を仮定した。
• 一般相対性理論補正を1つのモデルに実装し、その輸送結果への影響を評価した。
• 物理的条件を同一にした状態で、結果を高精度モンテカルロ輸送コードとベンチマーク比較した。

実験結果

リサーチクエスチョン

• RQ1 現実的な超新星大気において、ボルツマン解法の結果はモンテカルロ法と定量的にどの程度一致するか？
• RQ2 角度グリッド点数の増減が、特にニュートリノの非等方性が顕著な領域におけるボルツマン解法の精度にどの程度影響を与えるか？
• RQ3 空間的移流離散化法の選択が、ボルツマン解法の解の安定性および精度に与える影響は何か？
• RQ4 可変角メッシュの使用が、半透明領域における解像度向上にどの程度寄与するか？
• RQ5 ボルツマン解法は、フラックスおよびエディントン因子の不正確さにより、ホットバブル領域における局所的エネルギー付加率を過大評価しているか？

主な発見

• ボルツマン解法は、超新星モデルの大部分の領域でモンテカルロ法と良好な一致を示した。
• 半透明領域および外側領域では、角度分解能が不十分なため、フラックス因子およびエディントン因子が低く見積もられた。
• この低く見積もられた結果は、非等方性が強いホットバブル領域における局所的ニュートリノ加熱率の過大評価を示唆する。
• この誤差の傾向は、多群フラックス制限拡散近似と逆であり、後者では加熱が過大評価されがちである。
• 可変角メッシュを用いることで、半透明領域における角度分解能が著しく向上し、フラックスおよびエディントン因子の誤差が低減された。
• 本研究は、ボルツマン解法における角度分解能が、超新星におけるニュートリノ輸送の精度にとって重要な制限要因であることを確認した。