| 查看: 6789 | 回复: 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 | 第一性原理 |
» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有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
- 金币: 3681.7
- 散金: 5802
- 红花: 79
- 沙发: 35
- 帖子: 4560
- 在线: 916小时
- 虫号: 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
- 金币: 3681.7
- 散金: 5802
- 红花: 79
- 沙发: 35
- 帖子: 4560
- 在线: 916小时
- 虫号: 966216
- 注册: 2010-03-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 第一性原理
7楼2011-10-13 09:07:22
8楼2011-10-16 16:29:54
