Hi Cefele,
The structure you want is not what you define in your script, you can see
something is wrong from the surprisingly small number of bands outside the
light cone (usually characterizing low index contrast structures in z, not your
case here). I have extended the loweps material to infinity and have it
centered at z=0 and also extended the air-holes to infinity in z (although this
last step should not affect the result).
I obtain a band diagram which looks fine for all the bands you mention with at
least all points between M and K outside the light cone.
I hope this is what you want,
Regards,
Frederic
I have modified the script so that the air just extends to infinity
(define-param h 0.5) ; the thickness of the slab
(define-param eps 12.0) ; the dielectric constant of the slab
(define-param loweps 1.0) ; the dielectric constant of the substrate
(define-param r 0.3) ; the radius of the holes
(define-param supercell-h 4) ; height of the supercell
; triangular lattice with vertical supercell:
(set! geometry-lattice (make lattice (size 1 1 supercell-h)
(basis1 (/ (sqrt 3) 2) 0.5)
(basis2 (/ (sqrt 3) 2) -0.5)))
(set! geometry
(list (make block (material (make dielectric (epsilon loweps)))
(center 0)
(size infinity infinity infinity))
(make block (material (make dielectric (epsilon eps)))
(center 0) (size infinity infinity h))
(make cylinder (material air)
(center 0) (radius r) (height infinity))))
; 1st Brillouin zone of a triangular lattice:
(define Gamma (vector3 0 0 0))
(define M (vector3 0 0.5 0))
(define K (vector3 (/ -3) (/ 3) 0))
(define-param only-K false) ; run with only-K=true to only do this k-point
(define-param k-interp 4) ; the number of k points to interpolate
(if only-K
(set! k-points (list K))
(set! k-points (interpolate k-interp (list Gamma M K Gamma))))
(set-param! resolution (vector3 32 32 32))
(set-param! num-bands 9)
; Run even and odd bands, outputting fields only at the K point:
(run-zeven (output-at-kpoint K output-hfield-z))
On Apr 20 2009, [email protected] wrote:
(Message body was not text: suppressed)
_______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
