Hi,
is there anyone out there who successfully has set up a simulation with
aluminum using the Drude model? I want to compute transmission through a
grid made of aluminum, the Drude model parameters are:
omega_p = 1.8x10^16 rad/s
omega_tau = 1.1x10^15 rad/s
The meep manual suggests to choose some small omega and set
delta-epsilon accordingly.
What I do is:
Define a block:
(make block (center x x1 0)
(size thick B L)
(material (dielectric_drude w1 g1 de))
)
)
The function dielectric_drude looks like:
(define (dielectric_drude w1 g1 de)
(make dielectric (epsilon 1)
(polarizations
(make polarizability
(omega w1)
(gamma g1)
(delta-epsilon de)
)
)
)
)
The parameters are:
(define c 299762458.0)
(define a 500e-9)
(define w0 (/ c a))
; Drude-model parameter
(define w1 (* (/ 1 w0) 1.0e15))
(define de -324)
(define g1 1.8346)
This gives me NaNs, i.e. the output looks like:
#.#
#.#
#.#
#.#
#.#
#.#
#.#
#.#
Changing the parameters to e.g.
; (define w1 1.1)
; (define de 0.5)
; (define g1 1e-5)
give at least numbers (I didn´t check for correctness so far).
Any ideas?
best wishes,
andreas
--
Dr. Andreas Wilde
Gruppenleiter Physikalische Modellierung
Fraunhofer-Institut für Integrierte Schaltungen
Institutsteil Entwurfsautomatisierung
Zeunerstr. 38, 01069 Dresden, Germany
Tel: +49 351/4640-852
mailto:[EMAIL PROTECTED]
_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
