应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】encut收敛测试精确度不够
61/1返回列表
查看: 564  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

405963512

铜虫 (正式写手)

[交流] 【求助】encut收敛测试精确度不够

大家好，我做encut参数的收敛测试。
#!bin/sh
rm WAVECAR
for a in 400 450 500 550 600 650 700 750
do
cat >INCAR < SYSTEM=napthalene
ENCUT=$a
ISTART=0;ICHARG=2
PREC=Accurate
LREAL= AUTO
!
echo "a=$a"
vasp
E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'`
echo $a $E >>comment
done

这是脚本，计算过程没错，但是结果的精确度貌似不够，看侯老师的教程说，能量差小于0.001eV就可以了，但是我的结果
500 -256.477829
550 -256.507997
600 -256.560073
650 -256.622108
700 -256.671911
750 -256.697650
最小的差值也在0.029eV，所以我想请教大家，这个问题该如何处理，谢谢。
回复此楼

» 猜你喜欢
1楼 2010-02-03 15:14:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hgye

金虫 (小有名气)
★
405963512(金币+1): 2010-02-03 18:07
ice_rain(金币+1):谢谢 2010-02-03 21:09
不知你算的系统多大，我记得是说平均到每个原子的能量差小于1meV就可以了，要给能量差除以原子个数。
赞一下(1人) 回复此楼
2楼2010-02-03 15:48:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

405963512

铜虫 (正式写手)
引用回帖:
Originally posted by hgye at 2010-02-03 15:48:49:
不知你算的系统多大，我记得是说平均到每个原子的能量差小于1meV就可以了，要给能量差除以原子个数。

是这样的吗，貌似不是这么说的吧？？有知道的吗，朋友们？
赞一下 回复此楼
3楼2010-02-03 18:07:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xbelinda

银虫 (小有名气)
★
ice_rain(金币+1):谢谢 2010-02-03 21:13
405963512(金币+1): 2010-02-04 11:20
可以这么说，一般文章会说收敛标准的。 有的是1meV/f.u. 或者1meV/atom。
你可以继续加大ENCUT.可能你算的结构的赝势很hard.
赞一下(1人) 回复此楼
4楼2010-02-03 19:20:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)
★ ★
405963512(金币+1): 2010-02-04 11:20
ice_rain(金币+1):xiexie~ 2010-02-05 08:51
ice_rain(金币+1):xiexie 2010-02-05 08:51
其实收敛到多少主要取决于你要计算的能量差多大。如果你计算的是FM和AFM的能量差（一般在几十meV到几百meV),那么显然你的收敛还不够，因为你的能量差肯定是用两个体系的总能相减而不是每个原子。
而如果是计算晶格常数，这样的能量差基本就足够了。
收敛标准是自己定的，不是书上定的，不要拘泥
赞一下(2人) 回复此楼
5楼2010-02-03 23:20:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liumiao05

至尊木虫 (知名作家)

秦武汉风
★
ice_rain(金币+1):xiexie~ 2010-02-05 08:51
的确有点大，虽然不一定要那么严格，不过我觉得怎么也要在10mev一下吧
最后两个有减小的趋势，再增大试试
赞一下(1人) 回复此楼
天下英雄出我辈，一入江湖岁月催。鸿图霸业谈笑间，不胜人生一场醉。提剑跨骑挥鬼雨，白骨如山鸟惊飞。尘世如潮人如水，只叹江湖几人回。
6楼2010-02-04 00:36:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 405963512 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭