中高Z元素原子、离子的电子碰撞电离与激发截面快速计算方法  

Fast computation approach of electron-impact ionization and excitation cross-sections for atoms and ions with medium-and high-Z elements

作　　者：周旭 王川[1,2,3] 胡荣豪 陶治豪[1,2,3] 邓小良 梁亦寒 李晓亚[4] 吕蒙 祝文军[4] Zhou Xu;Wang Chuan;Hu Rong-Hao;Tao Zhi-Hao;Deng Xiao-Liang;Liang Yi-Han;Li Xiao-Ya;Lu Meng;Zhu Wen-Jun(College of Physics,Sichuan University,Chengdu 610064,China;Key Laboratory of High Energy Density Physics and Technology,Ministry of Education,Chengdu 610065,China;Key Laboratory of Radiation Physics and Technology,Ministry of Education,Chengdu 610065,China;Institute of Fluid Physics,China Academy of Engineering Physics,Mianyang 621000,China)

机构地区：[1]四川大学物理学院,成都610064 [2]高能量密度物理及技术教育部重点实验室,成都610065 [3]辐射物理及技术教育部重点实验室,成都610065 [4]中国工程物理研究院流体物理研究所,绵阳621000

出　　处：《物理学报》2024年第10期105-115,共11页Acta Physica Sinica

摘　　要：中高Z元素的原子数据如碰撞电离和碰撞激发截面在聚变工程、X射线与物质相互作用等工程及研究领域有非常广泛的需求.高能量密度等离子体中存在从基态到激发态的原子和各价态离子,其碰撞电离和碰撞激发截面需要分别计算.本文以73号元素钽(Ta)为例,基于相对论性Dirac-Fock理论和扭曲波模型计算了基态Ta原子到Ta^(72+)离子在入射电子能量范围为1-150 keV的碰撞电离与碰撞激发截面,并与相关实验和理论模型符合较好.通过分析Ta的碰撞电离和碰撞激发截面数据中的规律,给出两种减少计算量的方法:对初态能级进行随机抽样、筛去贡献小的反应道,并对计算量的优化程度和误差做出评估.最终结果可在误差5%内将计算效率提高一个数量级,本文所提方法可推广至其他中高Z元素的计算.The atomic data of medium-and high-Z elements,such as electron-impact ionization and excitation cross-sections,possess extensive applications in fields such as fusion science and X-ray interactions with matter.There are atoms and ions in high energy density plasma,with different charge states and energy states ranging from ground states to highly excited states,and the cross-sections of each charge state and energy state need to be calculated.The bottlenecks limiting computational performance are the inevitable relativistic effects of medium-and high-Z elements and the extremely complex electronic configurations.Taking tantalum(Ta)for example,by using the relativistic Dirac-Fock theory and distorted wave model,we compute the electron-impact ionization and excitation cross-sections of Ta from the ground state atom up to Ta^(72+) with the incident electron energy range of 1-150 keV.The detailed configuration accounting(DCA)reaction channel cross-sections are derived by summing and weighting the original detailed level accounting(DLA)cross-sections.After examining the data,two regularities are found.In terms of DLA,the pre-averaging DCA cross-sections have varying initial DLA energy levels but are typically close to each other.There is not a straightforward function that can explain the discrepancies between them.In terms of DCA,inner subshells typically contribute very little to the total cross-section as their ionization and excitation cross-sections are orders of magnitude smaller than those of outer subshells.We provide two techniques to reduce the computational costs based on the regularities.To minimize the total number of DLA reaction channels used in the computation,the initial DLA energy levels can be randomly sampled.Through a Monte Carlo numerical experiment,we determine the appropriate number of sampling points that can reduce the total number of DLA channels by an order of magnitude while maintaining a 5%error margin.In terms of impact ionization,since small cross-section DCA channels are insignificant,

关 键 词：碰撞电离 碰撞激发 扭曲波模型 原子数据 

分 类 号：O562.5[理学—原子与分子物理]