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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】encut收敛测试精确度不够
61/1返回列表
查看: 655  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 考博 +4 Eternitf 2026-04-29 4/200 2026-05-02 16:10 by yejenny
[找工作] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 16:06 by zahslv263y
[公派出国] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 15:51 by zahslv263y
[考博] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 15:51 by zahslv263y
[教师之家] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 15:36 by zahslv263y
[论文投稿] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 15:21 by zahslv263y
[公派出国] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 3/150 2026-05-02 14:21 by zahslv263y
[博后之家] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 3/150 2026-05-02 14:06 by zahslv263y
[教师之家] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 3/150 2026-05-02 14:06 by zahslv263y
[博后之家] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 3/150 2026-05-02 10:21 by ovidyohpyz
[找工作] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 4/200 2026-05-02 10:06 by ovidyohpyz
[教师之家] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 3/150 2026-05-02 09:36 by dz7y1p1i98
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 5/250 2026-05-02 09:21 by dz7y1p1i98
[论文投稿] 一篇MDPI论文改变了学习工作和生活 +3 helsinglee 2026-04-29 4/200 2026-05-01 11:45 by bobvan
[有机交流] 一个化合物的合成路线：CAS:367929-02-0 名称：8β-乙烯基雌二醇 50+4 zhquanbing 2026-04-28 4/200 2026-04-30 22:45 by zyqchem
[文学芳草园] 太白金星有点烦 +3 yexuqing 2026-04-27 3/150 2026-04-30 21:56 by myrtle
[考博] 中国地质大学（北京）博士招生补录，数理学院材料科学与工程专业和材料与化工专业 +3 liuxh329 2026-04-29 3/150 2026-04-30 15:42 by 顺利毕业zs
[基金申请] 收到国自然专家邀请后几年才会有本子送过来评 +3 barley5 2026-04-30 3/150 2026-04-30 14:58 by aspect3000
[论文投稿] 有没有快的中文核心比较快录用的，纳米材料光催化 5+3 2914439618 2026-04-27 4/200 2026-04-29 08:59 by 北京莱茵润色
[论文投稿] 有人投过CCC中国控制会议吗？ 19+3 初九凡 2026-04-26 3/150 2026-04-27 08:58 by 北京莱茵润色
信息提示
关闭
请填处理意见
关闭