pictureofgap.m

通过计算晶体结构来确定其光谱分布

%-------------------------draw the figure of the gap----------------------
%----------------------------OH-------------------------------------------
xti=[1,30,60,90];
xtilabel={'Γ','X','M','Γ'};
hN=hN+1;
foh=figure(hN);
lw=dth;
plot(OH,'-');
xlabel('k');
ylabel('{\omega}a/2{\pi}c');
title('H-polarization');
set(gca,'xlim',[1,Nk]);
set(gca,'XTick',xti);
set(gca,'XTickLabel',xtilabel);
set(gcf,'Position',[221 123 558 498]);
set(gca,'Position',[0.13 0.2249 0.780394 0.7001]);
KK=text(5,-0.20,['带隙中心频率点：',num2str(wh0)]);
kk=text(5,-0.14,['带隙宽度:      ',num2str(dth),'     ','( ',num2str(ohgapmin),'--',num2str(ohgapmax),' )']);
set(KK,'color',[1 0 0],'fontsize',[12]);
set(kk,'color',[1 0 0],'fontsize',[12]);
hold on
ddd=ones(1,Nk);
dd=ddd';
d=dd.*wh0;
ffh=plot(d,'r-');
set(ffh,'LineWidth',lw/(0.06321/16));
text(30,0.05,'带隙        红色宽线','FontSize',8);
hold off
%----------------------------OE-------------------------------------------
hN=hN+1;
foe=figure(hN);
lw=dte;
plot(OE,':');
xlabel('k');
ylabel('{\omega}a/2{\pi}c');
title('E-polarization');
set(gca,'xlim',[1,Nk]);
set(gca,'XTick',xti);
set(gca,'XTickLabel',xtilabel);
set(gcf,'Position',[221 123 558 498]);
set(gca,'Position',[0.13 0.2249 0.780394 0.7001]);
KK=text(5,-0.20,['带隙中心频率点：',num2str(we0)]);
kk=text(5,-0.14,['带隙宽度:      ',num2str(dte),'     ','( ',num2str(oegapmin),'--',num2str(oegapmax),' )']);
set(KK,'color',[1 0 0],'fontsize',[12]);
set(kk,'color',[1 0 0],'fontsize',[12]);
hold on
ddd=ones(1,Nk);
dd=ddd';
d=dd.*we0;
ffe=plot(d,'r-');
set(ffe,'LineWidth',lw/(0.06321/16));
text(30,0.05,'带隙        红色宽线','FontSize',8);
hold off
%----------------------------OHE-------------------------------------------
hN=hN+1;
fohe=figure(hN);
lw=dthe;
plot(OH,'-');
hold on
plot(OE,':');
xlabel('k');
ylabel('{\omega}a/2{\pi}c');
title('HE-polarization');
set(gca,'xlim',[1,Nk]);
set(gca,'XTick',xti);
set(gca,'XTickLabel',xtilabel);
set(gcf,'Position',[221 123 558 498]);
set(gca,'Position',[0.13 0.2249 0.780394 0.7001]);
kk=text(5,-0.20,['带隙中心频率点：',num2str(whe0)]);
KK=text(5,-0.14,['带隙宽度:      ',num2str(dthe),'     ','( ',num2str(ohegapmin),'--',num2str(ohegapmax),' )']);
set(KK,'color',[1 0 0],'fontsize',[12]);
set(kk,'color',[1 0 0],'fontsize',[12]);

hold on
ddd=ones(1,Nk);
dd=ddd';
d=dd.*whe0;
ffhe=plot(d,'r-');
set(ffhe,'LineWidth',lw/(0.06321/16));
text(30,0.15,'H-polarization  实线','FontSize',8);
text(30,0.10,'E-polarization  虚线','FontSize',8);
text(30,0.05,'带隙        红色宽线','FontSize',8);
hold off