24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2909  |  回复: 15
本帖产生 1 个 1ST强帖 ,点击这里进行查看

Gina88

木虫 (正式写手)


[交流] 【其他】vasp_band.exe处理的数据怎么最后不一样啊?

这个软件处理数据的时候,要求输入费米能级。
k点             费米能级=0                 费米能级=0.5           费米能级=1
0                         -0.28858700             -0.78858697           -1.2885870
0.01767765        -0.28196099             -0.78196096           -1.2819610
0.0353553258    -0.25929800             -0.75929797           -1.2592980
0.0530330203    -0.22188200             -0.72188199           -1.2218820
0.0707106963    -0.17048199             -0.67048198           -1.1704820

大家看看,最后的尾数怎么不一样啊?

[ Last edited by Gina88 on 2011-4-4 at 09:57 ]
回复此楼

» 猜你喜欢

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

» 抢金币啦!回帖就可以得到:

查看全部散金贴

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

夕阳西下

金虫 (著名写手)


Gina88(金币+1): 2011-04-04 21:52:35
最后是什么?
2楼2011-04-04 11:58:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzy870720z(金币+1): 谢谢交流 2011-04-04 13:48:49
Gina88(金币+2): 2011-04-04 21:52:27
不懂你要这到底是要什么?你要利用这个可执行文件运行,那你要对其中的最起码原理了解。网上随便找个程序,自己一运行,就提问,让别人不好回答吧。至少这不是常用的小程序。
3楼2011-04-04 12:48:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gina88

木虫 (正式写手)



franch(金币+1): 鼓励交流 2011-04-05 15:49:20
引用回帖:
Originally posted by youzhizhe at 2011-04-04 12:48:29:
不懂你要这到底是要什么?你要利用这个可执行文件运行,那你要对其中的最起码原理了解。网上随便找个程序,自己一运行,就提问,让别人不好回答吧。至少这不是常用的小程序。

k点             费米能级=0                 费米能级=0.5           费米能级=1
0                         -0.28858700             -0.78858697           -1.2885870
上面是那个软件产生的数据。
按我自己的理解,数据应该是:
k点             费米能级=0                 费米能级=0.5           费米能级=1
0                         -0.28858700             -0.78858700           -1.2885870

就是费米能级选用0,0.5,1,我觉得得到的能量本征值末尾几位数应该是一样的,可事实是不一样的。也就是对于K=0的点,如果费米能级=0.5,那K=0点的这条能带的本征能量应该为-0.28858700-0.5=-0.78858700,而实际的却是-0.78858697。

而当选费米能级=1时,我觉得应该是-1.28858700,而实际是-1.2885870。也就是小数点后的位数及精度好像是受费米能级的实际取值影响的。
4楼2011-04-04 21:59:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enola

捐助贵宾 (正式写手)


Gina88(金币+1): 2011-04-06 14:29:12
引用回帖:
Originally posted by Gina88 at 2011-04-04 21:59:47:
k点             费米能级=0                 费米能级=0.5           费米能级=1
0                         -0.28858700             -0.78858697           -1.2885870
上面是那个软件产生的数据。
...

这个软件好像会把价带顶(VASP默认最高占据态)移到0
5楼2011-04-06 11:56:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gina88

木虫 (正式写手)


引用回帖:
Originally posted by enola at 2011-04-06 11:56:01:
这个软件好像会把价带顶(VASP默认最高占据态)移到0

嗯,您说的是对的。
但是要自己填写费米能级的。

初学VASP,就瞎填个数,看看到底怎么移动。
结果却发现,我瞎填的数是整数,却会导致小数最后几位不同。也就是导致精度不一样。
不知道怎么回事,所以就发这个傻傻的问题了。
6楼2011-04-06 14:32:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
zzy870720z(金币+2): 谢谢指教 2011-04-07 12:04:04
zzy870720z(1ST强帖+1): 很详细的解释,呵呵,期待你更好的表现 2011-04-07 12:04:56
Gina88(金币+1): 2011-04-07 17:04:58
引用回帖:
Originally posted by Gina88 at 2011-04-04 09:52:40:
这个软件处理数据的时候,要求输入费米能级。
k点             费米能级=0                 费米能级=0.5           费米能级=1
0                         -0.28858700             -0.78858697          ...

计算出来的能带结构不是以费米能级为零点的,而是以真空能级为零点,但习惯上的能级图都是以费米能级为零点,所以会有一个平移。对一个体系的费米能级在outcar,是个固定的值,而不是像你这样乱取
7楼2011-04-07 09:32:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enola

捐助贵宾 (正式写手)


引用回帖:
Originally posted by 贺仪 at 2011-04-07 09:32:00:
计算出来的能带结构不是以费米能级为零点的,而是以真空能级为零点,但习惯上的能级图都是以费米能级为零点,所以会有一个平移。对一个体系的费米能级在outcar,是个固定的值,而不是像你这样乱取

贺仪,你好。。对于“以真空能级为零点”我觉得不对啊,在无限周期性的DFT的计算中,似乎不是真空能级啊?你有参考文献嘛?
不然我们可以把价带顶可以认为电离能啦?
8楼2011-04-07 12:46:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
zzy870720z(金币+2): 谢谢指点 2011-04-07 16:20:07
Gina88(金币+1): 2011-04-07 17:12:31
引用回帖:
Originally posted by enola at 2011-04-07 12:46:37:
贺仪,你好。。对于“以真空能级为零点”我觉得不对啊,在无限周期性的DFT的计算中,似乎不是真空能级啊?你有参考文献嘛?
不然我们可以把价带顶可以认为电离能啦?

你说的有道理,指示看到过这样的说法,我也没有文献。
我是这么理解的:在K-S方程中,每一项都有能量的来源,并没有加一个常数项。如果就两个电子来说,他的能量就是相互排斥能和交换关联能,如果只有一个电子呢?这些能量不就都为零了吗?这与单个电子放在真空里的能量也就是真空能级一样吗?从这点来看,解K-S方程的得到的能量应该就是真空能级为零点的了,不是等于也是近似吧,当然这不严谨,欢迎大家一起讨论。
就这个帖子的内容来说,计算的时候不管是不是以真空能级为零点,但总存在这样的参考零点,所以平移还是要的
零点也不是价带顶,你多画几个就知道了

[ Last edited by 贺仪 on 2011-4-7 at 13:53 ]
9楼2011-04-07 13:49:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haowenping

木虫 (著名写手)


引用回帖:
Originally posted by Gina88 at 2011-04-04 09:52:40:
这个软件处理数据的时候,要求输入费米能级。
k点             费米能级=0                 费米能级=0.5           费米能级=1
0                         -0.28858700             -0.78858697          ...

你用的那个程序我看是没有问题的呀!只是你要在OUTCAR中找到费米能级,然后再在把费米能级输进去就OK了!
10楼2011-04-07 16:08:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gina88

木虫 (正式写手)


引用回帖:
Originally posted by 贺仪 at 2011-04-07 13:49:47:
你说的有道理,指示看到过这样的说法,我也没有文献。
我是这么理解的:在K-S方程中,每一项都有能量的来源,并没有加一个常数项。如果就两个电子来说,他的能量就是相互排斥能和交换关联能,如果只有一个电 ...

费米能级不在价带顶,那它是真实费米能级吗?
11楼2011-04-07 17:13:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
Gina88(金币+1):谢谢参与
zzy870720z(金币+1): 鼓励交流 2011-04-08 08:26:55
最后一个费米能及改为1.0,,, 有效数字问题

其实数值计算出现这种情况很正常的,,,我觉得
12楼2011-04-07 18:18:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enola

捐助贵宾 (正式写手)



Gina88(金币+1):谢谢参与
引用回帖:
Originally posted by 贺仪 at 2011-04-07 13:49:47:
你说的有道理,指示看到过这样的说法,我也没有文献。
我是这么理解的:在K-S方程中,每一项都有能量的来源,并没有加一个常数项。如果就两个电子来说,他的能量就是相互排斥能和交换关联能,如果只有一个电 ...

恩。。
但是零点不是价带顶我不知道阿?因为VASP不是把最高占据态认为是费米能级的位置吗?平移后不就是拉?
13楼2011-04-09 10:26:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
franch(金币+1): 鼓励继续交流 2011-04-11 20:12:32
引用回帖:
Originally posted by enola at 2011-04-09 10:26:14:
恩。。
但是零点不是价带顶我不知道阿?因为VASP不是把最高占据态认为是费米能级的位置吗?平移后不就是拉?

最高占据态确实是费米能级的位置,但是没有人说最高占据就是价带顶呀,半导体物理里面有,对于本征半导体,费米能级处在价带顶和导带底的中间
14楼2011-04-09 10:38:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enola

捐助贵宾 (正式写手)


★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
franch(金币+1): 鼓励继续交流 2011-04-11 20:12:50
引用回帖:
Originally posted by 贺仪 at 2011-04-09 10:38:55:
最高占据态确实是费米能级的位置,但是没有人说最高占据就是价带顶呀,半导体物理里面有,对于本征半导体,费米能级处在价带顶和导带底的中间

对,有点搞混啦,哈哈哈。。。
我有个疑问啊,OUTCAR里有列出某个能带的占据态的,或是2,或是1,这个讲的是能带中电子的占据吗?但这个估计和ISMEAR有关?

[ Last edited by enola on 2011-4-11 at 17:41 ]
15楼2011-04-11 17:37:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guolianshun

金虫 (著名写手)


★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
gzqdyouxia(金币+1): 谢谢交流 2011-09-02 15:16:07
引用回帖:
7楼: Originally posted by 贺仪 at 2011-04-07 09:32:00:
计算出来的能带结构不是以费米能级为零点的,而是以真空能级为零点,但习惯上的能级图都是以费米能级为零点,所以会有一个平移。对一个体系的费米能级在outcar,是个固定的值,而不是像你这样乱取

很对,输出文件outcar里有一个费米能级,最下面那个就是。但是貌似文章里面也是很难讲清楚。
16楼2011-09-02 10:50:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Gina88 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见