24小时热门版块排行榜    

查看: 7973  |  回复: 13

zhjh888888

银虫 (小有名气)


[交流] 【求助】vasp计算Partial Charge,能带交叉的如何确定IBAND?

请问要做(图1)这种分析,能带交叉的地方怎么处理?比如这里(图2)两条红色的交叉 ,IBAND怎么设?vasp数据用origin画出一般这样 (图3)
origin里一条曲线并不一定是在一条能带上!红色48和红色49为一条能带,却在origin里的至少两条曲线上!先谢谢了!

图1

本图来自参考文献 APPLIED PHYSICS LETTERS 92, 173118 (2008 ),致谢!
图2

图3


[ Last edited by zhjh888888 on 2011-3-30 at 19:22 ]
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理相关文档 VASP算例 vasp Mr Wang收藏的投稿交流
待办事宜

» 猜你喜欢

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

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

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

lfhuang

木虫 (著名写手)



小木虫: 金币+0.5, 给个红包,谢谢回帖
这个问题可以解决的,只要楼主先把能带进行分类,具体程序请看:
http://blog.sciencenet.cn/home.php?mod=space&uid=345795&do=blog&quickforward=1&id=641013
12楼2014-01-09 20:33:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

bingmou

金虫 (著名写手)


★ ★
sunyang1988(金币+2): 谢谢交流 2011-03-30 22:59:40
zhjh888888(金币+5): 谢谢解答! 2011-03-31 10:55:07
一般是GAMMA点的进行编号吧,而且应该是指定K点的
2楼2011-03-30 21:43:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhjh888888

银虫 (小有名气)


引用回帖:
Originally posted by bingmou at 2011-03-30 21:43:09:
一般是GAMMA点的进行编号吧,而且应该是指定K点的

非常感谢回复!你的意思是不是:在输出文件PARCHG.0047.ALLK里,47是能带在GAMMA点编号,不管Origin里其它K点这条能带的号码是多少。要计算哪条能带对应的Partial Charge,只要确定其在GAMMA点的编号就可以了,不管Origin里这条能带在其它K点的编号是不是相同?
3楼2011-03-31 10:50:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)



zzy870720z(金币+1): 谢谢提示 2011-03-31 17:24:16
引用回帖:
Originally posted by zhjh888888 at 2011-03-31 10:50:21:
非常感谢回复!你的意思是不是:在输出文件PARCHG.0047.ALLK里,47是能带在GAMMA点编号,不管Origin里其它K点这条能带的号码是多少。要计算哪条能带对应的Partial Charge,只要确定其在GAMMA点的编号就可以了, ...

在VASP中是有一个KPUSE来指定K点的,手册的85页有详细的介绍
4楼2011-03-31 14:42:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhjh888888

银虫 (小有名气)


引用回帖:
Originally posted by bingmou at 2011-03-31 14:42:55:
在VASP中是有一个KPUSE来指定K点的,手册的85页有详细的介绍

谢谢,这个参数我用了。
.....................................................................
# ! plot band decomposed chargedensity
LPARD = .TRUE.
IBAND =   49
KPUSE = 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 !X to GAMMA
LSEPB = .TRUE.
LSEPK = .FALSE.
.....................................................................
以上是INCAR里的部分内容
问题的关键是IBAND的确定,所设的49(G点考察能带的标号,在X点这条能带的标号不是49了),是否就是要考察的边缘态能带的!origin里一条曲线有的不是一条能带!或者说同一条能带不同K点时曲线标号有时不一样(见一楼图3)!某条能带的标号应该按哪个K点来确定?GAMMA?
下边是KPOINTS
k-points along high symmetry lines,system:ZGNR
21 ! 21 intersections
Line-mode
rec
  0.5   0.0   0.0   ! X
  0.0   0.0   0.0   ! gamma
5楼2011-03-31 16:25:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)



zzy870720z(金币+1): 谢谢指教 2011-03-31 21:16:00
引用回帖:
Originally posted by zhjh888888 at 2011-03-31 16:25:54:
谢谢,这个参数我用了。
.....................................................................
# ! plot band decomposed chargedensity
LPARD = .TRUE.
IBAND =   49
KPUSE = 1 2 3 4 5 6 7 8 9 10 1 ...

为什么你要综合那么多个K点的呢?我个人觉得在不同k点同一条能带的性质并不一样,电荷密度分布也不一样,所以不能把整条能带的所有K点的电荷密度叠加进行分析。
当然你可以试试gamma点49和X点的47,也就是你关注的那一条能带,看看两者的电荷分布是否相同或者类似

[ Last edited by bingmou on 2011-3-31 at 19:27 ]
6楼2011-03-31 19:25:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhjh888888

银虫 (小有名气)


引用回帖:
Originally posted by bingmou at 2011-03-31 19:25:13:
为什么你要综合那么多个K点的呢?我个人觉得在不同k点同一条能带的性质并不一样,电荷密度分布也不一样,所以不能把整条能带的所有K点的电荷密度叠加进行分析。
当然你可以试试gamma点49和X点的47,也就是你关 ...

很感谢你的回复!谢谢!
我是参考上面那篇文献和下边的总结《用VASP进行Partial Charge分析实例》http://hi.baidu.com/datavault/bl ... 678cc537d3ca50.html
VASP可以得到某个能量范围、某个K点或者某个特定的态所对应的电荷密度。在文献中最常见的是价带顶部,导带底部,表面态或者局域态所对应的Partial Charge。通过分析这些态所对应的Partial Charge,可以得到体系的一些性质,比如局域态具体的是局域在哪个原子上等。我是想看边缘态的电荷分布。究竟是纳米带上哪个部位的电子激发的某个带。是不是应该把整条能带的所有计算的K点的电荷密度叠加?

[ Last edited by zhjh888888 on 2011-3-31 at 22:05 ]
7楼2011-03-31 22:01:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)


★ ★
zzy870720z(金币+2): 谢谢指教 2011-04-01 16:36:00
zhjh888888(金币+10): 2011-04-01 18:11:52
引用回帖:
Originally posted by zhjh888888 at 2011-03-31 22:01:36:
很感谢你的回复!谢谢!
我是参考上面那篇文献和下边的总结《用VASP进行Partial Charge分析实例》http://hi.baidu.com/datavault/bl ... 678cc537d3ca50.html
VASP可以得到某个能量 ...

我觉得直接计算gamma点的VBM的电荷密度或者CBM的电荷密度就可以了,很多时候激发出来应该是这里的电子
8楼2011-04-01 16:23:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhjh888888

银虫 (小有名气)


引用回帖:
Originally posted by bingmou at 2011-04-01 16:23:38:
我觉得直接计算gamma点的VBM的电荷密度或者CBM的电荷密度就可以了,很多时候激发出来应该是这里的电子

好的,我再结合一下文献!谢谢你了!讨论的过程学会了很多东西!谢谢!
9楼2011-04-01 18:14:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

国花1989

铜虫 (小有名气)



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
9楼: Originally posted by zhjh888888 at 2011-04-01 18:14:04
好的,我再结合一下文献!谢谢你了!讨论的过程学会了很多东西!谢谢!...

楼主,你好,你的问题解决了吗??我想问,边缘态看出来的啊??那条能带是边缘态,看编号是多少怎么看啊??谢谢喽,我真心求解啊~~~
10楼2013-12-03 22:28:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmumu012

铁虫 (初入文坛)



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
5楼: Originally posted by zhjh888888 at 2011-03-31 16:25:54
谢谢,这个参数我用了。
.....................................................................
# ! plot band decomposed chargedensity
LPARD = .TRUE.
IBAND =   49
KPUSE = 1 2 3 4 5 6 7 8 9 10 11 1 ...

请问,导带底不在Gamma点的电荷密度怎么设置啊,IBAND参数还是G点的排序吗?KPOINTS文件怎么设置啊
11楼2014-01-09 17:04:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

refors

木虫 (小有名气)



小木虫: 金币+0.5, 给个红包,谢谢回帖
只看Gamma点的话,3图楼主的问题仍然没有解决,这个带在Gamma点上的两个值差距为何如此之大?而且origin的结果一个横坐标只对应一个能量值,楼主如何确定两个是一个带的?
13楼2016-10-25 16:28:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sin81

金虫 (正式写手)


学习了。
14楼2018-12-12 20:43:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhjh888888 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工 +4 安全上岸! 2026-03-16 4/200 2026-03-17 14:02 by 勇敢太监王公公
[考研] 08工科 320总分 求调剂 +4 梨花珞晚风 2026-03-17 4/200 2026-03-17 13:38 by houyaoxu
[考研] 268求调剂 +7 一定有学上- 2026-03-14 8/400 2026-03-17 13:10 by zz820
[考博] 26申博 +4 八旬速览 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 341求调剂 +4 番茄头--- 2026-03-10 4/200 2026-03-13 23:12 by JourneyLucky
[考研] 332求调剂 +3 zjy101327 2026-03-11 6/300 2026-03-13 22:48 by JourneyLucky
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 333求调剂 +3 球球古力 2026-03-11 3/150 2026-03-13 21:27 by JourneyLucky
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[考研] 274求调剂0856材料化工 +12 z2839474511 2026-03-11 13/650 2026-03-13 10:39 by peike
[考研] 求调剂 资源与环境 285 +3 未名考生 2026-03-10 3/150 2026-03-13 10:31 by houyaoxu
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
信息提示
请填处理意见