| 查看: 5447 | 回复: 6 | ||
| 【奖励】 本帖被评价3次，作者byin增加金币 2 个 | ||
| 本帖产生 1 个 1ST强帖 ，点击这里进行查看 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[资源]
用matlab脚本绘制vasp中的能带图
|
||
|
把EIGENVAL和.m脚本放一起，修改脚本里的一些参数“手动修改区”，然后运行脚本即可。见笑了。 %%%%%%%%%%%%%%%%%%%% 此脚本用于VASP的能带绘制 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clc clear %%%%%%%%%%%%%%%%%%%%%%%%%%% 手动修改区 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% EFermi=6.1748; %设置费米能级 NBands=8; %一共计算的能带数目 inter=100; %每段PATH的插值点数目 NPath=4; %PATH数目 Ktotal=inter*NPath; %特殊K点的总数目 xmin=1; %x，y坐标轴的范围 xmax=Ktotal; ymin=-6; ymax=7; kpath={'L','\Gamma' ,'X','K','\Gamma'}; %设置K空间的路径 %%%%%%%%%%%%%%%%%%%%%%%%%% 图像边框 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% A4 = figure('PaperType','a4letter','PaperSize',[20.98 29.68]); axes1 = axes('Parent',A4,'Position',[0.2 0.1 0.6 0.8],'LineWidth',1.5,'FontSize',15.0,'FontWeight','bold'); %Position [left bottom width height] box(axes1,'on'); hold(axes1,'all'); %%%%%%%%%%%%%%%%%%%%%%%%%%% 数据读取 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% data=textread('EIGENVAL','','headerlines',7); for i=1:1:NBands, for j=1:1:Ktotal, P(j,: )=[j,data((NBands+1)*(j-1)+1+i,2)-EFermi]; %读出每条能带的数据（每个i中的各个j），做出曲线 end plot(P(:,1),P(:,2),'-r.'); hold on end %%%%%%%%%%%%%%%%%%%%%%%%%%%% 坐标轴设置 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% axis([xmin xmax ymin ymax]); set(gca,'XTick',[1,inter:inter:Ktotal]) set(gca,'XTickLabel','') j=1; %高对称点符号 for i=[1,inter:inter:Ktotal] text(i,ymin-0.7,kpath(j),'FontSize',18, 'horizontalAlignment', 'center') j=j+1; end ylabel('E-E_F (eV)'); X=[xmin xmax]; %辅助线 - 水平 Y=[0 0]; line(X,Y,'LineWidth',1.5,'Color',[0 0 0]) for i=1:1: (NPath-1) %辅助线 - 垂直 X=[i*inter i*inter]; Y=[ymin ymax]; line(X,Y,'LineWidth',1.5,'Color',[0 0 0]); end print -dpng Band_Structure_of_ %输出图像文件 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% END %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% [ Last edited by byin on 2011-5-30 at 21:27 ] |
回复此楼» 猜你喜欢
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有7人回复
-
面上意见出来了
已经有23人回复
-
能否申诉？
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何用matlab画出石墨烯高对称点的能带结构
已经有3人回复
- 在VASP中的脚本
已经有4人回复
- vasp 算能带k点数目问题
已经有10人回复
- 推荐几个vasp处理计算结果的python脚本
已经有66人回复
- 计算模拟需要的服务器 硬件 MS,VASP,matlab等等吧
已经有7人回复
- NBANDS 的设置问题
已经有6人回复
- 电机原理与设计的MATLAB分析
已经有140人回复
- 【求助】请问如何用P4VASP在LINUX系统下画能带图，或者导出数据
已经有5人回复
- vasp_bnd的Matlab改进版，更方便~
已经有3人回复
- 请教VASP CI-NBE分析脚本的使用问题
已经有6人回复
- matlab多维插值原理是什么？
已经有5人回复
- VASP有几个脚本文件？
已经有8人回复
- DOSCAR画态密度
已经有16人回复
- VASP能带横坐标怎么画
已经有6人回复
- 【活动or资源】vasp晶体结构优化
已经有142人回复
- 【求助】vasp作业提交脚本
已经有5人回复
- 【求助】求将vasp的输出文件转化成xyz坐标的脚本
已经有11人回复
- 【求助】Matlab 四维图形的绘制。
已经有12人回复
重新试过了，可以了，但图不好看，还需要不断改进才行。谢谢 
4楼2011-06-03 11:57:02
2楼2011-06-01 19:52:25
3楼2011-06-01 21:15:17
5楼2011-11-05 21:46:17
30