Home > SPlaC v1_0 > BookFigures > MakeFig6_13.m

MakeFig6_13

PURPOSE ^

Function to reproduce Fig. 6.13

SYNOPSIS ^

function MakeFig6_13()

DESCRIPTION ^

 Function to reproduce Fig. 6.13
 (a) Modification of dipole emission for a perpendicular
 dipole at a distance d=1nm from a silver sphere
 of radius a=10, 25, or 50 nm in water.
 (b)Same for a=25 nm, for both perpendicular and parallel dipoles
 at a distance d=1, 2, or 5 nm.

 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} for Mie series
0010 nNmax=50; 
0011 % N_{max} for Mtot in the ESA (necessary fo good convergence of M_{Tot}
0012 nNmaxESA=500; 
0013 
0014 % wavelengths (column)
0015 lambda=transpose(250:1:650);
0016 
0017 % dielectric constant of outside medium (real positive)
0018 epsilonM=1.77; % water
0019 
0020 % wavelength-dependent dielectric function of sphere material
0021 % vector column, same size as lambda [L x 1]
0022 epsilonIn=epsAg(lambda); % silver
0023 
0024 % Solve problems for all wavelength using Mie theory
0025 % for a=10nm, 25nm, and 50nm for d=1, d=2, and d=5 nm
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 % all necessary results are in structures stMdip
0032 
0033 %%%%%%%%%%%%%%%%%%
0034 % Drawing
0035 %%%%%%%%%%%%%%%%%%
0036 
0037 % create, position, and resize figure
0038 scrsz = get(0,'ScreenSize'); % scrsz(4) contains screen height
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 % plot Eta_{Rad}^{EM} for perp dipole (d=1nm)
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 % plot M_{Rad} and M_{Tot} for perp dipole (d=1nm)
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 % plot for a=25nm, Eta_{Rad}^{EM} for perp and para dipole (d=1,2,5 nm)
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 % plot M_{Rad} and M_{Tot} for perp dipole (d=1nm)
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}');

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)