hex.PhC k-points

Özgür Önder KARAKILINÇ <[email protected]> Mon, 25 Apr 2011 16:20:32 +0300
Newsgroups gmane.comp.science.photonic-bands
Message-ID <[email protected]>
Dear mpb users

I examine the tri-rods.ctl example in mpb wiki page.
This real space structure is triangular and corresponded reciprocal
structure is hexagonal, isn&#8217;t it?
When I change the real space basis for hexagonal structure as the
following code, what will the k-points become?
Does the k-point remain still same?
best...
Onder



(define-param npoly 1.46)
(define-param epspoly (* npoly npoly))
(define-param epsair 1)

(define-param a 1)
(define-param r (* 0.2 a))

(set! num-bands 8)

(set! default-material (make dielectric (epsilon epspoly)))

(set! geometry-lattice (make lattice (size 1 1 no-size)
                         (basis1 0.5 (/ (sqrt 3) 2))
                         (basis2 1 0)))
(set! geometry (list (make cylinder
                       (center 0 0 0) (radius r) (height infinity)
                       (material (make dielectric (epsilon epsair))))))

(set! k-points (list (vector3 0 0 0)          ; Gamma
                     (vector3 0 0.5 0)        ; M
                     (vector3 (/ -3) (/ 3) 0) ; K
                     (vector3 0 0 0)))        ; Gamma
(set! k-points (interpolate 4 k-points))

(set! resolution 32)

(run-tm (output-at-kpoint (vector3 (/ -3) (/ 3) 0)
                          fix-efield-phase output-efield-z))
(run-te)