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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】请教vasp如何使本征值按每条能带输出
191/2返回列表12下一页
查看: 3539  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 jRl3mE6ddZGq 2026-08-03 8/400 2026-08-04 06:17 by k7OM8YghWbkC
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 jRl3mE6ddZGq 2026-08-03 8/400 2026-08-04 06:15 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 6/300 2026-08-04 06:10 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:46 by k7OM8YghWbkC
[有机交流] 一个有机合成实验室都需要哪些设备？ 50+3 kf2781974 2026-07-31 9/450 2026-08-04 03:56 by kf2781974
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 12/600 2026-08-04 03:17 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀？ +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上提前没消息，有中的吗 +14 archvillain 2026-08-02 16/800 2026-08-03 12:23 by fuweiguochen
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[教师之家] 基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 12/600 2026-08-03 08:22 by syl200707
[基金申请] 微信指数没变化，科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液，是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了，估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
关闭
请填处理意见
关闭