Dear Questioner: Thr triangle lattice with air holes in dielectric material. I just applied that one and I got the same figure-10 in page 76 Regards
________________________________ From: "mpb-discuss-requ...@ab-initio.mit.edu" <mpb-discuss-requ...@ab-initio.mit.edu> To: mpb-discuss@ab-initio.mit.edu Sent: Thu, December 2, 2010 6:00:01 PM Subject: mpb-discuss Digest, Vol 58, Issue 1 Send mpb-discuss mailing list submissions to mpb-discuss@ab-initio.mit.edu To subscribe or unsubscribe via the World Wide Web, visit http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss or, via email, send a message with subject or body 'help' to mpb-discuss-requ...@ab-initio.mit.edu You can reach the person managing the list at mpb-discuss-ow...@ab-initio.mit.edu When replying, please edit your Subject line so it is more specific than "Re: Contents of mpb-discuss digest..." Today's Topics: 1. photonic crystal book fig.10 (?zg?r ?nder KARAKILIN?) 2. Re: photonic crystal book fig.10 (Frederic Brossard) 3. photonic crystal book fig.10 (?zg?r ?nder KARAKILIN?) ---------------------------------------------------------------------- Message: 1 Date: Wed, 1 Dec 2010 20:27:01 +0200 From: ?zg?r ?nder KARAKILIN? <okarakil...@pau.edu.tr> Subject: [MPB-discuss] photonic crystal book fig.10 To: mpb-discuss@ab-initio.mit.edu Message-ID: <145ed135a6c358283b5f74ea503b0395.squir...@eposta.pau.edu.tr> Content-Type: text/plain; charset=iso-8859-9 Hi everyone I am trying to page 76, figure 10 of "photonic crystals molding the flow of light" book. But I did not obtain same gap result. my ctl code is below. could you show me where is my mistake? thanks. ------------------------------ ;;phc book page 76 fig.10 (define-param epssub 13) ; dielectric substrate (define-param epsair 1) (define-param r 0.48) ; the radius of the holes (define-param supercell-y 31) ; the (odd) number of lateral supercell periods (set! default-material (make dielectric (epsilon epssub))) ; triangular lattice with vertical supercell: (set! geometry-lattice (make lattice (size 1 supercell-y no-size) (basis1 (/ (sqrt 3) 2) 0.5) (basis2 (/ (sqrt 3) 2) -0.5))) ;;;;;;;;;;;;;;;;;;;;;;;; Geometry (define-param a 1) (define-param d (* (/ (sqrt 3) 2) a)) (set! geometry (append (geometric-objects-lattice-duplicates (list (make cylinder (center 0 0) (radius r) (height infinity) (material (make dielectric (epsilon epsair))) ) (make cylinder (center (* .5 a) (* -1 d)) (radius r) (height infinity) (material (make dielectric (epsilon epsair))) ) ) 1 (* a (sqrt 3) a)) )) ; 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 9) ; the number of k points to interpolate (if only-K (set! k-points (list K)) (set! k-points (interpolate k-interp (list Gamma K')))) (define-param res 32) (set! grid-size (vector3 res (* res supercell-y) 1)) (define-param extra-bands 21) (set! num-bands (+ supercell-y extra-bands)) (set-param! mesh-size 7) (run-tm) ------------------------------ Message: 2 Date: Wed, 01 Dec 2010 18:37:51 +0000 From: Frederic Brossard <fs...@cam.ac.uk> Subject: Re: [MPB-discuss] photonic crystal book fig.10 To: mpb-discuss@ab-initio.mit.edu Message-ID: <4cf695ff.1090...@cam.ac.uk> Content-Type: text/plain; charset=ISO-8859-1; format=flowed Hi Ozgur, If you are referring to the second edition, then this corresponds to a structure without cavity or PhC waveguide, so no need for a supercell, (size 1 1) should do. Regards, Frederic Brossard On 01/12/2010 18:27, ?zg?r ?nder KARAKILIN? wrote: > Hi everyone > > I am trying to page 76, figure 10 of "photonic crystals molding the flow > of light" book. > But I did not obtain same gap result. > my ctl code is below. > > could you show me where is my mistake? > thanks. > > ------------------------------ > ;;phc book page 76 fig.10 > > (define-param epssub 13) ; dielectric substrate > (define-param epsair 1) > (define-param r 0.48) ; the radius of the holes > (define-param supercell-y 31) ; the (odd) number of lateral supercell periods > > (set! default-material (make dielectric (epsilon epssub))) > > ; triangular lattice with vertical supercell: > (set! geometry-lattice (make lattice (size 1 supercell-y no-size) > (basis1 (/ (sqrt 3) 2) 0.5) > (basis2 (/ (sqrt 3) 2) -0.5))) > > ;;;;;;;;;;;;;;;;;;;;;;;; Geometry > > (define-param a 1) > (define-param d (* (/ (sqrt 3) 2) a)) > > (set! geometry > (append > (geometric-objects-lattice-duplicates > (list > (make cylinder (center 0 0) (radius r) (height infinity) > (material (make dielectric (epsilon epsair))) > ) > (make cylinder (center (* .5 a) (* -1 d)) (radius r) (height >infinity) > (material (make dielectric (epsilon epsair))) > ) > ) > 1 (* a (sqrt 3) a)) > > )) > > > ; 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 9) ; the number of k points to interpolate > (if only-K > (set! k-points (list K)) > (set! k-points (interpolate k-interp (list Gamma K')))) > > > (define-param res 32) > (set! grid-size (vector3 res (* res supercell-y) 1)) > > (define-param extra-bands 21) > > (set! num-bands (+ supercell-y extra-bands)) > > (set-param! mesh-size 7) > > (run-tm) > > > > _______________________________________________ > mpb-discuss mailing list > mpb-discuss@ab-initio.mit.edu > http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss ------------------------------ Message: 3 Date: Thu, 2 Dec 2010 00:12:05 +0200 From: ?zg?r ?nder KARAKILIN? <okarakil...@pau.edu.tr> Subject: [MPB-discuss] photonic crystal book fig.10 To: mpb-discuss@ab-initio.mit.edu Message-ID: <7d6321f7411bcc74f0e6eebcc095d06b.squir...@eposta.pau.edu.tr> Content-Type: text/plain; charset=iso-8859-9 thanks for your reply Frederic, I am referring to second edition book. I forgot to say that I intentionally used supercell. I know that I do not see same graphs but it is expected to give same frequency gap region. If I am sure our ctl file used supercell, I am going to calculate cavity and other structure with modified ctl file. I used similar ctl file which is modified with lattice vector for square lattice in example figure 2 in the same book page 194. I obtained same graphs with supercell. But for hexagonal structure, although modifying to lattice vector, I do not obtain correct band gap frequency. if you show me my mistake, I am very appreciate to you. Re: [MPB-discuss] photonic crystal book fig.10 Frederic Brossard Wed, 01 Dec 2010 10:39:58 -0800 Hi Ozgur, If you are referring to the second edition, then this corresponds to a structure without cavity or PhC waveguide, so no need for a supercell, (size 1 1) should do. Regards, Frederic Brossard On 01/12/2010 18:27, ?zg?r ?nder KARAKILIN? wrote: Hi everyone I am trying to page 76, figure 10 of "photonic crystals molding the flow of light" book. But I did not obtain same gap result. my ctl code is below. could you show me where is my mistake? thanks. ------------------------------ ;;phc book page 76 fig.10 (define-param epssub 13) ; dielectric substrate (define-param epsair 1) (define-param r 0.48) ; the radius of the holes (define-param supercell-y 31) ; the (odd) number of lateral supercell periods (set! default-material (make dielectric (epsilon epssub))) ; triangular lattice with vertical supercell: (set! geometry-lattice (make lattice (size 1 supercell-y no-size) (basis1 (/ (sqrt 3) 2) 0.5) (basis2 (/ (sqrt 3) 2) -0.5))) ;;;;;;;;;;;;;;;;;;;;;;;; Geometry (define-param a 1) (define-param d (* (/ (sqrt 3) 2) a)) (set! geometry (append (geometric-objects-lattice-duplicates (list (make cylinder (center 0 0) (radius r) (height infinity) (material (make dielectric (epsilon epsair))) ) (make cylinder (center (* .5 a) (* -1 d)) (radius r) (height infinity) (material (make dielectric (epsilon epsair))) ) ) 1 (* a (sqrt 3) a)) )) ; 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 9) ; the number of k points to interpolate (if only-K (set! k-points (list K)) (set! k-points (interpolate k-interp (list Gamma K')))) (define-param res 32) (set! grid-size (vector3 res (* res supercell-y) 1)) (define-param extra-bands 21) (set! num-bands (+ supercell-y extra-bands)) (set-param! mesh-size 7) (run-tm) ------------------------------ _______________________________________________ mpb-discuss mailing list mpb-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss End of mpb-discuss Digest, Vol 58, Issue 1 ****************************************** To get in shape as in page 76 you refer to tutorials that avialable in MIT MPB site
_______________________________________________ mpb-discuss mailing list mpb-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss