Scattering properties of PT-symmetric chiral metamaterials

The combination of gain and loss in optical systems that respect parity-time (PT)-symmetry has pointed recently to a variety of novel optical phenomena and possibilities. Many of them can be realized by combining the PT-symmetry concepts with metamaterials. Here we investigate the case of chiral metamaterials, showing that combination of chiral metamaterials with PT-symmetric gain-loss enables a very rich variety of phenomena and functionalities. Examining a simple one-dimensional chiral PT-symmetric system, we show that with normally incident waves the PT-symmetric and the chirality-related characteristics can be tuned independently and superimposed almost at will. On the other hand, under oblique incidence, chirality affects all the PT-related characteristics, leading also to novel and uncommon wave propagation features, such as asymmetric transmission and asymmetric optical activity and ellipticity. All these features are highly controllable both by chirality and by the angle of incidence, making PT-symmetric chiral metamaterials valuable in a large range of polarization-control-seeking applications.


Introduction
A quantum mechanical Hamiltonian is called PT-symmetric, if and only if, it commutes with the combined action of parity, ̂, and time-reversal, ̂, operators, i.e. [̂̂,̂] = 0, with the brackets denoting the commutator. Since the action of parity on the momentum and position operators, ̂ and ̂ respectively, results to ̂→ −, ̂→ −̂ , and the time-reversal to ̂→ −, ̂→̂ and → − [1][2][3][4], the requirement for PT-symmetry in a one dimensional quantum Hamiltonian results to the symmetry condition for the potential * (− ) = ( ). In other words, the real part of the complex potential must be a symmetric function of position, while the imaginary part should be antisymmetric. In fact a characteristic feature in PT-symmetric Hamiltonians is the reality of the eigenvalues below a threshold value of the potential (region where also eigenstates are PT-symmetric -usually called PT-phase), while above that threshold (in the so-called PT-broken phase) the PT-symmetry of the eigenstates is not preserved anymore and the eigenvalues become complex. At the transition point between PT and PTbroken phase, the so-called exceptional point (EP), two or more eigenvalues and the corresponding eigenvectors coincide, making EPs to be associated with exotic wave propagation behavior [5].
Although the PT-symmetry concept started from quantum mechanics, it was very quickly transferred to optics, where in recent years is being studied very actively [6][7][8][9][10][11][12][13][14][15][16]. In optical systems the PT-symmetry requirement results in symmetry conditions for the complex material parameters [6], i.e. permittivity and permeability, realizable by proper combination of loss and gain in a system. Since both gain and loss can be introduced and adjusted externally, it is much more straightforward to realize PTsymmetry in optical systems than in quantum ones, and this has been the main reason for the increasingly growing studies on optical PT-symmetric systems . These studies revealed a series of novel and uncommon optical phenomena, such as unidirectional reflectionless perfect transmission [17], unidirectional invisibility [11,15], PT-breaking transitions [19], simultaneous coherent perfect absorption (CPA) and lasing [18][19][20] and others.
Many of the above phenomena have been studied and demonstrated in scattering rather than paraxial propagation configuration. In such configurations the identification of the different PT-related phases is usually done through the eigenvalues of the scattering matrix [15,17]. In the PT-symmetric phase the scattering matrix eigenvalues are unimodular, while in the PT-broken phase they are not unimodular but they form pairs of inverse moduli. At the EP two or more eigenvalues and eigenvectors coincide.
A major part of the existing studies regarding PT-symmetry in scattering configurations has been devoted also to the combination of PT-symmetry with metamaterials [21][22][23][24][25], and mainly with zero permittivity or permeability metamaterials [22,23,25], due to their unusual properties, such as constant phase throughout the whole metamaterial volume. Among the interesting relevant effects reported are nonreciprocal modes [25] at the interface between two mu-near-zero PT-symmetric metamaterials, exponentially bound interface modes characterized by zero attenuation in epsilon-near-zero metamaterials [23], etc.
A category of metamaterials that were until recently highly unexplored under the concept of PTsymmetry is chiral metamaterials [26][27][28], which are the subject of the present work. A chiral medium is such that it cannot be superimposed with its mirror image [29,30]. This asymmetry results in a strong magneto-electric coupling, i.e. excitation of a magnetic polarization by an applied electric field and of an electric polarization by a magnetic field. As a result, the medium responds differently to left-and right-handed circularly polarized waves (LCP/RCP), which are the eigenwaves in chiral media. Among the interesting and useful effects resulting from this different response are the circular dichroism, i.e. different absorption of LCP and RCP waves, and the optical activity, i.e. the polarization rotation of an initially linearly polarized electromagnetic wave passing through the medium. (Both of these effects are due to the different refractive index (real and imaginary part) that LCP and RCP waves experience in chiral media). In chiral metamaterials the magneto-electric coupling is orders of magnitude stronger than in natural chiral media [30][31][32][33][34][35][36], because of their unit cell being much larger than atomic size; as a result, circular dichroism and optical activity are also much stronger in metamaterials, a feature which can be exploited in a variety of polarization control applications, such as polarization converters, polarization filters [35] etc.
Based on the above, by combining PT-symmetry (gain-loss and all the effects coming from their interplay, i.e. coherent perfect absorption-lasing, PT phase transition etc.) with chiral metamaterials (polarization related effects), two different worlds of propagation characteristics can be merged, thus promising unique propagation and polarization control possibilities. This paper aims to reveal part of these possibilities, some of which were also pointed out by our recent studies [26,27].
More specifically, in this work we start by showing that in chiral media, despite the lack of mirror (and thus parity) symmetry, it is still possible to achieve full PT-symmetry and thus real eigenvalues; moreover we derive the necessary material parameter conditions for full PT-symmetry. Next we exploit these conditions in scattering configurations. In particular, we consider a simple model-system composed of a chiral PT-symmetric bilayer and we explore its scattering response under illumination with circularly polarized (CP) waves. First we investigate the case of normally incident waves. We find that the eigenvalues of the scattering matrix in this case are independent of chirality, i.e. the PTsymmetric phase is the same for chiral and non-chiral systems of the same permittivity and permeability values (obeying PT-symmetry). This is an important finding, as it allows us to tune the PT-related phenomena independently from the chirality-related phenomena. Next, we expand the analysis by allowing obliquely incident CP waves and we find that, now, the eigenvalues of the scattering matrix do depend on the chirality. Additionally, we find that any deviation from normal incidence opens up a rich phase diagram, including three distinct phases: full PT-symmetric phase, mixed PT-symmetric phase and broken PT-symmetric phase. As we demonstrate, both the boundaries of these phases and the polarization features of the waves transmitted through the system can be highly controlled by tuning the chirality parameter and/or the angle of incidence. Thus, as we demonstrate here, in PT-symmetric chiral metamaterials a wide range of advanced polarization control capabilities is possible.

PT-symmetry conditions in chiral systems
In order to examine if there is a possibility to achieve PT-symmetry in chiral media and to derive the associated conditions we transform the Maxwell's curl equations, ∇ × = and ∇ × = − with the chiral constitutive relations [30,37] = 0 + ( / ) and = 0 − ( / ) into an eigenvalue problem analogous to that of the Schrödinger's equation. (In the above relations , , refer to the relative permittivity, permeability and the chirality parameter respectively, 0 , 0 are the vacuum permittivity and permeability respectively and is the vacuum speed of light). Solving the resulting equations in terms of and , the magneto-electric coupling in the constitutive relations does not allow for the problem to be cast as a simple eigenvalue problem, because the eigenvalue ω appears in the generalized Hamiltonian as well. A way to overcome the problem is to solve the equations for the fields , , instead. In that case we obtain an eigenvalue problem of the form where ̂ is a tensor-operator analogous to the Hamiltonian in the Schrödinger problem, expressed as In. Eq. (2) 1 ( ) = − /( − 2 ), 2 ( ) = 0 2 /( − 2 ), 3 ( ) = 0 2 /( − 2 ) and 4 ( ) = /( − 2 ) (for simplicity we have omitted the space dependence in , and ). By requiring ̂ to commute with the PT operator (see e.g. Ref. [21,26]), i.e. [̂̂,̂] = 0, and taking into account that ̂̂( ∇ ×) = ∇ ×, we derive the material parameter conditions as to have a PT-symmetric chiral system. These conditions read as follows: Note that the first two of the conditions (3) are the ones that have been derived also in the case of nonchiral metamaterials [21]. As has been already widely discussed though [5][6][7][8][9][10][11][12][13][14][15][16][17][18][19], the full PT-symmetry, i.e. the existence of real eigenvalues, requires that, besides the PT-symmetry of the Hamiltonian, ̂, (ensured by the condition [̂̂,̂]=0), also the eigenstates of ̂ are PT-symmetric. In chiral media, where the eigenstates are left-and right-handed circularly polarized waves (denoted below by the subscript -and + respectively), the action of parity operator, ̂, and time-reversal operator, ̂, reveals that the eigenstates cannot be fully PT-symmetric. In particular, for eigenstates of the form ± = 0 (̂±̂) ( ± − ) , where 0 is the amplitude being a real number and ± = ( √ ± )/ are the wavenumbers for RCP(+)/LCP(-) waves, the action of time-reversal operator, ̂ ( → − , → − ), in ± gives ̂± = 0 (̂∓̂) − ( ± * + ) (4) where the star (*) denotes complex conjugation. The action of parity operator, ̂, in (4) flips the space ( → − , → − , → − ), i.e. ̂̂± = − 0 (̂∓̂) ( ± * − ) (5) Employing the conditions (3) to Eq. (5), to associate the wavenumbers ± with their PT-symmetric counterparts, we conclude that Thus, the action of ̂̂ in our case does not leave the eigenwaves unchanged; it rather transforms an initial RCP wave to LCP and vice versa -except for a phase factor. (Note that using the conditions (3) the above PT transformations lead to the same conclusion also for the fields H, B and D.) Nevertheless, despite the lack of PT-symmetry in the eigenwaves, there is still the possibility for real eigenvalues (PTsymmetric phase). This possibility is ensured by the degeneracy of RCP/LCP eigenwaves with respect to the frequency, as can be easily seen by applying the condition (3) into Eq. (1) [26]. Hence, our analysis demonstrates the possibility to obtain real eigenvalues in PT-symmetric chiral systems and provides the associated conditions. However, as numerical results in scattering systems have shown, the conditions (3) for the chirality parameter are not the only ones that lead to fully PT-symmetric phase, at least in scattering configurations [26]. There, there are cases beyond condition (3) for κ which yield similar features and similarly useful properties [27].

Chiral systems in scattering configurations under normal incidence: The scattering matrix
Having derived the necessary conditions for PT-symmetry in chiral systems [Eq. (3)], in the rest of the paper we turn our attention to scattering configurations, applying those conditions and investigating the scattering features and response of the resulting systems. Our model-system is a homogeneous onedimensional (1D) chiral bilayer (of thickness = 2 along z and infinite extent along the , directions), as shown in Fig. 1. The complex material parameters , , and of the two media involved in our system are the relative permittivity, relative permeability and chirality parameter, respectively. = * , = − * , with ε, μ, κ being the relative permittivity, permeability and chirality parameter. The amplitudes of the incident ( ± , ± ) and scattered ( ± , ± ) waves are also shown, where the subscripts +/-account for RCP/LCP waves.
We consider in our system normal incidence with four input and four output ports, each one allowing either left circularly polarized (LCP) or right circularly polarized (RCP) waves as illustrated in Fig. 1. We start with the most general form for the electric field outside the PT-symmetric structure, as where ± , ± and ± , ± are amplitudes of the ingoing and outgoing, respectively, LCP (-) and RCP (+) waves (as observed from the source point), ê + = (̂+̂)/√2 and ê − = (̂−̂)/√2, and 0 = / . Similar expression can be used for the electric field inside the chiral bilayer, with 0 replaced by the relevant chiral wavenumbers ± = ( √ ± )/ [30]. Due to the two possible circular polarizations at each of the two sides of the structure, the system should be described by a 4x4 scattering matrix 0 , defined by (see Fig. 1) where the last equation shows the S0 matrix as is simplified in our case. In Eq. To investigate when the system lies in PT-symmetric or PT-broken phase and to identify the position of exceptional points the first step is to evaluate the eigenvalues and eigenvectors of the scattering matrix 0 [17,38,39]. As was mentioned also in the Introduction, in the PT-symmetric phase the eigenvalues of the scattering matrix are unimodular, while in the PT-broken phase they are not. Therefore by calculating the scattering matrix eigenvalues one can identify the different possible PTrelated phases and the EP. The eigenproblem solution for our scattering matrix 0 yields two degenerate pairs of eigenvalues, σ, which are given by The corresponding eigenvectors are given by From Eqs. (8), (10) and Fig. 1 one can observe that the eigenvectors correspond to pure circularly polarized waves (either incoming or outgoing) of opposite handedness at the two opposite sides of the slab.
Note that depending on the arrangement of the input and output ports in the column matrices of Eq. (8), we can build our scattering matrix formalism in several ways, which characterize different physical processes. Here, we have chosen the case related with the position of exceptional points (EPs), in which the system passes from PT-symmetric phase to PT-broken phase [17][18][19]39] (similarly to nonchiral PT symmetric systems) [19]. In this case, the scattering matrix 0 satisfies the relation [17,19], which is a fundamental condition obeyed by PT-symmetry.
Although for PT-symmetry we should impose conditions (3) to the material parameters of our system, to obtain more general expressions we start our analysis with arbitrary parameters , , , , , . For arbitrary material parameters the reflection and transmission coefficients for circularly polarized waves are expressed as and ± = 0 (√ ± ), = or , are the wavevectors while = √ 0 / 0 , = or , are the wave impedances in the two chiral layers. Investigating further the above equations, one can observe that the reflection coefficients are totally independent of the chirality parameters of the two media, since all the terms , , ( ) , ( ) , ( ) , ( ) are chirality independent (note that they depend only on the impedances , and terms of the form ( + + − ) = 2 0√ and ( + + − ) = 2 0√ , all independent of chirality). Moreover, the product ++ −− is also independent of chirality. (Note that ++ = − ℎ 0 ( + ) and is the transmission for the system of the same permittivities and permeabilties but with zero chirality for both media.) Therefore the eigenvalues of the scattering matrix are independent of chirality, demonstrating that the different PT-symmetry-related phases of our system and the exceptional point are totally chirality independent. This is an important feature, which will be demonstrated and discussed further in the next sections.

Generalized unitarity relations for chiral PT-symmetric systems
By exploiting the PT-symmetry properties of our system following the procedures of [17], we can formulate the generalized unitarity relation for the four-port case. This can be done using the transfer matrix formalism, which is commonly used in electromagnetic theory [17,18]. We start from the general 4x4 scattering matrix configuration (Eq. 8), expressing the elements of the transfer matrix in terms of the elements of the scattering matrix 0 as Therefore, the transfer matrix can be written as It is straightforward to show that the transfer matrix satisfies the quite general condition det[ ] = 1 [17,18] and, additionally, the symmetry relation = ( −1 ) * which is a direct consequence of PTsymmetry. By taking into account the above two conditions, resulting to * = (with I the identity 4x4 matrix), we obtain the generalized unitarity relation for a general chiral PT-symmetric system, which is expressed as where ++ ≡ | ++ | 2 , −− ≡ | −− | 2 and ( ) ≡ | ( ) | 2 , ( ) ≡ | ( ) | 2 are the transmission and reflection power coefficients, respectively. It should be noted that the relation (15) is valid for all 1D chiral PTsymmetric structures and it holds in both PT-symmetric and PT-broken phase. It does not correspond though to flux conservation, which can occur only under strict conditions regarding the incident wave configurations.
Since, as was pointed out earlier, both the reflection coefficients and the product ++ −− are independent of chirality, all the terms of Eq. (15) are chirality-independent, although the coefficients ++ and −− alone are chirality-dependent. This shows that if we add chirality in a nonchiral PT-symmetric bilayer which is at an anisotropic transmission resonance (ATR) point (i.e. = 1 and one of the ( ) , ( ) is zero), which corresponds to unidirectional flux conservation, we will always be at the condition ++ =1/ −− . In that case, though, a flux-conserving process for waves incident from one side is realized only if ++ = −− = 1, a condition incompatible with the PT-symmetry requirement for κ (for non-zero Im(κ)) and normal incidence [see Eqs. (11) and paragraph after them].

PT-symmetric chiral bilayer under normal incidence: Different PT phases and scattering characteristics
In non-Hermitian systems the amplitude of the scattering matrix eigenvalues is either below (decay) or above (growth) unity. In the case of PT-symmetry, where gain and loss are in balance, the eigenvalues of the scattering matrix can be unimodular and are in fact unimodular below the exceptional point; thus their calculation offers the possibility to identify the exceptional point and the different PT-related phases. To investigate this possibility and its associated effects in the case of chiral systems, we investigated the system of Fig. 1 9)) as a function of frequency are plotted in Fig. 2(e), demonstrating the existence of the two different phases, the PT-symmetric one characterized by unimodular eigenvalues and the PT-broken one, characterized by eigenvalues of inverse moduli. At the exceptional point (at / = 22.25) all eigenvalues coincide. As has been discussed already in Section 2, the eigenvalues in our case are independent of chirality, since both the reflection coefficients and the transmission product ++ −− are chirality independent. This is also verified numerically in Fig. 2, where besides the eigenvalues of the chiral PT system we have also plotted the eigenvalues of a system with the above PT-symmetric permittivity and permeability values but with κ=0 in both media (see Fig. 2(a)) and for a system without any symmetry in κ, i.e. (− ) = −0.06 − 0.02 , ( ) = −0.04 − 0.01 (see Fig. 2(i)). For all three systems the scattering matrix eigenvalues shown in Fig. 2 are identical. All features are plotted as a function of the dimensionless frequency / , with c the vacuum speed of light and L=2d the total system length/thickness.
For the above-mentioned three systems we have calculated also the amplitudes of the scattering matrix eigenvectors (second row of Fig. 2) and the transmission (T) and reflection (R) power coefficients for circularly polarized incident waves (third row of Fig. 2). From the eigenvector plots (depicting the non-zero, non-unity components of Eqs. (10)) it can be observed that while in the case of non-chiral PTsymmetric systems (Fig. 2(b)) the ratio of the two non-vanishing components for each eigenvector is unimodular below the exceptional point, in chiral PT-symmetric structures (Fig. 2(f)) it is not. Considering a system excitation configuration of the form of an eigenvector, the corresponding scattered waves, which are pure circularly polarized waves of opposite handedness at the left and right side of the slab, above the exceptional point will be either exponentially growing or attenuating (depending on the incident wave configuration). At particular frequencies above EP these two modes (growing-attenuating) are expected to give simultaneous CPA and lasing of circularly polarized waves, analogous with the CPA-lasing modes of the non-chiral PT-symmetric systems [18][19][. The particular circular polarization favored in our foreseen CPA-laser modes depends mainly on the sign of the Im(κ) of the two chiral slabs. We have to note here that the above described eigenvector behavior seem to survive qualitatively even if the chiral system is beyond PT-symmetry (see Fig. 2(j)).
Regarding the transmission and reflection coefficients, for our isotropic chiral medium, as expected, ++ and −− do not depend on the side of incidence (since the system is reciprocal), while a similar condition does not hold for the reflection coefficients, as can be seen in Fig. 2(g). Comparing Fig. 2(g) with Figs. 2(c) and 2(k) (i.e. the case of non-chiral PT-symmetric structure and the case beyond PT-symmetry in κ, respectively), we can see that the chirality parameter has strong influence on the transmission coefficient and none on the reflection. Note that the chirality independence of the reflection coefficient was also observed in the corresponding analytical expressions (even for non-PT systems -see Section 2), while from the same expressions (see Eqs. (11), (12) and paragraph after them) one can conclude that the transmission (power) coefficient depends on the chirality parameter exponentially (in particular T depends on the Im(κΑ+κΒ) -see Fig. 1). Such a dependence, although having strong influence on the transmission magnitude, does not affect the frequency position of the zeros and the resonances of the transmission, something observable also from Fig. 2.
In the rest of the current subsection we investigate the impact of our system on the polarization of an incoming wave. For that, we examine first the optical activity, θ, of our structure, which for absence of cross-polarized transmission terms (e.g. +− ) is given by = 0.5[Arg( ++ ) − Arg( −− )]. Taking into account Eqs. (11a) and (11b), we can find that for a general chiral bilayer as the one of Fig. 1 (i.e. prior application of PT-symmetry conditions) Due to the symmetries imposed by applying PT [Eq. (3)], the real part of the chirality changes sign across the bilayer. Hence, the polarization rotation occurring in the gain slab is subsequently canceled out when the wave passes through the loss slab, resulting to zero optical rotation, as is demonstrated also in Fig. 2(h). Next we consider the other important polarization-related property, which is the circular dichroism (CD). A quantity directly related with the CD response is the transmitted wave ellipticty, η, for a linearly polarized incident wave, which is given by (The second r.h.s. of Eq. (17) has been obtained employing equations (11), concerning general material parameters in the two slabs). In a chiral PT-symmetric system, the imaginary part of chirality preserves its sign across the system implying a very strong circular dichroism (CD), resulting in a large degree of the transmitted wave ellipticity, as shown/confirmed also in Fig. 2(h).
Summarizing, as shown in Fig. 2 and discussed also in connection with analytical calculations, the PT-symmetric features of a chiral PT-symmetric bilayer (i.e. different PT-phases, exceptional point, structure resonances) are totally independent of chirality. On the other hand, the chirality strongly affects the polarization state of the wave passing through the bilayer (through its effect to its ellipticity). This shows that if one has the ability to control separately the system permittivity/permeability and the system chirality (as is to a large extent possible in chiral metamaterials) one can combine or superimpose almost at will the PT-related and the chirality-related features, achieving fascinating or important in applications effects, e.g. CPA-lasing for circularly polarized waves.

PT-symmetric chiral bilayer under oblique incidence: Controlling the PT-symmetry phase
As has been already discussed, in the case of waves normally incident on the structure of Fig. 1 the position of EP, which characterizes the transition from the PT-symmetric to the PT-broken phase, is totally independent of chirality. As we discuss in this section, the situation changes for oblique incidence, where chirality strongly affects the different PT-related phases and the EPs.
Since, as was shown recently [40], transverse electric (TE) and transverse magnetic (TM) polarized waves are associated with different exceptional points, it follows as a result that for circularly polarized light under oblique incidence a mixed phase is realizable. Therefore, it is interesting to examine the different PT-related phases and the phase transitions for CP light obliquely incident on our chiral bilayer structure. In this case, the transmission and reflection coefficients be considered [27]. By numerically calculating the eigenvalues of this scattering matrix the attainable PT-related phases of the bilayer can be identified.
As an example, we investigate here a system with permittivity and permeability values the same as in Fig. 2  To analyze further the impact of chirality on the different PT-related phases of a chiral bilayer we investigate the different PT-related phases for a system with permittivities and permeabilities as in Fig.  3 by scanning the PT-obeying chirality parameter (both real and imaginary part) at a fixed frequency, / = 15.5, and incidence angle = 45 . The results for the eigenvalues of the scattering matrix as a function of chirality are illustrated in Fig. 4(a). In Fig. 4(a) we can see that as we increase the chirality the system passes from PT-symmetric phase to mixed PT-symmetric phase (light grey) and, with further increase of the chirality, to PT-broken phase, indicating a possibility of a full control of exceptional points (EPs) and the associated PΤ phases by changing the chirality. Surprisingly, it is possible to achieve analogous control by changing only the real or only the imaginary part of the chirality, as illustrated in Figs. 4(b) and 4(c), respectively. All the above reveal a very rich behavior and possibility for different phases and phase re-entries as one changes the system chirality. Although for oblique incidence the PT-related features of a chiral bilayer can be highly controlled by chirality, however chirality can not offer an external dynamic control. Such a control can be offered by the angle of incidence as we show in Fig. 5. In Fig. 5 we plot the phase diagrams of the same systems as in Figs. 3(a) (see panels (a), (b) of Fig. 5) and 3(b) (see panels (c) and (d)) as we change the incidence angle. As can be observed there, the frequency positions of the exceptional points and, consequently, the frequency extent of the different PT-phases are highly dependent on the incidence angle. Moreover, incident angle controllable mixed phases and phase re-entries (see panel (d)) can be achieved. The possibility of PT-feature tuning by the incidence angle offers a great and practical way for a dynamic control of chiral PT-symmetric systems. The vertical dashed lines corresponds to the exceptional points of plots (a) and (c).
As was mentioned already, among the interesting and particularly useful characteristics associated with chiral media are the optical activity and the circular dichroism (CD), which for our double-layer chiral slab and for normal incidence depend exclusively on the chirality parameter of the two layers (see Eqs. (16), 17)). In the case of oblique incidence, where in the chiral bilayer there are cross-polarized transmission terms even for circularly polarized incident waves, the optical activity and the transmitted wave ellipticity (a measure of the CD response) cannot be obtained anymore by the simple relations (16) and (17). They can be calculated though through the Stokes parameters [27,29], which describe completely the polarization state of a wave. The four Stokes parameters for the transmitted wave in our case are defined by 0 = ⊥ ⊥ * + ∥ ∥ * , 1 = ⊥ ⊥ * − ∥ ∥ * , 2 = 2Re[ ⊥ ∥ * ] and 3 = 2Im[ ⊥ ∥ * ], where the star (*) denotes complex conjugation, the subscript ∥ indicates the transmitted electric field component that lies on the plane of incidence while the subscript ⊥ indicates the perpendicular component. Through Stokes parameters the optical activity is given by and the ellipticity by To be able to calculate θ and η employing the above formulas and to have a full picture for our system potential for polarization manipulation, one needs to calculate the scattering coefficients (transmission and reflection) for linearly polarized incident waves. These coefficients can be directly obtained from the circularly polarized reflection and transmission data according to the equations [30] In Eqs. (20), (21), as in the Stokes parameters, the subscripts ∥ and ⊥ indicate the components parallel and perpendicular to the plane of incidence, respectively, while, as in the circular polarization case, the first subscript refers to the transmitted (or reflected) component and the second to the incident one. In Fig. 6 we show the transmitted and reflected power coefficients (T=|t| 2 , R=|r| 2 respectively) for linearly polarized waves as well as the corresponding optical activity and ellipticity for the chiral PTsymmetric system examined in Fig. 2 under oblique incidence. In particular, we present those quantities at / = 15.5 as a function of chirality (left two columns) and as a function of the incidence angle (right two columns) keeping all the other parameters constant. Our results reveal a very rich behavior of propagation characteristics including asymmetric transmission, and asymmetric ellipticity and optical activity. Moreover, all these asymmetric effects are angle dependent, offering an additional degree of freedom for controlling the scattering and polarization properties of electromagnetic waves.

Conclusions
We investigated the effects of combining in the same system PT-symmetry (gain-loss media) and chirality. Starting by showing that such a combination is possible (despite the mirror-asymmetry of the chiral media) and by deriving the necessary conditions, we applied those conditions in a simple chiral bilayer and we investigated its different PT-related phases and its scattering-and polarization-related properties. In the case of normally incident electromagnetic waves we found that the PT-symmetric features, i.e. exceptional points, PT-phases etc., are not affected at all by chirality; as a result PTsymmetric and chiral effects can be controlled simultaneously and independently. In the case of oblique incidence, though, chirality highly affects the PT-symmetry related properties. In that case we showed that a full tuning of PT phases and scattering properties of our systems can be realized both by scanning the chirality parameter and/or the angle of incidence of the incoming waves. Moreover, we showed that by combining PT-symmetry with chirality we can achieve exotic propagation and scattering functionalities, such as asymmetric (side-dependent) transmission, ellipticity and optical activity. All the above effects can be exploited in a wide range of applications, in particular in applications where advanced polarization control is required (realizing, e.g., chiral unidirectional laser-absorbers, controllable polarization isolators etc.).
Author Contributions: M.K. and S.D. conceived the idea and designed the parametric study. C.M.S., E.N.E., and M.K. supervised the project. I.K. and S.D. carried out the analytical and numerical calculations. All authors discussed the results and commented on the manuscript, which was written by I. K..