magneto-optical material

Dear Steven and Others:
  I need to calculate the band diagram for the triangular lattice of
dielectric rods made of magneto-optical material in air,  but the diagram I
get makes no difference with lattice of ordinary dielectric rods.
  I have compiled MPB with "--with-hermitian-eps" flag, can some give me a
test case or figure out the errors in my input file.

(define-param CuO_eps 12)
(define-param CuO_gamma 0.4)
(define CuO (make dielectric-anisotropic
     (epsilon-diag CuO_eps CuO_eps CuO_eps)
     (epsilon-offdiag 0+0.4i 0 0)
    )
)
.
(set! num-bands 8)
(set! geometry-lattice (make lattice (size 1 1 no-size)
                         (basis1 (/ (sqrt 3) 2) 0.5)
                         (basis2 (/ (sqrt 3) 2) -0.5)))
(set! geometry (list (make cylinder
                       (center 0 0 0) (radius 0.2) (height infinity)
                       (material CuO))))
(set! k-points (list (vector3 0 -0.5 0)       ;-M
                     (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 4 k-points))
(set! resolution 32)
(run-tm)
(run-te)

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