Hi Steven, MPB users, Still with the same H1 cavity, one can find another issue, this time concerning the E field components. Attached are pictures of Hz, Ex, Ey. Hz gives the correct field profile whereas Ex, Ey are erroneous if compared with published results. The same problem can be found with other types of cavities, also for different cell size, resolution, k-point, phase fixed or not, even 2D calculation (hence z-even or TE)...
Again, any help would be appreciated.
Thanks,
Fred
; A triangular lattice of air holes with H1 defect
(define-param sx 5)
(define-param sy 5)
(define-param sz 2)
(define-param resxy 16)
(define-param resz 10)
(define-param bands 50)
(define r 0.36)
(define h 0.75)
(define eps 8.4)
(set! mesh-size 6)
(set! num-bands bands)
(define Gamma (vector3 0))
(define M0 (vector3 0 0.5 0))
(define K0 (vector3 (/ -3) (/ 3) 0))
(set! geometry-lattice (make lattice (size sx sy sz)
(basis1 (/ (sqrt 3) 2) 0.5)
(basis2 (/ (sqrt 3) 2) -0.5)))
(set! geometry (list (make block (material air)
(center 0)
(size infinity infinity infinity))
(make block (material (make dielectric (epsilon eps)))
(center 0)
(size infinity infinity h))))
(set! geometry (append geometry (geometric-objects-lattice-duplicates (list
(make cylinder
(center 0) (radius r) (height infinity)
(material air))))))
(set! geometry (append geometry
(list (make cylinder
(center 0) (radius r) (height h)
(material (make dielectric (epsilon eps)))))))
(set! k-points (list M0))
(set! resolution (vector3 resxy resxy resz))
(set! tolerance 1e-5)
(set! filename-prefix "3DH1-")
(run-zeven)
(fix-efield-phase bands)
(output-efield-x bands)
(output-efield-y bands)
(fix-hfield-phase bands)
(output-hfield-z bands)
(output-dpwr bands)
--
Frederic Brossard
Hitachi Cambridge Laboratory
Cavendish Laboratory
J J Thomson Avenue
Cambridge CB3 0HE
<<attachment: 3DH1hz.jpg>>
<<attachment: 3DH1Ex.jpg>>
<<attachment: 3DH1Ey.jpg>>
_______________________________________________ mpb-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
