| 查看: 1855 | 回复: 6 | ||
[求助]
【关于VASP计算流程的一些问题】 已有3人参与
|
|
最近才学vasp，看了很多教程，但是始终对于vasp计算流程还是很模糊，我把我要问的问题总结在下面，希望高手能指点一下。 1. 关于优化流程方面的 看有的说先优化ENCUT，K点，SIGMA，再优化晶包参数；有的说先优化晶包参数，在优化化ENCUT，K点，SIGMA。 在这里我说一下我的疑惑，一般优化ENCUT，K点，SIGMA的标准是能量变化小于1meV吧， 如果晶包参数和理想晶包偏差很大的话，那能量应该也会很大，可能的一种情况是ENCUT取很大的值也不能是能量变化在1meV。 如果是先优化晶包的话，那ENCUT，K点怎么取值？ ENCUT取ENMAX的1.3倍？那K点呢？ 2. 关于K点的选取 本人才接触第一性，没有学过的理论知识，什么布里渊区，K点之类的概念不知道是神马意思，也看过一些理论知识，实在没勇气看下去，所以像问大侠K点是怎么选出来的？ 如果是上面问题1里面说的，优化选出来的，那怎么变化K点？ 是和ENCUT一样慢慢变化K点吗？不过ENCUT只有一个变量，K点有三个变量，如果一个一个改变那计算量不是很大？ 3.计算DOS的时候，要增大K点，怎么个增大法？有没有标准？ 比如我的K点是 6 6 6 增大的话 增大到多少 9 9 9 还是多少？还是要根据DOS的精度要求，DOS图不好的话在把K点继续增高？ 4.最后说一下计算的流程，望高手指点一下看对不对 第一步：就是问题一提到的，优化各参数，晶包 第二步：将晶包优化的CONTCAR cp 成POSCAR 进行一下静态计算 第三步：利用静态计算的CHG,CHGCAR 进行DOS计算 第四步：利用静态计算的CHG,CHGCAR 进行band计算 不知道大概步骤是不是这个样子的，希望高手指点一下，方便像我一样的新手。 |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有12人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有12人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有12人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有9人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有22人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
2026年国自然面上资助率
已经有19人回复
|
本帖内容被屏蔽 |
2楼2014-01-05 21:27:55
谢谢支持！ 3楼2014-01-05 21:29:16
4楼2014-01-05 23:16:49
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
x-cool: 金币+5, ★★★★★最佳答案, 谢谢帮助 2014-01-06 14:57:18
感谢参与，应助指数 +1
x-cool: 金币+5, ★★★★★最佳答案, 谢谢帮助 2014-01-06 14:57:18
|
第一个问题，论坛里有个很好的帖子讲优化的。我只加一句，其实没必要1meV，几个都已经不错了，你可以做测试，如果用比较小的Encut就收敛的挺好了，何必用大的？如果收敛的不是特别好，1.3Enmax足够了。 第二个，我觉得既然要做计算，基本的知识太多了，也没说要全部看完，但我觉得起码要知道怎么一回事，这样万一结算结果很离谱，你也会大概知道哪出了问题。 k点，这个帖子不错：http://muchong.com/bbs/viewthread.php?tid=3180318&fpage=1 第三和第四，步骤都是没问题的，这方面教程也多的很，侯老师那个一般就够用了。问题是过程中会不会有问题。 |
5楼2014-01-05 23:31:46
【答案】应助回帖
★ ★ ★
感谢参与，应助指数 +1
x-cool: 金币+3, ★★★很有帮助, 谢谢帮助 2014-01-06 14:57:31
感谢参与，应助指数 +1
x-cool: 金币+3, ★★★很有帮助, 谢谢帮助 2014-01-06 14:57:31
|
1. 关于优化流程方面的 看有的说先优化ENCUT，K点，SIGMA，再优化晶包参数；有的说先优化晶包参数，在优化化ENCUT，K点，SIGMA。 在这里我说一下我的疑惑，一般优化ENCUT，K点，SIGMA的标准是能量变化小于1meV吧， 如果晶包参数和理想晶包偏差很大的话，那能量应该也会很大，可能的一种情况是ENCUT取很大的值也不能是能量变化在1meV。 如果是先优化晶包的话，那ENCUT，K点怎么取值？ ENCUT取ENMAX的1.3倍？那K点呢？ 问题看的有些晕，如果只回答问好的话，那就1.3*ENMAX吧，K点使劲取 2. 关于K点的选取 本人才接触第一性，没有学过的理论知识，什么布里渊区，K点之类的概念不知道是神马意思，也看过一些理论知识，实在没勇气看下去，所以像问大侠K点是怎么选出来的？ 如果是上面问题1里面说的，优化选出来的，那怎么变化K点？ 是和ENCUT一样慢慢变化K点吗？不过ENCUT只有一个变量，K点有三个变量，如果一个一个改变那计算量不是很大？ 请看教程 3.计算DOS的时候，要增大K点，怎么个增大法？有没有标准？ 比如我的K点是 6 6 6 增大的话 增大到多少 9 9 9 还是多少？还是要根据DOS的精度要求，DOS图不好的话在把K点继续增高？ 是有精度要求的，看你想要个什么样的精度，这个不是随便加大的 4.最后说一下计算的流程，望高手指点一下看对不对 第一步：就是问题一提到的，优化各参数，晶包 第二步：将晶包优化的CONTCAR cp 成POSCAR 进行一下静态计算 第三步：利用静态计算的CHG,CHGCAR 进行DOS计算 第四步：利用静态计算的CHG,CHGCAR 进行band计算 第三第四步我一般将WAVECAR也拷进去进行计算，另外，你要注意，这里面的INCAR和K点文件容易不一样，具体请参见侯老师的教程吧，写的很详细的 |
6楼2014-01-06 13:24:10
guoxu3822650
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 164.5
- 帖子: 86
- 在线: 26.3小时
- 虫号: 3877483
- 注册: 2015-05-18
- 专业: 凝聚态物性 II ：电子结构
7楼2016-07-06 00:55:18
15