0001
0002
0003
0004
0005 global noCheckSum;
0006 noCheckSum=true;
0007
0008
0009
0010 nNmax=50;
0011
0012 nNmaxESA=500;
0013
0014
0015 lambda=transpose(250:1:650);
0016
0017
0018 epsilonM=1.77;
0019
0020
0021
0022 epsilonIn=epsAg(lambda);
0023
0024
0025
0026 stMdipA10d1=DipSolveSingleSphere(nNmax,10,lambda,epsilonM,epsilonIn,1,nNmaxESA,'noplot','nocoeff');
0027 stMdipA25d1=DipSolveSingleSphere(nNmax,25,lambda,epsilonM,epsilonIn,1,nNmaxESA,'noplot','nocoeff');
0028 stMdipA25d2=DipSolveSingleSphere(nNmax,25,lambda,epsilonM,epsilonIn,2,nNmaxESA,'noplot','nocoeff');
0029 stMdipA25d5=DipSolveSingleSphere(nNmax,25,lambda,epsilonM,epsilonIn,5,nNmaxESA,'noplot','nocoeff');
0030 stMdipA50d1=DipSolveSingleSphere(nNmax,50,lambda,epsilonM,epsilonIn,1,nNmaxESA,'noplot','nocoeff');
0031
0032
0033
0034
0035
0036
0037
0038 scrsz = get(0,'ScreenSize');
0039 figAspectRatio=1;
0040 figHeight=scrsz(4)-150;
0041 figWidth=figAspectRatio*figHeight;
0042 figure('Name','Fig. 6.13', ...
0043 'Position',[(scrsz(3)-figWidth)/2 50 figWidth figHeight]);
0044
0045
0046 subplot(2,2,1);
0047 semilogy(lambda,stMdipA10d1.MRadPerp./stMdipA10d1.MTotPerp, '-', ...
0048 lambda,stMdipA25d1.MRadPerp./stMdipA25d1.MTotPerp, '-', ...
0049 lambda,stMdipA50d1.MRadPerp./stMdipA50d1.MTotPerp, '-');
0050 axis([250 650 2e-4 0.2]);
0051 title('Modified quantum yield \eta_{Rad}^{Perp} for d=1 nm');
0052 legend({'a=10nm','a=25nm','a=50nm'},'Location','NorthEast');
0053 xlabel('Wavelength [nm]');
0054 ylabel('\eta_{Rad}^{Perp}');
0055
0056
0057 subplot(2,2,3);
0058 h=semilogy(lambda,stMdipA10d1.MTotPerp, '-', ...
0059 lambda,stMdipA25d1.MTotPerp, '-', ...
0060 lambda,stMdipA50d1.MTotPerp, '-');
0061 axis([250 650 0.1 1e6]);
0062 title('Decay rate EFs for d=1 nm (Perp dipole)');
0063 xlabel('Wavelength [nm]');
0064 ylabel('Decay Rate EFs');
0065 legend({'a=10nm','a=25nm','a=50nm'},'Location','NorthEast');
0066 line(lambda,stMdipA10d1.MRadPerp, 'LineStyle','-','Color',get(h(1),'Color'));
0067 line(lambda,stMdipA25d1.MRadPerp, 'LineStyle','-','Color',get(h(2),'Color'));
0068 line(lambda,stMdipA50d1.MRadPerp, 'LineStyle','-','Color',get(h(3),'Color'));
0069 text(280,1e5,0,'M_{Tot}^{Perp}');
0070 text(280,10,0,'M_{Rad}^{Perp}');
0071
0072
0073 subplot(2,2,2);
0074 h=semilogy(lambda,stMdipA25d1.MRadPerp./stMdipA25d1.MTotPerp, '-', ...
0075 lambda,stMdipA25d2.MRadPerp./stMdipA25d2.MTotPerp, '-', ...
0076 lambda,stMdipA25d5.MRadPerp./stMdipA25d5.MTotPerp, '-');
0077 line(lambda,stMdipA25d1.MRadPara./stMdipA25d1.MTotPara, 'LineStyle','--','Color',get(h(1),'Color'));
0078 line(lambda,stMdipA25d2.MRadPara./stMdipA25d2.MTotPara, 'LineStyle','--','Color',get(h(2),'Color'));
0079 line(lambda,stMdipA25d5.MRadPara./stMdipA25d5.MTotPara, 'LineStyle','--','Color',get(h(3),'Color'));
0080 axis([250 650 2e-4 0.99]);
0081 title('Modified quantum yield for a=25 nm');
0082 legend({'d=1nm - Perp','d=2nm - Perp','d=5nm - Perp', ...
0083 'd=1nm - Para','d=2nm - Para','d=5nm - Para'},'Location','South');
0084 xlabel('Wavelength [nm]');
0085 ylabel('\eta_{Rad}^{Perp}');
0086
0087
0088
0089
0090 subplot(2,2,4);
0091 h=semilogy(lambda,stMdipA25d1.MTotPerp, '-', ...
0092 lambda,stMdipA25d2.MTotPerp, '-');
0093 line(lambda,stMdipA25d1.MTotPara, 'LineStyle','--','Color',get(h(1),'Color'));
0094 line(lambda,stMdipA25d2.MTotPara, 'LineStyle','--','Color',get(h(2),'Color'));
0095 line(lambda,stMdipA25d1.MRadPerp, 'LineStyle','-','Color',get(h(1),'Color'));
0096 line(lambda,stMdipA25d2.MRadPerp, 'LineStyle','-','Color',get(h(2),'Color'));
0097 line(lambda,stMdipA25d1.MRadPara, 'LineStyle','--','Color',get(h(1),'Color'));
0098 line(lambda,stMdipA25d2.MRadPara, 'LineStyle','--','Color',get(h(2),'Color'));
0099 axis([250 650 0.1 1e6]);
0100 title('Decay rate EFs for a=25 nm');
0101 xlabel('Wavelength [nm]');
0102 ylabel('Decay Rate EFs');
0103 legend({'d=1nm - Perp','d=2nm - Perp','d=1nm - Para','d=2nm - Para'}, ...
0104 'Location','NorthEast');
0105 text(280,1e5,0,'M_{Tot}');
0106 text(280,10,0,'M_{Rad}');