The output-epsilon HDF5 only outputs the non-dispersive dielectric constant (i.e. the \epsilon_\infty term). If you want the frequency-dependendent \epsilon you must compute it yourself from the formula in the manual.
Note that regardless of what is outputted by output-epsilon, the actual simulation uses the full dispersive dielectric function (effectively, via the auxiliary differential equation method). Regards, Steven G. Johnson On Sat, 8 Sep 2007, Jun She wrote: > How to get the value of the epsilon from the dispersive material? > > The refraction index of Ag I used is 0.158+2.365i. I used *Drude model* > and the * polarizability: omega=1e-20, gamma=0.00259181, > deps=3.45229e+37, esat=0 *Then I used get-epsilon-point to get the > epsilon of Ag. However, the value is 1 which is equal to the the > *instantaneous* dielectric function I used in my code. I set the > force-complex-fields as true. The result is the same. _______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss

