Dear MEEP users,
I am new to meep. I have developed a code for a silver
nano sphere. When i set the PML value to 0.01 or 0.001 then i got Quadrapole
mode in the nano sphere but at the same value of PML i got no value for the
Qaulity factor. And when i set the PML to 0.1 or 1, then i got no mode
inside the sphere and also no Q-factor. I have add a code with the email.
Any help would be greatly appreciated.
(set-param! resolution 76)
(set! Courant 0.09)
(define-param n 0.1534) ; for 447nm wavelength
(define-param k 2.44420)
(set! geometry-lattice (make lattice (size 0.6 0.6 no-size)))
(set! eps-averaging? false)
(set! geometry
(append
(list
(make block (center 0 0 0)(size 0.6 0.6 0.6)
(material air))
(make sphere (center 0 0 0) (radius 0.075)
(material (make medium (index n)
(E-polarizations
(make polarizability
(omega 1e-20) (gamma 0.2416) (sigma 3.512e41)))
)))
)))
(set! pml-layers (list (make pml (thickness 0.01))))
(define-param nfreq 500)
(define-param fcen 2.24)
(define-param wave 0.44733)
(define-param df 3)
(set! sources
(list (make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
;(src (make gaussian-src (wavelength wave) (width df)))
(component Hz)
(center -0.15 0 0) (size 0 0.48 0)
)))
;(run-until 200
;(at-beginning output-epsilon)
;(to-appended "hz" (at-every 1 output-hfield-z)))
(run-sources+ 2000
(at-beginning output-epsilon)
(after-sources (harminv Hz (vector3 0 0 0) wave df)))
Best Regards,
Khan
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
