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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】截断能优化该选择曲线的最低点还是继续优化以致稳定？
252/3返回列表上一页123下一页
查看: 2004  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭