24小时热门版块排行榜    

查看: 1940  |  回复: 13
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

nxc0907

[交流] 【求助】vasp计算单个原子能量总能是正值,怎么办?

我计算zn原子的能量,用的是paw_GGA势,可计算出的能量怎么都是正的,这有可能么?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nxc0907

INCAR:
ISTART=0;ICHARG=2
ISMEAR=0
PREC=Accurate


kpoints:
0
G
1 1 1
0 0 0
原来还做了自旋极化,也是正的。
3楼2009-09-03 11:53:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nxc0907

后来用paw_PBE做就没问题了,其他都没变,难道是赝势有问题?
4楼2009-09-03 11:54:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nxc0907


zdhlover(金币+1,VIP+0):欢迎常来交流讨论 9-3 23:39
引用回帖:
Originally posted by fanchen021 at 2009-9-3 14:49:
算单个原子能量,必然要加 自旋

这要看外层电子的排布了,都配对了就不用加了吧。
11楼2009-09-04 13:56:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nxc0907


zdhlover(金币+1,VIP+0):欢迎常来交流讨论 9-3 23:39
引用回帖:
Originally posted by martin7 at 2009-9-3 22:25:
paw_pbe
Paw_GGa?
这两个你怎么区分的。。。。。
另建议把结果贴出来
再则单原子vasp直接有实例,自己看看

我们服务器上的赝势文件直接有这种区分啊,结果就是用PAW_GGA势的话OSZCAR中最后的能量就是正的,但pbe就是负的很正常,另外不管用哪个势收敛的都很快,迭代几步就出来了。
12楼2009-09-04 14:00:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nxc0907


zdhlover(金币+1,VIP+0):欢迎常来交流讨论 9-3 23:39
引用回帖:
Originally posted by wuchenwf at 2009-9-3 23:17:
加自旋会精确一些,但是一般不会出现由正值到负值的这种情况。

我的这个不是加不加自选的问题,是势的问题,用paw—GGA加不加自选都是正的,用pbe则是负的能量。
13楼2009-09-04 14:04:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nxc0907

另外对于OSZCAR中的能量E究竟是什么含义?我看电子的迭代过程中这个能量一会正一会负,但一般最后收敛时都是负的,师兄说正的能量说明计算的结果就是错的,但为什么是这样呢?那中间迭代的过程中怎么有正的能量呢?
14楼2009-09-04 14:12:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 nxc0907 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见