用变分法计算基态氦原子能量.DOC

用变分法计算基态氦原子能量 第 PAGE 1页 共 NUMPAGES 10页 用变分法计算基态氦原子能量 李勇 (安庆师范学院物理与电气工程学院 安徽 安庆 246011) 指导教师：余春日 摘要：介绍了物理学中变分法的基本思想及其应用思路。分别运用微扰法和母函数及S 函数的变分法、带松弛坐标参数的变分法、坐标尺度伸缩的变分法和四参数基态波函数的变分法计算了基态氦原子能量，总结和比较了不同计算方法的优缺点。研究表明：根据不同的具体情况运用恰当的变分计算可使计算结果接近实验值。 关键词：变分法，氦原子，基态能量， 波函数 引言：关于变分法在物理学中的应用，人们早已开始研究和重视。近几十年来，变分法出现了许 多新的形式。例如，应用带松弛坐标参数法［１］，四参数法［２］等。与此同时，随着计算机技术 的不断进步，人们对计算基态氦原子能量的精度提出了更高的要求，文献［３］介绍了用变分法计算 氦原子基态能量但直接给出了结果，并没有计算过程和具体的计算方法。因此，本文探讨如何准确计 算基态氦原子能量，提高计算的精确度，对原子波函数和能量的精确计算在分子系统结合能的计算、 原子的超精细结的计算和对定态之间的跃迁几率的计算三个方面有具有重要的意义。 1. 物理学中的变分法 1.1变分法思想 设体系哈密顿算符的本征值由小到大的顺序列为与这些本征值对应的 本征函数是，和是基态能量和基态波函数。为简单起见，我们假定的本 征值是分立的，本征函数组成正交归一系。于是有 设是任意一个归一化的波函数，将按展开： 在上面所描写的状态中，体系能量的平均值是 这个不等式说明，用任意波函数算出的平均值总是一般大于体系基态能量，而且只有当 是体系的基态波函数时，的平均值才等于基态能量。 1.2变分法的应用思路 根据波函数算出的的平均值一般不小于，而我们可以选取足够多并算出的平均值， 这些平均值中最小的一个值最接近于基态能量。用变分法求体系基态能量的关键是：选取含有参量 的尝试波函数，并算出平均能量。然后由，求出的最小值。所得结果就是 的近似值。实际的应用中，常设计出一个含参数的试探波函数，算出平均值，然后 用求极值的方法求出函数的极小值，把它作为基态能量的近似值。这个参数可以是一个也可以 是一组。因此变分法处理问题结果的好坏，主要看试探函数的选择。对于函数选择的试探函数要满足 以下的几个条件： （1）关于束缚态态函数的一些基本要求，包括态函数的连续性和平方可积要求，而平方可积意味 着态函数在无限边界处必定是收敛的。 （2）对于基态，态函数是没有节点的，而第几个激发态就有几个节点。 （3）所取的试探态函数要使粒子被束缚住的主要区域的概率尽可能的大。 2. 变分法与微扰法的比较 2.1 微扰法计算基态氦原子能量 为了以后便于把用变分法计算的结果与微扰论算得的结果加以比较, 现先介绍用微扰论计基态氦 原子能量的方法和结果。若不考虑原子核运动所引起的影响, 氦原子的哈密顿算符可以写为: （１） 其中,和分别代表第一、第二电子的拉普拉斯算符，Z为原子序数。为简化起见，本文采用原子 单位。原子单位是以氢原子的第一波尔轨道半径为其长度单位，以氢原子离化能 作为能量单位。使用这个单位时，氦原子的哈密顿算符可以写为： （２） 现在用微扰论来解这个问题。令 其中为未受扰动时的哈密顿算符, 为微扰哈密顿算符。这个问题的零级近似波函数是下列方程的解。 若令 则方程（１）可分离为下列2 个方程 这两个方程正好与类氢原子的波方程相同。当原子处于最低态时可知: （３） 现在来计算一级微扰能量 （４） 上式右边第二积分可视为电荷在点1上的电势,并且此电荷以密度球对称地分布在核的周围。在 带电球外一点的位于球的全部电荷集中在球心时是一样的；而在球内一点的电势与该点的位置无关,都 等于球心的电势。应用此定理计算第二积分，可得: 故准确到一级近似, 氦原子基态能量为: （５） 2.2 变分法计算基态氦原子能量 若把式（３）