Hi, Dr. Steven suggest mpb for computing the band structures, also your structure consist a hole slab with square lattice, i don't know why you set your geometry to (10 10 10), you have to set it only as (1 1 10). 2. in the case when you want to introduce a defect (point defect or linear defect), you have to set your geometry to (1 6 10). i can send you a ctl.file example if you want. hope this help you with regards
--- On Wed, 7/7/10, [email protected] <[email protected]> wrote: From: [email protected] <[email protected]> Subject: [Meep-discuss] simple square slab To: [email protected] Date: Wednesday, July 7, 2010, 1:04 AM Dear Dr. Steven and meep users, I want to simulate a slab structure (without defect for the moment) with square lattice. I try to obtain the band diagram (I am particularly interested of the part above the light line). ;my file (define-param r 0.3) (define-param eps 11.36) (define-param loweps 1.0) (define-param dpml 1) (define-param supercell-h 10) ;dimensions of my slab: (define-param h 0.31) (define-param sx 5) (define-param sy 5) (set! eps-averaging? false) (set! default-material (make dielectric (epsilon loweps))) (set! geometry-lattice (make lattice (size 10 10 supercell-h))) (set! geometry (list (make block (material (make dielectric (epsilon eps))) (center 0) (size sx sy h)) (make cylinder (center 0 0 0) (radius r) (height h) (material air)))) (set! geometry (geometric-objects-lattice-duplicates (list (make cylinder (center 0 0 0) (radius r) (height h) (material air))) 1 1 supercell-h)) ;(set! pml-layers (list (make pml (direction Z) (thickness dpml)))) (set! pml-layers (list (make pml (thickness dpml)))) (define-param fcen 0.29) (define-param df 0.1) (set! sources (list (make source (src (make gaussian-src (frequency fcen)(fwidth df))) (component Ez) (center 0 0 (+ dpml (* -0.5 supercell-h))) ))) (set! resolution 8) (define-param k-interp 9) (run-k-points 150 (interpolate k-interp (list (vector3 0 0 0) (vector3 0.5 0 0)))) When I put pml layers in the three directions, the result is zero (w=0). When I calculate with only a pml layer in z direction, I obtain some value of w but if I executed h5topng -0z0 I got only black long rectangle. I don't find what is wrong in my file. I don't understand why it's not working when I put pml layers in x and y direction even far from my structure. thanks Regards, Elise _______________________________________________ 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
