Dear Steven,
Thank you for replying so fast.
The half of the modes I means that by using the rectangular lattice, the number
of the waveguide modes is smaller than that by using triangular lattices. Could
you please help me checking the following codes? It uses the rectangular
lattice. It should be very simple for you. For W1, the results are different
for the triangular and rectangular lattices, too.
Many thanks!
Best,
Jiangjun
; Here, Cartesian Coordination is used
;------------------------------------------------
(define a 1)
(define r (* 0.29 a))
(define h (* 0.6 a))
(define s1 (* 1 a))
(define s2 (* 1 a))
(define s3 (* 4 a))
(define d1 0.5)
(define d2 (* 0.5 (sqrt 3)))
; Define the positions of holes
(define wg-ratio 1) ; 1 for W1 defects
(define shift-y (* (- 1 wg-ratio) d2))
(define supercell-2 (* d2 (+ 5 (* 2 wg-ratio))))
(define p1 (vector3 0 (- (* d2 1) shift-y) 0))
(define p2 (vector3 0 (- (* d2 3) shift-y) 0))
(define p3 (vector3 0.5 (- (* d2 2) shift-y) 0))
(define p4 (vector3 -0.5 (- (* d2 2) shift-y) 0))
(define m1 (vector3 0 (- shift-y (* d2 1)) 0))
(define m2 (vector3 0 (- shift-y (* d2 3)) 0))
(define m3 (vector3 0.5 (- shift-y (* d2 2)) 0))
(define m4 (vector3 -0.5 (- shift-y (* d2 2)) 0))
;------------------------------------------------
;------------------------------------------------
(define eps-si 12.11)
(define eps-eff-si (sqr 2.77))
(define mat-si (make dielectric (epsilon eps-si)))
(define eff-si (make dielectric (epsilon eps-eff-si)))
;------------------------------------------------
;------------------------------------------------
(define cylinder-si (make cylinder
(center 0) (material mat-si)
(radius r) (height h)))
(define cylinder-air (make cylinder
(center 0) (material air)
(radius r) (height h)))
(define slab-si (make block
(center 0) (material mat-si)
(size infinity infinity h)))
;------------------------------------------------
;------------------------------------------------
(define-param k-interp 19)
(set-param! num-bands 15)
;------------------------------------------------
; Main body of ctl file
(set! default-material air) ; Air-bridged
(set! geometry-lattice (make lattice
(size 1 supercell-2 s3)))
(set! geometry (list
slab-si ; make slab
(shift-geometric-object cylinder-air p1)
(shift-geometric-object cylinder-air p2)
(shift-geometric-object cylinder-air p3)
(shift-geometric-object cylinder-air p4)
(shift-geometric-object cylinder-air m1)
(shift-geometric-object cylinder-air m2)
(shift-geometric-object cylinder-air m3)
(shift-geometric-object cylinder-air m4)))
(set-param! resolution (vector3 32 32 16))
(set-param! mesh-size 7)
(define kz 0) ; In-plane
(define Gamma (vector3 0 0 0))
(define K' (vector3 0.5 0 kz))
(set! k-points (interpolate k-interp (list Gamma K')))
(run-zeven
(output-at-kpoint K' fix-hfield-phase output-hfield-z)
(output-at-kpoint K' fix-efield-phase output-efield-y))
(display-eigensolver-stats)
(set! mesh-size 3)
好消息，您的网易邮箱容量永远用不完了！_______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
