24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

fangyongxinxi

新虫 (正式写手)

应助: 0 (幼儿园)
金币: 41
散金: 69
红花: 3
帖子: 876
在线: 138.5小时
虫号: 437418
注册: 2007-08-26
性别: MM
专业: 凝聚态物性I:结构、力学和

[求助] 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 ]

回复此楼

1楼 2011-10-12 10:44:41

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

liuhuimeiwoo

金虫 (小有名气)

应助: 1 (幼儿园)
金币: 1469.1
散金: 10
红花: 3
帖子: 235
在线: 40.2小时
虫号: 911880
注册: 2009-11-24
性别: MM
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★
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不一定要一样大的

赞一下(1人)

回复此楼

paper我要好多好的的paper

3楼2011-10-12 11:10:45

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 8 个回答

uuv2010

荣誉版主 (职业作家)

1ST强帖: 6
应助: 160 (高中生)
贵宾: 3.14
金币: 3679.7
散金: 5802
红花: 79
沙发: 35
帖子: 4560
在线: 915.2小时
虫号: 966216
注册: 2010-03-09
性别: GG
专业: 凝聚态物性 II ：电子结构
管辖: 第一性原理

【答案】应助回帖

★
fangyongxinxi(金币+5): 2011-10-12 21:21:07
youzhizhe(金币+1): 谢谢交流。 2011-10-13 09:37:27

可以将INCAR贴上来让大家一起分析

赞一下(1人)

回复此楼

2楼2011-10-12 10:55:47

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

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

>为什么在vasp在对结构优化过程中K点设置越大，能量反而越高呢？如果是根据体系能量的大小来判断体系
>的稳定性，看这晶格常数的差别，叫小生情何以堪。
参数的收敛测试可不是说能量逐渐降低啊，而是能量逐渐趋进于恒定值，这个值却不一定是最低的。而选择参数时候要选择收敛于恒定值的那组参数。
就比如你测试k点，假设发现在k>10的时候，得到的能量收敛于1eV，本应选择k=10作为参数，但发现当k=5时的能量值是0eV。那么这时候难道能说：“k=5的结果更稳定，所以选择作为参数”吗？其实这只能说明k=5的时候由于精度不够，使计算得到的能量误差达到了1eV之多，如果反而选择它作为参数，岂不奇怪。