Just use (set! default-material Ag) Regards,
Charalampos On Mon, Apr 14, 2014 at 3:08 PM, Allen Spencer <falcon_dai...@hotmail.com>wrote: > Hi Charlampos, > > > The version I use is 1.2.1 which just released few days ago. I know how to > define dispersive material in MEEP but what I want is to get the epsilon vs > frequency as calculate in > http://ab-initio.mit.edu/wiki/index.php/Meep_Tutorial/Material_dispersion. > Thanks a lot for helping! > > Regards, > Allen > > ------------------------------ > Date: Mon, 14 Apr 2014 13:24:44 +0300 > Subject: Re: [Meep-discuss] How to get material dispersion for metal > From: iph...@gmail.com > To: falcon_dai...@hotmail.com > CC: meep-discuss@ab-initio.mit.edu > > > First of all, which version of meep are you using? > > Here is my code for dispersive Ag: > ; Parameters > (define a 100e-9) ; define the unit of length in meters > (define cc (/ (* 2 (* pi 3e+8)) a)) ; unit of frequency in meep > ; eV->radpersec > (define (ev->radpersec l) > (/ l 6.5821e-16) ; divide by the reduced Planck's slab > ) > > ; Ag > ; Parameters from http://arxiv.org/pdf/0911.2737.pdf > (define Ag_epsb 3.7) > (define Ag_wp (/ (ev->radpersec 9.2) cc)) > (define Ag_gamma (/ (ev->radpersec 0.02) cc)) > > (define Ag > (make dielectric (epsilon Ag_epsb) > (E-susceptibilities > (make drude-susceptibility (sigma 1.0) (omega Ag_wp) (gamma Ag_gamma) > ) > ) > )) > > Then, write > (material Ag) > > in your geometry definition.. > > Regards, > > Charalampos > > > On Mon, Apr 14, 2014 at 11:12 AM, Allen Spencer <falcon_dai...@hotmail.com > > wrote: > > 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 > > >
_______________________________________________ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss