氦原子基态能级的探讨毕业论文正文(16页).doc

-氦原子基态能级的探讨毕业论文正文学号：毕业论文题 目：氦原子基态能级的探讨作 者届 别学 院专 业指导老师职 称完成时间-第 11 页毕业论文（设计）原创性声明本人所呈交的毕业论文（设计）是我在导师的指导下进行的研究工作及取得的研究成果。据我所知，除文中已经注明引用的内容外，本论文（设计）不包含其他个人已经发表或撰写过的研究成果。对本论文（设计）的研究做出重要贡献的个人和集体，均已在文中作了明确说明并表示谢意。 作者签名： 日期： 毕业论文（设计）授权使用说明本论文（设计）作者完全了解*学院有关保留、使用毕业论文（设计）的规定，学校有权保留论文（设计）并向相关部门送交论文（设计）的电子版和纸质版。有权将论文（设计）用于非赢利目的的少量复制并允许论文（设计）进入学校图书馆被查阅。学校可以公布论文（设计）的全部或部分内容。保密的论文（设计）在解密后适用本规定。 作者签名： 指导教师签名： 日期： 日期： 注 意 事 项1.设计（论文）的内容包括：1）封面（按教务处制定的标准封面格式制作）2）原创性声明3）中文摘要（300字左右）、关键词4）外文摘要、关键词 5）目次页（附件不统一编入）6）论文主体部分：引言（或绪论）、正文、结论7）参考文献8）致谢9）附录（对论文支持必要时）2.论文字数要求：理工类设计（论文）正文字数不少于1万字（不包括图纸、程序清单等），文科类论文正文字数不少于1.2万字。3.附件包括：任务书、开题报告、外文译文、译文原文（复印件）。4.文字、图表要求：1）文字通顺，语言流畅，书写字迹工整，打印字体及大小符合要求，无错别字，不准请他人代写2）工程设计类题目的图纸，要求部分用尺规绘制，部分用计算机绘制，所有图纸应符合国家技术标准规范。图表整洁，布局合理，文字注释必须使用工程字书写，不准用徒手画3）毕业论文须用A4单面打印，论文50页以上的双面打印4）图表应绘制于无格子的页面上5）软件工程类课题应有程序清单，并提供电子文档5.装订顺序1）设计（论文）2）附件：按照任务书、开题报告、外文译文、译文原文（复印件）次序装订3）其它摘 要本文主要研究氦原子的基态能级，通过应用双参数变分法，选择适当的试探波函数，求出氦原子基态能级的能量，将其计算值与实验值进行比较并分析了误差，然后比较了用微扰法和变分法求出的氦原子基态能量，通过比较，发现双参数变分法的优势是显而易见的。关键词：变分法；氦原子；基态能量AbstractIn this paper, we mainly study the ground-state energy of the helium atom. By using the double parameters variational method, and selecting an appropriate test wave function, we calculate the helium atom ground-state energy and compare it with the experimental values. Then we compare the perturbation method result with the double parameters variational method result, and we find that the superiority of the double parameters variational method is obvious.Keyword: The variational method；Helium atoms；Ground-state energy目 录摘 要IAbstractII目 录III第一章 绪 论11.1引言11.2 选题的依据和意义11.2.1 选题的依据11.2.2 选题的意义11.3 本文的主要研究内容2第二章 变分法介绍32.1变分法原理32.2变分法求体系基态能量的步骤52.2.1 选取一个参量的尝试波函数52.2.2 选取两个参量的尝试波函数5第三章 氦原子基态能量的变分计算63.1尝试波函数的选择63.2氦原子能量平均值的计算7第四章 结果比较174.1氦原子基态能量的计算值与实验值比较174.2一参变分法与二参变分法比较174.3变分法与微扰法比较18参考文献19致 谢20第一章 绪 论1.1引言 对氦原子基态能级的探讨一般选用微扰法及变分法，本文重点讨论变分法对氦原子基态能级的求解。变分法是解决氦原子和类氦离子基态问题的强有力工具，到目前为止，国内外为追求高精度所选取的变分参数个数已由数百增至数千，在忽略核质量的情况下，它们的非相对论基态波函数和能量的不确定度分别达到1010和10，这对于计算高精度的相对论修正和辐射修正具有非常重要的意义。在量子力学教科书中，一般介绍的近似求解法是微扰法和变分法，而变分法中选择的尝试波函数一般是一个参数型的，例如周世勋编量子力学、曾谨言著量子力学教程等介绍的便是用一个参数型的尝试波函数变分法求氦原子体系基态能级。l.2 选题的依据和意义1.2.1 选题的依据在量子力学中，对于具体物理问题的薛定谔方程，可以精确求解的问题是很少的。在经常遇见的许多问题中，由于体系的哈密顿算符比较复杂，往往不能求得精确的解，而只能求近似解。微扰法和变分法都是用来求问题的近似解的方法。用微扰法求氢原子和类氢离子是比较适合的，但是遇到比氢原子稍微复杂一点的氦原子时，微扰法就不及变分法容易和求解精确。用一参变分法即选用含一个参数的尝试波函数，这种波函数形式简单，其物理意义清晰，物理模型简单，适用于教育教学，但精确度比较低。选用含二参数的尝试波函数，这样的模型相对于更多参数的波函数要简单，又比一参变分法求解精确度高很多，这样既有利于理解怎样用变分法求基态能级，可适用于教学,又能求得比一参法更为精确的数值，因而具有重要的物理意义。1.2.2 选题的意义氦原子是比类氢离子这种单粒子体系复杂但是相对于其他粒子要简单的粒子，研究氦原子这种简单的多粒子体系，对于研究更复杂的多粒子体系具有重要的意义。变分法是解决氦原子和类氦原子的强有力工具，只要选择合适的试探波函数，对于提高求解能级近似值有很大的帮助。l.3 本文的主要研究内容本文主要研究氦原子的基态能级，通过应用双参数变分法，选择适当的试探波函数，求出氦原子基态能级的能量，并将计算值与试验值进行比较，再与用微扰法求出的氦原子基态能量结果进行对比，通过对比体现出用变分法求氦原子基态能级的优越性。第二章 变分法介绍2.1变分法原理 已知量子力学中用微扰法求解问题的条件是体系的哈密顿算符可以分为和两部分： 其中的本征值与本征函数是已知的，而很小。如果这些条件不能满足，微扰法就不能应用。因而在遇到不是很小的情况下，就需要寻找另外的求解方法，量子力学中求解问题的又一种简单方便的方法变分法的应用不受上述条件的限制。设体系哈密顿算符的本征值由小到大的顺序排列为：, ,，, (1)与这些本征值对应的本征函数是：, (2)和是基态能量和基态波函数。为简单起见，我们假定的本征值E是分立的，本征函数组成正交归一系。于是有 (3)设是任意一个归一化的波函数，将按展开： (4)在所描写的状态中，体系能量的平均值是 (5)将(4)式代入(5)式，得应用(3)式有 (6)由于是基态能量，所以(n=1,2, )，在上式中用代替，则 (7)最后一步用了的归一化条件。(6)式和(7)式给出： (8)这个不等式说明，用任意波函数算出的平均值总是大于体系基态能量，而只有当恰好是体系的基态波函数时，的平均值才等于能量。上面讨论中曾假定是归一化的，如果不是归一化的，那么(5)式应该写为： (9)(8)式应写为 (10)根据波函数算出的平均值总是不小于，我们可以选取很多的并算出的平均值，这些平均值中最小的一个最接近于。2.2变分法求体系基态能量的步骤 2.2.1 选取一个参量的尝试波函数选取含有一个参量的尝试波函数代入(5)式和 (9)式，算出平均能量，然后由 (11)求出的最小值。所得结果就是的近似值。2.2.2 选取两个参量的尝试波函数选取含有两个参量、的尝试波函数代入(5)式和 (9)式，算出平均能量，然后由 =0，=0 (12)求出的最小值，所得结果就是的近似值。第三章 氦原子基态能量的变分计算根据第二章变分法原理可知，氦原子的基态能量可由下式计算得到： (13)其中，为哈密顿算符, 为归一化的试探波函数。若波函数未归一化，氦原子的基态能量可由下式计算得到： (13.1)氦原子有一个原子核和两个电子，它们都处于运动状态。由于原子核的质量相对于电子非常大,核的运动比电子的运动要慢的多，因此近似把核看成是固定的，氦原子的哈密顿算符可写为：() (14)其中，是电子质量，与分别代表第一、第二个电子到核的距离；为两个电子之间的距离。为方便计算，采用原子单位(其中能量单位为哈特利,即)，氦原子体系的非相对论哈密顿算符可写为： =()+ (14.1)3.1尝试波函数的选择用变分法求近似解的关键在于选择合适的尝试波函数，尝试波函数选得好，可以在很大程度上提高计算结果的精度。对类氢离子,电子的基态波函数为(取原子单位,以下暂不考虑归一化常数),它所表示的状态有这样一个特点:电子的最概然半径为,电子在半径为的球壳附近分布的概率较大。对于类氦离子,如果不考虑两个电子相互排斥作用,则基态波函数为,二电子的最概然半径皆为,但若考虑两电子的相互作用,利用单参数变分法求得基态近似波函数为： (15)(15)式中=.在这种模型下,由于排斥作用,两个电子的最概然半径增长到.我们认为这个模型仍然比较粗糙,二电子之间的关系考虑得不够,现将试探波函数改进为如下形式：假如二电子可以区分即存在某种差别,我们设想由于二电子的径向排斥作用,电子1的最概然半径为且处于形式为的波函数描述的状态,电子2的最概然半径为且处于形式为的波函数描述的状态,则两个电子的基态波函数可写作： (16),为变分参数。显然,上述说法不满足全同粒子波函数交换对称性的要求,事实上我们不能说清某个电子处于哪个状态.由于基态为态,波函数的自旋部分是反对称的,空间部分要求是对称的,因此,满足交换对称性要求的波函数为: (17)值得强调,如果,则式(17)退化为单参数的.它们究竟会不会相等不是由想象而是由变分原理能量平均值取极小值条件来决定.在计算时，需注意以上波函数是未归一化的。3.2氦原子能量平均值的计算利用(14.1)式和(16)式，我们先计算，然后计算。有： (18)= (18.1) (18.2)+ (18.3)首先，求解(18)式中的(18.1)式，有= (18.1-1) (18.1-2)先求解(18.1)式中的(18.1-1)式，有= (18.1-1.1)接着求解(18.1)式中的(18.1-2)式，有= (18.1-2.1)然后，求解(18)式中的(18.2)式，有= (18.2-1) (18.2-2)先求解(18.2)式中的(18.2-1)式，有= (18.2-1.1)接着求解(18.2)式中的(18.2-2)式，有= (18.2-2.1)最后，求解(18)式中的(18.3)式，这里根据参考文献，有= (18.3-1)计算完了，接着计算。有：= (19)由于波函数（，）是未归一化的波函数，因此由(9)式、(18)式各个计算结果、(19)式计算结果可知： (20)将(20)式右边分子分母同乘可得： - (21)根据(12)式=0，=0，以及(21)式的结果，有： (22) (23)联立(22)、(23)可解得有用的，有：， (24)将(24)式的结果代入(20)式可解得： (25)上面所求采用的是原子单位制，因此在国际单位制中，最后求得的，亦即氦原子的基态能量为： (26) 第四章 结果比较4.1氦原子基态能量的计算值与实验值比较由参考文献可知用实验方法得出氦原子基态能量为，而用变分法算出来的结果由(26)式可知为。其相对误差为： (27)4.2一参变分法与二参变分法比较由参考文献可知用一参变分法计算能量，准确到第一级近似的结果为，而用实验方法得出氦原子基态能量为，可算出其相对误差为： (28)从(27)式和(28)式可以看出，用二参变分法计算氦原子基态能量比用一参变分法计算氦原子基态能量又要精确很多。由此可知，求多粒子体系的基态能量时，变分法的参数的多少对求解结果的影响也是很大的。4.3变分法与微扰法比较由参考文献可知用微扰法计算能量，准确到第一级近似的结果为，而用实验方法得出氦原子基态能量为，可算出其相对误差为： (29)比较（27）式、(28)式和(29)式可以看出，用变分法计算氦原子基态能量比用微扰法计算氦原子基态能量要精确很多。由此可知，求多粒子体系的基态能量时，变分法有着微扰法无法比拟的优越性。参考文献1 周世勋.量子力学M. 北京：高等教育出版社,1979.2 余凤军.氦原子及类氦原子基态的二参数变分法研究J,大学物理,2008,27(5)3 马二俊.类氦原子体系基态能量变分法的数值研究J,大学物理,2003,23(6)4 陈小波.应用双参数法对氦原子基态能级的研究J,太原师范学院学报, 2006,5(4)5 杨汉嵩,李元杰.氦原子组态能级的变分法计算J.黄河科技大学学报.2007,9(4)6陈玉红,李延龙. He原子基态能量与波函数的变分计算J,甘肃工业大学学报,2003,29(4)7 胡先权,许杰,马勇,殷霖.四参数法计算氦原子基态能级J.云南大学学报（自然科学版）,2005,27(5A)8李青仁.用变分法求解氦原子的基态能级J.松辽学刊（自然科学版），1998，(4)9 杨铜锁.氦原子及类氦离子基态能级与波函数研究J.陕西师范大学学报（自然科学版）2002,30(2)10孟现柱.氦原子及类氦原子基态能量的经验公式J.光谱实验室,2004,21(5)11 刘玉孝, 陈玉红,张丽杰. He原子非相对论基态能量的变分计算_英文_J,原子与分子物理学报,2005，22(2)12C.Le Sech,G.Hadinger,M.Auber-Frécon. Variational calculations in a halfspace for the interaction of ground and excited state helium atom with its mirror imageJ.ZEITSCHRIFT FüR PHYSIKD,1994.