dear meep user Now, I am working on simulation of photonic crystal. I am doing study about anomalous refraction in photonic crystal with triangular lattice. I am using plane wave as source. I have question about how to set the angle of incidence of light. I hope you can help me. I attach my ctl file. I modify it from code by Karel Van Acoleyen from meep discuss digest volume 22 issue 17. Please give me your suggestion about anything if you like. Thank you
Best regard
(define silicon (make dielectric (index 12.96)))
(define silica (make dielectric (index 1.45)))
(define SiN (make dielectric (index 2)))
(define-param h 0.5)
(define-param r 0.35)
(define-param h_air 2)
(define-param h_Si 0.255)
(define-param h_SiO2 0.61576)
(define-param n_bragg 9)
(define-param h_SiN 0.44643)
(set! geometry-lattice (make lattice (size 1 (* 13 (sqrt 3)) 3))) ;(+ (* h_air
2) h (* 2 n_bragg (+ h_SiN h_SiO2)))))
(set! default-material air) ; everything is air
(set! geometry
(append
(list (make block (center 0 0 0) (size infinity infinity h)
(material air)))
(geometric-object-duplicates (vector3 0 2 0) 0 1 (make cylinder
(center (/ (sqrt 3) 2) 0 0) (radius r) (height h) (material silicon)))
(geometric-object-duplicates (vector3 0 2 0) 0 1 (make cylinder
(center (/ (sqrt 3) -2) 0 0) (radius r) (height h) (material silicon)))
(geometric-object-duplicates (vector3 0 2 0) 0 1 (make cylinder
(center (* (/ (sqrt 3) 2) 3) 0 0) (radius r) (height h) (material silicon)))
(geometric-object-duplicates (vector3 0 2 0) 0 1 (make cylinder
(center (* (/ (sqrt 3) 2) -3) 0 0) (radius r) (height h) (material silicon)))
(geometric-object-duplicates (vector3 0 -2 0) 0 1 (make cylinder
(center (/ (sqrt 3) 2) 0 0) (radius r) (height h) (material silicon)))
(geometric-object-duplicates (vector3 0 -2 0) 0 1 (make cylinder
(center (/ (sqrt 3) -2) 0 0) (radius r) (height h) (material silicon)))
(geometric-object-duplicates (vector3 0 -2 0) 0 1 (make cylinder
(center (* (/ (sqrt 3) 2) 3) 0 0) (radius r) (height h) (material silicon)))
(geometric-object-duplicates (vector3 0 -2 0) 0 1 (make cylinder
(center (* (/ (sqrt 3) 2) -3) 0 0) (radius r) (height h) (material silicon)))
(geometric-object-duplicates (vector3 0 2 0) 0 1 (make cylinder
(center 0 -1 0) (radius r) (height h) (material silicon)))
(geometric-object-duplicates (vector3 0 2 0) 0 1 (make cylinder
(center (sqrt 3) -1 0) (radius r) (height h) (material silicon)))
(geometric-object-duplicates (vector3 0 2 0) 0 1 (make cylinder
(center (* (sqrt 3) -1) -1 0) (radius r) (height h) (material silicon)))))
(set! symmetries (list (make mirror-sym (direction Z) (phase -1)) (make
mirror-sym (direction Y) (phase -1))))
(set! pml-layers (list (make pml (direction Y) (thickness 1)) (make pml
(direction Z) (thickness 1))))
(define-param fcen 0.535) ; pulse center frequency
(define randomnumber 0.13)
(define-param k-interp 10)
(set-param! resolution 16)
(set! sources (list
(make source
(src (make continuous-src (frequency fcen)))
(component Hz ) (center 0 0 0)
)
)
)
(run-until 200
(at-beginning output-epsilon)
(at-end output-hfield-z))
____________________________________________________________________________________
Be a better friend, newshound, and
know-it-all with Yahoo! Mobile. Try it now.
http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ
TA1.ctl
Description: Binary data
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
