在三种理论水平下计算的排斥原子间势

Repulsive interatomic potentials calculated at three levels of theory

摘要 Abstract

核间距远小于平衡键长时的高能排斥相互作用是决定核阻止本领和千电子伏特(keV)及兆电子伏特(MeV)辐射事件中原子散射的关键量。这种相互作用传统上在冻结原子电子密度的轨道自由密度泛函理论框架内建模,遵循齐格勒-比尔斯马克-利特马克(ZBL)模型。本文中,我们利用ZBL模型计算了特定原子对的排斥原子间势,并将其与两种量子化学计算方法进行比较——即基于灵活高斯基组的第二阶莫耳-普莱舍特微扰理论以及采用数值原子轨道基组的密度泛函理论——这些方法超越了ZBL模型的局限性,允许在计算中放松密度。我们表明,这两种量子化学模型预测的排斥原子间势对于超过30电子伏特的势能值一致在约±1%以内,而ZBL特定原子对势和通用ZBL势则与这些计算结果存在较大差异。我们为所有原子对$Z_1$-$Z_2$($1 \leq Z_i \leq 92$)提供了新的基于三个指数的屏蔽函数拟合,与计算结果一致在约±2%以内。我们利用新的势能模拟了单晶硅中的离子注入深度分布,并显示出与实验的良好一致性。然而,我们也展示了在通道效应条件下,势能的吸引部分可以影响深度分布。所有计算得到的原子间势的完整数据集及其数据分析拟合均作为开放获取资源共享。

The high-energy repulsive interaction between nuclei at distances much smaller than the equilibrium bond length is the key quantity determining the nuclear stopping power and atom scattering in keV and MeV radiation events. This interaction is traditionally modeled within orbital-free density functional theory with frozen atomic electron densities, following the Ziegler-Biersack-Littmark (ZBL) model. In this work, we calculate atom pair specific repulsive interatomic potentials with the ZBL model, and compare them to two kinds of quantum chemical calculations - second-order M{\o}ller-Plesset perturbation theory in flexible Gaussian basis sets as well as density functional theory with numerical atomic orbital basis sets - which go well beyond the limitations in the ZBL model, allowing the density to relax in the calculations. We show that the repulsive interatomic potentials predicted by the two quantum chemical models agree within $\sim$ 1% for potential energies above 30 eV, while the ZBL pair-specific potentials and universal ZBL potentials differ much more from either of these calculations. We provide new pair-specific fits of the screening functions in terms of 3 exponentials to the calculations for all pairs $Z_1$-$Z_2$ for $1 \leq Z_i \leq 92$, and show that they agree within $\sim 2$% with the raw data. We use the new potentials to simulate ion implantation depth profiles in single crystalline Si and show very good agreement with experiment. However, we also show that under channeling conditions, the attractive part of the potential can affect the depth profiles. The full data sets of all the calculated interatomic potentials as well as analytic fits to the data are shared as open access.

在三种理论水平下计算的排斥原子间势 - arXiv