comparison of Guided Mode Expansion, MPB

Hello Steven and MPB Users,

I am trying to simulate the slab photonic crystal structure in your Phys
Rev. B 60 (1999) paper using two methods, MPB and a code by L.C.
Andreani called GME
(http://fisicavolta.unipv.it/dipartimento/ricerca/Fotonici/Index.htm).
Attached is my MPB .ctl file.  The band structure plots can be found at:
http://www.chriskang.net/incoming/r0.4_GME.pdf
http://www.chriskang.net/incoming/r0.4_MPB.pdf

My question pertains to the appearance of the bands in the light cone.
While the gap seems to remain intact under the light cone, there seems
to be some states that MPB is calculating that do not appear on the GME
band plot.  Is this merely a difference in how the modes are calculated,
or is there something wrong with how my structure is set up in MPB?  Any
help would be appreciated.

Best regards,
Chris

_______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss

; try a pc slab in air

(set! num-bands 12)

(set! k-points (list (vector3 0 0 0) 	 ; Gamma
		     (vector3 0 0.5 0) 	 ; M
		     (vector3 (/ -3) (/ 3) 0) ; K
		     (vector3 0 0 0) ))  ; Gamma

(set! k-points (interpolate 10 k-points))

(set! default-material (make dielectric (epsilon 12)))

(set! geometry (list (make cylinder (center 0 0 0)(radius 0.3)(height 0.6)(material air))
		     (make block (center 0 0 1.4)(size infinity infinity 2.2)(material air))
 	             (make block (center 0 0 -1.4)(size infinity infinity 2.2)(material air))      
		)
)

(set! geometry-lattice (make lattice (size 1 1 5)(basis1 (/ (sqrt 3) 2) 0.5)
                                                 (basis2 (/ (sqrt 3) 2) -0.5)         ) )

(set! resolution 32)


(run-te)
(run-tm)