photonic crystal slab

Dear Steven,

      I  have studied MPB for several days. I repeated Fig.3 (a) of  Fan's
article (PRB 65, 235112). The following is my result and Fan's figure. I
don't understand there are  some lines above the light line in the figure.
There are no these lines in Fan's figure. Can you help me about this? Thank
you very much! I can't paste the figures because the mail can't be larger
than 40k.



This is the control file:

(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.2) ; the radius of the holes

(define-param supercell-h 4) ; height of the supercell



; square lattice with vertical supercell:

(set! geometry-lattice (make lattice (size 1 1 supercell-h)))



(set! geometry

      (list (make block (material (make dielectric (epsilon loweps)))

                          (center 0 0 (* 0.25 supercell-h))

                          (size infinity infinity (* 0.5 supercell-h)))

                (make block (material (make dielectric (epsilon eps)))

                          (center 0) (size infinity infinity h))

                (make cylinder (material air)

                          (center 0) (radius r) (height supercell-h))))



; 1st Brillouin zone of a square lattice:

(define Gamma (vector3 0 0 0))

(define X (vector3 0.5 0 0))

(define M (vector3 0.5 0.5 0))



(define-param only-M false) ; run with only-K=true to only do this k-point

(define-param k-interp 12)   ; the number of k points to interpolate



;(set! k-points (interpolate k-interp (list Gamma X M Gamma)))

(if only-M

    (set! k-points (list M))

    (set! k-points (interpolate k-interp (list Gamma X M Gamma))))



(set-param! resolution (vector3 32 32 16))

(set-param! num-bands 9)



; Run even and odd bands, outputting fields only at the M point:

(run-zeven (output-at-kpoint M output-hfield-z))

(run-zodd (output-at-kpoint M output-dfield-z))



(display-eigensolver-stats)





Yours Sincerely,

Lina

_______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss