Hi Karel, I don't understand what you mean by a slab with a radius in your description below. Could you elaborate? Then I could go through your .ctl file.
cheers Manoj Karel Van Acoleyen wrote: > 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 al je foto's met je vrienden! Maak je eigen online plek! Windows > Live Spaces <http://spaces.live.com> > > > ------------------------------------------------------------------------ > > _______________________________________________ > meep-discuss mailing list > [email protected] > http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss _______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
