当前位置:文档之家› 第一性计算原理

第一性计算原理

第一性计算原理
第一性计算原理

Vasp

我所用第一原理是基于密度泛函(DFT)的从头计算,是以电子密度作为基本变量(HK定理),通过求解kohn-sham方程,迭代自洽得到体系的基态电子密度,然后求体系的基态性质。还有一种是基于hartree-fock自洽计算,通过自洽求解HF方程,获得体系的波函数,求基态性质。KS方程的计算水平达到了HF水平,同时还考虑了电子间的交换关联作用。关于DFT中密度泛函的Function其实是交换关联泛函,包括LDA,GGA,杂化泛函等等。一般LDA为局域密度近似,在空间某点用均匀电子气密度作为交换关联泛函的唯一变量,多数为参数化的CA-PZ方案;GGA为广义梯度近似,不仅将电子密度作为交换关联泛函的变量,也考虑了密度的梯度为变量,包括PBE,PE.RPBE等方案。

在处理计算体系中原子的电子态时有两种方法,一种是考虑所有电子叫做全电子法,比如WIEN2K中的FLAPW方法(线性缀加平面波);另一种是只考虑价电子而把芯电子和原子核构成离子实放在一起考虑即赝势法,一般贋势法是选取一个截断半径,截断半径以内波函数变化较平滑,和真实的不同,截断半径以外则和真实情况相同,而且贋势法得到的本征值和全电子法应该相同。贋势的测试标准应是贋势与全电子法计算结果的匹配度,而不是贋势与实验结果的匹配度,因为和实验结果的匹配可能是偶然的。

关于Ecut的收敛测试。一般情况下,总能相对于不同Ecut做计算,当截断能增大时总能变化不明显即可。但是在需要考虑体系应力时,还需要对应力进行收敛测试,而且应力相对于截断能要比总能更为苛刻。也就是某个截断能下总能已经收敛了,但应力未必收敛。(力的计算是在能量的基础上进行的,能量对坐标的一阶导数得到力。计算量的增大和误差的传递导致力收敛慢。)

K点也是需要经过测试的。

何时需要考虑自旋?例如BaTiO3中,三个元素分别为=+2,+4,-2价,离子全部为各个轨道满壳层的结构,此时就不必考虑自旋了。对于BaMnO3中,由于Mn+4价时d轨道还有电子但未满,因此需要考虑Mn(4s23d5)的自旋,Ba和O就不必考虑。其实设定自旋就是给定一个原子磁矩的初始值,只在刚开始计算时作为初始值使用。

几何优化包括晶格常数和原子位置的优化,一般情况下也有不优化几何结构直接计算电子结构的,但是对于缺陷形成的计算则往往要优化。

软件大致分为基于平面波的软件,如CASTEP,PWSCF.ABINIT等,计算量大概和体系原子数目的三次方相关;还有基于原子轨道线性组合的软件,比如openmx等,计算量和体系原子数目相关,一般可模拟较多原子数目的体系。

V ASP是使用贋势和平面波基组,进行从头量子力学分子动力学计算的软件包。V ASP中的方法基于有限温度下的局域密度近似(用自由能作为变量)以及对每一MD步骤用有效矩阵对角方案和有效Pulay混合求解瞬时电子基态。这些技术可以避免元氏的Car-Parrinello 方法存在的一切问题,而后者是基于电子、离子运动方程同时积分的方法。离子和电子的相互作用超缓Vanderbilt贋势(US-PP)或投影扩充波(PAW)方法描述。两种技术都可以相当程度地减少过度金属或第一行元素的每个原子所必须的平面波数量。V ASP可以很容易地计算力与张力,用于把原子衰减到其瞬时基态中。!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

V ASP程序亮点:

1、使用PAW方法或超软贋势,因此基组尺寸非常小,描述材料一般需要原子不超过100

个平面波,大多数情况下甚至每原子50个平面波就能得到可靠结果。

2、2. 在平面波程序中,某些部分代码的执行是三次标度。在VASP中,三次标度部分的前

因子足可忽略,导致关于体系尺寸的高效标度。因此可以在实空间求解势的非局域贡献,

并使正交化的次数最少。当体系具有大约2000个电子能带时,三次标度部分与其它部分可比,因此V ASP可用于直到4000个价电子的体系。

3、V ASP使用传统的自洽场循环计算电子基态。这一方案与数值方法组合会实现有效、稳

定、快速的Kohn-Sham方程自洽求解方案。程序使用的迭代矩阵对角化方案(RMM-DISS 和分块Davidson)可能是目前最快的方案。

4、V ASP包含全功能的对称性代码,可以自动确定任意构型的对称性。

5. 对称性代码还用于设定Monkhorst-Pack特殊点,可以有效计算体材料和对称的团簇。

Brillouin区的积分使用模糊方法或四面体方法。四面体方法可以用Bl?chl校正去掉线性四面体方法的二次误差,实现更快的k点收敛速度。

1、VASP能够进行哪些过程的计算?怎样设置?

我们平时最常用的研究方法是做单点能计算,结构优化、从头计算的分子动力学和电子结构相关性质的计算。

一般我们的研究可以按照这样的过程来进行

如果要研究一个体系的最优化构型问题可以首先进行结构弛豫优化,然后对优化后的结构进行性质计算或者单点能计算。

如果要研究一个体系的热力学变化过程可以首先进行分子动力学过程模拟,然后在某个温度或压强下进行性质计算或者单点能计算。

如果要研究一个体系的热力学结构变化可以首先在初始温度下进行NVT计算,然后进行分子动力学退火,然后在结束温度下进行性质计算研究。

2、什么是单点能计算(single point energy)?如何计算?

跟其它软件类似,VASP具有单点能计算的功能。也就是说,对一个给定的固定不变的结构(包括原子、分子、表面或体材料)能够计算其总能,即静态计算功能。

单点能计算需要的参数最少,最多只要在KPOINTS文件中设置一下合适的K点或者在INCAR文件中给定一个截断能ENCUT就可以了。还有一个参数就是电子步的收敛标准的设置EDIFF,默认值为EDIFF=1E-4,一般不需要修改这个值。

具体来说要计算单点能,只要在INCAR中设置IBRION=-1也就是让离子不移动就可以了。

3、什么是结构优化(structure optimization)?如何计算?

结构优化又叫结构弛豫(structure relax),是指通过对体系的坐标进行调整,使得其能量或内力达到最小的过程,与动力学退火不同,它是一种在0K下用原子间静力进行优化的方法。可以认为结构优化后的结构是相对稳定的基态结构,能够在实验之中获得的几率要大些(当然这只是理论计算的结果,必须由实验来验证)。

一般要做弛豫计算,需要设置弛豫收敛标准,也就是告诉系统收敛达成的判据(convergence break condition),当系统检测到能量变化减小到一个确定值时例如EDIFFG=1E-3时视为收敛中断计算,移动离子位置尝试进行下一步计算。EDIFFG这个值可以为负,例如EDIFFG=-0.02,这时的收敛标准是当系统发现所有离子间作用力都小于给定的数值,如0.02eV/A时视为收敛而中断。

弛豫计算主要有两种方式:准牛顿方法(quasi-Newton RMM-DIIS)和共轭梯度法(CG)两种。准牛顿方法计算速度较快,适合于初始结构与平衡结构(势能面上全局最小值)比较接近的情况,而CG方法慢一些,找到全局最小的可能性也要大一些。选择方法为IBRION=1时为准牛顿方法而IBRION=2时为CG方法。

具体来说要做弛豫计算,设置IBRION=1或者2就可以了,其它参数根据需要来设置。NSW是进行弛豫的最大步数,例如设置NSW=100,当计算在100步之内达到收敛时计算自动中断,而100步内没有达到收敛的话系统将在第100步后强制中止(平常计算步数不会超过100步,超过100步可能是计算的体系出了问题)。参数通常可以从文献中发现,例如收敛标准EDIFFG等。

有的时候我们需要一些带限制条件的弛豫计算,例如冻结部分原子、限制自旋的计算等等。冻结部分原子可以在POSCAR文件中设置selective dynamic 来实现。自旋多重度限制可以在INCAR中以NUPDOWN选项来设置。另外ISIF 选项可以控制弛豫时的晶胞变化情况,例如晶胞的形状和体积等。

费米面附近能级电子分布的smearing是一种促进收敛的有效方法,可能产生物理意义不明确的分数占据态情况,不过问题不大。在INCAR文件中以ISMEAR 来设置。一般来说K点只有一两个的时候采用ISMEAR=0,金属体材料用ISMEAR=1或2,半导体材料用ISMEAR=-5等等。不过有时电子步收敛速度依然很慢,还需要设置一些算法控制选项,例如设置ALGO=Very_Fast,减小真空层厚度,减少K点数目等。

弛豫是一种非常有效的分析计算手段,虽然是静力学计算但是往往获得一些动力学得不到的结果。

INCAR:

EDIFF 一般来说,用1E-4或者1E-5都可以,这个参数只是对第一个离子步的自洽影响大一些,对于长时间的分子动力学的模拟,精度小一点也无所谓,但不能太小。

IBRION=0分子动力学模拟

IALGO=48一般用48,对于原子数较多,这个优化方式较好。

NSW=1000 多少个时间步长。

POTIM=3 时间步长,单位fs,通常1到3.

ISIF=2 计算外界的压力.

NBLOCK= 1 多少个时间步长,写一次CONTCAR,CHG和CHGCAR,PCDAT. KBLOCK=50 NBLOCK*KBLOCK个步长写一次XDATCAR.(个离子步写一次PCDAT.)

ISMEAR=-1 费米迪拉克分布.

SIGMA =0.05单位:电子伏

NELMIN=8 一般用6到8,最小的电子scf数.太少的话,收敛的不好.

LREAL=A

APACO=10径向分布函数距离,单位是埃.

NPACO=200 径向分布函数插的点数.

LCHARG=F尽量不写电荷密度,否则CHG文件太大.

TEBEG=300 初始温度.

TEEND=300终态温度。不设的话,等于TEBEG.

SMASS=-3 NVE ensemble;-1用来做模拟退火。大于0 NVT系综。正确:SMASS=1,2,3是没有区别的。都是NVT ensemble。SMASS只要是大于0就是NVT系综。

CONTCAR是每个离子步之后都会写出来的,但是会用新的把老的覆盖

CHG是在每10个离子步写一次,不会覆盖

CHGCAR是在任务正常结束之后才写的。

5、收敛判据的选择

结构弛豫的判据一般有两中选择:能量和力。这两者是相关的,理想情况下,能量收敛到基态,力也应该是收敛到平衡态的。但是数值计算过程上的差异导致以二者为判据的收敛速度差异很大,力收敛速度绝大部分情况下都慢于能量收敛速度。这是因为力的计算是在能量的基础上进行的,能量对坐标的一阶导数得到力。计算量的增大和误差的传递导致力收敛慢。

到底是以能量为收敛判据,还是以力为收敛判据呢?关心能量的人,觉得以能量为判据就够了;关心力相关量的人,没有选择,只能用力作为收敛标准。对于超胞体系的结构优化,文献大部分采用Gamma点做单点优化。这个时候即使采用力为判据(EDIFFG=-0.02),在做静态自洽计算能量的时候,会发现,原本已经收敛得好好的力在不少敏感位置还是超过了结构优化时设置的标准。这个时候,是不是该怀疑对超胞仅做Gamma点结构优化的合理性呢?是不是要提高K点密度再做结构优化呢。

在我看来,这取决于所研究的问题的复杂程度。我们的计算从原胞开始,到超胞,到掺杂结构,到吸附结构,到反应和解离。每一步都在增加复杂程度。结构优化终点与初始结构是有关的,如果遇到对初始结构敏感的优化,那就头疼了。而且,还要注意到,催化反应不仅与原子本身及其化学环境有关,还会与几何构

型有关。气固催化反应过程是电子的传递过程,也是分子拆分与重新组合的过程。如果优化终点的构型不同,可能会导致化学反应的途径上的差异。仅从这一点来看,第一性原理计算的复杂性,结果上的合理性判断都不是手册上写的那么简单。

对于涉及构型敏感性的结构优化过程,我觉得,以力作为收敛判据更合适。而且需要在Gamma点优化的基础上再提高K点密度继续优化,直到静态自洽计算时力达到收敛标准的。

6、结构优化参数设置

结构优化,或者叫弛豫,是后续计算的基础。其收敛性受两个主要因素影响:初始结构的合理性和弛豫参数的设置

初始结构

初始结构包括原子堆积方式,和自旋、磁性、电荷、偶极等具有明确物理意义的模型相关参数。比如掺杂,表面吸附,空位等结构,初始原子的距离,角度等的设置需要有一定的经验积累。DFT计算短程强相互作用(相对于范德华力),如果初始距离设置过远(如超过4埃),则明显导致收敛很慢甚至得到不合理的结果。

比较好的设置方法可以参照键长。比如CO在O顶位的吸附,可以参照CO2中C-O键长来设置(如增长20%)。也可以参照文献。记住一些常见键长,典型晶体中原子间距离等参数,有助于提高初始结构设置的合理性。实在不行,可以先在小体系上测试,然后再放到大体系中算。

弛豫参数

弛豫参数对收敛速度影响很大,这一点在计算工作没有全部铺开时可能不会觉察到有什么不妥,反正就给NSW设置个“无穷大”的数,最后总会有结果的。但是,时间是宝贵的,恰当的设置3小时就收敛的结果,不恰当的设置可能要一个白天加一个黑夜。如果你赶文章或者赶着毕业,你就知道这意味这什么。

结构优化分电子迭代和离子弛豫两个嵌套的过程。电子迭代自洽的速度,有四个响很大的因素:初始结构的合理性,k点密度,是否考虑自旋和高斯展宽(SIGMA);离子弛豫的收敛速度,有三个很大的影响因素:弛豫方法(IBRION),步长(POTIM)和收敛判据(EDIFFG).

一般来说,针对理论催化的计算,初始结构都是不太合理的。因此一开始采用很粗糙的优化(EDIFF=0.001,EDIFFG=-0.2),很低的K点密度(Gamma),不考虑自旋就可以了,这样NSW<60的设置就比较好。其它参数可以默认。

经过第一轮优化,就可以进入下一步细致的优化了。就我的经验,

EDIFF=1E-4,EDIFFG=-0.05,不考虑自旋,IBRION=2,其它默认,NSW=100;

跑完后可以设置IBRION = 1,减小OPTIM(默认为0.5,可以设置0.2)继续优化。

优化的时候让它自己闷头跑是不对的,经常看看中间过程,根据情况调节优化参数是可以很好的提高优化速度。这个时候,提交两个以上的任务排队是好的方式,一个在调整的时候,下一个可以接着运行,不会因为停下当前任务导致机器空闲。

无论结构优化还是静态自洽,电子步的收敛也常常让新手头痛。如果电子步不能在40步内收敛,要么是参数设置的问题,要么是初始模型太糟糕(糟糕的不是一点点)。

静态自洽过程电子步不收敛一般是参数设置有问题。这个时候,改变迭代算法(ALGO),提高高斯展宽(SIGMA增加),设置自洽延迟(NELMDL)都是不错的方法。对于大体系比较难收敛的话,可以先调节AMIN,BMIX跑十多步,得到电荷密度和波函数,再重新计算。实在没办法了,可以先放任它跑40步,没有收敛的迹象的话,停下来,得到电荷密度和波函数后重新计算。一般都能在40步内收敛。

对于离子弛豫过程,不调节关系也不大。开始两个离子步可能要跑满60步(默认的),后面就会越来越快了。

总的说来,一般入门者,多看手册,多想多理解,多上机实践总结,比较容易提高到一个熟练操作工的水平。

如果要想做到“精确打击”,做到能在问题始发的时候就立刻采取有效措施来解决,就需要回归基础理论和计算方法上来了。

7、优化结果对初始结构和“优化路径”的依赖

原子吸附问题不大,但是小分子吸附,存在初始构型上的差异。slab上水平放置,还是垂直放置,可能导致收敛结果上的差异。根据H-K理论,理想情况下,优化得到的应该是全局最小,但在数值计算的时候可能经常碰到不是全局最小的情况。实际操作中发现,多个不同初始结构优化收敛后在能量和结构上存在一定差异。

为了加快收敛速度,特别是对于表面-分子吸附结构,初始放松约束,比如EDIFF=1E-3,EDIFFG=-0.3,NSW=30可能是很好的设置。但是下面的情况应当慎重:

EDIFF=1E-3;

EDIFFG=-0.1;!或者更小

NSW=500;!或者更大

电子步收敛约束较小,而离子步约束偏大,离子步数又很多,这种情况下,可能导致的结果是结构弛豫到严重未知的区间。

再在这个基础上提高约束来优化,可能就是徒劳的了——结果不可逆转的偏向不正常的区间。

好的做法,是对初始结构做比较松弛的约束,弛豫离子步NSW应该限制在一个较小的数值内。EDIFF=1E-3的话,EDIFFG也最好是偏大一些,如-0.3而不是-0.1.这样可以在较少的步数内达到初步收敛。

对于远离基态的初始结构,一开始在非常松弛的约束下跑若干离子步,时间上带来的好处是很大的。对于100个原子的体系用vasp做Gamma点优化,如果一开始就是正常优化(EDIFF=1E-4,EDIFFG=-0.02)设置,开始十个离子步可能都要花上几个小时。如果这个时候才发现输入文件有错误,那下午的时间就白费了,顺便带上晚上机器空转。

所以,我习惯的做法,是在初始几步优化后,会用xcrysden检查一下XDATCAR中的数据,用xdat2xyz.pl生成movie.xyz,然后看看弛豫过程是不是按照设想的那样。后续过程跑完一个收敛过程,就再检查一下movie.xyz。如此这般,才放心的展开后续计算。

8、目的导向的结构优化

结构优化到这个阶段,是高级的了。为了得到特定结构,或者为了验证某些猜想,需要设计合理的初始结构,然后在这个基础上小心优化,比如POTIM=0.1跑几步看看,然后修改优化参数。

我遇到过的一件跟结构优化关系很大的算例是CeO2氧空位结构电子局域的问题(https://www.doczj.com/doc/6710309853.html,/bbs/viewthread.php?tid=2954558)。按照一般方式(从优化好的bulk建slab模型,然后优化)得到一个O空位留下的两个电子均匀局域到O 次外层三个Ce原子上,得到空位形成能2.34eV.经高人指点后,调节空位附近O 原子位置,打破对称性后重新优化,两个电子完美的局域到两个Ce原子上了。并且空位形成能降低到2.0X eV。从这个例子可以看到,结构优化存在不少技巧的,这些技巧建立在研究者对模拟对象的物理意义的理解上。对物理图像的直观深入理解,才能做好模型预设,在此引导下才可能有目的的优化出不比寻常的结果。

目前第一性原理理论中的交换关联泛函部分包含经验参数。考虑这一点对优化结果的影响也很有意思。比如有专家提到,DFT+U参数对某些结构的收敛终态构型有影响。构型的变化可能影响表面反应过程。基于这一点,一个好的计算研究可能就出来了。

真实过程总是复杂多变的。无论何种模拟,估计都可以找到一些试验现象来验证。但是到底应该如何评判模拟结果,如何从第一性原理研究中得出有意义的结论需要很好的洞察力。这样的模拟不见得就必须建立的试验的基础上,完全凭空设计的模型有可能更能优美的解释本质。

第一性原理计算原理和方法

第二章 计算方法及其基本原理介绍 化学反应的本质就是旧键的断裂与新建的形成,参与成键原子的电子壳层重新组合就是导致生成稳定多原子化学键的明显特征。因此阐述化学键的理论应当描写电子壳层的相互作用与重排,借助求解满足适当的Schrodinger 方程的波函数描写分子中电子分布的量子力学,为解决这一问题提供了一般的方法,然而,对于一些实际的体系,不引入一些近似,就不可能求解其Schrodinger 方程。这些近似使一般量子力学方程简化为现代电子计算机可以求解的方程。这些近似与关于分子波函数的方程形成计算量子化学的数学基础。 2、1 SCF-MO 方法的基本原理 分子轨道的自洽场计算方法 (SCF-MO)就是各种计算方法的理论基础与核心部分,因此在介绍本文计算工作所用方法之前,有必要对其关键的部分作一简要阐述。 2、1、1 Schrodinger 方程及一些基本近似 为了后面介绍各种具体在自洽场分子轨道(SCF MO)方法方便,这里将主要阐明用于本文量子化学计算的一些重要的基本近似,给出SCF MO 方法的一些基本方程,并对这些方程作简略说明,因为在大量的文献与教材中对这些方程已有系统的推导与阐述[1-5]。 确定任何一个分子的可能稳定状态的电子结构与性质,在非相对论近似下,须求解 R AB =R 图2-1分子体系的坐标

定态Schrodinger 方程 ''12121212122ψψT p B A q p A p pA A pq AB B A p A A A E R Z r R Z Z M =??????? ?-++?-?-∑∑∑∑∑∑≠≠ (2、1) 其中分子波函数依赖于电子与原子核的坐标,Hamilton 算符包含了电子p 的动能与电子p 与q 的静电排斥算符, ∑∑≠+?-=p q p pq p e r H 12121?2 (2、2) 以及原子核的动能 ∑?-=A A A N M H 2121? (2、3) 与电子与核的相互作用及核排斥能 ∑∑≠+-=p A B A AB B A pA A eN R Z Z r Z H ,21? (2、4) 式中Z A 与M A 就是原子核A 的电荷与质量,r pq =|r p -r q |,r pA =|r p -R A |与R AB =|R A -R B |分别就是电子p 与q 、核A 与电子p 及核A 与B 间的距离(均以原子单位表示之)。上述分子坐标系如图2、1所示。可以用V(R,r)代表(2、2)-(2、4)式中所有位能项之与 ∑∑∑-+=≠≠p A pA A B A q p pq AB B A r Z r R Z Z r R V ,1 2121),( (2、5) 原子单位 上述的Schrodinger 方程与Hamilton 算符就是以原子单位表示的,这样表示的优点在于简化书写型式与避免不必要的常数重复计算。在原子单位的表示中,长度的原子单位就是Bohr 半径

第一性原理简介

第一性原理是什么 第一性原理怎么用 1什么是第一性原理 根据原子核和电子互相作用的原理及其基本运动规律,运用,从具体要求出发,经过一些近似处理后直接求解的算法,称为第一性原理。广义 的第一原理包括两大类,以Hartree-Fock自洽场计算为基础的从头算和 (DFT计算。 从定义可以看出第一性原理涉及到量子力学、、Hartree-Fock自洽场、等许多对我来说很陌生的物理化学定义。因此我通过向师兄请教和上网查资料一点点的了解并学习这些知识。 2第一性原理的作用 以密度泛函理论(DFT)为基础以及在此基础上发展起来的简单而具有一定精度的局域密度近似(LDA)和广义梯度近似(GGA)的第一性原理电子结构计算方法,与传统的解析方法一样,不但能够给出描述体系微观电子特性的物理量如波函数、态密度、费米面、电子间互作用势等,以及在此基础上所得到的体现体系宏观物理特性的参量如结合能、电离能、比热、电导、光电子谱、穆斯堡尔谱等等,而且它还可以帮助人们预言许多新的

物理现象和物理规律。密度泛函计算的一些结果能够与实验直接进行比较一些应用程序的发展乃至商业软件的发布,导致了基于密度泛函理论的第 一原理计算方法的广泛应用。 密度泛函理论(DFT)为第一性原理中的一类,在物理系、化学、材料科学以及其他工程领域中,密度泛函理论(DFT及其计算已经快速发展成 为材料建模模拟的一种“标准工具”。 密度泛函理论可以计算预测固体的晶体结构、晶格参数、能带结构、态密度(DOS、光学性能、磁性能以及原子集合的总能等等。 3第一性原理怎么用 目前我所学到的利用第一性原理的软件为Material Studio 、VASP软件。其中Materials Studio (简称MS是专门为材料科学领域研究者幵发的一款可运行在PC上的模拟软件。使化学及材料科学的研究者们能更方便地建立三维结构模型,并对各种晶体、无定型以及高分子材料的性质及相关过程进行深入的研究。模拟的内容包括了催化剂、聚合物、固体及表面、晶体与衍射、化学反应等材料和化学研究领域的主要课题。 模块简介 Materials Studio 采用了大家非常熟悉的Microsoft标准用户界面, 允许用户通过各种控制面板直接对计算参数和计算结果进行设置和分析。 目前,Materials Studio 软件包括如下功能模块: Materials Visualizer: 提供了搭建分子、晶体及高分子材料结构模型所需要的所有工具,可以操作、观察及分析结构模型,处理图表、表格或文本等形式的数据,并提供软件的基本环境和分析工具以及支持Materials Studio 的其他产品。是Materials Studio 产品系列的核心模块。 Discover: Materials Studio 的分子力学计算引擎。使用多种分子力学和动力学 方法,以仔细推导的力场作为基础,可准确地计算出最低能量构型、分子体系的结构和动力学轨迹等。

第一性原理计算方法讲义

第一性原理计算方法讲 义 标准化管理部编码-[99968T-6889628-J68568-1689N]

第一性原理计算方法 引言 前面讲述的有限元和有限差分等数值计算方法中,求解的过程中需要知道一些物理参量,如温度场方程中的热传导系数和浓度场方程中的扩散系数等,这些参量随着材料的不同而改变,需要通过实验或经验来确定,所以这些方法也叫做经验或者半经验方法。而第一性原理计算方法只需要知道几个基本的物理参量如电子质量、电子的电量、原子的质量、原子的核电荷数、布朗克常数、波尔半径等,而不需要知道那些经验或半经验的参数。第一性原理计算方法的理论基础是量子力学,即对体系薛定额方程的求解。 量子力学是反映微观粒子运动规律的理论。量子力学的出现,使得人们对于物质微观结构的认识日益深入。原则上,量子力学完全可以解释原子之间是如何相互作用从而构成固体的。量子力学在物理、化学、材料、生物以及许多现代技术中得到了广泛的应用。以量子力学为基础而发展起来的固体物理学,使人们搞清了“为什么物质有半导体、导体、绝缘体的区别”等一系列基本问题,引发了通讯技术和计算机技术的重大变革。目前,结合高速发展的计算机技术建立起来的计算材料科学已经在材料设计、物性研究方面发挥着越来越重要的作用。 但是固体是具有~1023数量级粒子的多粒子系统,具体应用量子理论时会导致物理方程过于复杂以至于无法求解,所以将量子理论应用于固体系统必须采用一些近似和简化。绝热近似(Born-Oppenheimei近似)将电子的运动和原子核的运动分开,从而将多粒子系统简化为多电子系统。Hartree-Fock近似将多电子问题简化为仅与以单电子波函数(分子轨道)为基本变量的单粒子问题。但是其中波函数的行列式表示使得求解需要非常大的计算量;对于研究分子体系,他可以作为一个很好的出发点,但是不适于研究固态体系。1964年,Hohenberg和Kohn提出了严格的密度泛函理论(Density Functional Theory, DFT)。它建立在非均匀电子气理论基础之上,以粒子数密度()r 作为基本变量。1965年,Kohn和Sham提出Kohn-Sham方程将复杂的多电子问题及其对应的薛定谔方程转化为相对简单的单电子问题及单电子Kohn-Sham方程。将精确的密度泛函理论应用到实际,需要对电子间的交换关联作用进行近似。局域密度近似(LDA)、广义梯度近似(GGA)等的提出,以及以密度泛函理论为基础的计算方法(赝

第一性原理计算原理和方法精编

第一性原理计算原理和 方法精编 Document number:WTT-LKK-GBB-08921-EIGG-22986

第二章 计算方法及其基本原理介绍 化学反应的本质是旧键的断裂和新建的形成,参与成键原子的电子壳层重新组合是导致生成稳定多原子化学键的明显特征。因此阐述化学键的理论应当描写电子壳层的相互作用与重排,借助求解满足适当的Schrodinger 方程的波函数描写分子中电子分布的量子力学,为解决这一问题提供了一般的方法,然而,对于一些实际的体系,不引入一些近似,就不可能求解其Schrodinger 方程。这些近似使一般量子力学方程简化为现代电子计算机可以求解的方程。这些近似和关于分子波函数的方程形成计算量子化学的数学基础。 SCF-MO 方法的基本原理 分子轨道的自洽场计算方 法(SCF-MO)是各种计算方法的理论基础和核心部分,因此在介绍本文计算工作所用方法之 前,有必要对其关键的部分作 一简要阐述。 Schrodinger 方程及一些基本近似 为了后面介绍各种具体在自洽场分子轨道(SCF MO)方法方便,这里将主要阐明用于本文量子化学计算的一些重要的基本 R AB =R 图2-1分子体系的坐标

近似,给出SCF MO 方法的一些基本方程,并对这些方程作简略说明,因为在大量的文献和教材中对这些方程已有系统的推导和阐述[1-5]。 确定任何一个分子的可能稳定状态的电子结构和性质,在非相对论近似下,须求解定态Schrodinger 方程 ''12121212122ψψT p B A q p A p pA A pq AB B A p A A A E R Z r R Z Z M =??????? ?-++?-?-∑∑∑∑∑∑≠≠ () 其中分子波函数依赖于电子和原子核的坐标,Hamilton 算符包含了电子p 的动能和电子p 与q 的静电排斥算符, ∑∑≠+?-=p q p pq p e r H 12121?2 以及原子核的动能 ∑?-=A A A N M H 2121? 和电子与核的相互作用及核排斥能 ∑∑≠+-=p A B A AB B A pA A eN R Z Z r Z H ,21? 式中Z A 和M A 是原子核A 的电荷和质量,r pq =|r p -r q |,r pA =|r p -R A |和R AB =|R A -R B |分别是电子p 和q 、核A 和电子p 及核A 和B 间的距离(均以原子单位表示之)。上述分子坐标系如图所示。可以用V(R,r)代表-式中所有位能项之和 ∑∑∑-+=≠≠p A pA A B A q p pq AB B A r Z r R Z Z r R V ,12121),( 原子单位

如何分析能带图及第一性原理的计算

分析能带图 能带结构是目前采用第一性原理(从头abinitio)计算所得到的常用信息,可用来结合解释金属、半导体和绝缘体的区别。能带可分为价带、禁带和导带三部分,倒带和价带之间的空隙称为能隙,基本概念如图所示: 如何能隙很小或为0 ,则固体为金属材料,在室温下电子很容易获得能量而跳跃至传倒带而导电;而绝缘材料则因为能隙很大(通常大于9电子伏特),电子很难跳跃至传导带,所以无法导电。一般半导体材料的能隙约为1至3电子伏特,介于导体和绝缘体之间。因此只要给予适当条件的能量激发,或是改变其能隙之间距,此材料距能导电。 能带用来定性地阐明了晶体中电子运动的普遍特点。价带(valence band),或称价电带,通常指绝对零度时,固体材料里电子的最高能量。在导带(conduction band)中,电子的能量范围高于价带,而所有在传导带中的电子均可经由外在的电

场加速而形成电流。对与半导体以及绝缘体而言,价带的上方有一个能隙(band gap),能隙上方的能带则是传导带,电子进入传导带后才能在固体材料内自由移动,形成电流。对金属而言,则没有能隙介于价带与传导带之间,因此价带是特指半导体与绝缘体的状况。 费米能级(fermi level)是绝对零度下的最高能级。根据泡利不相容原理,一个量 子态不能容纳两个或两个以上的费米子(电子),所以在绝度零度下,电子将从低到高依次填充各能级,除最高能级外均被填满,形成电子态的“费米海”。“费米海” 中每个电子的平均能量为(绝对零度下)为费米能级的3/5。海平面即是费米能级。一般来说,费米能级对应态密度为0的地方,但对于绝缘体而言,费米能级就位于价带顶。成为优良电子导体的先决条件是费米能级与一个或更多的能带相交。 能量色散(dispersion of energy)。同一个能带内之所以会有不同能量的量子态, 原因是能带的电子具有不同波向量(wave vector),或是k-向量。在量子力学中, k-向量即为粒子的动量,不同的材料会有不同的能量-动量关系(E-K relationship)。能量色散决定了半导体材料的能隙是直接能隙还是间接能隙。如导带最低点与价带最高点的K值相同,则为直接能隙,否则为间接能隙。 能带的宽度。能带的宽度或三度,即能带最高和最低能级之间的能量差,是一个非常重要的特征,它是由相互作用的轨道之间的重叠来决定的,因而反应出轨道之间的重叠情况,相邻的轨道之间重叠越大,带宽就越大。

第一性原理计算

实验一、第一性原理计算 1. 实验目的 (1) 掌握第一性原理和密度泛涵的计算方法; (2) 学会使用Visualizer 的各种建模和可视化工具; (3) 熟悉CASTEP 模块的功能。 2. 实验原理 CASTEP 是基于密度泛涵理论平面波赝势基础上的量子力学计算。 密度泛涵理论的基本思想是原子、分子和固体的基本物理性质可以用粒子密度函数进行描述。可以归纳为两个基本定理: 定理1:粒子数密度函数是一个决定系统基态物理性质的基本参量。 定理2:在粒子数不变的条件下能量对密度函数变分得到系统基态的能量。不计自旋的全同费米子的哈密顿量为:H T U V =++ 其中动能项为:()()T dr r r ψψ+=??? 库仑作用项为:11'()(')()(')2 ' U drdr r r r r r r ψψψψ++=-? V 为对所有粒子均相同的局域势u(r)表示的外场影响:()()()V dru r r r ψψ+=?粒子数密度函数为: ()()()r r r ρψψ+=ΦΦ 对于给定的()r υ,能量泛函[]E ρ定义为: []()()E dr r r T U ρυρ=+Φ+Φ ?;[]F T U ρ=Φ+Φ系统基态的能量: ' ''''[]''''[][]()()[][]()()[] E T U V G E F dr r r E G G F dr r r E G ρρυρφρυρρΦ=Φ+Φ+ΦΦ==+>?=+=? 3. 实验内容 材料的电子结构计算; 4. 实验设备和仪器 (1) 硬件:多台PC 机和一台高性能计算服务器。 软件:主要利用Materials studio 软件包里的Materials Visualizer 和CASTEP 模块 5. 实验步骤

第一节第一性原理计算方法综述

第一性原理计算的理论方法 随着科技的发展,计算机性能也得到了飞速的提高,人们对物理理论的认识也更加的深入,利用计算机模拟对材料进行设计已经成为现代科学研究不可缺少的研究手段。这主要是因为在许多情况下计算机模拟比实验更快、更省,还得意于计算机模拟可以预测一些当前实验水平难以达到的情况。然而在众多的模拟方法中,第一性原理计算凭借其独特的精度和无需经验参数而得到众多研究人员的青睐,成为计算材料学的重要基础和核心计算。本章将介绍第一性原理计算的理论基础,研究方法和ABINIT软件包。 1.1 第一性原理 第一性原理计算( 简称从头计算,the abinitio calculation) ,指 从所要研究的材料的原子组分出发,运用量子力学及其它物理规律,通过自洽计算来确定指定材料的几何结构、电子结构、热力学性质和光学性质等材料物性的方法。基本思想是将多原子构成的实际体系理解成为只有电子和原子核组成的多粒子系统,运用量子力学等最基本的物理原理最大限度的对问题进行”非经验”处理。【1】第一性原理计算就只需要用到五个最基本的物理常量即( m o.e.h.c.k b ) 和元素周期表中各组分元素的电子结构,就可以合理地预测材料的许多物理性质。用第一性原理计算的晶胞大小和实验值相比误差只有几个百分点,其他性质也和实验结果比较吻合,体现了该理论的正确性。

第一性原理计算按照如下三个基本假设把问题简化: 1.利用Born-Oppenheimer 绝热近似把包含原子核和电子的多粒子问题转化为多电子问题。 2.利用密度泛函理论的单电子近似把多电子薛定谔方程简化为比较容易求解的单电子方程。 3.利用自洽迭代法求解单电子方程得到系统基态和其他性质。以下我将简单介绍这些第一性原理计算的理论基础和实现方法:绝热近似、密度泛函理论、局域密度近似(LDA)和广义梯度近似(GGA)、平面波及赝势方法、密度泛函的微扰理论、热力学计算方法和第一性原理计算程序包ABINIT。 1.2量子力学与Born-Oppenheimer 近似固体是由原子核和核外的电子组成的,在原子核与电子之间,电子与电子之间,原子核与原子核之间都存在着相互作用。从物理学的角度来看,固体是一个多体的量子力学体系【2】,相应的体系哈密顿量可以写成如下形式: H (r,R) E H(r ,R) (1-1) 其中r,R 分别代表所有电子坐标的集合、所有原子核坐标的集合。在不计外场作用下,体系的哈密顿量日包括体系所有粒子( 原子核和电子) 的动能和粒子之间的相互作用能,即 H H e H N H e N (1-2) 其中,以是电子部分的哈密顿量,形式为: 22 1 e2 H e(r) r2i 1 e(1-3)

第一性原理计算原理和方法

第二章 计算方法及其基本原理介绍 化学反应的本质是旧键的断裂和新建的形成,参与成键原子的电子壳层重新组合是导致生成稳定多原子化学键的明显特征。因此阐述化学键的理论应当描写电子壳层的相互作用与重排,借助求解满足适当的Schrodinger 方程的波函数描写分子中电子分布的量子力学,为解决这一问题提供了一般的方法,然而,对于一些实际的体系,不引入一些近似, 确定任何一个分子的可能稳定状态的电子结构和性质,在非相对论近似下,须求解定态Schrodinger 方程 ''12121212122 ψψT p B A q p A p pA A pq AB B A p A A A E R Z r R Z Z M =??? ?????-++?-?-∑∑∑∑∑∑≠≠ (2.1) 其中分子波函数依赖于电子和原子核的坐标,Hamilton 算符包含了电子p 的动能和电子p

与q 的静电排斥算符, ∑∑≠+?-=p q p pq p e r H 12121?2 (2.2) 以及原子核的动能 ∑?-=A A A M H 2? (2.3) 和电子与核的相互作用及核排斥能 ∑∑≠+-=p A B A AB B A pA A eN R Z Z r Z H ,21? (2.4) 式中Z A 和M A 是原子核A 的电荷和质量,r pq =|r p -r q |,r pA =|r p -R A |和R AB =|R A -R B |分别是电子p 和q 、核A 和电子p 及核A 和B 间的距离(均以原子单位表示之)。上述分子坐标系如图2.1所示。可以用V(R,r)代表(2.2)-(2.4)式中所有位能项之和 ∑∑∑-+= ≠≠p A pA A B A q p pq AB B A r Z r R Z Z r R V ,1 2121),( (2.5) 原子单位 上述的Schrodinger 方程和Hamilton 算符是以原子单位表示的,这样表示的优点在于简化书写型式和避免不必要的常数重复计算。在原子单位的表示中,长度的原子单位是Bohr 半径 能量是以Hartree 为单位,它定义为相距1Bohr 的两个电子间的库仑排斥作用能 质量则以电子制单位表示之,即定义m e =1 。

第一性原理计算

钙钛矿型PbZrO3电子能带结构的第一性原理计算 班级:s1467 姓名:学号:201421801014 锆酸铅(PbZrO3)是最早发现的反铁电体之一,在工业上的一个重要应用是其固溶物Pb(Zr,Ti)O3。由于反铁电材料在相开关、电荷存储、电流源、电容、微电子及微型机电设备等方面有重要应用,其电子结构和物理特性一直为人们所关注。PbZrO3的有三个不同的相,在233℃以上为立方顺电相,具有钙钛矿结构,所属的空间群为Pm3m;当晶体处于233℃以下,将发生氧八面体的扭曲畸变和阳离子相对于O的移动,形成结构相变;230~233℃为正交铁电相,而230℃以下的基态为正交晶系,空间群为Pbam。基态正交相中离子移动主要由Pb、O之间的相对位移提供,由于相邻晶格之间Pb-O的位移相反,因此其为反铁电体。 1、原理及计算 采用第一性原理局域密度近似下的投影缀加平面波方法精确计算并比较了钙钛矿材料PbZrO3低温正交相(反铁电相)、高温立方相(顺电相)的电子能带结构,计算了PbZrO3材料正交相、立方相的电子结构。PbZrO3立方相的空间群为Pm3m,计算采用实验得到的晶格常量为a=4.11nm,Wyckoff坐标为Pb:(0,0,0),Zr:(0.5,0.5,0.5),O:(0.5,0.5,0)。正交相的空间群为Pmam,采用的晶格常数a=5.9411nm,b=11.8024nm,c=8.2564nm,各原子坐标见表1。正交相和立方相的多面体结构模型如图1所示。平面波截断能取为500eV,布里渊区积分分别采用5×5×5及7×3×5的K点网格,高斯展宽因子为0.1eV。 表1 正交相PbZrO3原胞内的原子位置

第一性原理计算方法讲义

第一性原理计算方法 引言 前面讲述的有限元和有限差分等数值计算方法中,求解的过程中需要知道一些物理参量,如温度场方程中的热传导系数和浓度场方程中的扩散系数等,这些参量随着材料的不同而改变,需要通过实验或经验来确定,所以这些方法也叫做经验或者半经验方法。而第一性原理计算方法只需要知道几个基本的物理参量如电子质量、电子的电量、原子的质量、原子的核电荷数、布朗克常数、波尔半径等,而不需要知道那些经验或半经验的参数。第一性原理计算方法的理论基础是量子力学,即对体系薛定额方程的求解。 量子力学是反映微观粒子运动规律的理论。量子力学的出现,使得人们对于物质微观结构的认识日益深入。原则上,量子力学完全可以解释原子之间是如何相互作用从而构成固体的。量子力学在物理、化学、材料、生物以及许多现代技术中得到了广泛的应用。以量子力学为基础而发展起来的固体物理学,使人们搞清了“为什么物质有半导体、导体、绝缘体的区别”等一系列基本问题,引发了通讯技术和计算机技术的重大变革。目前,结合高速发展的计算机技术建立起来的计算材料科学已经在材料设计、物性研究方面发挥着越来越重要的作用。 但是固体是具有?1023数量级粒子的多粒子系统,具体应用量子理论时会导致物理方程过于复杂以至于无法求解,所以将量子理论应用于固体系统必须采用一些近似和简化。绝热近似(Born-Oppenheimei 近似)将电子的运动和原子核的运动分开,从而将多粒子系统简化为多电子系统。Hartree-Fock 近似将多电子问题简化为仅与以单电子波函数(分子轨道)为基本变量的单粒子问题。但是其中波函数的行列式表示使得求解需要非常大的计算量;对于研究分子体系,他可以作为一个很好的出发点,但是不适于研究固态体系。1964年,Hohenberg和Kohn提出了严格的 密度泛函理论(Density Functional Theory, DFT )。它建立在非均匀电子气理论基础之上,以粒子数密度(『)作为基本变量。1965年,Kohn和Sham提出Kohn-Sham方程将复杂的多电子问题及其对应的薛定谔方程转化为相对简单的单电子问题及单电子Kohn-Sham方程。将精确的密度泛函理 论应用到实际,需要对电子间的交换关联作用进行近似。局域密度近似(LDA、广义梯度近似(GGA 等的提出,以及以密度泛函理论为基础的计算方法(赝势方法、全电子线形缀加平面波方法(FLAPW)等、的提出,使得密度泛函理论在化学和固体物理中的电子结构计算取得了广泛的应用,从而使得固体材料的研究取得长足的进步。 第一性原理计算方法的应用 1、体系的能量

第一节第一性原理计算方法综述

第一性原理计算的理论方法 随着科技的发展,计算机性能也得到了飞速的提高,人们对物理理论的认识也更加的深入,利用计算机模拟对材料进行设计已经成为现代科学研究不可缺少的研究手段。这主要是因为在许多情况下计算机模拟比实验更快、更省,还得意于计算机模拟可以预测一些当前实验水平难以达到的情况。然而在众多的模拟方法中,第一性原理计算凭借其独特的精度和无需经验参数而得到众多研究人员的青睐,成为计算材料学的重要基础和核心计算。本章将介绍第一性原理计算的理论基础,研究方法和ABINIT 软件包。 1.1第一性原理 第一性原理计算(简称从头计算,the abinitio calculation),指从所要研究的材料的原子组分出发,运用量子力学及其它物理规律,通过自洽计算来确定指定材料的几何结构、电子结构、热力学性质和光学性质等材料物性的方法。基本思想是将多原子构成的实际体系理解成为只有电子和原子核组成的多粒子系统,运用量子力学等最基本的物理原理最大限度的对问题进行”非经验”处理。【1】第一性原理计算就只需要用到五个最基本的物理常量即(b o k c h e m ....)和元素周期表中各组分元素的电子结构,就可以合理地预测材料的许多物理性质。用第一性原理计算的晶胞大小和实验值相比误差只有几个百分点,其他性质也和实验结果比较吻合,体现了该理论的正确性。

第一性原理计算按照如下三个基本假设把问题简化: 1.利用Born-Oppenheimer 绝热近似把包含原子核和电子的多粒子问题转化为多电子问题。 2.利用密度泛函理论的单电子近似把多电子薛定谔方程简化为比较容易求解的单电子方程。 3.利用自洽迭代法求解单电子方程得到系统基态和其他性质。 以下我将简单介绍这些第一性原理计算的理论基础和实现方法:绝热近似、密度泛函理论、局域密度近似(LDA)和广义梯度近似(GGA)、平面波及赝势方法、密度泛函的微扰理论、热力学计算方法和第一性原理计算程序包ABINIT 。 1.2量子力学与Born-Oppenheimer 近似 固体是由原子核和核外的电子组成的,在原子核与电子之间,电子与电子之间,原子核与原子核之间都存在着相互作用。从物理学的角度来看,固体是一个多体的量子力学体系【2】,相应的体系哈密顿量可以写成如下形式: ),(),(R r E R r H H ψψ= (1-1) 其中r,R 分别代表所有电子坐标的集合、所有原子核坐标的集合。在不计外场作用下,体系的哈密顿量日包括体系所有粒子(原子核和电子)的动能和粒子之间的相互作用能,即 N e N e H H H H -++= (1-2) 其中,以是电子部分的哈密顿量,形式为:

第一性原理简介

1什么是第一性原理? 根据原子核和电子互相作用的原理及其基本运动规律,运用量子力学原理,从具体要求出发,经过一些近似处理后直接求解薛定谔方程的算法,称为第一性原理。广义的第一原理包括两大类,以Hartree-Fock自洽场计算为基础的从头算和密度泛函理论(DFT)计算。 从定义可以看出第一性原理涉及到量子力学、薛定谔方程、Hartree-Fock自洽场、密度泛函理论等许多对我来说很陌生的物理化学定义。因此我通过向师兄请教和上网查资料一点点的了解并学习这些知识。 2第一性原理的作用 以密度泛函理论(DFT)为基础以及在此基础上发展起来的简单而具有一定精度的局域密度近似(LDA)和广义梯度近似(GGA)的第一性原理电子结构计算方法,与传统的解析方法一样,不但能够给出描述体系微观电子特性的物理量如波函数、态密度、费米面、电子间互作用势等,以及在此基础上所得到的体现体系宏观物理特性的参量如结合能、电离能、比热、电导、光电子谱、穆斯堡尔谱等等,而且它还可以帮助人们预言许多新的物理现象和物理规律。密度泛函计算的一些

结果能够与实验直接进行比较,一些应用程序的发展乃至商业软件的发布,导致了基于密度泛函理论的第一原理计算方法的广泛应用。 密度泛函理论(DFT)为第一性原理中的一类,在物理系、化学、材料科学以及其他工程领域中,密度泛函理论(DFT)及其计算已经快速发展成为材料建模模拟的一种“标准工具”。 密度泛函理论可以计算预测固体的晶体结构、晶格参数、能带结构、态密度(DOS)、光学性能、磁性能以及原子集合的总能等等。 3第一性原理怎么用? 目前我所学到的利用第一性原理的软件为Material Studio、V ASP软件。其中Materials Studio(简称MS)是专门为材料科学领域研究者开发的一款可运行在PC上的模拟软件。使化学及材料科学的研究者们能更方便地建立三维结构模型,并对各种晶体、无定型以及高分子材料的性质及相关过程进行深入的研究。模拟的内容包括了催化剂、聚合物、固体及表面、晶体与衍射、化学反应等材料和化学研究领域的主要课题。 模块简介 Materials Studio采用了大家非常熟悉的Microsoft标准用户界面,允许用户通过各种控制面板直接对计算参数和计算结果进行设置和分析。目前,Materials Studio软件包括如下功能模块: Materials Visualizer: 提供了搭建分子、晶体及高分子材料结构模型所需要的所有工具,可以操作、观察及分析结构模型,处理图表、表格或文本等形式的数据,并提供软件的基本环境和分析工具以及支持Materials Studio的其他产品。是Materials Studio产品系列的核心模块。 Discover: Materials Studio的分子力学计算引擎。使用多种分子力学和动力学方法,以仔细推导的力场作为基础,可准确地计算出最低能量构型、分子体系的结构和动力学轨迹等。

第一性计算原理

Vasp 我所用第一原理是基于密度泛函(DFT)的从头计算,是以电子密度作为基本变量(HK定理),通过求解kohn-sham方程,迭代自洽得到体系的基态电子密度,然后求体系的基态性质。还有一种是基于hartree-fock自洽计算,通过自洽求解HF方程,获得体系的波函数,求基态性质。KS方程的计算水平达到了HF水平,同时还考虑了电子间的交换关联作用。关于DFT中密度泛函的Function其实是交换关联泛函,包括LDA,GGA,杂化泛函等等。一般LDA为局域密度近似,在空间某点用均匀电子气密度作为交换关联泛函的唯一变量,多数为参数化的CA-PZ方案;GGA为广义梯度近似,不仅将电子密度作为交换关联泛函的变量,也考虑了密度的梯度为变量,包括PBE,PE.RPBE等方案。 在处理计算体系中原子的电子态时有两种方法,一种是考虑所有电子叫做全电子法,比如WIEN2K中的FLAPW方法(线性缀加平面波);另一种是只考虑价电子而把芯电子和原子核构成离子实放在一起考虑即赝势法,一般贋势法是选取一个截断半径,截断半径以内波函数变化较平滑,和真实的不同,截断半径以外则和真实情况相同,而且贋势法得到的本征值和全电子法应该相同。贋势的测试标准应是贋势与全电子法计算结果的匹配度,而不是贋势与实验结果的匹配度,因为和实验结果的匹配可能是偶然的。 关于Ecut的收敛测试。一般情况下,总能相对于不同Ecut做计算,当截断能增大时总能变化不明显即可。但是在需要考虑体系应力时,还需要对应力进行收敛测试,而且应力相对于截断能要比总能更为苛刻。也就是某个截断能下总能已经收敛了,但应力未必收敛。(力的计算是在能量的基础上进行的,能量对坐标的一阶导数得到力。计算量的增大和误差的传递导致力收敛慢。) K点也是需要经过测试的。 何时需要考虑自旋?例如BaTiO3中,三个元素分别为=+2,+4,-2价,离子全部为各个轨道满壳层的结构,此时就不必考虑自旋了。对于BaMnO3中,由于Mn+4价时d轨道还有电子但未满,因此需要考虑Mn(4s23d5)的自旋,Ba和O就不必考虑。其实设定自旋就是给定一个原子磁矩的初始值,只在刚开始计算时作为初始值使用。 几何优化包括晶格常数和原子位置的优化,一般情况下也有不优化几何结构直接计算电子结构的,但是对于缺陷形成的计算则往往要优化。 软件大致分为基于平面波的软件,如CASTEP,PWSCF.ABINIT等,计算量大概和体系原子数目的三次方相关;还有基于原子轨道线性组合的软件,比如openmx等,计算量和体系原子数目相关,一般可模拟较多原子数目的体系。 V ASP是使用贋势和平面波基组,进行从头量子力学分子动力学计算的软件包。V ASP中的方法基于有限温度下的局域密度近似(用自由能作为变量)以及对每一MD步骤用有效矩阵对角方案和有效Pulay混合求解瞬时电子基态。这些技术可以避免元氏的Car-Parrinello 方法存在的一切问题,而后者是基于电子、离子运动方程同时积分的方法。离子和电子的相互作用超缓Vanderbilt贋势(US-PP)或投影扩充波(PAW)方法描述。两种技术都可以相当程度地减少过度金属或第一行元素的每个原子所必须的平面波数量。V ASP可以很容易地计算力与张力,用于把原子衰减到其瞬时基态中。!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! V ASP程序亮点: 1、使用PAW方法或超软贋势,因此基组尺寸非常小,描述材料一般需要原子不超过100 个平面波,大多数情况下甚至每原子50个平面波就能得到可靠结果。 2、2. 在平面波程序中,某些部分代码的执行是三次标度。在VASP中,三次标度部分的前 因子足可忽略,导致关于体系尺寸的高效标度。因此可以在实空间求解势的非局域贡献,

第三节第一性原理计算简介

第一性原理计算简介 在物理学中,第一性原理计算或称从头计算是指,基于构建物理学的基础定理,不作任何假设,例如:经验模型和拟合参数,所进行的计算研究。特别地,在凝聚态物理中,指的是运用薛定愕方程在一定的近似情况下,但不包括拟合实验数据所得到的参数和模型,对物质的电子结构进行计算r 从而得到所研究物质的性质的一种研究方法。近些年,随着计算机技术的飞速发展,其运算能力越来越强大,使得人们可以处理更庞大更繁杂的物质结构体系,同时也使得计算物理成为了现代物理学,尤其是在凝聚态物理领域的一个重要分支。众所周知,固体是由相对重且带正电的粒子——原子核,以及相对轻且带负电的粒子——电子聚集在一起构成的。如果有个原子,需要处理的问题是包含有N+ZN(Z 为原子核所含的质子的个数)个粒子的电磁相互作用,是一个多体问题。另一方面,由于处理的是微观粒子的运动,所以需要运用量子力学来描述其基本的运动规律和相互作用。对于该系统,精确的多粒子哈密顿量可以写作: i 2i i i 1122R H M ?=--∑∑ Fuuuuuuuuj 其中位于為处的原子核的质量为M,.,位于巧处的电子的质量为m 一第一项是原子核的动能算符,第二项是电子的动能算符。后三项分别是描述电子与原子核,单个电子与其它电子以及单个原子核与其它原子核之间的库伦相互作用。很显然,直接精确求解(1.64)式几乎是不可能的。为了在合理的近似条件下得到体系的本征值,需要作不同层次的近似。 1.3.1波恩-奥本海默(Bom-Oppenheimer)近似 由于原子核的质量远大于电子质量,所以,原子核的运动速度远小于电子。因此,可以将原子“冻结”在固定的位置,并假设电子在瞬时与原子核是平衡的。或者说,只有电子在这个多体问题中是考察对象,原子核仅仅被当作一个带正电的外源场,相对于电子云是外在独立的。该近似被称为波恩-奥本海默(Bom-Oppenheimer)近似。原来的多体问题被简化成在原子的静电势下,瓜个带负电的粒子的相互作用。波恩-奥本海默认为,原子核不再运动,其动能为零,因此,(1.64)式的第一项被消除,最后一项退化为常数。(1.64)式简化为只含有电子气的动能,电 子与电子之间的相互作用所产生的势能,以及电子在可看作外源的原子核的势中的势能。(1.64)式可重写为: H = f + V + V^, 值得注意的是,(1.65)式中的动能以及电子与电子间的相互作用只取决于所处理的是系统是多电子系统,而不是多质子系统中强的原子内部作用力,并不依赖于特定的多电子系统本身,例如,Br2或者水分子,Cu 还是Fe, bcc-Fe 还是fcc-Fe,等等。因此,前两项是普适的,包含特定系统信息的部分均在第三项中。 1.3.2 密度泛函理论(Density functional theory) 在波恩-奥本海默近似后,该量子多体问题得到了极大的简化,但是,依然很难直接求解。存在许多方法将方程(1.65)进一步近似变为易于处理的形式,历史上非常重要的是Hartree-Fock 方法。该方法在处理原子以及分子时效果很好,因此在量子化学中被广泛使用。但对于处理固体问题,其精度不够高。本文中使用的是更为现代且可能更强大的方法:密度泛函理论。 密度泛函理论的建立可以追溯到1964年Hohenberg 和Kohn[7]提出的两条定理。 1.3. 2.1 Hohenberg-Kohn 定理 两条定理的原始表述如下: 第一定理:多电子体系(原子,分子,固体)基态时的电荷密度pOO 与外源的势之间存在着一一

第一性原理计算及相关理论方法

第二章第一性原理计算及相关理论方法 固体能带理论是凝聚态物理中最成功的理论之一,是固体电子论的支柱。原则上固体能带理论能够阐明和解释固体的许多基本物理性质,如力学,电学,光学及磁学等性质。固体能带理论的主要任务是确定固体电子能级,也就是能带。而要确定固体电子能带,其出发点便是求解组成固体的多粒子系统的定态薛定愕方程。但是对固体这样有1029/m3数量级个原子核和电子的复杂的多粒子系统而言,其薛定愕方程是无法精确求解的,必须作近似简化求解:通过绝热近似将核运动和电子运动分开;通过Hartree-Fock 自洽场方法或更严格、更精确的密度泛函理论(DFT) 将多电子问题简化为单电子问题;通过将固体抽象为平移对称性的理想晶体,将多体问题归结为单电子在周期性势场中的运动。而上述的这些近似简化最终把求解薛定愕方程的全部复杂性都归入了所谓的交换关联泛函,可见交换关联泛函在密度泛函理论中占有重要地位。目前密度泛函理论已经成为探索具有重要应用背景的功能材料和结构材料的重要理论方法。基于密度泛函理论,根据基函数选取的不同有多种具体的计算方法,通常都称为第一性原理计算(ab initio calculation)。所谓第一性原理,即从最基本的物理规律出发,求解体系的薛定愕方程以获取材料性能的信息,从而理解材料中出现的一些现象,预测材料的性能。 第一性原理计算方法[79, 80]有着半经验(HF)方法不可比拟的优势,因为它只需要知道构成微观体系各元素的原子序数,而不需要任何可调(经验或拟合)参数,就可以基于量子力学来处理体系中的电子运动,来计算出该微观体系电子的波函数和对应的本征能量,从而求得系统的总能量、电子结构以及成键、弹性、稳定性等性质。它们被广泛应用于原子,分子,固体,固体表面,界面,超晶格材料,低维材料的电子结构和物理性质的计算,并取得了惊人的成功。随着最近几十年计算机技术的飞速发展,第一性原理的规模和效率都有了极大的提高。目前计算物理已经从解析理论物理和实验物理中完全独立了出来,成为物理学中不可缺少的一个独立分支– 计算物理学[81]。 能带理论的周期性表象为布洛赫(Bloch)表象,同时还存在着局域化的瓦尼(Wannier)表象,Marzari N.等人[82, 83]发展了一套构造最局域Wannier 函数的数值方法,使得Wannier 轨道与第一性原理很好的结合,在传统的化学键分析和最近比较热的拓扑体系的非平庸的表面态计算上都有很大的应用。 传统的半导体模型分析k · p 方法[84],对于分析有效质量和光学性质有很大的帮助,同时它也一直是很好的能带近似手段,能准确的抓住物理实质,便于解析理论分析。 本章的主要任务是对密度泛函理论[85, 86],以及相应的数值的Wannier 方法 [82, 83]和解析的k · p 能带[84]分析方法作简要介绍。 2.1 密度泛函理论(DFT) 对于简单的晶体模型,描述多粒子体系的多体薛定谔方程可以写成如下形式: (),(,)H r R E r R ψ=ψ (2.1)

第一性原理简介精选文档

第一性原理简介精选文 档 TTMS system office room 【TTMS16H-TTMS2A-TTMS8Q8-

1什么是第一性原理? 根据原子核和电子互相作用的原理及其基本运动规律,运用,从具体要求出发,经过一些近似处理后直接求解的算法,称为第一性原理 。广义的第一原理包括两大类,以Hartree-Fock 自洽场计算为基础的从头算和(DFT )计算。 从定义可以看出第一性原理涉及到量子力学、、 Hartree-Fock 自洽场、等许多对我来说很陌生的物理化学定义。因此我通过向师兄请教和上网查资料一点点的了解并学习这些知识。

2第一性原理的作用 以密度泛函理论(DFT)为基础以及在此基础上发展起来的简单而具有一定精度的局域密度近似(LDA)和广义梯度近似(GGA)的第一性原理电子结构计算方法,与传统的解析方法一样,不但能够给出描述体系微观电子特性的物理量如波函数、态密度、费米面、电子间互作用势等,以及在此基础上所得到的体现体系宏观物理特性的参量如结合能、电离能、比热、电导、光电子谱、穆斯堡尔谱等等,而且它还可以帮助人们预言许多新的物理现象和物理规律。密度泛函计算的一些结果能够与实验直接进行比较,一些应用程序的发展乃至商业软件的发布,导致了基于密度泛函理论的第一原理计算方法的广泛应用。 密度泛函理论(DFT)为第一性原理中的一类,在物理系、化学、材料科学以及其他工程领域中,密度泛函理论(DFT)及其计算已经快速发展成为材料建模模拟的一种“标准工具”。 密度泛函理论可以计算预测固体的晶体结构、晶格参数、能带结构、态密度(DOS)、光学性能、磁性能以及原子集合的总能等等。 3第一性原理怎么用? 目前我所学到的利用第一性原理的软件为Material Studio、VASP软件。其中Materials Studio(简称MS)是专门为材料科学领域研究者开发的一款可运行在PC上的模拟软件。使化学及材料科学的研究者们能更方便地建立三维结构模型,并对各种晶体、无定型以及高分子材料的性质及相关过程进行深入的研究。模拟的内容包括了催化剂、聚合物、固体及表面、晶体与衍射、化学反应等材料和化学研究领域的主要课题。 模块简介

相关主题
相关文档 最新文档