Hello, It appears that the flux plane in the sample code is too close to the monolayer structure. Figure that the sphere radius is 0.3, making the diameter and thus the monolayer thickness 0.6, and the flux plane is at 1.0. I think there's also some confusion as the example is titled "Transmission" but with the flux plane being positioned above the monolayer (i.e. same side as the source) it is really accumulating reflected flux. The example states "The source plane and flux plane are positioned at the top and bottom of the computational cell respectively." However, the bottom of the cell (0,0,0) is along the bisection of the monolayer. If tranmission flux were to be collected, the flux plane would need to have a negative z (il.e. under the monolayer). In the example code, I changed the size of the computational cell to 12, and z coordinate of the source and flux planes to 10 and 8, respectively, and the results were normalized to 1 as they should be. Note the values come out to be negative, because the reflected flux is in the opposite direction. Just negate the values to get the reflection, or add 1 to the negative values to get the transmission. Both should be normalized have values between 0 and 1, with the exception of a very spurious points (can be removed with optimization).
- Ian _______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
