RESEARCH ARTICLE

Singular PT-symmetry broken point with infinite transmittance and reflectance----a classical analytical demonstration

  • Yingxin JIANG
Expand
  • Department of Physics, Jiangsu University, Zhenjiang 212013, China

Received date: 29 Sep 2019

Accepted date: 18 Nov 2019

Published date: 15 Dec 2021

Copyright

2020 Higher Education Press and Springer-Verlag GmbH Germany, part of Springer Nature

Abstract

To demonstrate the existence of singular parity-time symmetry (PT-symmetry) broken point in optics system, we designed a one-dimensional PT symmetric structure including N unit-cell with loss and gain materials in half. We performed an analytical deduction to obtain the transmittance and reflectance of the structure basing on Maxwell’s equations. We found that with the exact structure unit-cell number and the imaginary part of refraction index, the transmittance and reflectance are both close to infinite. Such strict condition is called the singular point in this study. At the singular point position, both the transmission and reflection are direction-independent. Away from the singular point, the transmittance and reflectance become finite. In light of classical wave optics, the single unit and total structure both become the resonance units. The infinite transmittance and reflectance result from the resonance matching of single unit and total structure. In light of quantum theory, the singular point corresponds to the single eigenvalue of electromagnetic scattering matrix. The infinite transmittance and reflectance mean a huge energy transformation from pumping source to light waves. Numerical calculation and software simulation both demonstrate the result.

Cite this article

Yingxin JIANG . Singular PT-symmetry broken point with infinite transmittance and reflectance----a classical analytical demonstration[J]. Frontiers of Optoelectronics, 2021 , 14(4) : 438 -444 . DOI: 10.1007/s12200-020-0969-3

1 Introduction

In designing optical transmission system, the material refraction indexes take the leading position. Diverse distributions of refractive index have led to great progress in optics ranging from the design and fabrication of photonic crystals [1] and photonic crystal fibers [2] to the exploration of nanoplasmonics [3] and metamaterials with negative refraction effect [4]. For most cases, people consider the ideal dielectric with real refractive index. Loss is a trivial problem that is usually solved by the loss correction or compensation gain [5]. Gain is also used to overcome the loss. For long time, loss and gain have taken a back seat in optics design. However, such a situation has been changed since PT-symmetric structure was used in optics [6,7].
The PT-symmetric structure comes from quantum mechanics. In quantum mechanics, the parity operator p^ leads to x^ x^ and p^ p^, and the operator T^ leads to x^ x^, ii and p^ p^, where x^ and p^ represent the position and momentum operators, respectively. For a non-Hermitian Hamiltonian H^=p ^2 /2+V(x^) associated with a complex potential V(x^), it is parity-time symmetric provided that it commutes with the parity-time operator, i.e., H^ p^ T^ = p^ T^ H^. In this case, H^ and p^ T^ have a common set of eigenfunctions and the non-Hermitian Hamiltonians can actually possess entirely real spectra as the Hermitian Hamiltonians possess. The lack of completeness or the breaking of PT-symmetry is associated with the presence of exceptional points for the discrete spectrum or the spectral singularity for the continuous spectrum [8].
Matter wave phenomena under the framework of non-Hermitian Hamiltonians are often taken as an analogy study in optical systems with balanced gain and loss, because optics can provide a productive test bed where the notions of parity time symmetry can be experimentally explored. In optics, PT-symmetry demands that the complex refractive index obeys the condition n(r) = n*(−r), in other words the real part of the refractive index should be an even function of position, whereas the imaginary part must be odd. In these optical systems, gain and loss have played the important role in determining optical properties. Many intriguing physical phenomena with PT-symmetric systems have been predicted and observed, such as loss induced transparency [9], power oscillations and nonreciprocity of light propagation [7], coherent perfect absorber-lasers [1012], optical switches [13], on-chip optical isolation [1416]. The one-dimensional sinusoidal [17,18] and layered PT symmetric complex crystals [19,20] with balanced gain/loss modulations may be highlighted since they show one-way enhanced reflection, while remains total transmission (unidirectional invisibility). The unique phenomena occur at the exception point or PT-symmetry breaking threshold. At the special position, the sum of reflectance and transmittance are larger than unit, which means that the extra energy is obtained from pumping light to the gain regions of the structure. Another important PT-symmetry phenomenon is that the spectral singularities in complex scattering potentials are achieved through an electromagnetic waveguide [8]. At the spectral singularities unbelievable infinite reflectance and transmittance occur. Other researchers also pointed out the existence of novel singular points in the broken-symmetry phase [1012,21,22]. The singular points are associated to the single eigenvector and eigenvalue of electromagnetic scattering matrix (S matrix) [21]. Although the concept of spectral singularities has been introduced into optics system, the analyzed method is also based on quantum theory [8]. In this paper. we demonstrate the existence of the singular broken PT-symmetry points through a classical analytical method from Maxwell equations. Thus our study can provide a more direct demonstration about the singular point. Different from the structures including one unit-cell with loss and gain materials in half embedded in the rectangular waveguide [1012,21,22], the structure in this study is a one-dimensional model including N unit-cell with loss and gain materials in half. We verify that the singular broken PT-symmetry point has infinite reflectance and transmittance. In this paper. we find that the direct physical reason for the singular point is the exact match between the unit resonance and the system resonance. We also find the singular point just corresponds to the single eigenvalue of the electromagnetic S matrix. The results are also demonstrated through numerical calculation and software simulation.

2 Model and analytical deduction

Optical PT-symmetry in one-dimensional structure requires the complex refractive index distributed in x direction in the form of n(x) = n*(−x). In this paper we construct such PT symmetric structure by placing two layers A and B with the same thickness d alternately in the x-axis. The refraction indexes for layers A and B are chosen as nA = n0 + iρ and nB = n0−iρ, respectively, to satisfy the PT-symmetry condition. The structure is schematized in Fig. 1.
Fig.1 Schematic diagram of PT symmetric layered structure

Full size|PPT slide

For such structure with N period units, we can obtain its transmission properties through strict analytical deduction. For layered structures, the transmittance and reflectance are dependent on the transfer matrixes [23]. The matrixes associated with layers A and B are
mA= [cos δA i/nAsinδA inAsinδA cosδA],
and
mB =[cosδ Bi/n Bsin δ Bin Bsinδ Bcos δB],
where δA =nAk0d=n0 k0d+iρk0d and δB=nBk0d=n0 k0diρk0d denote the phase change for normal incidence in layers A and B, respectively. For left incidence (from left to right), the matrix for a unit-cell is
m= m A mB = [m 11 m12 m21 m22]=[cos δAcosδB nB/nA sinδA sinδBi/nB cosδA sinδB i/ nAsin δAcos δB i nBcos δAsin δBinAsinδAcosδB cosδA cosδB n A/ nBsin δAsin δB ].
For the whole periodic structure with N unit-cells, we have the total matrix M
M=mN=[ M11M12M 21 M22]=[m11 uN1 (a)u N2(a ) m12u N1(a ) m21uN1(a) m22uN1(a) uN2 (a)],
where a=(m11+m 22)/2and uN are the Chebyshev Polynomials of the second kind from [23]
uN(a)= sin[(N +1)cos 1(a)] 1a2.
For background material with refraction index n0, the transmittance and reflectance of left incidence are respectively
Tleft=| 2n0( M11+M 22)n0 +(M 12n02+M21)|,
Rleft=| (M 11 M22)n0+( M12 n02 M21)( M11 +M22)n0+ (M12n 02+M 21)|.
To obtain T leftand Rleft, we first calculate the important parameter a. For simplicity, we denote δA =δ+iΔand δB=δiΔ, in which δ=n0k0d and Δ=ρk0d. Thus we have
a=(m11 +m22)/2=cosδA cosδB 12 nB nA+ nAnBsinδA sinδB = e i(δ+iΔ)+ei(δ+iΔ)2 ei(δiΔ )+ ei(δiΔ )2 12 n02ρ2n 02+ ρ2 ei(δ+iΔ)ei(δ+iΔ)2i ei(δi Δ)ei(δi Δ)2 i =cos(2δ+κ)=cos( α),
where α =2δ +κ, and κis a small quantity that is used to correct the approximation of n02ρ2n02+ρ2 =1 for n0>>ρ used in the conduction of Eq. (8). Here, we make the first constraint condition
2δ=jπ(j=1,2,3 ) .
Fromδ= n0 k0d, we have
d =jλ/(4 n0)(j=1 ,2,3 )
where λ is the incidence wavelength. Thus we have
u N(a)=sin[( N+1)α]sinα,
and
M11= m11 sin(Nα)sinαsin[( N1 )α]sinα,
M12= m12 sin(Nα)sinα,
M21= m21 sin(Nα)sinα,
M22= m22 sin(Nα)sinαsin[( N1 )α]sinα.
Next basing onα =2δ +κ=jπ+κ, we make following operations.
M11+ M22=(m11 +m22)sin (Nα)sin α 2 sin[(N 1)α]sinα= 2cosα sin(Nα)sinα2 sin[ (N1)α]sinα=2sinα[cosαsin(N α) sin(Nαα )]=2cos(Nα) =2cos(Nj π+Nκ).
For large N, we apply the second constraint condition
Nκ=(N 1+ 1/2)π ( N1=0 ,±1,± 2) .
That leads to
M11+ M22=0,
M11M22 = ±( m11m22 )sinα =± 4iρ n0sin δAsin δB(n02+ρ2)sinα ,
M12=± m12sinα ,
M21=± m21sinα ,
M12n 02= ± m12sinα n02=±[i n 02 nBcosδ Asinδ Bi n 02 nAsinδ Acosδ B]/sin α =±[ inAcosδAsinδB inBsinδAcosδB] n02n0 2+ρ2/sinα=± [i nAcos δAsin δBi nBsin δAcos δB]/ sinα.
In Eq. (22), we have made the only approximation in this study for the large contrast of n0 andρ
n02n0 2+ρ21.
Under this approximation, we further have
M12n 02+ M21=±[icosδA sinδB (n A+nB)isinδ Acosδ B(n A+ nB)]/sinα =i 2n0(cosδA sinδB +sinδAcosδB)/sinα,
M12n 02 M21=±[icosδA sinδB (n AnB)icosδ Bsinδ A(n A nB)]/sinα =±2 ρ(cosδAsinδB+cosδ Bsinδ A)/sinα ,
cos δ Asin δ B+ cosδB sinδA= ei(δ+iΔ )+ ei(δ+iΔ)2ei(δiΔ)ei(δi Δ)2 i+ e i(δiΔ )+ ei(δiΔ )2 ei(δ+iΔ)ei(δ+iΔ)2 i=sin (2δ).
Taking Eqs. (18), (19), (24) and (25) into Eqs. (6) and (7), we obtain
Tleft= |2n0 (M 11+ M22)n0+( M12 n02+ M21)| =| 2n 0sin α0 2in0sin (2δ )| =| sinαsin( 2δ)|,
Rleft=| (M11M22) n0+(M12 n02M21)(M11+ M22 )n0+( M12 n02+ M21)| =| 2ρn 02sinδ Asinδ B (n 02+ρ 2)±iρsin(2δ)sin α±n0sin(2δ)|.
From Eq. (9), we have sin (2δ) = 0. Sincesin α, sin δA and sinδB are all finite, T leftand Rleftare both close to infinite. The infinite values are based on a special status. Thus Eqs. (9) and (28) state the resonance condition and the resonance amplification effect, respectively. We will demonstrate in the next section that it is a special PT-symmetry broken point. However, we should notice that the infinite T leftand Rleftare based on the constraint conditions of Eqs. (9) and (17). Such conditions lead that such PT-symmetry broken points should be singular and only dependent on the two special structure parameters N and ρ if Eq. (9) holds.
Up to here, one may ask whether the infinite T leftand Rleft occur for only left incidence. In this study, we can demonstrate for right incidence (from right to left) infinite transmittance and reflectance occur at the same parameters N and ρ. In fact, for right incidence the matrix of a unit-cell is
m=mBmA=[ m11m12m 21m 22]=[ cosδB cosδ A nA/nB sinδB sinδ A i/nAcosδ BsinδA i/n BcosδAsin δB inB cosδ Asin δB inA cosδ Bsin δAcosδBcos δAnB/n A sinδ Bsin δA].
The total matrix M is
m=m N=[ m11m12m 21m 22]=[ m11 uN1( a) uN 2( a )m12u N1 (a )m 21uN 1(a) m22 uN1(a)u N2(a )].
From Eq. (30), we easily find that a= (m 11+m 22)/2=(m 11+ m22)/2=a. In this condition, we can demonstrate M 11+M 22=M 11+ M22, M11M22=(M11M22), M12'=M12 and M21'=M21. Comparing Eqs. (6) and (7) with
Tleft=| 2n0( M 11+M 22) n0 +( M12 n02+ M 21)|,
Rleft= |( M11 M22)n 0+( M 12n02 M21 )( M11+ M22)n 0+ (M12 n02+ M21)|,
we find T right is always equal to T left, thus the transmission is always reciprocal. However, without the constraint condition of Eq. (9), M12n 02M21 is not equal to zero. Because of M12 M22= (M11 M22), Rright is not equal to Rleft. Therefore, only at the isolated PT-symmetry broken point with M12n 02M21 = 0, the reflection is reciprocal.

3 Demonstration by numerical calculations and simulation

Here we consider a specific example. Taking arbitrary values of λ = 150 nm and n0 = 3, the value of d is taken as d=5λ/(4 n0) = 62.5 nm according to Eq. (10). Thus T left and Rleft are only dependent on the two parameters N and ρ. Based on Eqs. (6) and (7) and a Matlab program, the 3D plots of Tleft(N,ρ) and R left(N, ρ) are plotted in Fig. 2. We can see from the figure that there is an isolated point at N = 149 and ρ = 0.0158 at which both T leftand Rleftare close to infinite. Such result just demonstrates the above analytical formulations. It also demonstrates the existence of the singular PT-symmetry broken point aforementioned. About the point, two features should be mentioned. First, it is achieved through a finite period number and small modulation depth of index. Second, according to the analytical formulations in above section, if d =5λ /(4n0) holds, the wavelength λ can be arbitrary changed and the PT-symmetry broken point keep its position invariant. For example, if we take λ = 15 nm and d = 6.25 nm, the PT-symmetry broken point still occurs at N = 149 and ρ = 0.0158. Such feature provides us much flexibility in designing PT-symmetry dependent device. For simplicity, we do not show the result.
Fig.2 (a) Transmittance and (b) reflectance for left incidence as a function of N and ρ

Full size|PPT slide

The PT-symmetry breaking can be seen from the eigenvalues of scattering matrix S [21,22] that can be deduced from the matrix M in Eq. (4). For PT-symmetric systems, the eigenvalues of the S matrix either form pairs that are unimodular in the exact PT phase, or are with reciprocal moduli in the broken symmetry phase [22]. In Fig. 3, we plot them with different N for ρ = 0.0158. It is seen from the figure that out of the range of 145<N<153, all the two eigenvalues are unimodular. The PT-symmetry broken occurs in the range 145<N<153 in which the two eigenvalues form reciprocal relation. However, in the range 145<N<153, the two eigenvalues take a jumping just at N = 149 at which they are close to infinite and zero, respectively. The singularity of eigenvalues just corresponds to the singular PT-symmetry broken point. The dependence of the singular PT-symmetry broken point on N means that the structure length plays an important role in determining the optical properties of layered PT-symmetric structure. Some similar results have been reported [18,19]. In Ref. [18], authors proposed sinusoidal PT-symmetric complex crystals of finite thickness and found three regimes of lengths in which the structure takes on different properties. In Ref. [19], authors proposed periodic multilayer structures with PT that exhibits anisotropic reflection oscillation patterns as the number of unit-cells is increasing. Different from those in Refs. [18,19], the reflection and transmission in our study are close to infinite and occur at the isolated PT-symmetry broken point.
Fig.3 Two eigenvalues of the S matrix with ρ = 0.0158 as a function of N

Full size|PPT slide

To explore the transmission mechanism of infinite Tleft and Rleft, we perform a simulation through Comsol software based on finite element method. We use a 2D structure in the xy plane instead of infinite one-dimensional structure in x direction and periodic boundary condition in the y direction. We also use two ports at the two interfaces of the 2D structure to excite plane wave with Ez field. The source field value is taken as 0.02 V/m. We perform a frequency domain simulation with incidence wavelength λ = 150 nm and the structure parameters as given in Section 2. The results are shown in Fig. 4(a), through which we find huge field amplification in the PT-symmetric structure. Figure 4(b) is the localized region of Fig. 4(a). The dashed lines denote the interfaces of each unit cell. We clearly find that every interface is exactly at the localized maximum value or minimum value, which means an intense resonance effect occurring at the interfaces. Physically, besides each unit-cell can be function as a resonance unit, the whole structure forms a resonance system. Only if the unit resonance and the system resonance are exactly matched, can the infinite transmittance and reflectance occur. The matched condition is dependent on the values of N and ρ. Because this is an open system, the coupled resonance leads to a huge energy transformation from pumping energy source into the gain regions of the structure. The energy transformation will bring about the infinite transmittance and reflectance.
Fig.4 (a) Ez field distribution inside the PT-symmetric structure at the singular PT-symmetry broken point; (b) localized amplification of (a)

Full size|PPT slide

4 Conclusion

In this paper, we constructed a simple one-dimensional PT-symmetric structure in which the unit cell includes two layers with the conjugated refraction indexes. We found a singular PT-symmetry point corresponding to a special unit cell number and parameter ρ. At the singular PT-symmetry point, infinite transmittance and reflectance simultaneously occur. At the singular point position, both the transmission and reflection are direction-independent. Away from the singular point, the transmittance and reflectance will become finite. Such result means a huge energy transformation from pumping source to light waves. This study demonstrates the quantum theory for the singular point of PT-symmetry system is well in agreement with the classical analysis from the wave optics. Undoubtedly, the strange optical phenomena can be applied in nonlinear optics, energy transformation or other functional optical devices.
1
Joannopoulos J D, Villeneuve P R, Fan S. Photonic crystals: putting new twist on light. Nature, 1997, 386(6621): 143–149

DOI

2
Knight J C, Broeng J, Birks T A, Russell P S J. Photonic band gap guidance in optical fibers. Science, 1998, 282(5393): 1476–1478

DOI PMID

3
Barnes W L, Dereux A, Ebbesen T W. Surface plasmon subwavelength optics. Nature, 2003, 424(6950): 824–830

DOI PMID

4
Shelby R A, Smith D R, Schultz S. Experimental verification of a negative index of refraction. Science, 2001, 292(5514): 77–79

DOI PMID

5
Marani R, D’Orazio A, Petruzzelli V, Rodrigo S G, Martin-Moreno L, Garcia-Vidal F J, Bravo-Abad J. Gain-assisted extraordinary optical transmission through periodic arrays of subwavelength apertures. New Journal of Physics, 2012, 14(1): 013020

6
Feng L, Xu Y L, Fegadolli W S, Lu M H, Oliveira J E, Almeida V R, Chen Y F, Scherer A. Experimental demonstration of a unidirectional reflectionless parity-time metamaterial at optical frequencies. Nature Materials, 2013, 12(2): 108–113

DOI PMID

7
Rüter C E, Makris K G, El-Ganainy R, Christodoulides D N, Segev M, Kip D. Observation of parity–time symmetry in optics. Nature Physics, 2010, 6(3): 192–195

DOI

8
Mostafazadeh A. Spectral singularities of complex scattering potentials and infinite reflection and transmission coefficients at real energies. Physical Review Letters, 2009, 102(22): 220402

DOI PMID

9
Guo A, Salamo G J, Duchesne D, Morandotti R, Volatier-Ravat M, Aimez V, Siviloglou G A, Christodoulides D N. Observation of PT-symmetry breaking in complex optical potentials. Physical Review Letters, 2009, 103(9): 093902

DOI PMID

10
Longhi S. PT-symmetric laser absorber. Physical Review A, 2010, 82(3): 031801

DOI

11
Chong Y D, Ge L, Stone A D. PT-symmetry breaking and laser-absorber modes in optical scattering systems. Physical Review Letters, 2011, 106(9): 093902

DOI PMID

12
Ge L, Chong Y D, Rotter S, Tureci H E, Stone A D. Unconventional modes in lasers with spatially varying gain and loss. Physical Review A, 2011, 84(2): 023820

DOI

13
Nazari F, Nazari M, Moravvej-Farshi M K. A 2×2 spatial optical switch based on PT-symmetry. Optics Letters, 2011, 36(22): 4368–4370

DOI PMID

14
Bender N, Factor S, Bodyfelt J D, Ramezani H, Christodoulides D N, Ellis F M, Kottos T. Observation of asymmetric transport in structures with active nonlinearities. Physical Review Letters, 2013, 110(23): 234101

DOI PMID

15
Nazari F, Bender N, Ramezani H, Moravvej-Farshi M K, Christodoulides D N, Kottos T. Optical isolation via PT-symmetric nonlinear Fano resonances. Optics Express, 2014, 22(8): 9574–9584

DOI PMID

16
Peng B, Özdemir S K, Lei F, Monifi F, Gianfreda M, Long G L, Fan S H, Nori F, Bender C M, Yang L. Parity-time-symmetric whispering-gallery microcavities. Nature Physics, 2014, 10(5): 394–398

DOI

17
Lin Z, Ramezani H, Eichelkraut T, Kottos T, Cao H, Christodoulides D N. Unidirectional invisibility induced by PT-symmetric periodic structures. Physical Review Letters, 2011, 106(21): 213901

DOI PMID

18
Longhi S. Invisibility in PT-symmetric complex crystals. Journal of Physics A, Mathematical and Theoretical, 2011, 44(48): 485302

DOI

19
Zhu X F, Peng Y G, Zhao D G. Anisotropic reflection oscillation in periodic multilayer structures of parity-time symmetry. Optics Express, 2014, 22(15): 18401–18411

DOI PMID

20
Ding S, Wang G P. Extraordinary reflection and transmission with direction dependent wavelength selectivity based on parity-time-symmetric multilayers. Journal of Applied Physics, 2015, 117(2): 023104

DOI

21
Born M, Wolf E. Principles of Optics: Electromagnetic Theory of Propagation, Interference and Diffraction of Light. Elsevier, Cambridge University, 1997

22
Ge L, Chong Y D, Stone A D. Conservation relations and anisotropic transmission resonances in one-dimensional PT-symmetric photonic heterostructures. Physical Review A, 2012, 85(2): 023802

DOI

23
Schomerus H. Quantum noise and self-sustained radiation of PT-symmetric systems. Physical Review Letters, 2010, 104(23): 233601

DOI PMID

Outlines

/