2002固体物理讲义g0302_第1页
2002固体物理讲义g0302_第2页
2002固体物理讲义g0302_第3页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

3.2 晶体总能量和结合能由密度泛函研究晶体总能量和结合能晶体总能量及其相对于结构自由度参数的总能量计算一直是个重要的课题。晶体总能量（不包括核的动能部分）可分成两部分：一部分是原子核与芯电子组成的离子实能量，这部分能量基本上与晶体结构（）无关，上一个常数，约为Ry/原子；另一部分是总能量与离子实能量之差，即离子与价电子的相互作用、离子间的相互作用以及价电子间的相互作用，约为Ry/原子。赝势方法中常把总能量中不变的常数部分（即离子实能量）设为零。晶体的结合能定义为这部分总能加上核的动能（通常是零点振动能）与孤立原子能之差。下面讨论赝势方法中晶体总能量计算，在密度泛函理论中它是晶格电子能量与离子实排斥能之和： (3.2.1)其中电子与外场相互作用能量为 (3.2.2)电子间库仑相互作用能为 (3.2.3)而交换关联能在LDA框架下为 (3.2.4)离子间库仑相互作用能为 (3.2.5)这里表示价电子数，表示晶体格矢，表示原胞内原子的相对位矢。动能泛函可通过 Kohn-Sham方程 (3.2.7) 用单电子能量表示为： (3.2.6)于是可得到 (3.2.7)在晶体总能量计算中，如果利用晶格的周期性，将上面实空间（r空间）总能量表示转换到动量空间（k空间），将使计算变得十分便利， (3.2.8)这里，和分别是电子间库仑相互作用势，电子交换关联能，交换关联势和电子数密度的傅里叶分量，为原胞体积。其中电子间库仑相互作用势的傅里叶分量可以从泊松方程 (3.2.9)用傅里叶展开而得到 (3.2.10)而对于交换关联能和交换关联势，一般根据所采用的交换关联能近似在实空间求得和后，用傅里叶变换至动量空间，得到和。电子数密度的傅里叶分量既可以在实空间求得后变换至傅里叶空间，也可以直接在动量空间求得。在实际计算中还需要作些数学处理，因为离子间库仑相互作用能的求和 (3.2.5) 包括无穷项，是发散的，而在不考虑其他外场时，一般只考虑离子与电子的库仑相互作用， (3.2.11)它的傅里叶分量在也是发散的。该项在总能量计算式 (3.2.8) 中不出现，而在Kohn-Sham方程 (3.1.68) 中出现，既在能带计算时用到。上述三项单独都是发散的。但因为整个系统是电中性的，这些发散项相互抵消，为一常数。因此，在解Kohn-Sham方程 (3.1.68) 时，可先将和同时置零，这相当于势能作一平移，或说重新定义真空能级，而在总能量计算中补偿这一平移。记这些发散项的和为 (3.2.12)下面来确定和。对于形如的外场，将它的傅里叶分量在附近展开： (3.2.13)同样展开，有 (3.2.14)将 (3.2.13) 和 (3.2.14) 两式代入式 (3.2.12)，去掉的高次项后，可得 (3.2.15) (3.2.16)这里是误差函数，参数原则上是任意的，如果取得足够多，上述求和是与无关的。一般选取使在正格子和倒格子空间收
人人文库> 全部分类> 教育资料 > 课件下载

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论