Home > SPlaC v1_0 > Mie > MUL > MulDipScriptFullMonty.m

MulDipScriptFullMonty

PURPOSE ^

Example script to run the function MulDipFullMonty

SYNOPSIS ^

This is a script file.

DESCRIPTION ^

 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

EXAMPLE OF OUTPUT ^

Example figure output

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

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

This web page is part of the SPlaC package © 2008. Contact: Eric Le Ru
Generated on Wed 03-Dec-2008 11:10:14 by m2html © 2003 (adapted)