[论文解读] 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

一句话总结

本文提出了一种简化方法，用于计算早期宇宙中大质量冷暗物质遗迹的湮灭截面的低速展开，通过幅值平方的导数提取系数 $a_{\rm rel}$、$b_{\rm rel}$ 和 $c_{\rm rel}$，形式为 $\sigma v_{\rm rel} = a_{\rm rel} + b_{\rm rel}v_{\rm rel}^2 + c_{\rm rel}v_{\rm rel}^4$。该方法通过利用零速极限下角度无关性，避免了数值积分，从而实现快速、精确的热平均截面计算，且出错风险极低，经验证与狄拉克中微子及中性子湮灭的已知结果一致。

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} = 0$ 时Mandelstam变量 $s$ 和 $t$ 无角度依赖性的事实，将微分截面简化为仅含 $v_{\rm rel}^2$ 的级数，且系数含三角函数。
• 使用幅值平方 $|M|^2$ 的精确导数，计算幅值展开中的系数 $J_0$、$J_2$、$J_4$ 等。
• 对每个 $v_{\rm rel}^2$ 的幂次进行 $\cos\theta$ 的解析积分，避免对复杂相空间积分进行数值积分。
• 通过 $v_{\rm rel}^2$ 的泰勒展开计算相空间因子 $K(s)$，得到 $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$，该差异可通过不同的热平均处理方式解释。
• 当与符号计算软件结合时，该方法将计算时间缩短至数秒，显著提升了遗迹密度研究中的可靠性与效率。
• 即使在存在窄共振的情况下，该方法也具有鲁棒性，但当低速展开失效时需谨慎处理。