hi everyone!
below is the .ctl file I use to calculate the dispersion relation of
surface plasmon
I can get the dispersion relation when the source component set as Hz and
ky run from kmin to kmax.
the queer problem is that
when I set the source component as Hy, and kz run from kmin to kmax, I can
not get the dispersion relation.
Anyway,both configurations are the same for two infinite metal and air
layer.
anyone have any idea about that?
thank!
(define-param sx 4)
(set! geometry-lattice (make lattice (size sx no-size no-size)))
(set-param! resolution 40)
;(set! eps-averaging? false)
(define omega_eps 0.3)
(define Ag51 (make dielectric (epsilon 6.8)
(polarizations
(make polarizability
(omega 1e-20) (gamma 0.003226) (sigma (* 6.8 omega_eps omega_eps 1e+40)))
)))
(set! geometry (list
(make block (center 0 0 0) (size sx infinity infinity)
(material Ag51))
(make block (center 1 0 0) (size 2 infinity infinity)
(material air))
))
(define-param fcen 0.2) ; pulse center frequency
(define-param df 0.4) ; pulse width (in frequency)
(define-param runtime 200)
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Hz)
(center 0.0 0 0))
))
;(use-output-directory)
(define-param kmin 0.0)
(define-param kmax 3.0)
(define-param k-interp 10)
(define kpts (interpolate k-interp (list (vector3 0 kmin 0) (vector3 0
kmax 0 ))))
(define all-freqs (run-k-points runtime kpts)) ; a list of lists of
frequencies
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
