MNPBEM.pptVIP

二、模拟由下面几部分组成 1、定义介电函数 2、定义粒子边界 3、指定粒子如何嵌入介电环境 4、建立边界元方程 5、指定激励方法 6、由给定的激励通过计算电流或电荷来解边界基元方程 7、计算等离子纳米粒子对外部激励的响应 1、通过语句定义介电函数和球面 % table of dielectric functions epstab={epsconst(1.33^2),epstable(’gold.dat’)}; % nanosphere with 144 vertices and 10 nanometers diameter p = trisphere(144,10); comparticle命令定义了粒子边界是如何嵌入介电环境中的。 p=comparticle({eps1,eps2,...},{p1,p2,...},inout,closed); 由图1我们可以得到： % single sphere of Fig. 1(a) p=comparticle({eps1,eps2},{p},[2,1],1); % coupled spheres of Fig. 1(b) p=comparticle({eps1,eps2,eps3},{p1,p2},[2,1;3,1],1,2); % coated particles of Fig. 1(c) p=comparticle({eps1,eps2,eps3},{p1,p2},[2,1;3,2],1,2); 3、通过语句建立BEM方程 % quasistatic BEM solver bem = bemstat(p); 4、通过语句指定外部激励 % plane wave excitation for given light polarizations exc = planewavestat([1,0,0;0,1,0]); % planewave excitation e = exc(p,enei); 5、通过语句 % compute surface charge sig = bem \ e; 计算便面的电流来解BEM方程。 图3 1、如图3所示，模拟一个金属椭球的光散射谱的准静态近似。 % table of dielectric functions epsb = epsconst( 1.33 ^ 2 ); epsau = epstable( gold.dat ); epstab = { epsb, epsau }; % initialize sphere sph = trisphere( 256, 10 ); % plane wave excitation, polarization along long axis of ellipsoid exc = planewavestat( [ 0, 0, 1 ] ); % light wavelength in vacuum enei = linspace( 400, 800, 200 ); % scaling factor for long axis of ellipsoid scaling = [ 2, 3 ]; % scattering spectra sca =zeros( length( enei ), length( scaling ), %set up BEM solver in quasistatic approximation with eigenmode expansion bem = bemstateig( p, 20 ); for ien = 1 : length( enei ) sig = bem \ exc( p, enei( ien ) ); sca( ien, is, ic ) = exc.sca( sig ); end 2、 图4使用了完整的麦克斯韦方程，用边界基元法模拟了极化方向沿x轴和z轴的偶极子的总衰减率和辐射衰减率。 图 4 % table of dielectric functions epsb = epsconst( 1.33 ^ 2 ); epsau = epstable( gold.dat ); epstab = { epsb, epsau }; % diameter of sphere diameter = 60; % initialize sphere % use a higher number of vertices to