This article is about simulating localized plasmon resonances in metal nanospheres using MEEP package. Generally, I am interested in solving three problems in LSPR systems:

**Calculate the extinction, scattering, absorption spectra of metal nanoparticle**

The procedure for doing this is very similar to the method I mentioned here.**Calculating the electric field enhancement spatially as function of wavelength**

This involves taking electric field distributions with a particle in time domain and taking FFT of them. Also to be noted is that the electric fields near the particle should be normalized with electric fields with no nanoparticle. This has to be done by an external program like octave/matlab. One has to do this for three components of electric field in two planes. This was the tough part to pull off.**Calculate electric field at plasmon resonance as function of time**This was very easy. Get the electric fields at certain wavelength and then multiply it with a time varying sin function.

Below are some results for a silver sphere of 25 nm radius. I am happy with these results. They seem to match the analytical results. Some of the animations are large in size and may take time to load on your browser.

Update1: You can find my project file HERE. Enjoy!

(please remember that the project file is designed to run on ubuntu like system with meep and octave installed. The shell script does everything needed to get the results, make sure to “chmod” the shell file, so that you can run the file. In addition, also make sure to increase the resl (resolution) from 10 to 50 for decent results. This means it will take longer time to run the simulation. The results you see in this post were generated by resl=100 and it took 2 days to run on my computer. If you find any bugs, please let me know.)

Update2: This method also works well for gold nano-particles of 50 nm radius. See the last figure in above album.

Nice work! How long did you run the simulation of e. g. the 5th image? If it is full 3D, it has to consume really a lot of CPU time in such a good resolution!

It took approximately a day and half on my Intel Xeon(R) CPU E5620 @ 2.40 GHz*8, 60 Gb ram computer. For half the resolution it took around 5 hours.

Very interesting post. I have been looking for a similar example for some time. Could you please send me the code for these calculations?

Hi Dear,

I would like to reproduce the result in 5 image (for gold nanoparticle), for this aim i replaced only the Ag particle in your project file with the Au one but i’m not getting the right result.

I appreciate any help.

With regards

Beautiful! I have been working on MEEP calculated extinction, scatering and Absorption spectra of gold and silver nanoparticles using the same setup and a gaussian source. Nevertheless, in the long wavelength region I got instability, i.e., the intensity of each spectra oscillates and increase with the wavelength. Could you (or someone) please help me to know how to fix it? I will be very grateful for your help.

With regards

Your work is inspirational! You’v made a massive contribution to science! Thanks alot!

So what are some reasonable and extreme real EM field values.. like lets say from “sun light” on “good dots” to “extreme UV laser” on ” perfectly resonant dipole-dipole Q dots”? familiar units would help us alot.. like Volts/m3 or even better gauss… have u read about the insane “pseudo” magnetism in graphene – bent.. = 200-600 tesla!!!!!!!