Hi Steven, MPB users, Below is a simple test to view the z-component of the H field and magnetic energy density in a triangular lattice of air-holes. I have attached the pics for the first and second band. This looks a bit odd, should the field antinodes not be located in the dielectric for the first band? it almost looks like the opposite of what we would expect (first and second band). I am sure this must have been discussed elsewhere but couldn't find this particular topic.
Thanks,
Fred.
(define-param r 0.3)
(define K (vector3 (/ -3) (/ 3) 0))
(set! num-bands 2)
(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 r) (height infinity)
(material air))))
(set! default-material (make dielectric (epsilon 11.09)))
(define Brill (list K))
(set! k-points Brill)
(set! resolution 32)
(run-te fix-hfield-phase output-hfield-z output-hpwr)
--
Frederic Brossard
Hitachi Cambridge Laboratory
Cavendish Laboratory
J J Thomson Avenue
Cambridge CB3 0HE
h.k01.b01.z.te.png
Description: h.k01.b01.z.te.png
h.k01.b02.z.te.png
Description: h.k01.b02.z.te.png
hpwr.k01.b01.te.png
Description: hpwr.k01.b01.te.png
hpwr.k01.b02.te.png
Description: hpwr.k01.b02.te.png
_______________________________________________ mpb-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
