Dear Steven and Meep users, It is a while since I started to find the dispersion of a simple Metal/Insulator structure with Meep. BTW, I've reviewed the mailing list, but unfortunately nobody gives a complete solution to this open question yet. Anyway, My instruction is as follows: 1. Place a metal in Y+ region and dielectric in Y- region 2. Insert PML in Y direction, and Bloch in X direction; in C++: use_bloch(X, beta); where beta is my wavenumber 3. Put 10 Gaussian sources just beneath the M/I interface in dielectric region. Each source has appropriate phase shift. in C++: f.add_point_source(Hz, src, src_loc, polar(1., 2. * pi * beta * x)); //"Hz" is used since I am interested in plasmonic modes 4. Run the simulation till sources are dead and do "Hz" sampling in place of centered source 5. Do harminv peak detection Any idea about why I cant get correct dispersion with this method. I really appreciate if anybody can give me some hints or maybe a sample code which can calculate the foregoing dispersion correctly
Thanks Joe
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
