Dear Steven,
My codes is based on an example which is used to solve the bandstructure of
holey waveguide
http://ab-initio.mit.edu/wiki/index.php/Meep_Tutorial/Band_diagram%2C_resonant_modes%2C_and_transmission_in_a_holey_waveguide
What I do is changing Si to Ag. The codes is as follows,
-------------------------------------------------------------------------------------------------
; Some parameters to describe the geometry:
(define-param eps 13) ; dielectric constant of waveguide
(define-param w 1.2) ; width of waveguide
(define-param r 0.36) ; radius of holes
; The cell dimensions
(define-param sy 12) ; size of cell in y direction (perpendicular to wvg.)
(define-param dpml 1) ; PML thickness (y direction only!)
(set! eps-averaging? false)
(set! geometry-lattice (make lattice (size 1 sy no-size)))
(define-param Ag
(make dielectric (epsilon 1)
(polarizations
(make polarizability
(omega 1e-20) (gamma 0) (delta-epsilon 1e40))
)))
(set! geometry
(list (make block (center 0 0) (size infinity w infinity)
(material Ag))
(make cylinder (center 0 0) (radius r) (height infinity) (material
air))))
(set! pml-layers (list (make pml (direction Y) (thickness dpml))))
(set-param! resolution 40)
(define-param fcen 0.25) ; pulse center frequency
(define-param df 1.5) ; pulse freq. width: large df = short impulse
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Hz) (center 0.1234 0))))
(set! symmetries (list (make mirror-sym (direction Y) (phase -1))))
(define-param kx false) ; if true, do run at specified kx and get fields
(define-param k-interp 19) ; # k-points to interpolate, otherwise
(if kx
(begin
(set! k-point (vector3 kx))
(run-sources+
300 (at-beginning output-epsilon)
(after-sources (harminv Hz (vector3 0.1234 0) fcen df)))
(run-until (/ 1 fcen) (at-every (/ 1 fcen 20) output-hfield-z)))
(run-k-points 300 (interpolate k-interp (list (vector3 0) (vector3 0.5)))))
------------------------------------------------------------------------------------------------------------
Why holey-wvg-bands-eps-000000.00.png is black? If I use PEC, I can get the
structure image.
Thank you very much!
mario
__________________________________________________
Do You Yahoo!?
Tired of spam? Yahoo! Mail has the best spam protection around
http://mail.yahoo.com _______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
