[論文レビュー] Annihilation Cross Sections for Relic Densities in the Low Velocity Limit

James D. Wells|arXiv (Cornell University)|Apr 6, 1994

Computational Fluid Dynamics and Aerodynamics被引用数 15

ひとこと要約

本稿では、初期宇宙における質量のあるコールドレリクスの消失断面積の低速度展開を計算する簡素化された手法を提示する。相対速度 $v_{\rm rel}$ の平方振幅の導関数を用いて、$\sigma v_{\rm rel} = a_{\rm rel} + b_{\rm rel}v_{\rm rel}^2 + c_{\rm rel}v_{\rm rel}^4$ の係数 $a_{\rm rel}$, $b_{\rm rel}$, $c_{\rm rel}$ を抽出する。零速度極限における角度依存性の不在を活用することで、数値積分を回避し、熱平均断面積を高速かつ高精度に計算可能となり、誤差リスクが最小限に抑えられる。この手法は、ディラックニュートリノおよびニュートリノの消失断面積に対する既知の結果と整合しており、妥当性が確認されている。

ABSTRACT

A simple method to calculate the total annihilation cross section of massive cold relics is presented. Exact expressions for $\arel$, $\brel$, and $\crel$ of $σ\vrel=\arel+\brel \vrel^2+\crel \vrel^4$ are obtained from simple derivatives of the squared annihilation amplitude. Three examples are given. In one example, the computed annihilation cross section of massive Dirac neutrinos confirms recent correct calculations.

研究の動機と目的

• 初期宇宙における質量のあるコールドレリクスの熱平均消失断面積の計算を簡素化すること。
• 出版された消失断面積計算における頻発する不一致を解消するため、体系的かつ誤差に強くない手法を提供すること。
• 低速度展開 $\sigma v_{\rm rel} = a_{\rm rel} + b_{\rm rel}v_{\rm rel}^2 + c_{\rm rel}v_{\rm rel}^4$ の $a_{\rm rel}$, $b_{\rm rel}$, $c_{\rm rel}$ 係数を効率的に計算できること。
• ディラックニュートリノおよびニュートリノの消失断面積に対する既知の結果と照合し、手法の妥当性を検証すること。

提案手法

• 零速度極限における相対速度 $v_{\rm rel}$ のテイラー展開を用いて、微分断面積を $v_{\rm rel}$ のべき級数に展開する。
• Mandelstam 変数 $s$ および $t$ が $v_{\rm rel} = 0$ で角度依存性を示さないという事実を活用し、微分断面積を $v_{\rm rel}^2$ の級数に簡略化し、三角関数係数を含む形に変換する。
• 振幅の二乗 $|M|^2$ の正確な導関数を用いて、振幅展開の係数 $J_0$, $J_2$, $J_4$ などを計算する。
• 各 $v_{\rm rel}^2$ のべき乗に対して $\cos\theta$ について解析的積分を実行し、複雑な位相空間積分の数値積分を回避する。
• $K(s)$、すなわち位相空間因子を $v_{\rm rel}^2$ についてテイラー展開することで $K_0$, $K_2$, $K_4$ を得る。
• $J(s,t)/4\pi$ と $K(s)$ の展開を組み合わせ、立体角について統合することで、直接的に $\sigma v_{\rm rel}$ の係数を算出する。

実験結果

リサーチクエスチョン

• RQ1 数値積分を回避しつつ、消失断面積 $\sigma v_{\rm rel}$ の低速度展開を効率的かつ正確に計算する方法は何か？
• RQ2 振幅の二乗 $|M|^2$ の導関数を用いて、$a_{\rm rel}$, $b_{\rm rel}$, $c_{\rm rel}$ を体系的に抽出する手法は何か？
• RQ3 本手法は、ディラックニュートリノおよびニュートリノの消失断面積に対して、既存の計算とどのように比較できるか？
• RQ4 本手法は、ダークマター候補のレリック密度計算における計算作業と誤差率を低下させることができるか？

主な発見

• 本手法は、$|M|^2$ の導関数から直接 $a_{\rm rel}$, $b_{\rm rel}$, $c_{\rm rel}$ を計算可能であり、$\sigma v_{\rm rel}$ の数値積分を一切不要にした。
• ディラックニュートリノの場合、得られた消失断面積は Dolgov と Rothstein [9] の結果と一致し、正しさが確認され、以前の研究 [6,7] との不一致も解消された。
• ニュートリノのZボソン交換によるsチャネル消失の場合、$a_x$ は Ellis ら [8] の結果と一致するが、$b_x$ は $3/2 a_x$ の因子違いとなる。これは、熱平均化の処理方法の違いに起因する。
• 記号計算ソフトウェアと組み合わせることで、計算時間は数秒にまで短縮され、レリック密度研究における信頼性と効率性が著しく向上した。
• 狭い共鳴がある場合でも本手法は頑健であるが、低速度展開が破綻する場合には注意が必要である。