Dear Steven and MPB users, I tried to reproduce the band-diagram results of a famous paper from *M*. * Notomi* : *Phys*. *Rev*. *B* 62, 10696 (*2000*) In this paper figure 7(b) shows a band-diagram of hexagonal pillars of GaAs (n = 3.6) rods of radius 0.35a in air (n = 1).
It was very easy to modify the mpb example tri-rods.ctl to do the same diagram. I tried but could not match the band diagram as given in the paper ! Appearance wise both results look very similar on first hand but closer inspection shows the difference. In the paper e.g. band 5 goes from frequency value 0.56 to *0.635* (as explicitly mentioned in explanation of Fig. 6 in the same paper) but in my calculated values using mpb i get them in the range 0.56 - *0.615* (the upper frequency value doesn't match). As far as i know the results of said paper were not calculated using mpb then is this discrepancy arising because of that ? If so, then how to make sure which result is more accurate ? I'm just putting my .ctl file below for all to have a look and suggest if i made any mistake in setting up .ctl file itself or i'm missing out on something which may improve my result. Any help in this will be much appreciated. ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; (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.35) (height infinity) (material (make dielectric (epsilon 12.96)))))) (set! k-points (list (vector3 (/ -3) (/ 3) 0) ; K (vector3 0 0.5 0) ; M (vector3 0 0 0) ; Gamma (vector3 (/ -3) (/ 3) 0) ; K )) (set! k-points (interpolate 4 k-points)) (set! resolution 64) (run-te) ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -- ................................................................... Thanks n Regards, manish <http://web.iitd.ernet.in/%7Ephz118136/> kumar India.
_______________________________________________ mpb-discuss mailing list mpb-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss