Dear Steven. J
I have some questions to ask about you..
Firstly: I have simulated the PRL. 78. 3294. Fig1, but I found that my result is not correct My control file as follow:
(define-param n 3.5) ;dielectric constant
(define-param eps (* n n))
(set! geometry-lattice (make lattice (size 5 5 3))) ；the computional cell
;Geometry
(set! geometry
(list
(make block
(center 0 0 0)
(size 5 5 0.5) ; the slab size 5*5*0.5
(material (make dielectric (epsilon eps))))))
;Source
(set! sources
(list
(make source
(src (make gaussian-src (frequency 0.55) (fwidth 0.7)))
(component Ez) ; Is the result different if I use ‘Hz’ ’Ex’ or ‘Hx’, I do not know How it
; inffluence the result.
(center 0 0 0))))
(set! pml-layers (list (make pml (thickness 1.0))))
(set-param! resolution 20)
(define trans
(add-flux fcen df nfreq
(make flux-region (center 0 0 0.3) (size 5 5)))) ;Is the flux region correct?
(run-sources+ (stop-when-fields-decayed 50 Ez
(vector3 0 0 0.3) 1e-2))
(display-fluxes trans)
I computed the transmitted flux at the situation n=1 and n=3.5. But the result is wrong.
Secondly: I want to compute a transmission spectrum around a block. The control file as follow:
(set! geometry-lattice (make lattice (size 5 7 7))) ；the computional cell
;Geometry
(set! geometry
(list
(make block
(center 0 0 0)
(size 1 7 7) ; the block size 1*7*7
(material (make dielectric (epsilon eps))))))
;Source
(set! sources (list (make source (src (make continuous-src (frequency 0.15))) (component Ez) (center -1 0))))(set! resolution 10)(set! pml-layers (list (make pml (thickness 1.0))))
(define trans
(add-flux
(make flux-region (center 1 0 0) (size 0 7 7))))
(run-sources+ (stop-when-fields-decayed 50 Ez
(vector3 1 0 0))
(at-beginning output-epsilon))
(display-fluxes trans)
Is it correct?When I compute transmission and reflection spectrum, I can not set the source and the computed flux region at a correct location. Could you give me some help?
Thank you very much!
你 不 想 试 试 今 夏 最 “酷” 的 邮 箱 吗 ？
蕴 涵 中 华 传 统 文 化 于 世 界 一 流 科 技 之 中，创 新 Ajax 技 术，126 “D 计 划”火 热 体 验 中 ！
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
