FCPC方法在计算类锂原子体系时的应用

2014-06-28张晓娟郑维民

物理通报 2014年10期

张晓娟郑维民

(甘肃民族师范学院甘肃甘南 747000)

1 引言

近年来，高离化原子体系结构的实验和理论研究成为发展最为迅猛的领域之一.高离化原子一般指失去多余一个或两个电子的原子.高离化态原子光谱数据在X-ray激光、天体等离子体物理、磁约束聚变(MCF)和惯性约束聚变(ICF)研究中有特别重要的意义和广泛的应用价值.而高离化原子在地球上相当稀少，地球表面大气层对X射线有较高的吸收率，使得对高离化原子的研究和应用受到了很大限制.近年来由于诸多高新科技领域需求的推动和实验技术的发展，出现了以光电离和电子碰撞电离为主要机制的产生高离化原子的实验装置，使高核电荷数的高离化原子的能级结构与光谱特性越来越受到有关研究人员的重视[1～5].

2 FCPC方法简介

具有1s2-原子实的类锂原子体系能级结构的高精度计算是目前原子结构理论富有挑战性的课题之一[5].1991年Chung建立了全实加关联(Full-core plus correlation，简称FCPC)理论方法[6].FCPC方法的中心思想是,采用平方可积的slater型基函数,预先确定足够好的1s2-原子实波函数，将其作为三电子体系波函数中的单独一项，价电子效应通过将其乘以单电子slater轨道的线性组合予以考虑;原子实的弛豫和壳层间的电子关联则通过附加另一较大的CI波函数加以描述.由于所选取的体系波函数的特殊性，该方法既充分发挥了CI方法在描述电子关联效应方面的潜力，同时又较好地克服了传统的CI方法所受到的数值收敛问题的困扰.

3 FCPC方法的理论计算

下面主要介绍一下FCPC方法在计算类锂原子的激发能和电离能及能级的精细结构中的应用.类锂原子的1s2-原子实和1s2nl(l=d,f)态的非相对论能量均用Rayleigh-Ritz变分法(通过求H0的期待值的极值)计算

(1)

式中H0是类锂体系的非相对论哈密顿算符

(2)

上式中ri是第i电子与原子核的距离，rij是第i电子与第j电子之间的距离，对i

为了得到精确的结果，必须考虑相对论和质量极化修正.来自这两个效应的修正ΔE用一级微扰理论计算

ΔE=〈Ψ|H′|Ψ〉

(3)

H′=H1+H2+H3+H4+H5

(4)

其中，各项的具体形式如下

(5)

是电子动能的相对论修正，它是由于电子质量随运动速度的改变而对能量值的修正;

(6)

是Darwin项，可以认为该项来源于电子的相对论感生电矩，或者来源于电子的相对论非定域性;

(7)

是电子-电子接触(electorn-electrnoocntact)项，来源于自旋磁矩和自旋磁矩之间的相互作用;

(8)

是轨道-轨道相互作用，即电子轨道磁矩之间的相互作用.

以上4项即为相对论效应引起的能量修正.式(5)～(8)中的c即

(9)

为了得到更精确的结果需考虑QED效应(质量极化效应)，因为随着核电荷数的增大该项的贡献也越来越大，因此在计算组态能量时必须考虑QED效应.由于体系的电离能等于体系总能量与原子实能量之差，所以我们假定，在计算体系电离能时，1s2-原子实的QED效应基本上抵消,仅需考虑价电子的QED修正即可.对于类锂原子的1s2nl态来说，该项修正为