24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1602  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

kelly612

银虫 (小有名气)

[求助] 参数测试能量怎么都不稳定,求救!

最近在做一个体系的参数测试,搞了好久了,但是一直没有得出结果,虽然可以用文献上的参数,但是想自己验证一下,掌握方法。
脚本是这样
#!/bin/sh
for  i  in 350  400  450  500  550  600 650 700
do
cat > INCAR << !
SYSTEM = LMP
ENCUT = $i
ISTART = 0 ; ICHARG = 2
LDAU = .TRUE
LDAUTYPE = 1
LDAUL = -1 2 2 -1 -1
LDAUU = 0 5.5 5.5 0 0
LDAUJ = 0 1 1 0 0
ISPIN = 2
MAGMOM = 4*0 1*7 2*-7 1*7 4*0 16*0
ISMEAR = 0; SIGMA = 0.05
NELM = 30
LMIXMIX =4
LREAL = Auto
PREC = High
!
echo "ENCUT = $i eV" ;  mpirun -np 16 vasp5
E=$(grep "energy without entropy" OUTCAR | tail -1| awk '{printf "%12.6f \n", $5}')
echo $i  $E >>comment
done
KPOINTS 用的 3  6  9
结果是这样:
350        -200.25397
400        -199.72309
450        -199.49069
500        -199.42421
550        -199.44941
600        -199.49795
650        -199.53697
700        -199.57734
K点也取过很大的 趋势也是这样  sigma 也试了不同的值  能量始终都到不了稳定的值  另外我这个POSCAR是用的优化好的结构,我这个体系是反铁磁结构,所以加自旋加U,但我也试过不加的,能量趋势也是这样,然后用文献的值520 来测K点  已经增大到7 15 19了   能量还是一直降低  根本看不到能量变化小与1mev的希望啊……  各位大神出来指点一二吧,我实在没招了……   感激不尽啊!!!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

???????:
6?: Originally posted by kelly612 at 2013-11-13 20:54:05
лл????????  
???????????LiMnPO4, ????28?????????????10.5  6.18  4.8???????????18????????...

?????????????ж??????
?????^????
?]????????

???O????????
1) ?z??Y???Л]????e
    lattice constant and atom position
    ?????Л]????? key in ?e?`

2) POTCAR ?Л]????e
    (ex: element ?????e???????e)
   
3) ?????? VASP ???????????
   (?Еr?? vasp , ?o?e lib ??????????~?????????????e??
    ??????_?J?F??????VASP?????????????)

???_?J????????
??? input files ?G????
?п??????????????
???????F????F??
(script, POSCAR, POTCAR, KPOINTS)

» 本帖已获得的红花(最新10朵)

7楼2013-11-14 00:23:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

啸宇天涯

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
kelly612: 金币+5 2013-11-11 15:05:53
测试时是否收敛了,可否增大离子步,比如NSW=200等等,或者将每个能量的输出结果保存下来试试。
2楼2013-11-10 22:53:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kelly612

银虫 (小有名气)

引用回帖:
2楼: Originally posted by 啸宇天涯 at 2013-11-10 22:53:45
测试时是否收敛了,可否增大离子步,比如NSW=200等等,或者将每个能量的输出结果保存下来试试。

跑每个encut 时  是收敛了的,测试不是静态计算么  就一个离子步啊  没有NSW呢
3楼2013-11-11 15:07:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kelly612

银虫 (小有名气)

咋木有人支招呢    自己顶一个……
4楼2013-11-11 20:03:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 297,工科调剂? +3 河南农业大学-能 2026-04-14 3/150 2026-04-15 22:02 by noqvsozv
[考研] 生物学调剂 +9 纸扇zhishan 2026-04-13 9/450 2026-04-15 18:28 by AN流800
[考研] 211本科材料化工求调剂 +19 YHLAH 2026-04-11 23/1150 2026-04-14 22:25 by fenglj492
[考研] 材料专业344求调剂 +17 hualkop 2026-04-10 22/1100 2026-04-14 16:21 by sxdj2
[考研] 271求调剂 +35 2261744733 2026-04-11 41/2050 2026-04-14 15:36 by zs92450
[考研] 一志愿沪9,326求生物学调剂 +10 刘墨墨 2026-04-13 10/500 2026-04-14 15:16 by zs92450
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +6 zju2000 2026-04-14 17/850 2026-04-14 14:34 by jurkat.1640
[考研] 求调剂 +16 雪逢冬 2026-04-10 16/800 2026-04-14 14:27 by 逆水乘风
[考研] 农学0904 312求调剂 +4 Say Never 2026-04-11 4/200 2026-04-14 09:10 by zs92450
[考研] 085404 298分求调剂 +11 呼啦呼啦呼呼呼 2026-04-10 12/600 2026-04-14 08:38 by wfj257
[考研] 考研英一数一338分 +9 长江大学东校区 2026-04-13 10/500 2026-04-14 00:41 by 王珺璞
[考研] 一志愿中南大学 0855 机械 286 求调剂 +11 不会吃肉 2026-04-12 11/550 2026-04-13 21:59 by bljnqdcc
[教师之家] 山东双非院校考核超级无底线,领导幸灾乐祸,教师遭殃恐 +3 qut2026 2026-04-11 7/350 2026-04-12 20:24 by qut2026
[考研] 2本,初试303,0860求调剂 +6 floriea 2026-04-12 8/400 2026-04-12 18:13 by zhouxiaoyu
[考研] 电气工程专硕320求调剂 +5 小麻子111 2026-04-10 5/250 2026-04-12 10:47 by zhouyuwinner
[考研] 材料工程日语考生求调剂 +7 0856?调剂 2026-04-10 7/350 2026-04-11 21:33 by 蓝云思雨
[考研] 346,工科求调剂 +3 moser233 2026-04-09 3/150 2026-04-11 10:04 by zhq0425
[考研] 281求调剂 +11 觉得好的吧 2026-04-10 11/550 2026-04-11 09:35 by 逆水乘风
[考研] 调剂 +19 小张ZA 2026-04-10 20/1000 2026-04-10 22:08 by 猪会飞
[考研] 083200 初试305分 求调剂 暂不考虑跨专业 +15 Claireyyyy 2026-04-09 15/750 2026-04-09 16:11 by zhuimr
信息提示
请填处理意见