24小时热门版块排行榜    

查看: 3517  |  回复: 13

gleerat

木虫 (正式写手)

[交流] 【求助】如何使用VASP画费米面。 已有6人参与

不知道如何用VASP画费米面,请教用VASP做过费米面计算的同仁来给指一条明路。
1、VASP计算费米面的一般步骤;
2、对于计算结果的处理的一般步骤;
3、用来画费米面的工具。

我知道PWScf和Wine2K可以画费米面,而且Wine2K画费米面还相当简单,跟XCrysden的结合很紧密,不知道它是如何处理的。XCrysden对于VASP的支持就没有Wine2K那么好了,画费米面的时候还需要一个4×4×4的矩阵,而这个矩阵如何得来的根本就不知道。
所以,如果用XCrysden画图的话,请麻烦多多指教一下得到4×4×4的矩阵的方法。

先行谢谢各位了。
回复此楼
明察、慎思、笃行
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)


gavinliu7390(金币+1):谢谢交流! 2010-06-30 21:54:09
VASP不能画费米面。
2楼2010-06-30 20:54:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)


zxzj05(金币+1):鼓励有价值交流信息! 2010-06-30 22:52:09
引用回帖:
Originally posted by stractor at 2010-06-30 20:54:48:
VASP不能画费米面。

从官方论坛http://cms.mpi.univie.ac.at/vasp-forum/forum_viewtopic.php?4.4801的情况看来应该是可以的。
只是不知道如何画。

[ Last edited by gleerat on 2010-6-30 at 22:23 ]
明察、慎思、笃行
3楼2010-06-30 22:18:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beyondstar

木虫 (正式写手)

★ ★
lvjian8596(金币+1):感谢交流! 2010-06-30 22:38:26
aylayl08(金币+1):谢谢提供 2010-07-02 16:52:24
看这里,如果非要想画,需要自己做一个后处理程序;不过应该不难
http://bbs.zixia.net/disparticle ... ID=37950&pos=13
4楼2010-06-30 22:20:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Prettyswan

木虫 (小有名气)

★ ★
lvjian8596(金币+1):感谢交流! 2010-06-30 22:38:43
aylayl08(金币+1):谢谢提供 2010-07-02 16:52:35
在官方论坛上我看到admin是这样说的“you need a dense k-grid to interpolate the Fermi surface from the EIGENVAL file”,下面另外一个人又说“I think what you want to do is to plot the isosurface of the electron charge density, is that right? ”。
我也听一些人说,电荷密度文件是可以来画的。
5楼2010-06-30 22:29:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)

引用回帖:
Originally posted by beyondstar at 2010-06-30 22:20:57:
看这里,如果非要想画,需要自己做一个后处理程序;不过应该不难
http://bbs.zixia.net/disparticle ... ID=37950&pos=13

关于这个我也知道。这里也有一篇用PWScf做费米面的文章,原理都能看懂,就是那个4×4×4的矩阵不好弄。http://zhaoyunawuli.blog.163.com ... 762008111471211689/
现在关键的问题是,如果用本征值文件来画的话,如何将EIGENVAL转化为XCrysden所用的矩阵文件.bxsf。(我不会编写程序或者脚本)
明察、慎思、笃行
6楼2010-06-30 22:39:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)

★ ★
zxzj05(金币+1):感谢交流 2010-06-30 22:52:40
gavinliu7390(金币+1):谢谢交流! 2010-07-01 12:03:39
vasp貌似不能画费米面费米面要在静态计算时取三维的K点,KPOINTS里没有提供这一方案,只有画高对称K点能带图的Line-mode,本人新手不当之处敬请指教
心之所向无惧无悔
7楼2010-06-30 22:39:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Prettyswan

木虫 (小有名气)


zzy870720z(金币+1):费米面的定义就是那个,呵呵,后面的不太理解,希望能进一步解释,谢谢 2010-07-01 08:04:19
引用回帖:
Originally posted by sunyang1988 at 2010-06-30 22:39:54:
vasp貌似不能画费米面费米面要在静态计算时取三维的K点,KPOINTS里没有提供这一方案,只有画高对称K点能带图的Line-mode,本人新手不当之处敬请指教

我认为,K空间中能量等于费米能的等能面就是费米面,所以用把本征能量相同且等于费米能的K点做成的面在K空间中画出来就行了,不必非要能带计算的本征值的。
如有错误,敬请指正。
8楼2010-06-30 22:53:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)


gavinliu7390(金币+1):谢谢交流! 2010-07-01 12:03:51
引用回帖:
Originally posted by Prettyswan at 2010-06-30 22:53:09:

我认为,K空间中能量等于费米能的等能面就是费米面,所以用把本征能量相同且等于费米能的K点做成的面在K空间中画出来就行了,不必非要能带计算的本征值的。
如有错误,敬请指正。

是这样啊,但是要画出三维的等能面,K点应该很密。而且并不是直接算费米面上的K点能量,应该是画出三维的K点能量图与费米能相交的曲面得到等能面吧。呵呵,我也不太懂。欢迎批评指正
心之所向无惧无悔
9楼2010-07-01 09:34:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Prettyswan

木虫 (小有名气)

★ ★
zzy870720z(金币+2):谢谢交流讨论 2010-07-01 19:16:14
引用回帖:
Originally posted by sunyang1988 at 2010-07-01 09:34:30:

是这样啊,但是要画出三维的等能面,K点应该很密。而且并不是直接算费米面上的K点能量,应该是画出三维的K点能量图与费米能相交的曲面得到等能面吧。呵呵,我也不太懂。欢迎批评指正

在Wine2K中,画费米面用的是穿过费米能的带的K点。我想用VASP也该是同样的道理,计算密集的K点的能量本征值,然后找出穿过费米能的带的K点的坐标,在三维的K空间中画出这些K点组成的面来,这应该就是费米面。
10楼2010-07-01 17:10:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gleerat 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见