应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】vasp计算Partial Charge，能带交叉的如何确定IBAND？
141/1返回列表
查看: 8593  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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收藏的投稿交流
待办事宜

» 猜你喜欢

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

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

查看全部散金贴
1楼 2011-03-30 16:51:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

lfhuang

木虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
这个问题可以解决的，只要楼主先把能带进行分类，具体程序请看：
http://blog.sciencenet.cn/home.php?mod=space&uid=345795&do=blog&quickforward=1&id=641013
赞一下(1人) 回复此楼
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点的
赞一下(1人) 回复此楼
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页有详细的介绍
赞一下(1人) 回复此楼
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 ]
赞一下(1人) 回复此楼
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的电荷密度就可以了，很多时候激发出来应该是这里的电子
赞一下(1人) 回复此楼
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 的主题更新
141/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +9 Ldrop2023 2026-08-13 9/450 2026-08-16 13:24 by wukun_007
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +11 家与远方 2026-08-10 16/800 2026-08-16 10:28 by ray43
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] FileCode能看出啥？ +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
关闭
请填处理意见
关闭