Your code setting the k points looks fine to me, so it must be
something wrong with the rest of your .ctl file.
Here is the ctl file I used for that figure (with kx of either 0 or
0.5 for the different curves).
(define-param eps 13)
(define-param eps2 1)
(define-param w 0.2)
(set! default-material (make dielectric (epsilon eps2)))
(set! geometry-lattice (make lattice (size 1 no-size no-size))) ; 1d
cell
(set! geometry
(list
(make cylinder
(material (make dielectric (epsilon eps)))
(center 0 0 0) (axis 1 0 0)
(radius infinity) (height w))))
(define-param kx 0.5)
(define-param ky-interp 20)
(define-param ky-max 1)
(set! k-points (interpolate ky-interp (list (vector3 kx 0 0)
(vector3 kx ky-max 0))))
(set-param! resolution 64)
(set-param! num-bands 4)
(run-tm randomize-fields)
_______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
