24小时热门版块排行榜    

查看: 4346  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 明天应该可查了!? +4 chengyan1220 2026-08-23 4/200 2026-08-24 08:47 by gloomy6159
[基金申请] 能否退出参与的面上项目解除限项 +6 koalala 2026-08-24 6/300 2026-08-24 08:21 by 正能量1斤
[基金申请] 范进中举一文的中心思想 +5 炎黄贵胄 2026-08-22 6/300 2026-08-24 06:34 by sincosx
[基金申请] 建议基金发布提前给出明确的时间点 +11 kulium 2026-08-21 14/700 2026-08-24 02:12 by 流流伤
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见