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
