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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】能带图如何画出呢？
北京石油化工学院2026年研究生招生接收调剂公告
91/1返回列表
查看: 2868  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 307分材料专业求调剂 +5 Hll胡 2026-04-05 5/250 2026-04-05 15:37 by 今天星期八LQ
[考研] 材料与化工371求调剂 +10 陪琳看海 2026-04-04 11/550 2026-04-05 15:26 by luoyongfeng
[考研] 085410人工智能 初试316分 求调剂 +6 残星拂曙 2026-03-31 6/300 2026-04-05 12:15 by rainbow11
[考研] 专硕310求调剂 +5 捞捞我…. 2026-04-04 6/300 2026-04-04 23:33 by barlinike
[考研] 考研调剂 +5 四川王涛 2026-04-04 5/250 2026-04-04 22:18 by 啵啵啵0119
[考研] 085602 找调剂 +4 逆时针快乐 2026-04-02 4/200 2026-04-04 19:32 by 蓝云思雨
[考研] 材料383求调剂 +5 郭阳阳阳成 2026-04-04 5/250 2026-04-04 19:06 by dongzh2009
[考研] 求调剂：085600材料与化工，考材科基，总分319 +21 678lucky 2026-03-31 26/1300 2026-04-04 16:22 by dongzh2009
[考研] 320调剂 +4 农业工程与信息� 2026-04-03 4/200 2026-04-03 21:40 by lbsjt
[考研] 295求调剂 +3 尚偌呀 2026-04-03 4/200 2026-04-03 21:23 by zhq0425
[考研] 372分材料与化工（085600）一志愿湖南大学求调剂 +3 蓝笺片 2026-04-03 4/200 2026-04-03 17:58 by Jimmyandyou
[考研] 材料专硕322分 +13 哈哈哈吼吼吼哈 2026-04-01 13/650 2026-04-03 16:08 by 哦哦123
[考研] 机械专硕297 +3 Afksy 2026-04-03 3/150 2026-04-03 14:24 by 1753564080
[考研] 08工科，295，接受跨专业调剂 +8 lmnlzy 2026-03-30 8/400 2026-04-03 13:08 by nalakaiqi
[考研] 275求调剂 +14 waltzh 2026-04-01 14/700 2026-04-03 11:24 by qqq850810
[考博] 26年申博 +3 staryer 2026-03-30 4/200 2026-04-01 23:21 by ai4pharm
[考研] 0817化工学硕调剂 +11 努力上岸中！ 2026-03-31 11/550 2026-04-01 20:30 by 赖春艳
[考研] 0855机械初试280求调剂 +3 kazenotori 2026-03-31 3/150 2026-04-01 10:08 by JourneyLucky
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +5 双马尾痞老板2 2026-03-31 5/250 2026-04-01 09:04 by oooqiao
[考研] 调剂 +4 GK72 2026-03-30 4/200 2026-03-30 20:32 by dick_runner
信息提示
关闭
请填处理意见
关闭