24小时热门版块排行榜    

查看: 7352  |  回复: 15
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

xianglannihao

铜虫 (小有名气)

[交流] 关于VASP收敛判断

关于VASP收敛判断:
假定设置如下:
EDIFF    = 1E-05      (SCF convergence criteria)
EDIFFG   =    -0.03 (Ionic convergence criteria, eV/AA)

我的理解收敛条件:EDIFF 是电子步收敛的条件,电子步dE达到1E-5就认为收敛 。
EDIFFG是离子步收敛的标准,为负数时表示力的收敛判据。是否收敛看OUTCAR中的 reached required accuracy - stopping structural energy minimisation提示和最后一轮的如下两列:
d Force = 0.1423097E-03[ 0.108E-03, 0.177E-03]  d Energy = 0.1531149E-03-0.108E-04
d Force =-0.6844216E+00[-0.684E+00,-0.684E+00]  d Ewald  =-0.6844216E+00-0.177E-07


问题:1. 是否有这个提示就表示收敛了?
2. OUTCAR中为什么有两列d Force, 都表示什么意思?d Force d Energy表示的都是一个范围?
3. 有时在每一轮的离子步计算时,其电子步dE没达到也显示该轮收敛。
4. 优化好的构型,以相同的设置再重新优化一次,我发现其能量会发生变化,有时会大于1eV,但构型变化倒不是很大?难道我们要一直优化,直至能量不变,是否这样才能得到可靠的计算结果?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wmy8802217

木虫 (正式写手)

★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
xianglannihao: 金币+2 2015-10-07 10:30:44
引用回帖:
5楼: Originally posted by xianglannihao at 2015-10-05 16:21:40
只是很好奇,为什么优化后的构型,再次以相同设置优化按道理应该不变化才对吧?你试过拿着优化好的构型再优化吗?...

再优化跑了几个离子步?
凝聚态理论与计算物理
6楼2015-10-05 19:20:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答

wmy8802217

木虫 (正式写手)

★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
xianglannihao: 金币+3 2015-10-05 12:19:58
EDIFFG应该说的是每个原子的受力收敛吧,应该看OUTCAR里的Total-Force;
EDIFF我觉得应该是D eps和dE的收敛;
至于第四个问题,你应该加大NSW,增大离子步,因为你之前的结果应该是没收敛的。
凝聚态理论与计算物理
2楼2015-10-05 12:02:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xianglannihao

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by wmy8802217 at 2015-10-05 12:02:52
EDIFFG应该说的是每个原子的受力收敛吧,应该看OUTCAR里的Total-Force;
EDIFF我觉得应该是D eps和dE的收敛;
至于第四个问题,你应该加大NSW,增大离子步,因为你之前的结果应该是没收敛的。

谢谢回复。
1. OUTCAR里的每个原子在三个文献上Total-Force我看了,除了设置了固定原子位置的原子,其它原子的力都小于0.03,已达到力收敛标准。
2. EDIFF只表示其中一个吧?
3. NSW      =      500,已经足够大了吧,之前的结果都提示已达收敛。除了设置了固定原子位置的原子,其它原子的力都小于0.03,已达到力收敛标准。
3楼2015-10-05 12:19:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wmy8802217

木虫 (正式写手)

★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
xianglannihao: 金币+3 2015-10-05 16:19:06
引用回帖:
3楼: Originally posted by xianglannihao at 2015-10-05 12:19:47
谢谢回复。
1. OUTCAR里的每个原子在三个文献上Total-Force我看了,除了设置了固定原子位置的原子,其它原子的力都小于0.03,已达到力收敛标准。
2. EDIFF只表示其中一个吧?
3. NSW      =      500,已经足够大 ...

我觉得你可以把标准再调低一点看看,既然已经到了收敛标准,没必要再来一次;(也许对于你这个力收敛精度,对于总能的影响可能就在1eV左右)
比较能量你可以用更大的K点和截断能计算更精确的总能来比较。
凝聚态理论与计算物理
4楼2015-10-05 13:11:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +11 吃吃吃才有意义 2026-03-19 12/600 2026-03-24 09:18 by 吃吃吃才有意义
[考研] 材料292调剂 +4 橘颂思美人 2026-03-23 4/200 2026-03-23 23:16 by peike
[考研] 070300化学求调剂 +8 苑豆豆 2026-03-20 8/400 2026-03-23 20:57 by baobaoye
[考研] 328求调剂,英语六级551,有科研经历 +7 生物工程调剂 2026-03-17 12/600 2026-03-23 18:18 by YMU施老师
[考研] 求调剂一志愿武汉理工大学材料工程(085601) +3 WW.' 2026-03-23 5/250 2026-03-23 17:18 by 枫翼ljj
[考研] 350求调剂 +6 weudhdk 2026-03-19 6/300 2026-03-23 15:47 by tangyuan0840221
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 材料与化工085600,总分304,本科有两篇sci参与,求调剂 +4 幸运的酱酱 2026-03-22 5/250 2026-03-22 20:15 by edmund7
[考研] 308求调剂 +3 墨墨漠 2026-03-21 3/150 2026-03-22 16:54 by i_cooler
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 265求调剂 +3 Jack?k?y 2026-03-17 3/150 2026-03-21 03:17 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 材料学求调剂 +4 Stella_Yao 2026-03-20 4/200 2026-03-20 20:28 by ms629
[考研] 一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 求调剂 +3 eation27 2026-03-20 3/150 2026-03-20 19:32 by JourneyLucky
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
信息提示
请填处理意见