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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】请教vasp如何使本征值按每条能带输出
191/2返回列表12下一页
查看: 3196  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

TIGERYZZ

金虫 (小有名气)

[交流] 【求助】请教vasp如何使本征值按每条能带输出

vasp在计算能带后，EIGENVALUE文件输出的本质值是按能量大小排列的，请问如何确定哪些本征值是一条能带上的，如pi 带。在计算能带时是否可以按不同能带输出本征值？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-08-21 15:49:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

valenhou001

至尊木虫 (职业作家)
TIGERYZZ(金币+3):谢谢指导~~ 2010-08-25 08:52:41
这样的图，通过vasp也是可以弄出来了的。我前面提到了有个fat band的概念。操作的方式，可参考abinit里面。这个fat band 的概念好像是O.K. Anderson他们组最早提出来的，用来标准每个带的轨道贡献。在vasp中输出了PROCAR，里面给出了每个原子的每个轨道对每个k点每个本征值的贡献。在计算能带结构时，也加上LORBIT=11(针对PAW)方法，然后在能带结构计算完了得到PROCAR，想办法处理PROCAR里面的数据.

VASP弄的fat band图的例子，可参考这篇文献：
http://arxiv.org/PS_cache/arxiv/pdf/0908/0908.3502v1.pdf
赞一下(1人) 回复此楼
15楼2010-08-24 21:44:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meinvren

禁虫 (小有名气)
本帖内容被屏蔽
17楼2010-08-24 23:09:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

锐利的碎片

木虫 (正式写手)

star watcher
★
zzy870720z(金币+1):谢谢提示，最好给出链接，呵呵 2010-08-21 16:54:33
TIGERYZZ(金币+2): 2010-08-21 19:24:39
默认是按给出的k点输出能量。处理一下EIGENVAL就可以了。可以自己写，以前也有人发过。
赞一下(1人) 回复此楼
2楼2010-08-21 16:48:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

TIGERYZZ

金虫 (小有名气)
引用回帖:
Originally posted by 锐利的碎片 at 2010-08-21 16:48:26:
默认是按给出的k点输出能量。处理一下EIGENVAL就可以了。可以自己写，以前也有人发过。

我可能说的不清楚，默认按K点输出能量不错，然后不同k点的本征值输出是：
1 Energy
2 Energy
3 Energy
……
这个是按本征值能量从小到大排的。
那就有可能Gamma点的3 是 pi 带 而其他k点处的3并不是pi 带。
那么这个时候怎么才能确认哪些本征值文件属于一条能带？
赞一下 回复此楼
3楼2010-08-21 19:39:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
★ ★
qasd(金币+2):谢谢 2010-08-21 20:43:35
TIGERYZZ(金币+2):谢谢 2010-08-22 09:00:13
TIGERYZZ(金币+2):侯老师能够详细介绍一下吗？谢谢 2010-08-22 09:03:34
类似abinit做一个fat-band band structure。
赞一下(1人) 回复此楼
4楼2010-08-21 20:29:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zengchangan

铁虫 (初入文坛)

呵呵 也遇到相关问题

★
zzy870720z(金币+1):提议不错，祝你们尽快解决 2010-08-21 22:09:58
TIGERYZZ(金币+2):谢谢 2010-08-22 09:00:38
好像要找出高对称的K点 按K点画 可以通过EIGENVALUE文件读出每个K点的能量 进而分析K点走势 再用Origin画出能带图
赞一下(1人) 回复此楼
苦心专研奋力拼搏，天道酬勤，好好加油吧！
5楼2010-08-21 21:12:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)
★ ★
zzy870720z(金币+2):谢谢专家指导，期待你更好的表现 2010-08-22 13:11:20
TIGERYZZ(金币+3):您的意思是说将第n个本征值合在一起就是一条能带吗？但是这条能带会在不同k点处表现不同的特征，譬如pi或sigma。这样是一条能带吗？ 2010-08-22 14:08:40
看EIGENVAL这个输出文件。它里面给出的就是自己所选定的各个K点的一些列能带的本征值。正如你所说的，EIGENVAL里面是给定一个K点，然后就是各个能带的本征值。那么你要画出第n条，比如第5条能带，你就将每个K点的第5条能带值取出来。要确定哪些能带时pi带，我猜想要做原子的态密度分析，确定哪些原子轨道会形成pi键。
“在计算能带时是否可以按不同能带输出本征值？”。当然能，你已知道EIGENVAL保存数据规则，这只需要对EIGENVAL文件写个程序处理一下。

[ Last edited by stractor on 2010-8-22 at 12:33 ]
赞一下(1人) 回复此楼
6楼2010-08-22 12:30:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

TIGERYZZ

金虫 (小有名气)
引用回帖:
Originally posted by stractor at 2010-08-22 12:30:07:
看EIGENVAL这个输出文件。它里面给出的就是自己所选定的各个K点的一些列能带的本征值。正如你所说的，EIGENVAL里面是给定一个K点，然后就是各个能带的本征值。那么你要画出第n条，比如第5条能带，你就将每个K点的 ...

楼上还有2楼提出的处理本征值文件，这是个容易的事情。
但是如楼上说的，如果我将不同k点的第5个本征值合在一起，这是一条能带吗？毕竟这些本征值在不同k点可能表征处不同的pi或sigma特性。我不知道vasp可不可以做这样的是，有没有可能输出的本征值文件，那就是不同k点的第n个本征值合在一起，表现相同的能带特性（譬如pi），我觉得这样才可以称为一条能带。

[ Last edited by TIGERYZZ on 2010-8-22 at 14:17 ]
赞一下 回复此楼
7楼2010-08-22 14:11:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)
TIGERYZZ(金币+2):给您回复了短信息，这样的一条能带并不是只表现pi或sigma，在不同k点可能表现不同特征，所以我依然觉得不是一条能带。 2010-08-23 10:45:30
引用回帖:
Originally posted by TIGERYZZ at 2010-08-22 14:11:06:

楼上还有2楼提出的处理本征值文件，这是个容易的事情。
但是如楼上说的，如果我将不同k点的第5个本征值合在一起，这是一条能带吗？毕竟这些本征值在不同k点可能表征处不同的pi或sigma特性。我不知道vasp可不可 ...

是的同一条能带。对不同K点的第m个点画图，就是同一条能带！

[ Last edited by stractor on 2010-8-23 at 10:30 ]
赞一下 回复此楼
8楼2010-08-23 10:23:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

TIGERYZZ

金虫 (小有名气)

传个图片说下（来自 PRL 96 026402）
拿图中的最高价带来说明。在K点处是Pi带，而G点处是sigma带。
而vasp的本征值输出中，这两个点的肯定是一个number的。
如果按楼上的说法，这两个点就应该是一条带，那么一条带就表现两种特征了。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

嘻嘻哈哈yuan
9楼2010-08-23 10:56:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hgye

金虫 (小有名气)
TIGERYZZ(金币+2):谢谢 2010-08-23 15:14:03
不同K点第m个点连线 不一定 是同一条能带。
楼主想要的效果vasp不可能做到，软件所得本征值只是数学计算结果，不可能从成键机理上分析哪些点属于一条能带。可以画出能带图通过各能带的伸展自己分析。
赞一下 回复此楼
10楼2010-08-23 11:03:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 TIGERYZZ 的主题更新
191/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[考研] 本人考085602 化学工程 专硕 +7 不知道叫什么！ 2026-03-15 8/400 2026-03-15 20:11 by 棒棒球手
[考研] 326求调剂 +3 上岸的小葡 2026-03-15 4/200 2026-03-15 18:50 by 无际的草原
[考研] 085600材料与化工 求调剂 +11 enenenhui 2026-03-13 12/600 2026-03-14 21:48 by laoshidan
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 求调剂，药学 +3 归零lbm 2026-03-09 5/250 2026-03-14 02:21 by JourneyLucky
[考研] 材料与化工 一志愿山大 321分 求调剂 +7 每天散步 2026-03-09 8/400 2026-03-14 02:18 by JourneyLucky
[考研] 求调剂！ +4 朔朔话 2026-03-09 4/200 2026-03-14 01:38 by JourneyLucky
[考研] 295复试调剂 +5 简木ChuFront 2026-03-09 5/250 2026-03-14 01:29 by JourneyLucky
[考研] 0703求调剂 +7 jtyq001 2026-03-10 7/350 2026-03-14 01:06 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-10 3/150 2026-03-14 00:30 by JourneyLucky
[考研] 0703化学调剂 +4 快乐的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 泣血叩求调剂恩，愿以丹心报师恩 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] （081700）化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 274求调剂0856材料化工 +12 z2839474511 2026-03-11 13/650 2026-03-13 10:39 by peike
[考研] 290求调剂 +3 柯淮然 2026-03-10 8/400 2026-03-11 13:48 by 柯淮然
[考研] 一志愿山东大学，总分327，英语二79，有论文，有竞赛，已过四六级 +3 木木目目1 2026-03-09 3/150 2026-03-09 19:52 by yuningshan
信息提示
关闭
请填处理意见
关闭