-
Notifications
You must be signed in to change notification settings - Fork 28
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Check field computations #24
Comments
I will start with the grating, as that will be the simplest case and I already have a COMSOL model for it. |
Yeah so specifically the z-dependence of the permittivity is not captured well in the E-field computation, and so it causes a discontinuity in a place where there is no interface. This is in particular a problem for field localized in air, as you can see above. I think it could be worth trying to compute E from H(r) and eps(r) by just doing the curl numerically in real space, but I don't have time to experiment with this right now. However, if you do try it out and it looks good compared to other solvers, I will work on adding it to legume. |
Would be good to check the computed fields of the eigenmodes vs. some first-principle simulation.
Note that the H-field should be the most "correct" as all components are continuous across all interfaces. The D and especially the E fields, because of the way the expansion works, will have some spurious discontinuities even in cases when they're not supposed to. For example if a patterned and un-patterned slab of the same material sit on top of each other, there is no "interface" in the region where they have the same permittivity, but because the effective permittivity of the patterned slab is lower, there's a discontinuity in the guided modes, and so there will be one in the PhC modes too.
The text was updated successfully, but these errors were encountered: