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.
tri-defect.ctl
Description: Binary data
_______________________________________________ mpb-discuss mailing list mpb-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss