A generalized “ cut and projection ” algorithm for the generation of quasiperiodic plasmonic concentrators for high efficiency ultra-thin film photovoltaics

This report will present a generalized two-dimensional quasiperiodic (QP) tiling algorithm based on de Bruijn’s “cut and projection” method for use in plasmonic concentrator (PC) / photovoltaic hybrid devices to produce wide-angle, polarization-insensitive, and broadband light absorption enhancement. This algorithm can be employed with any PC consisting of point-like scattering objects, and can be finetuned to achieve a high spatial density of points and high orders of local and long-range rotational symmetry. Simulations and experimental data demonstrate this enhancement in ultra-thin layers of organic photovoltaic materials resting on metallic films etched with arrays of shallow subwavelength nanoholes. These devices work by coupling the incident light to surface plasmon polariton (SPP) modes that propagate along the dielectric / metal interface. This effectively increases the scale of light-matter interaction, and can also result in constructive interference between propagating SPP waves. By comparing PCs made with random, periodic, and QP arrangements, it is clear that QP is superior in intensifying the local fields and enhancing absorption in the active layer. ©2013 Optical Society of America OCIS codes: (350.6050) Solar energy; (240.6680) Surface plasmons; (310.6628) Subwavelength structures, nanostructures; (160.4890) Organic materials. References and links 1. H. J. Snaith, L. Schmidt-Mende, M. Chiesa, and M. Grätzel, “Light intensity, temperature, and thickness dependence of the open-circuit voltage in solid-state dye-sensitized solar cells,” Phys. Rev. B 74(4), 045306 (2006). 2. H. A. Atwater and A. Polman, “Plasmonics for improved photovoltaic devices,” Nat. Mater. 9(3), 205–213 (2010). 3. S. Pillai and M. A. Green, “Plasmonics for photovoltaics applications,” Sol. Energy Mater. Sol. Cells 94(9), 1481–1486 (2010). 4. K. Nakayama, K. Tanabe, and H. A. Atwater, “Plasmonic nanoparticle enhanced light absorption in GaAs solar cells,” Appl. Phys. Lett. 93(12), 121904 (2008). 5. V. E. Ferry, M. A. Vershuuren, H. B. T. Li, R. E. I. Schropp, H. A. Atwater, and A. Polman, “Improved redresponse in thin film a-Si:H solar cells with soft-imprinted plasmonic back reflectors,” Appl. Phys. Lett. 95(18), 183503 (2009). 6. R. A. Pala, J. White, E. Barnard, J. Liu, and M. L. Brongersma, “Design of plasmonic thin-film solar cells with broadband absorption enhancements,” Adv. Mater. (Deerfield Beach Fla.) 21(34), 3504–3509 (2009). 7. 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(S2 Suppl 2), A237–A245 (2010). 8. L. Dal Negro and S. V. Boriskina. “Deterministic aperiodic nanostructures for photonics and plasmonics applications.” Laser & Photonics Rev. 1–41 (2011). #178295 $15.00 USD Received 22 Oct 2012; revised 11 Dec 2012; accepted 5 Jan 2013; published 29 Jan 2013 (C) 2013 OSA 11 February 2013 / Vol. 21, No. 3 / OPTICS EXPRESS 2757 9. Y. A. Akimov and W. S. Koh, “Design of plasmonic nanoparticles for efficient subwavelength light trapping in thin-film solar cells,” Photonics 6, 155–161 (2011). 10. V. E. Ferry, M. A. Verschuuren, M. C. Lare, R. E. I. Schropp, H. A. Atwater, and A. Polman, “Optimized spatial correlations for broadband light trapping nanopatterns in high efficiency ultrathin film a-Si:H solar cells,” Nano Lett. 11(10), 4239–4245 (2011). 11. A. E. Ostfeld and D. Pacifici, “Plasmonic concentrators for enhanced light absorption in ultra-thin film organic photovoltaics,” Appl. Phys. Lett. 98(11), 113112 (2011). 12. S. Mei, T. Jie, L. Zhi-Yuan, C. Bing-Ying, Z. Dao-Zhong, J. Ai-Zi, and Y. Hai-Fang, “The role of periodicity in enhanced transmission through subwavelength hole arrays,” Chin. Phys. Lett. 23(2), 486–488 (2006). 13. F. Przybilla, C. Genet, and T. W. Ebbesen, “Enhanced transmission through Penrose subwavelength hole arrays,” Appl. Phys. Lett. 89(12), 121115 (2006). 14. T. Matsui, A. Agrawal, A. Nahata, and Z. V. Vardeny, “Transmission resonances through aperiodic arrays of subwavelength apertures,” Nature 446(7135), 517–521 (2007). 15. R. Penrose, “The role of aesthetics in pure and applied mathematical research,” Bull. Inst. Math. Appl. 10, 266– 271 (1974). 16. N. G. de Bruijn, “Algebraic theory of Penrose's non-periodic tilings of the plane, Pt. I & II,” Kon. Nederl. Akad. Wetensch. Proc. Ser. A 84, 39–66 (1981). 17. H. Raether, Surface Plasmons on Smooth and Rough Surfaces and on Gratings (Springer-Verlag, 1988). 18. J. A. Dionne, E. Verhagen, A. Polman, and H. A. Atwater, “Are negative index materials achievable with surface plasmon waveguides? A case study of three plasmonic geometries,” Opt. Express 16(23), 19001–19017 (2008). 19. D. Pacifici, H. J. Lezec, L. A. Sweatlock, R. J. Walters, and H. A. Atwater, “Universal optical transmission features in periodic and quasiperiodic hole arrays,” Opt. Express 16(12), 9222–9238 (2008). 20. D. Pacifici, H. J. Lezec, and H. A. Atwater, “All-optical modulation by plasmonic excitation of CdSe quantum dots,” Nat. Photonics 1(7), 402–406 (2007). 21. G. Dennler, M. C. Scharber, and C. J. Brabec, “Polymer-fullerene bulk-heterojunction solar cells,” Adv. Mater. (Deerfield Beach Fla.) 21(13), 1323–1338 (2009). 22. M. Reyes-Reyes, K. Kim, and D. L. Carroll, “High-efficiency photovoltaic devices based on annealed poly(3hexylthiophene) and 1-(3-methoxycarbonyl)-propyl-1-phenyl-(6,6)C61 blends,” Appl. Phys. Lett. 87(8), 083506 (2005). 23. W. Ma, C. Yang, X. Gong, K. Lee, and A. J. Heeger, “Thermally stable, efficient polymer solar cells with nanoscale control of the interpenetrating network morphology,” Adv. Funct. Mater. 15(10), 1617–1622 (2005). 24. M. A. Green, K. Emery, Y. Hishikawa, W. Warta, and E. D. Dunlop, “Solar cell efficiency tables (version 39),” Prog. Photovolt. Res. Appl. 20(1), 12–20 (2012). 25. L. Lüer, H.-J. Egelhaaf, D. Oelkrug, G. Cerullo, G. Lanzani, B.-H. Huisman, and D. de Leeuw, “Oxygeninduced quenching of photoexcited states in polythiophene films,” Org. Electron. 5(1-3), 83–89 (2004). 26. P. E. Shaw, A. Ruseckas, and I. D. W. Samuel, “Exciton diffusion measurements in poly(3-hexylthiophene),” Adv. Mater. (Deerfield Beach Fla.) 20(18), 3516–3520 (2008). 27. A. J. Moule and K. Meerholz, “Interference method for the determination of the complex refractive index of thin polymer layers,” Appl. Phys. Lett. 91(6), 061901 (2007). 28. E. Yablonovitch and G. D. Cody, “Intensity enhancement in textured optical sheets for solar cells,” IEEE Trans. Electron. Dev. 29(2), 300–305 (1982). 29. E. Yablonovitch, “Statistical ray optics,” J. Opt. Soc. Am. 72(7), 899–907 (1982). 30. T. Tiedje, E. Yablonovitch, G. D. Cody, and B. G. Brooks, “Limiting efficiency of silicon solar cells,” IEEE Trans. Electron. Dev. 31(5), 711–716 (1984). 31. J. Feng, V. S. Siu, A. Roelke, V. Mehta, S. Y. Rhieu, G. T. R. Palmore, and D. Pacifici, “Nanoscale plasmonic interferometers for multispectral, high-throughput biochemical sensing,” Nano Lett. 12(2), 602–609 (2012). 32. J. Gilot, I. Barbu, M. M. Wienk, and R. A. J. Janssen, “The use of ZnO as optical spacer in polymer solar cells: theoretical and experimental study,” Appl. Phys. Lett. 91(11), 113520 (2007).


Introduction
Thin film solar cells offer several advantages over traditional thick film devices-such as better electron transport properties and higher open circuit voltages [1]-but thicker films have the inherent advantage of being better absorbers of sunlight.To compensate for this, recently there has been a push to incorporate plasmonic concentrators (PCs) into thin film solar cells to improve the absorption properties of such devices, and thus improve the overall power conversion efficiency [2,3].There are numerous examples of this unification of the fields of plasmonics and photovoltaics with researchers using various arrays (e.g., periodic, random, or spiral) and nanoscale scatterers (e.g., disks, spheres, ridges, grooves, or holes) [4][5][6][7][8][9][10] as a means of improving optical absorption.The authors have recently demonstrated a particularly successful example of PCs arranged in two-dimensional quasiperiodic (QP) layouts [11].These deterministic aperiodic patterns were originally developed to tessellate planes with regular polygons, and can be understood as intermediaries between random and periodic: they lack long-range translational symmetry, but can have very high levels of local and long-range rotational symmetry, and can pack more holes into smaller areas.The points of the array correspond to the locations of the nanoscale light scatterers.
This report introduces a novel algorithm for generating QP patterns across a wide spectrum of point densities and orders of rotational symmetry (both local and long-range), and can also be used to reproduce certain periodic arrays.The results of this algorithm can be employed with PCs comprised of any type of symmetric scattering object (e.g., holes, spheres, or disks) to produce spectrally broad and angle-and polarization-insensitive enhancement.The scatterers used in this report are shallow holes with sub-wavelength diameters etched on the back metal contact, underneath the active layer.This particular type of PC is also known as a nanohole array (NHA).Nanoholes have two advantages in this context: in addition to having a relatively simple fabrication process -typically either via Nanoimprint Lithography (NIL), or Focused Ion Beam (FIB) milling, neither of which requires depositing extra material-nanoholes will not reflect or transmit a significant fraction of the incident light.It should also be noted that there have been studies on aperiodic NHAs for extraordinary transmission [12][13][14].However, those structures employ through-holes, which are not compatible with solar cells -the nanoholes of this paper are very shallow compared to the thickness of the metal substrate to prevent light from being transmitted.Further, those papers only consider a small handful of aperiodic patterns, and in some cases the patterns were not truly QP, in that they maintained translational symmetry through the use of repeating unit cells; this report will cover a wide variety of purely QP geometries.
When employed in a solar cell, an NHA works as follows: the holes scatter the incident light field and reradiate a fraction of it along the dielectric / metal interfacial plane in the form of surface plasmon polaritons (SPPs).This structure enhances absorption in the active layer through two methods: (1) having the light redirected along the interfacial plane allows for increased light-matter interaction compared to single-pass absorption through the thin film; (2) if two or more holes are close to each other, the waves can constructively interfere and create strong EM fields, both locally and averaged over the entire array area.One of the goals of this report is to fine-tune the QP generation algorithm to see how the array's properties (order of base polygon and inter-hole distance) affect absorption enhancement.Additionally, it will be argued that QP arrays are superior to comparable periodic and random arrays in absorption enhancement.It is believed that this concept should enhance absorption in any ultra-thin layer of absorbing material; an organic bulk heterojunction was chosen here as a proof of concept.
This report is broken into three sections.Section 1 will introduce the QP generation algorithm, provide a protocol on how to implement it, and explain some of the key properties as functions of its two tunable parameters.Section 2 delves into the underlying physics of plasmonic concentrators and thin film solar cells.In particular, the SPP propagation constants in a three layer system will be calculated, and are used to simulate the interference properties of several SPPs generated by an array of many nanoholes serving as localized sources of SPPs at various frequencies.In section 3, EM field enhancement and absorption enhancement are studied as functions of the two QP parameters using both simulations and experimental data.

Generalized construction algorithm for periodic and quasiperiodic arrays
To understand how quasiperiodic arrays were originally devised, imagine trying to tile a 2-D plane with regular N-sided polygons such that there are no empty spaces.This can be done easily with triangles (N = 3), squares (N = 4), and hexagons (N = 6), but there is a peculiar gap at N = 5.In 1974, Penrose developed a QP tiling method that employs two types of rhombi to tessellate the plane while maintaining five-fold rotational symmetry [15].In 1981, de Bruijn proposed an algebraic theory for the generation of Penrose's non-periodic tiling of the plane.His algorithm, known as the "cut and projection" method, involves taking a 5-D hyperspace, dividing it into unit cubes whose vertices are at integer coordinates, intersecting the hyperspace with a 2-D plane, then projecting the centers of the bisected unit cubes orthogonally onto the plane.De Bruijn demonstrated that the ensemble of the projected points can be turned by a similarity transformation into the full set of vertices of Penrose's tessellation [16].
This report introduces a generalization of the cut and projection method that can be used with a hypercube of any dimension (N).What follows is a step-by-step description of the algorithm for creating such patterns: 1. First, select a dimension number for the hypercube, N ≥ 3. The axes of the hypercube are also called grids.A system made with five grids is called "penta-grid" (also called "Penrose-grid"), one made with seven grids is called "hepta-grid," and so on.The grids exist in a complex space ("grid space") defined by z = u + iv, where u and v are real numbers.If N = 3, 4 or 6, the resulting array will be periodic; for all other values, the array is quasiperiodic.
2. Define a set of N real irrational numbers {γ j }, where j is the integer grid index (with 1 ≤ j ≤ N). γ j is the point on the jth grid where it is intersected by the 2-D plane.The user is free to pick any values for the set {γ j }, but in addition to being real and irrational, there are other conditions: the set must sum to zero, and must be selected such that every point in the complex plane belongs to no more than two of the N grids (in practice, this can be achieved by making each γ j value within a particular set unique).Different sets of {γ j } can produce slightly different real space arrays, but the reciprocal spaces will be identical.
3. The 2-D real space that contains the points of the tiling is called "array space," with real coordinates x and y.The function ( ) → maps the location of a vertex (i.e., a point of the QP array) from the N-hypercube onto a 2-D plane.It is defined by , and q     is the "roof" of q, defined as the smallest integer n such that n ≥ q.
4. Set a scaling factor a (this can be understood as an equivalent lattice constant) to define the characteristic hole-to-hole separation.For simple periodic arrays (like square and hexagonal), this is identical to the array "pitch" or "lattice constant" (the distance between nearest neighbors).Depending on the context, a can be expressed in terms of pixels or meters; for this report, a will be on the order of ~10 2 nm.
5. For a given value of z, the corresponding point (i.e., the location of the scatterer) in the real space array is located at ( ) 6.The process is repeated by rasterizing u and v in order to pick up every possible point in the array space (for the best results, the raster should be repeated multiple times, with the step size decreasing each time until a convergence in the number of points is observed).Typically, there are regions of neighboring points in grid space that are mapped to the same point in array space, so some array points will be found multiple times.These regions are called domains.The resultant arrays can be identified by just two tunable parameters: the scaling factor a, and the grid number N.Although the set {γ j } can be altered (within the rules laid out in step #2), that is not considered a tunable parameter because changing {γ j } usually does not produce significant differences in key properties like planar density of points or in the reciprocal spaces.
Figure 1 shows the major steps of this process for several representative periodic and quasiperiodic arrays.For each pattern, the left column shows its grid space and the middle column shows its array space -each domain-point pair (as determined by the mapping function f(z)) is drawn in the same color in both columns.The colors are chosen based on the distance of the hole from the central point of symmetry in each array space; that is, every hole with the same color is equidistant from a central point (18 colors in all).The array of points is then transformed to reciprocal space (right column -note that this is not on the same color scheme as the other two columns) using a discrete fast Fourier transform (DFFT).This is often necessary when studying patterns with high N values, as some of the relevant parameters (such as degree of rotational symmetry) are difficult to infer using the naked eye.For a given scaling factor, the area of the first Brillouin zone (FBZ) of every QP pattern is identical to the area of the FBZ of the corresponding tri-grid pattern.Since that quantity is well defined (a hexagon of area One of the most important properties is the order of rotational symmetry, which is found by counting the number of points that circumscribe the FBZ in the reciprocal space images.For example, the penta-grid example in Fig. 1 (3rd row, right column) displays 10-fold symmetry about the center.The order of local rotational symmetry is closely related to the number of nearest neighbors that a typical hole will have.For instance, the square array in Fig. 1 (second row) is defined by four-fold rotational symmetry, which means that every hole in a square array has four nearest neighbors (ignoring holes at the boundary).For the more complicated QP patterns, it is not accurate to say that each point will be surrounded by the same number of nearest neighbors, but rotational symmetry is still important because it indicates that the points in the array are spread out evenly.As will be discussed below, it is essential in this device that each hole and the available space are being used as efficiently as possible in the light concentration process.For example, a highly asymmetric pattern may have regions with few or no holes, thus diminishing the beneficial effects of the plasmonic concentrator -quantifying the pattern's rotational symmetry in such a manner helps the design process.
A quantitative overview of the predominant trends in these QP arrays is provided in Fig. 2. Figure 2(a) indicates that when N is even, the rotational symmetry of the array-as determined by analyzing its Fourier power spectrum -equals N, but when N is odd, the rotational symmetry of the array equals 2N. Figure 2(b) plots the density of holes as a function of grid number for certain a = 400 nm patterns.All QP arrays (open diamonds) have higher hole densities than the periodic arrays (full diamonds), which further supports the hypothesis that QP arrays are better than periodic ones; indeed, if rotational symmetry was of no concern, QP arrays would still be better than periodic arrays due to this fact.Finally, note that the density of holes increases as the scaling factor decreases (see Fig. 2(c), which reports data for a heptadeca-grid pattern, with N = 17 and varying a; the insets show the array and reciprocal space images for an N = 17, a = 250 nm QP pattern).As will be shown below, it is possible to create patterns that will produce much stronger absorption enhancement than a basic periodic array using only the two tunable parameters.For example, the data from Fig. 2 allows for a direct comparison of a = 400 nm heptadeca-grid and tri-grid patterns: although their FBZs are identical in size, the former has an order of rotational symmetry 5.67x greater than the latter (34 vs. 6), and a spatial density of points that is 1.66x greater (7.91 holes / μm 2 vs. 4.76 holes / μm 2 ).

The physics of plasmonic concentrators & thin film solar cells
Consider a semi-infinite dielectric / metal interface; when light (free space wavelength λ) is normally incident on this interface, most of the radiation will be reflected back in the direction it came (Fig. 3(a)).Now imagine that there is a shallow circular hole (diameter D λ << ; this is known as the sub-wavelength regime) etched on the metallic surface (Fig. 3(b)).When light hits the hole, the incident radiation is coupled to the electron plasma in the metal to form a surface plasmon polariton (SPP).SPP modes can only form if some sort of corrugation or prism is present on the interface to compensate for the mismatch in the incident wave's momentum and the SPP's momentum.Since an SPP is an electromagnetic wave, Maxwell's equations can be used to describe its propagation behavior.For the case of a semi-infinite dielectric / metal system, the parallel component of the wavevector is given by ( ) while the perpendicular component on the dielectric (i = 1) or metal (i = 2) side is ( ) ( ) ε ε ω = are the dimensionless frequency-dependent complex dielectric functions of the dielectric and metal materials, respectively [17].Other useful SPP parameters include the wavelength / Re[ ] , and skin depth , , . In the case of a three-layer system (vacuum / dielectric / metal), the dispersion relation will be calculated following the calculations developed by Dionne, et al. [18].When generalizing to a 2-D interface, the in-plane SPP wavevector will be referred to as SPP k .Although Fig. 3(b) shows only two holes, the figure does help demonstrate how the nanohole array (NHA) plasmonic concentrator works when scaled up to hundreds or thousands of holes, in which case the two tunable parameters of the QP array generator can profoundly affect the absorption enhancement process.Recall that the scaling factor (a) is directly related to the density of holes (Fig. 2(c)).Generally speaking, a smaller value for a is better, as it means there are more locations where the incident light has a chance to scatter into SPP modes, and having the holes closer together means that the SPP waves do not have to travel as far to meet a counter-propagating wave and constructively interfere with it.The grid number (N) is related to the rotational symmetry of the array, and thus how many nearest neighbors a typical hole will have.As the incident sunlight is randomly polarized, it is expected that, on average, the SPP modes will emerge from the holes omnidirectionally, so having as many nearest neighbors as possible should lead to better constructive interference patterns.
It has been shown experimentally that SPP modes do propagate on the interfacial plane [19] and can exhibit bulk-like absorption within an optically-thin superstrate dielectric material (e.g., by a 20nm-thick layer of CdSe quantum dots [20]).The absorbing dielectric material here will be the organic bulk heterojunction formed by the polymer poly(3hexylthiophene) (P3HT) and the fullerene (6,6)-phenyl C 61 butyric acid methyl ester (PCBM); the metal will be silver.P3HT:PCBM is one of the most widely studied active layers in organic thin film photovoltaics [21]; its highest recorded power conversion efficiency is 5% [22,23].By comparison, the current state of the art for all single-junction organic solar cells is 10% [24].One of the main setbacks hindering effective photovoltaic conversion in organic materials is the low exciton diffusion length (e.g., 4-12 nm in P3HT:PCBM [25,26]).Making the film thickness comparable to the diffusion length provides a way of avoiding this issue; incorporating a PC (such as an NHA) into such a device can enhance absorption in an otherwise optically thin layer.In this manner, the benefits of a thin film solar cell (improved carrier transport behavior) are merged with the benefits of a thick film solar cell (improved light absorption) to create a high-efficiency device.defines the scale of the extra light-matter interaction that occurs when the incident light is redirected perpendicular to the direction of incidence.Second, note that local constructive interference cannot occur if the SPP modes die out shortly after leaving each hole.Generally speaking, the device performs better when the SPP propagation length is as large as possible (this is why silver was chosen as the metallic substrate in this paper -its dielectric properties are such that Im[n SPP ] is relatively low).The inset in Fig. 4(b) shows normalized plots of the E x component of the SPP wave parameterized by d, with the incident light at λ = 700 nm.This shows the transition of the SPP mode from vacuum, to ultra-thin, to thin, and to bulk.The fact that the field amplitude rapidly decays on the metal side of the interface (quantified by the skin depth, δ z,2 ) further demonstrates that Ag is an appropriate choice for the substrate material -with most of the SPP's energy on the dielectric side, there is a better chance for absorption and photocurrent enhancement, whereas energy on the metal side is essentially unusable.Although the propagation length increases significantly at higher wavelengths, this should be tempered by the fact that the band gap of P3HT is 670 nm ( = 1.85 eV), so absorption beyond that wavelength will probably not produce any significant amount of photocurrent.
It appears as though the system reaches the "bulk limit" before d = 150 nm, as the data for d = 150 nm and d = ∞ look quite similar.Aside from the benefits inherent in using a 24-nmthick layer vis-à-vis the low exciton diffusion length, there are several reasons why a plasmonic concentrator etched on the back contact will have a greater impact when the active layer is thinner: i.The propagation length increases as d decreases, as seen in Fig. 4(b).This means that the plasmonic concentrator (PC) will be more effective if coated with a thinner absorbing layer.
ii.Having a thinner layer means that more of the incident light will reach the NHA in the first place; therefore, an increased fraction of the incident light is scattered and turned into usable SPP modes.
iii.SPP modes are confined to the dielectric / metal interface and have finite skin depths; thus, most of the SPP's energy will not reach the upper parts of thicker layers.
All this information, combined with previous research on SPP generation and absorption, suggests that the NHA paradigm should work to help improve the efficiency of a solar cell.Finally, note that the absorption properties of the devices discussed in this paper cannot be analyzed in terms of the well-known "ergodic limit," which claims that the maximum absorption enhancement factor a thin film can exhibit is 4n 2 , where n is the index of refraction of the material [28][29][30].The ergodic limit applies to films with textured / inhomogeneous surfaces -in this model, light can be trapped in the thin film in a waveguide-like mode, but reflections off the corrugated surface will randomize the paths of the rays, thus increasing the effective path length by 4n 2 when averaged over all angles, and when the film is adjacent to a perfectly reflecting surface on its back side.The derivation assumes that bulk absorption is weak enough throughout the layer such that the internal light intensity is uniform (that is, αd << 1, where α = 4πκ / λ is the absorption coefficient of the material, κ is the wavelengthdependent imaginary part of the index of refraction, and d is the thickness of the layer), while at the same time assuming that the layer is much thicker than the incident wavelength λ (that is, d >> λ).The first assumption is needed to give the direction of the light rays a chance to be sufficiently randomized by the corrugated surface; the second assumption is required when employing classical ray tracing arguments.In ultra-thin films, constructive interference between reflected waves within the layer is key in understanding absorption phenomena; optical ray tracing does not account for that.The second assumption clearly does not apply in this case, where the layer thickness is typically d = 24 nm, and the incident wavelength is on the order of ~10 2 nm.When trying to apply both assumptions simultaneously, it requires that 1 << d / λ << 1/(4πκ); since κ for most materials (including P3HT:PCBM) is large when λ is small, this can present a contradiction, or will only allow the approximation to work in a narrow wavelength range.

Description of simulation program
One of the hypotheses of this section is that NHAs made with higher N values will provide better absorption enhancement than ones with low N values, as high N patterns have higher orders of rotational symmetry (Fig. 2(a)) and higher densities of holes (Fig. 2(b)).The next set of simulations model the interference behavior of SPP modes in two different NHAs with the same scaling factor: one periodic (low N), and one quasiperiodic (high N).In theory, the latter should have stronger EM fields across the entire plane, which implies stronger absorption in the active layer.
Although there are commercially available programs to model such a system (e.g., Lumerical's Finite Difference Time Domain Maxwell Solver), they tend to be slow and computationally intensive.The authors present a quicker and easier simulation code, where an inherently 3-D problem is reduced to a simpler 2-D problem by confining the problem to the P3HT:PCBM / Ag interface, then using knowledge of the vertical field profiles (see the inset of Fig. 4(b) for examples) and the skin depths to complete the picture.In this model, each hole is assumed to be a source of planar SPP modes described by the three-layer k SPP wavevector (which is understood to have been created by a monochromatic, normally incident wave of wavelength λ).It is also assumed that the incident light is randomly polarized, so the SPP modes will propagate omnidirectionally and radially outward from each hole.To express this quantitatively, the field at any point on the plane (located at r) due to the jth hole (located at r j ) is given by ( ) ( ) where H 0 is the amplitude of the incident wave, β j is the scattering efficiency of the jth hole (defined as the fraction of incident light that is turned into an SPP, originally introduced in Fig. 3(b)), and φ j is the relative phase between the incident wave and the SPP wave emanating from the jth hole.The program renormalizes at the location of the hole itself (r = r j ) to avoid singularities.Each field spreads out evenly in a circle on the planar interface, which accounts for the inverse square root dependence.Since SPPs are transverse magnetic (TM) and confined to the plane, this is the only component of the magnetic field present [17].The fields due to every hole are summed together at every point in space to produce a single map (specific for each incident wavelength).A map can be identified by the function where the total SPP field intensity is normalized against the incident field to produce a quantity in arbitrary units.To simplify things further, any map can be described with a single number using the spatial average integral ( ) , also in arbitrary units.It is important to note that ( ) , I λ r only includes the field intensity of the interfacial SPP modes; it does not include any information about EM fields present in the absorbing layer due to classical light reflection and transmission behavior.In Eq. ( 3), the two parameters β j and φ j can be altered as needed to account for different devices and conditions.In general, β j and φ j can depend on the materials used, the incident wavelength, the angle of incidence, the size and shape of the hole, etc.In the current model, all the holes are identical in size and shape and are being illuminated by the same light, so it is not unreasonable to assume that β j and φ j are constants here.To obtain usable values for these parameters, data from an extraordinary transmission experiment was used, wherein light (λ = 514.5 nm) was normally incident on two holes (D = 50 nm) bored through a layer of Ag, with SiO 2 on top [19].It was found that the phase difference between the incident wave and the SPP wave was π / 2, so φ j is set accordingly here.The SPP wave amplitude was found to be 35% of the incident amplitude, so β j is set to 0.35 μm 1/2 in Eq. ( 3).There is reason to believe that the holes scatter SPP modes reasonably well at all of the relevant incident wavelengthsthat is, β j may not be exactly 35% in all cases, but it should not deviate by a large amount, either.A systematic study is still needed, but it is known from other experiments that the scattering coefficient of a 200nm-wide shallow groove at a vacuum / Ag interfaces varies between 17% and 27% for incident light in the visible range [31].I λ r (the highest intensity value observed) for each map.In the low wavelength limit, it is expected that the relative advantages and disadvantages of the different patterns should not matter as much because the SPP propagation length is so low (see Fig. 4(b)), so the SPP's energy will quickly dissipate and almost no constructive interference will occur.This is what is observed, but there still may be certain regions were two or more holes happen to be extremely close together by chance (this is more likely to happen in the heptadeca-grid case, due to its inherent non-periodicity), which will produce a large value for

Absorption enhancement as a function of grid number
even though the average intensity ( ) I λ will be roughly the same in the two patterns (see Table 1).As λ and SPP Λ increase, constructive interference becomes more prominent, so it is necessary to closely examine the mechanisms underlying this phenomenon.It was stated before that the heptadeca-grid pattern should produce higher SPP fields on average due to its higher density of holes (which implies more scattering centers and more chances for constructive interference), and because of its higher order of rotational symmetry.To better understand the importance of rotational symmetry, consider the fact that any given hole in the square array only has four directions along which other holes are spaced at distances that are integer multiples of a.The existence of these favored angles is confirmed by looking at the reciprocal space image of the square pattern (or indeed any periodic pattern), which shows discrete points as opposed to diffuse rings.Along every other direction, the hole distance are not integer multiples of a, which leads to the possibility of destructive interference.Indeed, that is what is observed in the top row of Fig. 5: as the wavelength increases, out of phase SPPs emanating from holes not along the preferred line of sight will reach this "test" hole.
Ergo, the amount destructive interference increases, which in turn causes ( ) I λ to decrease (see Table 1).Looking at the top row more closely (in particular, λ = 580 nm), it is clear that the local maxima are stronger near the edges of the array than in the center -this could be explained by the fact that holes on the edge of the pattern are less likely to have an off-angle neighbor, so there is less destructive interference.On the other hand, there is no favored angle in a quasiperiodic array, as the diffuse rings surrounding the first Brillouin zone can attest to.The diffusiveness of this ring (or, equivalently, the average number of nearest neighbors a typical hole will have) is a function of N. In practice, this means a typical hole in a QP pattern is more likely to have a nearest neighbor with the correct distance for constructive interference to occur, which is why increasing the wavelength (and increasing SPP Λ ) has a net positive effect, as seen in Table 1.This model can offer significant insight into SPP behavior in thin films, but it must first be tested against experimental data.Ostfeld and Pacifici studied absorption enhancement in layers of P3HT:PCBM resting on Ag films that had been etched with periodic and quasiperiodic NHAs [11].Although the data from that report clearly demonstrated that having an NHA will enhance absorption in organic thin films across a wide wavelength range, it was not well understood why some patterns worked better than others.To prepare the samples, P3HT and PCBM (1:1 weight ratio) were dissolved in dichlorobenzene at 50°C, then spincoated (at 5000 rpm for 60 sec to create a 24-nm-thick P3HT:PCBM layer, or at 1000 rpm for 60 sec to create a 150-nm-thick layer) onto 300-nm-thick Ag films that had been etched (via Focused Ion Beam milling) with periodic or QP NHAs, all of which was followed by annealing at 110°C for 10 minutes.The nanoholes were 70 nm deep, which is shallow enough to assume that no light is transmitted though the Ag layer.Each NHA had a total area of 100 μm × 100 μm, and the hole diameter was D = 80 nm in all cases.The incident beam was sent through a 40x objective lens with a numerical aperture of NA = 0.6, then collected by the same lens to record the combined specular and diffuse reflectance (R).The spectral absorptance (A) was obtained using the formula A = 1 -R (since transmittance T = 0).Note that this is in contrast to the simulations above, where the incident light is assumed to be normal.[11].The d = 0 case refers to just a vacuum / Ag interface.The equivalent lattice constant is a = 400 nm in all cases.(b) Spectral absorption enhancement when comparing a corrugated structure to an uncorrugated structure.The curves shown here are found by taking the "corrugated" curves from part (a) and dividing by an "uncorrugated" curve of the same thickness.(c) Comparison of heptadeca-grid to square in two different contexts (both at d = 24 nm).Left axis (solid magenta line): ratio of experimental absorptance (from panel (a)).Right axis (orange spheres): ratio of simulated normalized spatially-averaged field intensities (from simulations similar to those reported in Fig. 5).
Figure 6(a) plots the spectral absorptance in the P3HT:PCBM layer in a variety of devices: different active layer thicknesses were used, and back metal surfaces could have square, heptadeca-grid, or no NHA present (if an NHA is present, the scaling factor is a = 400 nm).The case plotted with the thin dotted line (d = 0, no NHA) is identical to a piece of pure, uncorrugated silver.This figure shows that the gap absorptance between a thin 150 nm layer and an ultra-thin 24 nm layer can be partially filled by incorporating an NHA into the device.This means that a 24-nm-thick P3HT:PCBM solar cell can benefit from the high absorption of a bulk-like layer without sacrificing its higher internal quantum efficiency (IQE).Of the two NHAs considered from this graph, the heptadeca-grid NHA outperforms the square NHA across the entire recorded spectrum, further confirming the original hypothesis.In Fig. 6(b), the "corrugated" (with NHA) curves of Fig. 6(a) are divided by the corresponding "uncorrugated" (without NHA) curves of the same thickness.The data shows that the NHA has very little effect in the d = 150 nm case, as the enhancement level is close to one for the entire recorded spectrum.On the other hand, the ultra-thin d = 24 nm cases show very strong absorptance enhancement, peaking at a factor of 6 near the band gap for the N = 17 case.These results confirm the hypothesis laid out in the end of section 2, which argued that the NHA will have a stronger effect if placed under a thinner absorbing layer.
Figure 6(c) compares heptadeca-grid to square in terms of both simulated and experimental data (everything in Fig. 6(b) has a = 400 nm and d = 24 nm).The left axis (solid magenta line) plots the ratio the experimental absorptance curves (from panel (a)), while the right axis (orange spheres) plots the ratio of the ( ) I λ values (from Fig. 5 and Table 1).
Absorptance is proportional to the square of the electromagnetic field, which is exactly what ( ) I λ measures.Since the only difference between the corresponding N = 4 and N = 17 curves is the nature of the extra SPP field present on the P3HT:PCBM / Ag interface, it is expected that these quantities are directly proportional to each other.Indeed, that is what is seen here, as the data follow the same trend across the entire recorded spectrum.Notice that both curves are at ~1 (meaning no relative enhancement) up until about 580 nm, which is also where the SPP Λ curve starts to increase (see Fig. 4(b)).In the high SPP Λ regime, heptadecagrid should dominate due to its higher density of holes and higher order of rotational symmetry, which is observed here.The fact that the experiments and simulations match up so well gives credence to the model developed above; further improvements to the model will be discussed in future reports.The rest of the data from [11] is summarized in Table 2.The experimental spectral absorptance was used to calculate the short circuit current density (J sc ) via the formula where e is the charge of an electron, ( ) IQE η λ is the internal quantum efficiency, and is the standard AM1.5G solar intensity spectrum (in units of photons / (s m 2 nm)).As a conservative estimate, it will be assumed that no significant exciton generation occurs beyond the band edge; thus, the range of integration was 400 nm < λ < 670 nm.As for the IQE, it is expected to be very close to 100% for a layer as thin as this one; one recent report put it at ~95% for a similar P3HT:PCBM device [32]; for simplicity, it will be assumed that IQE is not a function of wavelength.The data shows that not only does the presence of an NHA enhance absorption in all cases (when compared to an uncorrugated device whose J sc is 5.204 mA / cm 2 ), but also NHAs with higher values of N work better than ones with low N values.

Absorption enhancement as a function of array scaling factor
Although it is expected that a higher density of holes will produce better absorption in the active layer, there is limit on how dense the holes can be before the surface of the metal is completely etched away (at which point the plasmonic concentrator would essentially be a nanopillar array).There needs to be enough material surrounding the hole for the electron plasma to reside.Considering the chosen hole diameter (D = 80 nm) and the known surface plasmon polariton (SPP) propagation lengths (Fig. 4(b)), it is expected that the ideal scaling factor should be on the order of 10 2 nm, which is why a = 400 nm was the primary value used in this report.The next data sets look at how the device performs when the grid number is held constant but the scaling factor changes. Figure 7 shows simulations of ( ) I λ r for each map is provided in Table 3.As with Fig. 5, the results can largely be understood in terms of the SPP propagation length: when the wavelength is low, the SPP propagation length is also low, so the overall enhancement is weak.However, the 200 nm case produces better enhancement from the sheer volume of holes -there are more scattering centers, thus more SPPs mode present on the interface.There is only a moderate amount of constructive interference occurring, and is only visible in a few specific locations were there happen to be two or more holes very close together by chance.By the same token, the value of ( )

max[ , ]
I λ r for the 850 nm case probably would have been even lower had it not been for one small region where two holes are very close together; ultimately, the value for ( ) I λ in the 850 nm case is extremely small in the low wavelength regime, as expected.It is interesting to note that the 850 nm case is so sparse that increasing SPP Λ by a factor of 2.5 (from ( 450 nm) The experiments from section 3.2 were repeated for this group of devices, but the maximum incident wavelength was reduced to 650 nm. Figure 8(a) shows the experimental absorption spectrum as a function of scaling factor.It can be seen that the absorptance maximum redshifts as a increases -this is probably a combination of varying interference conditions by increasing the characteristic hole-hole spacing, and of longer propagation lengths needed as the holes move farther apart, which occurs if the wavelength increases.Figure 8(b) takes the ratio of the absorptance curve against that of an uncorrugated device (recall that the a = 200 and 850 nm patterns were shown in Fig. 7; the 400 nm pattern was shown in Fig. 5, lower row).Throughout almost the entire reported spectrum, the absorptance in a corrugated device equals or exceeds the absorptance in an uncorrugated device, as evidenced by the fact that the three normalized curves in Fig. 8(b) are greater than or equal 1.Again, SPP Λ provides the insight needed to understand this figure: at low wavelengths, the enhancement in the 400 nm and 850 nm cases is negligible, as the holes are too far apart and the propagation length is too low, but the 200 nm case is dense enough to see some enhancement.As SPP Λ increases, the enhancement factor starts to improve, but much more slowly for the very sparse 850 nm case.
Using Eq. ( 5), a plot of short-circuit current density versus scaling factor is presented in Fig. 9 for a heptadeca-grid and a random NHA; Scanning Electron Microscope images of the a = 200 nm and a = 850 nm heptadeca-grid NHAs (fabricated on Ag) are also provided.At each data point on the a axis, the two curves have the same density of holes, but the former was made with the standard QP algorithm, and the latter was made by randomizing the location of the points (while still avoiding accidental overlap).In the available data, the QP pattern peaks at λ = 300 nm, while the random pattern peaks at 200 nm.In other words, the QP pattern does not work best when the hole density is as high as possible, but the random one does.The latter point can be more easily explained: because there is no symmetry to speak of, hole density is the only relevant parameter.There appears to be something akin to a ~a-2 relationship between J sc and a for the random curve.This is not too surprising, as the density of holes is closely related to a, and the density of holes will have a strong effect on J sc .In fact, the density of holes in a square array is precisely defined by this a −2 expression, but more complicated patterns may have a slightly different relationship; this a −2 trend was also present in Fig. 2(c).The other key point in Fig. 9-that the smallest studied value for a does not produce the best enhancement-cannot be explained as easily.One phenomenon that has not yet been considered is the possibility of an SPP being scattered back up into the vacuum when it hits another hole.Although the odds of this happening are generally small (and only a fraction of the energy will be lost in such a process), it is more likely to occur in the very dense 200 nm case.This could help explain why a "compromise" value of a = 300 nm produces the best J sc for the heptadeca-grid case.Fig. 9. Calculated short-circuit current (J sc ) as a function of NHA scaling factor, as well as the percent increase with respect to an uncorrugated device (J sc = 4.09 mA / cm 2 ), as calculated using the data in Fig. 8 and assuming AM1.5G illumination conditions.In all cases, d = 24 nm, and IQE = 95%.The two curves refer to patterns with the same planar density of holes, but one is QP (N = 17) and the other is purely random.Insets: Scanning Electron Microscope images of heptadeca-grid NHAs (one at a = 200 nm, the other at a = 850 nm) etched on a 300-nmthick Ag film.
Looking at the general trend, Fig. 9 shows that the calculated J sc in the presence of either a random or heptadeca-grid NHA will produce a higher short-circuit current density when compared to the same solar cell without a PC present on the back metal contact, no matter what the scaling factor is.It is interesting to note that in both the case of extremely small a ( = 200 nm) and the extremely large a ( = 850 nm), the two patterns produce roughly the same J sc .The latter situation is easily explained, as the holes are so far apart that next to no constructive interference will occur at all, and any enhancement is purely due to increasing the scale of light-matter interaction.In the former case, it is possible that the density of points is so high that every hole will have a relatively close by neighbor along every line of sight, even in the random case.Thus, the otherwise beneficial effects of high rotational symmetry will be neutralized in this ultra-high density regime.

Conclusion
This report has introduced a generalized tiling algorithm that can generate 2-D quasiperiodic patterns; both the planar density of points and the orders of local and long-range rotational symmetry can be fine-tuned across a wide range of values.While this may prove to be interesting and insightful from a purely mathematical point of view, this algorithm was developed specially for use with in nanoscale plasmonic concentrators, in particular with plasmonic concentrators that have been incorporated into thin film solar cells.This merging of the fields of plasmonics and photovoltaics is needed to compensate for the poor absorption properties of thin film photovoltaic materials.By corrugating the back metal contact with a quasiperiodic array of nanoscale holes, it was found that absorption enhancement is a function of the two tunable parameters used in the QP array generation algorithm (N and a).This can be understood by examining the underlying physics of PCs, where incident light is coupled to surface plasmons in the metallic substrate then propagate along the interfacial plane, thus increasing the scale of light-matter interaction and increasing the chances of constructive interference.The layout of the nanoholes was found to be vital in determining the level of constructive interference-the level of enhancement was greater when the hole-hole spacing in the pattern is low, and when the order of local rotational symmetry was high.There are still numerous questions that need to be resolved in future research, such as how high the grid number can go before the enhancement level saturates, and if different arrays with the same symmetry behave similarly.Additionally, the simulations can be improved by further exploring the scattering mechanism of the nanoholes-specifically, extracting more information on the scattering efficiency (β) and launching phase (φ) of a nanohole as functions of intrinsic and extrinsic parameters.In any case, this novel paradigm should significantly expand the fields of plamonics and thin film photovoltaics, as the generation algorithm presented here can be employed not only for hole arrays, but also with any other planar arrangement of sub-wavelength scatterers (such as disks, nanoparticles, and more).

Fig. 1 .
Fig. 1.Visual representation of the generalized cut and projection algorithm, starting in complex "grid space", then moving to the real "array space" (which contains the actual pattern), then ending in reciprocal space (for analytical purposes).The algorithm works by mapping -via the function f(z) -domains in the complex plane to points in real space, thus determining the position of points in the array.Corresponding domain-hole pairs are identified by the same color (colors are based on the distance from the central point of symmetry; 18 shown here).The examples of planar tilings provided have the same scaling factor (a) and hole diameter (D), but are made with different grid numbers (N).The top two examples are periodic: tri-grid (N = 3, a.k.a.honeycomb) and quad-grid (N = 4, a.k.a.square).The bottom two examples are quasiperiodic: penta-grid (N = 5, a.k.a.Penrose) and octo-grid (N = 8).The figures in the right column were obtained by performing a 2-D discrete fast Fourier transform (DFFT) on the real space array; the scale bar in those figures is 2π / a.The insets in the right column show the first Brillouin zones of the Fourier power spectrum of each pattern; the scale bar in those inset figures is π / a.

Fig. 2 .
Fig. 2. (a) Order of rotational symmetry as a function of "grid number" N (equivalent to the dimension of the hypercube used in the "cut and projection" algorithm).If a pattern is made with an even number of grids, it has N-fold symmetry; if a pattern is made with an odd number of grids, if has 2N-fold symmetry.(b) Planar density of holes as a function of N, with the scaling factor held constant at a = 400 nm.(c) Planar density of holes as a function of a, with the grid number set at N = 17 (heptadeca-grid).
frequency, c is the speed of light, SPP n is the effective index of refraction of an SPP mode, and

Fig. 3 .
Fig. 3. Schematic detailing the process of absorption enhancement in a plasmonic concentrator (PC) realized using a nanohole array (NHA).(a) Light normally incident on a flat metal surface is mostly reflected back.(b) If light is normally incident on a nanocorrugated metal / dielectric interface, some fraction of the incident radiation can be converted to a propagating SPP mode (with amplitude scattering efficiency β).Constructive interference between SPP modes originating from neighboring holes can significantly increase the light intensity at the metal / dielectric interface, thus contributing to enhanced light absorption if the dielectric material is an optical absorber.

Fig. 4 .
Fig. 4. (a) Calculated SPP mode dispersion relation (incident wavelength vs. real part of the inplane SPP wavevector), shown for four different P3HT:PCBM thicknesses.The dielectric function for silver was determined via spectroscopic ellipsometry; the values for P3HT:PCBM came from Moule and Meerholz [27].(b) SPP propagation length as a function of incident wavelength, parameterized by d.Inset: component of the SPP electric field parallel to the metal / dielectric interface (E x ) for the same four P3HT:PCBM thicknesses, calculated at an incident wavelength of λ = 700 nm.The field profiles are to scale and normalized to the highest value in each curve.

Figure 4
Figure 4 plots the SPP dispersion relation (the real part of k x vs. incident wavelength) and SPP propagation length in four different configurations parameterized by d, the thickness of the middle P3HT:PCBM layer in a three-layer vacuum / P3HT:PCBM / Ag system.The dispersion relations have been numerically calculated using the experimentally determined complex dielectric constants of Ag and P3HT:PCBM [27].The d = 0 (no organic layer) and d = ∞ (no vacuum layer) cases were modeled as two-layer systems; the d = 24 nm and d = 150 nm cases were modeled as three-layer systems.The SPP propagation length graph provides significant insight in how an NHA can enhance absorption in the active layer.First, SPP Λ

Fig. 5 .
Fig. 5. Simulations of the normalized field intensity I(r,λ) in square (periodic, N = 4) and heptadeca-grid (QP, N = 17) NHAs at four different incident wavelengths.In all cases, the array scaling factor is a = 400 nm and the organic layer is 24 nm thick.The scale bars seen in the square, λ = 500 nm case (both in the main figure and the inset) apply to all other images.Note that the color bar has been rescaled from 0 to 10 for clarity; the color seen on the map does not necessarily correspond to the actual value of I(r,λ) (the zero value stays the same after the rescale).The calculated SPP propagation length is included for each of the four wavelengths studied (obtained from Fig. 4(b)).

Figure 5
Figure 5 presents a set of simulations of ( ) , I λ r at the P3HT:PCBM (d = 24 nm) / Ag interface for the square (N = 4) and heptadeca-grid (N = 17) NHAs at certain wavelengths, all with scaling factor a = 400 nm.The inset in each panel is a zoomed-in image from the same location.Table 1 lists the calculated values for

Fig. 6 .
Fig. 6.(a) Experimental spectral absorptance (A) in a P3HT:PCBM film for different thicknesses (d) and NHAs [11].The d = 0 case refers to just a vacuum / Ag interface.The equivalent lattice constant is a = 400 nm in all cases.(b) Spectral absorption enhancement when comparing a corrugated structure to an uncorrugated structure.The curves shown here are found by taking the "corrugated" curves from part (a) and dividing by an "uncorrugated" curve of the same thickness.(c) Comparison of heptadeca-grid to square in two different contexts (both at d = 24 nm).Left axis (solid magenta line): ratio of experimental absorptance (from panel (a)).Right axis (orange spheres): ratio of simulated normalized spatially-averaged field intensities (from simulations similar to those reported in Fig.5).

Fig. 7 .
Fig. 7. Simulations of the normalized field intensity (I(r,λ)) in PCs consisting of various NHAs with a = 200 or 850 nm, at three different incident wavelengths.In all cases, the grid number is N = 17 and the organic layer is d = 24 nm.The scale bars seen in the a = 200 nm, λ = 450 nm case (both in the main figure and the inset) apply to all other images.As with Fig. 5, the color bar has been rescaled from 0 to 10, with the zero point same after rescale).The calculated SPP propagation length is included for each of the three wavelengths studied (obtained from Fig. 4(b)).
and scaling factors, with N = 17 and d = 24 nm.A summary of ( ) λ .It is not until the SPP propagation length exceeds the characteristic distance of 850 nm does any sort of constructive interference start to occur.

Fig. 8 .
Fig. 8. (a) Experimental spectral absorptance of a 24-nm-thick P3HT:PCBM film in the presence of several PCs consisting of various quasiperiodic NHAs with the same order of rotational symmetry (heptadeca-grid, N = 17) but different equivalent lattice constants (a, in nm), as well as an uncorrugated structure (to serve as a reference for the absorption enhancement).(b) Absorptance enhancement in the specific NHAs, defined as that device's absorptance curve divided by the uncorrugated curve.

Table 1
lists the calculated values for

Table 1 . Summary of maps from Fig. 5 (a comparison of N = 17 vs. N = 4 when a = 400 nm): space-averaged SPP field intensity (<I(λ)>) and highest intensity that is observed on a given map
(max[I(r,λ)]).