0001
0002 scrsz = get(0,'ScreenSize');
0003 figAspectRatio=1;
0004 figHeight=scrsz(4)-150;
0005 figWidth=figAspectRatio*figHeight;
0006 figure('Name','Fig. 6.1', ...
0007 'Position',[(scrsz(3)-figWidth)/2 50 figWidth figHeight]);
0008
0009
0010
0011 epsilonM=1.0;
0012
0013
0014 for nn=1:2
0015
0016 if nn==1
0017
0018 lambda=transpose(200:2:800);
0019 epsilonIn=epsAg(lambda);
0020 sMetal='Ag';
0021 else
0022
0023 lambda=transpose(400:10:1000);
0024 epsilonIn=epsAu(lambda);
0025 sMetal='Au';
0026 end
0027
0028
0029 nNbStep=5000;
0030 disp(['Computations for ' sMetal]);
0031 stAllMd1=PlnSinglePlaneDecayRateEFs(epsilonM,lambda,epsilonIn,1,nNbStep);
0032 stAllMd5=PlnSinglePlaneDecayRateEFs(epsilonM,lambda,epsilonIn,5,nNbStep);
0033 stAllMd10=PlnSinglePlaneDecayRateEFs(epsilonM,lambda,epsilonIn,10,nNbStep);
0034 stAllMd20=PlnSinglePlaneDecayRateEFs(epsilonM,lambda,epsilonIn,20,nNbStep);
0035
0036
0037
0038 subplot(2,2,nn);
0039
0040 h=semilogy(lambda,stAllMd1.MTotPerp, '-', ...
0041 lambda,stAllMd5.MTotPerp, '-', ...
0042 lambda,stAllMd10.MTotPerp, '-');
0043 line(lambda,stAllMd1.MTotPerpESA,'LineStyle',':','Color',get(h(1),'Color'));
0044 line(lambda,stAllMd5.MTotPerpESA,'LineStyle',':','Color',get(h(2),'Color'));
0045 line(lambda,stAllMd10.MTotPerpESA,'LineStyle',':','Color',get(h(3),'Color'));
0046
0047 set(gca,'YLim',[1 1.2e6]);
0048 title(['Total M_{Tot} - Perp dipole in air on ' sMetal]);
0049 xlabel('Wavelength [nm]');
0050 ylabel('M_{Tot}');
0051 legend({'d=1nm','d=5mn','d=10nm','d=1nm (ESA)', ...
0052 'd=5mn (ESA)','d=10nm (ESA)'},'Location','NorthEast');
0053
0054
0055 subplot(2,2,2+nn);
0056
0057 h=semilogy(lambda,stAllMd1.MRadPerp, '-', ...
0058 lambda,stAllMd5.MRadPerp, '--', ...
0059 lambda,stAllMd10.MRadPerp, '-', ...
0060 lambda,stAllMd20.MRadPerp, '--');
0061 legend({'d=1nm','d=5mn','d=10nm','d=20nm'},'Location','South');
0062 line(lambda,stAllMd1.MRadPara,'LineStyle','-','Color',get(h(1),'Color'));
0063 line(lambda,stAllMd5.MRadPara,'LineStyle','--','Color',get(h(2),'Color'));
0064 line(lambda,stAllMd10.MRadPara,'LineStyle','-','Color',get(h(3),'Color'));
0065 line(lambda,stAllMd20.MRadPara,'LineStyle','--','Color',get(h(4),'Color'));
0066
0067 set(gca,'YLim',[0.01 4]);
0068 title(['Radiative M_{Rad} - dipole in air on ' sMetal]);
0069 xlabel('Wavelength [nm]');
0070 ylabel('M_{Rad}');
0071 text(600,0.08,'Para');
0072 text(600,2,'Perp');
0073 end