Retrieval of Effective Parameters of Subwavelength Periodic Photonic Structures

: We revisit the standard Nicolson–Ross–Weir method of effective permittivity and permeability restoration of photonic structures for the case of subwavelength metal-dielectric multilayers. We show that the direct application of the standard method yields a false zero-epsilon point and an associated spurious permeability resonance. We show how this artifact can be worked around by the use of the cycle shift operator to the periodic multilayer in question.


Introduction
Determining the effective macroscopic parameters of a slab containing a photonic structure is often both highly relevant both from a practical and from a fundamental point of view, providing a link between the actual microscopic composition of the structure and its macroscopic optical response.The most common ways to accomplish this task are various types of ellipsometry [1,2] or reflectance and transmittance measurements [3].
However, there are two major issues with such parameter restoration techniques.One is the fundamental applicability limit of the effective medium approximation approach.It is normally mitigated by considering photonic structures with feature sizes much smaller than the wavelength (deeply subwavelength structures), although an exception exist in special cases, such as hyperbolic metamaterials [4].The other issue is related to the artifacts of the restoration procedure itself.For example, it is known that the constitutive parameters of a slab are hard to obtain nearby the total transmission frequency or when we approach the magnetic or electric resonances of a structure [5].These difficulties can be overcome with a modified extraction procedure, based upon a generalized homogenization technique [6].
In this work, we address another difficulty arising when extraction is to be carried out at the point of total reflection.It is shown that application of the widely used standard Nicolson-Ross-Weir procedure [7] results in spurious features near this point, such as false zero-epsilon points and permeability poles.In the optics of subwavelength structures, total reflection is often associated with the true ε-near-zero regime [8].Therefore, being able to make a correction for artifacts in such a regime becomes of particular importance to determine the correct operational wavelengths of epsilon-near-zero behavior and the optical mechanism of subwavelength photonic structures.
The rest of the paper is organized as follows.In Section 2, we demonstrate the artifacts associated with the application of the standard restoration procedure to subwavelength metal-dielectric periodic photonic multilayers.In Section 3, we trace these artifacts to the total reflection from a photonic metamaterial slab at the point where artifacts are observed.In Section 4, we show how to amend the artifact behavior and ensure correct parameter restoration by applying the cycle shift operator to the multilayer in question.Finally, in Section 5, we summarize the results.

Spurious Features in the Nicolson-Ross-Weir Restoration Procedure
We consider periodic layered structures formed by alternating layers of metal and dielectric.It was shown previously that the nonlocal permittivity tensor of such plasmonic multilayers generally possesses off-diagonal components [9,10], which only vanish when light propagates normal to layer interfaces.Furthermore, in this case, no more than one eigenmode corresponds to each frequency.Hence, we will consider the normal incidence of light throughout the paper and will aim at the restoration of the longitudinal component of the permittivity tensor from the reflection r and transmission t coefficients.
We apply the standard Nicolson-Ross-Weir extraction method [7] to a finite multilayer slab of thickness D and calculate its effective permittivity ε and permeability µ as follows.First, the impedance Z and the refraction index n are given by: where for the given wavelength λ, one should choose the signs and the free parameter m in a proper manner (see [11,12]).Then, since Z = √ µ/ε, n = √ εµ, the effective permittivity and permeability can be determined as ε = n/Z, µ = Zn.
It can be seen that there is a singular point in Equation (1) where r = |r|e iϕr = 1.Therefore: If we suppose that R = 1 and T = 0, then Equation ( 1) results in Z = ±∞.Since Z = ± √ µ/ε, it follows that infinite Z implies either that µ must approach infinity or that ε must approach zero.On the other hand, the refractive index n = √ εµ remains finite as per Equation ( 1); therefore, we see both µ → ∞ and ε → 0 at the same time.
Therefore, at first glance, there should always be an epsilon-near-zero, mu-near-pole response for the case of total reflection.Figure 1 shows an example of such a situation, where we observe a spurious permeability resonance.We also see a poor match between the correct effective parameters ε eff , µ eff and extracted ε.Re(µ) Because of this apparent extraction artifact, which can be easily misinterpreted as optical magnetism in plasmonic multilayers, special care should be taken in applying the simple extraction procedure described above.In what follows, we address the subject in detail and show how the artifacts can be worked around using a cycle shift geometrical transformation on the multilayer structure.

Detailed Analysis of the Total-Reflection Point in Multilayers
We can characterize an arbitrary multilayer by the 2 × 2 transfer matrix M = m jk that connects the tangential components of propagating electric and magnetic fields.Let us place one unit cell of multilayer having thickness L = l m + l d (l d,m are the thicknesses of the dielectric and metal layers) between input and output host media having impedances of ξ i,o = ε −1/2 i,o .Knowing M and assuming normal incidence, the reflection coefficient from the single period is defined as the ratio of amplitudes of reflected (H − ) and incident (H + ) magnetic fields in the case of normal incidence [13]: Eigenmodes e ikxL of M provide us with the Bloch wave numbers x of the structure.Based on the knowledge of M 's eigenmodes, we can calculate the transfer matrix for a slab composed of N unit cells (with the total thickness D = N L) simply by raising M to the N -th power.Taking into account that M is a unimodular matrix, we have [14]: where are Chebyshev's polynomials.Now, it is possible to write down the reflection coefficient for the slab as follows: Now, we are ready to write down the reflection coefficient for a semi-infinite multilayer.It can be noticed that for complex k x : In the case of a semi-infinite structure, we are not interested in the output interface, and it can be assumed that We see that the total reflection occurs when: If the host medium is free space, then: i.e., m 11 + m 12 being the eigenvalue of M means that r ∞ = 1.
Here, we are concerned with multilayers formed by periodically alternating layers of two types, one made of a dielectric and the other made of metal.For this case, M is of the following form: where ) is the wave number in the respective layers.As is the usual meaning of metal that ε m < 0, β m is purely imaginary: resulting in: Based on Equations ( 13) and ( 14), we see that m 11 is always real, while m 12 is imaginary for a wave propagating across the layers.Separating the real and imaginary parts in Equation ( 11), we immediately obtain that ℜ[e ikxL ] = m 11 .As a result, defines the frequency where we observe total reflection.We note that since |r| 2 + |t| 2 = 1 in the case of total reflection, there is no transmission through the structure, regardless of whether the structure is infinite or finite.Consequently, we will obtain false material parameters for a finite structure at this point, as can be seen from Figure 1.

Obtaining Correct ε and µ Using the Cycle Shift Operator
However, there is a way to correct permittivity ε eff , and permeability µ eff can be obtained from reflection and transmission coefficients.We recall that multilayered metal-dielectric structures have a strong inherent spatial dispersion and are typically described by a nonlocal dielectric permittivity tensor ε(ω, k) [9].However, as we only consider the normal incidence of light, we can move from the spatial dispersion framework to a ε(ω) − µ(ω) description of the system [15,16].In such a paradigm, we can write the following expressions for the effective permittivity and permeability [17], valid for the case of normal incidence: where ε || (ω, 0) is the longitudinal component of ε(ω, 0), and Lagrange's notation is used for the derivative.Since we are obtaining ε eff and µ eff from the nonlocal effective model, they are much more accurate than the values yielded in the quasi-static limit [10].
We see that the two artifacts in the parameter retrieval seen in Figure 1 are connected.It is the mismatch in the retrieved ε that causes the faulty point where ε = 0 appears in the spectrum, which is where the pole in the effective µ appears.Calculations show that it can be traced to the mismatch in the retrieved impedance, Z, with respect to the Bloch impedance in the structure [18].
The method to avoid the appearance of the artifact features is based on a cycle shift operator applied to the multilayer structure.Taking a finite multilayer composed of N periods, this operator acts like a frame of length N L, moving across an infinite multilayer and cutting out the region that it covers, as depicted in Figure 2a.As can be seen in Figure 2b, different values of frame shift cause a significant difference in the phase of the reflection coefficient for a semi-infinite structure, corresponding to changing of the retrieved impedance to match the Bloch impedance.This, in turn, makes it possible to pick the value of the phase that would result in correct values of the extracted parameters.
For our example with l d = 30 nm and l m = 60 nm and for the structure containing five periods, we start from the dielectric layer and apply the cycle shift operator, shifting the frame with a step of 1 nm.Calculating r and t at each step and extracting material parameters, we find that the best correspondence of the extracted ε and µ with the accurate values ε eff and µ eff is achieved when the frame is shifted by 15 nm or 60 nm, making the unit cell symmetric.Figure 3 demonstrates the extraction procedure for a 15-nm frame shift, showing a much better matching with the theoretically calculated values, as well as the disappearance of the artifact features in comparison with the extraction without the frame shift operator application (Figure 1). Figure 3. Retrieved permittivity and permeability when the shift operator with a shift value of 15 nm applied in order to move the point of total reflection to the edge of the forbidden zone.The rest is as in Figure 1.Physically, the "magic" value of the shift corresponds to positioning the point where r = 1 (implying that ϕ r = ±2πm, m ∈ Z) at the edge of the stop band shown in Figure 2c.One can also explain why the cycle shift operator restores the validity of the extraction procedure by noting that the impedance mismatch (the difference between the wave impedance, which enters the extraction procedure, and the Bloch impedance, which is an attribute of the infinite periodic multilayer) can result from a poorly chosen termination of the multilayer.Namely, poor termination causes high impedance mismatch between the outer space and the multilayer and, accordingly, the associated strong reflection, leading to artifacts.Applying the cycle shift operator effectively cycles through all possible terminations, and the best correspondence in parameter extraction happens when a good impedance matching between the multilayer and the outer space is achieved.Our results agree with previous accounts [19,20], where it was shown that a symmetric unit cell facilitates a good matching between the effective medium and extracted parameters.

Conclusions
We have demonstrated that the singularity in impedance originating from r = 1 and the zero value of its phase leads to artifacts in the retrieval of constitutive parameters of subwavelength photonic metal-dielectric multilayers, namely spurious magnetic resonance and a large mismatch of effective permittivity.We have shown that these artifacts can be worked around by applying the cycle shift operator to the multilayer.The obtained results are useful for accurate experimental determination of effective macroscopic parameters of photonic multilayers and metamaterials, especially those that work in the epsilon-near-zero regime and those that are designed to exhibit artificial optical magnetism.The proposed procedure is applicable to lossless and lossy systems within the limits assessed by the standard Nicolson-Ross-Weir technique.These results can also be used to better understand why some metamaterial geometries lend themselves to homogenization and extraction-based parameter retrieval much more readily than some others [21].

Figure 1 .
Figure 1.A case of spurious constitutive parameter retrieval.(a) The reflection coefficient from a slab formed by ten alternating layers of dielectric and metal having thicknesses of l d = 30 nm and l m = 60 nm, respectively; (b) reflection coefficient from a semi-infinite structure formed by the same layers; (c) retrieved permittivity and (d) permeability, along with theoretical predictions.The material parameters are the same as those used throughout the paper.

Figure 2 .
Figure 2. (a) Diagram of how the cycle shift operator acts.A frame of chosen length N L, N ∈ N >0 , is cutting out the finite fragment from the infinite multilayer; (b) The reflection coefficient from the semi-infinite structure with (c) the corresponding band diagram of the multilayer.The allowed band is highlighted in (b,c).