- 1、本文档共8页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
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等,计算量和体系原子数目相关,一般可模拟较多原子数目的体系。
VASP是使用贋势和平面波基组,进行从头量子力学分子动力学计算的软件包。VASP中的方法基于有限温度下的局域密度近似(用自由能作为变量)以及对每一MD步骤用有效矩阵对角方案和有效Pulay混合求解瞬时电子基态。这些技术可以避免元氏的Car-Parrinello方法存在的一切问题,而后者是基于电子、离子运动方程同时积分的方法。离子和电子的相互作用超缓Vanderbilt贋势(US-PP)或投影扩充波(PAW)方法描述。两种技术都可以相当程度地减少过度金属或第一行元素的每个原子所必须的平面波数量。VASP可以很容易地计算力与张力,用于把原子衰减到其瞬时基态中。!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
VASP程序亮点:
使用PAW方法或超软贋势,因此基组尺寸非常小,描述材料一般需要原子不超过100个平面波,大多数情况下甚至每原子50个平面波就能得到可靠结果。
2. 在平面波程序中,某些部分代码的执行是三次标度。在VASP中,三次标度部分的前因子足可忽略,导致关于体系尺寸的高效标度。因此可以在实空间求解势的非局域贡献,并使正交化的次数最少。当体系具有大约2000个电子能带时,三次标度部分与其它部分可比,因此VASP可用于直到4000个价电子的体系。VASP使用传统的自洽场循环计算电子基态。这一方案与数值方法组合会实现有效、稳定、快速的Kohn-Sham方程自洽求解方案。程序使用的迭代矩阵对角化方案(RMM-DISS和分块Davidson)可能是目前最快的方案。VASP包含全功能的对称性代码,可以自动确定任意构型的对称性。5. 对称性代码还用于设定Monkhorst-Pack特殊点,可以有效计算体材料和对称的团簇。Brillouin区的积分使用模糊方法或四面体方法。四面体方法可以用Bl?chl校正去掉线性四面体方法的二次误差,实现更快的k点收敛速度。1、VASP
我们平时最常用的研究方法是做单点能计算,结构优化、从头计算的分子动力学和电子结构相关性质的计算。
一般我们的研究可以按照这样的过程来进行
如果要研究一个体系的最优化构型问题可
您可能关注的文档
最近下载
- 2024-2030年中国六氟磷酸锂行业市场现状供需分析及投资评估规划分析研究报告.docx
- 防护工程之三维网植被护坡(附施工图解、工艺、造价编制).pdf VIP
- 《中铁一局集团有限公司工程项目管理绩效考核办法》(中铁一成本[2015]623号).doc
- 机械打眼开炸石方(附施工图解、造价编制).pdf VIP
- 路基土石方(借土填方)(附施工图解、工艺、造价编制).pdf VIP
- 2023扬州龙川控股集团有限责任公司招聘试题及答案解析.docx
- 交通安全设施之交通标志(附施工图解、工艺、造价编制).pdf VIP
- 路基土石方之填挖交界(附施工图解、工艺、造价编制).pdf VIP
- 2纵横造价实战--抛石挤淤附施工图解、工艺、造价编制.pdf
- 2024《我的阿勒泰》作品简介PPT课件(精品).pptx
文档评论(0)