Hi all,
I want to calculate the dispersion of metal-dielectric-metal structure, i.e.,
Fig.1 of PRL 96, 073907 (2006) by H Shin & SH Fan.
My ctl file attached below [Drude model eps(w)=1-1.0/w/w assumed]. I can only
get the first and third band. Please tell me what's wrong. Thanks in advance.
==========================================
(define-param sz 5) ; size of cell in z direction
(define-param dpml 0.5) ; PML thickness
(set! geometry-lattice (make lattice (size no-size no-size sz)))
(set! eps-averaging? false)
(set! geometry
(list (make block (center 0 0 0) (size infinity infinity sz)
(material (make dielectric (epsilon 1.0)
(polarizations (make polarizability
(omega 1e-20) (gamma 0) (delta-epsilon 1e40))))))
(make block (center 0 0 0) (size infinity infinity 0.1)
(material (make dielectric (epsilon 4.0))))
))
(set! pml-layers (list (make pml (direction Z) (thickness dpml))))
(set-param! resolution 80)
(define-param fcen 0.25)
(define-param df 2.0)
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ez) (center 0 0 0))))
(define-param k-interp 8)
;square lattice
(define-param k-points (list (vector3 0 0 0) ; Gamma
(vector3 2 0 0) ; X
))
(set! k-points (interpolate k-interp k-points))
(run-k-points 300 k-points)
---------------------------------
Don't be flakey. Get Yahoo! Mail for Mobile and
always stay connected to friends.
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
