Hi All,
I am having some issues getting the expected reflectivity spectra from a
simple gold film in air. Below is my ctl file. Essentially, even for a
normalizing run with just air I get unrealistic featues in my specta,
notably a large peak at a frequency of 2 (corresponding to 0.5um in my
base units). I have attempted to increas ethe resolution (upto 100) and
have enlarged the cell size in the Z direction to increase the pml size
(to the level now). This also increased the distance between my source and
structure, but to no avail. It appears from the from the field profile
that I may be getting some form of standing wave formed accros the
periodic boundries?
Any advice would be gratefully recieved.
Thanks in advance,
Kenny
; Meep simulation to compare reflectivity of Gold
(define-param sx 1)
(define-param sy 1)
(define-param sz 8)
(define-param pml-t 1)
(set! ensure-periodicity true)
(set! geometry-lattice (make lattice (size sx sy sz)))
(set! k-point (vector3 0 0 0))
(define-param gold? false)
(define omega_d 7.045)
(define Au
(make dielectric (epsilon 5.967)
(polarizations
(make polarizability
(omega 1e-20) (gamma 5.310E-02)
(sigma (* (* omega_d omega_d) 1e+40)))
(make polarizability
(omega 2.168) (gamma 3.498E-01)
(sigma 1.09)))))
(set! geometry (if gold?
(list (make block (center 0 0 -2) (size sx sy 1)
(material Au)))
(list)))
(set! pml-layers (list (make pml (direction Z) (thickness pml-t))))
(define-param fcen 1.75)
(define-param df 1.5)
(set! sources (list
(make source (src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ex)
(center 0 0 (- (* sz 0.5) (+ pml-t 1))
(size sx sy 0)))))
(set! resolution 40)
(define-param nfreq 100)
(define refl
(add-flux fcen df nfreq
(make flux-region
(center 0 0 0)
(size sx sy 0)
(direction Z))))
(if gold? (load-minus-flux "refl-flux" refl))
(run-sources+
(stop-when-fields-decayed 50 Ex
(vector3 0 0 (* 0.25 sz))1e-3)
(at-beginning output-epsilon)
(to-appended "ex" (at-every 5 output-efield-x)))
(if (not gold?)
(save-flux "refl-flux" refl))
(display-fluxes refl)
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
