Dear Meep users,
I want to simulate a spatially varying dispersive material of which
the epsilon and sigma should be functions of position. But after searching
all related topics to apply to my .ctl file, I found it doesn't work. For a
simple example, I want to simulate a dispersive material as a thin film, and
there is a slit as vacuum on it. The material I made looks more like a
dielectric by checking the output results, meaning that sigma does not take
effect. The output epsilon is fine, but unfortunately, I can't check the
sigma. Can anyone tell me how to set it actually?
Jason
My code:
; The parameters================================================
(define-param Al_flm_Ly 0.10) ; the thickness of the metal film, t
(define-param Al_flm_Lx 0.10) ; the thickness of the metal film, w
(define-param Al_flm_Lz 0.10) ; the thickness of the metal film, L
(define-param Al_flm_cntr 0.51) ; the center of the metal film
(define-param slt_Lx 0.21) ; the width of the slit, w
(define-param slt_Lz 0.21) ; the length of the slit, L
; The spatially varying function of epsilon and sigma=====================
(define (eps p)
(if (and (< (vector3-x p) (/ slt_Lx 2)) (> (vector3-x p) (/ (- 0
slt_Lx) 2)) )
(if (and (< (vector3-z p) (/ slt_Lz 2)) (> (vector3-z p) (/
(- 0 slt_Lz) 2)) ) 1 3)
3)
)
(define (sgm w)
(if (and (< (vector3-x w) (/ slt_Lx 2)) (> (vector3-x w) (/ (- 0
slt_Lx) 2)) )
(if (and (< (vector3-z w) (/ slt_Lz 2)) (> (vector3-z w) (/
(- 0 slt_Lz) 2)) ) 0 1.6775468e+042)
1.6775468e+042)
)
; The dispersive material ==============================================
(define (Al p1) (make medium (epsilon (eps p1))
(E-polarizations
(make polarizability
(omega 1.000000e-20)
(gamma 0.462362658)
(sigma (sgm p1))))))
; Making the material===================================================
(set! geometry (list
(make block (center 0.0 Al_flm_cntr 0.0) (size Lx_phy Al_flm_Ly Lz_phy)
(material (make material-function (material-func Al))))
)
)
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
