| 查看: 696 | 回复: 2 | ||
cared1000铁杆木虫 (著名写手)
|
[求助]
ddscat 求助
|
|
程序是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 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  [ Last edited by cared1000 on 2012-6-22 at 22:22 ] |
回复此楼» 猜你喜欢
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
物理学I论文润色/翻译怎么收费?
已经有188人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
» 本主题相关商家推荐: (我也要在这里推广)
cared1000
铁杆木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 6037.4
- 红花: 2
- 帖子: 1671
- 在线: 251.6小时
- 虫号: 1172583
- 注册: 2010-12-18
- 性别: GG
- 专业: 光学
2楼2012-06-23 12:29:24
cared1000
铁杆木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 6037.4
- 红花: 2
- 帖子: 1671
- 在线: 251.6小时
- 虫号: 1172583
- 注册: 2010-12-18
- 性别: GG
- 专业: 光学
3楼2012-06-23 20:24:21
50