Home > SPlaC v1_0 > BookFigures > MakeFig3_10.m

MakeFig3_10

PURPOSE ^

Function to reproduce Fig. 3.10

SYNOPSIS ^

function MakeFig3_10()

DESCRIPTION ^

 Function to reproduce Fig. 3.10
 Propagation/decay lengths (normalized to the wavelength) for surface
 modes for a silver/air (left) and a gold/air (right) planar interface. 

 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 % wavelengths (column)
0006 lambda=transpose(200:2:1500); 
0007 
0008 % epsilon
0009 % metals: silver and gold
0010 eAg=epsAg(lambda); % column
0011 eAu=epsAu(lambda); % column
0012 % dielectric (air here)
0013 eM=1.0; % scalar
0014 
0015 % Get kx from dispersion relation (Eq. 3.15)
0016 k0=2*pi./lambda; % column
0017 kxAg=k0.*sqrt(eM*eAg./(eM+eAg));
0018 kxAu=k0.*sqrt(eM*eAu./(eM+eAu));
0019 
0020 % Get propagation length LSPP from Eq. 3.16
0021 LSPPAg=0.5./abs(imag(kxAg));
0022 LSPPAu=0.5./abs(imag(kxAu));
0023 
0024 % Get k1z outside (in air) and k2z inside (the metal)
0025 % using Eq. F.9/F.10
0026 % These are only valid up to a sign, which is irrelevant here
0027 % see Fig. 3.7 and Sec. F.2.5 for more detail on this point
0028 k1zOutAg=sqrt(k0.^2*eM-kxAg.^2);
0029 k2zInAg=sqrt(k0.^2.*eAg-kxAg.^2);
0030 
0031 k1zOutAu=sqrt(k0.^2*eM-kxAu.^2);
0032 k2zInAu=sqrt(k0.^2.*eAu-kxAu.^2);
0033 
0034 % Get decay lengths L1z outside (in air) and L2z inside (the metal)
0035 L1zOutAg=0.5./abs(imag(k1zOutAg));
0036 L2zInAg=0.5./abs(imag(k2zInAg));
0037 
0038 L1zOutAu=0.5./abs(imag(k1zOutAu));
0039 L2zInAu=0.5./abs(imag(k2zInAu));
0040 
0041 
0042 %%%%%%%%%%%%%%%%%%
0043 % Drawing
0044 %%%%%%%%%%%%%%%%%%
0045 
0046 % create, position, and resize figure
0047 scrsz = get(0,'ScreenSize'); % scrsz(3) contains screen width
0048 figAspectRatio=3;
0049 figWidth=3/4*scrsz(3);
0050 figure('Name','Fig. 3.10', ...
0051     'Position',[(scrsz(3)-figWidth)/2 scrsz(4)-150-figWidth/figAspectRatio figWidth figWidth/figAspectRatio]);
0052 
0053 % plot for Ag/air
0054 subplot(1,2,1);
0055 semilogy(lambda,LSPPAg./lambda,'-',lambda,L1zOutAg./lambda,':',lambda,L2zInAg./lambda,'--');
0056 % format graph
0057 axis([200 1500 0.005 500]);
0058 xlabel('Wavelength [nm]');
0059 ylabel('L/\lambda');
0060 legend('L_{SPP}/\lambda','L_{1z}/\lambda','L_{2z}/\lambda','Location','Best');
0061 title('Propagation/decay lengths - PSPPs on air/silver interface');
0062 set(gca,'YGrid','on','YMinorGrid','off','GridLineStyle','-');
0063 
0064 % plot for Ag/air
0065 subplot(1,2,2);
0066 semilogy(lambda,LSPPAu./lambda,'-',lambda,L1zOutAu./lambda,':',lambda,L2zInAu./lambda,'--');
0067 % format graph
0068 axis([200 1500 0.005 500]);
0069 xlabel('Wavelength [nm]');
0070 ylabel('L/\lambda');
0071 legend('L_{SPP}/\lambda','L_{1z}/\lambda','L_{2z}/\lambda','Location','Best');
0072 title('Propagation/decay lengths - PSPPs on air/gold interface');
0073 set(gca,'YGrid','on','YMinorGrid','off','GridLineStyle','-');
0074 
0075 
0076 
0077

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)