24小时热门版块排行榜    

查看: 3256  |  回复: 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

木虫 (正式写手)

引用回帖:
11楼: Originally posted by hs0218 at 2016-08-19 23:12:19
依照我的经验 你测试车模型ENCUT=450 K=15就足够了
我不知道你所谓的0.001eV的误差接受度是哪里来的,而且不同原子能量大小不一样,同一个原子所建立的不同大小的胞能量也不一样,那么0.001eV/atom 或者0.001eV/ce ...

感谢应助!我所说的0.001eV的误差接受度是从网上流传的一些资料里看来的,看来是我没理解准确,现在明白了。
17楼2016-08-21 11:25:16
已阅   回复此楼   关注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

木虫 (正式写手)

引用回帖:
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的回帖
信息提示
请填处理意见