> Hi, the ctl file that you need is as follow :(define-param h 0.6) ; the thickness of the slab(define-param eps 12) ; the dielectric constant of the slab(define-param loweps 1.0) ; the dielectric constant of the substrate(define-param r 0.3) ; the radius of the holes(define-param supercell-h 7) ; height of the supercell(define-param supercell-y 7) ; the (odd) number of lateral supercell periods; triangular lattice with vertical supercell:(set! geometry-lattice (make lattice (size 1 supercell-y supercell-h) (basis1 (/ (sqrt 3) 2) > 0.5) (basis2 (/ (sqrt 3) 2) -0.5)))(set! geometry (list (make block (material (make dielectric (epsilon loweps))) (center 0 0 0) (size infinity infinity supercell-h)) (make block (material (make dielectric (epsilon eps))) (center 0 0 0) (size infinity infinity h)) (make cylinder (material air) (center 0) (radius r) (height supercell- h))))(set! geometry (append ; duplicate the bulk crystal rods > over the supercell: (geometric-objects-lattice-duplicates geometry 1 1 supercell-h) ; add a rod of air, to erase a row of rods and form a waveguide: (list (make cylinder (center 0) (radius r) (height supercell-h) (material (make dielectric (epsilon eps))))) )); 1st Brillouin zone of a triangular lattice:(define Gamma (vector3 0 0 0))(define K' (lattice->reciprocal (vector3 0.5 0 0))) ; edge of Brillouin zone.(define-param only-K false) ; run with only-K=true to only do this k- point(define-param k-interp 4) ; the number of k points to interpolate(if only-K (set! k-points (list K)) (set! k-points (interpolate > k-interp (list Gamma K'))))(set-param! resolution (vector3 32 20 20)) (set-param! mesh-size 7)(set-param! num-bands 11); Run even and odd bands, outputting fields only at the K point:(run-zeven (output-at-kpoint K' output-hfield-z))(display-eigensolver-stats) Good Luck--- On Fri, 3/18/11,
Hi, Mahmoud, I have tried with the ctl file attached by you in this discussion. I just copy and paste and run your ctl file, however, when i plot the output (kx vs freq bands), the result is not the same as reported by steven in his book (Molding the Flow of Light, chapter 8, fig.8). May i know any mistake in my simulation? Thank you regards, DC _______________________________________________ mpb-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
