# Broadband mixing of -symmetric and -broken phases in photonic heterostructures with a one-dimensional loss/gain bilayer

## Abstract

Combining loss and gain components in one photonic heterostructure opens a new route to efficient manipulation by radiation, transmission, absorption, and scattering of electromagnetic waves. Therefore, loss/gain structures enabling -symmetric and -broken phases for eigenvalues have extensively been studied in the last decade. In particular, translation from one phase to another, which occurs at the critical point in the two-channel structures with one-dimensional loss/gain components, is often associated with one-way transmission. In this report, broadband mixing of the -symmetric and -broken phases for eigenvalues is theoretically demonstrated in heterostructures with four channels obtained by combining a one-dimensional loss/gain bilayer and one or two thin polarization-converting components (PCCs). The broadband phase mixing in the four-channel case is expected to yield advanced transmission and absorption regimes. Various configurations are analyzed, which are distinguished in symmetry properties and polarization conversion regime of PCCs. The conditions necessary for phase mixing are discussed. The simplest two-component configurations with broadband mixing are found, as well as the more complex three-component configurations wherein symmetric and broken sets are not yet mixed and appear in the neighbouring frequency ranges. Peculiarities of eigenvalue behaviour are considered for different permittivity ranges of loss/gain medium, i.e., from epsilon-near-zero to high-epsilon regime.

###### pacs:

Being a weaker condition compared to hermiticity, -symmetry can still yield real positive eigenvalues for Schrödinger’s equation (1); (2); (3). The investigations of -symmetry have not been restricted to the area of quantum mechanics. The realization of -symmetric structures in optics has been suggested in the 2000s. In Ref. dccclxxx((4)), a parallel-plate waveguide with -symmetric interior has been proposed. A large portion of the optics relevant theoretical studies has been focused on one-dimensional lattices with gradual variation of the refractive index (5); (6); (7). Moreover, -symmetry has been experimentally observed in optical lattices (8); (9). Later, conservation relations and -transition properties of one-dimensional photonic heterostructures (10); (11) have been studied. In this case, the eigenvalues of the -matrix are unimodular and flux conserving in -symmetric phase. In -broken phase, they have different magnitudes, one of which is larger and the remaining one is smaller than unity (11). They are known to correspond to amplification and attenuation, respectively, and satisfy the generalized conservation condition (10); (11). The two above-mentioned phases are separated by the spontaneous -symmetry breaking point (called critical, or exceptional, or phase transition point), at which -matrix has only one eigenvalue instead of the two ones. In this regime, the system can be fully transparent for the light incident from one side of the structure while reflectivity is enhanced for the opposite-side incidence (12). The role of the critical point for achieving unidirectional transmission has also been highlighted in many other theoretical and experimental studies, e.g., see Refs. dccclxxx((13); (14); (15); (16)). Unidirectional invisibility and cloaking in the structures having -symmetric components have been studied in detail (12); (17); (18). Moreover, lenses and cavities with -symmetry should be mentioned (19); (20); (21). While bilayers, multilayers, and waveguides with a finite-extent cross section were commonly considered in the studies of optical structures with -symmetry, the attention has also been attracted by two-dimensional and quasiplanar one-dimentional arrays (22); (23); (24); (25). Specifics of the epsilon-near-zero (ENZ) range of permittivity has been intestigated for the bilayer structures (26); (27). In particular, conditions of tunneling through -symmetric ENZ bilayers have been considered (26). Some exotic regimes like loss-induced superscattering and gain-induced absorption (24) are also worth mentioning. At the same time, new opportunities are promised by recent advances in full control of electromagnetic waves with the aid of metasurfaces (28); (29); (30); (31); (32). In particular, some extreme regimes of polarization conversion have recently been demonstrated (28); (33); (34); (35), whereas some others are expected to be achieved soon. This opens a new route to multichannel structures with -related properties, which do not need two-dimensional loss/gain components. Indeed, although the principal possibility of co-existence of -symmetric and -broken phases is a known feature for the case of physical dimensions higher than one (and, hence, more than two input and output channels may exist)(11), there is high demand in simpler structures that would enable easily realizable scenarios.

In this report, we theoretically demonstrate the co-existence of -symmetric and -broken sets of eigenvalues of the -matrix in a four-channel photonic heterostructure with a one-dimensional bilayer, while additional two channels required for the phase mixing are created due to thin polarization converting componenets (PCCs) that are placed at the end-faces. Our main goal is to show the principal possibility and find the basic features of this regime. Therefore, in the contrast with a two-dimensional case of loss/gain media, the function of the creation of additional transmission channels and that of combining gain and loss are in space, so far as they are performed by the different components of the entire structure. There are two input and two output channels (ports) in the studied heterostructures, whereas each of the PCCs is assumed to be passive, lossless, and quasiplanar. To the best of our knowledge, the possibility of the simultaneous existence of -symmetric and -broken phases in a photonic system with one-dimensional loss/gain component has not yet been discussed in the literature. We refer to the co-existence of -symmetric and -broken phases in the same frequency range as the phase mixing. The goals of this study include the validation of the suggested approach, finding and comparison of different configurations enabling the mixing of symmetric and symmetry-broken phases, and the formulation of the conditions providing the mixing at the minimal number of the structural components. These goals can be achieved by using the standard -matrix formalism, while advanced techniques are not required. We investigate a wide range of variation of permittivity of the loss and gain components. Among others, it includes the ENZ region, which is known for such exotic properties as squeezing, supercoupling, stretching of wavelength, enhancing nonreciprocity, and time-reversal symmetry breaking (36). Here, the study of the ENZ range and consideration of PCCs are restricted according to the goals of this paper. Design of PCCs is beyond its scope.

## Results

General Model.

The studied photonic heterostructures are assumed to have in the general case the following three components: The loss/gain bilayer and two PCCs, one at each of its end-faces. Capability in polarization conversion is directly connected with the possibility of the creation of additional channels required for obtaining a four-channel configuration. Different configurations can be accessed through adding or removing one or two PCCs. The general schematic of the structure is shown in Fig. 1(a). The region is the loss/gain medium that satisfies the symmetry condition . It is noteworthy that this condition is necessary but not sufficient for -symmetry. For the purposes of our study, we assume that PCCs are infinitesimally thin and capable of converting right/left circular polarization (CP) to left/right CP perfectly. Here, we consider three different types of PCCs. The first-type PCCs denoted by and change the polarization state for both reflected and transmitted waves. The second-type PCCs (denoted by and ) only change the transmitted waves’ polarization, whereas the third-type PCCs (denoted and ) only change the reflected waves’ polarization. The properties of PCCs of the three types are illustrated in Fig. 1(b). Polarization state cannot be changed by using only bilayer. All of the studied structures are free-standing structures adjusted with the vacuum half-spaces.

We build our formalism for a general 1D photonic heterostructure that has four channels, two input and two output, one allowing right circularly polarized light and the other allowing left circularly polarized light to pass through. We start from casting the most general expression for the electric field along the -direction that consists of right and left polarized waves. For the left-side half-space, , and normal incidence, it is given by

(1) |

where , are scattering amplitudes, and , with and , , is angular frequency, is velocity of electromagnetic wave, is index of refraction. To obtain for the right-side half-space, i.e., at , we replace and in Eq. (1). Similar expressions can be used for the bilayer region.

By exploiting the symmetry properties of the structure, we can cast the for the four-channel case. It is expected to be and really identical to such a well-known relation in the two-channel case, which is valid everywhere, except for coherent-perfect-absorption laser points, see Ref. dccclxxx((11)). It is connected with another key relation known as the generalized unitarity relation. Starting from the condition , which also implies , and using the standard -matrix formalism, which is commonly used in quantum mechanics and electromagnetic theory, we obtain the conservation relation in the studied four-channel case as follows:

(2) |

where is the transmittance, and are the reflectances
at left-side and right-side illumination, with no restrictions imposed on them; is transmission coefficient, and are reflection coefficients at left-side and right-side illumination. In line with the -matrix formalism, these coefficients are fully determined by the scattering amplitudes.
To achieve the purposes of this study, the conventional -matrix approach can be used for the arbitrary values
of , , and .

Various Configurations.
There are two general cases that are distinguished in terms of the physics they offer. The first one is the case, in which -symmetric and -broken eigenvalues may not exist simultaneously, and one phase is translated into another at the critical frequency (Case 1). In the second case, the
co-existence of the symmetric and symmetry-broken eigenvalues is possible, i.e., a may occur
(Case 2). This case is the focus of our study. It will be shown that there are various configurations
that differ in the end-faces of the region , which can be utilized for accessing Case 1 and Case 2.
We start from the simplest configuration, and then consider
more complex configurations by placing one or both of the components and at the end-faces.

Case 1: The simplest configuration to access this case is a bilayer enabling -symmetry (like the region ), which has been studied in detail in Refs. dccclxxx((10); (11)). For the configuration , -matrix has two eigenvalues. Let us generalize it by formally adding the polarization related degree of freedom to the system to obtain the -matrix, i.e., present it in the same form as used throughout the paper for more complex configurations. Since there are no PCCs in this case, {ceqn}

(3) |

yields the same set of eigenvalues twice, with the two eigenvalues in a set, which are given by {ceqn}

(4) |

Next, we add the PCCs to the left and right end-faces, which are assumed to be capable of changing the polarization for both reflected and transmitted waves ( and ). The configuration that we have now is , meaning that if the wave is transmitted it retains the initial polarization because of passing through two PCCs. In turn, if it is reflected, right/left CP is changed to left/right CP. Now, the two diagonal blocks of the -matrix in Eq. (3) are coupled due to the added PCCs. After some algebra, we obtain {ceqn}

(5) |

The four eigenvalues for this configuration are given as {ceqn}

(6) |

Thus, configurations and show the same moduli for the eigenvalues and the same basic physics. It is well-known that a preserves the eigenvalues. Hence, all other configurations, for which -matrices are connected by a unitary transformation with -matrix of one of two above discussed configurations, will have the same eigenvalues. Therefore, these configurations also belong to Case 1, in which the mixed phase for the eigenvalues is not possible. Indeed, they are either -symmetric or -broken at any fixed frequency, except at the critical point. We can gain access to such configurations by including different combinations of PCCs at one or both end-faces, so -symmetry is spontaneously broken at the critical point. We showed in our study that they include the configurations , , , , , and . Here, we do not discuss each of these cases separately, because eigenvalues in all the cases are the same as for the configurations and , i.e., they are given by Eq. (6). It is interesting that the configurations and do not yield the mixed phase despite the fact that polarization conversion is possible at one of the end-faces. Since the components and only change the polarization of transmitted waves, the ability of polarization conversion of the reflected wave is expected to be the (but not sufficient) condition of existence of the mixed phase. It is noticeable that the same -matrices and the same eigenvalues can be obtained in Case 1 for the structures with PCCs at two end-faced, with a PCC at one end-face, and without PCCs.

An example is presented in Fig. 2. At , all four eigenvalues are (log10, =1,2,3,4), being in -symmetric phase and, hence, are flux conserving. At , the eigenvalues become reciprocal in two sets, which are in the symmetry-broken phase and, hence, satisfy the generalized conservation relation. These properties are identical to those well-known for one-dimensional structures with two channels. One can see that it is not possible to simultaneously obtain the symmetric and symmetry-broken phases at any fixed value of .

Case 2: The simplest intuition for accessing this case says that a PCC should be added only to one of the end-faces of the component . However, it cannot ensure the phase mixing, as follows from the study of Case 1. Thus, let us first add a PCC, specifically to the left end-face, i.e., we now have configuration . For the light incident from the left side, polarization of both the reflected and the transmitted wave is assumed to be changeable due to this PCC. If it is incident from the right side, the transmitted wave’s polarization is still changed, while the reflected wave retains its polarization. These properties are described by the following -matrix: {ceqn}

(7) |

This configuration yields the mixed phase for the eigenvalues of the -matrix, so we can obtain symmetric and symmetry-broken sets of eigenvalues at fixed . The four eigenvalues corresponding to Eq. (7) are given by

(8a) | |||

(8b) |

As one set of eigenvalues () preserves the symmetry and unimodularity at , the other set of eigenvalues () is translated into a symmetry-broken phase in a wide -range, even near . Once is reached, this symmetry-broken set of the eigenvalues starts to overlap with the first set, whose eigenvalues are also in the symmetry-broken phase at that time. Thus, all of the eigenvalues are in the broken phase at . So, is the boundary between the mixing phase and the all-broken phase cases.

Fig. 3 and Fig. 4 show the maps of the magnitudes of eigenvalues , , in log10 scale, which are obtained for Case 2 from Eqs. (8a) and (8b), in wide ranges of variation in and . It is clearly seen that the symmetry of eigenvalues is broken even at small frequencies. Thus, there is a large region in -plane, where -symmetric and -broken sets of eigenvalues may co-exist. Starting from , the set of eigenvalues also experiences spontaneous symmetry breaking, so the mixing does not exist anymore.

As shown in Fig. 3(a,b), the value of and the extent to which the symmetry is broken strongly depend on . For instance, we obtain in the vicinity of . For the set of eigenvalues , a strong deviation from the unimodular case is observed in Fig. 3(c,d) even at very small values of . In particular, a strong anomaly of occurs nearly at and . Hence, -symmetric and -broken eigenvalues may co-exist also in ENZ regime, i.e., in the close vicinity of . It is noteworthy that the behaviours of eigenvalues in Case 1 and Case 2 at are very different. In Fig. 4, the boundary between the regions with the mixed phase for eigenvalues (at smaller ) and with all symmetry-broken eigenvalues (at larger ) is clearly seen. Its location can be controlled by variations in , while is fixed. A detailed investigation of these scenarios will be a subject of our future research.

Fig. 5 and Fig. 6 display the magnitudes for eigenvalues , , and in log10 scale at varying for and , respectively, for two different values of . The location of the critical frequency for and , , and, hence, width of the -range, in which the mixed phase is achieved, are strongly affected by . As the value of increases, the critical frequency is redshifted, whereas an increase of results in blueshift of . At the same time, weakly affects the width of mixing range. On the other hand, can strongly affect the magnitudes of and . -symmetric and -broken sets of eigenvalues can co-exist even at , although the difference between them is rather weak in this case. Since breaking symmetry at is a general property of and , it also occurs in the ENZ regime, e.g. at (not shown). However, in this regime, the difference between and at is weak, whereas and may significantly deviate from the unimodular case in a wide range of variation.

Next, let us add the component to the right end-face of the bilayer, so we now have the configuration . Then, we obtain exactly the same physical scenario for configuration , while the roles of right and left side illuminations are interchanged. The -matrix is now given as: {ceqn}

(9) |

It has four eigenvalues, which may yield the mixed phase:

(10a) | |||

(10b) |

eigenvalues | components | mixed phase | |
---|---|---|---|

2 | 1 | ✷ | |

4 | 3 | ✷ | |

or | 4 | 2 | ✳ |

4 | 3 | ✷ | |

4 | 3 | ✷ | |

or | 4 | 2 | ✷ |

or | 4 | 2 | ✳ |

or | 4 | 3 | ✳ |

or | 4 | 3 | ✷ |

or | 4 | 3 | ✳ |

Similarly to Case 1, we can obtain the sets of eigenvalues given by Eqs. (8a), (8b) and by Eqs. (10a), (10b), in different configurations. The -matrices of these configurations should be equivalent to the -matrices of the and cases, respectively, from which they must be obtainable via unitary transformations to yield the same eigenvalues. From the obtained results, the first set of eigenvalues [Eqs. (8a), (8b)], can be accessed in the configurations , , and , whereas the second set of eigenvalues [Eqs. (10a), 10b)] does in , , and . The simplest configurations for Case 2 are , , , and . They only need two components, i.e., a loss/gain bilayer and one PCC.

The properties of the studied configurations are summarized in Table 1. One can see therein which features are required for obtaining the mixing phase for eigenvalues, and which are for keeping only -symmetric eigenvalues within a certain frequency range. Structures with two components - one loss/gain bilayer and one PCC with peculiar properties - can be sufficient for obtaining the mixing phase, whereas those with three components do not always lead to it. Whether the mixing of the phase is achieved or not depends on the polarization conversion scenario at the end-faces, so the role of PCC(s) is evident. If it is achieved, it occurs in a rather wide frequency range. It is noticeable that despite having different end-faces of the studied configurations still enable the sets of symmetric eigenvalues. On the other hand, the effect exerted by one PCC can be compensated by that of the other, so phase mixing is not achieved in some three-component configurations. This case can be useful for the separation of two processes in one structure and, therefore, is promising for multifunctional operation.

## Discussion

The main goal of this study was to show that broadband mixing is possible, and, moreover, it may be achieved in photonic heterostructues with a one-dimensional gain/loss bilayer. We demonstrated a way to broadband mixing of -symmetric and -broken phases for eigenvalues in the structures with four channels and with a one-dimensional loss/gain component. As far as symmetric and broken phases show different properties related to flux conservation, amplification, and attenuation, it is expected that the broadband mixing of these phases may open new routes to efficient selective manipulation by electromagnetic radiation, including advanced regimes of directional selectivity, enhancement, and absorption. Earlier, such a mixing has been expected to occur for physical dimensions higher than one. The obtained results show that a two-dimensional loss/gain medium and, moreover, anisotropy are not required for the mixing. We have analytically derived the eigenvalues for different configurations of photonic heterostructures consisting of a bilayer of one-dimensional loss/gain medium and PCC(s), and proved the possibility of obtaining the mixed phase for eigenvalues in a wide but limited frequency range, while its width depends on both the real and the imaginary part of the permittivity of the loss/gain component. Therefore, the wideband phase mixing is a very general effect whose existence does not need any special adjustment of the parameters, although its appearance can be strongly sensitive to the parameter choice. While the principal possibility of obtaining four channels by replacing a two-dimensional loss/gain medium with a one-dimensional medium combined with a PCC could be expected, the obtained results indicate that this case can be achieved in different permittivity ranges for multiple configurations, which are distinguished due to the presence and properties of the PCC(s). The utilized model based on -matrix formalism properly describes a wide variety of photonic heterostructures with a 1D loss/gain component and their -related properties. It is noteworthy that the -symmetric eigenvalues set may exist for all of the considered configurations, regardless of whether the loss/gain component is end-faced with one or two PCCs, or not end-faced at all. According to the goals of this study, we clarified, based on the obtained results, which properties of PCCs are necessary and which are sufficient in different transmission/reflection scenarios. We showed that a two-component heterostructure can be sufficient to obtain wideband phase mixing, provided that the PCC(s) show the suitable properties. On the other hand, we detected such three-component heterostructures that keep immunity against phase mixing, in spite of containing a PCC that may lead to the mixing in other, even simpler configurations. Design of PCCs enabling the desired polarization manipulation will be one of the next steps. A further study of the peculiarities of different ranges of material parameters is planned. The difference between the cases with and without phase mixing for eigenvalues is especially intriguing for ultralow-permittivity regime, in which the mixing is possible even in close vicinity of . Thus, the consideration of the studied physical features in connection with recent advances in theory and applications of ENZ materials and specific coupling and transmission regimes realized with their use is a promising topic for future research.

## Methods

The standard -matrix formalism, which is commonly used in quantum mechanics, has been used to mathematically express transmission and reflection properties of different configurations. The eigenvalues of the -matrix were analyzed analytically. The values of , and were calculated by applying the conditions of continuity for the tangential components of electromagnetic field at the boundaries of the structural components, so the unknown coefficients in the general formulas for the field components are introduced unambiguously. The eigenvalues were calculated for various configurations via a simple custom-made Fortran code and the results are plotted via GNUPLOT. Both of these procedures were performed on a standard laptop that works on an Ubuntu operating system.

## Acknowledgements

This work is supported by the projects DPT-HAMIT and NATO-SET-193. One of the authors (E.O.) also acknowledges partial support from the Turkish Academy of Sciences. A.E.S. thanks the National Science Centre of Poland for support under the project MetaSel DEC-2015/17/B/ST3/00118. Work at Ames Laboratory was partially supported by the U.S. Department of Energy, Office of Basic Energy Science, Division of Materials Sciences and Engineering (Ames Laboratory is operated for the U.S. Department of Energy by Iowa State University under Contract No. DE-AC02-07CH11358). The European Research Council under ERC Advanced Grant No. 320081 (PHOTOMETA) supported work at FORTH. This work is origally published in Scientific Reports [Sci Rep. 2017; 7: 15504.].

### References

- Bender, C. M. & Boettcher, S. Real Spectra in Non-Hermitian Hamiltonians Having PT Symmetry. Phys. Rev. Lett. 80, 5243-5246 (1998).
- Bender, C. M., Boettcher, S., & Meisinger, P. N. PT-symmetric quantum mechanics. J. Math. Phys. 40, 2201-2229 (1999).
- Bender, C. M. Making sense of non-Hermitian Hamiltonians. Rep. Prog. Phys. 70, 947-1018 (2007).
- Ruschhaupt, A., Delgado, F., & Muga J. G., Physical realization of PT-symmetric potential scattering in a planar slab waveguide, J. Phys. A: Math. Gen. 38, L171-L176 (2005).
- Makris, K. G., El-Ganainy, R., Christodoulides, D. N., & Musslimani, Z. H. Beam Dynamics in PT Symmetric Optical Lattices. Phys. Rev. Lett. 100, 103904 (2008).
- Makris, K. G., El-Ganainy, R., Christodoulides, D. N., & Musslimani, Z. H. PT-symmetric optical lattices. Phys. Rev. A 81, 063807 (2010).
- Musslimani, Z. H., Makris, K. G., El-Ganainy, R., & Christodoulides, D. N. Optical Solitons in PT Periodic Potentials. Phys. Rev. Lett. 100, 030402 (2008).
- Guo, A. et al. Observation of PT-Symmetry Breaking in Complex Optical Potentials. Phys. Rev. Lett. 103, 093902 (2009).
- Ruter, C. E. et al. Observation of parityâtime symmetry in optics. Nature Phys. 6, 192-195 (2010).
- Chong, Y. D., Ge, L., & Stone, A. D. PT-Symmetry Breaking and Laser-Absorber Modes in Optical Scattering Systems. Phys. Rev. Lett. 106, 093902 (2011).
- Ge, L., Chong, Y. D., & Stone, A. D. Conservation relations and anisotropic transmission resonances in one-dimensional PT-symmetric photonic heterostructures. Phys. Rev. A 85, 023802 (2012).
- Lin, Z., Ramezani, H., Eichelkraut, T., Kottos, T., Cao, H., & Christodoulides, D. N. Unidirectional invisibility induced by PT-symmetric periodic structures, Phys. Rev. Lett. 106, 213901 (2011).
- Xu, Y.-L., Feng, L., Lu, M.-H., & Chen, Y.-F. Asymmetric optical mode conversion and transmission by breaking PT-symmetry on silicon photonic circuits, Phys. Lett. A 376, 886-890 (2012).
- Hahn, C., et al. Unidirectional Bragg gratings using Parity-Time symmetry breaking in plasmonic systems, IEEE J. Quant. Electron. 52, 4600712 (2016).
- Yin, X. & Zhang, X. Unidirectional light propagation at exceptional points, Nat. Mat. 12, 175-177 (2013).
- Feng, L., et al. Experimental demonstration of a unidirectional reflectionless parity-time metamaterial at optical frequencies. Nat. Mat. 12, 108-113 (2013).
- Sounas, D. L., Fleury, R., & Alù, A. Unidirectional cloaking based on metasurfaces with balanced loss and gain. Phys. Rev. Appl. 4, 014005 (2015).
- Zhu, X., Feng, L., Zhang, P., Yin, X., & Zhang, X. One-way invisible cloak using parity-time symmetric transformational optics, Opt. Lett. 38, 2821-2824 (2013).
- Valagiannopoulos, C. A., Monticone, F., & Alù, A. PT-symmetric planar devices for field transformation and imaging, J. Opt. 18, 044028 (2016).
- Peng, B., et al. Parity-time-symmetric whispering-gallery microcavities, Nat. Phys. 10, 394-398 (2014).
- Kulishov, M., Kress, B., & Slavik, R. Resonant cavities based on Parity-Time-symmetric diffraction gratings, Opt. Express 21, 9473-9483 (2013).
- Turduev, M., Botey, M., Giden, I., Herrero, R., Kurt, H., Ozbay, E., & Staliunas, K. Two-dimensional complex parity-time-symmetric photonic structures, Phys. Rev. A 91, 023825 (2015).
- Poo, Y. et al. Manipulating one-way space wave and its refraction by time-reversal and parity symmetry breaking. Sci. Rep. 6, 29380 (2016).
- Feng, S. Loss-induced super scattering and gain-induced absorption, Opt. Express 24, 1291-1304 (2016).
- Y.-L. Xu, et al. Experimental realization of Bloch oscillations in a parity-time synthetic silicon photonic lattice, Nat. Commun. 7, 11319 (2016).
- Savoia, S., Castaldi, G., Galdi, V., Alù, A., & Engheta, N. Tunneling of obliquely incident waves through PT-symmetric epsilon-near-zero bilayers. Phys. Rev. B 89, 085105 (2014).
- Savoia, S., Castaldi, G., Galdi, V., Alù, A., & Engheta, N. PT-symmetry-induced wave confinement and guiding in -near-zero metamaterials, Phys. Rev. B 91, 115114 (2015).
- Zhang, L., Mei, S., Huang, K., & Qiu, C.-W. Advances in full control of electromagnetic waves with metasurfaces. Adv. Opt. Mat. 4, 818-833 (2016).
- Yu, N. & Capasso, F. Flat optics with designer metasurfaces. Nat. Mat. 13, 139-150 (2014).
- Jain, A., Moitra, P., Koschny, T., Valentine, J., & Soukoulis, C. M. Electric and magnetic response in dielectric dark states for low loss subwavelength optical meta atoms. Adv. Opt. Mat. 10, 1431-1438 (2015).
- Monticone, F., Estakhri, N. M., & Alù, A. Full Control of nanoscale optical transmission with a composite metascreen. Phys. Rev. Lett. 110, 203903 (2013).
- Delvin, R. C., Khorasaninejad, M., Chen, W. T., Oh, J., & Capasso, F. Broadband high-efficiency dielectric metasurfaces for the visible spectrum. PNAS 113, 10473-10478 (2016).
- Singh, R. et al. Terahertz metamaterial with asymmetric transmission. Phys. Rev. B 80, 153104 (2009).
- Mutlu, M. & Ozbay, E. A transparent 90 polarization rotator by combining chirality and electromagnetic wave tunneling. Appl. Phys. Lett. 100, 051909 (2012).
- Serebryannikov, A. E., Mutlu, M., & Ozbay, E. Dielectric inspired scaling of polarization conversion subwavelength resonances in open ultrathin chiral structures. Appl. Phys. Lett. 107, 221907 (2015).
- Engheta, N. Pursuing Near-Zero Response. Science 340, 286-287 (2013).