Broadband mixing of PT\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathscr{P}}{\mathscr{T}}$$\end{document}-symmetric and PT\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathscr{P}}{\mathscr{T}}$$\end{document}-broken phases in photonic heterostructures with a one-dimensional loss/gain bilayer

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 PT\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathscr{P}}{\mathscr{T}}$$\end{document}-symmetric and PT\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathscr{P}}{\mathscr{T}}$$\end{document}-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 PT\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathscr{P}}{\mathscr{T}}$$\end{document}-symmetric and PT\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathscr{P}}{\mathscr{T}}$$\end{document}-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.

and cavities with PT -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 PT -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 PT -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 PT -related properties, which do not need two-dimensional loss/gain components. Indeed, although the principal possibility of co-existence of PT -symmetric and PT -broken phases is a known feature for the case of physical dimensions higher than one (i.e., when more than two input and output channels may exist) 11 , there is high demand in simpler structures that would enable similar and easily realizable scenarios.
In this report, we theoretically demonstrate the co-existence of PT -symmetric and PT -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 the function of combining gain and loss are separated in space, so far as they are performed by the different components of the entire heterostructure. 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 PT -symmetric and PT -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 PT -symmetric and PT -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, and 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 M is the loss/gain medium that satisfies the symmetry condition ε ε = − ⁎ z z ( ) ( ). This condition is necessary but not sufficient for PT -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 L 0 and R 0 change the polarization state for both reflected and transmitted waves. The second-type PCCs (denoted by L 1 and R 1 ) only change the transmitted waves' polarization, whereas the third-type PCCs (denoted L 2 and R 2 ) only change the reflected waves' polarization. The properties of PCCs of the three types There are four channels in total, which contain two input and two output channels corresponding to the right and left circular polarizations (CPs); they are described by eight coefficients for the scattering amplitudes, ± a , ± b , ± c , ± d . The middle region labeled by M consists of gain and loss components satisfying ε ε = − ⁎ z z ( ) ( ). Different types of PCCs (labeled as L i and R i ) can be added at the end-faces of region M to extend a variety of achievable scenarios. (b) Three types of PCCs are considered: L 0 component allows changing the polarization of both the transmitted and reflected waves at left-side incidence (perfect polarization conversion from right/left CP to left/right CP is assumed), whereas L 1 and L 2 do so only for polarization of the transmitted and only for polarization of the reflected waves, respectively. For right-side incidence, the components R 0 , R 1 and R 2 show the same properties as their left-side counterparts, i.e., L 0 , L 1 , and L 2 .
Scientific RepoRts | 7: 15504 | DOI:10.1038/s41598-017-14982-w 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 z-direction that consists of right (+) and left (−) polarized waves. For the left-side half-space, < − z L/2, and normal incidence, it is given by To obtain E for the right-side half-space, i.e., at > z L/2, we replace → a c and → b d in equation (1). Similar expressions can be used for the bilayer region.
By exploiting the symmetry properties of the structure, we can cast the generalized conservation relation 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 11 . It is connected with another key relation known as the generalized unitarity relation. Starting from the condition ( ), 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: are the reflectances at left-side and right-side illumination, with no restrictions imposed on them; t is transmission coefficient, r L and r R are reflection coefficients at left-side and right-side illumination. In line with the -matrix formalism, all 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 t, r L , and r R .

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 PT -symmetric and PT -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 mixed phase 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 M, 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 L i and R i at the end-faces.
Case 1: The simplest configuration to access this case is a bilayer enabling PT -symmetry (like the region M) 10,11 . For the configuration M, -matrix has two eigenvalues. Let us generalize it by formally adding the polarization related degree of freedom to the system to obtain the X 4 4 -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, yields the same set of eigenvalues twice, with the two eigenvalues in a set, which are given by 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 (L 0 and R 0 ). The configuration that we have now is L MR 0 0 , 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 equation (3) are coupled due to the added PCCs. After some algebra, we obtain The four eigenvalues for this configuration are given as Scientific RepoRts | 7: 15504 | DOI:10.1038/s41598-017-14982-w Thus, configurations M and L MR 0 0 show the same moduli for the eigenvalues and the same basic physics. It is well-known that a unitary transformation 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 PT -symmetric or PT -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 that PT -symmetry is spontaneously broken at the critical point. We showed in our study that they include the configurations L MR and MR 1 do not yield the mixed phase despite the fact that polarization conversion is possible at one of the end-faces. Since the components L 1 and R 1 only change the polarization of transmitted waves, the ability of polarization conversion of the reflected wave is expected to be the necessary (but not sufficient) condition of existence of the mixed phase. It is noticeable that the same -matrix 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 ω ω < c , all four eigenvalues are unimodular (log10 λ = 0 m , m = 1, 2, 3,4), being in PT -symmetric phase and, hence, are flux conserving. At ω ω > c , 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 ωL c / . 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 M. However, it cannot ensure the phase mixing, as follows from the study of Case 1. Thus, let us first add a PCC, specifically L 0 to the left end-face, i.e., we now have configuration L M 0 . 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 changed while the reflected wave is assumed to retain its polarization. These properties are described by the following -matrix: 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 equation (7) are given by As one set of eigenvalues (λ 1,2 ) preserves the symmetry and unimodularity at ω ω < c , the other set of eigenvalues (λ 3,4 ) is translated into a symmetry-broken phase in a wide ω-range, even near ω = L c / 0. Once ω ω = c 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 ω ω > c . So, ω ω = c is the boundary between the mixing phase and the all-broken phase cases. Figures 3 and 4 show the maps of the magnitudes of eigenvalues λ m , = m 1, 2, 3, 4, in log10 scale, which are obtained for Case 2 from equations (8a) and (8b), in wide ranges of variation in ε Re( ) and ω. It is clearly seen that the symmetry of eigenvalues λ 3,4 is broken even at small frequencies. Thus, there is a large region in ω ε L c Re ( / , ( )) -plane, where PT -symmetric and PT -broken sets of eigenvalues may co-exist. Starting from ω ω = c , the set of eigenvalues λ 1,2 also experiences spontaneous symmetry breaking, so the mixing does not exist anymore.
As shown in Fig. 3(a,b), the value of ω ω = c and the extent to which the symmetry is broken strongly depend on ε Re( ). For instance, we obtain ω = min L c ( / ) 5 c in the vicinity of ε = . Re( ) 0 5. For the set of eigenvalues λ 3,4 , a strong deviation from the unimodular case is observed in Fig. 3(c,d) even at very small values of ε Re( ). In particular, a strong anomaly of λ 3,4 occurs nearly at ε < < . Re 0 ( ) 0 5 and ω < < L c 2 / 4. Hence, PT -symmetric and PT -broken eigenvalues may co-exist also in ENZ regime, i.e., in the close vicinity of ε = Re( ) 0. It is noteworthy that the behaviours of eigenvalues in Case 1 and Case 2 at ε < < Re 0 ( ) 1 are very different. In Fig. 4, the boundary between the regions with the mixed phase for eigenvalues (at smaller ωL c / ) and with all symmetry-broken eigenvalues (at larger ωL c / ) is clearly seen. Its location can be controlled by variations in ε Re( ), while ε Im( ) is fixed. A detailed investigation of these scenarios will be a subject of our future research. Figures 5 and 6 display the magnitudes for eigenvalues λ 1 , λ 2 , λ 3 and λ 4 in log10 scale at varying ωL c / for ε = Re( ) 1 and ε = Re( ) 2, respectively, for two different values of ε Im( ). The location of the critical frequency for λ 1 and λ 2 , ω ω = c , and, hence, width of the ω-range, in which the mixed phase is achieved, are strongly affected by ε Im( ). As the value of ε Im( ) increases, the critical frequency is redshifted, whereas an increase of ε Re( ) results in blueshift of ω c . At the same time, ε Re( ) moderately affects the width of mixing range in the considered case. On the other hand, ε Re( ) can strongly affect the magnitudes of λ 3 and λ 4 . PT -symmetric and PT -broken sets of eigenvalues can co-exist even at ε = Re( ) 1, although the difference between the magnitudes is rather weak in this case. Since breaking symmetry at ω ω = c is a general property of λ 1 and λ 2 , it also occurs in the ENZ regime, e.g. at ε = .
Re( ) 0 02 (not shown). However, in this regime, the difference between λ 1 and λ 2 at ω ω > c is weak, whereas λ 3 and λ 4 may significantly deviate from the unimodular case in a wide range of ε Im( ) variation. ; ε = . Im( ) 0 2. The upper panels display λ 1 (a) and λ 2 (b) that are symmetric (unimodular) at ω ω < c ; the lower ones display eigenvalues λ 3 (c) and λ 4 (d), which are symmetry-broken at any ω. Deviation from unimodular case (here -blue color corresponding to 0 at the scale bar) indicates the extent to which the symmetry is broken. Until the value of ω c is reached, there is a wide range of ω, where symmetric (unimodular) and symmetry-broken sets of eigenvalues co-exist. After hitting ω c , the symmetric set also experiences a spontaneous symmetry breaking, and the two eigenvalue sets with broken symmetry start to overlap.
Next, let us add the component R 0 to the right end-face of the bilayer, so we now have the configuration MR 0 . Then, we obtain exactly the same physical scenario as for configuration L M 0 , while the roles of right and left side illuminations are interchanged. The -matrix is now given as: Figure 4. Same as Fig. 3 but for ε = Re( ) [1,15]; λ 1 (a), λ 2 (b), λ 3 (c), and λ 4 (d). As the ε Re( ) increases, the critical frequency ω c is blueshifted that results in a larger region of simultaneous symmetric (unimodular) and symmetry-broken sets of eigenvalues.
It has four eigenvalues, which may yield the mixed phase: Similarly to Case 1, we can obtain the sets of eigenvalues given by equations (8a)  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 PT -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 very important. It is noticeable that despite having different end-faces all 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 second PCC, so phase mixing is not achieved in some of 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 mixing of PT -symmetric and PT -broken phases for eigenvalues in the structures with four channels and with a one-dimensional bilayer. 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 bilayer and, moreover, its 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 PT -related properties. It is noteworthy that the PT -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 demonstrated 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 shown to be possible even in close vicinity of ε = Re( ) 0. 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 r L , r R and t 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.