24小时热门版块排行榜    

查看: 3248  |  回复: 19
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

vigaryang

木虫 (正式写手)

[求助] k值和encut值测试遇到的问题 已有2人参与

利用侯博士教程里提供的脚本,对Ni bulk进行了k值的测试(静态计算),结果如下图:
k值和encut值测试遇到的问题
根据该教程及其它很多资料里的说法,两个相邻k值对应的总能变化在0.001eV左右才OK。那么我都已经把k值取到19了,不论encut取值是多少,相邻k值对应的总能变化仍然大于0.001eV。不知道是哪里出了问题?

对ENCUT的测试也遇到了类似问题,当ENCUT大于等于550eV时,两个相邻ENCUT值对应的总能变化才小于0.001eV。而文献中利用VASP对该体系进行计算时,选取的ENCUT通常在390~450 eV之间,小于我所得到的550eV。也不知道为什么会存在这种差异。

求解惑,谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vigaryang

木虫 (正式写手)

引用回帖:
3楼: Originally posted by mylifecg at 2016-08-17 20:26:56
你确定是静态下测试这些参数的吗?

确定是的。INCAR如下:
SYSTEM = Ni
ENCUT = 550 eV
ISTART = 0
ICHARG = 2
ISMEAR = -5
PREC=Accurate
4楼2016-08-17 21:26:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答

qwangar

金虫 (小有名气)


漫天飘雪: 金币+1, 谢谢交流 2016-08-17 19:10:30
试试选定一个encut 继续增大k值。对小的原胞,k值可以是很大的

发自小木虫IOS客户端

» 本帖已获得的红花(最新10朵)

I know nothing.
2楼2016-08-17 17:31:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mylifecg

木虫 (正式写手)

你确定是静态下测试这些参数的吗?
3楼2016-08-17 20:26:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vigaryang

木虫 (正式写手)

送红花一朵
引用回帖:
2楼: Originally posted by qwangar at 2016-08-17 17:31:53
试试选定一个encut 继续增大k值。对小的原胞,k值可以是很大的

按照您说的试了一下,设定ENCUT=550eV,发现只有当K值大于等于23的时候,相邻K值对应的TOTEN之差才小于0.001eV。我的原胞是含有4个Ni原子的fcc,不知这个K值是否太大了。
5楼2016-08-17 21:28:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见