help

MPB users:

I am trying to model a photonic crystal slab structure (air in dielectric) with one defect (circular hole of different size) in the center.  Slab is surrounded by water cladding regions on top and bottom.  I am having a hard time adding a single defect in the center of the lattice pattern.  Perhaps I am missing a command or I have something incorrect.  The code that I have does work, however I have looked at the geometry of the structure and it does not appear to have a defect in the center.  Another concern is I do not want the defect to extend into the water clad regions.  The defect could be filled with water or air either one.   I am new to MPB and would greatly appreciate any help you experts could provide.  I have attached my code below.

Thank You

(define-param h 0.9)
(define-param eps 5.76)
(define-param loweps 1.77)
(define-param r 0.28)
(define-param supercell-h 12)
(set! geometry-lattice (make lattice (size 1 1 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.25 supercell-h))
(size infinity infinity (* 0.5 supercell-h)))
(make block (material (make dielectric (epsilon loweps)))
(center 0 0 (* -0.25 supercell-h))
(size infinity infinity (* 0.5 supercell-h)))
(make block (material (make dielectric (epsilon eps)))
(center 0) (size infinity infinity h))
(make cylinder (material air) (center 0) (radius r) (height h))))
(set! geometry (append geometry
                (list (make cylinder (center 0 0 0)
                (radius 0.15) (height h) 
                (material air)))))
(define Gamma (vector3 0 0 0))
(define M (vector3 0 0.5 0))
(define K (vector3 (/ -3) (/ 3) 0))
(define-param k-interp 4)
(set! k-points (interpolate k-interp (list Gamma M K Gamma)))
(set-param! resolution (vector3 24 24 24))
(set-param! num-bands 9)
(run-zeven)
(run-zodd)
(display-eigensolver-stats)


Randy Tompkins

Graduate Research Assistant

West Virginia University
Physics Dept.

Phone:  (304)293-3422 x1436

Home:  (304)-906-9821

Fax: (304)-293-5732


_________________________________________________________________
Insert movie times and more without leaving Hotmail®. 
http://windowslive.com/Tutorial/Hotmail/QuickAdd?ocid=TXT_TAGLM_WL_HM_Tutorial_QuickAdd_062009

_______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss