Abstract: Using Mie theory and Rigorous Coupled Wave Analysis (RCWA) we compare the properties of dielectric particle and void resonators. We show that void resonators—low refractive index inclusions within a high index embedding medium—exhibit larger bandwidth resonances, reduced peak scattering intensity, different polarization anisotropies, and enhanced forward scattering when compared to their particle (high index inclusions in a low index medium) counterparts. We evaluate amorphous silicon solar cell textures comprising either arrays of voids or particles. Both designs support substantial absorption enhancements (up to 45%) relative to a flat cell with anti-reflection coating, over a large range of cell thicknesses. By leveraging void-based textures 90% of above-bandgap photons are absorbed in cells with maximal vertical dimension of 100 nm.
©2011 Optical Society of America
Solar cells utilize textured surfaces to increase light absorption by reducing reflections and increasing the effective optical path length of an incident light ray. In thick solar cells, increased absorption is realized through, for example, deterministic texturing of wavelength-scale pyramidal structures which can be described with geometric optics and ray-tracing techniques [1,2]. Due to texture sizes comparable to or larger than a wavelength, however, these techniques are not applicable to thin film solar cells, which are instead often fabricated on roughened surfaces [3,4]. More recently, researchers have investigated the possibility of coupling photovoltaic (PV) cells to plasmonic [5–7] or dielectric [8,9] resonators. By leveraging similar resonances within the photovoltaic active material, it is also possible to achieve nanoparticle-based textures without coupling to external resonators [10–12].
Typically, the individual scattering element in such designs is a high-index dielectric or negative permittivity plasmonic nanoparticle surrounded by a low-index embedding medium, typically air. In this particle resonator geometry, scattering occurs as a result of the dielectric contrast between air and a dielectric or plasmonic nanoparticle (Fig. 1(a) ). For a given choice of material, however, one can also construct a void resonator, wherein light scatters off nanoscale voids within a high-index dielectric (Fig. 1(b)) or plasmonic medium. Such resonators can form the constituent subunits in photonic crystals  and metamaterials . Recently, solar cell textures based on inclusions of low index materials within dielectric  or plasmonic [16,17] nanovoids have been proposed and demonstrated .
In this paper we compare the generalized scattering properties of single particle and void resonators. Next we propose and evaluate solar cell textures comprising periodic arrays of particle (Fig. 1(c)) or void (Fig. 1(d)) resonators. In comparison with previous investigations regarding anti-reflection effects in dielectric gratings [19–21] we will demonstrate how properties of single resonators are manifested as properties of periodic resonator arrays. In particular, we exploit the scattering properties of particle and void resonators in thin film amorphous silicon (a-Si) solar cell textures. Ultimately, we show that these active material textures can enable significant absorption enhancements over equivalent thin film cells with anti-reflection coatings. We demonstrate absorption of 90% of incident solar photons in 100 nm thick a-Si solar cells that take advantage of void resonator geometries.
2. Scattering properties of single resonators
In both particle and void geometries the scattering properties of individual spherical and cylindrical resonators can be described with analytical Mie theory . While this work will concentrate on cylindrical scatterers, results can be applied to spherical resonators as well. In Mie theory the incident, internal, and scattered fields are decomposed into a set of vector harmonics and which are indexed by their azimuthal phase dependence eimφ. At normal incidence the scattered fields are purely transverse electric (TE) or transverse magnetic (TM) with excitation coefficients am and bm, respectively:
Here nemb and nint are the refractive index of the embedding medium and the resonator respectively, k is the wavevector, r0 is the cylinder radius, and λ is the free space wavelength.
The scattering and extinction properties of single resonators can be determined easily once the Mie coefficients have been calculated. For particle resonators nint > nemb and nr > 1, whereas for void resonators nint < nemb and nr < 1. In Fig. 2 we plot the lowest order (m = 0 or 1) Mie coefficients for particle (red) and void (blue) resonators comprising air and a material with refractive index 4—a value that approximates the refractive index of e.g. silicon or germanium (for a discussion of Mie resonances in spherical Si or Ge particles see e.g [23,24].). A unique feature of the void configuration is a redshift in the TE1 resonance relative to the TE0 resonance. In all other cases resonances shift to higher frequency with increasing mode order (larger m).
At resonance the Mie coefficients equal 1 independent of their geometry, which can be proven analytically . Although the peak values are identical, the void geometry supports larger bandwidth (lower Q) resonances. While small bandwidths are required for applications such as filters and sensors, the large bandwidths available in the inverse resonator geometry may be useful for scattering radiation across the solar spectrum. A more detailed analysis of the resonance bandwidths is given in Appendix A.
2.1 Simplified resonance conditions
In either geometry the various resonance frequencies can be determined analytically by setting the Mie coefficients equal to one , and this condition can be further simplified by taking Mie theory in the limit of large index contrast. In the particle configuration a large index contrast implies nr ≫ 1 and the lowest-order resonances occur for values of x ≪ 1. In the void geometry a large index contrast implies nr ≪ 1 and the lowest-order resonances occur for values of xnr ≪ 1. By taking limits of the Mie coefficients in these respective regimes (see Appendix B for details) we derive simplified resonance conditions, shown in Table 1 .
In both configurations the resonance conditions are functions of the variable 2πr0nH/λ, where nH is the refractive index of the high-index medium, and the relative resonance frequencies are entirely governed by the location of Bessel function zeroes. For instance, the first zeroes of the Bessel functions J0, J1, Y0, Y1, and Y2 occur when the arguments equal 2.40, 3.83, 0.89, 2.20, and 3.38, respectively. From a comparison of these values the red-shift of the TE1 void resonance relative to the TE0 void follows directly. Similarly, the derived resonance conditions capture a feature that is evident in comparing the particle and void Mie coefficients in Fig. 2: for a given choice of two materials the void geometry will exhibit a red-shifted TE1 but blue-shifted TE0, TM0, and TM1 resonance frequency (r0 λ−1).
2.2 Scattering cross section
The scattering cross section (per unit length) for a TE illuminated cylindrical resonator is Fig. 3(a) we plot the scattering cross sections for cylindrical particles and voids. Scattering from particles results from a sequence of narrow resonances, whereas voids have significantly smaller variations with normalized frequency due to strongly overlapping broadband resonant modes. Additionally, very small voids scatter more TE radiation than their particle counterparts, whereas the opposite is true for TM polarization. Regardless, the larger bandwidths of void resonators is offset by the reduction in scattering due to large nemb and the particle resonators have a slightly larger unpolarized scattering cross section over most normalized frequencies.
2.3 Forward scattering properties
Although the total cross section may be slightly weaker for voids, for the solar cell textures in Fig. 1 the principle interest is in forward scattering, where the effect of an interface on scattering properties is very important. Previously, researchers noted that classical dipoles at an interface emit preferentially into the high index material and demonstrated the importance of this phenomenon in plasmon enhanced solar cells . In Fig. 3(b) we plot the fraction of light emitted into a substrate for both parallel (blue) and perpendicular (red) dipoles embedded within a low-index (n = 1, dashed) or high-index (n = 4, solid) medium as a function of distance from the interface. The fraction of forward (into the substrate) emitted radiation decreases rapidly with distance for dipoles sitting above a high-index substrate. This situation approximates scattering in the particle geometry and necessitates careful design of particle resonators that are located within the interface near-field  in order to suppress back-scattering. In contrast, a dipole embedded within the high-index substrate emits more than 98% of its radiation in the forward direction regardless of orientation or distance from an interface. This situation is akin to the void geometry, and suggests that void resonators may exhibit substantial preference for forward scattering.
2.4 Single particle summary
We have demonstrated the possibility of designing resonators in either a void or particle configuration for any given choice of two materials. In addition to deriving resonance conditions for the lowest order Mie resonances we have demonstrated a systematic comparison of the two scattering geometries. Both resonator types display an increase in scattering with normalized frequency up to the first Mie resonance, after which the scattering is mostly constant for the void geometry while the scattering by particles is dominated by narrow resonance peaks. Compared to particles, ultra-small voids scatter TE radiation more efficiently, while the opposite is true for TM polarization. Voids tend to exhibit larger bandwidths, have smaller unpolarized cross-sections, but are likely to scatter a larger portion of light in the forward direction in the presence of an interface. In the subsequent section we demonstrate how these various scattering properties are evident in dielectric resonator array solar cell textures and evaluate the use of both void and particle resonators for enhancing light absorption in solar cells.
3. Resonator arrays
In Fig. 1 we present two different solar cell textures comprising arrays of rod-shaped particle (Fig. 1(c)) or void (Fig. 1(d)) resonators. Despite significant differences between single scatterers and resonator arrays—square rather than circular cross sections, the presence of an interface, and periodicity dependent coupling between resonators —we will demonstrate that single resonator properties are preserved in solar cell textures. Based on previous investigations of absorption in particle rods, we anticipate that the essential properties of lowest order Mie resonances are unchanged despite the differing cross-sections . The electromagnetic response of arrays was calculated with a rigorous coupled wave algorithm (RCWA) implemented via DiffractMod from RSoft . In RCWA reflected and transmitted light is decomposed into a basis set of diffracted plane-wave components :Fig. 4(a) we plot the scattering properties of textures based on the designs in Fig. 1 for a substrate index of n = 4 but without an ARC. Simulations are presented for a periodicity 40 times larger than the resonator diameter such that inter-particle coupling is weak and there is a near continuum of diffracted orders kxj. The sum of the power contained in all diffracted orders is a good approximation to the single resonator scattering cross-sections, and allows for a direct comparison with Fig. 3(a).
The general results derived previously for a single resonator can be observed in the case of resonator arrays as well. Both resonator geometries display an increasing scattering cross-section with normalized frequency up to the first resonance, after which the scattering cross-section plateaus. Ultra-small void-based arrays scatter TE radiation more efficiently than the particle counterparts, while the reverse is true for TM polarization. Void arrays have larger bandwidths but less efficient scattering, especially for TM polarization.
By dividing the power in the transmitted orders by the total scattering we can determine the proportion of forward scattering (Fig. 4(b), bottom) and compare with dipole calculations (Fig. 3(b)). The comparison is most valid at small normalized frequencies, where the scattering arises from single dipolar resonances. As particles get larger, and their center moves away from the interface, the proportion of forward scattered light is reduced dramatically.
Furthermore, voids are superior forward scatterers, as expected from our dipole calculations. In the top of Fig. 4(b), however, we show that back-scattering can be reduced for particles by adding an anti-reflection coating; light scattered into oblique angles is totally internally reflected at the boundary of the air-ARC boundary and is redirected into the substrate.
While this analysis was presented for large periodicity arrays a practical solar cell texture requires a substantial fill fraction to maximize light scattering. In Fig. 5 we plot the normalized scattering intensity for arrays of 100nm diameter particles, as a function of fill fraction (i.e. periodicity). While the scattering amplitude depends on fill fraction (Figs. 5(c) and 5(d)), the frequency dependence (Figs. 5(a) and 5(b)) is nearly independent of periodicity for fill fractions less than approximately 45%. These results indicate that the properties of dielectric textures are primarily governed by the Mie resonances of individual resonators located within an inhomogeneous environment.
4. Thin film absorption enhancements
In previous sections we detailed the optical response of single particle and void resonators and showed that the response of resonator arrays is primarily determined by single resonator properties (at fill fractions less than 50%). In real solar cells, with finite thickness, scattering off surface texture enhances light absorption through light trapping. Light trapping effects are particularly important in thin-film cells where the single-pass absorption of a normal incident light ray is weak. In this section we will demonstrate significant absorption enhancements in resonator arrays atop very thin a-Si layers. Absorption enhancements may arise from anti-reflection effects, excitation of lossy Mie resonances, and coupling to waveguide modes. As these mechanisms are interrelated, finding an optimal cell design is non-trivial.
We use a numerical optimization algorithm, Covariance Matrix Adaptation Evolution Strategy (CMA-ES) , to efficiently search large areas of parameter space for the global maximum in absorbed photons. We calculate the total fraction of absorbed photons by integrating the product of the solar spectrum, bAM1.5(λ), and the simulated absorption spectrum of the cell, A(λ), for photons with energy larger than the a-Si mobility gap (730nm):
We calculate the absorption in the amorphous silicon (a-Si) layer using a dielectric function from , and we model the back contact as a perfect electrical conductor (PEC) to avoid complexities associated with small but finite absorption effects in real metals. For the ARC we assume an index of 2, which approximates the refractive index of ZnO, ITO and Si3N4 in the optical range. For every value of the maximum a-Si thickness t, we allow the square rod diameter d, the periodicity a, and the ARC coating thickness tARC to vary (Fig. 6(c) ). The optimized absorption values as a function of t are compared to an equivalent maximum thickness flat cell with optimized ARC and plotted in Figs. 6(a) and 6(b).
Both particle and void resonator arrays absorb significantly more solar photons than the untextured cell for all thicknesses, with enhancement factors as large as 45%. For thin cells between 0 and 50 nm and 75-112 nm the void geometry has superior absorption properties. In particular, 100 nm thick solar cells with an array of void resonators already absorb over 90% of the incoming solar photons, more than the absorption of a 300 nm thick flat cell. For cells thicker than ~140 nm the particle geometry absorbs more light than the void geometry.
For these thin cells the maximal possible rod diameter is small and we expect the absorption properties to be similar to the scattering properties for small resonators described in sections 2 and 3. As expected, in this regime the void (particle) resonator geometries absorb TE (TM) light more efficiently (Fig. 6(d)), especially at long wavelengths. While future investigations will focus on detailing the interplay of Mie resonances and waveguide modes in these structures, these results clearly demonstrate that leveraging Mie resonance effects within patterned grating structures can enhance light absorption in thin solar cells of varying thicknesses.
In this study we demonstrate a number of significant differences in the optical properties of dielectric particle and void resonators. Using Mie theory we derive simplified resonance conditions for the lowest order resonances in both configurations in the limit of large index contrast. Voids resonate at larger normalized frequencies than particles, with the exception of the TE1 resonance, which shows an unusual redshift, and have higher bandwidth resonances except for the TM0 mode. Due to the early onset of the TE1 mode ultrasmall voids are more efficient TE scatterers than ultrasmall particles, while for TM polarization the opposite is true. In general, however, the advantages from larger bandwidth resonances in the void configuration are partially offset by an inverse relationship between peak scattering cross-section and embedding medium refractive index.
We suggest two distinct active medium solar cell textures which rely on void or particle resonances to enhance light absorption. Using RCWA calculations we show that void geometries have preferential scattering into the active medium as expected from calculations of single dipole emitters. We show that the particle and void solar cell textures preserve many of the general features evident from Mie theory and calculate optimized absorption values for a-Si solar cells of varying thicknesses which exploit these textures. Enhancement factors as large as 45% are observed, and void-based textures allow for 90% absorption of incident above-bandgap solar photons in cells whose maximal thickness is 100 nm, a three-fold reduction in active layer thickness compared to untextured cells with equivalent absorption. We anticipate that single particle calculations will bring new understanding of basic properties of void resonators and resonator arrays. The results of this study may lead to novel solar cell textures and may be used to refine understanding of other void-based solar cell textures  (see Appendix C).
Appendix A: Bandwidths
For solar cell applications, where the illuminating radiation is broadband, the bandwidth of resonators is a very important quantity. An examination of the Mie coefficients in Fig. 2 suggests that the peak values of the Mie coefficients are identical, but the void geometry supports larger bandwidth (lower Q) resonances. This trend becomes more evident as the index contrast between the particle and embedding medium is increased.
In the particle geometry, as the index contrast increases the resonator size gets smaller and Q increases exponentially . With voids, however, the Q factor is independent of nr. This fact can be proven by deriving simplified Mie coefficients in the limit of small nr, as we do in Appendix B. The Mie coefficients depend only on the normalized frequency x, thus the resonance lineshape of a void resonator looks identical for all values of index contrast where the approximation holds. The possibility to have very small voids with large bandwidths is markedly different than what occurs in the particle geometry for the same choice of materials. In both cases the resonator size relative to the free space wavelength gets smaller with increasing index contrast. However, in the void geometry the size relative to the embedding medium wavelength is constant and it is this normalized size that determines the bandwidth.
Differences in bandwidth between the two scattering geometries also increase markedly with higher azimuthal mode index n. In the void configuration the Q factor is nearly constant with m whereas in the particle configuration the Q factor increases exponentially with m (Fig. 7 ) a scaling which has been observed in low order whispering gallery microcavities . In the particle geometry a light-ray within a high index medium encounters a low-index medium, and at high m the cylindrical boundary looks approximately planar leading to near-total internal reflection and high Qs. In voids, the light-ray approaches the boundary from the low-index side, there is finite transmission even for a planar interface, and bandwidths remain large.
Appendix B: Simplified Mie coefficients and resonance conditions
Below, we derive the large index contrast limits for particle and void resonators presented in Table 1. For void resonators we can take the limit of the Mie coefficients as nr → 0:
The Mie coefficients become functions solely of the variable x = 2πr0nemb/λ and can be solved directly to determine the resonance values in Table 1. The simplified Mie coefficients are reasonably accurate starting from index contrast of approximately 3, and are accurate even up to very high values of the mode index m. These simplifications highlight the fact that in the high index contrast limit the void resonance lineshapes, and thus bandwidths, are independent of the specific value of nint and remain quite large even for resonant voids that are significantly smaller than the freespace wavelength.
These properties of voids are in stark contrast with particle resonators, where the bandwidth depends strongly on the particle refractive index nint. For particles, we were unable to find a simplified limit of the Mie coefficient as nr → ∞. To determine the resonance conditions we set the Mie coefficients equal to unity , which implies that
The Bessel functions have arguments of both x and nrx and the resonances cannot be solved for arbitrary values of nr. However in the large index contrast limit (nr → ∞) the resonance condition for TE polarization simplifies to Jm(nrx) = 0. For the lowest order TM resonances, the limit of large index contrast also implies that x → 0 and with suitable Taylor expansions about x near 0 we can derive simplified resonance conditions:
We would like to refer the reader to chapters 10 and 15 in  for an interesting discussion of the resonant condition in the limit of very large nr.
Appendix C: Electric field intensity within void resonators
In reference  Yu et al. suggest enhancing light absorption in low index materials by embedding them within a high index medium. For deeply subwavelength spherical voids of low refractive index (nL) surrounded by a high refractive index medium (nH) they derive a maximal absorption enhancement factor (Eq. (20)) of
This equation can be rewritten as
Note that this generalized form applies to planar inclusions as well (Eq. (17)). The enhancement factor depends directly on the average electric field intensity within the void. For a deeply subwavelength inclusion the internal field can be derived from electrostatic considerations. In Fig. 8 we show the average squared electric field enhancement for a spherical void of index nint=1 surrounded by a medium of index nemb=4 as a function of its normalized frequency (its size relative to the incident wavelength). In addition to recovering the electrostatic result, , this result highlights the fact that the internal field drops off quickly as the particle size increases. To fully exploit the spherical void geometries it is important to construct inclusions that are significantly smaller than the wavelength of the illuminating radiation.
We would like to thank R. Zia for help with dipole emitter calculations. This material is based upon work supported as part of the Center for Re-Defining Photovoltaic Efficiency through Molecule Scale Control, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences under Award Number DE-SC0001085.
References and links
1. P. Campbell, “Enhancement of light absorption from randomizing and geometric textures,” J. Opt. Soc. Am. B 10(12), 2410–2415 (1993). [CrossRef]
2. M. A. Green, “Two new efficient crystalline silicon light-trapping textures,” Prog. Photovolt. Res. Appl. 7(4), 317–320 (1999). [CrossRef]
3. O. Isabella, K. Jager, J. Krč, and M. Zeman, “Light scattering properties of surface-textured substrates for thin-film solar cells,” Proceedings of the 23rd European Photovoltaic Solar Energy Conference and Exhibition (EU PVSEC), (2008), Session 3AV 1, pp. 476–481.
4. J. Müller, B. Rech, J. Springer, and M. Vanecek, “TCO and light trapping in silicon thin film solar cells,” Sol. Energy 77(6), 917–930 (2004). [CrossRef]
6. V. E. Ferry, M. A. Verschuuren, H. B. T. Li, E. Verhagen, R. J. Walters, R. E. I. Schropp, H. A. Atwater, and A. Polman, “Light trapping in ultrathin plasmonic solar cells,” Opt. Express 18(S2Suppl 2), A237–A245 (2010). [CrossRef] [PubMed]
8. Y. A. Akimov, W. S. Koh, S. Y. Sian, and S. Ren, “Nanoparticle-enhanced thin film solar cells: Metallic or dielectric nanoparticles?” Appl. Phys. Lett. 96(7), 073111 (2010). [CrossRef]
9. J. Grandidier, D. M. Callahan, J. N. Munday, and H. A. Atwater, “Light absorption enhancement in thin-film solar cells using whispering gallery modes in dielectric nanospheres,” Adv. Mater. (Deerfield Beach Fla.) 23(10), 1272–1276 (2011). [CrossRef] [PubMed]
10. M. Kroll, S. Fahr, C. Helgert, C. Rockstuhl, F. Lederer, and T. Pertsch, “Employing dielectric diffractive structures in solar cells - a numerical study,” Phys. Status Solidi 205(12), 2777–2795 (2008). [CrossRef]
11. L. Cao, P. Fan, A. P. Vasudev, J. S. White, Z. Yu, W. Cai, J. A. Schuller, S. Fan, and M. L. Brongersma, “Semiconductor nanowire optical antenna solar absorbers,” Nano Lett. 10(2), 439–445 (2010). [CrossRef] [PubMed]
12. S. Bandiera, D. Jacob, T. Muller, F. Marquier, M. Laroche, and J. J. Greffet, “Enhanced absorption by nanostructured silicon,” Appl. Phys. Lett. 93(19), 193103 (2008). [CrossRef]
13. J. D. Joannopoulos, S. G. Johnson, J. N. Winn, and R. D. Meade, Photonic crystals: molding the flow of light, (Princeton Univ Press, 2008).
14. F. Falcone, T. Lopetegi, M. A. G. Laso, J. D. Baena, J. Bonache, M. Beruete, R. Marqués, F. Martín, and M. Sorolla, “Babinet principle applied to the design of metasurfaces and metamaterials,” Phys. Rev. Lett. 93(19), 197401 (2004). [CrossRef] [PubMed]
15. Q. G. Du, C. H. Kam, H. V. Demir, H. Y. Yu, and X. W. Sun, “Enhanced optical absorption in nanopatterned silicon thin films with a nano-cone-hole structure for photovoltaic applications,” Opt. Lett. 36(9), 1713–1715 (2011). [CrossRef] [PubMed]
17. T. V. Teperik, F. J. García de Abajo, A. G. Borisov, M. Abdelsalam, P. N. Bartlett, Y. Sugawara, and J. J. Baumberg, “Omnidirectional absorption in nanostructured metal surfaces,” Nat. Photonics 2(5), 299–301 (2008). [CrossRef]
18. N. N. Lal, B. F. Soares, J. K. Sinha, F. Huang, S. Mahajan, P. N. Bartlett, N. C. Greenham, and J. J. Baumberg, “Enhancing solar cells with localized plasmons in nanovoids,” Opt. Express 19(12), 11256–11263 (2011). [CrossRef] [PubMed]
19. T. K. Gaylord, W. E. Baird, and M. G. Moharam, “Zero-reflectivity high spatial-frequency rectangular-groove dielectric surface-relief gratings,” Appl. Opt. 25(24), 4562–4567 (1986). [CrossRef] [PubMed]
21. A. Gombert, K. Rose, A. Heinzel, W. Horbelt, C. Zanke, B. Bläsi, and V. Wittwer, “Antireflective submicrometer surface-relief gratings for solar applications,” Sol. Energy Mater. Sol. Cells 54(1-4), 333–342 (1998). [CrossRef]
22. C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley Inter-Science, 1998).
23. A. García-Etxarri, R. Gómez-Medina, L. S. Froufe-Pérez, C. López, L. Chantada, F. Scheffold, J. Aizpurua, M. Nieto-Vesperinas, and J. J. Sáenz, “Strong magnetic response of submicron silicon particles in the infrared,” Opt. Express 19(6), 4815–4826 (2011). [CrossRef] [PubMed]
24. R. Gómez-Medina, B. García-Cámara, I. Suárez-Lacalle, F. González, F. Moreno, M. Nieto-Vesperinas, and J. J. Sáenz, “Electric and magnetic dipolar response of germanium nanospheres: interference effects, scattering anisotropy, and optical forces,” J. Nanophotonics 5(1), 053512 (2011). [CrossRef]
26. K. R. Catchpole and A. Polman, “Design principles for particle plasmon enhanced solar cells,” Appl. Phys. Lett. 93(19), 191113 (2008). [CrossRef]
27. A.B. Evlyukhin, C. Reinhardt, A. Seidel, B.S. Luk’yanchuk, and B.N. Chichkov, “Optical response features of Si-nanoparticle arrays,” Phys. Rev. Lett. B 82, 045404 (2010).
28. RSoft Design Group, Inc., http://www.rsoftdesign.com.
29. M. G. Moharam, E. B. Grann, D. A. Pommet, and T. K. Gaylord, “Formulation for stable and efficient implementation of the rigorous coupled-wave analysis of binary gratings,” J. Opt. Soc. Am. A 12(5), 1068–1076 (1995). [CrossRef]
30. The algorithm can be found at http://www.lri.fr/~hansen/cmaes_inmatlab.html#python.
31. J. Leng, J. Opsal, H. Chu, M. Senko, and D. E. Aspnes, “Analytic representations of the dielectric functions of materials for device and structural modeling,” Thin Solid Films 313–314(1-2), 132–136 (1998). [CrossRef]
33. H. C. van der Hulst, Light scattering by Small Particles (Dover, 1981).