| 查看: 7346 | 回复: 7 | ||||||
fangyongxinxi新虫 (正式写手)
|
[求助]
vasp计算中，如何选择正确的K点
|
|
对于这个问题，自己一直比较郁闷，有两点，第一是如何选取K点，第二是再什么地方设置K点。 1. 自己的情况是这样的：计算一个含有四个原子的金属体系，两次优化结构，K点设置分别为：6*6*6,10*10*10，得到的晶格常数差别比较大0.04埃；利用grep entropy OUTCAR分别查看两次优化后体系的能量，得到6*6*6（-28eV）设置的能量比10*10*10(-27.54eV)能量反而更低。不知道为什么。 像自己在用ATK计算的时候，会对体系做一个K点的循环计算，比如从5*5*5到10*10*10，然后比较每次体系的能量，做一条能量-K点的曲线，在K点大到一定程度的时候，能量的变化会趋于一个很小的值，这样就算是找到了一个合适的K点。同样的方法用于寻找cut-off energy（见附图），这也是为什么常在文章中看到“。。。取更大的cut-off/K点，对体系能量影响小于，。。所以我们的参数是对的。。。，” so,自己的问题是：为什么在vasp在对结构优化过程中K点设置越大，能量反而越高呢？如果是根据体系能量的大小来判断体系的稳定性，看这晶格常数的差别，叫小生情何以堪。 2. 对K点的选取，有的软件有自动生成的，可以有相应的软件来完成，或根据体系手动自己设置（比如自己这个四个原子体系，咨询了一个高手，说K点要设置成27*13*13，我当时看了就懵了，这么大的K点数，计算起来一定非常慢），或者像上面计算一个循环。 so, 问题是，vasp能像上面那样，计算个循环，然后比较能量吗？  补充： 自己用MS计算了上面的体系，则没有出现上面的问题。K点大小的变化（5*5*5 10*10*10）对优化后结果（晶格常数）影响不大：0.002埃。so，自己再一次怀疑了VASP。 求高手解答。 [ Last edited by fangyongxinxi on 2011-10-12 at 10:52 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理 | 仿真建模与计算 | VASP | 第一性原理 |
» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
什么时候开奖？
已经有10人回复
-
跳槽后在研项目怎么办？
已经有9人回复
-
2026国自然函评费到账
已经有16人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
让我中一个面上吧！
已经有14人回复
-
今天放榜吗？
已经有16人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp计算中如何比较掺杂前后结合能
已经有13人回复
- vasp自洽计算K点不匹配问题，急！！
已经有15人回复
- vasp计算DOS,以及如何判断参数的设置是合理的，达到要求的精度
已经有7人回复
- 【求助】vasp计算Partial Charge，能带交叉的如何确定IBAND？
已经有13人回复
- 【求助】晶格常数大了以后 VASP的计算值准确吗
已经有5人回复
- 【求助】Vasp计算能带结构，高对称点如何在origin中表示
已经有17人回复
- 【求助】VASP怎样精确计算出费米能级处的态密度值？
已经有22人回复
- 【求助】第一原理的分子动力学的材料学计算，软件选择VASP还是MS？或者其他
已经有17人回复
- 【求助成功】如何用VASP计算反铁磁
已经有14人回复
- 【讨论】用vasp计算过弹性常数的请进！！更新中！！
已经有52人回复
uuv2010
荣誉版主 (职业作家)
- 1ST强帖: 6
- 应助: 160 (高中生)
- 贵宾: 3.14
- 金币: 3677.7
- 散金: 5807
- 红花: 79
- 沙发: 35
- 帖子: 4564
- 在线: 916.6小时
- 虫号: 966216
- 注册: 2010-03-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 第一性原理
2楼2011-10-12 10:55:47
liuhuimeiwoo
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1469.1
- 散金: 10
- 红花: 3
- 帖子: 235
- 在线: 40.2小时
- 虫号: 911880
- 注册: 2009-11-24
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★
uuv2010(金币+1): 多谢提示 2011-10-12 20:01:36
uuv2010(金币+1): 多谢提示 2011-10-12 20:01:36
| 有些赝势也会对晶格参数造成很大的影响的 这个看你计算的要求是什么了 一般去k点是要慢慢增加k点 然后看能量的变化 一般是先变小后变大的 你去能量最低值的k点应该是不会有什么问题的 如果你觉得这样不可信 你可以参考意见发出来的paper跟你一样或者差不多的体系都是去的什么样的k点 其次 你才四个原子 k点取27*13*13应该也不会很慢 看你自己取得和别人提供的k点 可能你的k点取得有点问题 应该晶格常数差不多是b=c=2a吧 你的k点也要按照比例去的 ka kb kc不一定要一样大的 |
3楼2011-10-12 11:10:45
acridine
木虫 (正式写手)
- 1ST强帖: 11
- 应助: 13 (小学生)
- 金币: 4901.6
- 红花: 13
- 帖子: 706
- 在线: 389.5小时
- 虫号: 305478
- 注册: 2006-12-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★
uuv2010(金币+1): 多谢提示 2011-10-12 20:02:04
fangyongxinxi(金币+10): 2011-10-12 21:20:56
uuv2010(金币+1): 多谢提示 2011-10-12 20:02:04
fangyongxinxi(金币+10): 2011-10-12 21:20:56
|
>为什么在vasp在对结构优化过程中K点设置越大，能量反而越高呢？如果是根据体系能量的大小来判断体系 >的稳定性，看这晶格常数的差别，叫小生情何以堪。 参数的收敛测试可不是说能量逐渐降低啊，而是能量逐渐趋进于恒定值，这个值却不一定是最低的。而选择参数时候要选择收敛于恒定值的那组参数。 就比如你测试k点，假设发现在k>10的时候，得到的能量收敛于1eV，本应选择k=10作为参数，但发现当k=5时的能量值是0eV。那么这时候难道能说：“k=5的结果更稳定，所以选择作为参数”吗？其实这只能说明k=5的时候由于精度不够，使计算得到的能量误差达到了1eV之多，如果反而选择它作为参数，岂不奇怪。  |
4楼2011-10-12 17:52:00
fangyongxinxi
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 41
- 散金: 69
- 红花: 3
- 帖子: 876
- 在线: 138.5小时
- 虫号: 437418
- 注册: 2007-08-26
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
5楼2011-10-12 21:22:43
fangyongxinxi
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 41
- 散金: 69
- 红花: 3
- 帖子: 876
- 在线: 138.5小时
- 虫号: 437418
- 注册: 2007-08-26
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
6楼2011-10-12 21:31:29
uuv2010
荣誉版主 (职业作家)
- 1ST强帖: 6
- 应助: 160 (高中生)
- 贵宾: 3.14
- 金币: 3677.7
- 散金: 5807
- 红花: 79
- 沙发: 35
- 帖子: 4564
- 在线: 916.6小时
- 虫号: 966216
- 注册: 2010-03-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 第一性原理
7楼2011-10-13 09:07:22
8楼2011-10-16 16:29:54
40