VASP四大输入文件——KPOINTS
本文转载于微信公众号,VASP学习交流,将会持续更新
今天这里讲的是常规计算使用的五行式的KPOINTS,一般计算要在第一布里渊区均匀撒点,能带计算在高对称点连线路径上取值,能带计算的以后再说。
五行式KPOINTS
K-POINTS
0
Gamma
1 1 10 0 0
解释
第一行随便写都行,不识别,但是一定不能什么都不写
第二行是数字0,表示格子自动生成,很多人会误写成氧符号o,
第三行G表示的是以gamma点为中心生成网格,还有一种是原始的Monkhorst-Pack网格,又可以叫M点,二者的区别是M在G的基础上在三个方向平移了1/(2N)个单位,所以G点是M点的特殊情况;VASP只认第一个字母,大小写均可,所以Gamma可以写为G,g,gamma,ga....,Monhorst可以写为M, m等。。。。
第四行表示在xyz三个方向上生成对应数目的K点,本例是1*1*1,意思是在倒格矢a,b,c方向上都只取一个k点,一共1个k点,4*4*1意思是在倒格矢a,b方向上都取4个k点,c方向上取一个,一共16个,K点的密度由KPOINTS决定,KPOINTS取点越多,包含到计算里的信息越多,但是经常有人拿不准K点数目怎么取,经常出错,所以在接下来会总结几个方法教大家怎么取K点
第五行表示shift的值,一般都写成 0 0 0,一般不用改,Gamma center的K点就相当于MP方法shift了0.5 0.5 0.5,
KPOINTS的选取
方法1,使用官网建议的取值,每个晶格矢量的长度(a, b, c)乘以这个方向上的K点数目(k)。需要做收敛性测试,比如要算一个绝缘体体系, a,b的晶胞长度为6Å,c方向真空层15Å,K点就是3 3 1
方法2,把结构拖到MS里,点击CASTEP Tools图标,查看不同K点下,K点密度,一般密度在0.03附近最佳,对于正交晶系,晶格常数取模,然后取倒数。这三个倒数的比值就是三个方向上k-points密度的取值,例如,fcc Cu的KPOINTS则是7 7 7,
方法3,对于非正交晶系,倒格矢长度和实空间晶格常数不满足反比关系,所以用以上方法就不行了,这个时候用vaspkit可以产生所有体系的K点,vaspkit--102--1(M点),2(G点)--选择倒格子中k点间距,单位是Angstrom-1,一般计算使用0.04,精确计算0.03或0.02
方法4,简单粗暴,一般晶格常数乘以K点数目大约在30~40之间即可,K点数目足够多的情况下,前面的问题都不是问题,但是不要过分大,导致撒点过于密集,浪费机时,过分小撒点过疏,计算不准;
不管做什么,正式计算前一定要做收敛性测试。
KPOINTS注意事项
真空层一定不要加K点,用1即可,因为这个方向都没有结构,还加K点一来没意义,二来纯粹搞笑,比如说二维体系在xy方向是周期性体系,z方向存在真空层,那么这个方向一定不要加K点,用1即可;
对于原子或者分子的计算,K点取一个gamma点就够了(1 1 1),也就是这个KPOINTS文件可以用于绝大多数的原子或者分子计算,不用再修改;但如果你要用ISMEAR = -5 来计算,那么就需要把 1 1 1 改成 2 2 2或者3 3 3
对于原子或者分子的计算,K点取一个gamma点(1 1 1)就够了,因为这些原子和分子不是周期性的,多的K点是能提高周期性镜像分子间的相互作用精度,但是这部分能量是我们不想要的。即:对于含有真空层的体系,在真空层的方向上永远只使用一个K点(真空层一定不要加K点,用1即可)。多余的K点只会增加真空层两边体系的相互作用的精度,而这一部分是我们不想要的。
对于六角晶系hexagonal的结构,官网建议用gamma centered,因为M平移之后,网格的对称性和晶胞的对称性会出现不匹配的情况,从而导致计算出错;
对于特别大的体系,也可以使用单Gamma点,只含有一个Gamma点的计算可以使用vasp_gam版加速计算
一般需要做收敛性测试,测试,测试
自动生成KPOINTS的脚本
五行式每次输入挺麻烦的,直接写一个小脚本,每次一键操作,直接调用,把脚本命名为kpoints.sh(或者别的名字)
echo K-POINTS >>KPOINTS
echo 0 >>KPOINTS
echo G >>KPOINTS
echo $1 $1 $1>>KPOINTS
echo 0 0 0>>KPOINTS
chmod u+x kpoints.sh 赋予脚本可执行权限
VASP四大输入文件——KPOINTS相关推荐
- vasp测试计算机,科学网—PWSCF 自洽计算、kpoints测试和ecut测试 - 叶小球的博文
关注: 1) 自洽计算的目的是为了check结构建立是否合理,程序能否顺利运行: 2) kpoints测试的重要性: 选择合适的参数,可在保证质量的情况下,节省计算时间 3) ecut 测试的重要 ...
- Vasp二维材料单胞基本计算范例
Vasp二维材料单胞基本计算范例,以InSe为例子. 总的来说,对于非磁性的简单二维材料体系,计算分为以下几步 1. 结构优化 2.静态计算 3.能带和态密度的计算 上面几个过程,实际需要改动的参数只 ...
- vasp 模拟退火_VASP 计算问题小结
本文转载自博主一个人就是一个叠加态,有部分删减修改,文中对相关概念方法做了详细的总结,留坑待填... 1. 第一原理计算的一些心得 1. 第一性原理 第一性原理其实是包括基于密度泛函的从头算和基于Ha ...
- VASP计算笔记_弹性常数计算
二阶弹性常数计算 主要采用能量–应变方法 AELAS软件计算 必要的文件 elas_strain.sh elas_cal.sh POTCAR INPOS 优化后的POSCAR.可在脚本设置或者单独给出 ...
- vasp能带结构计算
Si (band structure) VASP中能带结构计算流程与态密度相同,具体不同点如下: 1. scf自洽计算 与态密度dos计算中的scf自洽计算完全相同. 2. 能带结构非自洽计算 能带计 ...
- Vasp进行分子动力学模拟关键词解析及计算示例1
针对周期性体系的分子动力学模拟(包含计算所需输入文件) 通过vasp进行分子动力学的计算需要进行以下具体步骤: 1 选择一个足够大的晶胞制作成POSCAR,原子数越多越好(100个以上),一般k点较小 ...
- vasp可以做哪些计算及命令
VASP是目前材料微观反应机理和计算材料电子结构性质科学研究中最流行的一款软件,它可以处理金属及其氧化物.半导体.晶体.掺杂体系.纳米材料.分子.团簇.表面体系和界面体系等. 结构性质:晶格常数.原子 ...
- 用VASP求硅的电子态密度和能带
本文介绍用VASP如何求硅的电子态密度和能带,共分为如下5个部分: (1). 生成4个输入文件: POSCAR POTCAR INCAR KPOINTS (2). 优化晶格参数,求出能 ...
- vasp结构优化英语_VASP态密度、能带、DFT+U实例分析教程
今天主要记录强关联电子体系+U.DOS和能带计算过程中参数的设置,相信大家在读一些计算类文献时会看到文章中计算部分有The GGA + U calculations were performed us ...
- VASP(1)_参数测试_截断能(ENCUT)测试
必要输入文件: encut.sh POTCAR KPOINTS和POSCAR可以在encut.sh里设置,也可以独立给出. encut.sh脚本参考 #!/bin/bash#生成KPOINTS文件,K ...
最新文章
- python下载方法-python实现下载文件的三种方法_python
- 【Linux】5.linux下的export命令和环境变量
- 区间比较_Simulink(离散PIDamp;区间检测模块)+AURIX功能安全板过流比较电路+电机结构应用与工作原理学习网站...
- [POI2008] Poc (原名 Trians) Treap+Hash
- datagrid 重载本地数据_音视频系列3:使用ffmpeg + nginx搭建本地转发服务器
- PageObjects 设计模式
- android静态类保存context,Android:静态获取Context是一个好习惯吗?
- laravel项目白屏问题解决办法及原因
- cpython和jpython_在J内部使用Python
- 《社交网站界面设计(原书第2版)》——3.2 注册
- FFmpeg — 屏幕录制器材
- leetcode python3 简单题118. Pascal's Triangle
- HR人事管理系统源码
- google保存网页为图片
- 在本地存储中存入和取出对象
- 金庸感人的经典爱情对白
- 【oiClass 2085】马(排序,DP)
- 双核Linux内核的引导
- JAVA初学(七):银行账户演示程序【改进版】
- Ubuntu14.04和16.04官方默认更新源sources.list和第三方源推荐(干货!)