# Super Unit Cells in Aperture-Based Metamaterials

Authors: Dragan Tanasković; Zoran Jakšić; Marko Obradov; Olga Jakšić
Journal: Journal of Nanomaterials
 (2015)

---

## Abstract

An important class of electromagnetic metamaterials are aperture-based metasurfaces. Examples include extraordinary optical transmission arrays and double fishnets with negative refractive index. We analyze a generalization of such metamaterials where a simple aperture is now replaced by a compound object formed by superposition of two or more primitive objects (e.g., rectangles, circles, and ellipses). Thus obtained “super unit cell” shows far richer behavior than the subobjects that comprise it. We show that nonlocalities introduced by overlapping simple subobjects can be used to produce large deviations of spectral dispersion even for small additive modifications of the basic geometry. Technologically, some super cells may be fabricated by simple spatial shifting of the existing photolithographic masks. In our investigation we applied analytical calculations andab initio finite element modeling to prove the possibility to tailor the dispersion including resonances for plasmonic nanocomposites by adjusting the local geometry and exploiting localized interactions at a subwavelength level. Any desired form could be defined using simple primitive objects, making the situation a geometrical analog of the case of series expansion of a function. Thus an additional degree of tunability of metamaterials is obtained. The obtained designer structures can be applied in different fields like waveguiding and sensing.


---

## Body

## 1. Introduction

Nanocomposites based on two-dimensional arrays of subwavelength apertures in conductive material represent probably the most important class of metamaterials for the optical wavelength range. An example of such artificial structures is the fishnet metamaterials, which are the structures of choice for achieving negative value of effective refractive index in the optical range [1–3]. Another example is the extraordinary transmission arrays [4–6], where subwavelength apertures ensure, contraintuitively, an almost 100% optical transmission in a wavelength range determined by the surface plasmon resonance effects.The function of these aperture-based metamaterials is based on the electromagnetic interplay between metal and dielectric parts in their composite periodic structures, where the dimensions of the features and the sizes of their unit cells are subwavelength. The fact that the radiation in the operating wavelength range is much larger than the structural details ensures the possibility to describe them by the effective values of their optical parameters [7]. In the ranges where the effective value of relative magnetic permeability and the relative dielectric permittivity are simultaneously below zero, the effective refractive index, being equal to the square root of the permittivity-permeability product, will be below zero too [8].One can observe subwavelength metal-dielectric composites from a more general point of view, where it can be said that the possibility to tailor their geometry and composition results in engineering of frequency dispersion, where one is able to reach different desired values of effective refractive index in the targeted range. It is also possible to localize electromagnetic field at deep subwavelength level using plasmonic nanocomposites and to ensure extreme concentrations at very small volumes [9]. A typical example of this kind of applications is nanoantennas [10–12], where one modifies geometry and composition to obtain an almost limitless range of different modes. Actually a nanoaperture can be observed as a Babinet-type complementary analogue [13] to plasmonic nanoparticle-based antenna. Strong spatial nonlocalities affect the electromagnetic behavior of the plasmonic nanocomposites.Since different geometrical shapes result in different mode distributions and different subwavelength field concentration levels and spatial shapes, a question is posed: under which conditions one could apply effective medium approximation (EMA) [14, 15]? In the particular case of double fishnets, the question is whether it is possible to make relatively small modifications of the aperture shape in order to obtain large changes of the effective properties. One could expect that the effects of field concentration should be largest around the sharpest corners and thus one could tune and reshape the field distribution by introducing such forms.Changes of unit cell geometry within a fishnet structure have been reported in [16, 17]. An introduction of deeply subwavelength notches at the corners of a rectangular nanoaperture resulted in fine tuning of the effective optical properties [16].In this paper we consider the possibility to apply rectangular patch-like extensions to the fishnet apertures, thus redistributing the electromagnetic field. In other words, we utilize compound unit cells that represent a superposition of larger and smaller rectangles. Albeit the effective geometrical and material parameters of these structures (like their aperture area and metal-to-dielectric ratio) remain very similar to those of the unperturbed initial structures, the field distribution is sharply changed. We use finite element method to investigate the amount of change thus caused to the values of the effective refractive index and permittivity and to check if this offers a possibility to tailor the effective parameters in a wider range of values.


## 2. Method

### 2.1. Basic Concept

The basic element of a fishnet structure, whether single-metal layer (extraordinary optical transmission array) [18], double layer (double fishnet) [19], or multiple [20], is a single subwavelength aperture. This aperture may have different shapes, including circular, elliptical, and rectangular. One can define a unit cell for a fishnet by cutting four neighboring apertures through the middle and thus defining a cross-like metal structure where a quarter of aperture is located in each corner.In this work we modify the basic element by superposing additional apertures onto those already existing. We denote the basic aperture by the unperturbed structure, while the additional apertures are subobjects, substructures, or “patches.” Both the unperturbed structure and the added substructures are simple geometrical forms like circles, rectangles, and so forth, denoted by the primitive objects or primitives. Thus one obtains a compound unit cell (a “super” unit cell) consisting of the basic object and a number of subobjects.Figure1 shows four situations as characteristic examples, out of a literary infinite number of possible combinations. In Figure 1(a), both the basic object and the substructures are circles; in Figure 1(b), both are hexagons. A more general case is shown in Figure 1(c), where the main object is a square, while the subobjects are octagons. The most interesting example from the practical point of view is in Figure 1(d), where the same basic shape, a square, is replicated three times to form additional deeply subwavelength corners formed by simple overlap of the main object.Examples of different aperture forms consisting from “primitive” subobjects. (a) Object and subobjects are circles; (b) object and subobjects are hexagons; (c) object is square, subobjects octagons; (d) only main object is replicated, subobject features are generated by intersecting these replicas.
(a)
(b)
(c)
(d)One could continue to add further subobjects, making additive superstructures. Any desired form could be thus defined using simple primitive objects, making the observed situation a geometrical analog to the case of expansion of a function into series (in our case, a complex object that is expanded into primitives).The existence of the superposed object will change the resulting geometry and thus introduce nonlocalities in its electromagnetic response. These nonlocalities will be caused by several factors, among which is the geometrical shape itself, since larger field concentrations are expected on sharper corners. Another factor is coupling between different subobjects, since narrower gaps between the subobjects result higher field localizations.In this work we choose to analyze the case when rectangular subobjects are added to rectangular basic structures. Figure2, left side, shows the case of an unperturbed fishnet: this is a double fishnet with rectangular nanoapertures. Figure 2, right side, shows the case when smaller rectangular subobjects are added to each corner of the basic structure.Figure 2
Square apertures versus compound apertures.It has been already mentioned that the effective medium theory in principle does not take into account the nonlocal nature of plasmonic resonances in complex geometries like those in the compound unit cells described here. Still, it is possible to introduce the effective parameters even for such structures as descriptive figures of merit that give insight into the spectral dispersion of their electric and magnetic response. Thus regardless of the validity of the EMA the effective parameters have a physical meaning that relates with the measured or simulated scattering parameters of the structure and helps in predicting its spectral behavior in different situations.


### 2.2. Extraction of Effective Electromagnetic Parameters

Here we utilize direct extraction of effective constitutive parameters (complex dielectric permittivityε, permeability μ, and refractive index n) from the known complex transmittance (t) and reflectance (r) spectra, as proposed in [21] and described in [22]. If k
0 is the free space wave number and d is the total thickness of the fishnet (top and bottom layer plus dielectric), then the wave impedance that can be calculated as
(1)
z
=
±
1
+
r
+
t
exp
⁡
j
k
0
d
1
+
r
-
t
exp
⁡
j
k
0
d
1
-
r
-
t
exp
⁡
j
k
0
d
1
-
r
+
t
exp
⁡
j
k
0
d
,
while the effective refractive index is
(2)
n
=
Im
⁡
ln
⁡
⁡
e
j
n
k
0
d
+
2
m
π
k
0
d
-
j
Re
ln
⁡
⁡
e
j
n
k
0
d
k
0
d
,
where m is an integer that represents the branch index of the real part of n and the exponential part is
(3)
exp
⁡
j
n
k
0
d
=
t
exp
⁡
j
k
0
d
1
-
r
z
-
1
z
+
1
-
1
.The values ofε and μ are further calculated from the obtained wave impedance and refractive index. All expressions except the first have multiple branches and it is necessary to choose the right one with a proper sign. For lossy materials without active media the conditions are Re
{
z
}
>
0 and Im
⁡
{
n
}
>
0.The figure of merit of a metamaterial describes its absorptive losses and is defined as the ratio between the real and imaginary part of refractive index [18]:
(4)
FOM
=
Re
n
Im
⁡
(
n
)
.


### 2.3. Structures

We consider several cases. Figure3(a) shows the basic unit cell: a cruciform structure where a quarter of an aperture is located in each corner. Figure 3(b) shows the structure perturbed by four 30 nm × 30 nm subobjects, while Figure 3(c) shows 60 nm × 60 nm subobjects and Figure 3(d) shows 30 nm × 60 nm subobjects. Each rectangular subobject is placed so that the intersection of its diagonals coincides with a corner of the main object, while the sides of the main object and the subobjects are parallel.Types of unit cells considered; (a) basic (unperturbed) shape—simple rectangular apertures; (b) compound aperture (super unit cell) obtained by adding 30 nm × 30 nm square patches on each corner of the basic shape; (c) super cell obtained by adding four 60 nm × 60 nm squares; (d) supercell obtained by adding four 30 nm × 60 nm rectangles.
(a)
(b)
(c)
(d)The dimensions of the basic structure were as follows: the aperture dimensions are 60 nm × 80 nm, square unit cell with a 300 nm side, the top and the bottom metal layer are 30 nm thick, and the thickness of the dielectric part is 40 nm.It has been assumed that the metal part was gold, well described by Drude model. Material parameters for gold were taken as follows: plasma frequencyω
p
=
7.175 PHz and collision frequency γ
=
6.5 THz [23]. The dielectric part was aluminum oxide with a refractive index of 1.63.


### 2.4. Numerical

We utilized Comsol Multiphysics software package (RF module) to calculate the scattering parameters of our fishnets. The unit cell is treated as a waveguide surrounded by air. We applied the Perfect Electric Conductor and the Perfect Magnetic Conductor boundary conditions on the surfaces perpendicular to the respective incident fields. Thus we simulated infinite periodic structure. Two ports were added above and below the unit cell, each of them parallel to the surface. One port was active and incident light entered the waveguide through it. Parametric sweep of incident light wavelength was used to determine the spectral dependence.


## 2.1. Basic Concept

The basic element of a fishnet structure, whether single-metal layer (extraordinary optical transmission array) [18], double layer (double fishnet) [19], or multiple [20], is a single subwavelength aperture. This aperture may have different shapes, including circular, elliptical, and rectangular. One can define a unit cell for a fishnet by cutting four neighboring apertures through the middle and thus defining a cross-like metal structure where a quarter of aperture is located in each corner.In this work we modify the basic element by superposing additional apertures onto those already existing. We denote the basic aperture by the unperturbed structure, while the additional apertures are subobjects, substructures, or “patches.” Both the unperturbed structure and the added substructures are simple geometrical forms like circles, rectangles, and so forth, denoted by the primitive objects or primitives. Thus one obtains a compound unit cell (a “super” unit cell) consisting of the basic object and a number of subobjects.Figure1 shows four situations as characteristic examples, out of a literary infinite number of possible combinations. In Figure 1(a), both the basic object and the substructures are circles; in Figure 1(b), both are hexagons. A more general case is shown in Figure 1(c), where the main object is a square, while the subobjects are octagons. The most interesting example from the practical point of view is in Figure 1(d), where the same basic shape, a square, is replicated three times to form additional deeply subwavelength corners formed by simple overlap of the main object.Examples of different aperture forms consisting from “primitive” subobjects. (a) Object and subobjects are circles; (b) object and subobjects are hexagons; (c) object is square, subobjects octagons; (d) only main object is replicated, subobject features are generated by intersecting these replicas.
(a)
(b)
(c)
(d)One could continue to add further subobjects, making additive superstructures. Any desired form could be thus defined using simple primitive objects, making the observed situation a geometrical analog to the case of expansion of a function into series (in our case, a complex object that is expanded into primitives).The existence of the superposed object will change the resulting geometry and thus introduce nonlocalities in its electromagnetic response. These nonlocalities will be caused by several factors, among which is the geometrical shape itself, since larger field concentrations are expected on sharper corners. Another factor is coupling between different subobjects, since narrower gaps between the subobjects result higher field localizations.In this work we choose to analyze the case when rectangular subobjects are added to rectangular basic structures. Figure2, left side, shows the case of an unperturbed fishnet: this is a double fishnet with rectangular nanoapertures. Figure 2, right side, shows the case when smaller rectangular subobjects are added to each corner of the basic structure.Figure 2
Square apertures versus compound apertures.It has been already mentioned that the effective medium theory in principle does not take into account the nonlocal nature of plasmonic resonances in complex geometries like those in the compound unit cells described here. Still, it is possible to introduce the effective parameters even for such structures as descriptive figures of merit that give insight into the spectral dispersion of their electric and magnetic response. Thus regardless of the validity of the EMA the effective parameters have a physical meaning that relates with the measured or simulated scattering parameters of the structure and helps in predicting its spectral behavior in different situations.


## 2.2. Extraction of Effective Electromagnetic Parameters

Here we utilize direct extraction of effective constitutive parameters (complex dielectric permittivityε, permeability μ, and refractive index n) from the known complex transmittance (t) and reflectance (r) spectra, as proposed in [21] and described in [22]. If k
0 is the free space wave number and d is the total thickness of the fishnet (top and bottom layer plus dielectric), then the wave impedance that can be calculated as
(1)
z
=
±
1
+
r
+
t
exp
⁡
j
k
0
d
1
+
r
-
t
exp
⁡
j
k
0
d
1
-
r
-
t
exp
⁡
j
k
0
d
1
-
r
+
t
exp
⁡
j
k
0
d
,
while the effective refractive index is
(2)
n
=
Im
⁡
ln
⁡
⁡
e
j
n
k
0
d
+
2
m
π
k
0
d
-
j
Re
ln
⁡
⁡
e
j
n
k
0
d
k
0
d
,
where m is an integer that represents the branch index of the real part of n and the exponential part is
(3)
exp
⁡
j
n
k
0
d
=
t
exp
⁡
j
k
0
d
1
-
r
z
-
1
z
+
1
-
1
.The values ofε and μ are further calculated from the obtained wave impedance and refractive index. All expressions except the first have multiple branches and it is necessary to choose the right one with a proper sign. For lossy materials without active media the conditions are Re
{
z
}
>
0 and Im
⁡
{
n
}
>
0.The figure of merit of a metamaterial describes its absorptive losses and is defined as the ratio between the real and imaginary part of refractive index [18]:
(4)
FOM
=
Re
n
Im
⁡
(
n
)
.


## 2.3. Structures

We consider several cases. Figure3(a) shows the basic unit cell: a cruciform structure where a quarter of an aperture is located in each corner. Figure 3(b) shows the structure perturbed by four 30 nm × 30 nm subobjects, while Figure 3(c) shows 60 nm × 60 nm subobjects and Figure 3(d) shows 30 nm × 60 nm subobjects. Each rectangular subobject is placed so that the intersection of its diagonals coincides with a corner of the main object, while the sides of the main object and the subobjects are parallel.Types of unit cells considered; (a) basic (unperturbed) shape—simple rectangular apertures; (b) compound aperture (super unit cell) obtained by adding 30 nm × 30 nm square patches on each corner of the basic shape; (c) super cell obtained by adding four 60 nm × 60 nm squares; (d) supercell obtained by adding four 30 nm × 60 nm rectangles.
(a)
(b)
(c)
(d)The dimensions of the basic structure were as follows: the aperture dimensions are 60 nm × 80 nm, square unit cell with a 300 nm side, the top and the bottom metal layer are 30 nm thick, and the thickness of the dielectric part is 40 nm.It has been assumed that the metal part was gold, well described by Drude model. Material parameters for gold were taken as follows: plasma frequencyω
p
=
7.175 PHz and collision frequency γ
=
6.5 THz [23]. The dielectric part was aluminum oxide with a refractive index of 1.63.


## 2.4. Numerical

We utilized Comsol Multiphysics software package (RF module) to calculate the scattering parameters of our fishnets. The unit cell is treated as a waveguide surrounded by air. We applied the Perfect Electric Conductor and the Perfect Magnetic Conductor boundary conditions on the surfaces perpendicular to the respective incident fields. Thus we simulated infinite periodic structure. Two ports were added above and below the unit cell, each of them parallel to the surface. One port was active and incident light entered the waveguide through it. Parametric sweep of incident light wavelength was used to determine the spectral dependence.


## 3. Results and Discussion

Figure4 shows the transmittance and reflectance for the unperturbed aperture shape. A sharp reflectance dip is observed between 800 and 900 nm, as well as a high reflectance band between 600 nm and about 820 nm.Figure 4
Transmittancet and reflectance r for the unperturbed (basic) shape.In the case when a 30 nm × 60 nm subobject is added to each corner of the basic shape, a drastic change is observed in the spectra (Figure5). A dramatic rise of transmission and drop of reflection occur where a reflection band existed. Basically a new transmission window occurs within the plasmonic bandgap.Figure 5
Transmittancet and reflectance r for the 30 nm × 60 nm squares on the corners of the basic shape.Similar situation, but less drastic, is observed for square subobjects. For 30 nm × 30 nm subobject (Figure6) transmission is increased, but its overall level decreases with the wavelength. For 60 nm × 60 nm (Figure 7) the situation is very similar, but the transmission increases with the wavelength.Figure 6
Transmittancet and reflectance r for the 30 nm × 30 nm squares on the corners of the basic shape.Figure 7
Transmittancet and reflectance r for the 60 nm × 60 nm squares on the corners of the basic shape.Figure8 shows the spectral absorption for the unperturbed case and the 30 nm × 60 nm case. A distinct redshift is observed in the case with superposed subobjects, stronger for shorter wavelengths.Figure 8
Absorption coefficients for the unperturbed aperture (A
0, blue line) and for the one with 30 nm × 60 nm subobjects (A
m, red line).A comparison of Figure4 with Figures 5–7 shows that relatively small deviations from the unperturbed object cause large changes of the resulting spectral characteristics. Such sensitivity is convenient for the use of metasurfaces and 2D plasmonic crystals in sensing. However, it requires a careful design of the structures and accurate fabrication methods when aiming at specific values of the resulting optical parameters. It also points out that inaccurate fabrication and the appearance of spatial noise may completely wipe out the desired optical properties, although the effective values of optical parameters will be practically unchanged. The reason for this is again the appearance of nonlocalities that may become very strong in some geometries.We analyzed the dependence of the optical response of the structure to perturbations of the basic geometry of the unit cell and summarized our results in Figure9. The unperturbed structure is identical to that used to obtain Figures 4–8, but the dimensions of the corner patches are now varied. The starting corner patch used is 30 nm × 60 nm. We multiplied the sides of this corner patch with a scaling factor a that varied from 0.1 (furnishing a corner patch with the dimensions 3 nm × 6 nm) to 1.2. We can see that the highly reflective band between 600 nm and 800 nm is gradually replaced by a transmission band in the same wavelength range. At the same time, the wavelengths of its edges remain practically constant.Figure 9
Spectral reflection coefficient of the compound unit cell-based double fishnet for rectangular corner patches. The basic patch is 30 nm × 60 nm, and the scaling factor for its both side lengths isa.Figures10–12 show the optical constants extracted by using the procedures outlined in [21, 22]. Figure 10(a) shows the extracted relative dielectric permittivity for the basic structure as shown in Figure 3(a), while Figure 10(b) shows the extracted permittivity for the structure with 30 nm × 60 nm subobjects, shown in Figure 3(d). Figure 11 shows relative magnetic permeability and Figure 12 extracted effective refractive index, all for the same structures.Extracted effective relative dielectric permittivity for the basic structure (a) and for the one with 30 nm × 60 nm subobjects as in Figure3(d) (b).
(a)
(b)Extracted effective relative magnetic permeability for the basic structure (a) and for the one with 30 nm × 60 nm subobjects as in Figure3(d) (b).
(a)
(b)Extracted effective refractive index for the basic structure (a) and for the one with 30 nm × 60 nm subobjects as in Figure3(d) (b).
(a)
(b)Figure10(a) shows that in the unperturbed situation the real part of the relative dielectric permittivity (notation: Re) becomes negative for wavelengths above 600 nm. At the resonant frequencies between 800 nm and 900 nm there is the characteristic increase of the absorption losses, as seen from the peak in the imaginary component of the permittivity (notation Im). This behavior is actually similar to that of metallic structures, however contrary to metals; it is possible here to tune the spectral dispersion through choosing the aperture geometry. This is seen in Figure 10(b), where the peaks of both the real and imaginary component are shifted to lower wavelengths.Figure11 shows the real and imaginary component of effective permeability. This value should be understood as a parameter describing high-frequency magnetic response of the metal-dielectric nanocomposite under consideration. Here too one can observe negative values of the real part and an increase of the imaginary component connected with higher losses.Figure12 shows the effective refractive index as the synthetic parameter obtained from the effective permittivity and permeability. Since there is a range of overlapping negative values of the real parts of ε and μ, this means that there will be a range of negative effective refractive index values; that is, the structure will behave as negative index metamaterial. This behavior is qualitatively preserved in the modified structure as well. However, one should note a beneficial property of the modified structure: it has lower losses (i.e., smaller imaginary part of effective refractive index) compared to the unmodified one. This can be also seen in Figure 13 that shows the comparison of figures of merit between the unperturbed case and the case perturbed with 30 nm × 60 nm rectangles. The value of the obtained figure of merit is about 50% larger than in the unperturbed case.Figure 13
Comparison of figures of merit for the basic structure (dashed blue lines) and for the one with 30 nm × 60 nm subobjects (solid red line).It is important to mention here some technological aspects of the described approach. Since the main objects have subwavelength dimensions, the subobjects must be deeply subwavelength, and both are intended for fabrication by the top-down approach. With the main objects in the dimension range of several hundreds of nanometers, the subobjects have dimensions of tens of nanometers. Having in mind the limitations of the state-of-the-art lithographic equipment, it is obviously rather difficult to reproduce such shapes.We argue that one possible solution is to generate supercells using identical overlapping objects, an example being the one shown in Figure1(d). A square with the dimensions comparable to the wavelength and thus more easily delineated is replicated there three times, each time with a slight misplacement. The newly created compound structure contains new deeply subwavelength features created at the intersections of the overlapping “big” objects. These new features ensure new hotspots of electromagnetic field and thus cause reshaping of the optical spectra.The procedure obviously cannot create objects like those in Figures1(a) to 1(c), which contain smaller subobjects. However, it does create sharp corners as a result of addition of identical objects. Alternatively, it can replicate the main object twice at a subwavelength distance, thus ensuring the appearance of a field hotspot within the gap.The main object (which can also be circular, etc.) could be thus used to pattern a photolithographic mask at the maximum available resolution of the existing equipment. This resolution is obviously insufficient to get any deeply subwavelength features. As the next step, however, we use alignment equipment to simply shift the mask and repeat pattern definition on the substrate. In this way the shapes like those in Figure1(d) can be generated, ensuring the appearance of sharp corners. Thus we reduce the problem of resolution of lithographic masks to the problem of alignment of the existing lower resolution masks.


## 4. Conclusion

We constructed 2D arrays of nanoapertures in double fishnet plasmonic configuration with complex 2D shapes on deep subwavelength level. To define the aperture geometry we applied a simple additive approach, where an initial form is modified by overlapping it with additional objects/subobjects. In this manner one can use a set of primitives to generate additional nonlocalities without increasing the complexity of the system or requiring finer resolutions, which contributes to a lower fabrication cost. We described a possible approach to the fabrication of some of such objects, which is to simply shift the same photolithographic mask and repeat its pattern, where field localizations could be expected near the deeply subwavelength corners at the intersection points of the “large” objects.We show that the nonlocalities introduced by overlapping simple objects can be used to produce large deviations of spectral dispersion from the unperturbed case even in the case of relatively minor additive modifications of the basic geometry.Our approach points out to the possibility to tailor the dispersion and the resonance spectra for plasmonic nanocomposites by adjusting the local geometry and exploiting localized interactions at a subwavelength level. Additional resonances can be introduced and spectral dispersion can be shifted and tuned up to a desired level. The obtained freedom in tailorability could be used to fabricate structures for plasmonic sensing of chemical or biological analytes and generally for the fabrication of planar waveguides utilizing two-dimensional plasmonic crystals.



---

*Source: 312064-2015-01-31.xml*