I've been working on simulating a non-ideal metal, and it seems to run well with my calculated polarizations if I choose n>=1. However, there seems to be a limit on how low I can make n. At my wavelength of interest, Ag has n=0.226, and I can't go below about n=0.6 without the simulation "blowing up" (energy densities become inf and eventually nan, starting at the metal's outer surface). The lower I make n, the sooner the simulation blows up. I've tried a few geometries with the same effect.
Any ideas? Regards, Bryce _______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
