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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】能带图如何画出呢？
北京石油化工学院2026年研究生招生接收调剂公告
91/1返回列表
查看: 2867  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

匿名

用户注销 (正式写手)
本帖仅楼主可见

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-11-30 16:13:43
已阅   同方向广播   申请1ST强帖   回复此楼   编辑   删除   查看我的主页

youzhizhe

荣誉版主 (职业作家)

有志者

优秀版主优秀版主
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1):谢谢交流 2010-12-01 00:47:28
引用回帖:
Originally posted by 晨曦3999 at 2010-11-30 16:13:43:
采用VASP软件计算出相应的能带后，采用什么方法可以将能带画出来呢？希望各位大侠能给出详细的步骤！非常感谢！

你都计算好了 那就用origin，一切OK。
赞一下(2人) 回复此楼
微博http://weibo.com/10986069
2楼2010-11-30 16:48:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (正式写手)
本帖仅楼主可见
赞一下
3楼2010-11-30 16:53:07
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

gaogh

银虫 (初入文坛)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
youzhizhe(金币+2):鼓励交流讨论。 2010-11-30 17:01:26
好多列用不同的y就可以了啊，不仅是能带计算，带自旋的态密度计算也同样处理
赞一下(3人) 回复此楼
4楼2010-11-30 16:59:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

za5010

铁虫 (初入文坛)

一个处理能带数据的代码

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2):谢谢分享，呵呵 2010-12-01 00:48:02
一个处理能带数据的代码，先编译一下生成一个可执行程序，然后在你的vasp计算目录下面运行生成的可执行程序即可得到一个可以用origin处理的能带数据

如果没有考虑自旋转，生成的文件名是bnd.dat
考虑了自选，则生成两个文件bndup.dat和bnddn.dat，一个自旋向上，另一个是自旋向下

程序执行过程中会要你输入费米能级（单位是eV）。

bnd.dat文件中，第一列是x坐标（没有意义，画图的时候要把它的标度隐掉，并写上高对称点），它按KPOINTS中的高对点顺序依次排开。

程序从网上找的，自己相应改一下.......版权不属于我......呵呵


program prog
real, allocatable :: e(:, ,eup(:,,edn(:,
real, allocatable :: k(:, ,dk( ,x(
real, dimension(3) ::k0,a
real, dimension(6) ::xxxx
character(len=32):: xx, yy ,nosence



write(6,*) 'fermi level (eV)'
read(5,*) ef

open(10,file='EIGENVAL', status='old')


read(10,*) iii, iii, iii, ispin
read(10,*) (xxxx(i),i=1,5)
read(10,*) xxxx(6)
read(10,*) xx
read(10,*) yy
read(10,*) nn,nk,nbands

allocate(e(nk,nbands))
allocate(eup(nk,nbands))
allocate(edn(nk,nbands))
allocate(k(nk,3))
allocate(dk(0:nk))
allocate(x(1:nk))


wtk=0.2500000E-01

if(ispin.eq.2) then
      open(13,file='bndup.dat')
          open(14,file='bnddn.dat')
          do i=1,nk
           read(10,*)
           read(10,*) (k(i,j),j=1,3),wtk
          do j=1,nbands
               read(10,*) jj,eup(i,j),edn(i,j)
           enddo                                                          
      enddo
      
          do n=1,nk
      dk(0)=0
          if (n.eq.1) then
        k0=k(n,
      endif
      a=k(n,-k0
      dk(n)=dk(n-1)+sqrt(dot_product(a,a))
      x(n)=dk(n)     
      k0=k(n,
      enddo

      do i=1,nk
        write(13,9030) x(i), ((eup(i,j)-ef),j=1,nbands)
        write(14,9030) x(i), ((edn(i,j)-ef),j=1,nbands)
      enddo

else
      open(11,file='bnd.dat')
          do i=1,nk
           read(10,*)
           read(10,*) nosence,nosence,nosence,(k(i,j),j=1,3)
           read(10,*)
           do j=1,nbands
               read(10,*) jj,e(i,j)
           enddo
      enddo
          
          do n=1,nk
      dk(0)=0
          if (n.eq.1) then
        k0=k(n,
      endif
      a=k(n,-k0
      dk(n)=dk(n-1)+sqrt(dot_product(a,a))
      x(n)=dk(n)     
      k0=k(n,
      enddo

      do i=1,nk
        write(11,9030) x(i),((e(i,j)-ef),j=1,nbands)
      enddo
                              
endif
9030 format (1x,f9.4,2x ,2000(f12.6,2x))


stop
end program prog
赞一下(2人) 回复此楼
5楼2010-11-30 20:39:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liangpei
a)
6楼2012-05-18 13:08:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mxyxmc

金虫 (初入文坛)
学习一下～
回复此楼
7楼2012-10-26 15:26:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lu_wei07

银虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
421680楼: Originally posted by za5010 at 2010-11-30 20:39:25
一个处理能带数据的代码，先编译一下生成一个可执行程序，然后在你的vasp计算目录下面运行生成的可执行程序即可得到一个可以用origin处理的能带数据

如果没有考虑自旋转，生成的文件名是bnd.dat
考虑了自选，则 ...

你好，有这个程序倒出来的文件0KB 啊 ，谢谢了
赞一下 回复此楼
这只是一场游戏
8楼2012-11-14 13:17:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阳光的影子

新虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
你好，我想请教你用VASP计算能带时手动输入line模式的KPOINTS文件怎么弄？是根据自己体系的晶系类型找到高对称点坐标，然后根据文献中的路径去自己编辑的吗？这样做对不对？还是从MS里找？
赞一下 回复此楼
9楼2013-12-04 18:02:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 晨曦3999 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0832食品科学与工程学硕282调剂 +6 鱼在水中游a 2026-04-02 9/450 2026-04-05 11:45 by flysky1234
[考研] 270求调剂 +9 小杰pp 2026-03-31 11/550 2026-04-05 11:02 by 风雨无晴
[考研] 358求调剂 +6 秋gk 2026-04-04 6/300 2026-04-05 11:01 by xiayan13521
[考研] 工科求调剂 +14 11ggg 2026-04-03 14/700 2026-04-05 10:58 by 风雨无晴
[考研] 材料工程302分求调剂 +6 zyx上岸！ 2026-04-04 6/300 2026-04-05 07:57 by qlm5820
[考研] 323分（计算机视觉和大模型项目）能直接上手 +3 chaoxiicy 2026-04-01 3/150 2026-04-05 00:50 by chongya
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 22:08 by 啵啵啵0119
[考研] 材料调剂 +10 懒羊羊轻置玉臀 2026-04-02 11/550 2026-04-04 21:56 by laoshidan
[考研] 278求调剂 +3 依旧！ 2026-04-02 4/200 2026-04-04 20:27 by 蓝云思雨
[考研] 348分环境工程·调剂 +10 吴彦祖24k 2026-04-03 11/550 2026-04-04 14:19 by 无际的草原
[考研] 一志愿C9的化学工程（085602） 340分，感觉校内调剂无望，求调剂 +9 万事宜臻 2026-04-04 9/450 2026-04-04 11:49 by 啵啵啵0119
[考研] 085600材料与化工调剂 +26 kikiki7 2026-03-30 27/1350 2026-04-04 09:18 by qlm5820
[考研] 335求调剂 +7 沈清璃 2026-04-03 7/350 2026-04-03 18:55 by lijunpoly
[考研] 求调剂 +9 akdhjs 2026-03-31 11/550 2026-04-03 13:32 by akdhjs
[考研] 316求调剂 +14 舟自梗 2026-04-01 18/900 2026-04-03 10:28 by linyelide
[考研] 一志愿西安交大材料学硕（英一数二）347，求调剂到高分子/材料相关专业 +7 zju51 2026-03-31 9/450 2026-04-01 19:35 by CFQZAFU
[考研] 262求调剂 +9 励志一定发文章 2026-03-31 10/500 2026-04-01 12:22 by sunshine0013
[考研] 08工科，295，接受跨专业调剂 +6 lmnlzy 2026-03-31 6/300 2026-04-01 11:02 by 逆水乘风
[考研] 326求调剂 +4 崽崽仔 2026-03-31 4/200 2026-04-01 09:58 by 我的船我的海
[考研] 370求调剂 +3 080700调剂 2026-03-30 3/150 2026-03-31 01:09 by A_Zhe
信息提示
关闭
请填处理意见
关闭