Dear Meep users, I am a beginner of MEEP. I am wondering to know how to calculate the material dispersion for metallic material like gold and silver as shown in the Tutorial http://ab-initio.mit.edu/wiki/index.php/Meep_Tutorial/Material_dispersion. I tried to use the same method as given in the tutorial but it seems that the source doesn't work in a PEC like material. Here is my code: (set! geometry-lattice (make lattice (size no-size no-size no-size)))(set-param! resolution 200)(set! default-material(make dielectric (epsilon 1)(E-polarizations (make polarizability(omega 1e-20) (gamma 0.038715) (sigma 4.4625e+41))(make polarizability(omega 0.65815) (gamma 3.1343) (sigma 7.9247))(make polarizability(omega 3.6142) (gamma 0.36456) (sigma 0.50133))(make polarizability(omega 6.6017) (gamma 0.052426) (sigma 0.013329))(make polarizability(omega 7.3259) (gamma 0.7388) (sigma 0.82655))(make polarizability(omega 16.365) (gamma 1.9511) (sigma 1.1133)))))(define-param fcen 0.5)(define-param df 0.5)(set! sources (list (make source (src (make gaussian-src (frequency fcen) (fwidth df))) (component Ez) (center 0 0 0)))) (define-param kmin 0.1)(define-param kmax 1)(define-param k-interp 99)(define kpts (interpolate k-interp (list (vector3 kmin) (vector3 kmax)))) (define all-freqs (run-k-points 200 kpts)) ; a list of lists of frequencies 20) (map (lambda (kx fs) (map (lambda (f) (print "eps:, " (real-part f) ", " (imag-part f) ", " (real-part (sqr (/ kx f))) "," (imag-part (sqr (/ kx f))) "\n")) fs)) (map vector3-x kpts) all-freqs) In Aaron Webster's note, meep-fields-analtic-chi1 has been used but in the newest version, I didn't find it. Can anyone tell me how to solve this problem? Thanks in advance! Regards Spencer
_______________________________________________ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
