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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】encut收敛测试精确度不够
北京石油化工学院2026年研究生招生接收调剂公告
61/1返回列表
查看: 616  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 283求调剂 +5 baiiyu 2026-04-05 6/300 2026-04-05 20:35 by 啵啵啵0119
[考研] 材料与化工306分找调剂 +12 沧海轻舟e 2026-04-03 13/650 2026-04-04 23:45 by lqwchd
[考研] 一志愿上海大学生物学346 +3 上海大学346调剂 2026-04-03 3/150 2026-04-04 20:20 by dongzh2009
[考研] 求生物学调剂 +14 15172915737 2026-04-01 14/700 2026-04-04 20:13 by babysonlkd
[考研] 309分085801求调剂 +11 MY_angel 2026-03-31 11/550 2026-04-04 19:11 by 蓝云思雨
[考研] 求调剂：085600材料与化工，考材科基，总分319 +21 678lucky 2026-03-31 26/1300 2026-04-04 16:22 by dongzh2009
[考研] 本科211，专业085404，293分请求调剂 +5 莲菜就是藕吧 2026-04-04 5/250 2026-04-04 14:08 by 这是一个无聊的�
[考研] 材料专硕322分 +13 哈哈哈吼吼吼哈 2026-04-01 13/650 2026-04-03 16:08 by 哦哦123
[考研] 一志愿北京交通大学材料工程总分358 +4 cs0106 2026-04-03 4/200 2026-04-03 13:41 by 百灵童888
[考研] 交通运输考试264分求工科调剂 +4 jike777 2026-04-02 4/200 2026-04-02 21:53 by zllcz
[考研] 296求调剂 +4 sdhu 2026-04-02 4/200 2026-04-02 21:29 by baoball
[考研] 348求调剂 +11 zzzzyk123 2026-04-01 11/550 2026-04-02 16:52 by Wang200018
[考研] 282求调剂 +13 呼吸都是减肥 2026-04-01 13/650 2026-04-02 14:10 by baoball
[考研] 材料工程322分 +8 哈哈哈吼吼吼哈 2026-04-01 8/400 2026-04-02 11:53 by 3041
[考研] 一志愿北交大材料工程总分358 +3 cs0106 2026-04-02 5/250 2026-04-02 11:37 by olim
[考研] 272求调剂，接受跨专业调剂！ +4 闲鱼卢 2026-03-31 4/200 2026-04-02 11:18 by guyan1000
[考研] 354求调剂 +4 lxb598 2026-03-31 5/250 2026-04-02 09:55 by Jaylen.
[考研] 292求调剂 +17 木虫er12138 2026-04-01 17/850 2026-04-01 21:37 by 七度不信任
[考研] 379求调剂 +3 ?苦瓜不苦 2026-04-01 3/150 2026-04-01 20:09 by 暮云清寒
[考研] 考研材料工程351分调剂 +5 整个好的 2026-03-31 5/250 2026-04-01 09:36 by topgun2009
信息提示
关闭
请填处理意见
关闭