Example script to run the function MulDipFullMonty It produces 4 figures that illustrate the capabilities of the Mie codes for dipole emission close to a silver nanoshell. See the function MulDipFullMonty for more details. This file is part of the SPlaC v1.0 package (copyright 2008) Check the README file for further information
0001 global noCheckSum; 0002 noCheckSum=true; 0003 % use 'clear all' or 'clear global noChecksum' to erase noCheckSum 0004 % clear global noCheckSum; % same as noCheckSum=false; 0005 % noCheckSum=false; 0006 0007 % N_{max} 0008 nNmax=50; 0009 % N_{max} for ESA approximation in the evaluations of M_{Tot} 0010 nNmaxESA=500; 0011 0012 % sphere radii in nm (scalar [1 x 1]) 0013 a1=24; 0014 a2=30; 0015 0016 % wavelengths in nm (column vector [L x 1]) 0017 lambda=transpose(200:1:800); 0018 0019 % dielectric constant of outside medium (real positive scalar) 0020 epsilonM=1.77; % water 0021 0022 % wavelength-dependent dielectric function of sphere material 0023 % vector column, same size as lambda [L x 1] 0024 epsilonShell=epsAg(lambda); % silver 0025 %epsilonIn=1.5^2; % silica 0026 epsilonIn=epsilonM; 0027 0028 Ca={a1,a2}; 0029 Cepsilon={epsilonIn,epsilonShell,epsilonM}; 0030 0031 % distance of dipole from surface (in nanometer) 0032 % note that the dipole is by convention on the positive z-axis (theta=0) 0033 % and aligned in the (xOz) plane (py=0) 0034 d=1; 0035 0036 % Number of theta for radiation profile computations 0037 nNbTheta=361; 0038 % selected wavelengths for radiation profile plots 0039 lambda0=[618,487,301]; 0040 0041 [stMdip,stEfarAllP]=MulDipFullMonty(nNmax,Ca,lambda,Cepsilon,d,nNmaxESA,lambda0,nNbTheta); 0042 0043 % % uncomment the following to check MRadPerp and MRadPara 0044 % % from averages of radiated field 0045 % figure; 0046 % plot(lambda,stMdip.MRadPara,lambda,stEfarAllP.MRadPara, ... 0047 % lambda,stMdip.MRadPerp,lambda,stEfarAllP.MRadPerp); 0048