24小时热门版块排行榜    

CyRhmU.jpeg
查看: 7633  |  回复: 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的回帖
普通回帖
★ ★
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的回帖

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的回帖

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的回帖
★ ★
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 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见