Greetings MEEP users
I tried to simulate case for the Poisson's spot.
I made a gold slab and generated the plane wave source
But the problem is light is penetrating the gold slab.
I increased the thickness of gold to avoid skin depth problem,
but the result was same. could you help me to understand this?
Here is my code,
-----------------------------------------------------------------------------------------------
(define-param source-component Ez)
(define-param theta_deg 0)
(define-param fcen 0.8) ; pulse center frequency
(define-param df 0.02) ; turn-on bandwidth
(define wlength (/ 1 fcen)) ; pulse wavelength
(define-param sx 16) ; the size of the comp cell in X, not including PML
(define-param dpml 5) ; thickness of PML layers
; set periodicity in the X direction. The PML takes care of the
; non-desired Y periodicity that comes for free with the k-point declaration
(set! k-point (vector3 0 0 0))
(set! ensure-periodicity true)
(define gold (make dielectric
(epsilon 4.98)
(E-polarizations
(make polarizability (omega 1e-10) (gamma 6.13e-4) (sigma (* .142 .142
1e20)))
(make polarizability (omega 0.0592) (gamma 0.0211) (sigma 1.76))
(make polarizability (omega 0.0458) (gamma 0.0119) (sigma 0.952))
)
))
(define sy (+ sx (* 2 dpml))) ; cell size in Y direction, including PML
(set! geometry-lattice (make lattice (size sx sy no-size)))
(set! geometry(list(make block (center 0 0) (size 4 2
infinity)(material gold))))
; we'll only have PML in the Y direction (top and bottom)
(set! pml-layers (list (make pml (thickness dpml) (direction Y))))
(set-param! resolution 50)
; pw-amp is a function that returns the amplitude exp(ik(x+x0)) at a
; given point x. (We need the x0 because current amplitude functions
; in Meep are defined relative to the center of the current source,
; whereas we want a fixed origin.) Actually, it is a function of k
; and x0 that returns a function of x ...
(define ((pw-amp k x0) x)
(exp (* 0+1i (vector3-dot k (vector3+ x x0)))))
(define theta_rad (/ (* pi theta_deg) 180))
; direction of k (length is irrelevant)
(define-param kdir (vector3 (sin theta_rad) (cos theta_rad)))
; k with correct length
(define k (vector3-scale (* 2 pi fcen) (unit-vector3 kdir)))
(set! sources (list
(make source
(src (make continuous-src (frequency fcen) (fwidth df)))
(component source-component) (center 0 (* -0.5 sx)) (size sx 0)
(amp-func (pw-amp k (vector3 0 (* -0.5 sx)))))
))
(run-until 40
(at-every 1 (output-png source-component "-Zc bluered")))
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
