Dear Professor and MPB users
When I simulate photonic crystal, I have got a problem. It's so strange that I can't explain it.
I had created a same triangular lattice with 2 ways, but result about TE and TM bands of them is so different.
Case 1: Fist, I make a triangular lattice
with "size (1 1 no-size)":
(define-param r 0.36)
(define-param eps 13)
(define Si (make dielectric (epsilon eps)))
(set! geometry-lattice (make lattice (size 1 1 no-size)
(basis1 (/ (sqrt 3) 2) 0.5)
(basis2 (/ (sqrt 3) 2) -0.5)))
(set! default-material (make dielectric (epsilon 13)))
(set! geometry
(list (make cylinder (material air)
(center 0 0 0) (radius r) (height infinity))))
(set! geometry (geometric-objects-lattice-duplicates geometry))
(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-param! resolution 16)
(set-param! num-bands 10)
(run-tm (output-at-kpoint (vector3 (/ -3) (/ 3) 0)
fix-efield-phase output-efield-z))
(run-te)
Then, I run the calculation following the tutorial.
...unix% mpb-data -r -m 3 -n 32 epsilon.h5....
Case 2: I make a triangular lattice with (size 3 3 no-size). I only change this parameter in case 1.
Then I run the calculation with:
...unix% mpb-data -r -n 32 epsilon.h5 (without "-m 3").
Two case take same lattice, but the result is so different. Can you explain this? Please help me. I just learn MPB for 2 months, so my knowledge is not good.
I'm so excited that I can't wait any more.
I'm looking forward to hearing from you.
Thank you very much!
sincerly,
Hoa Nguyen An.