Home > SPlaC v1_0 > BookFigures > MakeFig3_14c.m

MakeFig3_14c

PURPOSE ^

Function to reproduce Fig. 3.14(c)

SYNOPSIS ^

function MakeFig3_14c()

DESCRIPTION ^

 Function to reproduce Fig. 3.14(c)
 Wavelength-dependence of the reflection coefficient at theta_inc=53º for three
 dielectric refractive index close to water (nM=1.32, 1.33, or 1.34),
 for coupling to PSPPs at a silver/dielectric 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 (column)
0006 lambda=transpose(500:0.5:800); % [601 x 1]
0007 
0008 % angles in degrees
0009 aideg=53; % [1 x 1]
0010 
0011 % metal thickness
0012 Lmetal=50; % 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 Krestschmann configuration, 2 interfaces, 3 layers
0025 % Prism/Metal/Dielectric
0026 % all epsilons are column (lambda-dependent)
0027 nNbSurf=2;
0028 Cepsilon{1}=nP^2+0*lambda; % Prism
0029 Cepsilon{2}=eAg; % Metal
0030 CL{1}=0; % interfaces positions
0031 CL{2}=Lmetal;
0032 
0033 % loop through each possible nM
0034 for ii=1:nNbnM
0035     Cepsilon{3}=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 [601 x 1]
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(c)', ...
0051     'Position',[(scrsz(3)-figWidth)/2 scrsz(4)-150-figWidth/figAspectRatio figWidth figWidth/figAspectRatio]);
0052 
0053 % create plot
0054 plot(lambda,Rp{1},lambda,Rp{2},lambda,Rp{3});
0055 title('Dielectric/Ag PSPP excitation in the Kretschmann configuration - L_{metal}=50 nm, \theta_{inc}=53^\circ');
0056 % format graph
0057 axis([500 800 0 1.15]);
0058 xlabel('Wavelength [nm]');
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)