NanoComp / meep

free finite-difference time-domain (FDTD) software for electromagnetic simulations
GNU General Public License v2.0
1.2k stars 612 forks source link

Documentation about materials question #386

Closed kitchenknif closed 6 years ago

kitchenknif commented 6 years ago

The documentation for how meep handles dispersive materials includes a formula for Drude-Lorentz dispersion, that is significantly different from models I've seen in FDTD-related literature.

Equation from your documentation: render2

Equation from, e.g., doi:10.1016/j.cplett.2007.08.027 render1

For the model used and data from dot:10.1364/AO.37.005271, these two models would be equivalent, because conductivity is set to zero. But for fits where the conductivity is not zero, these two models are vastly different.

Am I misunderstanding something here?

stevengj commented 6 years ago

They are inequivalent when you have both σ and Lorentzians, but since both of them have to be fit to experimental data and both have the same asymptotic behaviors it doesn't seem like a problem.

You do have to do a separate fit, though, if you have a fit from the other form.