0001 global noCheckSum;
0002 noCheckSum=true;
0003
0004
0005
0006 nNmax=50;
0007
0008
0009 epsilonM=1.77;
0010
0011
0012 scrsz = get(0,'ScreenSize');
0013 figAspectRatio=1.2;
0014 figHeight=scrsz(4)-150;
0015 figWidth=figAspectRatio*figHeight;
0016 figure('Name','Fig. 6.11', ...
0017 'Position',[(scrsz(3)-figWidth)/2 50 figWidth figHeight]);
0018
0019
0020 aArray=[10,25,50,50];
0021 sTitle={'Silver (a=10 nm) in water', ...
0022 'Silver (a=25 nm) in water', ...
0023 'Silver (a=50 nm) in water', ...
0024 'Gold (a=50 nm) in water'};
0025 CYscale={[0.5 4000], [0.5 500], ...
0026 [0.5 300], [2 100]};
0027
0028 lambdamin=[300,300,300,300];
0029 lambdamax=[600,600,800,800];
0030
0031
0032
0033 for nn=1:4
0034
0035 lambda=transpose(lambdamin(nn):1:lambdamax(nn));
0036
0037 if (nn<=3)
0038
0039
0040 epsilonIn=epsAg(lambda);
0041 else
0042
0043
0044 epsilonIn=epsAu(lambda);
0045 end
0046
0047
0048 a=aArray(nn);
0049
0050
0051 stM=PweSolveSingleSphere(nNmax,a,lambda,epsilonM,epsilonIn,'noplot','coeff');
0052
0053
0054
0055 nNbTheta=721;
0056 stEsurf=PweEsurf(lambda,epsilonM,stM,a,nNbTheta,'Outside');
0057
0058
0059 indtpiov2=(nNbTheta+1)/2;
0060 MLocA=abs(stEsurf.Ecr(:,indtpiov2)).^2+ ...
0061 abs(stEsurf.Ect(:,indtpiov2)).^2;
0062
0063
0064 subplot(2,2,nn);
0065 hg=semilogy(lambda,MLocA,'--', ...
0066 lambda,stEsurf.MLocAve,'-', ...
0067 lambda,sqrt(stEsurf.F0E4Ave),'-');
0068 set(hg(3),'LineWidth',1.5);
0069 axis([[lambdamin(nn) lambdamax(nn)] CYscale{nn}]);
0070 title(sTitle{nn});
0071 legend({'M_{Loc}(A)','<M_{Loc}>','(<F^0_{E4}>)^{1/2}'},'Location','SouthEast');
0072 xlabel('Wavelength [nm]');
0073 ylabel('Enhancement Factor');
0074 end