0001
0002
0003
0004
0005 global noCheckSum;
0006 noCheckSum=true;
0007
0008
0009
0010 nNmax=50;
0011 nNmaxESA=500;
0012
0013
0014 a=50;
0015
0016
0017 epsilonM=1.77;
0018
0019
0020 d=1;
0021
0022
0023
0024
0025 stMdip340Mie=DipSolveSingleSphere(nNmax,a,340,epsilonM,epsAg(340),d,nNmax,'noplot','coeffSum');
0026
0027 stMdip340ESA=DipSolveSingleSphere(0,a,340,epsilonM,epsAg(340),d,nNmaxESA,'noplot','coeffSum');
0028
0029 stMdip394Mie=DipSolveSingleSphere(nNmax,a,394,epsilonM,epsAg(394),d,nNmax,'noplot','coeffSum');
0030
0031 stMdip394ESA=DipSolveSingleSphere(0,a,394,epsilonM,epsAg(394),d,nNmaxESA,'noplot','coeffSum');
0032
0033
0034
0035 lambda=transpose(300:1:650);
0036
0037
0038 epsilonIn=epsAg(lambda);
0039
0040 stMdipN50=DipSolveSingleSphere(nNmax,a,lambda,epsilonM,epsilonIn,d,nNmax,'noplot','nocoeff');
0041
0042 stMdipN500=DipSolveSingleSphere(nNmax,a,lambda,epsilonM,epsilonIn,d,nNmaxESA,'noplot','nocoeff');
0043
0044
0045
0046
0047
0048
0049
0050
0051 scrsz = get(0,'ScreenSize');
0052 figAspectRatio=1.3;
0053 figHeight=scrsz(4)-150;
0054 figWidth=figAspectRatio*figHeight;
0055 figure('Name','Fig. H.1', ...
0056 'Position',[(scrsz(3)-figWidth)/2 50 figWidth figHeight]);
0057
0058
0059
0060 subplot(2,2,1);
0061 loglog(1:nNmax,stMdip340Mie.MRadPerpCoeffs,'c-',1:nNmax,stMdip340Mie.MRadParaCoeffs,'m-');
0062 axis([0.8 550 1e-3 9900]);
0063 xlabel('Angular momentum order l');
0064 ylabel('Series Term in EF');
0065 line(1:nNmax,stMdip340Mie.MTotPerpCoeffs,'LineStyle','-','Color','b');
0066 line(1:nNmax,stMdip340Mie.MTotParaCoeffs,'LineStyle','-','Color','r');
0067 line(1:nNmaxESA,stMdip340ESA.MTotPerpCoeffs,'LineStyle',':','Color','b');
0068 line(1:nNmaxESA,stMdip340ESA.MTotParaCoeffs,'LineStyle',':','Color','r');
0069 legend({'M_{Rad}^{Perp} (Mie)','M_{Rad}^{Para} (Mie)', ...
0070 'M_{Tot}^{Perp} (Mie)','M_{Tot}^{Para} (Mie)', ...
0071 'M_{Tot}^{Perp} (ESA)','M_{Tot}^{Para} (ESA)'},'Location','South');
0072 title('(a) Silver sphere in water, a=50nm, d=1nm, \lambda=340nm');
0073
0074
0075 subplot(2,2,3);
0076 loglog(1:nNmax,stMdip394Mie.MRadPerpCoeffs,'c-',1:nNmax,stMdip394Mie.MRadParaCoeffs,'m-');
0077 axis([0.8 550 1e-3 600]);
0078 xlabel('Angular momentum order l');
0079 ylabel('Series Term in EF');
0080 line(1:nNmax,stMdip394Mie.MTotPerpCoeffs,'LineStyle','-','Color','b');
0081 line(1:nNmax,stMdip394Mie.MTotParaCoeffs,'LineStyle','-','Color','r');
0082 line(1:nNmaxESA,stMdip394ESA.MTotPerpCoeffs,'LineStyle',':','Color','b');
0083 line(1:nNmaxESA,stMdip394ESA.MTotParaCoeffs,'LineStyle',':','Color','r');
0084 legend({'M_{Rad}^{Perp} (Mie)','M_{Rad}^{Para} (Mie)', ...
0085 'M_{Tot}^{Perp} (Mie)','M_{Tot}^{Para} (Mie)', ...
0086 'M_{Tot}^{Perp} (ESA)','M_{Tot}^{Para} (ESA)'},'Location','South');
0087 title('(b) Silver sphere in water, a=50nm, d=1nm, \lambda=394nm');
0088
0089
0090 axtmp=subplot(2,2,[2,4]);
0091 postmp=get(axtmp,'Position');
0092 delete(axtmp);
0093 ax1=axes('Position',postmp,...
0094 'XAxisLocation','bottom','XColor','k', ...
0095 'YAxisLocation','left','Ycolor','b', ...
0096 'YScale','Log', ...
0097 'XLim', [300 650], 'YLim', [1e-3 1e6]);
0098 set(get(ax1,'XLabel'),'String','Wavelength [nm]');
0099 set(get(ax1,'YLabel'),'String','Decay rate EF');
0100
0101 line(lambda,stMdipN500.MTotPerp,'Color','b','LineStyle','-');
0102 line(lambda,stMdipN50.MTotPerp,'Color','b','LineStyle',':');
0103 line(lambda,stMdipN50.MRadPerp,'Color','c','LineStyle','-');
0104 legend({'M_{Tot}^{Perp} (L_{Max}=500)','M_{Tot}^{Perp} (L_{Max}=50)', ...
0105 'M_{Rad}^{Perp}'},'Location','NorthEast');
0106 title('(c) Silver sphere in water, a=50nm, d=1nm');
0107
0108 ax2=axes('Position',get(ax1,'Position'),...
0109 'XAxisLocation','top','XColor','k', ...
0110 'XTickLabel',[], ...
0111 'YAxisLocation','right','Ycolor',[0,0.5,0], ...
0112 'Color','none', ...
0113 'YScale','Log', ...
0114 'YLim',[1e-3 10],'XLim',[300 650]);
0115 set(get(ax2,'Ylabel'),'String','Modified Quantum Yield');
0116
0117 line(lambda,stMdipN50.MRadPerp./stMdipN500.MTotPerp,'Color',[0,0.5,0],'LineStyle','-');
0118 line(lambda,stMdipN50.MRadPerp./stMdipN50.MTotPerp,'Color',[0,0.5,0],'LineStyle',':');
0119 legend({'\eta_{Rad}^{Perp} (L_{Max}=500)','\eta_{Rad}^{Perp} (L_{Max}=50)'}, ...
0120 'Location','SouthEast','Color','w');
0121
0122
0123
0124
0125