Dear Steven and other users,
As we all known Green’s function has a singularity at r=r’, so if meep
calculate LDOS from the trace of the Green's function, can I get vacuum LDOS
from following codes? The analytical solution of vacuum LDOS is omega^2/(pi^2*
c^3).The result of my code seems not right. I want to compare simulation
results and analytical solutions in order to simulate other model’s LDOS. How
can I get vacuum LDOS?
My code:
(set-param! resolution 30)
(set! Courant 0.1)
(define-param sxy 20)
(define-param sz 5)
(define-param dpml 1)
(define-param sXY (+ sxy (* 2 dpml)))
(define-param sZ (+ sz (* 2 dpml)) )
(define-param z (/ sz 2))
(set! geometry-lattice (make lattice (size no-size no-size sZ)))
(set! force-complex-fields? true)
(set! k-point (vector3 0 0 0))
(set! ensure-periodicity true)
(set! pml-layers (list (make pml (direction Z) (thickness dpml))))
(set! geometry (list
(make block (center 0 0 (/ z 2)) (size sxy sxy z) (material air))
))
(define-param fcen (/ 1 3))
(define-param df 0.1)
(set! sources (list (make source
(src (make gaussian-src (frequency fcen) (fwidth df))) (component Ez) (
center 0 0 1))))
(run-sources+ 5000 (dft-ldos fcen df 30))
Result:
ldos0:, 0.2833333333333333, 0.014037260731577704
ldos0:, 0.2867816091954023, -0.006127334388567867
ldos0:, 0.29022988505747127, 0.0017494695483672368
ldos0:, 0.29367816091954024, 1.0749165090993693e-4
ldos0:, 0.2971264367816092, -6.747295758148794e-4
ldos0:, 0.3005747126436782, 6.79867221417562e-4
ldos0:, 0.30402298850574716, -4.752211708063314e-4
ldos0:, 0.30747126436781613, 2.5390916705574356e-4
ldos0:, 0.3109195402298851, -8.54016430309763e-5
ldos0:, 0.3143678160919541, -2.6226179435950882e-5
ldos0:, 0.31781609195402305, 8.402467557645676e-5
ldos0:, 0.321264367816092, -9.851241374346014e-5
ldos0:, 0.324712643678161, 8.673338269835055e-5
ldos0:, 0.32816091954022997, -5.8283507567008535e-5
ldos0:, 0.33160919540229894, 1.9745369321159384e-5
ldos0:, 0.3350574712643679, 1.9375036796669827e-5
ldos0:, 0.3385057471264369, -5.337955212916454e-5
ldos0:, 0.34195402298850586, 7.8557271588279e-5
ldos0:, 0.34540229885057483, -8.619056960090422e-5
ldos0:, 0.3488505747126438, 6.883483910166926e-5
ldos0:, 0.3522988505747128, -2.238868790286368e-5
ldos0:, 0.35574712643678175, -6.246056334498924e-5
ldos0:, 0.3591954022988507, 1.8821578338003457e-4
ldos0:, 0.3626436781609197, -3.327395108874167e-4
ldos0:, 0.36609195402298866, 4.5442497355070765e-4
ldos0:, 0.36954022988505764, -4.509844620586398e-4
ldos0:, 0.3729885057471266, 7.372528219273671e-5
ldos0:, 0.3764367816091956, 0.0010519713257993227
ldos0:, 0.37988505747126455, -0.0034275621325830763
ldos0:, 0.3833333333333335, 0.0076687403827055065
I will be thankful if anybody helps me.
Mayawei
_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
