| 查看: 3053 | 回复: 7 | |||||
abdoman木虫 (著名写手)
|
[求助]
VASP计算流程讨论 已有1人参与
|
|
本人新手，学习VASP，主要是利用小木虫来学习，发现资料很多，良莠不齐，确实很伤神。 在此，写下自己的经验，与教训。只有大致流程，具体可以搜索论坛。 【注意】 vasp计算很多关键词，建议一定要以官方说明书为参考。而且是对应版本的参考教程。vasp不断在更新完善，所以很有必要将所学教程与所用计算的版本统一。 比如 计算用的势能文件是否支持。 我是采用的VASP5.12说明问题, 本人化学专业，非物理，这个不足，不准确的地方还望批评指正！ VASP 所需文件 INCAR 总调度，各项参数的设定 POTCAR 势能的选择。 目测现在PBE是最流行的，当然也会涉及到是否含GW，是否含P层，s层电子（这个我不清楚，也不懂，不深究，最简单方法就是查询相关体系，别人用啥，我用啥） KPOINTS 感觉就是周期性的设置，你考虑多大范围，当然越多，越大，就越好 POSCAR 这个应该就是分子坐标（也许是晶胞）文件 2. 在计算的时候，一般都是先对参数优化，在做精确计算和分析 在参数优化的时候，一般来说，有： A。 晶体常数 B。INCAR 中的ISMEAR，SIGMA, ENCUT, K 点数目 A 的优化 这个有教程， http://muchong.com/bbs/viewthread.php?tid=3234667 说明一点，我在优化是设置了 LWAVE = .FALSE.和 LCHARG = .FALSE. 因为没有用。 其他参数的选择，因为后续会再优化，所以不要太离谱就行。 I. ISMEAR 的选择，官网有建议 For semiconductors or insulators use the tetrahedron method (ISMEAR=-5） Avoid to use ISMEAR >0 for semiconductors and insulators http://cms.mpi.univie.ac.at/vasp ... _SMEARINGS_tag.html ENCUT 建议看看POTCAR 你们ENMAX 最大是多少，15%-20% 扩大一点就行，太小，不准，太大，太贵！ K 点的话 看你的体系大小，7 7 7 我的最爱，至于G还是M，涉及到布里渊中心位置，我不懂，貌似奇数比较靠谱。偶数要平移一下。 至于脚本，批处理，这个是简单linux 知识。 完成后，就对能量做图，是个双曲线，找到最低。就知其CONSCAR 为新的POSCAR。 B 的优化，最主要就是有ENCUT， SIGMA， KPOINT， 类似A中的优化，能量变化小到（1meV）你认为可以了就行。但不是所有都能达到这个标准。论坛有很多帖子介绍，可以看看。 3. band 计算 A。 计算前用最优的参数，进行SCF计算。 B。 然后拷贝CHGCAR 到新文件夹中，设置ISTART=1，NBANDS （增大，曾多少，我一般2-3倍，因为光线性质计算也需要这么多的）。 ICHARG=11 ， 再就是设计KPONITS 文件，走一个流程，主要是布里渊区的高对称点。 http://www.cryst.ehu.es/cgi-bin/cryst/programs/nph-table 这个里面有对应的点群。 进行计算，数据处理，搜索论坛，我也没搞明白 4. dos 计算 A。 同样的最优参数的SCF计算。也可以直接拷贝3中的CHGCAR。 B。 dos计算，增大k点，原因：http://cms.mpi.univie.ac.at/vasp ... e_calculations.html 设置ISTART=1， ICHARG=11 5. 进行光学性质计算 http://cms.mpi.univie.ac.at/wiki ... c_properties_of_SiC 需要3 中的WAVECAR。 增大的NBANDS 因为： Reasonable results are usually only obtained, if the parameter NBANDS is roughly doubled or tripled in the INCAR file with respect to the VASP default. http://cms.mpi.univie.ac.at/vasp ... lectric_matrix.html 数据处理，还不会。 最主要就是流程。 |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP | 第一性原理计算-tg | vasp |
» 猜你喜欢
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有5人回复
-
评审感受-评审感受-评审感受
已经有20人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于vasp计算共吸附
已经有14人回复
- 关于VASP的一些问题？
已经有3人回复
- VASP计算中total-force为什么等于零
已经有7人回复
- 请问vasp里面有Kubo-Greenwood公式的计算么
已经有7人回复
- VASP能带计算问题
已经有10人回复
- VASP计算neb的步骤
已经有14人回复
- vasp计算中，如何选择正确的K点
已经有7人回复
- vasp用neb计算迁移势垒
已经有11人回复
- vasp计算静态出现错误！
已经有5人回复
- vasp加电场的计算
已经有14人回复
- vasp计算结果处理态密度要用到的程序
已经有18人回复
- VASP中计算dos的相关设置
已经有6人回复
- 【活动】vasp计算分子动力学讨论合集及专贴（筹划中）
已经有33人回复
- 【求助】VASP加电场的具体计算步骤是怎么样的？
已经有31人回复
- 【原创】vasp单机双核计算流程
已经有8人回复
- 【其他】vasp neb 计算讨论
已经有20人回复
- 【讨论】用vasp计算过弹性常数的请进！！更新中！！
已经有52人回复
wujunlin1501
铜虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 111.5
- 帖子: 5
- 在线: 35小时
- 虫号: 4121151
- 注册: 2015-10-06
- 性别: GG
- 专业: 无机非金属类光电信息与功
2楼2016-03-29 14:35:58
lvxinyuabc
木虫 (正式写手)
材料大数据
- 应助: 3 (幼儿园)
- 金币: 3083.9
- 散金: 5
- 红花: 2
- 帖子: 520
- 在线: 211.4小时
- 虫号: 1043016
- 注册: 2010-06-17
- 专业: 金属基复合材料
3楼2016-03-30 01:22:03
guoxu3822650
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 164.5
- 帖子: 86
- 在线: 26.3小时
- 虫号: 3877483
- 注册: 2015-05-18
- 专业: 凝聚态物性 II ：电子结构
4楼2016-06-06 15:37:56
guoxu3822650
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 164.5
- 帖子: 86
- 在线: 26.3小时
- 虫号: 3877483
- 注册: 2015-05-18
- 专业: 凝聚态物性 II ：电子结构
5楼2016-06-06 20:40:00
*#？！
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3695.6
- 散金: 170
- 红花: 2
- 沙发: 1
- 帖子: 676
- 在线: 120小时
- 虫号: 4214821
- 注册: 2015-11-12
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
6楼2016-09-16 17:52:47
*#？！
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3695.6
- 散金: 170
- 红花: 2
- 沙发: 1
- 帖子: 676
- 在线: 120小时
- 虫号: 4214821
- 注册: 2015-11-12
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
7楼2016-09-16 18:02:32
8楼2016-09-17 22:41:39
30