| 查看: 2255 | 回复: 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 高斯学习 |
» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有262人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
 |
8楼2012-10-25 22:55:23
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 (幼儿园)
- 金币: 4097.3
- 散金: 249
- 帖子: 1051
- 在线: 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
投票:
