用于计算半轻子电弱相互作用中谐振子核矩阵元的 Mathematica 脚本

原子核中的半轻子电弱相互作用（如\beta衰变、\mu俘获、带电流和中性流中微子反应以及电子散射）由一组具有确定宇称和角动量的多极算符描述，这些算符是通过对底层核电荷算符和三维流算符进行投影得到的。如果将这些核算符近似为其单体形式，并按核子速度展开至|\vec{p}|/M 阶（其中\vec{p}和 M 分别为核子动量和质量），则可得到一组七个多极算符。核结构计算通常在由谐振子轨道构成的 Slater 行列式基中进行，这种选择允许保持平移不变性。多极算符的谐振子单粒子矩阵元可以解析计算，并表示为 q^2 的有限多项式，其中 q 是三维动量转移的大小。虽然此类矩阵元的结果以表格形式提供，但对量子数有一定限制，确定响应函数解析形式的任务仍然相当繁琐，需要将表格矩阵元与核密度矩阵折叠，并进行后续代数运算以评估算符乘积。在此，我们提供了一个用于生成这些矩阵元的 Mathematica 脚本，使用户能够通过符号操纵执行所有此类计算。这将消除手工计算可能伴随的错误，并加快电弱核截面和反应率的计算。我们通过计算^{12}C 中带电流和中性流中微子散射的截面来说明新脚本的使用。