Dear Steven and other Meep users,
I'm trying to calculate J*E at the source position to invetigate the
decay rate variation (or LDOS) in my cavity structure.
First, I located a point source having a fixed frequency and
polarization Ex and waited for some time until transients died away.
And I tried to get the magnitude of the electric field (Ex) at the
same source position. (I used complex field to get the amplitude of
the electric field at the chosen frequency - not time-varying electric
fields) Then, I ran the simulation again without my structure and
obtained the field value again at the source position to get the
normalized decay rate.
I repeated simulations with different point source positions
(pointsrc_y), but my result is quite different from what I expected
for my structure.. Do you have any ideas or comments for my code?
Thanks very much. Any comment will be really appreciated..
Best Regards,
Young
(set-param! force-complex-fields? true)
(define-param lambda 600)
(define wcen (/ 1 lambda))
(define-param pointsrc_y 100)
...
(set! sources (list
(make source
(src (make continuous-src (frequency wcen)))
(size 0 0)
(component Ex)
(center 0 pointsrc_y))))
(define (print-field-at-point) (print "At y=" pointsrc_y " : " (magnitude
(get-field-point Ex (vector3 0 pointsrc_y))) "\n"))
(run-until (* 10 lambda)
(at-end print-field-at-point))
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss