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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】vasp计算Partial Charge，能带交叉的如何确定IBAND？
141/1返回列表
查看: 8624  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 10/500 2026-08-17 23:43 by sytp
[基金申请] 93BebMhtakh前后11位开头都是大写 +3 且听虎啸 2026-08-17 4/200 2026-08-17 22:31 by ccec101
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +36 医学老男孩 2026-08-13 85/4250 2026-08-17 21:23 by lcmdlut
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 时间戳变了，能看出什么问题？ +4 基诺咪客 2026-08-17 4/200 2026-08-17 16:10 by Vivilian
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +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
[基金申请] 是这周出结果还是下周出结果？ +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
[基金申请] 重要来源：本周末出结果 +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
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
信息提示
关闭
请填处理意见
关闭