frequency conversion

Özgür Önder KARAKILINÇ <[email protected]>
Newsgroups gmane.comp.science.photonic-bands
Message-ID <[email protected]>
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)

;;;;;;;;;;;;;;;;;;;;;;;;
lmpx.com only provides a reader for public news (NNTP) servers. It is not affiliated with the servers or forums shown here and is not responsible for the content of articles, which is written by their respective authors.