Hi,everyone
I am a new MPB user. When I try to simulate the ctl files in examples, I have a
problem.
When I am trying the tri-rods ctl file, I could not understand why the K
wavevector in irreducible Brillouin zone is (vector3 -1/3, 1/3, 0). I chceked
the book,"Molding the flow of light",it seems that the K point should be
(vector3 1/3, 1/3,0), and I can not find out where I was wrong. Please help me!
Thanks a lot!
Bing
; A triangular lattice of dielectric rods in air. (This structure has
; a band-gap for TM fields.) This file is used in the "Data Analysis
; Tutorial" section of the MPB manual.
(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 (make dielectric (epsilon 12))))))
(set! k-points (list (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 (output-at-kpoint (vector3 (/ -3) (/ 3) 0)
fix-efield-phase output-efield-z))
(run-te)
---------------------------------
抢注雅虎免费邮箱3.5G容量，20M附件！ _______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
