Dear Steven and mpb users, I'm using MPB to simulate the bandgap of a silica-air photonic crystal. The pitch of the photonic crystal a=3um, air holes' radius r=1.35um, and the MPB CTL file is shown as follow. Basing on the code, gaps can be calculated with certain initial kz which means frequency of the modes.
My question is how to get the frequency(c/a)-wavevector k(2pi/a) plot (as
the attached Fig) with MPB? Is there an example?
Thank you ^_^
**************************************************
(define-param nsio2 1.457) ; the dielectric constant of SiO2
(define-param nair 1.0003) ; the dielectric constant of air
(define-param r 0.45) ; the hole radius
(set! geometry-lattice (make lattice (size 1 1 no-size)
(basis1 (/ (sqrt 3) 2) 0.5)
(basis2 (/ (sqrt 3) 2) -0.5)))
(define-param kz 5) ; pitch-3um wavelength-0.6um
(set! k-points (list (vector3 0 0 kz) ; Gamma
(vector3 0 0.5 kz) ; M
(vector3 (/ -3) (/ 3) kz) ; K
(vector3 0 0 kz))) ; Gamma
(set! k-points (interpolate 4 k-points))
(set-param! resolution 32)
(set-param! num-bands 10)
(set! default-material (make dielectric (index nsio2)))
(set! geometry (list (make cylinder
(center 0 0 0) (radius r) (height infinity)
(material (make dielectric (index nair))))))
(run)
**************************************************
Best regards,
He Yan
<<attachment: GAP.jpg>>
_______________________________________________ mpb-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
