Dear All,
Is there any one who used mpb for simulating the confined TM bands of 2D
planar waveguide? I wanted to do it but I couldn't see the second and the
third band in the frequency vs. wavevector graph. I include the ctl file in
the mail. Thank you in advance
Kadir.
(define-param eps 12) ; dielectric constant of
waveguide
; The cell
dimensions
(define-param sy 9) ; size of cell in y direction (perpendicular to
wvg.)
(define-param sx 1) ; size of cell in x direction
(set! geometry-lattice (make lattice (size sx sy no-size)))
(set! geometry
(list
(make block (center 0 0) (size 1 0.2 infinity) (material (make
dielectric (epsilon eps))))
)
)
(set-param! resolution 64)
(define Gamma (vector3 0 0 0))
(define K (lattice->reciprocal (vector3 0.5 0 0))) ; edge of Brillouin zone.
(define M (lattice->reciprocal (vector3 0.5 0.5 0))) ; edge of Brillouin
zone.
(set! k-points (interpolate 19 (list Gamma K ;M Gamma
)))
(define-param extra-bands 5) ; number of extra bands to compute above the
gap
(set! num-bands (+ sy extra-bands))
(run-tm fix-efield-phase output-efield-z)
_______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
