Example script to run the function MulPweFullMonty for a silver nanoshell It produces 6 figures that illustrate the capabilities of the Mie codes for a sphere multilayer under PWE. See the function MulPweFullMonty 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 0010 % sphere radii in nm (scalar [1 x 1]) 0011 a1=24; 0012 a2=30; 0013 0014 % wavelengths in nm (column vector [L x 1]) 0015 lambda=transpose(200:1:800); 0016 0017 % dielectric constant of outside medium (real positive scalar) 0018 epsilonM=1.77; % water 0019 0020 % wavelength-dependent dielectric function of sphere material 0021 % vector column, same size as lambda [L x 1] 0022 epsilonShell=epsAg(lambda); % silver 0023 %epsilonIn=1.5^2; % silica 0024 epsilonIn=epsilonM; 0025 0026 Ca={a1,a2}; 0027 Cepsilon={epsilonIn,epsilonShell,epsilonM}; 0028 0029 nNbTheta=361; 0030 % selected wavelengths for plots 0031 lambda0=[618,487,301]; 0032 % Parameters to compute electric field everywhere at one wavelength 0033 % (lambda0(1)) 0034 rMax=4*a2; 0035 nNbRout=300; 0036 CnNbRin={100,50}; 0037 0038 [stMMul,stEsurfMul,stEfarScaMul,CstEmapMul]= ... 0039 MulPweFullMonty(nNmax,Ca,lambda,Cepsilon,lambda0,nNbTheta,rMax,nNbRout,CnNbRin);