程序是ddscat7.2 计算金纳米棒消光横截面。为什么只有520纳米左右的横向吸收峰,纵向表面等离子体吸收峰算不出来呢?请懂的大侠明示,非常感谢!
附计算出来的曲线。
附 程序
'**** Initial Memory Allocation ****'
100 100 100 = dimensioning allowance for target generation
'**** Target Geometry and Composition ****'
'cylinder1' = CSHAPE*9 shape directive
50 20 1 = shape parameters 1 - 6
1 = NCOMP = number of dielectric materials
'gold.par' = file with refractive index 1
'**** Additional Nearfield calculation? ****'
0 = NRFLD (=0 to skip nearfield calc., =1 to calculate nearfield E)
0.0 0.0 0.0 0.0 0.0 0.0 (fract. extens. of calc. vol. in -x,+x,-y,+y,-z,+z)
'**** Error Tolerance ****'
1.00e-5 = TOL = MAX ALLOWED (NORM OF |G>=AC|E>-ACA|X>/(NORM OF AC|E>
'**** Maximum number of iterations ****'
1000 = MXITER
'**** Integration limiter for PBC calculations ****'
1.00e-2 = GAMMA (1e-2 is normal, 3e-3 for greater accuracy)
'**** Angular resolution for calculation of , etc. ****'
0.5 = ETASCA (number of angles is proportional to [(3+x)/ETASCA]^2 )
'**** Wavelengths (micron) ****'
0.98413 0.18788 44 'TAB' = wavelengths (1st,last,howmany,how=LIN,INV,LOG,TAB)
'**** Refractive index of ambient medium ****'
1.33 = NAMBIENT
'**** Effective Radii (micron) **** '
0.015 0.015 1 'LIN' = eff. radii (1st,last,howmany,how=LIN,INV,LOG,TAB)
'**** Define Incident Polarizations ****'
(0,0) (1,0) (0,0) = Polarization state e01 (k along x axis)
1 = IORTH (=1 to do only pol. state e01; =2 to also do orth. pol. state)
'**** Specify which output files to write ****'
1 = IWRKSC (=0 to suppress, =1 to write ".sca" file for each target orient.
'**** Specify Target Rotations ****'
0. 0. 1 = BETAMI, BETAMX, NBETA (beta=rotation around a1)
0. 0. 1 = THETMI, THETMX, NTHETA (theta=angle between a1 and k)
0. 0. 1 = PHIMIN, PHIMAX, NPHI (phi=rotation angle of a1 around k)
'**** Specify first IWAV, IRAD, IORI (normally 0 0 0) ****'
0 0 0 = first IWAV, first IRAD, first IORI (0 0 0 to begin fresh)
'**** Select Elements of S_ij Matrix to Print ****'
6 = NSMELTS = number of elements of S_ij to print (not more than 9)
11 12 21 22 31 41 = indices ij of elements to print
'**** Specify Scattered Directions ****'
'LFRAME' = CMDFRM (LFRAME, TFRAME for Lab Frame or Target Frame)
2 = NPLANES = number of scattering planes
0. 0. 180. 5 = phi, theta_min, theta_max, dtheta(deg) for plane 1
90. 0. 180. 5 = phi, theta_min, theta_max, dtheta(deg) for plane 2
|