Dear meep users,
A few weeks ago I sent a mail to try to simulate waveguide modes of line
defects in photonic crystal slabs. I now have a slab of thickness 0.5a, radius
0.3a and just one missing row of holes. I fixed the previous mentioned errors,
but the band diagram stays the same. My mode is lying between 0.2 and 0.3
normalized f, while it should be between 0.25 and 0.35 as mentioned in articles
and that is where my band gap is. I added my ctl file again.
Thanks in advance,
Karel
------------------------------------------------------------------
(define silicon (make dielectric (index 3.5)))
(define silica (make dielectric (index 1.45)))
(define SiN (make dielectric (index 2)))
(define-param h 0.5)
(define-param r 0.3)
(define-param h_air 2)
(define-param h_Si 0.255)
(define-param h_SiO2 0.61576)
(define-param n_bragg 9)
(define-param h_SiN 0.44643)
(set! geometry-lattice (make lattice (size 1 (* 13 (sqrt 3)) 3))) ;(+ (* h_air
2) h (* 2 n_bragg (+ h_SiN h_SiO2)))))
(set! default-material air) ; everything is air
(define bragg_top (list
(make block (center 0 0 (+ (/ h 2) (/ h_SiO2 2))) (size 25 25 h_SiO2) (material
silica))
(make block (center 0 0 (+ (/ h 2) h_SiO2 (/ h_SiN 2))) (size 25 25 h_SiN)
(material SiN))))
(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)))
;(geometric-objects-duplicates (vector3 0 0 (+ h_SiN h_SiO2)) 0 n_bragg
bragg_top)
))
(set! symmetries (list (make mirror-sym (direction Z) (phase -1)) (make
mirror-sym (direction Y) (phase -1))))
(set! pml-layers (list (make pml (direction Y) (thickness 1)) (make pml
(direction Z) (thickness 1))))
(define-param fcen 0.30) ; pulse center frequency
(define-param df 0.20) ; 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 0 0 0)
)
)
)
(run-k-points 300 (interpolate k-interp (list (vector3 0 0 0) (vector3 0.5 0 0)
)))
_________________________________________________________________
Deel je gedachten direct met je vrienden
http://messenger.live.com
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
