Hello everybody,
I am using Meep to do some simulation about the propagation of
light in metallic slit. The struture is every simple, just a silver
plate with a slit of air. A point source is put in front of the
entrance of the slit. The dielectric function of silver is modeled by
Drude model. The PML layers are put on the boundaries in all
directions. After some steps, a strange field occur at the joint of
the metal and the PML layer. A figure of the result is attached in
this message. I am confused and don't know what's wrong in my ctl
script.
The version of meep I use is 0.20.3 and all tests passed when "make check".
Thanks!
My ctl script:
(define a 100e-9) ; define the unit, 100 nm
(define sx 40) ; size of x, 4 um
(define sy 10) ; size of y
(define dpml 1)
(set! geometry-lattice
(make lattice (size sx sy no-size)))
; ------------------------------------
; define material silver
(define cc (/ (* 2 (* pi 3e+8)) a)) ; unit of frequency in meep
(define wp_ag (/ 1.3798e16 cc))
(define gamma_ag (/ 2.72293e13 cc))
(define Ag
(make dielectric (epsilon 3.7)
(polarizations (make polarizability
(omega 1e-20)
(gamma gamma_ag)
(delta-epsilon (* wp_ag wp_ag 1e40))))))
(define wid 1)
(set! geometry
(list
(make block
(center 0 0)
(size 20 sy infinity)
(material Ag))
(make block
(center 0 0)
(size 20 wid infinity)
(material air))))
; set source
(set! sources (list
(make source
(src (make continuous-src (wavelength 6.5) (width 20)))
(component Hz)
(center -15 0)
(size 0 0))))
; ---------------------------------------------
(set! pml-layers (list
(make pml (thickness dpml) (direction ALL))))
(set! resolution 10)
(run-until 50 (at-beginning output-epsilon)
(at-every 2 (output-png Hz "-Zc dkbluered")))JL Liu 2008-09-25
test-hz-000012.00.png
Description: Binary data
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
