Hi Meep user I am trying to simulate narrow band filter example in the PhC molding the flow of light book, page 196-197, figure 4.
First question: I have obtained the resonance frequency=0.3786 and Q=510. Although the dimension, pml and resolution are increased why I cannot get the exact value? Second question: According to the figure, 0.3765 frequency is not transmitted but resonance frequency 0.3803 is transmitted. I have observed the 0.3765 frequency is blocked in my simulation. But, neither the resonance frequency 0.3786 (according to my calculation) nor 0.3803 (according to book) is transmitted. What I missed the point? my ctl codes for Q and field transmittance are as the following best… Onder ------------------------------------ ;; Q calculation ;; PhC Molding the flow of light, page 196-197 , Figure 4 (define-param nbase 1) ;air (define-param epsbase (* nbase nbase)) (define-param epsrod 11.4) ;; rods (define-param a 1) (define-param r 0.2) ; rods radius (define-param sx 25) (define-param sy 15) (set! geometry-lattice (make lattice (size sx sy no-size))) (set! default-material (make dielectric (epsilon epsbase))) (set! geometry (list (make cylinder (center 0 0) (height infinity) (radius r) (material (make dielectric (epsilon epsrod)))) )) (set! geometry (geometric-objects-lattice-duplicates geometry)) (set! geometry (append geometry (geometric-object-duplicates (vector3 1 0) 0 10 (make cylinder (center 3 0) (radius r) (height infinity) (material (make dielectric (epsilon epsbase)))) ))) (set! geometry (append geometry (geometric-object-duplicates (vector3 -1 0) 0 10 (make cylinder (center -3 0) (radius r) (height infinity) (material (make dielectric (epsilon epsbase)))) ))) (set! geometry (append geometry (list (make cylinder (center 0 0) (radius r) (height infinity) (material (make dielectric (epsilon epsbase) )) ) ))) (define-param fcen 0.38) (define-param df 0.1) (set! sources (list (make source (src (make gaussian-src (frequency fcen) (fwidth df) )) (component Ez)(center 0 0) ) )) (set! pml-layers (list (make pml (thickness 2.0)))) (set! resolution 32) (run-sources+ 300 (at-beginning output-epsilon) (after-sources (harminv Ez (vector3 0) fcen df))) ---------------------- ---------------------- ;; field transmission ;; PhC Molding the flow of light, page 196-197 , Figure 4 (define-param nbase 1) ;air (define-param epsbase (* nbase nbase)) (define-param epsrod 11.4) ;; rods (define-param a 1) (define-param r 0.2) ; rods radius (define-param sx 25) (define-param sy 15) (set! geometry-lattice (make lattice (size sx sy no-size))) (set! default-material (make dielectric (epsilon epsbase))) (set! geometry (list (make cylinder (center 0 0) (height infinity) (radius r) (material (make dielectric (epsilon epsrod)))) )) (set! geometry (geometric-objects-lattice-duplicates geometry)) (set! geometry (append geometry (geometric-object-duplicates (vector3 1 0) 0 10 (make cylinder (center 3 0) (radius r) (height infinity) (material (make dielectric (epsilon epsbase)))) ))) (set! geometry (append geometry (geometric-object-duplicates (vector3 -1 0) 0 10 (make cylinder (center -3 0) (radius r) (height infinity) (material (make dielectric (epsilon epsbase)))) ))) ;;; hole at the center (set! geometry (append geometry (list (make cylinder (center 0 0) (radius r) (height infinity) (material (make dielectric (epsilon epsbase) )) ) ))) (set! sources (list (make source (src (make continuous-src (frequency 0.3803) )) (component Ez)(center -9 0) (size 0 1) (amplitude 1) ) )) (set! pml-layers (list (make pml (thickness 2.0)))) (set! resolution 32) (run-until 400 (at-beginning output-epsilon) (to-appended "ez" (at-every 5 output-efield-z)))
_______________________________________________ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss