0001
0002
0003
0004
0005 global noCheckSum;
0006 noCheckSum=true;
0007
0008
0009
0010 lambda=418;
0011
0012
0013 nNmax=50;
0014
0015 nNmaxESA=500;
0016
0017
0018 a=25;
0019
0020
0021 epsilonM=1.77;
0022
0023
0024
0025 epsilonIn=epsAg(lambda);
0026
0027 dArray=[0.1 0.2 0.25:0.25:15];
0028 nNbD=length(dArray);
0029 etaRadPerp=zeros(1,nNbD);
0030 etaRadPara=zeros(1,nNbD);
0031 stPtArray=[];
0032
0033
0034 for indd=1:nNbD
0035 d=dArray(indd);
0036 stMdip=DipSolveSingleSphere(nNmax,a,lambda,epsilonM,epsilonIn,d,nNmaxESA,'noplot','nocoeff');
0037
0038 etaRadPerp(indd)=stMdip.MRadPerp./stMdip.MTotPerp;
0039 etaRadPara(indd)=stMdip.MRadPara./stMdip.MTotPara;
0040
0041 stPtArray=[stPtArray, struct('x',a+d,'y',0,'z',0), struct('x',0,'y',0,'z',a+d)];
0042 end
0043
0044 stM=PweSolveSingleSphere(nNmax,a,lambda,epsilonM,epsilonIn,'noplot','coeff');
0045
0046 stEpts=PweEpts(lambda,epsilonM,stM,stPtArray,'Outside');
0047
0048 MLocPerpA=zeros(1,nNbD);
0049 MLocParaC=zeros(1,nNbD);
0050 for indd=1:nNbD
0051 MLocPerpA(indd)=abs(stEpts(2*indd-1).Er).^2;
0052 MLocParaC(indd)=abs(stEpts(2*indd).Etheta).^2+abs(stEpts(2*indd).Ephi).^2;
0053 end
0054 MFluoPerpA=MLocPerpA.*etaRadPerp;
0055 MFluoParaC=MLocParaC.*etaRadPara;
0056
0057
0058
0059 lambda=transpose(300:1:650);
0060 epsilonIn=epsAg(lambda);
0061
0062 stMdipD1=DipSolveSingleSphere(nNmax,a,lambda,epsilonM,epsilonIn,1,nNmaxESA,'noplot','nocoeff');
0063 stMdipD2=DipSolveSingleSphere(nNmax,a,lambda,epsilonM,epsilonIn,2,nNmaxESA,'noplot','nocoeff');
0064 stMdipD5=DipSolveSingleSphere(nNmax,a,lambda,epsilonM,epsilonIn,5,nNmaxESA,'noplot','nocoeff');
0065
0066
0067 stM=PweSolveSingleSphere(nNmax,a,lambda,epsilonM,epsilonIn,'noplot','coeff');
0068
0069
0070 nNbTheta=721;
0071 stEsurf=PweSurfProperties(stM,a+1,nNbTheta);
0072 F0E4PerpAveD1=stEsurf.F0E4PerpAve;
0073 F0E4ParaAveD1=stEsurf.F0E4ParaAve;
0074 stEsurf=PweSurfProperties(stM,a+2,nNbTheta);
0075 F0E4PerpAveD2=stEsurf.F0E4PerpAve;
0076 F0E4ParaAveD2=stEsurf.F0E4ParaAve;
0077 stEsurf=PweSurfProperties(stM,a+5,nNbTheta);
0078 F0E4PerpAveD5=stEsurf.F0E4PerpAve;
0079 F0E4ParaAveD5=stEsurf.F0E4ParaAve;
0080 clear stEsurf stM;
0081
0082
0083 MFluoPerpD1=F0E4PerpAveD1./stMdipD1.MTotPerp;
0084 MFluoParaD1=F0E4ParaAveD1./stMdipD1.MTotPara;
0085 MFluoPerpD2=F0E4PerpAveD2./stMdipD2.MTotPerp;
0086 MFluoParaD2=F0E4ParaAveD2./stMdipD2.MTotPara;
0087 MFluoPerpD5=F0E4PerpAveD5./stMdipD5.MTotPerp;
0088 MFluoParaD5=F0E4ParaAveD5./stMdipD5.MTotPara;
0089
0090
0091
0092
0093
0094
0095 scrsz = get(0,'ScreenSize');
0096 figAspectRatio=2.5;
0097 figWidth=7/8*scrsz(3);
0098 figure('Name','Fig. 6.14', ...
0099 'Position',[(scrsz(3)-figWidth)/2 scrsz(4)-150-figWidth/figAspectRatio figWidth figWidth/figAspectRatio]);
0100
0101
0102 axtmp=subplot(1,2,1);
0103 postmp=get(axtmp,'Position');
0104 delete(axtmp);
0105
0106 ax1=axes('Position',postmp,...
0107 'XAxisLocation','bottom','Xcolor','k', ...
0108 'YAxisLocation','left','Ycolor','b', ...
0109 'XLim',[-0.5 14.9],'YLim',[1 1000], ...
0110 'YScale','log');
0111 set(get(ax1,'XLabel'),'String','Distance from surface, d [nm]');
0112 set(get(ax1,'YLabel'),'String','Enhancement Factor');
0113
0114 line(dArray,MLocPerpA,'Color','c','LineStyle','-');
0115 line(dArray,MFluoPerpA,'Color','b','LineStyle','-');
0116 line(dArray,MLocParaC,'Color','c','LineStyle',':');
0117 line(dArray,MFluoParaC,'Color','b','LineStyle',':');
0118
0119 legend({'M_{Loc}^{Perp}(A)','M_{Fluo}^{Perp}(A)', ...
0120 'M_{Loc}^{Para}(C)','M_{Fluo}^{Para}(C)'},'Location','SouthEast');
0121
0122 title('EFs and modified quatum yield for a silver sphere (a=25nm) in water');
0123
0124 axes('Position',get(ax1,'Position'),...
0125 'XAxisLocation','top','XColor','k', ...
0126 'YAxisLocation','right','Ycolor','r', ...
0127 'Color','none', ...
0128 'YScale','log', ...
0129 'YLim',[2e-4 1],'XLim',[-0.5 14.9], ...
0130 'XTickLabel',[]);
0131 line(dArray,etaRadPerp,'Color','r','LineStyle','-');
0132 line(dArray,etaRadPara,'Color','r','LineStyle',':');
0133
0134 legend({'\eta_{Rad}^{Perp}','\eta_{Rad}^{Para}'},'Location','NorthEast','Color','w');
0135
0136
0137 subplot(1,2,2);
0138 semilogy(lambda,MFluoPerpD1,'b-', lambda,MFluoParaD1,'b:', ...
0139 lambda,MFluoPerpD2,'g-', lambda,MFluoParaD2,'g:', ...
0140 lambda,MFluoPerpD5,'r-', lambda,MFluoParaD5,'r:');
0141 axis([300 650 0.02 25]);
0142 xlabel('Wavelength [nm]');
0143 ylabel('Average Fluorescence EF');
0144 legend({'d=1nm - Perp','d=1nm - Para', ...
0145 'd=2nm - Perp','d=2nm - Para', ...
0146 'd=5nm - Perp','d=5nm - Para'}, ...
0147 'Location','NorthEast');
0148 title('Average Fluorescence EFs for a silver sphere (a=25nm) in water');
0149