Re: mpb-discuss Digest, Vol 58, Issue 1

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: "[email protected]" 
<[email protected]>
To: [email protected]
Sent: Thu, December 2, 2010 6:00:01 PM
Subject: mpb-discuss Digest, Vol 58, Issue 1

Send mpb-discuss mailing list submissions to
    [email protected]

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
    [email protected]

You can reach the person managing the list at
    [email protected]

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? <[email protected]>
Subject: [MPB-discuss] photonic crystal book fig.10
To: [email protected]
Message-ID:
    <[email protected]>
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 <[email protected]>
Subject: Re: [MPB-discuss] photonic crystal book fig.10
To: [email protected]
Message-ID: <[email protected]>
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
> [email protected]
> 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? <[email protected]>
Subject: [MPB-discuss] photonic crystal book fig.10
To: [email protected]
Message-ID:
    <[email protected]>
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
[email protected]
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
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss