24小时热门版块排行榜

返回列表

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

vigaryang

木虫 (正式写手)

应助: 16 (小学生)
金币: 3466.3
散金: 677
红花: 6
帖子: 471
在线: 495.7小时
虫号: 227097
注册: 2006-03-23
性别: GG
专业: 应用化学

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

利用侯博士教程里提供的脚本，对Ni bulk进行了k值的测试（静态计算），结果如下图：

根据该教程及其它很多资料里的说法，两个相邻k值对应的总能变化在0.001eV左右才OK。那么我都已经把k值取到19了，不论encut取值是多少，相邻k值对应的总能变化仍然大于0.001eV。不知道是哪里出了问题？

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

求解惑，谢谢！

回复此楼

» 猜你喜欢

Cas 72-43-5需要30g，定制合成，能接单的留言已经有8人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有6人回复
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？已经有8人回复
磺酰氟产物，毕不了业了！已经有5人回复
论文终于录用啦！满足毕业条件了已经有25人回复
2026年机械制造与材料应用国际会议（ICMMMA 2026）已经有3人回复
自荐读博已经有3人回复
不自信的我已经有5人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

1楼 2016-08-17 16:59:59

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

从心而觅

新虫 (小有名气)

应助: 34 (小学生)
金币: 115.6
红花: 17
帖子: 206
在线: 133小时
虫号: 4330822
注册: 2016-01-04
性别: GG
专业: 理论和计算化学

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
vigaryang: 金币+5, ★★★很有帮助, 谢谢！ 2016-08-21 11:30:27

ISMEAR决定了如何确定每个波函数的占有数，驰豫计算时一般不需要这个参数。
进行任何的静态计算或态密度计算，且K点数目（从IBZKPT文件中读取）大于4时，取ISMEAR＝－5；当由于原胞较大而K点数目较少（小于4个）时，取ISMEAR＝0，并设置一个合适的SIGMA值。另外对半导体或绝缘体的计算（不论是静态还是结构优化），取ISMEAR＝－5;当体系呈现金属性时，取ISMEAR＝1和2，以及设置一个合适的SIGMA值。在进行能带结构计算时，ISMEAR 和SIGMA用默认值就好。一般说来，无论是对何种体系，进行何种性质的计算，采用ISMEAR＝0 ，并选择一个合适的SIGMA值都能得到合理的结果。
另外收敛标准是0.001ev/atom，需要除以原子数。