| 查看: 6635 | 回复: 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 | 第一性原理 |
» 猜你喜欢
不自信的我
已经有5人回复
-
磺酰氟产物，毕不了业了！
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有16人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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计算过弹性常数的请进！！更新中！！
已经有51人回复
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