Dear MPB user and St.V.Johnson,
I wish to calculate the band diagram of a overlapping titania shell
hollow model ( inverse opal). I wonder my calculation is correct since
the result is different from the transmission spectra:
(set! default-material (make dielectric (epsilon 1)))
(set! geometry-lattice (make lattice
(basis-size (sqrt 0.5) (sqrt 0.5) (sqrt 0.5))
(basis1 0 1 1)
(basis2 1 0 1)
(basis3 1 1 0)))
(set! k-points (interpolate 4 (list
(vector3 0 0.5 0.5) ; X
(vector3 0 0.625 0.375) ; U
(vector3 0 0.5 0) ; L
(vector3 0 0 0) ; Gamma
(vector3 0 0.5 0.5) ; X
(vector3 0.25 0.75 0.5) ; W
(vector3 0.375 0.75 0.375)))) ; K
;(define-param rbig 0.380953778) ; sphere of shell
;(define-param rsmall 0.35355339) ; radius of the air spheres
(define diel1 (make dielectric (epsilon 1)))
(define diel2 (make dielectric (epsilon 6.25)))
(set! geometry (list (make sphere (center 0) (radius 0.380953778) ; make
titania shell
(material diel2))
(make sphere (center 0) (radius 0.35355339)
; make air sphere
(material diel1))
))
(set-param! resolution 16)
(set-param! mesh-size 10)
(set-param! num-bands 20)
(run)
thanks in advance
_______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss