Perhaps my answer will not be of any help to you...
But if you want to do the calculation for spherical nanoparticles, have you considered using Mie theory? Mie theory lets you calculate the electric field around the particles quite easily.If you are interested, I have developped a MATLAB program which uses Mie theory to compute the dielectric field around the nanoparticle for different excitation wavelengths. Contact me if you are interested.
But of course, if you want to consider non-spherical nanoparticles, you must rely on numerical methods.
Also, wether using Mie theory or FDTD, be careful which dielectric function you use for silver. There are many reported measurements and they do not all agree well with each other. One dataset could give you much higher field enhancement at the resonance wavelength compared to another. Some of the most used datasets for silver are the ones from Palik and the ones from Johnson and Christie... and they can both give quite different results.
For further readings about that:
- Babar & Weaver, Optical constants of Cu , Ag , and Au revisited, Applied optics, 2015- Jiang et al., Re-evaluation of literature values of silver optical constants, Optics Express, 2015
- McPeak et al., Plasmonic Films Can Easily Be Better: Rules and Recipes, ACS Photonics, 2015