24小时热门版块排行榜    

查看: 3011  |  回复: 18

TIGERYZZ

金虫 (小有名气)

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

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

» 猜你喜欢

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

已阅   回复此楼   关注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
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就可以了。可以自己写,以前也有人发过。
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。
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画出能带图
苦心专研奋力拼搏,天道酬勤,好好加油吧!
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 ]
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朵)

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 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见