24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

xianggui7895

银虫 (小有名气)

1ST强帖: 1
应助: 0 (幼儿园)
金币: 442.1
散金: 44
红花: 4
帖子: 190
在线: 53.4小时
虫号: 980738
注册: 2010-03-24
性别: GG
专业: 无机非金属类光电信息与功

[交流] 【原创】利用Matlab绘制Materials studio能带图已有7人参与

Materials Studio提供的能带图既呆板又难看，能不能把能带图画漂亮点呢？

答案当然是肯定的。MS提供了将能带图（实际上是所有chart图形）导出为csv文件的功能。能带图的CSV文件是逗号分割的数据位置，一共两列，第一列是K点位置坐标，第二列是各K点相应的能量值。傻瓜点的做法的是导入到Origin中重新作图，也有专门的软件对数据进行处理画图。不过，我在这里介绍使用Matlab对能带数据进行处理和绘图的方法，感兴趣的大家可以尝试一下，matlab强大的画图功能不会让你失望的：）

第一步：在MS中将能带图导出（file->export...）为csv文件；

第二步：打开Matlab，将csv文件保存到工作目录下；

第三步：新建m文件，内容如下，不同能带图进行相应修改：）

CODE:

% plot band structure from *.csv file clear;  

% general information 

Filename = 'ZnO Band Structure-dot.csv'; 

Num_Pts = 67;   % number of Band structure K points 

% Note: MS exported csv files of band structure have two columns 

%   first is K points positions, and second is corresponding energy 

%   actual points number is Num_Pts-1:)  

% computing process 

M = csvread(Filename); 

x = M(:,1); y = M(:,2); 

L = size(M,1); 

x = reshape(x,Num_Pts,L/Num_Pts); 

y = reshape(y,Num_Pts,L/Num_Pts); 

x(67,:) = []; y(67,:) = []; 

plot(x,y,'-+','LineWidth',2); 

ylabel('Energy(eV)'); 

title('Band Structure of ZnO'); 

line([0 1],[0 0],'Color','k','LineWidth',2,'LineStyle','--'); %Fermi energy line  

% special K-points marker id_kpts = [1 7 21 27 41 53 59 66];

label_kpts = 'G|A|H|K|G|M|L|H'; loc_kpts = x(id_kpts); set(gca,'XTick',loc_kpts,'XTickLabel',label_kpts); line([loc_kpts;loc_kpts],repmat(get(gca,'yLim')',1,size(id_kpts,2)),'Color','k');