; xy triangular lattice of cylindrical rods ; ; Rod radius is R*lattice size = R ; Rod permitivity is (n^2+k^2)*Eo ; Rod height is infinite ; Rod axis // z ; turned around x towards y by A ; turned around z towards y by B ; ; Plane wave basis MxMx1 ; Number of computed intermediate k-vector steps = S ; ; +------------------+ ; | Input parameters | ; +------------------+ ; (define-param A 0) (define-param B 0) (define-param k 0.00) (define-param n 1.87) (define-param R 0.30) (define-param M 32) (define-param S 8) ; ; +------------------+ ; | Relative epsilon | ; +------------------+ ; (define E (+ (* k k) (* n n))) (display "\n\n \n \n MIT Photonic-Bands output\n\n") (display "

Hexagonal crystal geometry

\n\n\n") (display " \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n \n
Background dielectric constantεb=1.0000ε0
Dielectric constant of the rodsεr=") (display E) (display "ε0

Rod axis [default: z-axis]
  turned around the x-axis towards the y-axisα=") (display A) (display "°
  turned around the z-axis towards the y-axisβ=") (display B) (display "°
Radius of the rodsr=") (display R) (display "lattice

The computation bases on a set of ") (display (* M M)) (display " plane waves.
A total of ") (display (+ (* S 3) 4)) (display " k-vectors will be examined.
\n") ; ; +------------------+ ; | Lattice geometry | ; +------------------+ ; (set! geometry-lattice (make lattice (basis1 (sqrt 0.75) -0.5 0) (basis2 (sqrt 0.75) 0.5 0) (basis3 0 0 1) ) ) ; ; +-----------------+ ; | Rod axis vector | ; +-----------------+ ; (define axe (rotate-vector3 (vector3 0 0 1) (deg->rad B) (rotate-vector3 (vector3 (sqrt 0.75) -0.5 0) (deg->rad (- A)) (vector3 0 0 1)))) ; ; +----------------+ ; | Setup geometry | ; +----------------+ ; (set! geometry (list (make cylinder (axis axe) (center 0 0 0) (height infinity) (radius R) (material (make dielectric (epsilon E))) ) )) ; ; +--------------+ ; | Quantization | ; +--------------+ ; (set! grid-size (vector3 M M 1)) (set! k-points (list (vector3 (/ -3) (/ 3) 0) ; K (vector3 0 0 0) ; G (vector3 0 0.5 0) ; M (vector3 (/ -3) (/ 3) 0) ; K )) (set! k-points (interpolate S k-points)) (set! num-bands 20) (display "

TE-Polarization

\n
")
(run-te display-group-velocities
;	(output-at-kpoint (vector3 0 0.5 0) fix-efield-phase output-efield)
)
(display "

TM-Polarization

\n
")
(run-tm display-group-velocities
;	(output-at-kpoint (vector3 0 0.5 0) fix-efield-phase output-efield)
)
(display "

Comments

\n\n

Created with the HTML generator written by Marcel Leutenegger.

\n\n")