Design a Stratiform Metamaterial with Precise Optical Property

: In this work, a stratiform metamaterial is arranged as multiple periods of metal-dielectric symmetrical film stack to provide precise equivalent refractive index and admittance. There are multiple solutions of equivalent refractive index retrieved from the characteristic matrix of the film stack. The correct refractive index is derived by connecting different branches of solution at different ranges of wavelength or thickness of the dielectric layer. The refractive index of an Ag-TiO 2 five-layered symmetrical film stack shown in previous work is demonstrated to be positive real instead of negative real. The associated type I iso-frequency curve supports negative refraction. In order to extend the operating wavelength of type I metamaterial, the number of the metal-dielectric symmetrical film stack is increased to reduce the thickness of the dielectric film to approach subwavelength requirement.


Introduction
A metamaterial with a negative index of refraction has been predicted to be effective for constructing a perfect lens [1]. Archetypical left-handed metamaterials comprise split rings as inductive-capacitive resonators, which support the simultaneous reversal of electric and magnetic responses to light at microwave frequencies [2]. Using advanced lithographic nanofabrication, the scales of resonators were shrunk to exhibit a left-handed response to light at higher frequencies, such as infrared or visible frequencies [3]. Recently, a stratiform metamaterial that comprises metal film (M) and dielectric film (D) was developed for transverse magnetic (TM) polarization mode [4]. Light incident on the side of this stack will couple to waveguide modes [5,6]. The anti-symmetry mode of the transverse electric field in a plasmonic waveguide composed of metal film (M)/dielectric film (D)/metal film (M) (MDM) causes a negative group velocity at frequencies greater than the surface plasmon resonance frequency [6]. A combination of a pair of typical MDM waveguides as a waveguide composed of metal film (M)/dielectric film (D)/metal film (M)/dielectric film (D)/metal film (M) (MDMDM) results in electric-field coupling between adjacent MDM waveguides. The symmetrical electric field distribution leads to a positive coupling constant, thus the stack exhibits backward wave phenomenon.
In 2014, a film stack comprising TiO2 and Ag films alternately was claimed to yield an index of negative unity over a wide range of incidence angles at a wavelength of 363.8 nm [7]. Later, whether the negative refraction in the TiO2-Ag multilayer came from a negative index of refraction was questioned and discussed. The negative refraction also occurs in a hyperbolic metamaterial. Hyperbolic metamaterials have been developed for more than 10 years. In most cases, hyperbolic metamaterials are metal-dielectric multilayers [8,9]. Its property is described by a diagonal permittivity tensor. The principal permittivities are roughly estimated using effective medium approximation [10]. The signs of the two tangential permittivities are the same, but the signs of its tangential and vertical permittivities are opposite. Accordingly, there are two kinds of hyperbolic shaped isofrequency curves (type I and type II) [11]. A previous work on the decomposition of the Bloch wave in the layered structure into a number of harmonics showed that the iso-frequency curve of the main harmonic mode corresponds to the type I hyperbolic metamatrerial [12]. The negative refraction of energy is invoked by the shape of the iso-frequency contour [11]. Therefore, the multilayer is a hyperbolic metamaterial, and light propagation is mainly dominated by right-handed harmonics with a positive phase index.
Although the relevant works have explained the wave phenomenon in a metal-dielectric multilayer, a comprehensive design method for negative refraction is still required. A multilayer with precise optical constants, including refractive index and admittance, are necessary for a metamaterial design because, in real applications, both transmission and reflection are affected by both optical constants. In this work, a stratiform metamaterial is designed by arranging a multilayer composed of multiple periods of metal-dielectric symmetrical film stack (SFS). Based on the calculation of characteristic film matrixes [13], the metal-dielectric symmetrical film stack is equivalent to a layer with the equivalent refractive index and admittance that can be derived precisely by choosing the correct solution from a multi-branch of solution [14,15]. The refractive index versus angle of incidence can be used to plot a precise iso-frequency curve that indicates whether the stratiform metamaterial is a hyperbolic metamaterial. A type I hyperbolic metamaterial supports negative refraction from normal incidence to a certain oblique angle of incidence. In order to apply to real applications, the transmission is taken into account by considering the extinction coefficient and admittance [16].

Methods
Exploiting thin-film optics, symmetrical film stacks have been widely designed for edge filters [17]. An all-dielectric symmetrical film stack is equivalent to a single layer with an equivalent refractive index that is a function of wavelength , thicknesses, and refractive indexes of all constituent films. The equivalent refractive index is pure real and pure imaginary in different wavelength ranges, providing passband and stopband characteristics, respectively. An edge filter was designed across the passband and stopband. The equivalent model of a metal-dielectric symmetrical film stack was still valid, thus the product of characteristic film matrixes of the constituent thin films yielded a resultant matrix that had the same form of that of a single layer.
where δ and η (j = A, B) are the phase and admittance of each layer of the stack. The equivalent phase thickness of the stack γ and the equivalent admittance Eeq were complex. At normal incidence, the equivalent refractive index Neq is derived from the phase thickness γ through the relationship γ = 2 (2 + )/ . Both Eeq and Neq can be tailored separately. Both Eeq and Neq are represented with the equivalent permittivity εeq and the equivalent permeability μ as = / and N = , respectively. According to effective medium approximation, the equivalent refractive index and admittance of a metal-dielectric composite are generally complex. In case of oblique incidence and TM polarization state, the equivalent phase thickness should be = 2 (2 + ) / . Figure 1 shows the equivalent scheme of a 3-layered symmetrical film stack ABA. The θeq is the equivalent angle of refraction in the film stack. The and are thicknesses of film A and film B, respectively. The and are angles of refraction in film A and film B, respectively. The equivalent admittance for TM polarization should be modified as = / . The product of the 3 characteristic film matrixes ABA is shown in Equation (2).
The permittivity εeq can be retried from the matrix elements firstly through Equation (3).
The equivalent angle of refraction can be derived through the law of refraction = (N0 and θ0 are the refractive indexes of the cover medium and angle of incidence, respectively), as shown in Equation (4).
Next, we can get the equivalent refractive index Neq shown in Equation (5).
Finally, we can have the equivalent permeability μeq and admittance In this work, the time-dependent phase factor of a propagating harmonic wave takes the form (ω is the angular frequency) and an absorbing material has a refractive index n-ik. Unlike an alldielectric symmetrical film stack, most metal-dielectric symmetrical film stacks have a different equivalent refractive index and equivalent admittance because the equivalent permeability is not unity for most metal-dielectric composites in subwavelength scale. Both and can be tailored separately for novel applications [16]. In a previous work, a 7-layered metal-dielectric symmetrical film stack was designed and fabricated to have an near unity and an with large imaginary part that was comparable to that of a metal. Such a layered metamaterial can couple most incident light into the upper interface and dissipate its energy with a very small thickness [18].
It is noted that there were multi-solutions for the phase thickness and refractive index. There were multiple branches of solutions of equivalent phase thickness retrieved from matrix element M11:  [19]. It is particular that a specific branch was unable to satisfy the criteria over the whole wavelengths. The continuity of the real part of the equivalent refractive index as a function of wavelength or any constitutive parameter required choosing different branches for regions separated by the discontinuous points. For a metal-dielectric SFS in the form of MDM…DM, a m-th branch was chosen thus that the at thickness of dielectric layer = 0 was equal to the refractive index of metal. The m-th branch was correct for the range form = 0 to its first discontinuous point at d = .
The at thickness d larger than relied on a branch next to the m-th branch to connect the m-th branch at d = to keep the continuous condition till its discontinuous point at d = .
The correct branch for the thickness d larger than was chosen to keep continuous at the = .

Results
A five-layered SFS of Ag (33 nm)/TiO2 (28 nm)/Ag (30 nm)/TiO2 (28 nm)/Ag (33 nm) in a previous work [7] was examined here. Such SFS was claimed to have a negative index around -1 at a wavelength of 363.8 nm. The refractive indexes of Ag and TiO2 were 0.0785-1.59i and 2.8-0.05i, respectively. The equivalent versus thickness of TiO2 from 0 nm to 60 nm at normal incidence with branch (m = 0) is shown in Figure 2a. At the thickness of d = 0 nm, the index of refraction was 0.0785, which was the real part of the refractive index of Ag. It means that the branch m = 0 was correct from d = 0 nm to its first discontinuous point at = 22 nm. At the thickness of d = 28 nm in Figure 2a, the index of refraction was −1.058, which was the proposed value of negative index. However, the index of refraction after = 22 nm needed to be connected with another branch (m = 1), as shown in Figure 1b. For the branch (m = 1), there was a discontinuous point at = 37 nm, thus the index after d = 37 nm was connected with branch (m = −1), as shown in Figure 1c. Figure 1d shows the correct real part of as a function of d from d = 0 nm to d = 60 nm.  Figure 3a shows both the real part and imaginary part of as a function of d. The at d = 28 nm was 1.34 and the versus the angle of incidence is shown in Figure 3b. With the correct as a function of angle of incidence. The iso-frequency curve (IFC) in Figure 3c represents the wave vector component versus the other component . The wave vector component along the interface is denoted as that is decided by the incident wave vector and continuous across each interface of a multilayer. The other wave vector component perpendicular to the interface is that can be derived through dispersion relation, as shown in Equation (8).
The IFC shows that the multilayer is a type I hyperbolic metamaterial. Compared with the IFC estimated by effective permittivity tensor, the IFC calculated using the SFS equivalent model is precise. The ray vector derived from the IFC demonstrates that the imaging using the multilayer as a flat lens was based on the negative refraction instead of the nature of the negative index of refraction. Here an SFS arrangement was used to design a type I hyperbolic metamaterial in order to provide a precise iso-frequency curve. High transmission at normal incidence was necessary for a type I metamaterial. The thicknesses of the metal and dielectric films needed to be within the subwavelength scale. However, the type I iso-frequency curve of an alternatively arranged metaldielectric multilayer usually occurs at short wavelengths around the plasmonic resonate wavelength of the metal at which the permittivity of the metal film lead to low loss [20,21]. Besides the requirement of refractive index for type I isofrequency curve, the admittance should match to environment (air in this work) from normal incidence to an oblique angle to support the function of imaging. In order to overcome the frequency restriction for the occurrence of negative refraction, the type I hyperbolic metamaterial composed of Ag and Ta2O5 films was designed at a wavelength of 600 nm. According to our previous work [22], a thin metal film can be inserted between two metal films to perform a huge admittance locus and admittance matching to achieve high transmission property. The method, called modified Fabry-Perot (FP) design, can be extended to any odd-numbered metaldielectric film stack. At a wavelength of 600 nm, the silver film usually takes 10 nm to 20 nm to achieve the modified FP design. Therefore, each silver film in the film stack was set to have the same thickness of 20 nm for the five-layered symmetrical film stack MDMDM = Ag/Ta2O5/Ag/Ta2O5/Ag. The optical constants of Ag and Ta2O5 were adopted from the commercial optical thin-film software (Essential Macleod) [23]. respectively. The ranges of for low and low | | considerably overlaped, thus a multilayer with low loss could be derived by tuning the dielectric layer within both ranges. The as a function of and angle of incidence is shown in Figure 4c. The maximum is 2.156 at = 90 and =126 nm. At the range of between 67 nm and 74 nm, the versus showed that the SFS was equivalent to a dielectric-like anisotropic medium. The versus showed near type I hyperbolic shape at the range from = 100 nm to = 113 nm. The optimum transmission for the metamaterial in air was considered here. The transmittance of Air/SFS/Air was greater than from = 0 to = 60 . The thickness was tuned in the range that corresponded to type I hyperbolic metamaterial to yield a maximum of 47.56% at =100 nm. Figure 4d plots the isofrequency curve of Ag (20 nm)/Ta2O5 (100 nm)/Ag (20 nm)/Ta2O5 (100 nm)/Ag (20 nm), where was the wave vector in free space. Unlike a smooth and typical hyperbolic curve, the ratio / kept around 1.2 within the range of / between +1.56 and −1.56. It means that the reflected or diffracted ray vector was perpendicular to the interface for the parallel-to-interface component of the wave vector during the aforementioned range. The procedure for constructing the aforementioned modified FP SFS can be executed for a sevenlayered MDMDMDM for a wavelength of 600 nm with 20 nm-thick Ag films and equally thick Ta2O5 films. Figure 5a shows and as functions of . Figure 5b shows and as functions of . There are three ranges of corresponding to low ( < 0.1), which were : (64 nm, 67 nm), (83 nm, 97 nm), and (108 nm, 117 nm). The thickness of the dielectric layer in the seven-layered SFS with a low extinction coefficient was less than that of MDMDM. At = 113 nm, had a minimum magnitude of 0.1327 and the was 4.2736. varied from 1.28 to 4.02 and 2.037 to 9.053 at a low-magnitude (| | ≤ 0.3), which ranges were ( =84 nm, =92 nm) and ( = 110 nm, = 116 nm), respectively. The as a function of and angle of incidence is shown in Figure  5c. The type I hyperbolic metamaterial occurred at and ranged from 75 nm to 80 nm. The maximum of 46.49 % occurred at = 85 nm, and the corresponding iso-frequency curve is shown in Figure 5d. The optimum thickness for type I iso-frequency curve was shifted from = 100 nm in MDMDM structure to = 80 nm here in MDMDMDM structure. Compared with the five-layered case, the ratio / keeps around 1.1 within a smaller range of / between +1.7 and −1.7.

Conclusions
In conclusion, the equivalent model of a symmetrical film stack was applied to calculate the precise equivalent refractive index and admittance. In order to satisfy the criteria of branch selection, the correct branch solution of the refractive index was only valid for a certain range of wavelength or thickness of the dielectric layer. We clarify that the previous stratiform metamaterial is a type I hyperbolic metamaterial instead of a negative index medium. Based on the modified Fabry-Perot design for a metal-dielectric multilayer, a low loss hyperbolic metamaterial can be designed for an odd-numbered metal-dielectric symmetrical film stack. It is demonstrated that the required thickness of dielectric film for a type I metamaterial can be reduced by increasing the number of films in a symmetrical film stack. A metal-dielectric multilayer should be arranged as a multiple symmetrical film stack to exhibit the exact optical property. The method we proposed here can be applied to design a hyperbolic metamaterial in a precise and more flexible way.