第六讲第一原理计算方法简介及castep使用_图文
- 格式:ppt
- 大小:2.70 MB
- 文档页数:5
CASTEP 计算理论总结XBAPRSCASTEP 特点是适合于计算周期性结构,对于非周期性结构一般要将特定的部分作为周期性结构,建立单位晶胞后方可进行计算。
CASTEP 计算步骤可以概括为三步:首先建立周期性的目标物质的晶体;其次对建立的结构进行优化,这包括体系电子能量的最小化和几何结构稳定化。
最后是计算要求的性质,如电子密度分布(Electron density distribution),能带结构(Band structure)、状态密度分布(Densityof states)、声子能谱(Phonon spectrum)、声子状态密度分布(DOS of phonon),轨道群分布(Orbitalpopulations)以及光学性质(Optical properties)等。
本文主要将就各个步骤中的计算原理进行阐述,并结合作者对计算实践经验,在文章最后给出了几个计算事例,以备参考。
CASTEP 计算总体上是基于DFT ,但实现运算具体理论有:离子实与价电子之间相互作用采用赝势来表示;超晶胞的周期性边界条件;平面波基组描述体系电子波函数;广泛采用快速fast Fourier transform (FFT) 对体系哈密顿量进行数值化计算;体系电子自恰能量最小化采用迭带计算的方式;采用最普遍使用的交换-相关泛函实现DFT 的计算,泛函含概了精确形式和屏蔽形式。
一, CASTEP 中周期性结构计算优点与MS 中其他计算包不同,非周期性结构在CASTEP 中不能进行计算。
将晶面或非周期性结构置于一个有限长度空间方盒中,按照周期性结构来处理,周期性空间方盒形状没有限制。
之所以采用周期性结构原因在于:依据Bloch 定理,周期性结构中每个电子波函数可以表示为一个波函数与晶体周期部分乘积的形式。
他们可以用以晶体倒易点阵矢量为波矢一系列分离平面波函数来展开。
这样每个电子波函数就是平面波和,但最主要的是可以极大简化Kohn-Sham 方程。
第一性原理计算方法引言前面讲述的有限元和有限差分等数值计算方法中,求解的过程中需要知道一些物理参量,如温度场方程中的热传导系数和浓度场方程中的扩散系数等,这些参量随着材料的不同而改变,需要通过实验或经验来确定,所以这些方法也叫做经验或者半经验方法。
而第一性原理计算方法只需要知道几个基本的物理参量如电子质量、电子的电量、原子的质量、原子的核电荷数、布朗克常数、波尔半径等,而不需要知道那些经验或半经验的参数。
第一性原理计算方法的理论基础是量子力学,即对体系薛定额方程的求解。
量子力学是反映微观粒子运动规律的理论。
量子力学的出现,使得人们对于物质微观结构的认识日益深入。
原则上,量子力学完全可以解释原子之间是如何相互作用从而构成固体的。
量子力学在物理、化学、材料、生物以及许多现代技术中得到了广泛的应用。
以量子力学为基础而发展起来的固体物理学,使人们搞清了“为什么物质有半导体、导体、绝缘体的区别”等一系列基本问题,引发了通讯技术和计算机技术的重大变革。
目前,结合高速发展的计算机技术建立起来的计算材料科学已经在材料设计、物性研究方面发挥着越来越重要的作用。
但是固体是具有~1023数量级粒子的多粒子系统,具体应用量子理论时会导致物理方程过于复杂以至于无法求解,所以将量子理论应用于固体系统必须采用一些近似和简化。
绝热近似(Born-Oppenheimei近似)将电子的运动和原子核的运动分开,从而将多粒子系统简化为多电子系统。
Hartree-Fock近似将多电子问题简化为仅与以单电子波函数(分子轨道)为基本变量的单粒子问题。
但是其中波函数的行列式表示使得求解需要非常大的计算量;对于研究分子体系,他可以作为一个很好的出发点,但是不适于研究固态体系。
1964年,Hohenberg 和Kohn 提出了严格的密度泛函理论(Density Functional Theory, DFT )。
它建立在非均匀电子气理论基础之上,以粒子数密度()r ρ作为基本变量。
CASTEP 实战守则模型选定我们在进行材料物理模拟所需要做的第一步 (也是很重要的一步) 是模型的选定或建构。
CASTEP 虽然内建了很多功能来预测晶胞参数 (边长,夹角) 与原子位置,但仍然仰赖使用者告诉它 "要进行计算的系统是什么"。
在选定模型时,我们需切记如果系统内原子太多或是超晶胞体积太大,则计算量都会以平面波数的 3 次方增加到计算机难以负荷或使用者难以等待的程度,]此,我时时应考虑设计出一个足以表现出我们所想要研究的物理,而却又能使所有采用的超晶胞越少越好的模型。
如果研究上涉及到一系列原子数不同的大小系统需要做计算,最好能先计算小结构,不要一开始就送入大结构到计算机中。
如果要仿真的系统是含有杂质,则单位晶胞必须进一步放大成超晶胞以便使化学成分里的分数变成整数,因此晶胞会变得很大。
在某些特殊的情况,相互取代的元素种类是很类似的,(即在化学行为上类似),则下一个版本的CASTEP会提供一种叫做虚拟晶体近似 (Virtural Crystal Approximation,VCA) 的方法,则模型里面的原子就可以指定成如0.3A元素加0.7B元素这种样子,因此总可以以最小晶胞来做计算的模型。
但这种方法的精确度通常只适用于合金材料,故要小心使用。
模型的选定有许多人为抉择会含在里面,例如表面计算的层数,因此有些情况也需要进行所谓的收敛性测试。
Vps (poseudopotential)选择Vps选单:MS接口的设定是选用USP优先于NCP,USP有加速计算与减少内存使用的效果,其精确度也与NCP (norm-conserving potential) 相当。
至于什么时候使用NPC 呢 (就是在poseudopotential 选项中那些延伸文件名是 .recpot者),使用到NCP的场合有:1.某些CASTEP计算的功能尚未支持到USP,因此需要选用NCP。
2.为了要与已经发展的文献比较或进行验证3.对计算的结果存疑者,能提供『多一种选择』( 注:至于延伸档名是*.psp者也是属于norm-conserving的一种,是TM potential,在文献上也常被使用,但所需的截止动能较高,因此计算代价较大。
castep动力学计算标题:探索材料科学领域的新篇章:CASTEP动力学计算导语:材料科学领域一直在追求新材料的发现和性能优化,而CASTEP动力学计算作为一种重要的计算方法,为材料科学家们提供了强大的工具。
本文将从人类的视角出发,探索CASTEP动力学计算的应用和意义,希望为读者带来全新的科学发现之旅。
1. 引言材料科学是一门关乎人类生活的重要学科,通过对材料微观结构与性能的研究,可以提升生活质量和推动社会发展。
而CASTEP动力学计算作为一种先进的计算方法,能够模拟和预测材料的结构和性质变化,为材料科学家们提供了强有力的支持。
2. CASTEP动力学计算的原理CASTEP动力学计算基于第一性原理,通过求解薛定谔方程来描述材料的电子结构和原子核的运动。
通过计算材料的势能面、原子间的相互作用和能量变化等信息,可以预测材料的结构演化、相变和响应性质等重要参数。
3. CASTEP动力学计算的应用3.1 新材料的发现CASTEP动力学计算能够通过模拟材料的结构和性质变化,帮助科学家们发现新材料。
通过对材料的能带结构、电子密度等进行计算,可以预测材料的电子性质,从而为新材料的设计和合成提供指导。
3.2 材料性能的优化除了新材料的发现,CASTEP动力学计算还能够帮助科学家们优化已有材料的性能。
通过调整材料的晶格结构、掺杂或添加其他元素,可以改变材料的光学、电子、磁性等性质,从而提高材料的应用性能。
4. CASTEP动力学计算的意义4.1 省时高效CASTEP动力学计算能够在计算机上模拟材料的结构和性质变化,避免了传统实验方法的耗时和耗费资源的缺点。
科学家们可以通过计算快速地获取材料的相关信息,从而指导实验工作的开展。
4.2 推动科学发展CASTEP动力学计算的应用不仅可以帮助科学家们解决材料科学中的问题,还能够推动科学的发展。
通过计算模拟和预测,科学家们可以更好地理解材料的性质和行为规律,为材料科学的理论研究提供新的思路和方法。
4CASTEP的使用方法及应用CASTEP是一个用于计算固体和分子材料的电子结构和晶体结构的第一性原理程序。
它采用了密度泛函理论(DFT)和平面波基组进行计算,可以模拟材料的电子结构、能带结构、振动态和晶体结构等属性。
CASTEP 广泛应用于材料科学、化学、物理学以及生物科学等领域。
下面将介绍CASTEP的使用方法以及其在材料科学研究中的应用。
使用方法:1.安装和准备输入文件:首先需要安装CASTEP程序,并准备好模拟的材料的晶体结构文件。
CASTEP采用标准的输入文件格式,用于描述模拟系统的原子坐标、晶格参数以及计算参数等。
2.执行计算:使用CASTEP提供的命令行界面或者图形化界面,加载输入文件,并选择所需的计算任务,如计算材料的电子结构、优化晶体结构、计算能带结构等。
在计算过程中,CASTEP会自动构建系统的电子密度和施加合适的周期性边界条件。
3.分析和理解计算结果:CASTEP会输出计算结果,如能带图、电子密度分布、振动频谱等。
用户可以通过分析这些结果来理解材料的性质,如能带结构揭示了材料的导电性质,电荷密度分布有助于理解材料的化学键等。
应用:1.材料的能带结构计算:通过CASTEP可以计算材料的能带结构,揭示材料的电子能级和导电性质。
这对于理解材料的禁带宽度、载流子的输运行为以及材料的光学性质等都非常重要。
能带结构的计算成果可以用于材料设计和功能调控。
2.晶体结构的优化:CASTEP可以通过结构优化算法,寻找材料的最稳定晶格参数和原子坐标。
通过晶体结构的优化,可以预测材料的物理和化学性质,指导实验合成设计以及改进材料性能。
晶体结构的优化可以用于研究材料的力学性质、相变过程等。
3.材料的电子密度分布计算:CASTEP可以计算材料的电子密度分布图,展示了材料中电子云的分布情况。
通过分析电子密度分布,可以了解材料之间的化学键类型和强度,理解分子的极性和分子间相互作用力。
电子密度分布的计算可以用于理解材料的反应性和化学稳定性。
CASTEP软件的主要功能及原理半导体、非线性光学材料、金属氧化物、玻璃、陶瓷等固体材料,对电子工业、航空航天以及石化、化工等工业领域有着非常重要的战略意义。
对这些材料而言,其电子的结构与性质,以及表面和界面的性质与行为都非常重要。
半导体和其他固体材料的许多性能由电子性质决定,而电子性质又由原子结构决定,特别是缺陷在改变电子结构上的作用对半导体性质尤为重要。
分子模拟,特别是量子物理技术,可用来预测原子和电子结构及分析缺陷对材料性能的影响。
CASTEP能有效的研究存在点缺陷、空位、替代杂质、位错等的半导体和其它材料中的的性能。
CASTEP的量子力学方法,为深入了解固体材料的这些性质并进而设计新的材料,提供了强有力的工具。
1.CASTEP软件的主要功能基于密度泛函平面波赝势方法的CASTEP软件可以对许多体系包括象半导体、陶瓷、金属、矿石、沸石等进行第一原理量子力学计算。
典型的功能包括研究表面化学、带结构、态密度、和光学性质。
它也能够研究体系电荷密度的空间分布和体系波函数。
CASTEP还可以用来计算晶体的弹性模量和相关的机械性能,如泊松系数等。
CASTEP中的过度态搜索工具提供了研究气相或者材料表面化学反应的技术。
总的来说,它可以实现:计算体系的总能;进行结构优化;执行动力学任务;在设置的温度和关联参数下,研究体系中原子的运动行为;计算周期体系的弹性常数;化学反应的过度态搜索等。
除此之外,计算一些晶体的性质,如能带结构、态密度、聚居数分析、声子色散关系、声子太密度、光学性质、应力等。
量子力学计算精确度高但计算密集。
直到最近,表征固体和表面所需的扩展体系的量子力学模拟对大多数研究者来说才切实可行。
然而,不断发展的计算机功能和算法的进步使这种计算越来越容易实现。
与许多该领域一流专家一起工作推动固体量子力学发展,通过提供可方便直接进入上述CASTEP计算方法中。
2.CASTEP软件的主要理论(1)密度泛函理论(DFT)CASTEP的理论基础是电荷密度泛函理论在局域电荷密度近似(LDA)或是广义梯度近似(GGA)的版本。
CASTEP 计算理论总结XBAPRSCASTEP 特点是适合于计算周期性结构,对于非周期性结构一般要将特定的部分作为周期性结构,建立单位晶胞后方可进行计算。
CASTEP 计算步骤可以概括为三步:首先建立周期性的目标物质的晶体;其次对建立的结构进行优化,这包括体系电子能量的最小化和几何结构稳定化。
最后是计算要求的性质,如电子密度分布(Electron density distribution),能带结构(Band structure)、状态密度分布(Densityof states)、声子能谱(Phonon spectrum)、声子状态密度分布(DOS of phonon),轨道群分布(Orbitalpopulations)以及光学性质(Optical properties)等。
本文主要将就各个步骤中的计算原理进行阐述,并结合作者对计算实践经验,在文章最后给出了几个计算事例,以备参考。
CASTEP 计算总体上是基于DFT ,但实现运算具体理论有:离子实与价电子之间相互作用采用赝势来表示;超晶胞的周期性边界条件;平面波基组描述体系电子波函数;广泛采用快速fast Fourier transform (FFT) 对体系哈密顿量进行数值化计算;体系电子自恰能量最小化采用迭带计算的方式;采用最普遍使用的交换-相关泛函实现DFT 的计算,泛函含概了精确形式和屏蔽形式。
一, CASTEP 中周期性结构计算优点与MS 中其他计算包不同,非周期性结构在CASTEP 中不能进行计算。
将晶面或非周期性结构置于一个有限长度空间方盒中,按照周期性结构来处理,周期性空间方盒形状没有限制。
之所以采用周期性结构原因在于:依据Bloch 定理,周期性结构中每个电子波函数可以表示为一个波函数与晶体周期部分乘积的形式。
他们可以用以晶体倒易点阵矢量为波矢一系列分离平面波函数来展开。
这样每个电子波函数就是平面波和,但最主要的是可以极大简化Kohn-Sham 方程。
CASTEP的原理及其应用1. 概述CASTEP(Computer Assisted Software for Theoretical and Experimental Phases)是一种用于计算材料性质的软件程序,基于第一性原理密度泛函理论(DFT)和平面波基组。
CASTEP具有广泛的应用领域,包括材料科学、固体物理、化学、地球科学等。
2. CASTEP的原理2.1 第一性原理密度泛函理论第一性原理密度泛函理论是一种计算量子力学的方法,基于电子的波函数和粒子的哈密顿量,通过求解波函数的动态方程来描述系统的行为。
在DFT中,体系的能量可以写为电子密度的函数,通过最小化总能量来得到系统的平衡态。
2.2 平面波基组在CASTEP中,平面波基组被用来展开电子波函数。
平面波基组具有周期性的性质,适用于固体材料的计算。
通过使用平面波基组,CASTEP可以对固体材料的电子结构进行准确的描述。
2.3 有效核势在DFT中,电子-电子相互作用被描述为有效核势。
CASTEP使用广义擦除近似(generalized gradient approximation,GGA)来近似描述核与电子之间的交换关联作用,以提高计算效率。
2.4 能带结构能带结构是材料中电子能量与波矢的关系。
CASTEP可以通过计算能带结构来研究材料的导电性、电子输运等性质。
2.5 势能面搜索CASTEP采用有效的势能面搜索算法,通过改变原子位置来寻找系统的平衡结构。
这对于计算晶体结构、力学性质、相变等有很大的意义。
3. CASTEP的应用3.1 材料设计与发现CASTEP在材料设计与发现方面有着广泛的应用。
它可以计算材料的晶体结构、电子结构、能带结构等,为材料研究提供基础数据。
通过计算不同元素、不同结构的材料的电子结构和能带结构,可以预测材料的性质,为新材料的设计提供指导。
3.2 材料性质预测除了材料设计与发现,CASTEP还可以用于预测材料的性质。
CASTEP 实战守则模型选定我们在进行材料物理模拟所需要做的第一步 (也是很重要的一步) 是模型的选定或建构。
CASTEP 虽然内建了很多功能来预测晶胞参数 (边长,夹角) 与原子位置,但仍然仰赖使用者告诉它 "要进行计算的系统是什么"。
在选定模型时,我们需切记如果系统内原子太多或是超晶胞体积太大,则计算量都会以平面波数的 3 次方增加到计算机难以负荷或使用者难以等待的程度,]此,我时时应考虑设计出一个足以表现出我们所想要研究的物理,而却又能使所有采用的超晶胞越少越好的模型。
如果研究上涉及到一系列原子数不同的大小系统需要做计算,最好能先计算小结构,不要一开始就送入大结构到计算机中。
如果要仿真的系统是含有杂质,则单位晶胞必须进一步放大成超晶胞以便使化学成分里的分数变成整数,因此晶胞会变得很大。
在某些特殊的情况,相互取代的元素种类是很类似的,(即在化学行为上类似),则下一个版本的CASTEP会提供一种叫做虚拟晶体近似 (Virtural Crystal Approximation,VCA) 的方法,则模型里面的原子就可以指定成如0.3A元素加0.7B元素这种样子,因此总可以以最小晶胞来做计算的模型。
但这种方法的精确度通常只适用于合金材料,故要小心使用。
模型的选定有许多人为抉择会含在里面,例如表面计算的层数,因此有些情况也需要进行所谓的收敛性测试。
Vps (poseudopotential)选择Vps选单:MS接口的设定是选用USP优先于NCP,USP有加速计算与减少内存使用的效果,其精确度也与NCP (norm-conserving potential) 相当。
至于什么时候使用NPC 呢 (就是在poseudopotential 选项中那些延伸文件名是 .recpot者),使用到NCP的场合有:1.某些CASTEP计算的功能尚未支持到USP,因此需要选用NCP。
2.为了要与已经发展的文献比较或进行验证3.对计算的结果存疑者,能提供『多一种选择』( 注:至于延伸档名是*.psp者也是属于norm-conserving的一种,是TM potential,在文献上也常被使用,但所需的截止动能较高,因此计算代价较大。