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 mpb-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss