hi, this is what i use. as steven notes in the manual (which you
should have read) the fcc unit cell has an edge length of sqrt(2)
times the diameter of your spheres, so you need to take that factor
into account if you convert the freqs output by mpb to real frequency
units.
anyway here is the ctl file.. cheers:
(set! num-bands 10)
(set! default-material (make dielectric (index 3.9)))
(set! resolution 16)
(set! mesh-size 5)
(set! geometry-lattice
(make lattice
(basis1 0 1 1)
(basis2 1 0 1)
(basis3 1 1 0)
(basis-size (sqrt 0.5) (sqrt 0.5) (sqrt 0.5))
(size 1 1 1)
)
)
(set! geometry
(list
(make sphere
(center 0 0 0)
(radius (/ (sqrt 0.5) 2))
(material
(make dielectric
(index 1.0)
)
)
)
)
)
(set! k-points
(interpolate 4
(list
(vector3 0 0.5 0.5) ; X
(vector3 0 0.625 0.375) ; U
(vector3 0 0.5 0) ; L
(vector3 0 0 0) ; Gamma
(vector3 0 0.5 0.5) ; X
(vector3 0.25 0.75 0.5) ; W
(vector3 0.375 0.75 0.375) ; K
)
)
)
(run)
maryam saberi <[EMAIL PROTECTED]> wrote:
> Hello
> I' d like simulate a fcc inverse opal structure with Mit Photonic Bands
> package, but I don't know how to write it's program.
> Please help me.
> thanks
>
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
nate lipkowitz
416 823 8057
[EMAIL PROTECTED]
_______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
