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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】截断能优化该选择曲线的最低点还是继续优化以致稳定？
252/3返回列表上一页123下一页
查看: 2360  |  回复: 24
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

ylli

木虫 (小有名气)

能量和体积间,可以有这种抛物线

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
mazuju028(金币+1):谢谢交流 2010-11-07 23:25:34
应该是你的计算不够准确. 另外,你注意到能量变化其实已经很小. 具体不是很清楚你这种原因. 没碰到过. 你固定一个合适的k点,然后再做下看看. 一定记住, 其他参数保持一致,不要用不同的参数去作.
赞一下(2人) 回复此楼
11楼2010-11-07 21:15:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

awmc2008

至尊木虫 (文坛精英)

学习者
奇怪
回复此楼
xuexizhe
12楼2010-11-07 21:27:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)
引用回帖:
Originally posted by ylli at 2010-11-07 21:15:03:
应该是你的计算不够准确. 另外,你注意到能量变化其实已经很小. 具体不是很清楚你这种原因. 没碰到过. 你固定一个合适的k点,然后再做下看看. 一定记住, 其他参数保持一致,不要用不同的参数去作.

这是我用4 4 4的Kmesh做到encut优化
引用回帖:
280 -489.207017
290 -489.279444
300 -489.313468
310 -489.317510
320 -489.309982
330 -489.285661
340 -489.261826
350 -489.234883

这是我的encut脚本，自认应该不会有什么问题的
引用回帖:
s2 afm1-dopSe:3.125%
PREC = Accurate
ISPIN = 2
MAGMOM = 16*0 16*3 16*-3 32*0
ISMEAR = -5
ISTART = 0
ICHARG = 2
ENCUT = $i
EDIFF = 1E-5
LORBIT = 11
LREAL = Auto
NPAR ＝ 1
LPLANE ＝ .TRUE.
LSCALU ＝ .FALSE.
NSIM ＝ 4
!
echo "ENCUT = $i eV"
/home/software/openmpi-1.2.2-intel9/bin/mpirun -np 4 /home/bin/vasp.openmpi
E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'`
echo $i $E >>comment
done
赞一下 回复此楼
明察、慎思、笃行
13楼2010-11-07 21:40:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ylli

木虫 (小有名气)

截断能含义

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2):谢谢指教 2010-11-08 14:48:47
gleerat(金币+2):谢谢解答，赠送7个金币 2010-11-08 21:32:51
说下截断能含义. 截断能决定着展开电子波函数的平面波的个数. 在平面波(PW)法求解电子波函数中, 平面波应该构成正交归一完备系, 所以精确的讲,需要无穷多的PW来展开电子波函数(对平面波的波矢求和). 然而,数值计算中,只能使用有限多的PW,这就要求在某个能量(平面波的波矢)进行截断. 截断能越大, 使用的平面波越多, 计算得到的电子总能也越接近于电子的基态. 所以,随着截断能升高, 电子的总能应该降低才合理. 你这个出现振荡, 是不合理的.
赞一下(2人) 回复此楼
14楼2010-11-08 12:57:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ym23

铁杆木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by ylli at 2010-11-08 12:57:49:
说下截断能含义. 截断能决定着展开电子波函数的平面波的个数. 在平面波(PW)法求解电子波函数中, 平面波应该构成正交归一完备系, 所以精确的讲,需要无穷多的PW来展开电子波函数(对平面波的波矢求和). 然而,数值计算 ...

那既然能量应该一直下降，我们到底应该取那个能量值作为计算的截断能呢？
赞一下(1人) 回复此楼
博观约取、厚积薄发---壁立千仞、无欲则刚
15楼2010-11-08 16:26:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)
★ ★
bingmou(金币+2):很详细，谢谢 2010-11-08 18:20:10
引用回帖:
Originally posted by ylli at 2010-11-08 12:57:49:
说下截断能含义. 截断能决定着展开电子波函数的平面波的个数. 在平面波(PW)法求解电子波函数中, 平面波应该构成正交归一完备系, 所以精确的讲,需要无穷多的PW来展开电子波函数(对平面波的波矢求和). 然而,数值计算 ...

我纠正一下，POTCAR中的默认值是300eV。
也就是说，从300eV开始总能是随着波数的增加而增加的。那么，会不会存在由于VASP本身的原因（比如说算法、赝势等）造成CaFe2As2（Se摻杂）体系的计算结果出现曲线的走向呢？因为我还从来没有见到类似的结果（我也找不出我错的地方）。
赞一下(1人) 回复此楼
明察、慎思、笃行
16楼2010-11-08 17:28:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ylli

木虫 (小有名气)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
sunyang1988(金币+1):谢谢交流 2010-11-08 19:31:06
gleerat(金币+1): 2010-11-13 13:06:39
引用回帖:
Originally posted by ym23 at 2010-11-08 16:26:44:


那既然能量应该一直下降，我们到底应该取那个能量值作为计算的截断能呢？

呵呵,这就是进行截断能测试的目的. 当两次能量差别很小(meV量级), 这时候的截断能基本就可以了. 第一性原理的计算只能准确到meV.
赞一下(2人) 回复此楼
17楼2010-11-08 18:08:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ylli

木虫 (小有名气)
★
bingmou(金币+1):感谢参与 2010-11-08 18:20:40
gleerat(金币+5): 2010-11-08 21:30:23
引用回帖:
Originally posted by gleerat at 2010-11-08 17:28:54:

我纠正一下，POTCAR中的默认值是300eV。
也就是说，从300eV开始总能是随着波数的增加而增加的。那么，会不会存在由于VASP本身的原因（比如说算法、赝势等）造成CaFe2As2（Se摻杂）体系的计算结果出现曲线的走向 ...

对于Fe基超导这类体系很难收敛, 我听说.
赞一下(1人) 回复此楼
18楼2010-11-08 18:10:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

limin2007

铁杆木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by ylli at 2010-11-08 18:08:45:

呵呵,这就是进行截断能测试的目的. 当两次能量差别很小(meV量级), 这时候的截断能基本就可以了. 第一性原理的计算只能准确到meV.

呵呵，谢谢，正好解决了我今天的问题了。
赞一下(1人) 回复此楼
19楼2010-11-08 19:22:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)
引用回帖:
Originally posted by ylli at 2010-11-08 18:10:00:

对于Fe基超导这类体系很难收敛, 我听说.

那么，如何看是否收敛呢？
还有，在无法收敛时，一般都是怎么做的呢？
赞一下 回复此楼
明察、慎思、笃行
20楼2010-11-13 17:10:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gleerat 的主题更新
252/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 听说今天filecode变了 +20 布布和一二 2026-08-06 40/2000 2026-08-07 09:49 by olivermiaoer
[基金申请] filecode +12 等待解的谜 2026-08-06 16/800 2026-08-07 09:01 by wlwhappy
[基金申请] 基金中了 +9 laoda193707 2026-08-06 9/450 2026-08-07 08:54 by ZJTJZ
[基金申请] 固定端突然变了，今天 +4 archvillain 2026-08-06 8/400 2026-08-07 08:06 by 带你看大海
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 6/300 2026-08-06 20:19 by 苏知砚
[有机交流] 一个有机合成实验室都需要哪些设备？ 50+3 kf2781974 2026-07-31 12/600 2026-08-06 15:11 by eddyin
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[教师之家] 咨询面上基金 +4 李长云 2026-07-31 7/350 2026-08-06 11:20 by 李长云
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的，举个手。玩一下，释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！ +22 低垂的野花 2026-07-31 30/1500 2026-08-05 18:03 by 低垂的野花
[基金申请] 好消息？这个有何含义？？？ +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的？有收到消息的吗？ +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[论文投稿] 十年后又回来了，论文投稿求助 +3 哈哈114477 2026-08-01 3/150 2026-08-04 15:40 by tegsgjy20
[基金申请] 纯娱乐，不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息，有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[基金申请] 什么时候能放榜呀？ +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
关闭
请填处理意见
关闭