24小时热门版块排行榜    

查看: 3194  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料284求调剂,一志愿郑州大学英一数二专硕 +15 想上岸的土拨鼠 2026-02-28 15/750 2026-03-02 20:13 by hypershenger
[考研] 中国科学技术大学材料与化工281求调剂,有科研和获奖经历 +5 wsxw 2026-03-02 5/250 2026-03-02 20:12 by hypershenger
[考研] 302材料工程求调剂 +5 Doleres 2026-03-01 6/300 2026-03-02 19:53 by 张晓芳0105
[考研] 材料工程274求调剂 +5 Lilithan 2026-03-01 5/250 2026-03-02 19:39 by caszguilin
[考研] 课题组接收材料类调剂研究生 +4 gaoxiaoniuma 2026-02-28 5/250 2026-03-02 19:33 by Xinyuu
[考研] 材料085601调剂 +5 多多子. 2026-03-02 5/250 2026-03-02 19:15 by zhukairuo
[考研] 085600 英一数二272求调剂 5+6 vida_a 2026-03-01 16/800 2026-03-02 19:13 by zhukairuo
[考研] 化工京区271求调剂 +6 11ing 2026-03-02 6/300 2026-03-02 18:52 by caszguilin
[考研] 一志愿山东大学材料与化工325求调剂 +5 半截的诗0927 2026-03-02 5/250 2026-03-02 18:37 by 明亮9527
[考研] 一志愿东北大学材料专硕328,求调剂 +3 shs1083 2026-03-02 3/150 2026-03-02 17:27 by houyaoxu
[考研] 材料化工调剂 +12 今夏不夏 2026-03-01 14/700 2026-03-02 16:09 by 今夏不夏
[考研] 265分求调剂不调专业和学校有行学上就 +6 礼堂丁真258 2026-02-28 9/450 2026-03-02 12:04 by 52hz~~
[考研] 322求调剂 +3 熊境喆 2026-03-01 3/150 2026-03-02 08:44 by houyaoxu
[考研] 0857调剂 +4 一ll半 2026-02-28 5/250 2026-03-02 02:33 by 908055542
[考研] 272求调剂 +6 田智友 2026-02-28 6/300 2026-03-01 21:40 by 公瑾逍遥
[基金申请] 刚录用,没有期刊号,但是在线可看的论文可以放为代表作吗 10+3 arang1 2026-03-01 3/150 2026-03-01 16:43 by babero
[考研] 304求调剂 +6 曼殊2266 2026-02-28 7/350 2026-03-01 15:14 by wjLi2017
[考研] 311求调剂 +9 南迦720 2026-02-28 10/500 2026-03-01 10:55 by sunny81
[硕博家园] 2025届双非化工硕士毕业,申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
[论文投稿] 求助coordination chemistry reviews 的写作模板 10+3 ljplijiapeng 2026-02-27 4/200 2026-03-01 09:07 by babero
信息提示
请填处理意见