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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】vasp计算Partial Charge，能带交叉的如何确定IBAND？
141/1返回列表
查看: 6885  |  回复: 13
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 zhjh888888 的 5 个金币

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返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件) 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 实验室太吵闹，无法安静学习，怎么办？ +7 utahh 2024-05-31 12/600 2024-06-02 01:23 by dking22
[教师之家] 双非本科毕业论文，气人 +7 河西夜郎 2024-05-27 12/600 2024-06-01 22:12 by iwdw2012
[基金申请] 化学口B0109(高分子合成)，拿青年基金一般需要怎样的文章水平？ +13 salmon95 2024-05-30 22/1100 2024-06-01 21:22 by salmon95
[硕博家园] 每天学术时间不能保证，能保证的只有： +10 hahamyid 2024-05-27 10/500 2024-06-01 21:11 by 小小芝麻影
[教师之家] 研究生在毕业答辩时挂了，遗憾 +10 spacr2022 2024-05-27 13/650 2024-06-01 17:45 by keyaner23
[考博] 24年博士招生 +8 abinit432 2024-05-27 10/500 2024-06-01 17:38 by czp97
[考博] 24or25材料专业申博 +4 农夫三拳有点痛 2024-05-30 11/550 2024-06-01 14:45 by Napoleonsky
[硕博家园] 讨厌酒桌文化 +6 hahamyid 2024-05-31 6/300 2024-06-01 09:03 by lizhiweiouc
[硕博家园] 哈工大硕博招生!博士每月入学! +4 nailooo 2024-05-30 5/250 2024-06-01 06:47 by anevay
[有机交流] 苯酚和炔溴反应成醚 6+3 TONGMEIMEI 2024-05-26 12/600 2024-05-31 14:41 by TONGMEIMEI
[文学芳草园] 对对子啊 +5 天若孤独 2024-05-29 7/350 2024-05-31 09:00 by wjykycg
[电化学] 有锂金属软包电池代工公司代做软包嘛 +5 yyywj995 2024-05-26 6/300 2024-05-30 08:33 by 那年的旧城
[硕博家园] 又想换工作 +15 brightmj 2024-05-27 20/1000 2024-05-29 23:25 by zyqchem
[博后之家] 2024公派博后申请 +4 326lhpqk 2024-05-27 5/250 2024-05-29 20:03 by @古月胡
[论文投稿] 高手朋友推荐比较容易投稿和录用的SCI期刊，不投稿收费SC，对分区没有要求 5+3 xintangren 2024-05-28 4/200 2024-05-29 10:46 by xintangren
[有机交流] 机理求助 120+4 15147165026 2024-05-26 10/500 2024-05-28 14:42 by 江东闲人
[基金申请] 面上基金会评专家，有回避机制吗？ +4 huang1991js 2024-05-27 4/200 2024-05-27 19:08 by 星火12
[硕博家园] 我是很理想化一人 +6 hahamyid 2024-05-26 6/300 2024-05-27 18:13 by 大飞鱼鱼鱼
[基金申请] 感觉自然基金限制通过比例就是有点扯，学学B口，化学学部，不限制比例。 +10 wsjing 2024-05-26 14/700 2024-05-27 11:57 by kanmiaolucky
[硕博家园] 周日 +6 1加油哦棒 2024-05-26 9/450 2024-05-27 10:30 by hahamyid
信息提示
关闭
请填处理意见
关闭