任意温度密度下的原子结构求解

本文采用非自洽TF统计模型,对具有较高实用价值的含温有界原子结构进行计算了。文中根据分子轨道理论、重点讨论了有界原子径向波函数的两类外边界条件的选取及其物理含义。本文结合类氢近似下的径向薛定谔方程和波函数,给出了径向波函数的中心边界条件。文中采用四级四阶龙格库塔差分格式,取代了沿用已久的罗曼诺夫差分格式,彻底更换了以势的拟合形式V(r)=-z/r+A+Br求解波函数初始值的方法,使波函数在r=0处的一阶导数值的选取不再依赖势的形式。本文在简单的含温有界TF势的基础上,进行了求解。鉴于简单TF势的不合理性与Shalitin电子自作用修正的某些缺陷,本文对电子的自作用修正作了进一步的改进,使平均原子模型的电子自作用修正在概念上更加合理、物理图像更加明确、计算结果得以较大的改善。本文针对Fe^(26),Rb这两种具有代表性的原子,采用改进后的电子自作用修正方法,详细计算了在不同温度不同物质密度下的能级、波函数、电子占据数、化学势等。为了较准确地计算重Z原子,我们还考虑了含电子自作用修正的相对论质量修正、自旋轨道耦合修正以及Darwin项修正,计算出了相对论三项效应对轻、重元素能级的贡献大小。本文把两类电子自作用修正后的计算结果同相对论自洽场HFS方法的计算结果作了比较,证实了?