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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】截断能优化该选择曲线的最低点还是继续优化以致稳定？
101/1返回列表
查看: 2018  |  回复: 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

木虫 (正式写手)
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的回帖

gleerat

木虫 (正式写手)
引用回帖:
Originally posted by sunyang1988 at 2010-11-06 14:34:50:
取一个K点做测试吧，原子这么多K点没必要这么密

我做过Kmesh的优化了，3、4、5、6、7都做了，5是合理的。4×4×4的Kmesh比5差一点，前者总能变化在10meV量级，后者在1meV量级。
赞一下 回复此楼
明察、慎思、笃行
8楼2010-11-07 20:55:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)
引用回帖:
Originally posted by ylli at 2010-11-07 21:05:42:
不应该再升高.

那，我这个东西是怎么回事？以前没见过这个V型的encut曲线。
赞一下 回复此楼
明察、慎思、笃行
10楼2010-11-07 21:08:07
已阅   回复此楼   关注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的回帖

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的回帖

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

木虫 (正式写手)
引用回帖:
Originally posted by ustbwgs at 2010-11-14 08:04:13:
你的纵坐标的取值太小了，太小了会使本来趋于稳定的让你看起来有很大的变化，也不知道你的体系有多大，一般纵坐标最小分度有0.1，0.05足够了，这样做出图来，就会使你对它的收敛性有更好的认识，不要老是取作图软 ...

一般来说做encut优化的时候都是以OUTCAR中encut的最大值作为优化取值的最小值，我OUTCAR中为300，所以就从300eV开始算了。
赞一下 回复此楼
明察、慎思、笃行
22楼2010-11-14 17:41:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)
引用回帖:
Originally posted by ustbwgs at 2010-11-14 18:36:55:
你算收敛时不要一开始就取那么小的间隔，可以先取比较大一点的间隔，比如50，当你发现在一个大的间隔内比较稳定时，再确定小范围。
for example:先取300，350，400，450，500，假如当你发现400开始能量收敛时，再 ...

由于这么大的kmesh，使用300eV本来就很大了，所以我就不敢使用大的能量间隔，而事实也证明了当encut取到370时计算无法进行。
赞一下 回复此楼
明察、慎思、笃行
25楼2010-11-16 10:38:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gleerat 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭