1D Primera zona de Brillouin

(set! num-bands 10)

(set! geometry-lattice 
 (make lattice 
  (size 1 no-size no-size)
 )
)
         
(set! geometry
  (list
   (make block
    (center 0 0 0)
    (size 1.0 1.0 0)
    (material (make dielectric (epsilon 5.52)))
  )
   (make block
    (center 0 0 0)
    (size 0.5 1 0)
    (material (make dielectric (epsilon 2.1316)))
  )
 )
)

(set! k-points 
 (list 
  (vector3 0.0 0 0)          ; Gamma
  (vector3 0.5 0 0)        ; M
 )
)
(set! k-points (interpolate 20 k-points))

(set! resolution 32)

(run-tm)
(run-te)

mpb a01.ctl > a01.out
grep tmfreqs a01.out > a01.tm.dat
grep tefreqs a01.out > a01.te.dat
gnuplot
plot 'a01.tm.dat' u 6:7 w l lc 1,'' u 6:8 w l lc 1,'' u 6:9 w l lc 1,'' u 6:10 w l lc 1

11