Dear mpb users,
I am trying to calculate band solve of following ctl file structure with
supercell.
Actually, frequency defination is different between meep and mpb.
Frequency is equivalent a/lambda (so units of a/c)
in meep but in contrast, frequency is unit of (c/a) in mpb.
So I do not use directly frequency value of mpb graphs in to the meep sim.
is it correct thinking so far?
I want to use frequency value of mpb calculation directly in to the meep sim.
As I undestand from the manual page, it needs to transform k points
reciprocal->cartesian
But at the beginning in mpb calculation, can I use
(define Gamma (vector3 0 0 0))
(define K' (lattice->cartesian (vector3 0.5 0 0)))
for compatible with meep. I am not sure this usage.
How can I modify this defination for compatible with meep.
I am waiting for your suggestion. Could you check my ctl file and show me
if there are wrong thinking.
Thanks in advance.
Ozgur Onder Karakilinc
Electrical-Electronics Eng.Dept.
20070 Kinikli-Denizli TURKEY
Tel: (+90 258) 296 3213
;;; hex.PhC with defect
(define-param epspoly 2.25)
(define-param epsair 1)
(define-param r 0.2)
;;;;;;;;;;; lattice
(define-param supercell-y 31)
(set! default-material (make dielectric (epsilon epspoly)))
(set! geometry-lattice (make lattice
(basis1 (/ (sqrt 3) 2) 0.5)
(basis2 (/ (sqrt 3) 2) -0.5)
(size 1 supercell-y 1)))
;;;;;;;;;;;;;;;;;;;;;;;; Geometry
(define-param a 1)
(define-param d (* (/ (sqrt 3) 2) a)) ; define d cell size
(set! geometry
(append
(geometric-objects-lattice-duplicates
(list
(make cylinder (center 0 0 0) (radius r) (height infinity)
(material (make dielectric (epsilon epsair)))
)
(make cylinder (center (* .5 a) (* -1 d) 0) (radius r)
(height
infinity) (material (make dielectric (epsilon epsair)))
)
)
1 (* a (sqrt 3) a))
))
(set! geometry (append geometry
(list (make cylinder (center 0 0 0)
(radius r) (height infinity)
(material (make dielectric (epsilon
epspoly))))
)
)
)
;;;;;;; how can I modify this defination for compatible with meep.
(define Gamma (vector3 0 0 0))
(define K' (lattice->reciprocal (vector3 0.5 0 0))) ; edge of Brillouinzone.
(set! k-points (interpolate 4 (list Gamma K')))
(define-param extra-bands 5)
(set! num-bands (+ supercell-y extra-bands))
(define-param res 32)
(set! grid-size (vector3 res (* res supercell-y) 1))
; for TM
(run-tm)
; for TE
;(run-te)
;;;;;;;;;;;;;;;;;;;;;;;;
_______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
