[论文解读] Mutual neutralization in low energy H^+ + H^- collisions

Michael Stenrup, Åsa Larson+1 more|ArXiv.org|Feb 11, 2009

High-Energy Particle Collisions Research被引用 3

一句话总结

本研究采用分子紧密耦合方法，对低能 H⁺ + H⁻ 碰撞中的互中性化过程进行了完全量子力学的从头算研究，全面考虑了核量子效应和核身份（自旋统计）。研究计算了 H(1) + H(n)，n=1,2,3 的态分辨及总中性化截面，结果与先前的理论研究吻合良好，但比 Moseley 等人（1970 年）的实验截面低两到三倍，表明该实验在低能区可能存在高估。

ABSTRACT

The mutual neutralization of H^+ and H^- ions at low collision energies is studied by means of a molecular close-coupling approach. All degrees of freedom are treated at the full quantum level taking into account also the identity of the nuclei. The relevant ^1Sigma_g,u^+ electronic states and their associated nonadiabatic radial couplings are calculated for internuclear distances between 0.5 and 50 a_0. Following a transformation into a strictly diabatic basis, these quantities enter into a set of coupled equations for the motion of the nuclei. Numerical solution of these equations allows the cross sections for scattering into the H(1)+H^*(n), n=1,2,3 channels to be calculated. In the present paper, results are reported for the collision energy region 0.001 to 100 eV, with special emphasis on the important energy region below 10 eV. The low temperature rate coefficient is obtained from a parametrization of the calculated neutralization cross section and is estimated to be valid in the range 10 to 10 000 K.

研究动机与目标

• 为最简单的离子-离子反应：H⁺ + H⁻ → H(1) + H(n) 提供基准质量的量子力学研究，以实现精确的理论验证。
• 解决低能实验截面（Moseley 等人，1970 年）与先前理论预测之间的差异，特别是在 0.001–10 eV 范围内，该范围与星际环境密切相关。
• 在散射计算中严格考虑核量子效应、电子非绝热耦合以及核身份（自旋统计）。
• 计算态特异及总中性化截面，并推导出适用于 10 至 10,000 K 温度范围的低温速率系数，以用于天体化学建模。

提出的方法

• 采用分子紧密耦合方法，对所有核和电子自由度（包括核自旋和身份）进行完全量子处理。
• 利用从头算电子结构方法，计算 1Σg⁺ 和 1Σu⁺ 电子态在 0.5 至 50 a₀ 间距范围内的绝热势能曲线及非绝热径向耦合矩阵元。
• 将绝热表象转换为严格绝热基，以实现对耦合核薛定谔方程的精确求解。
• 通过数值求解所得耦合方程，获得 H(1) + H(n)，n=1,2,3 中性化过程的散射振幅和截面。
• 利用适当的核自旋统计，结合 gerade 和 ungerade 对称性贡献，以考虑全同粒子交换效应。
• 通过拟合计算得到的截面，推导出低温速率系数，其有效范围为 10–10,000 K。

实验结果

研究问题

• RQ1 在低能区（0.001–10 eV）H⁺ + H⁻ 碰撞中，H(n) 最终态（n=1,2,3）的准确态分辨中性化截面是多少？
• RQ2 非绝热耦合与核量子效应如何影响这一典型离子-离子反应中的互中性化动力学？
• RQ3 为何 Moseley 等人（1970 年）的低能实验截面与理论预测存在差异？哪一方更可靠？
• RQ4 在低碰撞能量下，包含核身份（自旋统计）对总截面的影响有多大，特别是在 n=2 生成过程中？
• RQ5 该反应的低温速率系数是多少？与实验和理论估计相比如何？

主要发现

• 在 0.001–10 eV 范围内，计算得到的总中性化截面与 Bates 和 Lewis（1955 年）、Fussen 和 Kubach（1986 年）以及 Eerden 等人（1995 年）的先前理论研究结果吻合良好。
• 本研究结果在整体上比 Moseley 等人（1970 年）报告的实验截面低两到三倍，表明该实验测量在低能区可能存在高估。
• 在低于若干 eV 的碰撞能量下，大部分中性产物以 H(n=3) 最终态生成，与常规预期一致。
• 在低能区，gerade 和 ungerade 对称性贡献的加权仅在 n=2 生成时显著，此时截面较小。
• 低温速率系数已通过拟合获得，估计在 10 至 10,000 K 的天体物理相关温度范围内有效。
• 本研究证实，在低能区核身份效应对总截面的影响可忽略不计，仅在 n=2 生成的微小贡献中略有影响。