Dear Renaud, Could you please provide the contents of your "parameters" file?
Thanks -----Original Message----- From: [email protected] [mailto:[email protected]] On Behalf Of Renaud Vallée Sent: March 2, 2010 10:36 AM To: Nizamov Shawkat Cc: [email protected] Subject: Re: [Meep-discuss] problems with material dispersion of gold (metal) in meep C++ Hello Shavkat, thanks for your reply. Yes, I took care of this, via the lines 'vec xx = p - v.center();' inserted in each of the relevant functions used to define the geometry. Renaud Le mardi 02 mars 2010 à 13:59 +0100, Nizamov Shawkat a écrit : > 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 -- Dr. Renaud Vallée Centre de Recherche Paul Pascal 115, avenue du docteur Albert Schweitzer 33600 Pessac France Tél: +33556845612 _______________________________________________ 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

