第三章 密度泛函理论 3.2.doc

文档介绍：晶体总能量和结合能
由密度泛函研究晶体总能量和结合能
晶体总能量及其相对于结构自由度参数的总能量计算一直是个重要的课题。晶体总能量(不包括核的动能部分)可分成两部分:一部分是原子核与芯电子组成的离子实能量,这部分能量基本上与晶体结构()无关,上一个常数,约为Ry/原子;另一部分是总能量与离子实能量之差,即离子与价电子的相互作用、离子间的相互作用以及价电子间的相互作用,约为Ry/原子。赝势方法中常把总能量中不变的常数部分(即离子实能量)设为零。晶体的结合能定义为这部分总能加上核的动能(通常是零点振动能)与孤立原子能之差。下面讨论赝势方法中晶体总能量计算,在密度泛函理论中它是晶格电子能量与离子实排斥能之和:
()
其中电子与外场相互作用能量为
()
电子间库仑相互作用能为
()
而交换关联能在LDA框架下为
()
离子间库仑相互作用能为
()
这里表示价电子数,表示晶体格矢,表示原胞内原子的相对位矢。动能泛函可通过 Kohn-Sham方程() 用单电子能量表示为:
()
于是可得到
()
在晶体总能量计算中,如果利用晶格的周期性,将上面实空间(r空间)总能量表示转换到动量空间(k空间),将使计算变得十分便利,
()
这里,,和分别是电子间库仑相互作用势,电子交换关联能,交换关联势和电子数密度的傅里叶分量,为原胞体积。其中电子间库仑相互作用势的傅里叶分量可以从泊松方程
()
用傅里叶展开而得到
()
而对于交换关联能和交换关联势,一般根据所采用的交换关联能近似在实空间求得和后,用傅里叶变换至动量空间,得到和。电子数密度的傅里叶分量既可以在实空间求得后变换至傅里叶空间,也可以直接在动量空间求得。
在实际计算中还需要作些数学处理,因为离子间库仑相互作用能的求和() 包括无穷项,是发散的,而在不考虑其他外场时,一般只考虑离子与电子的库仑相互作用,
()
它的傅里叶分量在也是发散的。该项在总能量计算式() 中不出现,而在Kohn-Sham方程(.