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