On Fri, 24 Nov 2006, Sung Park LEE wrote:
According to the paper (J. Appl. Phys., Vol. 92, pp.654). When r' = 0.25a, They get monopole, dipole, quadrupole, and hexapole. In my calculation, I get same result when r' = r = 0.35a; I don't understand my results are not correct when reducing the radius of nearest neighbor holes. I think that my results are not correct because my symmetries are not coorect. Could you give me some suggestions about my mistakes?
You shouldn't be using any symmetry at all except for the z=0 mirror plane, if you want to see all of those modes at once. You should also output the epsilon.h5 file [ use (at-beginning output-epsilon) ] and visualize it e.g. with h5topng to make sure the structure looks like what you think it should be (it is easy to get this wrong!).
Most likely there is a mismatch between the structure you are simulating and the one on the paper; read the paper carefully. For example, you say that your slab thickness is 0.5a, but the one reported in the paper is 0.4a. You could also have a bug in your .ctl file that makes the geometry different than what you think, which is why it is important to look at the epsilon.h5 file.
Steven _______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
