Hello everybody
I have been trying to simulate plasmons using meep and a code like the
following for some time.
I got bound solutions and I found out that the decay length into the metal is
quite close to the analytically expected decay length.
But I have not been able to get a propagation length which is of the magnitude
of the theoretically expected one.
Does anybody know what's wrong about the code?
Thanks in advance
Johannes
Code for the materials:
_________________________________________________________________________________________________
(define Au
(make dielectric (epsilon 1)
(E-polarizations
(make polarizability (omega (* a 1e-13)) (gamma (* a 42717.8)) (sigma
4.0258e+39))
(make polarizability (omega (* a 334488.6)) (gamma (* a 194245.2)) (sigma
11.3629))
(make polarizability (omega (* a 668977.3)) (gamma (* a 278068.9)) (sigma
1.183639))
(make polarizability (omega (* a 2393004)) (gamma (* a 701217.2)) (sigma
0.65677))
(make polarizability (omega (* a 3469010)) (gamma (* a 2010155.9)) (sigma
2.645486))
(make polarizability (omega (* a 10735876.6)) (gamma (* a 1784476.8))
(sigma 2.01482623))
)
)
)
________________________________________________________________________________________________
Simulation:
(define-param a 1e-6) ; unit size (equals 1 um)
(define-param l 20)
(define-param w 6)
(define-param h 3)
(define-param dpml 1)
(define-param pad 4)
(define-param sx l)
(define-param sy w)
(define-param sz (* 2 h))
(define-param zcen (* -0.5 h))
(include "materials_Rakic.h")
(set-param! resolution 10)
(set-param! Courant 0.1)
(set! geometry-lattice (make lattice (size sx sy sz)))
(set! eps-averaging? false)
(set! geometry (list
(make block (center 0 0 (* -0.5 h)) (size l w h) (material Au))
))
(set! pml-layers (list (make pml (thickness dpml))))
(define-param fcen (* a 1200000))
(set! sources (list
(make source
(src (make continuous-src
(frequency fcen)
))
(component Ez)
(center (+ 2 (* -0.5 l)) 0 0)
(size 0 0 0)
)
))
(use-output-directory)
(run-until 200
(at-beginning output-epsilon)
(to-appended "hy" (at-every 1 output-hfield-y))
)
___________________________________________________________________________________
In order to find out about the decay length into the metal I used a much
smaller cell but a higher resolution.
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
