Hi everyone,
I am an undergraduate student and doing my final year project on Photonic
Crystal with Heterostructure. I am new user of MPB. i build an
increased-Index Waveguide with supercell.
the structure is like this:
...O...
...O...
...O...
...O...
...x...
...x...
...x...
...O...
...O...
...O...
...O...
it is air holes embedded in alumina. O and x are airholes with different
radius. the wave propagates along the gamma X direction. i want to find a
guided band in the bandgap. So i want to have a look at the H field of the
band24, 25. here i got a problem, i cant make MPB to put out the H field
at a k point of (0.1 0 0) (reciprocal). what i can get is the H field at k
point (0.5 0 0). Is that MPB can only output the field at edge of the
Brillouin zone? or something i have done wrong? Could somebody give me a
help? Thanks in advance.
=========
My CRL-file is like this:
=========
(define-param supercell-y 23)
(set! geometry-lattice (make lattice(size 1 supercell-y no-size)))
(define-param eps 8.9)
(define-param r1 0.4)
(define-param r2 0.3)
(set! default-material (make dielectric (epsilon eps)))
(set! geometry (list (make cylinder(center 0 0 0) (radius r1)(height
infinity)
(material (make dielectric (epsilon 1))))))
(set! geometry (geometric-objects-lattice-duplicates geometry))
(set! geometry (append geometry
(list (make block (center 0 0 0) (size 1 3 infinity)
(height infinity) (material(make dielectric(epsilon
8.9)))))))
(set! geometry (append geometry
(list (make cylinder(center 0 1 0) (radius r2)
(height infinity)(material (make dielectric(epsilon
1)))))))
(set! geometry (append geometry
(list (make cylinder(center 0 0 0) (radius r2)
(height infinity)(material (make dielectric(epsilon
1)))))))
(set! geometry (append geometry
(list (make cylinder(center 0 -1 0) (radius r2)
(height infinity)(material (make dielectric(epsilon
1)))))))
(define Gamma (vector3 0 0 0))
(define X (vector3 0.5 0 0))
(set! k-points (interpolate 4 (list Gamma X)))
(define-param extra-bands 5) ; number of extra bands to compute above the gap
(set! num-bands (+ supercell-y extra-bands))
(set-param! resolution 32)
;(use-output-directory)
(run-te (output-at-kpoint (vector3 0.1 0 0)
fix-hfield-phase output-hfield-z))
_______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
