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
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