Home > SPlaC v1_0 > BookFigures > MakeFig6_14.m

MakeFig6_14

PURPOSE ^

Function to reproduce Fig. 6.14

SYNOPSIS ^

function MakeFig6_14()

DESCRIPTION ^

 Function to reproduce Fig. 6.14
 (a) distance dependence of various EFs relevant to SEF near a silver
 sphere in water at the dipolar LSP resonance (418 nm)
 (b) wavelength dependence of the average approximate fluorescence
 EF for a perpendicular or parallel dipole at a distance d=1,2, and 5nm

 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 % wavelengths (column)
0010 lambda=418;
0011 
0012 % N_{max} for Mie series
0013 nNmax=50;
0014 % N_{max} for Mtot in the ESA (necessary fo good convergence of M_{Tot}
0015 nNmaxESA=500; 
0016 
0017 % sphere radius in nm
0018 a=25;
0019 
0020 % dielectric constant of outside medium (real positive)
0021 epsilonM=1.77; % water
0022 
0023 % wavelength-dependent dielectric function of sphere material
0024 % vector column, same size as lambda [L x 1]
0025 epsilonIn=epsAg(lambda); % silver
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 % Solve dipole emission problems for each distance using Mie theory
0033 % at lambda=418
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     % create point array for PWE calculations in the same loop
0041     stPtArray=[stPtArray, struct('x',a+d,'y',0,'z',0), struct('x',0,'y',0,'z',a+d)];
0042 end
0043 % Solve PWE excitation problem at lambda=418 mn
0044 stM=PweSolveSingleSphere(nNmax,a,lambda,epsilonM,epsilonIn,'noplot','coeff');
0045 % Calculate LFIEF at all relevant points
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 % Solve dipole emission problem for all wavelength for d=1, 2, and 5 nm
0059 lambda=transpose(300:1:650); 
0060 epsilonIn=epsAg(lambda); % silver
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 % Solve PWE problem for all wavelength
0067 stM=PweSolveSingleSphere(nNmax,a,lambda,epsilonM,epsilonIn,'noplot','coeff');
0068 % computes surface averaged properties (in particular <F>)
0069 % at various distances
0070 nNbTheta=721; % nb of theta for averages
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 % computes average Fluo EF for relevant cases
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 % Drawing
0092 %%%%%%%%%%%%%%%%%
0093 
0094 % create, position, and resize figure
0095 scrsz = get(0,'ScreenSize'); % scrsz(3) contains screen width
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 % plot distance dependence
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 % plot lambda dependence
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

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)