2010/3/2, Renaud Vallée <[email protected]>: > Dear Steven, everybody, > > I want to obtain the transmission/reflection spectra of a single gold > spherical nanoparticle located at the center of the simulation box. > Everything works fine except for the polarizability, which is not taken > into account. I use the Drude-Lorentz model described in Vial et al., > PHYSICAL REVIEW B 71, 085416 (2005͒) with units adapted to my case. The > functions sigma_1 and sigma_2 here below are defined accordingly, and > the add-polarizability terms are called in the program. However, they > are not taken into account in the calculation, i.e. if I remove the > add-polarizability terms in the program, everything works the same way. > Owing to the fact that the same program written with the libctl > interface works perfectly, i.e. it reproduces adequately the expected > plasmon resonance peak(s), I guess something is wrong with my C++ > implementation. Does anybody know why? The way the sigma_1 and sigma_2 > functions are defined is not correct or is it due to the call of the > add-polarizability that is not fine?
Check your coordinates origin - they are different for libctl (center of volume) and c++ (corner of volume). Hope it helps, Shavkat Nizamov _______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
