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’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)