This describes some of the potential reasons for your problem: https://meep.readthedocs.io/en/latest/FAQ/#why-are-the-fields-blowing-up-in-my-simulation
Try bumping the resolution (or drop the Courant factor) and be sure to use absorbing layers in place of pml if you use dispersive cladding. Make sure you really do need dispersive SiO2 though... usually you don’t since the dispersion is so small. I only see 1% relative differences in my designs... It all depends on your application. > On Sep 4, 2020, at 8:18 AM, [email protected] wrote: > > Dear all, > > first of all: Thanks for you kind support! > > I am trying to simulate a simple Bragg reflector made up of double layers of > Si3N4/SiO2using a Guassian source centered around 1550 nm (1/0.64). > (Meep/Python) > > When simply using fixed epsilon values (e.g.: material = mp.Medium(epsilon = > 1.989*1.989)) everything works fine. > > When instead - which is what I would prefer to do - referring to Meep's > material list for either SiO2 or Si3N4 (e.g. material = Si3N4_NIR) or both, > the simulation keeps aborting after some time steps with an error message > referring to the code line containing the "sim.run(...)"-command, and giving > the error message: "RuntimeError: meep: simulation fields are NaN or Inf". > > I have been scratching my head for some time now, so I was wondering if you > might give me a hint as to what might be going wrong here? > > Thanks in advance! > > Best > Sebastian > > P.S. I could imagine that this topic has been raised in the past already, > unfortunately I didn't succeed in accessing the meep-discuss > archive._______________________________________________ > meep-discuss mailing list > [email protected] > http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
