| 查看: 2083 | 回复: 8 | |||
bin2boy木虫 (正式写手)
|
[求助]
关于ddscat7.2模拟计算问题
|
请教一下关于用ddscat7.2模拟计算的问题，我在模拟一个有效半径为15nm的银纳米椭球（长轴=30nm，短轴=10nm），文献上分别给出光沿长轴、短轴方向偏振以及两个方向的平均结果，如图， 请问如何设置这个average方向会出现两个消光峰？我把我的代码贴上来请帮我看看如何设置： ' ========== Parameter file for v7.2 ===================' '**** Preliminaries ****' 'NOTORQ' = CMDTRQ*6 (DOTORQ, NOTORQ) -- either do or skip torque calculations 'PBCGS2' = CMDSOL*6 (PBCGS2, PBCGST, GPBICG, QMRCCG, PETRKP) -- solution method 'GPFAFT' = CMETHD*6 (GPFAFT, FFTMKL) -- FFT implementation 'GKDLDR' = CALPHA*6 (GKDLDR, LATTDR) -- recipe for dipole polarizabilities 'NOTBIN' = CBINFLAG (NOTBIN, ORIBIN, ALLBIN) -- binary output? '**** Initial Memory Allocation ****' 100 100 100 = dimensioning allowance for target generation '**** Target Geometry and Composition ****' 'ELLIPSOID' = CSHAPE*9 shape directive 60 20 22.5 = shape parameters 1 - 3 1 = NCOMP = number of dielectric materials 'Ag_evap' = 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 (fract. extension of 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 ****' 100000 = 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 ) '**** Vacuum Wavelengths (micron) ****' 0.3 0.8 40 'LIN' = wavelengths (1st,last,howmany,how=LIN,INV,LOG,TAB) '**** Refractive index of ambient medium ****' 1.0000 = NAMBIENT '**** Effective Radii (micron) **** ' 0.015 0.015 1 'INV' = a_eff (1st,last,howmany,how=LIN,INV,LOG,TAB) '**** Define Incident Polarizations ****' (0,0) (1.,0.) (0.,0.) = Polarization state e01 (k along x axis) 2 = 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) 1 = NPLANES = number of scattering planes 0. 0. 180. 1 = phi, theta_min, theta_max (deg) for plane A |
回复此楼» 收录本帖的淘帖专辑推荐
 高斯学习 |
» 猜你喜欢
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
-
大豆异黄酮分离
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有124人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有1人回复
-
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有4人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
tianwaiflyke
金虫 (正式写手)
- 应助: 26 (小学生)
- 金币: 3680.7
- 红花: 1
- 帖子: 832
- 在线: 259.7小时
- 虫号: 1505177
- 注册: 2011-11-22
- 专业: 凝聚态物性 II ：电子结构
2楼2012-10-24 06:41:35
tianwaiflyke
金虫 (正式写手)
- 应助: 26 (小学生)
- 金币: 3680.7
- 红花: 1
- 帖子: 832
- 在线: 259.7小时
- 虫号: 1505177
- 注册: 2011-11-22
- 专业: 凝聚态物性 II ：电子结构
★
fegg7502: 金币+1, 鼓励交流 2012-10-24 10:40:35
fegg7502: 金币+1, 鼓励交流 2012-10-24 10:40:35
|
我也有个问题求教楼主，望多多赐教： http://muchong.com/bbs/viewthread.php?tid=5096750&authorid=1505177 非常感谢！ |
3楼2012-10-24 06:43:42
zclian 
银虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4092.3
- 散金: 249
- 帖子: 1050
- 在线: 909.7小时
- 虫号: 1407064
- 注册: 2011-09-18
- 专业: 污染控制化学
【答案】应助回帖
感谢参与，应助指数 +1
|
我也需要此软件，请问你安装的Linux系统是哪个版本的？gfortran最新的编译器吗？或者用的是ifort编译器吗？ 应改变此行的(0,0) (1.,0.) (0.,0.) = Polarization state e01 (k along x axis)，极化方向有关。我以前和导师算过 我的邮箱zichao8806@163.com |
4楼2012-10-25 17:30:59
bin2boy
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 4457.8
- 散金: 42
- 红花: 4
- 沙发: 1
- 帖子: 608
- 在线: 39.1小时
- 虫号: 1615407
- 注册: 2012-02-13
- 性别: GG
- 专业: 光学
5楼2012-10-25 20:18:36
zclian
银虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4092.3
- 散金: 249
- 帖子: 1050
- 在线: 909.7小时
- 虫号: 1407064
- 注册: 2011-09-18
- 专业: 污染控制化学
6楼2012-10-25 22:18:32
zclian
银虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4092.3
- 散金: 249
- 帖子: 1050
- 在线: 909.7小时
- 虫号: 1407064
- 注册: 2011-09-18
- 专业: 污染控制化学
7楼2012-10-25 22:45:20
 |
8楼2012-10-25 22:55:23
bin2boy
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 4457.8
- 散金: 42
- 红花: 4
- 沙发: 1
- 帖子: 608
- 在线: 39.1小时
- 虫号: 1615407
- 注册: 2012-02-13
- 性别: GG
- 专业: 光学
9楼2012-10-25 23:10:32