Dear Steven,
It seems that Meep does not support slab structure with drude type metal or
dielectric.
I got the result about 3D photonic crystal (slab structure with the air hole)
dispersion curve. But, when I simulate the slab structure without the air hole,
the result is wrong.
Could you give me some suggestion to solve the problem?
Thank you very much!!
<3D metal slab case (.ctl file)>
(set! geometry-lattice (make lattice (size 1 1 10)))
(set! geometry
(list
(make block (center 0 0 0) (size 1 1 10)
(material (make dielectric (epsilon 1.0))))
(make block (center 0 0 0) (size 1 1 0.5)
(material (make dielectric (epsilon 1.0)
(polarizations (make polarizability (omega 1.0e-20) (gamma 0.15)
(delta-epsilon 10.0e+40 ) )))))))
(set! pml-layers (list (make pml (direction Z) (thickness 1))))
(set-param! resolution 40)
(define-param k-interp 15)
(define-param fcen 0.5)
(define-param df 1.0)
(set! sources
(list (make source (src (make gaussian-src (frequency fcen) (fwidth df)))
(component Hy) (center 0 0 0) )))
(run-k-points 300 (interpolate k-interp (list (vector3 0 0 0) (vector3 0.5 0 0)
)))
<3D dielectric slab case (.ctl file)>
(define-param r 0.3)
(define-param dpml 0.5)
(set! geometry-lattice (make lattice (size 1 1 7)))
(set! geometry
(list
(make block (center 0 0 0) (size 1 1 7)
(material (make dielectric (epsilon 1.0)))
)
(make block (center 0 0 0) (size 1 1 0.5)
(material (make dielectric (epsilon 12))))
))
(set! pml-layers (list (make pml (direction Z) (thickness dpml))))
(set-param! resolution 40)
(define-param fcen 0.3) ; pulse center frequency
(define-param df 0.6) ; pulse freq. width: large df = short impulse
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Hz) (center 0.123 0.1234 0))))
(define-param k-interp 15)
(run-k-points 300 (interpolate k-interp (list (vector3 0 0 0) (vector3 0.5 0 0)
))))
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
