24小时热门版块排行榜    

查看: 1551  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 一志愿北京化工大学 070300 学硕 336分 求调剂 +5 vv迷 2026-03-22 5/250 2026-03-23 07:36 by Iveryant
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
[考研] 求调剂 +5 Zhangbod 2026-03-21 7/350 2026-03-22 13:13 by Zhangbod
[考研] 285求调剂 +6 ytter 2026-03-22 6/300 2026-03-22 12:09 by 星空星月
[考研] 280求调剂 +11 咕噜晓晓 2026-03-18 12/600 2026-03-21 22:40 by ACS Nano——
[考研] 一志愿南大,0703化学,分数336,求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 材料学学硕080502 337求调剂-一志愿华中科技大学 +4 顺顺顺mr 2026-03-18 5/250 2026-03-21 10:22 by luoyongfeng
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 316求调剂 +6 梁茜雯 2026-03-19 6/300 2026-03-21 06:32 by Ecowxq666!
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 材料与化工 322求调剂 +4 然11 2026-03-19 4/200 2026-03-20 22:12 by luoyongfeng
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 350求调剂 +5 weudhdk 2026-03-19 5/250 2026-03-20 22:04 by luoyongfeng
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
信息提示
请填处理意见