PT-symmetric transport in non-PT-symmetric bi-layer optical arrays

We study transport properties of an array created by alternating $(a,b)$ layers with balanced loss/gain characterized by the key parameter $\gamma$. It is shown that for non-equal widths of $(a,b)$ layers, i.e., when the corresponding Hamiltonian is non-PT-symmetric, the system exhibits the scattering properties similar to those of truly PT-symmetric models provided that without loss/gain the structure presents the matched quarter stack. The inclusion of the loss/gain terms leads to an emergence of a finite number of spectral bands characterized by real values of the Bloch index. Each spectral band consists of a central region where the transmission coefficient $T_N \geq 1$, and two side regions with $T_N \leq 1$. At the borders between these regions the unidirectional reflectivity occurs. Also, the set of Fabry-Perrot resonances with $T_N=1$ are found in spite of the presence of loss/gain.


Introduction -
The possibility of deliberately intermixing loss and gain in optical structures as a means to attain physical properties that are otherwise out of reach in standard arrangements has recently been explored in a number of studies [1][2][3][4][5][6][7][8].In general, such nonconservative configurations can display surprising behavior resulting from the possible presence of exceptional points [9,10] and mode bi-orthogonality.A particular class of non-Hermitian systems that has recently attracted considerable attention is that respecting paritytime (PT ) symmetry [1][2][3][4][5].In this regard, if a Hamiltonian is PT -symmetric, then its ensuing spectrum happens to be entirely real in spite of the fact that the potential involved is complex [9,10].Once however a loss/gain parameter exceeds a critical value, this symmetry can spontaneously break and consequently the eigenvalue spectrum ceases to be real and hence enters the complex domain [9].
The availability of optical attenuation and amplification makes optics an ideal ground where PT -symmetric effects can be experimentally observed and investigated [3,4,11].In the paraxial regime, a necessary (but not sufficient) condition to impose this symmetry in photonics is that the complex refractive index distribution obeys n(x) = n * (−x) [1,2].This latter condition directly implies real part of the refractive index potential must be an even function of position while its imaginary component must be antisymmetric.
Thus far, the ramifications of parity-time symmetry on optical transport have been systematically considered in several theoretical and experimental investigations.Processes emerging from PT-symmetry, like double refraction, band merging, unidirectional invisibility, abrupt phase transitions and nonreciprocal wave propagation have been predicted and observed [1,3,[11][12][13][14].In addition this same symmetry was recently utilized to enforce single-mode behavior in laser micro-cavities as well as to judiciously control wave dynamics around exceptional points and photonic crystals [15].Finally, similar concepts are nowadays used in other fields, like for example opto-mechanics, acoustics, nonlinear optics, plasmonics, metamaterials, and imaging to mention a few [16].In view of these developments, the following question naturally arises.To which extent could a system depart from an exact parity-time-symmetry and still exhibit PT -symmetric characteristics?Recently, a possibility of synthesizing complex potentials that support entirely real spectra despite the fact that they violate the necessary condition for PT -symmetry has been discussed in Ref. [17] within the context of optical supersymmetry.
In this paper we show that it is indeed possible to observe PT -scattering behavior even in systems that strictly speaking lack this symmetry.This is explicitly demonstrated in perfectly matched multilayer non-Hermitian arrangements having unequal layer widths.In this case, the incorporation of loss/gain domains leads to a finite number of bands-all associated with a real Bloch index that is reminiscent of actual PT -symmetric lattices [12].Special attention has been paid to an emergence of unidirectional reflectivity.Our results indicate that a PT -symmetric response can be quite robust and therefore can be observed in more involved structures even in the absence of a strict PT -symmetry.
Here the parameter γ > 0 measures the strength of balanced loss and gain inside a and b layers, respectively, and ϕ = ϕ a + ϕ b is the purely real phase shift of the wave passing each of the unit cells.In the case of no loss/gain (γ = 0) the structure is known as the matched quarter stack for which the a and b layers are perfectly matched and have equal optic paths, n Within any layer an electromagnetic wave propagates along the x -direction (Fig. 1) perpendicular to the stratification according to the 1D Helmholtz equation.Its general solution inside every unit cell can be presented as a superposition of traveling forward and backward plane waves.By combining these solutions with the boundary conditions for the wave at the interfaces between a and b layers, one can obtain the unit-cell transfer matrix M relating the wave amplitudes of two adjacent unit cells (see, e.g., Ref. [18]).For our model the elements of the M -matrix read Note that det M = 1 in accordance with the general condition inherent for transfer matrices.On the other hand, the transfer matrix M (γ) has specific symmetry, where " * " stands for the complex conjugation.This symmetry can be compared with that emerging for various PT -symmetric Hamiltonians [7,[19][20][21].Clearly, it differs from the standard one, M 22 = M * 11 , M 21 = M * 12 , that originates from time-reversal symmetry.
Frequency band structure -The dispersion relation for the Bloch index µ B determined through the eigenvalues exp(±iµ B ) of the transfer M -matrix, is defined by the matrix trace, 2 cos µ B = M 11 + M 22 .Thus, one obtains This relation determines the dependence of the Bloch wave number κ = µ B /d on the frequency ω (or, on the phase shift ϕ ∝ ω) for the wave inside the scattering region.One of the important conclusions that follows from Eq. ( 3), is a finite number N band of spectral bands defined by the condition | cos µ B | < 1.This number N band depends on the parameter γ only, and for γ 1 can be estimated as In our model this determinative condition is provided not due to the PT -symmetry but merely due to the balance between loss and gain.Outside the bands, where | cos µ B | > 1, the index µ B is purely imaginary, thus creating spectral gaps.Here the waves are known as the evanescent Bloch states, attenuated on the scale of the order of |µ B | −1 .Therefore, for a sufficiently long structure, N |µ B | 1, the transmission is exponentially small.It should be noted that in the absence of loss/gain (γ = 0), there are no spectral gaps since all spectral bands are touching.
When γ exceeds the critical value γ cr = 1, there are no spectral bands since the Bloch index µ B (ω) becomes imaginary for any frequency ω.Two examples of the dependence µ B (ω) are shown in Fig. 2 for different values of γ.It is quite instructive that within any spectral band the maximal value of the Bloch index µ given γ from the interval 0 < γ < 1 the higher the band the smaller the value of µ Transmittance -With the standard transfer matrix approach we derived analytical expression for the transmittance T N of our model consisting of N unit (a, b) cells and connected to perfect c -leads by the impedance Z, Here F(−γ, ϕ) < 0 for any value of ϕ, unlike F(γ, ϕ) that can be either positive or negative as a function of ϕ.Two examples of dependence T N (ϕ) are given in Fig. 3.The analysis shows that a distinctive property of the transmittance is that every spectral band consists of the central region with T N 1 and those with T N 1, see Figs. 3 and 4. At the borders between these regions (circles in Fig. 3) the transmission is perfect, T N = 1, and this happens due to vanishing the function F(γ, ϕ).The corresponding values ϕ s (γ) for which F(γ, ϕ s ) = 0 can be compared with exceptional points emerging in the PT -symmetric models (see, e.g., Ref. [22] and references therein).However, in our case the Bloch index µ B (γ, ϕ s ) does not vanish, in contrast to truly exceptional points for which both conditions, T N = 1 and µ B = 0, are fulfilled.These conditions are known as that providing invisibility of a model in scattering process.The only situation for the invisibility to occur in our setup is when two conditions, F(γ, ϕ) = 0 and µ B (γ, ϕ) = 0, meet.As one can see in Fig. 4, the sequence of exceptional points ϕ = ϕ sp ≡ ϕ s (γ sp ) emerges for the corresponding sequence of γ = γ sp for which two curves ϕ = ϕ s (γ) and µ B (γ, ϕ) = 0 touch each other.Note that at every fixed value of γ = γ sp the corresponding exceptional point ϕ sp is located on the top of the highest band.Conclusions -We have analytically studied the model that reveals the PT -symmetric transport although its Hamiltonian, in general, is not the PT -symmetric one.Our analysis exhibits that such a situation emerges in the bi-layer optical setup which without loss/gain is known as the quarter stack, provided all layers are perfectly matched (no reflections due to boundary conditions).
In the absence of loss/gain terms, all spectral bands are touched thus creating a perfect transmission for any wave frequency.When the balanced loss and gain are alternatingly included in all layers, the spectral gaps emerge between the bands, and the total number of bands is defined by the loss/gain parameter γ only.
Inside the spectral bands the Bloch index µ B appears to be real as it happens in the known models with the PT -symmetric Hamiltonians.The analysis shows that each of the spectral bands consists of central region where the transmission coefficient is greater than one, T N 1, and two side regions with T N 1.At the borders between these regions the transmission is perfect, T N = 1, however, µ B does not typically vanish as happens in the known PT -symmetric models.These borders are of specific interest since one of the reflectances vanishes, thus leading to the so-called "unidirectional reflectivity", the effect which is important both from theoretical and experimental points of view.The ratio between left and right reflectances is specified by the loss/gain parameter γ and the wave frequency ω, being independent of the number N of cell units.It is important that a strong difference between left and right reflectances occurs in a relatively large wave-propagation region where the Bloch index µ B is real.
The analytical expressions display that inside the spectral bands the Fabry-Perot resonances emerge in spite of the presence of loss and gain.These resonances exist both in the regions with T N > 1 and T N < 1.It is also worthwhile that the invisibility (when both relations, T N = 1 and µ B = 0 hold) is realized only for a quite specific values γ sp and for the corresponding wave frequencies located in the highest spectral bands.
In conclusion, our study unexpectedly manifests that the PT -symmetric properties may occur in non-PTsymmetric systems, the fact that may be important in view of experimental realizations of matched quarter
0) b d b .As one can see, for γ = 0 and d a = d b the corresponding Hamiltonian has no the PT -symmetry.A particular PT -symmetric realization is achieved only when d a = d b with ε (0) a = ε (0) b , and µ a = µ b .
Within these bands the Bloch index µ B is real in spite of the presence of loss and gain.This property has been widely discussed in view of the nature of the PTsymmetric systems.As is known, the real-valued solutions µ B to the dispersion equation 2 cos µ B = M 11 +M 22 can exist only when the trace M 11 + M 22 of the transfer M -matrix is real-valued, i.e. if ImM 22 = −ImM 11 .

FIG. 2 .
FIG. 2. The dependence (3) of the Bloch index µB on the phase shift ϕ ∝ ω for two values of γ.Spectral bands and gaps are shown by shaded areas.

10 FIG. 3 .
FIG.3.The logarithm of TN defined by Eq. (6) as a function of ϕ for γ = 0.25 and different values N of unit cells.Shaded regions correspond to spectral bands, and circles stand for the borders between the regions with TN 1 and TN 1.

FIG. 4 .
FIG.4.Spectral bands and gaps in dependence on the loss/gain parameter γ.Full curve shows the border µB = 0 above which µB is imaginary (gaps).Below this curve µB is real.Between full and dashed curve µB is real and TN 1. Inside dashed regions µB is real and TN 1.The dependence ϕ = ϕs(γ) is depicted by the dashed curve on which TN = 1.The horizontal dash-dotted line is shown for γ = 0.25.
Apart from the specific points ϕ s (γ), the perfect transmission, T N = 1, occurs due to Fabry-Perot resonances located within spectral bands.They are defined by the condition sin(N µ B )/ sin µ B = 0 specifying the resonant valuesµ res B = mπ/N (m = 1, 2, 3, . . ., N −1) of the Bloch index µ B .Since the spectrum µ B (ω) is restricted by µ (max) B, the total number of such resonances in a given spectral band can vary from zero to 2(N − 1), depending on the band number and loss/gain parameter γ.Specifically, for a fixed γ the higher spectral band contains the smaller number of the Fabry-Perot resonances.This result is principally different from what happens in an array of bi-layers with real and positive optical parameters, where N − 1 number of the Fabry-Perot resonances emerge in any spectral band.Unidirectional reflectivity -The expressions for the 1 remain to be valid.Thus, the ratio R extremely small values of the order of 10 −7 .