Home > SPlaC v1_0 > BookFigures > MakeFig3_14b.m

MakeFig3_14b

PURPOSE ^

Function to reproduce Fig. 3.14(b)

SYNOPSIS ^

function MakeFig3_14b()

DESCRIPTION ^

 Function to reproduce Fig. 3.14(b)
 Angle dependence of the reflection coefficient at 633 nm for three
 dielectric refractive index close to water (nM=1.32, 1.33, or 1.34),
 for coupling to PSPPs at a silver/water interface with a sapphire prism
 using the Otto ATR configuration.
 Thickness of dielectric layer is 500 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=50:0.005:55; % row [1 x 1001]
0010 
0011 % gap thickness
0012 Lgap=500; % in nm
0013 
0014 % calculate epsilon for silver
0015 eAg=epsAg(lambda);
0016 % refractive indices
0017 nP=1.766; % Sapphire - Prism
0018 nM={1.32, 1.33, 1.34}; % cell of 3 scalars
0019 nNbnM=length(nM);
0020 
0021 % Initialize reflectivity cell Rp (one cell element for each nM)
0022 Rp=cell(1,nNbnM);
0023 
0024 % Defines Otto configuration, 2 interfaces, 3 layers
0025 % Prism/Dielectric/Metal
0026 % all epsilons are scalar here (one lambda only)
0027 nNbSurf=2;
0028 Cepsilon{1}=nP^2+0*lambda; % Prism
0029 Cepsilon{3}=eAg; % Metal
0030 CL{1}=0; % interfaces positions
0031 CL{2}=Lgap;
0032 
0033 % loop through each possible nM
0034 for ii=1:nNbnM
0035     Cepsilon{2}=nM{ii}^2+0*lambda; % Dielectric for PSPP
0036     % solves the problem for TM waves
0037     stResTM=PlnMultiRef ('TM',nNbSurf, lambda, Cepsilon, CL, aideg);
0038     % extract reflectivity for this nM: rP is [1 x 1001]
0039     Rp{ii}=abs(stResTM.rP).^2;
0040 end
0041 
0042 %%%%%%%%%%%%%%%%%%
0043 % Drawing
0044 %%%%%%%%%%%%%%%%%%
0045 
0046 % create, position, and resize figure
0047 scrsz = get(0,'ScreenSize'); % scrsz(3) contains screen width
0048 figAspectRatio=1.5;
0049 figWidth=3/4*scrsz(3);
0050 figure('Name','Fig. 3.14(b)', ...
0051     'Position',[(scrsz(3)-figWidth)/2 scrsz(4)-150-figWidth/figAspectRatio figWidth figWidth/figAspectRatio]);
0052 
0053 % create plot
0054 plot(aideg,Rp{1},aideg,Rp{2},aideg,Rp{3});
0055 title('Dielectric/Ag PSPP excitation in the Otto configuration - L_{gap}=500 nm, \lambda=633 nm');
0056 % format graph
0057 axis([50.5 55 0 1.15]);
0058 xlabel('angle of incidence [degrees]');
0059 ylabel('Reflectivity');
0060 legend('n_M=1.32','n_M=1.33','n_M=1.34','Location','East');
0061 
0062 
0063 
0064 
0065

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)