Impact of neutrino-electron scattering and an improved treatment of pair processes on binary neutron star mergers

技術要約：中性微子・電子散乱および改良された対生成プロセスの連星中性子星合体への影響

問題提起
マルチメッセンジャー観測による連星中性子星（BNS）合体は、高密度物質の性質や重元素核合成に関する独自の制約を与える。しかし、現在のシミュレーションにおける中性微子物理学の正確性は、結果として生じる電磁波信号（キロノバ）の信頼性の高いモデリングを制限している。最先端のコードは通常、エネルギー積分型の（「グレイ」）2モーメント・スキームを採用している。これらは効率的ではあるが、中性微子の圧力テンソルおよび、極めて重要な中性微子のエネルギースペクトルに対して近似的な解析的クローザーに依存している。これらの近似は、全分布関数に敏感な反応、例えば電子に対する非弾性散乱（$\nu + e^- \to \nu + e^-$）や中性微子・反中性微子の対消滅（$\nu + \bar{\nu} \leftrightarrow e^+ + e^-$）の相互作用率を正確に計算することを妨げる。さらに、膨大な計算コストをかけずに、グリッドセル内で中性微子の分布関数 $f_\nu$ を捉えることは、歴史的に角度依存性を無視するか、大きな統計的（ショット）ノイズを受け入れるかの選択を強いてきた。

手法
著者らは、SpECコード内での改良されたモンテカルロ（MC）輸送アルゴリズムを用いた一連のBNS合体シミュレーションを提示している。核心となる手法の進展は以下の通りである：

1. 適応型パケット重み付け（Adaptive Packet Weighting）: 反応率のために $f_\nu$ を計算できるように、著者らはMCパケットの重み付けを変更した。定数エネルギーまたは定数個数のパケットを用いる代わりに、パケットが利用可能な位相空間の一定の割合（$\Delta f_\nu$）を表すスキームを実装した。これにより、半透明領域および光学的に薄い領域での $f_\nu$ の推定が可能となった。このスキームは、以下の3つの処方箋の間で動的に遷移する：

   • ショットノイズを抑えるために、高エネルギーの中性微子に対しては定数エネルギーパケットを使用。
   • 平衡仮定が成立する光学的に厚い領域では、定数個数パケットを使用。
   • 関連するエネルギー領域（O(10 MeV)）において、分布関数を解像するために定数 $\Delta f_\nu$ パケットを使用。
     このアプローチにより、流体への結合におけるショットノイズの増加という代償はあるものの、単一のパケットが角度積分された分布関数の変化 $\Delta f_\nu < 1$ を表す状態に到達することが可能となった。

2. 改良された中性微子・物質結合: 中性微子パケットと流体の間の結合は、「期待値」法（低密度領域でのノイズ抑制に使用）と、運動量転送イベントの明示的な計上（保存則を確保するために高密度領域で使用）の加重平均によって処理される。

3. 新しい反応カーネルの導入:

   • 非弾性散乱: シミュレーションには、電子に対する中性微子の非弾性散測が含まれている。合体条件下での核に対する弾性散乱と比較して、このプロセスが副次的であるため、著者らは等方的な中性微子分布を仮定し、計算をルジャンドル展開の0次項に簡略化している。
   • 対プロセス: $\nu\bar{\nu} \leftrightarrow e^+e^-$ の扱いは、中性微子のエネルギー分布を仮定するのではなく、シミュレートされたエネルギー分布を利用することで改善された。角度依存性は、応力エネルギーテンソル・モーメントから導出された局所的な「フォームファクター」を用いて近似され、グリッドセル内の平均的な伝搬方向に対して等方的なレートを補正している。

4. シミュレーション設定: 著者らは、全質量（$M_{tot} \approx 2.7M_\odot$ から $2.9M_\odot$）は異なるが、質量比と状態方程式（SFHo）は一定である3つの連星構成をシミュレートした。これにより、即時崩壊（$\sim 4-6$ ms以内）を起こす系と、長寿命の残骸（$> 10$ ms）を形成する系を比較することができる。一つの構成（$M127-M118$）は、ベースモデル、新しいカーネルを持たない適応型重み付けモデル、およびフル版の高度なモデル（Adapt + Kernels）の3つの異なる中性微子処理を用いて発展させた。

主要な結果

• 分布関数の解像度: 適応型重み付けスキームは、種や時刻に応じて $\Delta f_\nu < 1$（$\sim 0.01$ から $0.6$の範囲）を、計算コストを劇的に増大させることなく達成することに成功した。これにより、ブロック因子や相互作用率を計算するために、シミュレートされた$f_\nu$ を直接使用することが可能となった。これは、これまでのグローバルな合体シミュレーションでは到達不可能であった能力である。
• 中性微子の特性への影響: 非弾性散乱と改良された対プロセスの導入は、重レプトン中性微子（$\nu_x$）の集団を著しく変化させる。
  • 不透明度: 非弾性散乱により、脱離領域および高温の潮汐腕における $\nu_x$ に対する不透明度は2〜5倍に増加する。
  • エネルギーと光度: $\nu_x$ の平均エネルギーは約10%減少し、残骸が高度なシミュレーションにおいてわずかに高温であるにもかかわらず、その総光度はベースモデルと比較して約30%低下した。
  • スペクトル形状: 重レプトン中性微子は、散乱が支配的な領域において、$T_\nu \ll T_{fluid}$ および正のスペクトル指数を特徴とする、高度に非熱的な分布を示す。
• 物質放出への影響: 改良された中性微子物理学は、$M127-M118$ 構成において、合体後最初の10 msの間に放出される非束縛物質の総質量を50%増加させる。これは、アウトフローが発生する領域における重レプトン中性微子によるエネルギー沈着の強化に起因する。放出物の組成も変化し、$Y_e \sim 0.2$ 付近の物質の顕著な増加が見られる。
• 連星質量依存性: 連星の総質量は、アウトフローの特性を強く決定する。高質量系（$M_{tot} \gtrsim 2.7M_\odot$）は、主にコアバウンス振動に関連した高温・高 $Y_e$（$> 0.25$）のアウトフローを生成する。低質量系は、スパイラルアームのダイナミクスによって駆動されると思われる、中間的な $Y_e$（$\sim 0.2$）の放出物を主成分として生成するが、これが中性微子物理学に最も敏感な成分である。
• 対消滅: 合体初期相（$t < 10$ ms）において、極領域における対消滅は、極領域が高温・高密度な物質（$T \sim 6.5$ MeV, $\rho \sim 10^9$ g/cm$^3$）で満たされているため、アウトフローを駆動する上で副次的である。なぜなら、逆反応（$e^+e^- \to \nu\bar{\nu}$）が順方向のプロセスを支配するか、あるいは均衡させているためである。

意義と主張
本論文は、非弾性散乱を含み、かつ中性微子のシミュレートされたエネルギー分布を考慮した対消滅の扱いを含む、初の合体シミュレーションを提示していると主張している。著者らは、合体の定性的な進化（例：崩壊時間、一般的なアウトフローの幾何学）は先行研究と同様であるが、放出物の定量的な詳細は著しく変化することを強調している。

主な意義は、グローバルなシミュレーション内で $\Delta f_\nu < 1$ を用いて中性微子の分布関数を捉えることが可能であることを示した点にあり、これにより、より正確な反応率の計算が可能となる。著者らは、中性微子の物理学の改良は、モンテカルロ法と2モーメント・スキームのどちらを選択するかよりも、初期の時間的なアウトフローの質量と組成に対してより実質的な影響を与えるが、磁場（ここではモデル化されていない）の導入ほどの影響はないと結論付けている。結果は、平均的なアウトフロー特性に頼ることは、核合成収量における決定的な違いを覆い隠してしまう可能性があることを示唆している。詳細な $Y_e$ の分布は、連星の質量および中性微子の物理的処理によって大きく変化するためである。論文は控えめな姿勢を維持しており、現在の $\Delta f_\nu$ 値は高精度な反応率計算のための有用性の限界付近にあり、ショットノイズがこの解像度のトレードオフであると述べている。