Issue 
EPJ Appl. Metamat.
Volume 7, 2020



Article Number  1  
Number of page(s)  14  
DOI  https://doi.org/10.1051/epjam/2020001  
Published online  01 September 2020 
https://doi.org/10.1051/epjam/2020001
Research Article
Features of transmission of electromagnetic waves through composite nanoresonators including epsilonnearzero metamaterials
Gomel State Technical University, October ave. 48, Gomel 246746, Belarus
^{*} email: starodub@tut.by
Received:
13
December
2019
Accepted:
13
April
2020
Published online: 1 September 2020
Transmission of electromagnetic waves through nanometric multilayers (nanoresonators) including a main composite layer made of two alternating strips of lowabsorbing dielectrics that is sandwiched between epsilonnearzero (ENZ) or metallic spacer layers has been modeled. Analytical models are based on exact solutions of electromagnetic boundary problems. The spacers with the definite properties lead to extreme dependences of amplitude transmission coefficients on the system parameters and drastic increase in phase difference of the transmitted waves. These effects are most pronounced for subwavelength multilayer thicknesses due to multibeam interference features in the nanoresonator, and they can be amplified when the main layer and (or) the whole system thicknesses decrease. The investigated transmission features take place under variations of the system parameters such as anisotropy of the main layer materials, nonideal realization of ENZ materials, oblique incidence of the exciting radiation (for small incidence angles). The obtained results can have applications in development of ultrathin nanophotonics devices using phase transformation of transmitted waves.
Key words: Epsilonnearzero metamaterial / nanoresonator / transmission / phase difference / multibeam interference
© E. Starodubtsev, published by EDP Sciences, 2020
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
The progress in theory and technologies of nanophotonics, plasmonics, optics of metamaterials (MMs) and metasurfaces provided possibilities to reduce dimensions of many optical elements such as modulators, filters, absorbers, phase shifters, holograms, nanoresonators and others down to nanometer values (see, e.g. recent reviews [1–12] and numerous references in these works). The main purpose of such ultrasmall devices is to control the radiation characteristics (amplitude, phase, polarization) at deeply subwavelength dimensions of a composite material or system. In particular, large phase differences of two or more waves at the optical system output are required for many applications. Traditionally, the approaches based on increasing optical path difference of the waves (usage of proper waves in anisotropic media, propagation through different isotropic materials, application of optically dense materials, etc.) are used in this case. However, these methods require, as a rule, increasing geometric dimensions of optical elements (to several wavelengths and larger), and application of high optical density materials can lead to an undesirable high reflection or small transmission due to the large mismatch of electromagnetic properties of the device components and ambient media.
In the recent papers [13,14], new possibilities to develop ultrathin transformers of reflected waves phase differences and binary holograms on the basis of perforated multilayers including epsilonnearzero (ENZ) MMs were revealed and investigated analytically. ENZ MMs are characterized by small values of the real and imaginary parts of permittivity, and these materials are promising for many optics applications [15–21]. In particular, ENZ components allowed implementation of new operation regimes of nanoresonators: alloptical dynamic control of transmittance [22], hybridization of ENZ modes in double nanocavities [23], engineering of ENZ resonances in metalinsulatormetal nanocavities [24]. According to the results [13,14], the usage of ENZ substrates for such multilayers can lead to the extreme values of amplitude reflection coefficients depending on the system parameters. The ENZ (or metallic for thicker systems) substrate has a significant impact on the transformation of phase difference of the reflected waves. According to the modeling results, the optimization of this system properties allows decreasing the holograms thicknesses down to the values of the order of 10–20 nm as well as obtaining nanometric thickness phase transformers of reflected waves. The analysis [13,14] pointed to the link between the features of interaction of electromagnetic waves with ENZ layers (for the reflection problems) and realization of regimes of effective phase transformation of the radiation using such lowdimensional systems. These results make it possible to consider the more general problem of application of ENZ MM multilayer components for phase transformation of not only reflected but also transmitted radiation.
The present paper is devoted to analytical modeling (using exact solutions of the corresponding boundary problems) of transmission of electromagnetic waves through composite multilayers of subwavelength thicknesses that include a main dielectric layer sandwiched between two ENZ MM (or metallic) spacer layers. The numerical and graphical analysis of such a nanoresonator model has the following main goals: (i) to specify the features of operation of the multilayers including ENZ MM spacers that determine new possibilities of obtaining ultrathin devices with the large phase difference and transmission of output waves; (ii) to optimize parameters of the system under consideration for typical materials of the layers. The obtained results confirm, in general, possibilities to apply such nanoresonators as effective ultracompact phase transformers of a new type that can be useful for nanophotonics applications.
The work has the following structure. The statement of the problem, solutions of electromagnetic boundary problems and some consequences of these solutions are given in Section 2. The numerical modeling and graphical analysis of the results (Sect. 3) explain the main features of amplification of transmitted waves phase difference using the ENZ or metallic spacers of the nanoresonators. In particular, the conditions are considered in detail when decreasing the main layer and (or) the whole system thickness down to significantly subwavelength values leads to increasing phase difference of the output waves and rather high transmission. Section 4 includes a brief discussion of the results and summarizes the paper. The additional calculations excluded from the main text are contained in Appendix.
2 Statement of the problem and analytical model
Let us investigate transmission of an electromagnetic plane monochromatic wave through a layered system under oblique incidence of the wave (Fig. 1). The multilayer scheme and geometry of the corresponding boundary problem are illustrated in Figure 1a and b. Layers of thicknesses d_{l} (here and below indexes α, l = 1,2,3, β correspond to the quantities characterizing the multilayer components in Fig. 1b) are, respectively, the upper, main, and bottom layers. The upper and bottom layers can play a role of spacers, protective or controlling layers, nanoresonator mirrors depending on applications of the system. Media in front of the upper layer (e.g. air) and behind the bottom layer (the substrate) are assumed to be semiinfinite (optically thick). So, the reflected wave inside of the substrate is absent. Coordinate axis Z is perpendicular to the layers boundaries and plane XZ is the incidence plane. Phase multiplier exp [ik_{0} (m_{1}x + m_{3}z) − iωt] is used in expressions for the fields below where m = (m_{1}, 0, m_{3}) = k/k_{0} is a complex refraction vector [25], k is a wave vector (vectors are denoted with boldface symbols), k_{0} = ω/c is a wave number for vacuum, i^{2} = − 1, r = (x, y, z) is a radiusvector. The transversal component of refraction vector is supposed to be real (m_{1} is constant at boundaries z = 0, d_{l}) and the longitudinal component m_{3} is complex in the general case. Electric field strengths of the incident, E_{in}, reflected, E_{r1,2}, and transmitted, E_{t1,2}, waves are considered at the boundaries z = 0 and z = z_{3} = d_{1} + d_{2} + d_{3} (Fig. 1).
Isotropic materials of the input nonabsorbing medium (superstrate), absorbing layers 1, 3 (spacers), and substrate are characterized by scalar dielectric permittivities ε_{α,1,3,β}, correspondingly. The main composite layer is formed from the alternating “strips” of two different materials that are assumed to be optically uniaxial in the general case (Fig. 1). The Strips 1 and 2 materials are characterized by the permittivity tensors and , where is double singular and is nonsingular complex eigenvalues of tensor (here and below the upper indexes (1,2) point to the quantities corresponding to Strips 1, 2). So, the optical axis (OA, Fig. 1b) for the Strip 1 (2) material is parallel to axis X (Y), respectively.
All the multilayer components are also assumed to be absorbing (besides the superstrate) and nonmagnetic. Permittivities ε_{α,1,3,β} or components of tensors can have the negative real parts, all the imaginary parts of these quantities are positive according to the causality condition. So, we consider absorbing linear materials of the layers and substrate including the possible case of the MM (or metallic) layers. The denotations a ′ = Re (a) , a ″ = Im (a) will be also used for the real and imaginary parts of scalar or vector quantities.
Let us also assume that the subwavelength thickness of the strips d_{2} (along axis Z) is much less than their linear sizes in plane XY: d_{2} ≪ l_{x} ≪ l_{y} where l_{x,y} are characteristic dimensions of the strips along axes X, Y (the value of d_{2} is significantly increased in Fig. 1b for illustrative purposes). The strips can be obtained with various methods and materials, e.g. they can be formed from liquid crystals with the mutually perpendicular director orientations. The similar schema of application of liquid crystals cells are used in components of display technologies for phase modulation of the transmitted radiation [26–28].
One can note that the considered geometry of inclusions (the strips) of two materials in the main composite layer is not critical for the further analysis. For example, this layer can be originated by the strips with different (unequal) values of parameters l_{x,y} or by the system of parallel to axis Z microholes in a homogeneous host material [13,14,29]. It is important that there are, mainly, two transmitted waves at the output of the system (in the substrate), and these waves correspond to interaction of the electromagnetic field with two different materials (Strips 1, 2) in the main layer.
We use the constitutive relations D = ε E, B = μ H (μ = 1) and Maxwell's equations for monochromatic waves [25](1)where ε ^{−1} is the inverse to ε tensor, and the dispersion equation, m ^{2} = ε, for isotropic materials of layers 1, 3 and media outside of the multilayer. Then the electric field strength of the waves in media v = α, 1, 3, β (Fig. 1b) can be written as(2) where (and below) upper indexes “+” and “–” correspond to the quantities characterizing transmitted and reflected proper waves in an isotropic medium (that decay in the directions along and oppositely to axis Z, respectively), and values of are chosen from the corresponding dispersion equations for the transmitted waves, according to the condition In equation (2), the common for all the waves factor exp(ik_{0}m_{1}x − iωt) is omitted and complex scalar coefficients are determined from the solutions of the boundary problems (for the substrate ). The magnetic fields strengths corresponding to equation (2) are determined by the second relation of equation (1). Coefficients in equation (2) are different for the waves corresponding to interaction of the incident wave with the Strips 1 and 2 materials (e.g. for fields E_{r1} and E_{r2}, E_{t1} and E_{t2}, Fig. 1b).
The dispersion equations for uniaxial materials of the strips and corresponding TM and TE proper waves in these materials take the form(3) (4) (5) (6)where equations (3), (4) and (5), (6) correspond to the materials of Strips 1 and 2, the values of m_{3} are determined by the corresponding dispersion equations according to the condition for the transmitted (reflected) waves, and the phase multiplayer exp [ik_{0} (mr) − iωt] is omitted for each of the proper waves. Complex scalar coefficients P, Q, T, R in equations (3)–(6) for the transmitted and reflected waves are determined from the solutions of the boundary problems.
To solve the electromagnetic boundary problems, we apply the standard approach using continuity of tangential components of the resulting fields E, H on the multilayer boundaries. Then the following expressions for the amplitude transmission coefficients of wave 1 (transmitted through Strip 1, Fig. 1) for TM and TE polarizations of the incident wave, can be obtained (7) (8)
where the denotations for equations (2)–(6) are also used, (according to the dispersion equations from Eqs. (3) and (4)).
In a similar manner, we have the expressions for the amplitude transmission coefficients of wave 2 (transmitted through Strip 2, Fig. 1) for TM and TE polarizations of the incident wave, (9) (10)
where (according to the dispersion equations from Eqs. (5), (6)). In equations (7)–(10), the expressions for quantities A_{±}, B_{±}^{ (1,2)}, C_{±}^{ (1,2)}, D_{±}, E_{±}^{ (1,2)}, F_{±}^{ (1,2)} are given in Appendix, values of η_{α,β,1,3}, η ^{ TM,TE } are chosen for the transmitted waves (Im (η_{α,β,1,3}) > 0, Im (η ^{ TM,TE }) > 0), and the parameters are also used: a = ε_{α}/η_{α}, b = ε_{1}/η_{1}, , d = ε_{3}/η_{3}, f = ε_{β}/η_{β}, g = η ^{ TE (1)}/η_{1}, h = η_{β}/η_{3}, , q = η ^{ TE (2)}/η_{1}.
Quantities , can be considered as the energy transmission coefficients at the boundary z = z_{3}, and , are the phase shifts arising in transmission of the waves corresponding to rays 1 and 2 (Fig. 1). So, the exact solutions of the boundary problems are determined by equations (7)–(10). Let us note that quantities in equations (7)–(10) characterize relations of the corresponding components of field E.
The phase shift difference (PSD) arising for transmitted (TM or TE) waves 1 and 2 at the boundary z = z_{3} and determining phase transformation effectiveness takes the form(11)
Besides the exact PSD (Eq. (11)) at the multilayer output, we will also use the phase differences arising between TM or TE waves (that is, for the TM or TE exciting wave) passing one time (for one pass): (i) only through the different strips of the main layer(12)(ii) only through the different strips of the main layer, but for the case when the thickness of these strips (d_{2}) is equal to the whole multilayer thickness (d_{1} + d_{2} + d_{3}), that is, when replacing layers 1, 2, 3 by only the Strips 1, 2 materials(13)
Then relations Δϕ/ΔΦ, Δϕ/ΔΦ^{∼} for the corresponding polarizations characterize the relative changes of waves 1, 2 PSD at the multilayer output (in presence of layers 1, 2, 3) in comparison with a “one pass” PSD of these waves arising for passing through: (i) only the main layer, equation (12); (ii) the whole multilayer filled by only materials of the strips, equation (13).
Using the obtained transmission characteristics (Eqs. (7)–(10)) and the boundary problems solutions, one can also determine quantities that serve as peculiar amplitude reflection coefficients characterizing the radiation parameters inside of the main layer relatively to ones for the incident wave. In particular, we have for the field inside of Strip 1(14) (15) (16) (17) where the denotations from equations (3), (4), (7), (8) are used. Quantities and characterize reflection of the proper waves of the corresponding polarizations from the upper (z = z_{1} = d_{1}) and bottom (z = z_{2} = d_{1} + d_{2}) boundaries of the main layer (Fig. 1b) inside of the layer. As above, we also use energy reflection coefficients at these boundaries that are determined by the relations , . In a strict sense, quantities and are not reflection coefficients, as the numerators and denominators in equations (14)–(17) (e.g. ) are taken for the different boundaries z = 0, z_{1,2} without regard for the corresponding phase multipliers. However, these quantities can be used to characterize qualitatively energy distribution of the proper waves in the main layer near boundaries z = z_{1,2}. The similar to equations (14)–(17) relations can be obtained for the field inside of the Strip 2 material.
Using equations (3), (14), (15) and the expression for the timeaveraged Pointing vector, , where c is the speed of light in vacuum and asterisk (here and below) denotes complex conjugation, one can obtain the normalized value of the resulting electromagnetic field intensity inside of the main layer (for Strip 1 and TM incident wave) (18)
where the denotations are used: ,, S_{in} (z = 0) is the timeaveraged Pointing vector and ψ is incidence angle of the incident wave on the boundary z = 0. Equation (18) takes into account multibeam interference of the opposite transmitted and reflected TM proper waves inside of the main layer.
The standard relations for phase incursion of the waves inside of the resonator are used usually to determine conditions when transmission maxima of the similar nanoresonators are realized (e.g. [22,30]). However, obtaining compact analytical expressions for phases of the transmitted waves is not trivial in the case under consideration. So, the detailed numerical analysis of phases, amplitude and energy transmission coefficients of the output waves will be executed below using equations (7)–(13).
Obviously, the obtained simplified model considers only two “ways” of the transmission in the system that correspond to “conventional” transmission of rays 1 and 2 (Fig. 1). Clearly, a part of the radiation also reflects from the strips vertical walls and propagates from one to other strip with the incidence angle ψ growth. Moreover, the analytical model does not account for the strips form, their location in plane XY, and possible excitation of additional nearfield modes including the effects of extraordinary interaction of radiation with the system (under conditions of surface plasmons excitation for the ordered inclusions arrays) [31–34]. However (similarly to the reflection models [13,14]), the multilayer (Fig. 1) includes dielectric strips with the deeply subwavelength thickness (d_{2} ≪ λ) and “macroscopic” dimensions in plane XY (λ < l_{x} ≪ l_{y}). For optical applications, the values of l_{x,y} (or other parameters characterizing sizes of the main layer inhomogeneities and distances between them in plane XY) can be estimated as several microns or larger. Exactly that and the condition of small incidence angles allow the use of transmission coefficients (Eqs. (7)–(10)) for the main part of the transmitted radiation, and without application of more complicated diffraction theory methods (e.g. [35–38]). The similar approach (based on exact analytical solutions of the boundary problems) was recently confirmed experimentally for calculation of the waves reflected from the nanometric binary holograms [29].
Fig. 1 (a) Scheme of the multilayer with the main composite layer formed from alternating “strips” of two different materials that can be optically uniaxial in the general case. The rays in the substrate transmitted through the multilayer are denoted by numbers 1 and 2. E_{in} is the incident wave field, the reflected rays in the superstrate are not shown. (b) Geometry of the corresponding electromagnetic boundary problems. The multilayer includes upper and bottom layers, which electromagnetic properties can correspond to isotropic ENZ MMs or metals, and the main layer. The superstrate and substrate materials are isotropic and optically thick. E_{r1}, E_{r2}, and E_{t1}, E_{t2} are the fields of the waves reflected from and transmitted through Strips 1 and 2, respectively, when the incident wave interacts with the multilayer. OA denotes the optical axis direction in the materials of Strips 1 or 2. The other denotations are given in the main text. 
3 Numerical and graphical analysis
3.1 Main assumptions and parameters
The main goal of the analysis below is to determine parameters of the investigated systems when rather large values of PSD and transmission of waves 1, 2 take place for the minimal main layer or whole multilayer thicknesses. According to equations (11)–(13), it means that we consider the cases when Δϕ/ΔΦ > 1, Δϕ/ΔΦ^{∼} > 1, that is, when it is possible to “amplify” PSD of waves 1, 2 in comparison with the “one pass” PSD corresponding to the onetime run of the radiation through the main layer or the whole system. Indexes TM, TE are further omitted when both types of the waves are considered or the choice of the waves is clear from the context.
We choose some model ENZ MMs and silver as a metallic component frequently used in photonics applications for the materials of layers 1, 3 (the upper and bottom layers, Fig. 1). The cases are considered when the main layer is formed from the strips of transparent (lowabsorbing) dielectrics (the numerical analysis shows that it is a rather optimal case). For all the figures below, air is assumed to be the superstrate and substrate medium (ε_{α, β} = 1), and the spacer layers thicknesses are equal (d_{1} = d_{3}). The effects under consideration also take place for the general (less “symmetric”) cases when ε_{α} ≠ ε_{β}, d_{1} ≠ d_{3}.
The main features of the nanoresonators occur for both oblique (for small incidence angles) and normal incidence of the exciting wave (when the difference between TE and TM waves disappears). Isotropic dielectrics are chosen as the main layer materials for many figures below (Figs. 2–10, Sect. 3.2). The impact of optical anisotropy of the strips materials (that is not a determining factor for the effects under investigation within the considered parameter ranges) is characterized by Figures 11, 12 data and numerical estimations given below (Sect. 3.3). For all the graphs in Figures 2–10, the parameters are chosen: ψ = 0 (besides the Fig. 8 data), , (the Strip 2 material is air). The other parameter values are given in the figures captions.
Phase shifts in the graphs further are calculated in relative units of radians/π. Phases ϕ_{1,2} are determined below using builtup functions of the computer algebra system (Mathcad, PTC) and values of t_{1,2} (Eqs. (7)–(10)). As these builtup functions return the principal argument of a complex number (between –π and π), so phase values “jumps” can be in graphs of dependences of quantities ϕ_{1,2} on the system parameters when phase takes values ±π.
3.2 Case of isotropic main layer materials
The data in Figure 2 illustrate phase and energy characteristics of the transmitted waves depending on the main layer thickness for the following choice of the parameters of layers 1, 3: “conventional” layers (Fig. 2a,d); one of the layers is conventional, and the other one is ENZ (Fig. 2b,e); both the layers are ENZ (Fig. 2c,f). When thickness d_{2} is varying, the phases and transmission coefficients modulation is significantly amplified in the transition from “conventional” to ENZ layers 1, 3. Figure 2b, e data are “intermediate“ in comparison with the cases of “conventional” and ENZ layers 1, 3, and the corresponding graphs do not change when the layers properties are rearranged: ε_{1} ↔ ε_{3}. When the values of d_{2} change, for the case of ENZ layers 1, 3, the additional modulation of quantities ϕ_{1,2} and periodic resonance changes of quantities T_{1,2} from 0 up to 0.9 are realized (the data in Fig. 2c,f in comparison with ones in Fig. 2a,d). It is seen that positions of characteristic points in the graphs (phase jumps, transmission extrema) depend strongly on the layers 1, 3 material. The distances between the neighboring abrupt changes of ϕ_{1,2} values from π to –π (Fig. 2a–c) and maxima of functions T_{1,2} (d_{2}) (Fig. 2d–f) are well described by constant values Δd_{2} and Δd_{2}/2, correspondingly, where . So, these properties for each of the nanoresonators corresponding to Strips 1 and 2 (for lowabsorbing main layer materials) are similar to ones for the “ideal” FabryPerot resonator [39] and have the pronounced interference nature. One can see that the presence of both ENZ layers 1, 3 leads to the significant amplification of the interference effects.
Changes of quantities in the transition from “conventional” to ENZ layers 1, 3 (for the data of Fig. 2a,d and c,f) are illustrated in Figure 3. For “conventional” layers 1, 3, the graphs of functions are rather smooth and periodic (Fig. 3a). It is seen, that the resonance character of the dependences in Figure 2c, f is due to the similar features of graphs (Fig. 3b). For the case of ENZ layers 1, 3, the trajectories of the tips of vectors on the complex plane change strongly and similarly when thickness d_{2} is varied (Fig. 3c–f). For quantities t_{1,2} characterizing waves 1, 2, the trajectory form changes cardinally (Fig. 3c,e in comparison with Fig. 3d,f). The data in Figure 3g,h illustrate the trajectories of vectors t_{1,2} tips for changes of parameters in the range (− 2 ;2) when d_{2} = 39 nm. For the values of the order the tip of vector t_{1} “moves” on the complex plane much “faster” than one for vector t_{2} providing larger values of Δϕ for ENZ layers 1, 3. However, the trajectories of the tips of vectors t_{1,2} for changing are similar (Fig. 3g,h). In Figure 3c–h graphs, there are denotations of the characteristic points (values of d_{2}, Fig. 3c–f, and parameters , Fig. 3g,h) describing the “motion” of vectors t_{1,2} tips for increasing d_{2} or .
Figure 4 illustrates phase differences (Δϕ, ΔΦ, ΔΦ^{∼}) and transmission coefficients (T_{1,2}) of waves 1, 2 for “conventional” and ENZ layers 1, 3 depending on the system parameters. For “conventional” layers 1, 3, quantities T_{1,2} depend weakly on parameters d_{2} and λ (when neglecting frequency dispersion of the multilayer components), and quantities Δϕ, ΔΦ have close values (Fig. 4a,b). So, in this case PSD of waves 1 and 2 at the boundary z = z_{3} coincides practically with “one pass” phase difference ΔΦ (Eq. (12)). The situation changes drastically when choosing ENZ layers 1, 3 for rather thin (subwavelength) main layers (Fig. 4c–h). With that, the ranges of parameters d_{2}, λ, d_{1,3}, values can be realized when the condition Δϕ ≫ ΔΦ is satisfied. So, ENZ layers 1, 3 (isotropic and the same ones for Strips 1, 2) lead to the drastic increase of waves 1, 2 PSD at the multilayer output in comparison with quantity ΔΦ. In this case, the strong dependences of quantities T_{1,2} on the system parameters also take place (Fig. 4c,d,f,g). It is important that large values of Δϕ and Δϕ/ΔΦ can be obtained for the rather high transmission. For example, the first crossing of the functions T_{1,2} (d_{2}) graphs in Figure 4c (for increasing d_{2}) corresponds to the main layer thickness d_{2} = 39 nm when Δϕ ≈ π/2, Δϕ/ΔΦ ≈ 7 (Fig. 4e), T_{1,2} ≈ 0.43. The corresponding point is also pointed in Figure 3c–f, and this value of d_{2} is used for the data in Figure 3g, h. The graphs of functions Δϕ/ΔΦ (d_{2}), Δϕ/ΔΦ (d_{1,3}), have the pronounced maxima near values d_{2} = 36 nm, d_{1,3} = 330 nm, , respectively (Fig. 4e–h). Let us note that condition Δϕ/ΔΦ^{∼} < 1 takes place for all the data in Figure 4 (that characterize the case of a rather large dielectric contrast between isotropic materials of the strips: ). So, the fulfilment of conditions ΔΦ ≪ Δϕ < ΔΦ^{∼} is possible for the considered parameters (for rather small values of d_{2} and optimally selected values of d_{1,3} for ENZ layers 1, 3). For these data and small dielectric contrasts (), the numerical analysis also showed the simultaneous fulfilment of more optimal conditions of PSD amplification: Δϕ/ΔΦ > 15, 1 < Δϕ/ΔΦ^{∼} < 2.
The considered effects are also realized for larger values of parameter (e.g. for ). However, in this case the satisfaction of condition Δϕ ≫ ΔΦ worsens (as the straight line ΔΦ(d_{2}) slope increases), and it is difficult to obtain simultaneously large values T_{1,2} (as the maxima of functions T_{1} (d_{2}) and T_{2} (d_{2}) move apart with the growth of parameter ). Increasing parameter is also not optimal, as it leads to decreasing waves 1, 2 transmission.
The graphs in Figure 5 characterize the multilayer corresponding to Figure 4 data but for “metallic” layers 1, 3 and rather optimal (for realization of large PSD and transmission of waves 1, 2) parameters d_{1,2,3}. Qualitatively, the graphs in Figure 5a,b are similar to ones in Figure 4c,e. In particular, with the growth of values d_{2}, the first crossing of the functions T_{1,2} (d_{2}) graphs in Figure 5a realizes for d_{2} = 156 nm, when Δϕ ≈ 0.7π, Δϕ/ΔΦ ≈ 2, Δϕ/ΔΦ^{∼} ≈ 2 (Fig. 5b), T_{1,2} ≈ 0.26. However, functions (Δϕ/ΔΦ (d_{1,3})) have no maxima and decrease (increase) monotonically for the considered parameters values (Fig. 5c,d). That explains the choice of rather small (negative) values of and thicknesses d_{1,3} for realization of large values of Δϕ, T_{1,2}. With that, we have the condition d_{2} ≪ d_{1,3} for ENZ layers 1, 3 (Fig. 4), and the opposite one d_{2} ≫ d_{1,3} for “metallic” layers 1, 3 (Fig. 5).
The transmission characteristics of waves 1, 2 depending on parameter d_{2} for the case of “metallic” layers 1, 3 are illustrated in Figure 6. The comparison of the graphs in Figures 2, 6 shows that functions ϕ_{1,2} (d_{2}), T_{1,2} (d_{2}) have a number of similar features for one or two layers 1, 3 of ENZ or “metalliclike” materials (comparing Fig. 2b, e, f and Fig. 6a, c, d). This leads to the qualitative features considered above for the data in Figure 5. As for the data in Figure 2b, the graphs in Figure 6a,c do not change for the layers properties rearrangement: ε_{1} ↔ ε_{3}.
The graphical analysis above illustrates two main “regimes” of operation of the multilayers depending on layers 1 and 3 properties: for ENZ and metallic materials of the layers. These regimes can be realized for both the different and the same material of layers 1, 3 (assuming these layers to be similar). In particular, one material for two regimes can be used for the different wavelengths of the exciting radiation taking into account frequency dispersion of the dielectric properties. The data in Figure 7 illustrate possibility to use silver as the layers 1, 3 material. The literary data contain a rather wide variation of the Ag permittivity values subject to fabrication methods, thicknesses of the layers, etc. For the numerical analysis, the following values were used: ε_{Ag} = 2.27 × 10^{−3} + 0.08i, − 19.74 + 0.85i for λ = 292, 650 nm, correspondingly [40]. The graphs in Figure 7a, b, c (the ENZ regime) and Figure 7d, e, f (the metallic regime) well correspond to ones in Figure 4c, e, g and Figure 5a, b, d, respectively. In particular, for increasing values of d_{2}, the first crossing of the graphs T_{1,2} (d_{2}) in Figure 7a and d corresponds to the following values of parameters: d_{2} = 19 nm, Δϕ ≈ 0.5π, Δϕ/ΔΦ ≈ 6, Δϕ/ΔΦ^{∼} ≈ 0.5, T_{1,2} ≈ 0.2 and d_{2} = 171 nm, Δϕ ≈ 0.6π, Δϕ/ΔΦ ≈ 2, Δϕ/ΔΦ^{∼} ≈ 2, T_{1,2} ≈ 0.3. So, silver layers 1, 3 can provide the realization of the system under investigation for two considered operation regimes.
Dependences of the main transmission characteristics of the multilayers on the incidence angle (under ψ = 0 ÷ π/10 and TE, TM polarizations of the incident wave) for the parameters corresponding to Figure 7 data are illustrated in Figure 8. The graphs in Figure 8a–d and e–h correspond to the ENZ and metallic regimes, respectively, for silver as the layers 1, 3 material. The thicknesses d_{1,2,3} for the data in Figure 8 correspond to the first crossing of the graphs T_{1,2} (d_{2}) in Figure 7a and d, when rather large values of parameters Δϕ, Δϕ/ΔΦ, Δϕ/ΔΦ^{∼}, T_{1,2} take place simultaneously. With growth of the incidence angle, phase characteristics Δϕ, Δϕ/ΔΦ have larger values and change (decrease) stronger in the case of ENZ layers 1, 3 and TM polarization of the incident wave (the left panels of Fig. 8 in comparison with the right ones). Parameters ΔΦ and Δϕ/ΔΦ^{∼} change rather weakly for all the data in Figure 8. For angles ψ = 0 ÷ π/20, quantities T_{1,2} have the comparable values for both ENZ (Fig. 8a,c) and metallic (Fig. 8e,g) layers 1, 3. Thus, Figure 8 data characterize the possibilities of simultaneous realization of rather large values of PSD and transmission for both the normal and oblique incidence of the incident radiation on the ultrathin multilayers (for rather small ψ values).
According to the similar analysis (the data are not shown), Figures 2–5 data do not change qualitatively for cases of anisotropic materials of the strips (for rather small anisotropy values of the strips material, e.g. for parameters ; ). These data changes are also very small for oblique incidence and TM or TE polarizations of the exciting wave up to incidence angles of the order of ψ = π/20 (the changes do not exceed several per cents, especially for TE waves and metallic layers 1, 3).
Electromagnetic field inside of the main layer for various parameters of the multilayer is characterized by the graphs in Figures 9 and 10. It is seen that the transition from “conventional” (Fig. 9a,c) to ENZ (Fig. 9b) or metallic (Fig. 9d) layers 1, 3 is accompanied by the significant amplification of the reflection (described by parameters R_{21}, R_{23}) of the waves on the boundaries z = z_{1,2} inside of the layer for the definite thicknesses d_{2}. With that, the values of d_{2} for the “peaks” of values of R_{21}, R_{23} and T_{1} for the corresponding graphs in Figure 9b and Figure 2f, Figure 9d and Figure 6d coincide. Near to the “peaks” of values of R_{21}, R_{23} and T_{1} (that is, when the optical interference effects are amplified in the main layer), the leapfrog change of the wave 1 phase takes place (Fig. 9b,d). As for Figure 2 data, the distances between the neighboring abrupt changes of ϕ_{1} values from π to –π and maxima of functions R_{21} (d_{2}), R_{23} (d_{2}) (Fig. 9b,d) are well described by values Δd_{2} and Δd_{2}/2, respectively. The values R_{21}, R_{23} > 1 in Figures 9, 10 data characterize qualitatively reflection from the corresponding boundaries without taking into account the opposite waves interference inside of the main layer. For all the data in Figures 9 and 10, the normalized intensity of the resulting field inside of the main layer for Strip 1 (determining by Eq. (18) with account of the opposite waves interference) has the values I ^{ TM (1)} < 1 (values of I ^{ TM (1)} (z = d_{1} + 0.5d_{2}) are given in Figs. 9, 10). This condition points to the conservation of the energy fluxes balance in the main layer. So, Figure 9 data confirm the interference nature of the considered effects of PSD amplification for waves 1, 2 at the multilayer (nanoresonator) output.
The graphs in Figure 10 illustrate the effect of parameters , d_{1,3} on the energy fluxes and field intensity inside of the main layer. For ENZ layers 1, 3, the pronounced maxima of the functions characterizing dependences of quantities R_{21}, R_{23} on parameters , d_{1,3} (Fig. 10a,b) take place. For metallic layers 1, 3, such maxima are absent (the data in Fig. 10с for , Fig. 10d). When varying parameters d_{1,3}, the abrupt phase ϕ_{1} changes take place for ENZ layers 1, 3 (Fig. 10b). For metallic layers 1, 3 (Fig. 10d), the values of R_{21},R_{23}, I ^{ TM (1)} decrease quickly with growth of d_{1,3}. The graphs of quantity I ^{ TM (1)} given in Figures 9, 10 depend very weakly on the choice of coordinate z of the point inside of the main layer. That is conditioned by the small value of parameter determining small absorption of the radiation inside of the main layer. According to the numerical analysis results, dependences I ^{ (1)} (z) enhance with the growth of absorption (for values of the order of 10^{−2} ÷ 10^{−1}). With that, values of I ^{ (1)} decrease with increasing z from d_{1} to d_{1} + d_{2}.
Fig. 2 Dependences of phases ϕ_{1,2} (a–c) and energy transmission coefficients T_{1,2} (d–f) of waves 1, 2 on the main layer thickness d_{2} for various materials of layers 1, 3. The values of parameters are: ε_{1} = ε_{3} = 2 + 10^{−3} i (a,d); ε_{1} = 2 + 10^{−3} i, ε_{3} = 0.01 + 0.01i (b,e); ε_{1} = ε_{3} = 0.01 + 0.01i (c,f). The values of parameters for all the graphs are: d_{1,3} = 0.25 μm, λ = 0.63 μm. Parameter d_{2} is used as the abscissa axis for all the graphs. 
Fig. 3 The effect of parameter d_{2} on the real and imaginary parts of amplitude transmission coefficients of waves 1, 2 (a, b) and trajectories of the tips of vectors t_{1,2} on the complex plane depending on parameters d_{2} (c–f) and (g,h). The values of parameters are: ε_{1} = ε_{3} = 2 + 10^{−3} i (a, c, e); ε_{1} = ε_{3} = 0.01 + 0.01i (b, d, f); , d_{2} = 39 nm (g, h). The values of parameters for all the graphs are: d_{1,3} = 0.25 μm, λ = 0.63 μm. 
Fig. 4 Dependences of PSD Δϕ, ΔΦ, relations Δϕ/ΔΦ, Δϕ/ΔΦ^{∼} and energy transmission coefficients T_{1,2} of waves 1, 2 on the main layer thickness d_{2} (a, c, e), wavelength λ (b, d), parameters (f, h), d_{1,3} (g) for various materials of layers 1, 3. The values of parameters are: ε_{1} = ε_{3} = 2 + 10^{−3} i (a, b); ε_{1} = ε_{3} = 0.01 + 0.01i (c, d, e, g); d_{2} = 39 nm (b, d, f, h); (f, h). The values of parameters for all the graphs are: d_{1,3} = 0.25 μm (besides g); λ = 0.63 μm (besides b, d). Figures e and h detail the data in Figures c and f, respectively. The black horizontal dashed line in Figures e, g, h characterize equality to 1 of the corresponding quantities. 
Fig. 5 Dependences of PSD Δϕ, ΔΦ, relations Δϕ/ΔΦ, Δϕ/ΔΦ^{∼} and energy transmission coefficients T_{1,2} of waves 1, 2 on the main layer thickness d_{2} (a, b), parameters (c), d_{1,3} (d) for a metallic material of layers 1, 3. The values of parameters are: ε_{1} = ε_{3} = − 13 + 0.5i (a, b, d); d_{2} = 156 nm (c, d); (c). The values of parameters for all the graphs are: d_{1,3} = 13 nm (besides d), λ = 0.6 μm. Figure b details the data in Figure a. The black horizontal dashed line in fig. b characterizes equality to 1 of the corresponding quantities. 
Fig. 6 Dependences of phase differences ϕ_{1,2} (a, b) and energy transmission coefficients T_{1,2} (c, d) of waves 1, 2 on the main layer thickness d_{2} for various materials of layers 1, 3. The values of parameters are: ε_{1} = 2 + 10^{−3} i, ε_{3} = − 13 + 0.5i (a, c); ε_{1} = ε_{3} = − 13 + 0.5i (b, d). The values of parameters for all the graphs are: d_{1,3} = 13 nm, λ = 0.6 μm. Parameter d_{2} is used as the abscissa axis for all the graphs. 
Fig. 7 Dependences of PSD Δϕ, ΔΦ, relations Δϕ/ΔΦ, Δϕ/ΔΦ^{∼} and energy transmission coefficients T_{1,2} of waves 1, 2 on thicknesses d_{2} (a, b, d, e), d_{1,3} (c, f) for silver layers 1, 3 and two different wavelengths of the incident wave. The values of the Ag permittivity for λ = 292 nm (a, b, c), 650 nm (d, e, f) are given in the text. The values of parameters are: d_{1,3} = 125 nm (a, b), 9 nm (d, e); d_{2} = 19 nm (c), 171 nm (f). Figures b and e detail the data in figures a and d, respectively. The black horizontal dashed line in figures b, c, e characterizes equality to 1 of the corresponding quantities. 
Fig. 8 Dependences of PSD Δϕ, ΔΦ, relations Δϕ/ΔΦ, Δϕ/ΔΦ^{∼} and energy transmission coefficients T_{1,2} of waves 1, 2 on incidence angle ψ for silver layers 1, 3, two different wavelengths and TM (a, b, e, f), TE (c, d, g, h) polarizations of the incident wave. The values of parameters are: (a–d) λ = 292 nm, d_{1,3} = 125 nm, d_{2} = 19 nm; (e–h) λ = 650 nm, d_{1,3} = 9 nm, d_{2} = 171 nm. Figures b, d, f, h detail the data in figures a, c, e, g, respectively. The black horizontal dashed line in figures b, d characterizes equality to 1 of the corresponding quantities. Parameter ψ is used as the abscissa axis for all the graphs. Values of ψ (as Δϕ, ΔΦ, ΔΦ^{∼}) are given in the relative units of radians/π. 
Fig. 9 Dependences of quantities ϕ_{1}, R_{21}, R_{23}, I ^{ (1)} on the main layer thickness d_{2} for various materials of layers 1, 3. The values of parameters are: ε_{1} = ε_{3} = 2 + 10^{−3} i (a, c), 0.01 + 0.01i (b), −13 + 0.5i (d); d_{1,3} = 250 nm (a, b), 13 nm (c, d); λ = 0.63 μm (a, b), 0.6 μm (c, d). Denotations of the curves in Figures b, d are given in Figures a, c. I ^{ (1)} = I ^{ TM (1)} (z = d_{1} + 0.5d_{2}) for all the graphs. Parameter d_{2} is used as the abscissa axis for all the graphs. 
Fig. 10 Dependences of quantities ϕ_{1}, R_{21}, R_{23}, I ^{ (1)} on parameters (a, c), d_{1,3} (b, d) for various materials of layers 1, 3. The values of parameters are: (a, b), 0.5 (c, d); (b), −13 (d); d_{2} = 39 nm (a, b), 156 nm (c, d); d_{1,3} = 250 nm (a, b), 13 nm (c, d); λ = 0.63 μm (a, b), 0.6 μm (c, d). I ^{ (1)} = I ^{ TM (1)} (z = d_{1} + 0.5d_{2}) for all the graphs. The black horizontal dashed line in figs. a, b, c characterizes equality to 1 of the corresponding quantities. 
3.3 Case of uniaxial main layer materials
The data in Figure 11 correspond to ones in Figure 4c,e,f,h but for the case of an uniaxial material (with the different orientations of the optical axis for the strips, Fig. 1) of Strips 1, 2. It is seen that the main features of the considered effects also take place for the anisotropic components of the main layer. For ENZ layers 1, 3, there are pronounced maxima of the functions characterizing dependences of quantities Δϕ, Δϕ/ΔΦ, Δϕ/ΔΦ^{∼}, T_{1,2} on parameters d_{2} and (Fig. 11). With that, conditions Δϕ/ΔΦ ≫ 1, Δϕ/ΔΦ^{∼} > 1 take place (Fig. 11b, d). Small anisotropy values of the strips material (ε_{o} − ε_{e} <0.5, as for the case of isotropic strips with small values ) lead to the growth of possible values Δϕ/ΔΦ, Δϕ/ΔΦ^{∼} (the data in Fig. 11b,d in comparison with ones in Fig. 4e,h). According to Figure 11 data, the parameters ranges are also possible for anisotropic materials of the strips when rather large transmission and amplification of waves 1, 2 PSD are realized simultaneously.
The data in Figure 12 illustrate impact of parameters (the main layer material is chosen similarly to the data in Fig. 11) on the properties of the considered systems for the choice of ENZ (Fig. 12a,b) and metallic (Fig. 12c,d) layers 1, 3. The graphs in Figure 12a,b and c,d are obtained for the values of d_{1,2,3} corresponding to the first maximum of functions Δϕ (d_{2}), Δϕ/ΔΦ (d_{2}), Δϕ/ΔΦ^{∼} (d_{2}) in Figures 11a,b and 5a,b, respectively. It is seen that the graphs in Figure 12a and b, c and d are “symmetric” relatively to the vertical axis determined by the condition (for the chosen TM polarization of the incident wave). Functions increase (Fig. 12a,с) or decrease (Fig. 12b,d) monotonically for the investigated ranges of values . For all the cases considered in Figure 12, one of functions is constant and the other one is characterized by one (Fig. 12a,b) or two (Fig. 12c,d) maxima. The condition Δϕ/ΔΦ ≫ 1 can be satisfied for rather wide ranges of values (with growth of ) only for ENZ layers 1, 3 (Fig. 12a,b). For metallic layers 1, 3, this condition is also possible but for the narrow ranges of values (for ) corresponding to the quick changes of functions (Fig. 12c, d). With that, for “metallic” parameters of the main layer materials (for ), the condition Δϕ ≈ ΔΦ takes place (Fig. 12c,d).
Fig. 11 Dependences of PSD Δϕ, ΔΦ, relations Δϕ/ΔΦ, Δϕ/ΔΦ^{∼} and energy transmission coefficients T_{1,2} of waves 1, 2 on the main layer thickness d_{2} (a, b), parameters (c, d) for a ENZ material of layers 1, 3 and anisotropic materials of the main layer. The values of parameters are: ε_{1} = ε_{3} = 0.01 + 0.01i (a, b); (c, d); d_{2} = 39 nm (c, d). The values of parameters for all the graphs are: ; ; d_{1,3} = 250 nm; λ = 0.63 μm; ε_{α,β} = 1; ψ = 0. Figures b and d detail the data in Figures a and c, respectively. The black horizontal dashed line in Figures b, d characterizes equality to 1 of the corresponding quantities. TM polarization of the incident wave is used. 
Fig. 12 Dependences of PSD Δϕ, ΔΦ, energy transmission coefficients T_{1,2} of waves 1, 2 on parameters (a, c) and (b, d) for various materials of layers 1, 3 and anisotropic materials of the main layer. The values of parameters are: (a) ε_{1} = ε_{3} = 0.01 + 0.01i, , , d_{1,3} = 250 nm, d_{2} = 28 nm, λ = 0.63 μm; (b) ε_{1} = ε_{3} = 0.01 + 0.01i, , , d_{1,3} = 250 nm, d_{2} = 28 nm, λ = 0.63 μm; (c) ε_{1} = ε_{3} = − 13 + 0.5i, , , d_{1,3} = 13 nm, d_{2} = 156 nm, λ = 0.6 μm; (d) ε_{1} = ε_{3} = − 13 + 0.5i, , , d_{1,3} = 13 nm, d_{2} = 156 nm, λ = 0.6 μm. The values of parameters for all the graphs are: ε_{α,β} = 1; ψ = 0. The black horizontal dashed line in all the figures characterizes equality to 0 of the corresponding quantities. TM polarization of the incident wave is used. 
4 Conclusion
The characteristic properties of transmission of electromagnetic waves through nanocomposite multilayers including ENZ MM or metallic spacer (upper and bottom) layers and the main layer made of two alternating strips of lowabsorbing dielectrics have been considered. The analytical modeling based on the exact solutions of electromagnetic boundary problems points to new possibilities of such nanoresonators usage as ultrathin phase transformers of the transmitted radiation. Counterintuitively, but the “addition” to the system of two isotropic spacer layers (that are the same for rays 1, 2, Fig. 1, and do not change the “one pass” phase difference of the transmitted waves) with the definite properties gives the opportunities to increase significantly the exact PSD (by an order of magnitude or more) of the transmitted waves. The numerical and graphical analysis shows that the considered effects are most pronounced under subwavelength multilayer thicknesses and accompanied by the significant amplification of multibeam interference inside of the main layer.
Two main regimes of the composite nanoresonator usage for amplification of the transmitted waves PSD are revealed: (i) for ENZ spacers when the conditions d_{2} ≪ d_{1,3}, ΔΦ ≪ Δϕ ≲ ΔΦ^{∼} are realized; (ii) for metallic spacers when d_{2} ≫ d_{1,3}, ΔΦ, ΔΦ^{∼} < Δϕ. These regimes were investigated for both model (ENZ and metallic) materials and silver (using two different wavelengths of the incident radiation) as spacer materials. For the considered ranges of the nanoresonator parameters, maximal values of PSD and transmission of the output waves are not realized simultaneously. But if the system is slightly tuned away from the maximal values, the parameter regions take place when PSD and transmission of the output waves are rather large for the whole multilayer thicknesses down to values of λ/2 and less (e.g. this is confirmed by the Figs. 4, 5, 7 data). With that, the main layer thicknesses can be less than λ/15 and λ/4 for the regimes (i) and (ii), respectively. So, the conditions are possible when decreasing the whole system thickness down to significantly subwavelength values can lead to substantial increasing PSD of the transmitted waves. The investigated transmission features also take place when changing the system parameters such as possible presence of absorption (losses) of the nanoresonator components, anisotropy of the main layer materials, nonideal realization of ENZ materials [41], oblique incidence of the exciting radiation (for small incidence angles).
Let us note that mechanisms of realization of large PSD for reflected and transmitted waves using the nanometric multilayers with ENZ components considered in [13,14] and the given paper are different. For reflection problems, the features of reflection of one of the waves from the ENZ substrate play a key role, and the interference effects are not main determining factors. For transmission, the effects of amplification of multibeam interference inside of the multilayer (that forms the nanoresonator) correspond to the conditions of obtaining large PSD of the transmitted waves. However, for both cases, the application of ENZ components for the considered layered systems leads to the quickly changing and extreme values of amplitude reflection or transmission coefficients depending on the system parameters. Namely that determines the main possibilities to amplify PSD of the output waves and decrease the whole system thickness.
The considered features of phase transformations of the transmitted waves can take place for various layered systems including the investigated nanoresonators or their modifications. The paper results also point to possibilities of exhibition of the investigated effects for the case when the main layer is originated by one optically anisotropic material. With that, the effective PSD amplification shall relate to the phase difference of proper waves (with different polarizations) in this material. One of the simplest methods to validate the obtained results is, apparently, a registration of the transmitted waves interference at the output of the considered systems, e.g. in the regimes of zeroorder phase plates of subwavelength thicknesses.
The obtained results can be applied to develop various ultrathin devices using phase transformation of transmitted waves, such as modulators, phase plates, filters, metasurfaces, nanoresonators, sensors, liquid crystals cells with reduced active material volumes, data processing system components, and for other nanophotonics applications. In general, together with the reflection problems analysis [13,14], the modeling results point to the new and perspective approaches (that can be simply realized within the frames of existing thin films technologies) for realization of multilayered nanometric optical elements including ENZ components for effective phase transformation of both reflected and transmitted radiation. Moreover, the considered “amplification” of the transmitted waves phase differences using pair of nanometric ENZ or metallic layers could be of interest for phasecontrast microscopy methods [42–44]. One can expect that addition of such suitable layers on the substrate and cover glass (between which there is a transparent sample, a “phase object”, e.g. biological cell) enhances weak phase contrasts of different regions of the sample and the image quality. That is of interest, in particular, for many biology and medicine applications of phasecontrast microscopy of transparent objects.
Acknowledgments
The author expresses the sincere appreciation to Dr Yuri Gritsai for the very fruitful discussions of the problem statement and results of this paper.
Appendix
The expressions for quantities A_{±}, B_{±}^{(1)}, C_{±}^{ (1)}, D_{±}, E_{±}^{ (1)}, F_{±}^{ (1)} in equations (7) and (8) take the form(A1) (A2) (A3) (A4) (A5) (A6)
In a similar manner, we have for quantities B_{±}^{ (2)}, C_{±}^{ (2)}, E_{±}^{ (2)}, F_{±}^{ (2)} in equations (9), (10) (A7) (A8) (A9) (A10)
The expressions for quantities A_{±}, D_{±} (Eqs. (A1), (A4)) are also true for equations (9) and (10). The relations for the parameters in equations (A1)(A10) are given in the main text.
References
 A.M. Urbas et al., J. Optics 18, 093005 (2016) [CrossRef] [Google Scholar]
 N.M. Estakhri, A. Alù, JOSA B 33, A21 (2016) [CrossRef] [Google Scholar]
 T. Cui, S. Liu, L. Zhang, J. Mater. Chem. C 5, 3644 (2017) [CrossRef] [Google Scholar]
 I. Liberal, N. Engheta, Nat. Photonics 11, 149 (2017) [CrossRef] [Google Scholar]
 S.A. Tretyakov, J. Optics 19, 013002 (2017) [CrossRef] [Google Scholar]
 P. Genevet et al., Optica 4, 139 (2017) [CrossRef] [Google Scholar]
 L. Huang, S. Zhang, T. Zentgraf, Nanophotonics 7, 1169 (2018) [CrossRef] [Google Scholar]
 VC. Su et al., Opt. Express 26, 13148 (2018) [CrossRef] [Google Scholar]
 V.S. Asadchy, A. DíazRubio, S.A. Tretyakov, Nanophotonics 7, 1069 (2018) [Google Scholar]
 M.A. Remnev, V.V. Klimov, Physics Uspekhi 61, 157 (2018) [CrossRef] [Google Scholar]
 M. Pu et al., Adv. Optical Mater. 7, 1801376 (2019) [CrossRef] [Google Scholar]
 S. Yoo, QH. Park, Nanophotonics 8, 249 (2019) [CrossRef] [Google Scholar]
 E.G. Starodubtsev, Mater. Res. Express 5, 126202 (2018) [CrossRef] [Google Scholar]
 E. Starodubtsev, EPJ Appl. Metamat. 6, 22 (2019) [CrossRef] [EDP Sciences] [Google Scholar]
 R.W. Ziolkowski, E. Heyman, Phys. Rev. E 64, 056625 (2001) [CrossRef] [Google Scholar]
 A. Alù et al., Phys. Rev. B 75, 155410 (2007) [CrossRef] [Google Scholar]
 I.V. Lindell, A.H. Sihvola, PIER 89, 213 (2009) [CrossRef] [Google Scholar]
 E.G. Starodubtsev, Metamaterials 4, 32 (2010) [CrossRef] [Google Scholar]
 R. Maas et al., Nat. Photonics 7, 907 (2013) [CrossRef] [Google Scholar]
 N. Kinsey et al., Nat. Rev. Materials 1 (2019) [Google Scholar]
 C. Valagiannopoulos et al., J. Appl. Phys. 121, 063102 (2017) [CrossRef] [Google Scholar]
 J. Kim et al., Nano Lett. 18, 740 (2018) [CrossRef] [Google Scholar]
 V. Caligiuri et al., Nanophotonics 8, 1505 (2019) [CrossRef] [Google Scholar]
 V. Caligiuri et al., Nano Lett. 19, 3151 (2019) [CrossRef] [Google Scholar]
 F.I. Fedorov, Optics of Anisotropic Media , 2nd edn. (Editorial URSS, Moscow, 2004 (in Russian)) [Google Scholar]
 M. Schadt, Annu. Rev. Mater. Sci. 27, 305 (1997) [CrossRef] [Google Scholar]
 DK. Yang, ST. Wu, Fundamentals of Liquid Crystal Devices , 2nd edn. (Wiley, 2014) [CrossRef] [Google Scholar]
 R. Häussler et al., Appl. Optics 56, F45 (2017) [CrossRef] [Google Scholar]
 Z. Yue et al., Nat. Commun. 8, 15354 (2017) [CrossRef] [Google Scholar]
 A.M. Shaltout et al., Nat. Commun. 9, 2673 (2018) [CrossRef] [Google Scholar]
 T.W. Ebbesen et al., Nature 391, 667 (1998) [CrossRef] [Google Scholar]
 E. Popov et al., Phys. Rev. B 62, 16100 (2000) [CrossRef] [Google Scholar]
 A. Degiron et al., Appl. Phys. Lett. 81, 4327 (2002) [CrossRef] [Google Scholar]
 D. Xiang et al., Opt. Commun. 284, 471 (2011) [CrossRef] [Google Scholar]
 A.S. Marathay, J.F. McCalmont, JOSA A 18, 2585 (2001) [CrossRef] [Google Scholar]
 V. Kivijärvi et al., New J. Phys. 17, 063019 (2015) [CrossRef] [Google Scholar]
 J.M. Coupland, N.I. Nikolaev, Opt. Express 27, 30380 (2019) [CrossRef] [Google Scholar]
 C. Valagiannopoulos, PIER 71, 59 (2007) [CrossRef] [Google Scholar]
 A. Yariv, P. Yeh, Optical Electronics in Modern Communications , 6th edn. (Oxford University Press, 2007) [Google Scholar]
 J. Smajic et al., J. Comput. Theor. Nanosci. 6, 1 (2009) [CrossRef] [Google Scholar]
 E. Starodubtsev, EPJ Appl. Metamat. 5, 1 (2018) [CrossRef] [EDP Sciences] [Google Scholar]
 F. Zernike, Science, New Series 121, 345 (1955) [Google Scholar]
 M. Teschke, S. Sinzinger, JOSA A 26, 1015 (2009) [CrossRef] [Google Scholar]
 K. Toda et al., Sci. Reports 9, 9957 (2019) [Google Scholar]
Cite this article as: Evgenii Starodubtsev, Features of transmission of electromagnetic waves through composite nanoresonators including epsilonnearzero metamaterials, EPJ Appl. Metamat. 7, 1 (2020)
All Figures
Fig. 1 (a) Scheme of the multilayer with the main composite layer formed from alternating “strips” of two different materials that can be optically uniaxial in the general case. The rays in the substrate transmitted through the multilayer are denoted by numbers 1 and 2. E_{in} is the incident wave field, the reflected rays in the superstrate are not shown. (b) Geometry of the corresponding electromagnetic boundary problems. The multilayer includes upper and bottom layers, which electromagnetic properties can correspond to isotropic ENZ MMs or metals, and the main layer. The superstrate and substrate materials are isotropic and optically thick. E_{r1}, E_{r2}, and E_{t1}, E_{t2} are the fields of the waves reflected from and transmitted through Strips 1 and 2, respectively, when the incident wave interacts with the multilayer. OA denotes the optical axis direction in the materials of Strips 1 or 2. The other denotations are given in the main text. 

In the text 
Fig. 2 Dependences of phases ϕ_{1,2} (a–c) and energy transmission coefficients T_{1,2} (d–f) of waves 1, 2 on the main layer thickness d_{2} for various materials of layers 1, 3. The values of parameters are: ε_{1} = ε_{3} = 2 + 10^{−3} i (a,d); ε_{1} = 2 + 10^{−3} i, ε_{3} = 0.01 + 0.01i (b,e); ε_{1} = ε_{3} = 0.01 + 0.01i (c,f). The values of parameters for all the graphs are: d_{1,3} = 0.25 μm, λ = 0.63 μm. Parameter d_{2} is used as the abscissa axis for all the graphs. 

In the text 
Fig. 3 The effect of parameter d_{2} on the real and imaginary parts of amplitude transmission coefficients of waves 1, 2 (a, b) and trajectories of the tips of vectors t_{1,2} on the complex plane depending on parameters d_{2} (c–f) and (g,h). The values of parameters are: ε_{1} = ε_{3} = 2 + 10^{−3} i (a, c, e); ε_{1} = ε_{3} = 0.01 + 0.01i (b, d, f); , d_{2} = 39 nm (g, h). The values of parameters for all the graphs are: d_{1,3} = 0.25 μm, λ = 0.63 μm. 

In the text 
Fig. 4 Dependences of PSD Δϕ, ΔΦ, relations Δϕ/ΔΦ, Δϕ/ΔΦ^{∼} and energy transmission coefficients T_{1,2} of waves 1, 2 on the main layer thickness d_{2} (a, c, e), wavelength λ (b, d), parameters (f, h), d_{1,3} (g) for various materials of layers 1, 3. The values of parameters are: ε_{1} = ε_{3} = 2 + 10^{−3} i (a, b); ε_{1} = ε_{3} = 0.01 + 0.01i (c, d, e, g); d_{2} = 39 nm (b, d, f, h); (f, h). The values of parameters for all the graphs are: d_{1,3} = 0.25 μm (besides g); λ = 0.63 μm (besides b, d). Figures e and h detail the data in Figures c and f, respectively. The black horizontal dashed line in Figures e, g, h characterize equality to 1 of the corresponding quantities. 

In the text 
Fig. 5 Dependences of PSD Δϕ, ΔΦ, relations Δϕ/ΔΦ, Δϕ/ΔΦ^{∼} and energy transmission coefficients T_{1,2} of waves 1, 2 on the main layer thickness d_{2} (a, b), parameters (c), d_{1,3} (d) for a metallic material of layers 1, 3. The values of parameters are: ε_{1} = ε_{3} = − 13 + 0.5i (a, b, d); d_{2} = 156 nm (c, d); (c). The values of parameters for all the graphs are: d_{1,3} = 13 nm (besides d), λ = 0.6 μm. Figure b details the data in Figure a. The black horizontal dashed line in fig. b characterizes equality to 1 of the corresponding quantities. 

In the text 
Fig. 6 Dependences of phase differences ϕ_{1,2} (a, b) and energy transmission coefficients T_{1,2} (c, d) of waves 1, 2 on the main layer thickness d_{2} for various materials of layers 1, 3. The values of parameters are: ε_{1} = 2 + 10^{−3} i, ε_{3} = − 13 + 0.5i (a, c); ε_{1} = ε_{3} = − 13 + 0.5i (b, d). The values of parameters for all the graphs are: d_{1,3} = 13 nm, λ = 0.6 μm. Parameter d_{2} is used as the abscissa axis for all the graphs. 

In the text 
Fig. 7 Dependences of PSD Δϕ, ΔΦ, relations Δϕ/ΔΦ, Δϕ/ΔΦ^{∼} and energy transmission coefficients T_{1,2} of waves 1, 2 on thicknesses d_{2} (a, b, d, e), d_{1,3} (c, f) for silver layers 1, 3 and two different wavelengths of the incident wave. The values of the Ag permittivity for λ = 292 nm (a, b, c), 650 nm (d, e, f) are given in the text. The values of parameters are: d_{1,3} = 125 nm (a, b), 9 nm (d, e); d_{2} = 19 nm (c), 171 nm (f). Figures b and e detail the data in figures a and d, respectively. The black horizontal dashed line in figures b, c, e characterizes equality to 1 of the corresponding quantities. 

In the text 
Fig. 8 Dependences of PSD Δϕ, ΔΦ, relations Δϕ/ΔΦ, Δϕ/ΔΦ^{∼} and energy transmission coefficients T_{1,2} of waves 1, 2 on incidence angle ψ for silver layers 1, 3, two different wavelengths and TM (a, b, e, f), TE (c, d, g, h) polarizations of the incident wave. The values of parameters are: (a–d) λ = 292 nm, d_{1,3} = 125 nm, d_{2} = 19 nm; (e–h) λ = 650 nm, d_{1,3} = 9 nm, d_{2} = 171 nm. Figures b, d, f, h detail the data in figures a, c, e, g, respectively. The black horizontal dashed line in figures b, d characterizes equality to 1 of the corresponding quantities. Parameter ψ is used as the abscissa axis for all the graphs. Values of ψ (as Δϕ, ΔΦ, ΔΦ^{∼}) are given in the relative units of radians/π. 

In the text 
Fig. 9 Dependences of quantities ϕ_{1}, R_{21}, R_{23}, I ^{ (1)} on the main layer thickness d_{2} for various materials of layers 1, 3. The values of parameters are: ε_{1} = ε_{3} = 2 + 10^{−3} i (a, c), 0.01 + 0.01i (b), −13 + 0.5i (d); d_{1,3} = 250 nm (a, b), 13 nm (c, d); λ = 0.63 μm (a, b), 0.6 μm (c, d). Denotations of the curves in Figures b, d are given in Figures a, c. I ^{ (1)} = I ^{ TM (1)} (z = d_{1} + 0.5d_{2}) for all the graphs. Parameter d_{2} is used as the abscissa axis for all the graphs. 

In the text 
Fig. 10 Dependences of quantities ϕ_{1}, R_{21}, R_{23}, I ^{ (1)} on parameters (a, c), d_{1,3} (b, d) for various materials of layers 1, 3. The values of parameters are: (a, b), 0.5 (c, d); (b), −13 (d); d_{2} = 39 nm (a, b), 156 nm (c, d); d_{1,3} = 250 nm (a, b), 13 nm (c, d); λ = 0.63 μm (a, b), 0.6 μm (c, d). I ^{ (1)} = I ^{ TM (1)} (z = d_{1} + 0.5d_{2}) for all the graphs. The black horizontal dashed line in figs. a, b, c characterizes equality to 1 of the corresponding quantities. 

In the text 
Fig. 11 Dependences of PSD Δϕ, ΔΦ, relations Δϕ/ΔΦ, Δϕ/ΔΦ^{∼} and energy transmission coefficients T_{1,2} of waves 1, 2 on the main layer thickness d_{2} (a, b), parameters (c, d) for a ENZ material of layers 1, 3 and anisotropic materials of the main layer. The values of parameters are: ε_{1} = ε_{3} = 0.01 + 0.01i (a, b); (c, d); d_{2} = 39 nm (c, d). The values of parameters for all the graphs are: ; ; d_{1,3} = 250 nm; λ = 0.63 μm; ε_{α,β} = 1; ψ = 0. Figures b and d detail the data in Figures a and c, respectively. The black horizontal dashed line in Figures b, d characterizes equality to 1 of the corresponding quantities. TM polarization of the incident wave is used. 

In the text 
Fig. 12 Dependences of PSD Δϕ, ΔΦ, energy transmission coefficients T_{1,2} of waves 1, 2 on parameters (a, c) and (b, d) for various materials of layers 1, 3 and anisotropic materials of the main layer. The values of parameters are: (a) ε_{1} = ε_{3} = 0.01 + 0.01i, , , d_{1,3} = 250 nm, d_{2} = 28 nm, λ = 0.63 μm; (b) ε_{1} = ε_{3} = 0.01 + 0.01i, , , d_{1,3} = 250 nm, d_{2} = 28 nm, λ = 0.63 μm; (c) ε_{1} = ε_{3} = − 13 + 0.5i, , , d_{1,3} = 13 nm, d_{2} = 156 nm, λ = 0.6 μm; (d) ε_{1} = ε_{3} = − 13 + 0.5i, , , d_{1,3} = 13 nm, d_{2} = 156 nm, λ = 0.6 μm. The values of parameters for all the graphs are: ε_{α,β} = 1; ψ = 0. The black horizontal dashed line in all the figures characterizes equality to 0 of the corresponding quantities. TM polarization of the incident wave is used. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.