;璁扮殑鏀筫psinf
(reset-meep)
(set! resolution 100)
;(set! k-point (vector3 0 0 0))
;(set! pml-layers (list (make pml (direction Z) (thickness 0.2))))
;(set! pml-layers (list (make pml (direction X) (thickness 0.1))))
(set! pml-layers (list (make pml (direction Y) (thickness 0.1))))
(define-param pi 3.14159265358979323846)
(define-param c 2.99792458e8)
(define-param aa 200e-9); aa is the lattice constant
(define-param pi2cqaa (/ (* 2 pi c) aa));
(define-param a 1); a is the lattice constant
(define-param epsrod 1)
;(define-param bc (/ 280e-9 aa)); bc is the side length of square holes
;(define-param bcy (/ (* 280e-9 (sqrt 2) (/ 1 3)) aa))
;(define-param bcx (* bcy 9 0.5))
(define-param r1 (/ 50e-9 aa)); r1 is the outer raduii of the nanotube;
(define-param r2 (/ 20e-9 aa)); r2 is the inner raduii of the nanotube;
;; prl,95,137404, 2005 Au , gold.
(define-param epsinf 1)
(define-param wp 1.37e16)
(define-param tau 4.08e13)
(define-param wpw (/ wp pi2cqaa))
(define-param tauw (/ tau pi2cqaa))
(define-param x (* 1 a))
(define-param y (* 7.5 a))
;(define-param z (* 1.5 a))