Hi dear users.
I want to calculate the transmission (Fig.3c) that is illustrated in this paper
(Phys. Rev. B, Vol. 54, Num. 11 , pages:7837-7842), but my result is different.
The program that I've used in Meep is as below, but the result is different
from Fig. 3c. Help me to correct the mistake please.
Best.
Hodjat
;some parameters to describe the geometry:
(define-param eps 1.843)
(define-param r 0.2) ; radius of holes
(define-param Nx 7) ; number of holes in X direction
(define-param Ny 11) ; number of holes in Y direction
(define-param rd 0.2) ; r of defect
; the cell dimensions:
(define-param pad 4) ; padding from last hole and pml edge
(define-param dpml 1) ; pml thickness
(define sx (+ (* 2 (+ pad dpml r)) Nx -1))
(define sy (+ (* 2 (+ dpml r)) Ny
-1))
(define wx (+ Nx -1 (* 2 r)))
(define wy (+ Ny -1 (* 2 r)))
;*******************************
(set! geometry-lattice ( make lattice (size sx sy no-size)))
(set! geometry
(list
(make cylinder (center 0 0 0) (radius r) (height infinity) (material (make
dielectric (epsilon eps))))
)
)
(set! geometry (geometric-objects-lattice-duplicates geometry))
(set! geometry
(append geometry
(list
(make cylinder (center 0 0 0) (radius rd) (height infinity) (material air))
)
(list
(make block (center (/ (+ sx wx) 4) 0) (size (* 0.5 (- sx wx)) infinity
infinity) (material air))
)
(list
(make block (center (/ (+ (* -1 wx) (* -1 sx)) 4) 0) (size (* 0.5 (- sx wx))
infinity infinity) (material air))
)
(list
(make block (center 0 (/ (+ sy wy) 4)) (size infinity (* 0.5 (- sy wy))
infinity) (material air))
)
(list
(make block (center 0 (/ (+ (* -1 wy)
(* -1 sy)) 4)) (size infinity (* 0.5 (- sy wy)) infinity) (material air))
)
)
)
;*********************************
(set! pml-layers (list (make pml (thickness dpml))))
(set-param! resolution 20)
;some parameters to describe the source:
(define-param fcen 0.35)
(define-param df 0.20)
(define-param nfreq 500)
(set! sources (list
(make
source (src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ez) (center (+ dpml (* -0.5 sx)) 0) (size 0 sy)) ;
))
(set! symmetries (list
(make mirror-sym (direction Y) (phase -1))
))
(define trans (add-flux fcen df nfreq
(make flux-region (center (- (* 0.5 sx) dpml 0.5) 0) (size 0 sy))
))
(run-sources+
(stop-when-fields-decayed 50 Ez (vector3 (- (* 0.5 sx) dpml 0.5) 0) 1e-3)
(at-beginning output-epsilon)
)
(display-fluxes trans)
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
