Home > SPlaC v1_0 > BookFigures > MakeFigH_1.m

MakeFigH_1

PURPOSE ^

Function to reproduce Fig. H.1

SYNOPSIS ^

function MakeFigH_1()

DESCRIPTION ^

 Function to reproduce Fig. H.1
 Convergence of Mie series for decay rate EF of dipoles
 Note that the first few terms for M_{Rad}^{Perp} in the
 original figure in (a) and (b) were slightly wrong.

 This file is part of the SPlaC v1.0 package (copyright 2008)
 Check the README file for further information

EXAMPLE OF OUTPUT ^

Example figure output

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 %%%%%%%%%%%%%%%%%%
0002 % Computations
0003 %%%%%%%%%%%%%%%%%%
0004 
0005 global noCheckSum;
0006 noCheckSum=true; % no check for numerical problems (much faster)
0007 % noCheckSum=false; % to check for numerical problems with Mie
0008 
0009 % N_{max}
0010 nNmax=50;
0011 nNmaxESA=500;
0012 
0013 % sphere radius in nm (scalar [1 x 1])
0014 a=50;
0015 
0016 % dielectric constant of outside medium (real positive scalar)
0017 epsilonM=1.77; % water
0018 
0019 % distance of dipole from surface (in nm)
0020 d=1;
0021 
0022 % solve decay rate EFs for perp and para dipoles
0023 % and return coeffs in the sum of Mie series
0024 % lambda=340, Mie only (nNMaxESA=nNmax);
0025 stMdip340Mie=DipSolveSingleSphere(nNmax,a,340,epsilonM,epsAg(340),d,nNmax,'noplot','coeffSum');
0026 % lambda=340, ESA (nNmax=0);
0027 stMdip340ESA=DipSolveSingleSphere(0,a,340,epsilonM,epsAg(340),d,nNmaxESA,'noplot','coeffSum');
0028 % lambda=340, Mie only (nNMaxESA=nNmax);
0029 stMdip394Mie=DipSolveSingleSphere(nNmax,a,394,epsilonM,epsAg(394),d,nNmax,'noplot','coeffSum');
0030 % lambda=340, ESA (nNmax=0);
0031 stMdip394ESA=DipSolveSingleSphere(0,a,394,epsilonM,epsAg(394),d,nNmaxESA,'noplot','coeffSum');
0032 
0033 % solve decay rate EFs for perp and para dipoles
0034 % for all wavelengths (lambda in nm)
0035 lambda=transpose(300:1:650);
0036 % wavelength-dependent dielectric function of sphere material
0037 % vector column, same size as lambda [1 x 1]
0038 epsilonIn=epsAg(lambda); % silver
0039 % nNmax = 50, no ESA (nNmaxESA=nNmax)
0040 stMdipN50=DipSolveSingleSphere(nNmax,a,lambda,epsilonM,epsilonIn,d,nNmax,'noplot','nocoeff');
0041 % nNmax = 50, with ESA up to nNmaxESA=500
0042 stMdipN500=DipSolveSingleSphere(nNmax,a,lambda,epsilonM,epsilonIn,d,nNmaxESA,'noplot','nocoeff');
0043 
0044 
0045 
0046 %%%%%%%%%%%%%%%%%%
0047 % Drawing
0048 %%%%%%%%%%%%%%%%%%
0049 
0050 % create, position, and resize figure
0051 scrsz = get(0,'ScreenSize'); % scrsz(4) contains screen height
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 % plot results for lambda=340 nm (Fig. H.1(a))
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 % plot results for lambda=394 nm (Fig. H.1(b))
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 % plot lambda dependence (Fig. H.1(c))
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

This web page is part of the SPlaC package © 2008. Contact: Eric Le Ru
Generated on Wed 03-Dec-2008 11:10:14 by m2html © 2003 (adapted)