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');