Hello,
When I run the following ctl file using MPBI vs MPB, I get different results
between the two, though I believe I should be getting the same output:
*******************************************************************
(set! default-material (make dielectric (index 1.0)) )
(define-param eps_air 1.0)
(define-param eps_sub 12)
(define-param thickness_air 10.0)
(define-param thickness_sub 1.0)
(set! geometry-lattice ( make lattice
(basis1 1 0 0)
(basis2 0 1 0)
(basis3 0 0 1)
(size 1 1 10)))
(set! num-bands 8)
(set! grid-size (vector3 16 16 160))
(set! k-points (list (vector3 0 0 0) ; Gamma
(vector3 0.5 0 0) ; X
(vector3 0.5 0.5 0) ; M
(vector3 0 0 0))) ; Gamma
(set! k-points (interpolate 4 k-points))
(set! geometry
(list
; PHOTONIC LAYER
(make cylinder
(material (make dielectric (epsilon eps_sub)))
(center 0 0 0)
(radius 2.0)
(height thickness_sub))
; AIR HOLES IN SUBSTRATE
(make cylinder
(material (make dielectric (epsilon eps_air)))
(center 0 0 0)
(radius 0.4)
(height thickness_air)
)
)
)
(run-even)
*****************************************************************
Any help that could be provided would be greatly appreciated.
Mark Herrera
_______________________________________________
mpb-discuss mailing list
mpb-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
