Hello,
I'm currenlty working on photonic crystal lenses. I am interested in using MPB
to develop a crystal structure, then use meep to do field analysis and examine
how the fields are focused by the crystal lense.
I have a crystal that has properties I want, now I'm trying to get meep to use
the geometry design in MPB. I would like to make a crystal slab that's in the
center of my simulation, then put a source and 'detector' in the air on either
side. I have a geometry almost working, but meep isn't setting up the crystal
right. In the code below it makes a square lattice, not a triangle lattice.
I'm pretty sure this is because of the basis vectors for the geometry-lattice
not being set, but when I leave those on, I can't get my source to be a plane
wave in the x,y, directions.
Can you set the geometry-lattice, then define a geometry, then change the
lattice and have it not affect the previosuly define geometry? This doesn't
sound like it would work.
Any sugestions or thoughts?
(set! geometry-lattice (make lattice
(size sx sy no-size))) ; uncomment for full simulation
space
;(size 1 1 no-size) ;
uncomment for bandgap analysis
;(basis1 (/ (sqrt 3) 2) 0.5)
;(basis2 (/ (sqrt 3) 2) -0.5)))
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;; Set up the simulation space. First define the photonic crystal by
repicated a unit cell over the entire
;;; simulation domain. Second overwrite the left and right edges of the
domain with the dielectric block.
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
(set! default-material (make dielectric (epsilon eps)))
(set! geometry (append
(geometric-objects-lattice-duplicates (list ;uncomment for
full simulation space
(make cylinder (center 0 0 0) (radius (/ r a))
(height infinity)
(e1 (/ (sqrt 3) 2) 0.5 0) (e2 (/ (sqrt 3) 2) -0.5)
(e3 0 0 1)
(material air))))))
(set! geometry (append geometry
(list
(make block (center (- (/ SlabWidth -2) (/ sx 2) 0)) (size (- sx (/
SlabWidth 2)) sy 0)
(e1 1 0 0) (e2 0 1 0) (e3 0 0 1)(material (make dielectric (epsilon air))))
(make block (center (+ (/ SlabWidth 2) (/ sx 2) 0)) (size (- sx (/
SlabWidth 2)) sy 0)
(e1 1 0 0) (e2 0 1 0) (e3 0 0 1)(material (make dielectric (epsilon air))))
)))
Nathan R. Huntoon, Ph.D.
Director Innovation Gymnasium
Caruth Institute for Engineering Education
[email protected]
214-768-1402 (office)
Sent from my HP TouchSmart tm2 running Linux Mint 11
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
