Home > SPlaC v1_0 > BookFigures > MakeFig3_15b.m

MakeFig3_15b

PURPOSE ^

Function to reproduce Fig. 3.15(b)

SYNOPSIS ^

function MakeFig3_15b()

DESCRIPTION ^

 Function to reproduce Fig. 3.15(b)
 Angle dependence of the reflection coefficient and local field intensity
 EF at 633 nm for coupling to PSPPs at a silver/air interface with a sapphire
 prism using the Kretschmann ATR configuration.
 Thickness of metal layer is 50 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 % lambda in nm
0006 lambda=633; % [1 x 1]
0007 
0008 % angles in degrees (row)
0009 aideg=34.5:0.005:37; % row [1 x 601]
0010 
0011 % metal thickness
0012 Lmetal=50; % in nm
0013 
0014 % calculate epsilon for silver
0015 eAg=epsAg(lambda); % [1 x 1]
0016 % refractive indices
0017 nP=1.766; % Sapphire - Prism
0018 nM=1.0; % Air - Dielectric for PSPPs
0019 
0020 % Defines Kretschmann configuration, 2 interfaces, 3 layers
0021 % Prism/Metal/Dielectric
0022 nNbSurf=2;
0023 Cepsilon{1}=nP^2+0*lambda; % Prism
0024 Cepsilon{2}=eAg; % Metal
0025 Cepsilon{3}=nM^2+0*lambda; % Dielectric for PSPP
0026 CL{1}=0; % interfaces positions
0027 CL{2}=Lmetal;
0028 
0029 % solves the problem for TM waves
0030 stResTM=PlnMultiRef ('TM',nNbSurf, lambda, Cepsilon, CL, aideg);
0031 
0032 
0033 % Reflectivity (reflection coefficient $R^p=|r^p|^2$)
0034 % Rp is row [1 x 601]
0035 Rp=abs(stResTM.rP).^2;
0036 % Local field EF on metal/air interface, i.e. interface 2, on the air side
0037 % MLoc is row [1 x 601]
0038 MLoc=stResTM.MinPerp{2}+stResTM.MinPara{2};
0039 
0040 % Local field EF from Eq. 3.27 [1 x 601]
0041 MLocEq3_27=2*abs(eAg)/(nM^2 * imag(eAg)) * imag(eAg*sqrt(1/(eAg+nM^2)))*nP*cos(aideg*pi/180).*(1-Rp);
0042 
0043 %%%%%%%%%%%%%%%%%%
0044 % Drawing
0045 %%%%%%%%%%%%%%%%%%
0046 
0047 % create, position, and resize figure
0048 scrsz = get(0,'ScreenSize'); % scrsz(3) contains screen width
0049 figAspectRatio=1.5;
0050 figWidth=3/4*scrsz(3);
0051 figure('Name','Fig. 3.15(b)', ...
0052     'Position',[(scrsz(3)-figWidth)/2 scrsz(4)-150-figWidth/figAspectRatio figWidth figWidth/figAspectRatio]);
0053 
0054 line(aideg,Rp,'Color','b');
0055 ax1 = gca;
0056 set(ax1,'XColor','k','YColor','b')
0057 set(ax1,'XLimMode','Manual','YLimMode','Manual');
0058 set(ax1,'XLim',[34.5 37],'YLim',[0 1.15]);
0059 set(get(ax1,'XLabel'),'String','Angle of incidence [degrees]');
0060 set(get(ax1,'YLabel'),'String','Reflectivity');
0061 legend('R','Location','NorthWest');
0062 title('Air/Ag PSPP excitation in the Kretschmann configuration - L_{metal}=50 nm, \lambda=633 nm');
0063 
0064 ax2=axes('Position',get(ax1,'Position'),...
0065            'XAxisLocation','top','XColor','k', ...
0066            'YAxisLocation','right','Ycolor','r', ...
0067            'Color','none', ...
0068            'YScale','Log', ...
0069            'XTickLabel',{}, ...
0070            'XLimMode','Manual','YLimMode','Manual', ... 
0071            'XLim',[34.5 37],'YLim',[0.2 800]);
0072 set(get(ax2,'Ylabel'),'String','Local field intensity EF');
0073 
0074 line(aideg,MLoc,'Color','r');
0075 indMarker=10:20:length(aideg);
0076 line(aideg(indMarker),MLocEq3_27(indMarker),'Color','r','Marker','o','LineStyle','none');
0077 legend({'M_{Loc}','M_{Loc} from Eq. 3.27'},'Location','SouthEast','Color','w');
0078        
0079

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)