I have some problems calculate a band diagram of metallic pc.
When I calculate my own ctl code, I couldn't get modes. What is the problem?
/* ctl file*/
(set! geometry-lattice (make lattice (size 1 1 no-size)))
(set! geometry
(list
(make cylinder (center 0 0 0) (radius 0.3) (height infinity)
(material (make dielectric (epsilon 1) (polarizations (make polarizability (omega 1e-20) (gamma 0.051603) (delta-epsilon 3.6372e+41) ))))
)
)
)
(set-param! resolution 20)
(define-param k-interp 4)
(define-param fcen 0.1)
(define-param df 5)
(set! sources
(list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ez) (center 0.48 0.1234)
)
)
)
(set! symmetries
(list
(make mirror-sym (direction X) (phase -1) )
(make mirror-sym (direction Y) (phase -1) )
)
)
(run-k-points 300 (interpolate k-interp (list (vector3 0 0 0) (vector3 0.5 0 0) (vector3 0.5 0.5 0) (vector3 0 0 0))))
/* result file(grep freqs:) */
freqs:, 1, 0, 0, 0
freqs:, 2, 0.1, 0.0, 0.0
freqs:, 3, 0.2, 0.0, 0.0
freqs:, 4, 0.3, 0.0, 0.0
freqs:, 5, 0.4, 0.0, 0.0
freqs:, 6, 0.5, 0, 0
freqs:, 7, 0.5, 0.1, 0.0
freqs:, 8, 0.5, 0.2, 0.0
freqs:, 9, 0.5, 0.3, 0.0
freqs:, 10, 0.5, 0.4, 0.0
freqs:, 11, 0.5, 0.5, 0
freqs:, 12, 0.4, 0.4, 0.0
freqs:, 13, 0.3, 0.3, 0.0
freqs:, 14, 0.2, 0.2, 0.0
freqs:, 15, 0.1, 0.1, 0.0
freqs:, 16, 0, 0, 0
--
Ki-young Choi
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
