Consideramos la banda proyectada
sean los siguientes programas
(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 11.56)))
)
(make block
(center 0 0 0)
(size 0.5 1 0)
(material (make dielectric (epsilon 1.96)))
)
)
)
(set! k-points
(list
(vector3 0.0 0.0 0) ; Gamma
(vector3 0.0 1.0 0) ; M
)
)
(set! k-points (interpolate 20 k-points))
(set! resolution 32)
(run-tm)
(run-te)
(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 11.56)))
)
(make block
(center 0 0 0)
(size 0.5 1 0)
(material (make dielectric (epsilon 1.96)))
)
)
)
(set! k-points
(list
(vector3 0.5 0.0 0) ; Gamma
(vector3 0.5 1.0 0) ; M
)
)
(set! k-points (interpolate 20 k-points))
(set! resolution 32)
(run-tm)
(run-te)
mpb a02_0p00.ctl > a02_0p00.out mpb a02_0p50.ctl > a02_0p50.out grep tmfreqs a02_0p00.out > a02_tm_0p00.dat grep tmfreqs a02_0p50.out > a02_tm_0p50.dat
luego en gnuplot
plot 'a02_tm_0p00.dat' u 2:7 w l lc 1,'' u 2:8 w l lc 1,'' u 2:9 w l lc 1,'' u 2:10 w l lc 1, 'a02_tm_0p50.dat' u 2:7 w l lc 1,'' u 2:8 w l lc 1,'' u 2:9 w l lc 1,'' u 2:10 w l lc 1
