Hi Steven and meep users!
I was trying to simulate a 3D Ag block in the visible and the problem is the
following.
I use a far-UV drude oscillator to make e' and e" almost constant at
frequencies around 2.
Parameters are epsilon=-11, omega=1000, gamma=0.001, delta-epsilon 0.001, so
that for frequency=2 e'=-10 e"=2, corresponding n=0.31 so I set the
Currant=0.05.
Even in the case I set the fwidth=0.05 it's OK without the block and blowup
with it.
Is it because if internal averaging or Currant is too high?..I've read all the
mailing list archive and it seems that I'm not the only one with this problem
but still I'm not satisfied with the eps=-1e10, a want a real metal)
I have tried eps-averaging?=false, and versions 0.10.1 and 0.10
Regards,
Max
the script I ran is the following:
(reset-meep)
(set! eps-averaging? false)
(define-param no-cube? false)
(set! geometry-lattice (make lattice (size 3 3 6)))
(set! geometry
(if no-cube?
(list
(make block
(center 100 100 100)
(size 1 1 1)
;daleko
(material (make dielectric (epsilon -5)
(polarizations
(make polarizability
(omega 1.11) (gamma 0.64) (delta-epsilon 1)))))))
(list
(make block
(center 0 0 0)
(size 1 1 1)
;blizko
(material (make dielectric (epsilon -11)
(polarizations
(make polarizability
(omega 1000) (gamma 1) (delta-epsilon 0.000001)))))))
))
(set! k-point (vector3 0 0 0))
(set! ensure-periodicity true)
(set! pml-layers (list (make pml (direction Z) (thickness 1.0))))
(set! resolution 10)
(set! sources (list
(make source
(src (make gaussian-src (frequency 2)(fwidth 0.05)))
(component Ex)
(center 0 0 -2)
(size 3 3 0))))
(define trans ; transmitted flux
(add-flux 2 1 200
(if no-cube?
(make flux-region
(center 0 0 2) (size 3 3 0))
(make flux-region
(center 0 0 2) (size 3 3 0)))
))
(if (not no-cube?) (load-minus-flux "zero-flux" trans))
(run-sources+ 1
(at-beginning output-epsilon)
(to-appended "ex" (at-every 0.5 output-efield-x)))
(display-fluxes trans)
(if no-cube? (save-flux "zero-flux" trans))
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
