Dear Meep users:
I am trying to use symmetry to reduce the computation effort.
My structure is the bowtie antenna and I used plane wave to impede it.
My original idea is cut computation domain into quarter. However, when I
calculated the electrical energy in the metal,
the value is going up to infinity. If I didn't apply symmetry, this wouldn't
happen. I'm wondering if my symmetry application is wrong.
Attach ctl file is below:
Thank you
Cheers
Adrian
(define lam 1550) ; free space wavelength (meters)
(define fcen (/ 1 lam))
; Drude model of Au
(define-param wp_au 6.7368e-3) ; resonant plasma frequency
(define-param gama_au 6.556e-5) ; gamma of Gold
(define-param e_au 12) ; epsilon of Gold 9.84
; Geometry
(define-param L 150 ) ; length of triangle
(define-param t_m 24 ) ; thickness of metal
(define-param t_ito 52 ) ; thickness of ito
(define-param t_inp 100 ) ; thickness of silica
(define-param g 20) ; gap between bowtie antenna
(define-param ang (/ (* 60 pi) 180)) ; angle of triangle
(define-param r 12) ; curvature of apex
(define-param time_turnon (* 0.5 lam)) ;
(define-param dpmlz 64 ) ; thickness of PML in z
; the distance between object and PML
(define-param padx 120) ; to fit the length equal to 800nm in the x direction
(define-param pady 120) ; to fit the length equal to 800nm in the y direction
(define-param padz 128) ; to fit the length equal to 176nm in the z direction
; lattice size
(define-param sy (+ (* L (sin (/ ang 2)) 2) pady pady) ) ; (+ (* L (sin (/ ang
2)) 2) padx padx)
(define-param sx (+ (* L (cos (/ ang 2)) 2) g padx padx) ) ;(+ (* L (cos (/ ang
2)) 2) g pady pady)
(define-param sz (+ t_m t_ito t_inp padz )) ;
(set! geometry-lattice (make lattice (size sx sy sz)))
(set-param! resolution 0.5)
(set! eps-averaging? false)
(define Au
(make dielectric (epsilon e_au)
(polarizations
(make polarizability
(omega 1e-20) (gamma gama_au) (delta-epsilon (* wp_au wp_au
1e+20 1e+20))))))
(define ito (make dielectric (epsilon 3.8025)))
(define inp (make dielectric (epsilon 10.03)))
(set! pml-layers (list (make pml (thickness dpmlz) (direction Z))))
(set! k-point (vector3 0 0 0))
(set! geometry
build bowtie structure
)
(define-param kdir (vector3 0 0 -1)) ; direction of k (length is irrelevant)
(define k (vector3-scale (* 2 pi fcen 3.17)
(unit-vector3 kdir))) ; k with correct length
(define ((pw-amp k x0) x)
(exp (* 0+1i (vector3-dot k (vector3+ x x0)))))
(define-param run_time 9 )
(set! sources
(list
(make source
(src (make continuous-src (frequency fcen) (end-time (* run_time (/ 1
fcen))) (width (* 0.8 lam)) (cutoff 3.424 )))
(component Ex) (center 0 0 (+ t_ito (* (- t_inp dpmlz) 0.5))) (size sx sy
0 )
(amp-func (pw-amp k (vector3 0 0 (+ t_ito (* (- t_inp dpmlz) 0.5))))))
))
(define (f1 r ez ey ex eps)
(if (> eps 10.03)
(+ (* ez (conj ez)) (* ey (conj ey)) (* ex (conj ex)))
0
))
(define Een_m 0)
(define (do-Een_m)
(set! Een_m (integrate-field-function (list Ez Ey Ex Dielectric) f1 )) )
(set! symmetries (list (make mirror-sym (direction Y) (phase 1)) (make
mirror-sym (direction X) (phase -1))))
(define-param execute_time 20 )
(run-until (* execute_time (/ 1 fcen))
(after-time (* (- execute_time 10) (/ 1 fcen)) do-Een_m )
)
______________________________________________________________________________________
Yahoo!奇摩迷你筆，英文單字一點就翻。http://tw.mini.yahoo.com/
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
