On Feb 11, 2008, at 1:19 PM, Seong Kyu Kim wrote: > (1) dwpr is supposed to write E*D/2 or epsilonE*E, whose time > average under vacuum should give 1/4. However, the time averaged > dwpr gives almost zero in vacuum area? (I get the same result if I > accumulated many h5 files and averaged them individually.) Isn't > this strange? > Why do you think this is strange? It's obviously a positive quantity; the absolute magnitude will depend on your source etc. > (2) The time-averaged dwpr shows negative values inside the gold > target, which is understandable since epsilon of gold > > = -22+1.48 i at the calculated frequency. Then, what is the value of > epsilon that calculates dwpr. Is it -22 in this case? > Yes, it should be if you have waited long enough to get a constant- frequency (harmonic) field (and you calculated things correctly).
> > (3) The output-epsilon shows the value of epsilon at the gold region > is 5.967. How did this value come out? > This is just the non-frequency-dependent part of the dielectric constant (the dielectric constant at frequency = infinity) that you specified. Steven _______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
