Dear Steven, dear MEEP user,
I am kindly asking you for help.
I have a problem with calculation of the band structure of metallic rods
arranged in a triangular lattice. I started with 2D problem with rods of
dielectric material to compare it with the result of MPB calculation. I
am wondering that results are not consistent. I guess some mistake in my
definition of the lattice in MEEP ctl file.
The listing is here
(define-param eps 14);
(define-param sx 1)
(define-param sy (/ 1 (sqrt 3) ))
(define-param a sy)
(define-param r (* a 0.1) ) ; radius of holes
(define-param K (/ 1 (* 2 (sqrt 3))))
(set! geometry-lattice (make lattice (size sx sy no-size)))
(set! geometry
(list
(make cylinder (center 0 0) (radius r) (height
infinity) (material(make dielectric (epsilon eps)))
)
(make cylinder (center (/ sx 2) (/ sy 2)) (radius r) (height
infinity) (material(make dielectric (epsilon eps)))
)
(make cylinder (center (/ sx -2) (/ sy 2)) (radius r)
(height infinity) (material(make dielectric (epsilon eps)))
)
(make cylinder (center (/ sx 2) (/ sy -2)) (radius r)
(height infinity) (material(make dielectric (epsilon eps)))
)
(make cylinder (center (/ sx -2) (/ sy -2)) (radius r)
(height infinity) (material(make dielectric (epsilon eps)))
)
)
)
(set-param! resolution 10)
(define-param fcen 0.31) ; pulse center
frequency
(define-param df 1.5) ; pulse freq. width: large df = short impulse
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ez) (center 0.1234 0))))
(define-param k-interp 20) ; # k-points to interpolate, otherwise
(set! k-point (vector3 0.5 0))
(set! ensure-periodicity true)
(at-beginning output-epsilon)
(run-k-points 300 (interpolate k-interp (list (vector3 0 0) (vector3 .5
0.0) (vector3 0.5 K) (vector3 0.0 0. )))
Best Regards,
Maksim
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
