首页时事民生政务教育文化科技财富体娱健康情感

旅行百科职场楼市企业乐活学术汽车时尚创业美食幽默美体文摘

VASP干货 | VASP表面计算步骤

2024-10-25 18:07 广东

一、概述

vasp用“slab” 模型来模拟表面体系结构。

vasp计算表面的大概步骤是：

材料体性质的计算；表面模型的构造；表面结构的优化；表面性质的计算。

二、分步介绍

1、材料体性质计算：

本步是为了确定表面计算时所需的一些重要参数：ENCUT、SIGMA、晶格参数。

在计算前，要明确：何种PP；ENCUT；KPOINTS ;SIGMA;PREC;EX-CO，这其实是准备proper input files。

a. 何种PP

选择的PP能使计算得到的单个原子能量值在1meV~10meV之间,

所求得的单原子能量（对称性破缺时）可用来提高结合能的精度。

b.ENCUT

选择的ENCUT应使得总能变化在0.001eV左右为宜。

注意：试探值最小为POTCAR中的ENMAX（多个时，取最大的），递增间隔50；

另外，在进行变体积的结构优化时，最好保证ENCUT＝1.3ENMAX，以得到合理精度。

c.PREC

控制计算精度的最重要参数，决定了（未指定时）ENCUT、FFT网格、ROPT取值。

一般计算取NORMAL；当要提高Stress tensor计算精度时，HIGH 或ACCURATE，并手动设置ENCUT。

d.EDIFF＆EDIFFG

EDIFF 判断电子结构部分自恰迭代时自恰与否，一般取默认值＝1E－4；

EDIFFG 控制离子部分驰豫

e.ISTART＆ICHARGE

ISTART = 1, ICHARG = 11：能带结构、电子态密度计算时；

ISTART =0, ICHARG = 2：其余计算

ISTART = 1，ICHARG = 1（其他所有不改变）：断点后续算设置

f.GGA&VOSKOWN

GGA＝91： Perdew -Wang 91；

GGA＝PE：Perdew-Burke-Ernzerhof

VOSKOWN＝1（ GGA＝91时）；VOSKOWN＝默认（其余情况）

g.ISIF

控制结构参数之优化。在对原胞进行变形状或者体积的优化时，ENCUT要取大（比如1.3ENMAX或PREC＝HIGH）,以消除Pulay Stress导致的误差。

h.ISMEAR&SIGMA

进行任何静态计算时，且K点数目大于4，ISMEAR＝－5；

当原胞太大，导致K点数目小于4时，ISMEAR = 0，并且要设置一个SIGMA；

对绝缘体和半导体，不论是静态计算还是结构优化，ISMEAR = -5；

对金属体系，SMEAR=1和 2，并且设置一个SIGMA

能带结构计算，用默认值：ISMEAR＝1，SIGMA＝0.2；

一般来说，对于任何体系，任何计算，采用ISMEAR＝0，并选择合适的SIGMA都会得到合理结果。

选择的SIGMA应使得entropy T*S EENTRO 绝对值最小。K 点数目变化后，SIGMA需再优化。

i.RWIGS

一般取POTCAR中以A为单问的RWIGS值

j.K points

选择的K点应使得总能变化在0.001eV左右即可。

k. 一些重要的参数在默认下的值NSW =0，IBRION＝－1，ISIF＝2：静态计算。

2：

a.体材料结合能修正

在OUTCAR中energy without entropy之后的那个能量值，就是修正值

b. 结构参数优化。简单情况（没有内部自由度如晶胞形状、原子位置）：静态计算，得出E～V关系，然后用Birch-Murnaghan状态方程拟合。

复杂情况：总思路是先“建立好房子”，再“放好桌子”。

先算一步结构优化（取ISIF＝5，只改变“房子”形貌，房间大小不变，家具不予考虑），接着算一步静态自恰计算，从而得到某结构参数下的能量，如此循环得到E～V关系。用状态方程拟合得到平衡体积。

在该体积下，重复1（取ISIF＝2,房子造好后，考虑的是如何放家具。此处一般是使得每个家具受力达到某中小即可认为达到稳定结构）、2两步，便得到了所有的晶格参数值，如离子坐标。

c. VASP得到的总能即是结合能，不过还要减去前面得到的修正值。

d.自恰的电荷密

优化得到晶格参数后，再进行静态的自恰计算，就得到了自恰的电荷密度。

此时的POSCAR为从优化晶格参数时可CONTCAR得到。

KPOINTS 不变

典型的INCAR设置是：

ENCUT = 250

ISTART = 0; ICHARG = 2

ISMEAR = -5

PREC = Accurate

计算完后，注意保存相关结果。

本文转自孙祥科学网博客。

链接地址：https://blog.sciencenet.cn/blog-3474450-1293286.html

点击阅读原文，提交计算需求！

http://mp.weixin.qq.com/s?__biz=MzkzOTIwMzI0Mg==&mid=2247532644&idx=1&sn=94130e201c0d98a43e02ee0a2b47b84b

DFT计算参谋长

华算科技朱参谋长公众号，十余年VASP重度使用经验，每天给大家带来期刊速览，学术新闻，欢迎关注！

最新文章

学习VASP一定要知道的七点！！K点、赝势、INCAR参数...

1天倒计时！理论计算全场8折，速来享华算科技年度钜惠！

急急急！最后2天！吹爆这个理论计算全场8折福利！快冲！！！

INCAR、POTCAR 、POSCAR、KPOINTS核心输入文件解析！

年终狂欢倒计时3天！理论计算8折起，全场钜惠不容错过，速来薅羊毛！

VASP简介及计算前期准备工作

预存高至30%增值，科研双11超级重磅加码，全年无忧，狂发顶刊！

单原子催化HER/OER/CO2RR/NO3RR反应计算

做计算不会选泛函？来看看这篇总结！

惊！华算科技计算全场8折，2000元+返利送送送，仅限本周！

如何用VASP做第一性原理分子动力学计算？

惊！华算科技计算全场8折，2000元+返利送送送，仅限本周！

惊！华算科技计算全场8折，2000元+返利送送送，仅限本周！

抢疯了！华算科技超级重磅加码，科研双11享高至30%增值经费预存，错过等一年！

预存高至30%增值，科研双11超级重磅加码，全年无忧，狂发顶刊！

抢疯了！华算科技超级重磅加码，科研双11享高至30%增值经费预存，错过等一年！

如何利用理论计算预测平衡电位

预存高至30%增值，科研双11超级重磅加码，全年无忧，狂发顶刊！

预存高至30%增值，科研双11超级重磅加码，全年无忧，狂发顶刊！

抢疯了！华算科技超级重磅加码，科研双11享高至30%增值经费预存，错过等一年！

VASP考虑范德瓦尔斯计算结构要点总结

抢爆了！经费预存享高至17%增值，更有 ¥ 8500+返利！科研人都惊了……

导师要求学点DFT计算，怎么入门更轻松？

居然免费算！快来薅华算科技的羊毛！！

炸裂！Nature 正刊又+1，审稿人惊呼：这计算，绝了！

抢爆了！经费预存享高至17%增值，更有 ¥ 8500+返利！科研人都惊了……

VASP计算光学性质INCAR

绝了！Nature正刊又双叒叕+1，审稿人连夸计算太给力了！

VASP赝势文件如何合并？

VASP非金属催化计算培训：HER、OER/ORR、CO2RR、NO3RR、NRR、d带中心、差分电荷密度、态密度、吸附能

VASP零基础入门培训：晶体与二维材料的结构、弹性、电子、光学、磁性、吸附、催化性质

VASP干货 | VASP表面计算步骤

【今晚19点】直播揭秘：VASP vs QE，全面展示过渡态计算流程！

【免费公开课】全流程实操解析，一堂课理清 QE与VASP 过渡态计算差异！

计算干货 | 半导体中载流子迁移率计算

VASP非金属催化计算培训：HER、OER/ORR、CO2RR、NO3RR、NRR、d带中心、差分电荷密度、态密度、吸附能

VASP钙钛矿/半导体计算专题培训：二维结构/异质结/掺杂/电子/吸附/缺陷性质

VASP锂电池计算专题培训：Li离子电池/燃料电池/电位/容量/离子扩散/析氢反应/氧还原反应/自由能/过电位

VASP计算声子谱

VASP零基础入门培训：晶体与二维材料的结构、弹性、电子、光学、磁性、吸附、催化性质

学习计算材料学，你一定要记住这10个常用的Linux命令!

今晚七点！免费分享！【分子动力学特训营】GROMACS：建模+MD+分析

明晚开讲！【分子动力学特训营】GROMACS：建模+MD+分析

VASP二维材料计算培训：二维/异质结构建模、能带、功函数、静电势、吸附性质、虚频消除、态密度、d带中心、CO2RR、电池性质

VASP单原子催化计算培训：单/双原子模型，电子性质，HER、OER/ORR、CO2RR、NO3RR

今晚七点开讲！从零开始追赶诺奖，如何在化学研究中使用AI？

明晚开讲！从零开始追赶诺奖，如何在化学研究中使用AI？

VASP钙钛矿/半导体计算专题培训：二维结构/异质结/掺杂/电子/吸附/缺陷性质

QE学习笔记2：输入参数CONTROL各项含义

VASP燃料电池计算培训：热门催化剂建模、d带中心、功函数、催化剂对HOR与ORR反应中间体的吸附、反应自由能、反应过渡态

分类

原创标签

时事社会财经军事教育体育科技汽车科学房产搞笑综艺明星音乐动漫游戏时尚健康旅游美食生活摄影宠物职场育儿情感小说曲艺文化历史三农文学娱乐电影视频图片新闻宗教电视剧纪录片广告创意壁纸头像心灵鸡汤星座命理教育培训艺术文化金融财经健康医疗美妆时尚餐饮美食母婴育儿社会新闻工业农业时事政治星座占卜幽默笑话独立短篇连载作品文化历史科技互联网

发布位置

广东北京山东江苏河南浙江山西福建河北上海四川陕西湖南安徽湖北内蒙古江西云南广西甘肃辽宁黑龙江贵州新疆重庆吉林天津海南青海宁夏西藏香港澳门台湾美国加拿大澳大利亚日本新加坡英国西班牙新西兰韩国泰国法国德国意大利缅甸菲律宾马来西亚越南荷兰柬埔寨俄罗斯巴西智利卢森堡芬兰瑞典比利时瑞士土耳其斐济挪威朝鲜尼日利亚阿根廷匈牙利爱尔兰印度老挝葡萄牙乌克兰印度尼西亚哈萨克斯坦塔吉克斯坦希腊南非蒙古奥地利肯尼亚加纳丹麦津巴布韦埃及坦桑尼亚捷克阿联酋安哥拉