0001
0002
0003
0004
0005
0006 lambda=transpose(200:2:1500);
0007
0008
0009
0010 eAg=epsAg(lambda);
0011 eAu=epsAu(lambda);
0012
0013 eM=1.0;
0014
0015
0016 k0=2*pi./lambda;
0017 kxAg=k0.*sqrt(eM*eAg./(eM+eAg));
0018 kxAu=k0.*sqrt(eM*eAu./(eM+eAu));
0019
0020
0021 LSPPAg=0.5./abs(imag(kxAg));
0022 LSPPAu=0.5./abs(imag(kxAu));
0023
0024
0025
0026
0027
0028 k1zOutAg=sqrt(k0.^2*eM-kxAg.^2);
0029 k2zInAg=sqrt(k0.^2.*eAg-kxAg.^2);
0030
0031 k1zOutAu=sqrt(k0.^2*eM-kxAu.^2);
0032 k2zInAu=sqrt(k0.^2.*eAu-kxAu.^2);
0033
0034
0035 L1zOutAg=0.5./abs(imag(k1zOutAg));
0036 L2zInAg=0.5./abs(imag(k2zInAg));
0037
0038 L1zOutAu=0.5./abs(imag(k1zOutAu));
0039 L2zInAu=0.5./abs(imag(k2zInAu));
0040
0041
0042
0043
0044
0045
0046
0047 scrsz = get(0,'ScreenSize');
0048 figAspectRatio=3;
0049 figWidth=3/4*scrsz(3);
0050 figure('Name','Fig. 3.10', ...
0051 'Position',[(scrsz(3)-figWidth)/2 scrsz(4)-150-figWidth/figAspectRatio figWidth figWidth/figAspectRatio]);
0052
0053
0054 subplot(1,2,1);
0055 semilogy(lambda,LSPPAg./lambda,'-',lambda,L1zOutAg./lambda,':',lambda,L2zInAg./lambda,'--');
0056
0057 axis([200 1500 0.005 500]);
0058 xlabel('Wavelength [nm]');
0059 ylabel('L/\lambda');
0060 legend('L_{SPP}/\lambda','L_{1z}/\lambda','L_{2z}/\lambda','Location','Best');
0061 title('Propagation/decay lengths - PSPPs on air/silver interface');
0062 set(gca,'YGrid','on','YMinorGrid','off','GridLineStyle','-');
0063
0064
0065 subplot(1,2,2);
0066 semilogy(lambda,LSPPAu./lambda,'-',lambda,L1zOutAu./lambda,':',lambda,L2zInAu./lambda,'--');
0067
0068 axis([200 1500 0.005 500]);
0069 xlabel('Wavelength [nm]');
0070 ylabel('L/\lambda');
0071 legend('L_{SPP}/\lambda','L_{1z}/\lambda','L_{2z}/\lambda','Location','Best');
0072 title('Propagation/decay lengths - PSPPs on air/gold interface');
0073 set(gca,'YGrid','on','YMinorGrid','off','GridLineStyle','-');
0074
0075
0076
0077