计算材料心得体会

《计算材料心得体会》由会员分享，可在线阅读，更多相关《计算材料心得体会（23页珍藏版）》请在金锄头文库上搜索。

1、1湖南工业大学课 程 设 计资 料 袋理 学院（系、部） 2011 2012 学年第 一 学期 课程名称 计算材料学 指导教师 雷军辉 职称 讲师 学生姓名 余晓燕 专业班级 应用物理 081 班 学号 08411200135 题 目 计算BN的弹性常数 成 绩 起止日期 2011 年 12 月 4 日 2011 年 12 月 12 日目 录 清 单序号 材 料 名 称 资料数量 备 注1 课程设计任务书2 课程设计说明书3 课程设计图 张4562湖南工业大学课程设计任务书20112012 学年第 1 学期理学院 学院（系、部） 应用物理学 专业 081 班级课 程 名 称： 计算材料学 一、

2、 设计题目： 计算 BN 的弹性常数二、 完成期限：自 2011 年 12 月 4 日至 2011 年 12 月 12 日共 2 周内容及任务1. DFT 基本理论，CASTEP 使用方法2. 晶体模型的建立与几何优化，相关性质的计算。3. 计算 BN 的弹性常数4. 结果分析5. 报告写作与修改起止日期 工作内容11-12-4-6 熟悉DTT理论，软件安装，认识界面，熟悉基本操作11-12-7 晶体模型建立，进行结构优化，计算物理性质11-12-8 物理性质，力学性质的计算11-12-9 计算 BN 的弹性常数进度安排11-12-10-12 写出课程设计的总结实验报告.，修改成文主要参考资料

3、1 Kohn W, Sham L J, Self-consistent equations including exchange and correlation effects J. Physical review, 1965, 140(4):A1133-A1338. 2 Hohenberg P, Kohn W. Inhomogeneous electron gas J. Physical review, 1964,136(3):B864- B871. 3 谢希德, 陆栋.固体能带理论M.上海：复旦大学出版，1998. 4 Perdew J P, Chevary J A, Vosko S H.

4、 Atoms, molecules, solids, and surfaces: Applications of the generalized gradient approximation for exchange and correlationJ. Physical review B, 1992, 46(11): 6671-6687. 指 导 教 师 （签字）： 年 月 日系 （ 教 研 室 ） 主 任 （ 签字）： 年 月 日3（计算材料）设 计 说 明 书计算 BN 的弹性常数 起止日期： 2011 年 12 月 4 日 至 2011 年 12 月 12 日学 生 姓 名 余 晓 燕班

5、 级 081学 号 08411200135成 绩指 导 教 师 (签 字 )理学院（部）2011 年 12 月 12 日4计算BN的弹性常数背景 ：近年来，随着材料、物理、计算机和数学等学科的发展，应用计算的方法研究材料的结构、能量和性能已成为一门迅速发展的新兴学科-计算材料学。这种方法不仅能进行材料的计算模拟，而且能进行材料的计算机设计和相关性能的预测。随着计算机技术的飞速发展，第一性原理计算的方法在材料的结构和性能等方面的研究已取得了巨大的成功，第一性原理的方法是基于量子力学理论，从电子运动的层次研究材料的结构和相关性能。目前，CASTEP软件的主要功能是对半导体、非线性光学材料、金属氧化

6、物、玻璃、陶瓷等固体材料，对电子工业、航空航天以及石化、化工等工业领域有着非常重要的战略意义。对这些材料而言，其电子的结构与性质，以及表面和界面的性质与行为都非常重要。CASTEP的量子力学方法，为深入了解固体材料的这些性质并进而设计新的材料，提供了强有力的工具。基于密度泛函平面波赝势方法的CASTEP软件可以对许多体系包括像半导体、陶瓷、金属、矿石、沸石等进行第一性原理量子力学计算。典型的功能包括研究表面化学、能带结构、态密度、热学性质和光学性质。它也能够研究体系电荷密度的空间分布和体系波函数。CASTEP还可以用来计算晶体的弹性模量和相关的机械性能，如泊松系数等。半导体和其他固体材料的许多

7、性能由电子性质决定，而电子性质又由原子结构决定，特别是缺陷在改变电子结构上的作用对半导体性质尤为重要。分子模拟，特别是量子物理技术，可用来预测原子和电子结构及分析缺陷对材料性能的影响。CASTEP能有效的研究存在点缺陷、空位、替代杂质、位错等的半导体和其它材料中的的性能。除此以外，它还可以被用来计算固体的振动性质，如声子色散关系、声子态密度等。这些计算结果可以用来分析表面吸附的振动性质，可以解释实验中的振动谱，可以研究在高温高压下的相稳定性等等。总的来说，它可以实现如下的功能：1.计算体系的总能；2.进行结构优化；3.执行动力学任务：在设置的温度和关联参数下，研究体系中原子的运动行为；4.计算

8、周期体系的弹性常数；5.化学反应的过度态搜索。除此之外，计算一些晶体的性质，如能带结构、态密度、声子色散关系、声子态密度、光学性质、应力等。下面介绍一下密度泛函理论、交换关联泛函近似、赝势方法和K-S方程迭代解法。一、基础理论：1. Hohenberg-Kohn 定理和密度泛函理论 ：5密度泛函理论(DFT)是用量子力学的理论求解多电子体系基态能量方法，其核心是用电子密度函数取代波函数作为研究的基本量，由Hohenberg 和Kohn 在1964 年创建1,2。根据量子力量的相关知识，大量电子和原子核相互作用的多粒子体系，在非相对论前提下，系统粒子运动的波函数可以由以下定态薛定谔方程来描述：(

9、1-1)哈密顿量 仅考虑电子-电子作用、电子-原子核作用、原子核-原子核作用以及各个粒子的动能，对其它外场的情况可忽略。因此其哈密顿量可以写成如下形式： (1-2)其中，(1-3)(1-4)(1-5)对于上述方程，是无法直接求解的，必须对多粒子系统的电子能级计算采用一些简化和近似。在实际的多粒子体系中，原子核的质量远远大约电子， 但是运动速度比电子小的多。因此考虑粒子运动时，将原子核的运动和电子的运动分开，考虑核的运动时忽略其电子分布，考虑电子运动时假定原子核处于相对静止的状态，这就是绝热近似3。通过近似，可以独立的处理原子核运动和电子的运动，因此可以将薛定谔方程写成电子运动方程和原子核运动方

10、程。其电子运动方程是：(1-6)原子核的运动方程：(1-7)通过绝热近似，得到了多电子的薛定谔方程，但不能实际求解，要求解上述方程， 必须将多电子问题简化为单电子问题。单电子近似理论的源于H.Thomas 和E.Fermi 在1927 年的工作，就是用粒子数密度表示多粒子的基态系统的能量。P.Hohenberg 和W.Kohn 根据的均匀电子气的理论提出著名的Hohenberg-Kohn 定理1，这个定理包含如下内容：不计自旋的情况下，将粒子数密度函数 表示成全同费米子系统的基态能量的唯一泛函；在粒子数不变的情况下，能量泛函 对正确的粒子数密度 取等于基态能量的极小值。因此，对于基态非间并多粒

11、子系统，不考虑自旋的条件下，其哈密顿算符为(1-8)式(1-8)中，外场作用看成原子核-电子作用，相同的局域势 对外场的作用用 表示。6对于给定的外场 ，多电子系统的能量表示成电子数密度 的泛函为：(1-9)(1-10)(1-11)(1-12)式中， 包括体系中电子之间的相互作用能和电子的动能， 是外场对电子的作用能， 是系统中原子核间的排斥能。在式(1-10) 中，前两项表示无相互作用粒子模型的动能和库仑排斥能，复杂的电子相互作用用交换关联能表示。根据Hohenberg-Kohn 定理，假设能得到能量泛函E( )，然后就能将电子数密度变分，就能确定系统的基态和基态所有的性质，因此确定E( )

12、成为问题的关键所在，而要确定能量泛函E()，必须要确定动能泛函T 、电子数密度 以及交换关联泛函 。为了解决上述问题，W.Kohn 和L.J.Sham 提出了如下假设：假定已知无相互作用的电子系统和未知的有相互作用的电子系统密度函数相同，未知的相互作用电子系统的动能泛函T 可用已知的无相互作用电子系统的动能泛函 来代替；假定密度函数 用N 个单电子波函数 构成，于是有：(1-13)则(1-14)对能量泛函 进行变分得到(1-15)7(1-16)式(1-13)、(1-15)和(1-16) 就是 Kohn-Sham 方程。这个方程的核心就是有相互作用动能泛函能否用未知的无相互作用的动能泛函来代替。

13、而将所有复杂问题都归入中，所以求解Kohn-Sham 方程的关键是找到准确的 ，这样密度泛函理论精确求解量子多体问题的中心是构造交换关联泛函。2交换关联能近似 ：根据密度泛函理论，能将多电子的基态特性问题转化成等效的单电子问题， 而其它所有复杂问题都归结到交换关联能泛函 ，但是交换关联泛函是未知的。因此得到可靠并准确的交换关联 ，成为求解Kohn-Sham 方程的关键。W.kohn 和L.J.Sham 提出了交换关联泛函局域密度近似(LDA，Local Density Approximation)，其基本思想是：在局域密度近似中，利用均匀电子气密度函数来获得非均匀电子气密度泛函 。对变化平缓的

14、密度函数，非均匀交换关联能密度用均匀电子气 代替，则可表示为：(1-17)相应的局域交换关联势可以表示为：(1-18)局域密度近似虽然在大多数的材料计算中显示出巨大的成功，但是由于点r处的交换关联作用仅依赖于点r处的近邻和近邻的电荷密度，因此，对于与均匀电子气或空间变化缓慢的电子气相差太远的体统，LDA 不适用。因此，人们对局域密度近似应用多种方法进行修正，应用较广的是广义梯度近似(GGA)，其泛函与局域密度和密度梯度都有关4，因此能更好的描述真实体系的电子密度的不匀性，其交换关联能密度泛函 可表示为(1-19)目前，在交换关联泛函GGA 的构建上有两个方向，一个是 Becke 为首的， 这类

15、泛函包含若干个实验参数，这些参数通过计算和实验数据来获得，这种形式的好坏由实践的工作所决定的。另一个是Perdew-Wang 91 的，这类泛函以物理规律为基础，不包含实验参数。随着研究的不断深入，不仅出现了非局域的相互作用交换关联泛函，还有密度高阶梯度的近似交换关联泛函，如Vaner Wals 和Meta-GGA 等。GGA 和LDA 相比在能量精确度和开放体系方面更有优势。83. 赝势法：在晶体的近自由电子能带计算中，计算量大而且收敛速度慢。对于固体而言，价电子的化学性质活泼，对于结构和性质的影响较大，而内层电子的能带较窄，较稳定，而且相邻原子的作用对内层电子的状态影响较小。因此，人们关注的是价电子，将原子核和内层电子近似看出粒子实。对于固体中的价电子波函数而言，在离子实的内部区域，变化剧烈，存在若干个节点；而在离子实之间的区域，变化平缓。离子实内部的这一特点要求价电子波函数与内层电子波函数正交，而价电子与内层电子