Dear meep users,
I am trying to calculate the band diagram for waveguide modes in line defects
in photonic crystals. I looked to find some examplary ctl files, but did not
find one, so I made one of my own. Now when I run this I get a waveguide band
diagram, but if I only want the even mode and specify the mirror-symmetry plane
around the y axis, harminv does not produce any results, although my source is
located at y = 0 so that is also symmetric. If I furthermore do not specify
this I get a band diagram but it seems like the bandgap has shifted (the gap
does not corresond with the bandgap, and I used another program and that gives
the expected better results, but I am trying to compare programs). The
structure is just one missing row of holes in an air slab. My ctl files is
given below. Does anyone has some advice?
Thanks in advance,
Kind regards,
Karel Van Acoleyen
*********************
ctl file:
(define silicon (make dielectric (index 3.5)))
(define silica (make dielectric (index 1.45)))
(define SiN (make dielectric (index 2)))
(define-param h 0.6)
(define-param r 0.3)
(define-param h_air 2)
(define-param h_Si 0.255)
(define-param h_SiO2 0.61576)
(define-param h_SiN 0.44643)
(set! geometry-lattice (make lattice (size 1 (* 13 (sqrt 3)) 3)))
(set! default-material air) ; everything is air
(set! geometry (append (list (make block (center 0 0 0) (size infinity infinity
h) (material silicon)))
(geometric-object-duplicates (vector3 0 (sqrt 3) 0) 0 5 (make cylinder (center
0 (sqrt 3) 0) (radius r) (height h) (material air)))
(geometric-object-duplicates (vector3 0 (sqrt 3) 0) 0 5 (make cylinder (center
-0.5 (/ (sqrt 3) 2) 0) (radius r) (height h) (material air)))
(geometric-object-duplicates (vector3 0 (sqrt 3) 0) 0 5 (make cylinder (center
0.5 (/ (sqrt 3) 2) 0) (radius r) (height h) (material air)))
(geometric-object-duplicates (vector3 0 (* -1 (sqrt 3)) 0) 0 5 (make cylinder
(center 0 (* -1 (sqrt 3)) 0) (radius r) (height h) (material air)))
(geometric-object-duplicates (vector3 0 (* -1 (sqrt 3)) 0) 0 5 (make cylinder
(center -0.5 (/ (* -1 (sqrt 3)) 2) 0) (radius r) (height h) (material air)))
(geometric-object-duplicates (vector3 0 (* -1 (sqrt 3)) 0) 0 5 (make cylinder
(center 0.5 (/ (* -1 (sqrt 3)) 2) 0) (radius r) (height h) (material air)))
))
(set! symmetries (list (make mirror-sym (direction Z) (phase -1))))
(set! symmetries (list (make mirror-sym (direction Y) (phase 1))))
(set! pml-layers (list (make pml (direction Y) (thickness 1))))
(set! pml-layers (list (make pml (direction Z) (thickness 1))))
(define-param fcen 0.25) ; pulse center frequency
(define-param df 0.15) ; pulse freq. width
(define randomnumber 0.13)
(define-param k-interp 10)
(set-param! resolution 16)
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Hz) (center randomnumber 0 0)
)
)
)
(run-k-points 300 (interpolate k-interp (list (vector3 0 0 0) (vector3 0.5 0 0)
)))
_________________________________________________________________
De snelle weg naar een praatje
http://messenger.live.com
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
