How to set k-points in piperiod pc waveguide?

Dear Steven and MPB's users:
  I design a PC waveguide , but the period of waveguide is not equal to the period of crystal .When i calculate the mode of  the waveguide , how to set the k-point ?
the structure of  PC  waveguide:
 *    *    *    *    *    *   *    *    *   *    *    *   *    *    *   *    *    *   *    *    *   *    *    *   *    *    *   *    *    *
 *    *    *    *    *    *   *    *    *   *    *    *   *    *    *   *    *    *   *    *    *   *    *    *   *    *    *   *    *    *
 * * * * * * * * * * * * * ** * * * * * ** * * * * * ** * * * * * ** * * * * * ** * * * * * ** * * * * * ** * * 
 * * * * * * * * * * * * * ** * * * * * ** * * * * * ** * * * * * ** * * * * * ** * * * * * ** * * * * * ** * *
 *    *    *    *    *    *   *    *    *   *    *    *   *    *    *   *    *    *   *    *    *   *    *    *   *    *    *   *    *    *
 *    *    *    *    *    *   *    *    *   *    *    *   *    *    *   *    *    *   *    *    *   *    *    *   *    *    *   *    *    *
 My ctl 
 (define-param h  0.5) ; widelength of  waveguide
 (define-param r 0.3)    ;
 (define-param rd 0.2)
 (define-param eps 1)
(define-param supercell-y (+ 8 h))
(set! default-material (make dielectric (epsilon 11.56)))
(set! geometry-lattice (make lattice 
    (basis1 1 0)
    (basis2 0 1)
    (size 1 supercell-y no-size)))
(set! geometry (append
                      (list (make block (center 0 0 0)(size a supercell-y no-size)(material (make dielectric (epsilon 11.56)))))
(list (make cylinder (center 0 (+ (/ h 2)  0.5 ))(radius rd)(height infinity) (material (make dielectric (epsilon 1)))))
(list (make cylinder (center 0 (+ (/ h 2)  -0.5 ))(radius rd)(height infinity) (material (make dielectric (epsilon 1)))))
(list (make cylinder (center 0.5 (+ (/ h 2)  0.5 ))(radius rd)(height infinity) (material (make dielectric (epsilon 1)))))
(list (make cylinder (center 0.5 (+ (/ h 2)  0.5 ))(radius rd)(height infinity) (material (make dielectric (epsilon 1)))))
(list (make cylinder (center 0 (+ (/ h 2)  1.5 ))(radius r)(height infinity) (material (make dielectric (epsilon 1)))))
(list (make cylinder (center 0 (+ (/ h 2)  -1.5 ))(radius r)(height infinity) (material (make dielectric (epsilon 1)))))
(list (make cylinder (center 0 (+ (/ h 2)  2.5 ))(radius r)(height infinity) (material (make dielectric (epsilon 1)))))
(list (make cylinder (center 0 (+ (/ h 2)  -2.5 ))(radius r)(height infinity) (material (make dielectric (epsilon 1)))))
(list (make cylinder (center 0 (+ (/ h 2)  3.5 ))(radius r)(height infinity) (material (make dielectric (epsilon 1)))))
(list (make cylinder (center 0 (+ (/ h 2)  -3.5 ))(radius r)(height infinity) (material (make dielectric (epsilon 1)))))))
 
(define Gamma (vector3 0 0 0))
(define K' (lattice->reciprocal (vector3  0.25 or 0.5 or 1 ? 0 0))) 
(set! k-points (interpolate 19 (list Gamma K')))
(define-param extra-bands 2)
(set! num-bands (+ 9 extra-bands))
(define-param res 32) 
(set! grid-size (vector3 res (* res supercell-y) 1))
(run-te (output-at-kpoint K' fix-hfield-phase output-hfield-z))

 
 
 
 

-------------------------------------------------------------------
新浪空间——与朋友开心分享网络新生活！(http://space.sina.com.cn/ )

_______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss