We will calculate the properties of chromatic plasmonic polarizers for color filtering and polarimetry, reproducing figures 2b, 2d, 3 of Ellenbogen et al., "Chromatic Plasmonic Polarizers for Active Visible Color Filtering and Polarimetry", Nano Lett.12(2):1026-31 (2012) to show the variation in the transmission spectrum as a function of polarization and the position of the dip in transmission as a function of the arm length of the cross-shaped structures of the array.
The file sp_filter.fsp can be used to simulate an array of chromatic plasmonic polarizers (CPPs). The simulation file contains one unit cell of the periodic structure. The polarization angle of the source can be set in the model analysis group which is at the top of the Objects Tree. If set to 0 or 90 degrees, the group will automatically use either symmetric or anti-symmetric boundaries to reduce the simulation volume.
To obtain the polarization information required for figure 2, two simulations are necessary, as the results for arbitrary linear polarizations can be obtained by summing the results from two orthogonal polarizations with appropriate weighting. The simulations are performed by running the parameter sweep named "figure2" in the Optimization and Sweeps window where the transmission spectrum is recorded as a result of the sweep.
To obtain information on the transmission minima as a function of the arm length of the cross structure to recreate figure 3, we use the "figure3" sweep object which will record the transmission minimum as a function of CPP horizontal length.
The results of the sweeps are already contained in the simulation file. However, you can re-run the parameter sweeps to refresh the results. The script file sp_filter.lsf will then analyze the results and plot them as they are displayed in figures 2b, 2d and 3 of Ellenbogen et al.
We can see that there is good agreement with the published results. To improve the results, it would be possible to use a smaller mesh size (here we use 5nm over the Al CPP), replace the ITO with a better material model (here we used a dielectric with n=1.9), and draw the structure in a way that is closer to what has been fabricated. Also, it should be noted that the polarization angle definition was reversed between figures 2b and 2d of the publication.