mpb in meep with dispersive materials

[ryouta1120]

I'm interested in photonic crystals, so I'm considering using mpb packaged in meep.
Here, I have a simple question.
Is it available to apply dispersive Drude-Lorentz metals on mpb?

According to the discription of mpb, I guess it only allows to specify as constant.
However, testing some materials in materials.py, I could run mpb.modesolver without any error.
My environment is built by anaconda.

[oskooi]

Is it available to apply dispersive Drude-Lorentz metals on mpb?

No. As stated in the "Features" list of https://mpb.readthedocs.io/en/latest/:

Support for arbitrary, anisotropic dielectrics including gyrotropic/magneto-optic materials 
and non-orthogonal unit cells. Lossy and wavelength-dependent ε and μ are not supported.

Note, however, that Meep (not MPB) has a frequency-domain eigensolver which you can use for structures involving metals: https://meep.readthedocs.io/en/latest/Eigensolver_Math/. This feature has not yet been thoroughly tested and so convergence may be an issue.

[ryouta1120]

Thank you for replying.

Now in my case, is it the best way to use Hatmienv following band diagram tutorial?

[smartalecH]

You could also try the frequency-domain shift-invert eigensolver within Meep. The manual description is here. There's an example script here.

You can use PML, Lorenz-Drude materials, etc.

If you just care about generating band diagrams, however, Harminv is probably the way to go. It will be significantly faster and more efficient.

[ryouta1120]

Thank you for your suggestion.

I will try both.