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
