Hello!
I tried to use meep to simulate the propagation of surface plasmon in
a 2d case waveguide. The unit of meep is set as 100nm. The metal is
assumed to be lossless, with a bulk plasmon frequency to be 1e16
rad/s. The whole control file is as follows:
*************************************************************************************************
(set! geometry-lattice (make lattice (size 8 3 no-size)))
(set! default-material
(make dielectric (epsilon 1)
(polarizations
(make polarizability
(omega 1e-20) (gamma 0) (delta-epsilon 0.53088e+040))
)))
(set! geometry (list
(make block (center 0 0) (size 6 1 infinity)
(material (make dielectric (epsilon 1.5))))))
(set! sources (list
(make source
(src (make continuous-src (frequency 0.15803))) ;free
space wavelength632.8nm
(component Ey)(amplitude 1)
(center -3 0))))
(set! pml-layers (list (make pml (thickness 1))))
(set! resolution 50)
(run-until 300
(at-beginning output-epsilon)
(at-end output-efield-y))
******************************************************************************************************
A snapshot of the electric field is shown as attached. It seems to be
rather quirky, though the field is concentrated as the surfaces. The
plasmon wavelength can not be so short! And,at the right end, where
the waveguide meets the pml tightly, still
there are fields existing.
I am puzzled with these. Any explanations?
Very much thanks!
2dspp-ey-000300.00.png
Description: PNG image
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
