You should only be able to activate the SPP for TM polarization.. Which appears to indeed be what you are observing! --Judson
On Mon, Dec 10, 2012 at 9:35 AM, huang bo <abohu...@gmail.com> wrote: > hi everyone! > below is the .ctl file I use to calculate the dispersion relation of surface > plasmon > I can get the dispersion relation when the source component set as Hz and > ky run from kmin to kmax. > the queer problem is that > when I set the source component as Hy, and kz run from kmin to kmax, I can > not get the dispersion relation. > Anyway,both configurations are the same for two infinite metal and air > layer. > > anyone have any idea about that? > thank! > > > (define-param sx 4) > (set! geometry-lattice (make lattice (size sx no-size no-size))) > (set-param! resolution 40) > ;(set! eps-averaging? false) > (define omega_eps 0.3) > (define Ag51 (make dielectric (epsilon 6.8) > (polarizations > (make polarizability > (omega 1e-20) (gamma 0.003226) (sigma (* 6.8 omega_eps omega_eps 1e+40))) > ))) > (set! geometry (list > (make block (center 0 0 0) (size sx infinity infinity) > (material Ag51)) > (make block (center 1 0 0) (size 2 infinity infinity) > (material air)) > )) > (define-param fcen 0.2) ; pulse center frequency > (define-param df 0.4) ; pulse width (in frequency) > (define-param runtime 200) > (set! sources (list > (make source > (src (make gaussian-src (frequency fcen) (fwidth df))) > (component Hz) > (center 0.0 0 0)) > )) > ;(use-output-directory) > (define-param kmin 0.0) > (define-param kmax 3.0) > (define-param k-interp 10) > (define kpts (interpolate k-interp (list (vector3 0 kmin 0) (vector3 0 > kmax 0 )))) > (define all-freqs (run-k-points runtime kpts)) ; a list of lists of > frequencies > _______________________________________________ > meep-discuss mailing list > meep-discuss@ab-initio.mit.edu > http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss _______________________________________________ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
