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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】截断能优化该选择曲线的最低点还是继续优化以致稳定？
52/1返回列表
查看: 2013  |  回复: 24
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

gleerat

木虫 (正式写手)

[交流] 【求助】截断能优化该选择曲线的最低点还是继续优化以致稳定？ 已有11人参与

我做个了截断能的优化
引用回帖:
280 -489.204150
290 -489.276800
300 -489.310773
310 -489.314845
320 -489.307293
330 -489.283012
340 -489.259010
350 -489.232337

在这里是应该选择最低点的310eV还是应该继续优化？
回复此楼

» 猜你喜欢
明察、慎思、笃行
1楼 2010-11-05 20:58:40
已阅   回复此楼   关注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的回帖
查看全部 25 个回答

lzl8181

铁杆木虫 (著名写手)

单身浏览汉
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1):谢谢交流 2010-11-06 08:25:49
貌似收敛时再增加截断能总能应该不变，最后应该是一条水平线，K点测试也该如此。
赞一下(2人) 回复此楼
版权乜有，盗用不究！
3楼2010-11-06 05:59:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
cenwanglai(金币+1):谢谢回复！ 2010-11-06 20:22:28
引用回帖:
Originally posted by gleerat at 2010-11-05 20:58:40:
我做个了截断能的优化
在这里是应该选择最低点的310eV还是应该继续优化？

你每次改变的能量太小了，试一试每次幅度为50eV。这个明显没有收敛
赞一下(2人) 回复此楼
心之所向无惧无悔
4楼2010-11-06 09:24:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)
cenwanglai:啥机子配置？ 2010-11-06 20:23:10
引用回帖:
Originally posted by sunyang1988 at 2010-11-06 09:24:57:

你每次改变的能量太小了，试一试每次幅度为50eV。这个明显没有收敛

我的结构是80个原子的，使用了5×5×5的K点，当encut的值取到370eV时已经算不动了，提示错误
引用回帖:
ENCUT = 370 eV
running on    3 nodes
distr:  one band on    1 nodes,    3 groups
vasp.4.6.28 25Jul05 complex
POSCAR found :  4 types and   80 ions
LDA part: xc-table for Pade appr. of Perdew
POSCAR, INCAR and KPOINTS ok, starting setup
forrtl: error (78): process killed (SIGTERM)
3 processes killed (possibly by Open MPI)
赞一下(1人) 回复此楼
明察、慎思、笃行
5楼2010-11-06 14:05:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 25 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭