| 查看: 1864 | 回复: 10 | ||
[求助]
参数测试能量怎么都不稳定，求救！
|
||
|
最近在做一个体系的参数测试，搞了好久了，但是一直没有得出结果，虽然可以用文献上的参数，但是想自己验证一下，掌握方法。 脚本是这样 #!/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是用的优化好的结构，我这个体系是反铁磁结构，所以加自旋加Ｕ，但我也试过不加的，能量趋势也是这样，然后用文献的值520 来测K点 已经增大到7 15 19了 能量还是一直降低 根本看不到能量变化小与1mev的希望啊…… 各位大神出来指点一二吧，我实在没招了…… 感激不尽啊！！！  |
回复此楼» 猜你喜欢
filecode
已经有13人回复
-
filecode
已经有11人回复
-
听说今天filecode变了
已经有26人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有6人回复
-
固定端突然变了，今天
已经有6人回复
-
基金中了
已经有4人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
-
求各位大神看下
已经有33人回复
-
咨询面上基金
已经有7人回复
-
影响面上的因素
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- NEB的计算出现负值，测试的其他能量都比初始的00和终止的10低
已经有17人回复
- ENCUT收敛测试时，能量差为什么总不能达到1meV呢？
已经有20人回复
- hartree 能量单位和KJ/mol单位换算
已经有10人回复
- 求助：北京低碳清洁能源研究所怎么样？硕士进去待遇如何？
已经有19人回复
- 求助：关于肋片效率问题以及换热系数测量误差分析
已经有14人回复
- 求助，如何使用Forster能量转移理论计算药物小分子与蛋白之间的能量转移啊？
已经有11人回复
- 求助castep计算中：晶胞总能量E 与晶胞体积V 的关系图如何得到，看内容
已经有24人回复
- 【求助】用DMOL3怎么算钛原子能量
已经有11人回复
- 【求助】vasp 算出的0k下的能量如何扩展到其他温度压力
已经有12人回复
- 【求助】反应路径能量图怎么画？
已经有7人回复
- 【求助】VASP里的几个能量怎么区别呢？？？？？？？？
已经有23人回复
- 【求助】请问分子最小能量结果怎么不一致？
已经有12人回复
- 【求助】怎么求N在最稳定单质形态下的能量
已经有5人回复
- 【求助】能量密度是怎么计算的
已经有7人回复
- 【求助】中科院青岛能源所怎么样
已经有32人回复
- 【求助】单原子的总能量怎么算？
已经有37人回复
|
?????????????ж?????? ?????^???? ?]???????? ???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朵）
kelly6127楼2013-11-14 00:23:03
2楼2013-11-10 22:53:45
3楼2013-11-11 15:07:28
4楼2013-11-11 20:03:45
5楼2013-11-13 17:50:09
6楼2013-11-13 20:54:05
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★
kelly612: 金币+8 2013-11-17 18:36:29
kelly612: 金币+8 2013-11-17 18:36:29
| LDAU中的U值应该指定4个就可以了，你指定了5个，大概是你吧自旋上和自旋下的Mn看成两个不同原子了。可以把自旋上和自旋下的Mn当成同一个原子指定U和POTCAR；然后指定U值。U值不知道你根据什么选的5.5-1的组合，可以参照文献指定；截断能也可以参考文献值--也就是你找到文献上用VASP计算的文章，采用和文献值相同的PAW--也就是电子组态要相同，然后直接拿来文献的截断能用即可，无需测试。另外，能量差小于1meV这个也不是严格的，能量-截断能的变化趋势趋于平滑即可。。。另外，截断能太大了也不好，我以前计算发现过，ECUT太大反而会让体系收敛缓慢，甚至scf不收敛。 |
8楼2013-11-14 10:24:12
|
非常感谢你的关注，U值是参考文献上的，赝势也和文献一样，不把Mn当成两类原子的话 也是一样的情况，用文献上的encut是没问题，我就是纳闷我怎么就测不出来，主要是我这个能量-截断能的变化趋势 没有趋于平滑，不知道跟我的磁构型有没有关系，文献是反铁磁构型，我就这么设置了初始磁矩，或者是我结构优化的不好，没有进入基态构型？但我是分步优化，过程都挺正常的呀。 350 -200.21188 400 -199.69089 450 -199.45863 500 -199.39341 550 -199.41708 600 -199.46676 650 -199.50511 700 -199.54608 |
9楼2013-11-17 18:26:15
10楼2013-11-17 18:29:15
50