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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】截断能优化该选择曲线的最低点还是继续优化以致稳定？
252/3返回列表上一页123下一页
查看: 2402  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然评审意见 +3 wangmingqi 2026-08-28 3/150 2026-08-28 17:43 by James-1988
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +16 爱看书的可乐 2026-08-26 18/900 2026-08-28 17:09 by winsaint
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +11 Kittylucky 2026-08-27 12/600 2026-08-28 16:26 by lyxiaomc
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
关闭
请填处理意见
关闭