Dear All,
Are there alternative ways for computing the material dispersion curve other
than the one suggested in the material_dispersion.ctl example file? The problem
I am ecountering is as follows:
I recently adapted the Lorentz-Drude model parameters (Rakic et al., Appl. Opt.
vol. 37, No.22, pp. 5271) for use in Meep. I then modified the material_example
code in order to check my adaptation. However, this example causes Meep to hang
on the first call to harminv. By removing the first and last
make-polarizability lines from the gold definition, I got the example to run,
but harminv only returned zeros. I believe this might be related to the values
of delta-epsilon being too large, but these values are set by the
parametrization. Any help will be appreciated. Below is the modified version of
the materials_example.ctl file.
Thanks,
Gustavo
(set! geometry-lattice (make lattice (size no-size no-size no-size)))
(set-param! resolution 20)
(define gold (make dielectric (epsilon 1)
(polarizations
(make polarizability (omega 1.00000e-6) (gamma 0.04272) (delta-epsilon
4.02583e13))
(make polarizability (omega 0.33449) (gamma 0.19425) (delta-epsilon 11.36290))
(make polarizability (omega 0.66898) (gamma 0.27807) (delta-epsilon 1.18364))
(make polarizability (omega 2.39300) (gamma 0.70122) (delta-epsilon 0.65677))
(make polarizability (omega 3.46901) (gamma 2.01016) (delta-epsilon 2.64549))
(make polarizability (omega 10.73590) (gamma 1.78448) (delta-epsilon 2.01483))
)))
(set! default-material gold)
(define-param fcen 1.72)
(define-param df 6.4)
(set! sources (list (make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ez) (center 0 0 0))))
(define-param kmin 0.3)
(define-param kmax 8)
(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
(map (lambda (kx fs)
(map (lambda (f)
(print "eps:, " (real-part f) ", " (imag-part f)
", " (sqr (/ kx f)) "\n"))
fs))
(map vector3-x kpts) all-freqs)
_________________________________________________________________
It’s the same Hotmail®. If by “same” you mean up to 70% faster.
http://windowslive.com/online/hotmail?ocid=TXT_TAGLM_WL_hotmail_acq_broad1_122008_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
