Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Singular properties generated by finite periodic PT-symmetric optical waveguide network

Open Access Open Access

Abstract

In this work we investigate the extraordinary characteristics of one-dimensional (1D) finite periodic parity-time (PT) symmetric network. On the basis of the transfer matrix method, three simple expressions are analytically obtained for transmission, left reflection and right reflection coefficients. For this periodic structure, we provide new criteria for the PT-symmetry breaking transition in terms of the elements of the transfer matrix and the scattering matrix. These criteria indicate that the exceptional points are related only to the cell structure, but not to the cell number. Utilizing these criteria and expressions, the relationships between the transmittances (reflectances) and the cell number are considered in detail. Furthermore, the conditions for ultrastrong transmission are analytically derived. We also show how a PT-symmetric network can become unidirectionally and bidirectionally transparent at specific frequencies. The conditions and related properties of unidirectional and bidirectional transparencies are also examined. Finally, we find that the finite periodic PT-symmetric network with certain cell number can be viewed as a unidirectionally invisible structure at the exceptional points. Our work may pave the way for designing a diversefamily of optical structures and networks with new properties and functionalities.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Designing and fabricating man-made optical materials or structures that can control the propagation of light waves is an unremitting goal all the time. After decades of efforts, people have made many remarkable achievements, such as photonic crystals [1], photonic crystal fibers [2], negative-index metamaterials [3], and so on. However, it is noteworthy that these remarkable works mainly focus on the modulation of the real part of optical parameters to obtain optical properties not found in nature. In recent years, a new family of artificial optical structures, which aims to explore novel optical properties in the complex domain, has attracted considerable attention. This class of optical structures takes advantage of the concept of PT-symmetry in quantum mechanics [4, 5] to attain new exotic optical functionalities and properties. These PT-synthetic structures rely on judiciously balanced gain and loss regions such that the complex refractive index is required to satisfy the condition: n(r)=n*(r). Under this condition, the optical modes can be created and absorbed in a controlled manner. Researchers have always been interested in these PT-symmetric optical systems because it can exhibit some intriguing optical phenomena that are not available in the Hermitian system. These phenomena include double refraction [6], power oscillations [6–8], absorption enhanced transmission [9, 10], unidirectional invisibility [11–14], and coherent perfect absorbers and lasers (CPA-Lasers) [15, 16], etc. Moreover, within the framework of PT-symmetry, other properties also intrigue people’s interests, such as Bloch oscillations [6, 12, 17–20], the Talbot effect [21], andthe realization of optical isolators and circulator [22–25].

In many practical applications, such as metamaterials and plasmonics [3, 26, 27], the loss is an inevitable problem and impedes the progress of applied physics. The concept of PT-symmetry may provide an effective way to overcome losses. However, in practice, it is a very challenging task to fabricate an artificial optical structure with delicately balanced gain and loss elements due to technical difficulties. As a result, many exotic phenomena, like the PT-phase transition or CPA-Lasers, can be experimentally observed only in elemental two-component systems [7, 9, 28, 29]. On the other hand, although there exist many studies on the PT-symmetric periodic systems, most of them mainly focused on layered periodic structures. Therefore, it will be very important to implement PT-symmetry concepts on a variety of new platforms or optical structures. It will help to design new optical devices and materials, such as unidirectional on-chip devices and laser oscillators.

In this paper, we construct a 1D PT-symmetric network which is connected by 1D gain or loss waveguide segments of different length in a periodic manner. This periodic PT-symmetric structure can be described by a simple network equation. Utilizing the network equation and the transfer matrix method, we obtain three expressions for total transmission, left reflection and right reflection coefficients, which depends on the single cell transmission and reflection coefficients, the dimensionless Bloch wavevector, and the total number of unit cells. For this periodic PT-symmetric system, new criteria for the spontaneous PT-symmetry breaking transition are analytically derived. These criteria indicate that the exceptional point is only related to the unit cell structure, and the cell number has no impact on the frequency position and number of the exceptional points. According to these criteria and expressions, the relationships between the transmittances (reflectances) and the cell number are systematically studied. Moreover, the conditions for producing ultrastrong transmission are found for any finite periodic PT-symmetric network. We also demonstrate that our designed PT-symmetric network can display singular scattering characteristics such as unidirectionalor bidirectional transparency. The conditions and related properties of unidirectional and bidirectional transparencies are introduced in detail. Finally, we find that the way of adjusting cell number is an effective route to design a one-way invisible network near the exceptional point. Our work may contribute to the design and manufacture of large-scale PT-symmetric optical lattices or networks, such as two-dimensional or three-dimensional optical networks. This kind of PT-symmetric structure can be arranged more flexible. Therefore, it can be used to implement various PT-synthetic devices or networks with new characteristics and functionalities.

This paper is organized as follows: The model and main equations of the finite periodic PT-symmetric network is introduced in Section 2. In Section 3, we discuss the extraordinary properties of the finite periodic PT-symmetric network and show the results of numerical simulations. Finally, the conclusion is drawn in Section 4.

2. Model and main equations

The designed 1D PT-symmetric periodic optical waveguide network is schematically shown in Fig. 1. This periodic structure contains N cells, and is composed of three kinds of waveguide segments: passive (involving no gain or loss) segments (black lines), loss segments (red lines), and gain segments (green lines). The refractive indices of them are, respectively, n0=2, n1=2+0.005ι and n2=20.005ι, and the refractive index distribution of the whole structure satisfies the condition:n(x)=n*(x). The lengths of the two passive segments connected to nodes 0 and 2N are d. The lengths of the two loss or gain segments between adjacent nodes (blue dots) are d1

(d1=d) and d2, and their length ratio is d1:d2=1:2.

 figure: Fig. 1

Fig. 1 Schematic of the 1D PT-symmetric periodic optical waveguide network containing N cells. Red lines labeled n1 are the loss segments (n1=2+0.005ι), green lines labeled n2 are the gain segments (n2=20.005ι), and black lines labeled n0 are the passive segments (n0=2). The lengths of the two loss (gain) segments between adjacent nodes (blue dots) are d1 and d2, and their length ratio is d1:d2=1:2. The lengths of the two passive segments are d.

Download Full Size | PDF

In this PT-symmetric optical network, the wave function φij(x) in any segment between nodes i and j obeys the Helmholtz equation:

2φij(x)x2+(nijωc)2φij(x)=0,
where nij is the refractive index of the dielectric medium between nodes i and j, ω is the optical wave frequency, and c is the light speed in vacuum. For passive segments 10 and 2N2N+1, Eq. (1) admits the solution:
{φ1,0(x)=α1eιk(xd)+β1eιk(xd),φ2N,2N+1(x)=α2eιkx+β2eιkx,
where the wave vector k=n0ω/c, α1,2 and β1,2 are, respectively, the amplitudes of the forward and backward propagating waves in corresponding passive segment. These amplitudes are related through the transfer matrix MN for N cells:
(α2β2)=MN(α1β1)=(mN(11)mN(12)mN(21)mN(22))(α1β1).

When the optical waves are, separately, incident from the left or right side, the transmission and reflection coefficients for left (L) and right (R) incidence waves are defined as tN(L)=α2/α1, rN(L)=β1/α1, tN(R)=β1/β2 and rN(R)=α2/β2, and we have the boundary conditions: β2=0 for left incidence and α1=0 for right incidence. In terms of the transfer matrix (MN) elements, these coefficients can be expressed as follows [30]:

tN(L)=tN(R)=tN=1mN(22),rN(L)=mN(21)mN(22),rN(R)=mN(12)mN(22).

In order to derive the transfer matrix MN, we will analyze the PT-symmetric waveguide network. The wave function φij(x) in any passive, loss or gain segments can be decomposed into two opposite traveling waves:

φij(x)=αijezijx+βijezijx,
where zij=ιnijω/c, The wave function φij(x) is continuous at the nodes of the PT-symmetric network, and has the boundary conditions: φij(x)|x=0=φi and φij(x)|x=lij=φj. By means of this condition, Eq. (5) can be rewritten as follows:
φij(x)=φisinh [zij(lijx)]sinh (zijlij)+φjsinh (zijx)sinh (zijlij),
where lij is the length of the corresponding waveguide segment. At each node i, the wave function is continuous and the derivative of the wave function at the node i gives the following flux conservation condition [31, 32]:
jφij(x)x|x=0=0.
Substituting Eq. (6) into Eq. (7), we obtain the network equation of the PT-symmetric network:
φiεi+φi1κi1,i+φi+1κi,i+1=0,
where
{εi=qcoth (zi1,idq)+qcoth (zi,i+1dq),κi1,i=qcsch (zi1,idq),κi,i+1=qcsch (zi,i+1dq),
here q is the number of the waveguide segments between adjacent nodes. The network eqation [Eq. (8)] provides the following recursion relation:
(φi+1φi)=(εiκi,i+1κi1,iκi,i+110)(φiφi1)=Γi(φiφi1).

Then, from Eq. (2), we have(φ0φ1)=(11eιkdeιkd)(α1β1)=QL(α1β1), and

(φ2N+1φ2N)=(eιkdeιkd11)(α2β2)=QR(α2β2).

Combining Eqs. (3), (10) and (11), we obtain the transfer matrix MN:

MN=QR1Γ2NΓ2N1Γ2N2Γ3Γ2Γ1Γ0QL,

Thus, the transfer matrix for one cell is M1=QR1Γ2Γ1Γ0QL. Because the PT-symmetric optical network is periodic, these recursion matrices satisfy the following relation:

Γ2N1Γ2N2=Γ2N3Γ2N4==Γ3Γ2=Γ=(Γ11Γ12Γ21Γ22),
where Γ is a unimodular transfer matrix for one period, and its Nth power can be calculated analytically [33]:
ΓN1=(Γ11UN2(u)UN3(u)Γ12UN2(u)Γ21UN2(u)Γ22UN2(u)UN3(u)),
where UN(u) is the Chebyshev Polynomials of the second kind with argument u=1/2(Γ11+Γ22):
UN(u)=sin [(N+1)arccos u]1u2.

By use of the generalized Floquet-Bloch theorem [31], one can derive the dispersion relation: cosK=u, where K is the dimensionless Bloch wavevector. At present, combining Eqs. (12), (14), and (15), we can rewrite the transfer matrix MN as:

MN=(m11sin NKsin Ksin (N1)Ksin Km12sin NKsin K  m21sin NKsin Km22sin NKsin Ksin (N1)Ksin K),
where m11, m12, m21 and m22 are the elements of the transfer matrix M1. In this case, the transmission and reflection coefficients in Eq. (4) become
{tN=sin K(1t1sin NKsin (N1)K)1, rN(L)=r1(L)(1t1sin (N1)Ksin NK)1, rN(R)=r1(R)(1t1sin (N1)Ksin NK)1.

It is found that the transmission coefficient tN for N cells is obtained as a function of t1, K, and N, while the left (right) reflection coefficient rN(L) (rN(R)) for N cells is derived as a function of r1(L) (r1(R)), t1, K, and N. These expressions will help us accurately analyze the extraordinary optical properties and numerically calculate the transmittance T=|tN|2 and reflectances RL=|rN(L)|2 and RR=|rN(R)|2.

3. Discussions and numerical results

3.1. PT-symmetry breaking transition

For the PT-symmetric optical network, the amplitudes of the forward and backward traveling optical waves in passive segments are also related through the scattering matrix S(ω):

(α2β1)=S(ω)(α1β2)=(tNrN(R)rN(L)tN)(α1β2),
and the eigenvalues of S(ω) are s1,2=tN(1±ι1/TN1). The related study of PT-symmetric systems [34] reveal that (i) when TN<1, the eigenvalues are unimodular, and the system is in the PT-symmetric phase; (ii) when TN>1, the eigenvalues form pairs with reciprocal moduli, and the system is in the broken phase; (iii) at the exceptional (EP) point (TN=1), the eigenvalues are degenrate, and the PT-symmetry breaking transition happens. For our designed periodic PT-symmetric network, we can derive more explicit criteria for the transition between the PT-symmetric phase and the broken phase in terms of the elements of the transfer matrix and the scattering matrix. Using Eqs. (16) and (18), we have
s1,2=sin K1/t1×sin NKsin (N1)K(1±ιsin NKsin K1T11).

From this rigorous result, we can conclude that for a finite periodic PT-symmetric network, (i) when T1<1, the eigenvalues are unimodular, and the system is in the PT-symmetric phase; (ii) when T1>1, the eigenvalues are nonunimodular (except the point sinNK=0), and the system is in the broken phase. (iii) when T1=1 or sinNK=0, the eigenvalues are degenerate. But the symmetry-breaking transition can only happen at T1=1; (iv) At the exceptional point, T1 and TN are equal to 1. T1=TN=1 is a necessary condition for the spontaneous PT-symmetry breaking transition. These criteria indicate that for a finite periodic PT-symmetric optical waveguide network, the exceptional point is determined by its unit cell structure, and the cell number N has no effect on the exceptional point.

 figure: Fig. 2

Fig. 2 Logarithm of the absolute values of the eigenvalues of the corresponding scattering matrix for the finite periodic PT-symmetric network with different cell numbers. Red solid lines correspond to N = 1, blue dashed lines correspond to N = 4, and black dashed lines correspond to N = 8.

Download Full Size | PDF

In order to understand these criteria more intuitively and elucidate the influence of the cell number on the exceptional points and the eigenvalues of the scattering matrix, we plot the logarithm of the absolute values of the eigenvalues of the scattering matrix for the finite periodic PT-symmetric network with three different cell numbers (N= 1, 4, 8) in Fig. 2. From Fig. 2, we observe that no matter what the cell number N is, there only exists two exceptional points between frequency 62.12πc/d and 62.20πc/d, and these two exceptional points are, respectively, located at frequency 62.1287πc/d and 62.1895πc/d. Between these two exceptional points, the eigenvalues are generally nonunimodular, and the system is in the broken phase. However, when N>1, the eigenvalues become degenerate and are unimodular at the point sinNK=0. On both sides of these two exceptional points, the eigenvalues are unimodular, and the system is in the PT-symmetric phase. Similarly, in the PT-symmetric phase, the eigenvalues can be also degenerate when sinNK=0. Although the eigenvalues become degenerate at the point sinNK=0, the symmetry-breaking transition does not happen. Thus, we can conclude that the cell number has no effect on the number and frequency location of the exceptional points. From Fig. 2(a), one can also observe that the cell number N has great influence on the eigenvalues of the scattering matrix. When N = 4, the absolute value of one eigenvalue of the scattering matrix is very large, while the other approaches zero at frequency 62.177πc/d. However, when N is equal to 1 or 8, the absolute value of one eigenvalue is relatively small and the other is relatively big at this frequency.

Next, we will show the influence of the cell structure on the PT-symmetry breaking transition. In Fig. 3(b), we plot the logarithm of the absolute values of the eigenvalues of the scattering matrix for the PT-symmetric network with different waveguide length ratios. Here, the PT-symmetric network contains only one cell. From Fig. 3(a), we observe that when the length ratio is d1:d2=1:1, the eigenvalues are unimodular, and the system is in the PT-symmetric phase. When we keep the length d1 unchanged and extend the length d2, the frequency region where the system is in the broken phase becomes larger and larger. From Fig. 3(b), we observe that when the length ratio is d1:d2=1:5, the system is completely in the broken phase. when the length ratio is d1:d2=1:8, the system is still in the broken phase, and the eigenvalues of these two cases are degenerate at the same frequency. If we further extend the length d2, the absolute values of the eigenvalues are approximately equal to the case: d1:d2=1:8. (In the following investigation, we still focus on the PT-symmetric network with d1:d2=1:2)

 figure: Fig. 3

Fig. 3 Logarithm of the absolute values of the eigenvalues of the corresponding scattering matrix for the PT-symmetric network with different length ratios. Here the PT-symmetric network contains only one cell. (a) Green solid lines correspond to d1:d2=1:1, magenta solid lines correspond to d1:d2=1:3, and blue solid lines correspond to d1:d2=1:4. (b) Cyan solid lines correspond to d1:d2=1:5, black solid lines correspond to d1:d2=1:8, and red dashed lines correspond to d1:d2=1:17.

Download Full Size | PDF

3.2. The influence of cell number on transmission and reflection

After introducing the PT-symmetry breaking transition of the finite periodic PT-symmetric optical network, it is necessary to investigate the transmission and reflection properties related to the number of unit cells. From Eq. (17), the total transmission and reflection for a finite periodic PT-symmetric network are obtained as a function of the cell number N.

Thus, we may argue that the transmission and reflection curves should oscillate as the cell number N increases. However, in fact, the results may be different from our intuition. According to the discussions in Subsection 3.1, we have known that TN and T1 are equal to 1 at the exceptional point. Combining this condition and Eq. (17), we can deduce that when T1=1, TN will be equal to 1 whether the point T1=1 is the exceptional point or not. Therefore, the transmission curve does not oscillate with the cell number N when T1=1. When T11, the transmission curve will oscillate as the cell number N is increasing. Moreover, according to Fig. 2, we can also deduce that the transmission TN for N cells satisfies the following relations: (i) when T1<1, TN1 ; (ii) when T1>1, TN1. Both cases are equal to 1 at the point sinNK=0. In Fig. 4, we have, respectively, calculated the transmittance, left reflectance and right reflectance of the finite periodic PT-symmetric network with the cell number N increasing from 1 to 30 around the point T1=1. These results in Fig. 4(a) are in excellent agreement with our prediction. From Fig. 4(b) and 4(c), we observe that if the point T1=1 is the exceptional point, the left reflection is always zero regardless of the increment of the cell number N whereas the right reflection curve oscillates with the cell number N; if the point T1=1 is not the exceptional point, both left and right reflections are always zero no matter what the cell number N is. In the PT-symmetric phase (T1<1) or the broken phase (T1>1), both left and right reflection curves will oscillate with the cell number N, and they will vanish at sinNK=0.

 figure: Fig. 4

Fig. 4 (a) Transmittance (TN), (b) left reflectance (RN(L)), and (c) right reflectance (RN(R)) for the finite periodic PT-symmetric network as a function of cell number N. The black dot-solid lines correspond to the PT-symmetric phase (T1<1), the red circle-solid lines correspond to the exceptional point (T1=1), the blue dot-solid lines correspond to the broken phase (T1>1), and the cyan dot-dashed lines correspond to the degeneracy point (T1=1, not the exceptional point).

Download Full Size | PDF

3.3. Conditions for generating ultrastrong transmission

In Fig. 2, when N = 4, the absolute value of one eigenvalue of the scattering matrix is very large at frequency 62.177πc/d. It implies that the transmission for N = 4 is very strong at this frequency. However, when N=1 or 8, the transmission is relatively weak at frequency 62.177πc/d. Here, we will analyse this phenomenon and derive the conditions for generating ultrastrong transmission. For 1D PT-symmetric system, the elements of the transfer matrix MN have the property: mN(22)(ω)=mN(11)*(ω*) [14]. Using this property and Eq. (4), we can rewrite the transfer matrix M1, which has the form: M1=(1/t1*,rR/t1;rL/t1,1/t1). Meanwhile, the eigenvalues of M1 are written as λ1=1/λ2=exp(ιK). According to the theory of matrices, we can deduce 1/t1=cosK+ιIm(1/t1). Substituting this equation into Eq. (17), the transmission coefficient tN can be rewritten as

tN=1cos NK+ιIm(1/t1)sin NKsin K.

Since cosNK and sinNK cannot be equal to zero at the same time, from this equation, the transmission coefficient tN goes to infinity only when cosNK=0 (sinNK=±1) and Im(1/t1)=0. In Fig. 5(a) and 5(b), we show the imaginary part of 1/t1 and the values of sinNK with N=34 in the frequency range 62.10πc/d62.40πc/d. The transmittances of the periodic PT-symmetric network for N = 1 and N=34 are displayed in Fig. 5(c). From Fig. 5(a) and 5(b), we observe that at the frequency 62.351πc/d, the imaginary part of 1/t1 is zero and the value of sin34K is approaching zero. The ultrastrong transmission (red line) can be clearly observed near the frequency 62.351πc/d in Fig. 5(c). However, the transmission is relatively weak at frequency 62.177πc/d due to sin34K0.4. For a finite periodic PT-symmetric network, the condition Im(1/t1)=0 designates the frequency region where the ultrastrong transmission can be created. Then, we can obtain ultrastrong transmission by adjusting the cell number to satisfy the condition: sinNK=±1 (cosNK=0). Under these two conditions (Im(1/t1)=0, sinNK=±1/cosNK=0), the ultrastrong transmission can be created in a controlled manner. It is noteworthy that these two stong transmissions in Fig. 5(c) are realized in the broken phase (T1>1).

 figure: Fig. 5

Fig. 5 (a) The imaginary part of 1/t1 between frequency 62.10πc/d and 62.40πc/d, (b) the values of sinNK with N=34 in the same frequency range, and (c) the transmittance of the periodic PT-symmetric network with two different cell numbers, where the green line corresponds to N = 1, and the red line corresponds to N=34

Download Full Size | PDF

3.4. Unidirectional and bidirectional transparencies

In PT-symmetric system, we have S*(ω)=PTS(ω)PT=S1(ω). From this relation, one can deduce that rN(L)rN(R)*=1|tN|2, and this equation implies that RN(L)RN(R)=|TN1| [35]. When we have perfect transmission TN = 1, the geometric mean of the two reflectances RN(L)RN(R) is required to vanish. This constraint can be satisfied by the following two cases: the first case is that only one reflection becomes zero, and the second case is that both reflections are all zero. From Eqs. (4) and (17), the left reflectance RN(L) is generally different from the right reflectance RN(R) due to r1(L)r1(R), and the left and right reflectances can be equivalent at some frequencies (i.e. sinNK=0). Therefore, these two cases can be realized in the PT-symmetric system. When only one reflection vanishes, we have reflectionless perfect transmission in one direction but not the other, and this phenomenon is the unidirectional transparency. When both reflections vanish, we have reflectionless perfect transmission in each direction, and this phenomenon is the bidirectional transparency.

We first examine the unidirectional transparency properties of the PT-symmetric network. This phenomenon is closely related to the PT-symmetric structures, and it can only happen at the exceptional point. Based on Eq. (19), we have calculated the logarithm of the absolute values of the eigenvalues of the scattering matrix for the finite periodic PT-symmetric network with four cells in Fig. 6(a). The transmittance, left reflectance and right reflectance of optical waves propagating through the PT-symmetric network are shown in Fig. 6(b). The results in Fig. 6(a) and 6(b) show that the exceptional point is at 62.1895πc/d, and the transmittance is unitary at the frequency ω=62.1895πc/d. Below the exceptional point, the eigenvalues are nonunimodular and the system is in the PT-broken phase. The transmittance is superunitary. While above the exceptional point, the eigenvalues are unimodular and the system is in the PT-symmetric phase. The transmittance is subunitary. At the exceptional point, the reflectance for left incident waves vanishes whereas the reflectance for right incident waves is RN(R)15. This phenomenon is the so-called unidirectional transparency mentioned above. From Eq. (17), the left and right reflection coefficients for N cells have the same term: 1[t1sin(N1)K]/sinNK. Therefore, the conditions for unidirectional transparency can be expressed as |t1|=1, r1(L)=0 and r1(R)0, or |t1|=1, r1(R)=0 and r1(L)0.

 figure: Fig. 6

Fig. 6 Unidirectional transparency characteristics of the finite periodic PT-symmetric network with four cells. (a) Logarithm of the absolute values of the eigenvalues of the corresponding scattering matrix. (b) Transmittance, left reflectance and right reflectance. Phases and delay times of transmitted, left-reflected, and right-reflected waves of the finite periodic PT-symmetric are shown in (c) and (d).

Download Full Size | PDF

Unidirectional transparency is a singularity with highly directional responses. The physical mechanism of this singularity is the resonance trapping [36]. The occurrence of resonance trapping for transmission and reflection can be judged from the divergence of the corresponding delay time. The delay time is associated with the time experienced by a wave packet inside the potential before it exits from the scattering region, and it is defined as τt,r=dϕt,r/dω, where ϕt,r is the phase of the transmission or reflection coefficient. In Fig. 6(c), we plot the phases of transmitted, left-reflected, and right-reflected waves. At the exceptional point (ω=62.1895πc/d), there is an abrupt change of π in the phase of left-reflected waves. It indicates that the delay time for the left reflection diverges. In fact, in Fig. 6(d), the delay time for left reflection also behaves like a δ(ω) function due to the step-function behavior of the phase. As a result, the light reflecting from the left side will be trapped in the lossy waveguide segment for a long time and get absorbed completely by the loss. While the light reflecting from the right side will delay in the gain waveguide segment for a relatively short time and is enhanced. It is also necessary to note that the phase relation between transmitted and reflected waves. From Fig. 6(c), in the broken phase, there is a π phase difference between the left-reflected and right-reflected waves, while in the PT-symmetric phase, the phase of the left-reflected waves is the same as that of the right-reflected waves. Moreover, there is a ±π/2 phase difference between transmitted waves and reflected waves either in the broken phase or in the PT-symmetric phase. These results can be also derived from the above equation, that is, rN(L)rN(R)*=1|tN|2. Since transmitted and reflected waves have such phase relations, the delay times for transmission and reflections are exactly equivalent in the broken phase and the PT-symmetric phase [Fig. 6(d)].

 figure: Fig. 7

Fig. 7 Bidirectional transparency characteristics of the finite periodic PT-symmetric network with four cells in the PT-symmetric phase. (a) Logarithm of the absolute values of the eigenvalues of the corresponding scattering matrix. (b) Transmittance, left reflectance and right reflectance. Phases and delay times of transmitted, left-reflected, and right-reflected waves of the PT-symmetric are shown in (c) and (d).

Download Full Size | PDF

 figure: Fig. 8

Fig. 8 Bidirectional transparency characteristics of the finite periodic PT-symmetric network with four cells in the broken phase. (a) Logarithm of the absolute values of the eigenvalues of the corresponding scattering matrix. (b) Transmittance, left reflectance and right reflectance. Phases and delay times of the transmitted, left-reflected, and right-reflected waves of the PT-symmetric are shown in (c) and (d).

Download Full Size | PDF

Next, we analyze the bidirectional transparency properties. From Eq. (17), the conditions for this phenomenon are written as r1(L)=r1(R)=0 or sinNK=0. The first case is related to the unit cell structure of the finite periodic PT-symmetric network, and the second case is related to the periodicity of the PT-symmetric network. The number of points corresponding to bidirectional reflectionless increases with the number of cells. According to Subsection 3.1, the eigenvalues of the scattering matrix can be degenerate in the PT-symmetric phase and the broken phase, thus bidirectional transparency can happen in the PT-symmetric phase and the broken phase. In this paper, we show the bidirectional transparency properties of the finite periodicPT-symmetric network with four cells in the PT-symmetric phase (Fig. 7) and the broken phase (Fig. 8). From Fig. 7, one can observe that the eigenvalues of the corresponding scattering matrix are unimodular and degenerate at ω=62.902πc/d, as marked by a darkened dot. At the degeneracy point ω=62.902πc/d, the transmittance TN is unitary. Away from this point, the transmittance TN are subunitary. Although the left and right reflectances are quite different inthe frequency domain [62.88πc/d,62.902πc/d) and (62.902πc/d,62.92πc/d], they both vanish at the degeneracy point 62.902πc/d. Since the system is in the PT-symmetric phase, the phases of left-reflected and right-reflected waves must be equivalent. the phase difference between transmitted waves and reflected waves is still ±π/2. At the degeneracy point, the phases of two reflected waves experience an abrupt change of π simultaneously. Therefore, the delay time for both reflections all behaves as a delta function. Both left-reflected and right-reflected waves are trapped in the PT-symmetric network for a long time and get absorbed. Another case for the bidirectional reflectionless is in the broken phase. From Fig. 8, one can observe that the absolute values of a scattering matrix is unimodular and degenerate at the frequency 62.1474πc/d. On both sides of this frequency 62.1474πc/d, the absolute values of the eigenvalues are nonunimodular, and the transmission is superunitary. The phases of the two reflected waves all have an abrupt change of π, but the directions of two abrupt changes are opposite. Thus the delay time of two reflected waves satisfies the relation: δL(ω)=δR(ω) [Fig. 8(d)]. It is necessary to note that all bidirectional transparencies are caused by Bragg resonances, which are mainly attributed to the increase of the cell number.

3.5. Unidirectional invisibility

Over the past decades, optical cloaking has attracted considerable attention, specifically in connection to transformation optics [37, 38]. In recent years, a new notion of unidirectional invisibility, which results from a completely different physical process, was proposed [11]. This kind of invisibility is based on the concept of PT-symmetry and arises because of the spontaneous PT-symmetry breaking. It is well known that at the exceptional point, the PT-symmetric network can display unidirectional transparency when the light is incident from different sides of the structure. However, unidirectional transparency is not in general equivalent to unidirectional invisibility. This is due to the fact that the phase of transmitted waves might be different from the phase of the waves after passing through a homogeneous background medium with the same length. For example, in Fig. 6(c), the transmission is exactly unitary at the exceptional point and the phase of the transmitted wave is 0.26π. However, for a homogeneous background medium, the phase of transmitted waves is 0.5π. In this case, simple time-of-flight measurements can sense the existence of the PT-symmetric optical network, and the PT-symmetric optical network cannot be viewed as an invisible structure. Therefore, if we want to design an invisible structure, it is necessary to examine the phase of the transmitted wave and compare it with the phase of a wave propagating in a homogeneous waveguide network. Using Eq. (17), we deduce that the phase ϕtN for N cells is

ϕtN=1ιln (|1/t1sin NKsin (N1)K|1/t1sin NKsin (N1)K).

One can observe that ϕtN is closely related to the number of cells, and the cell number has no impact on the frequency positon and the number of the exceptional points. It signifies that we can design a unidirectionally invisible PT-symmetric network by adjusting the cell number of the structure at the exceptional points. In order to compare the phases of the transmitted waves traveling through the PT-symmetric network (ϕPT) with those of the transmitted waves traveling through a homogeneous background medium (ϕ), we plot the absolute value of the phase difference of these two phases (|ϕPTϕ|) as a function of the cell number at the exceptional point (ω=62.1895πc/d) in Fig. 9. From this figure, one can observe that the phase difference has an approximate periodicity and approaches zero when the cell number N is 1, 25, 53 and 79. Therefore, the finite periodic PT-symmetric network with such cell number can be seen as a unidirectionally invisible structure. In summary, adjusting cell number is a valid way to realize optical cloaking in PT-symmetric networks at given frequencies.

 figure: Fig. 9

Fig. 9 The absolute value of the phase differences for transmitted waves traveling through the PT-symmetric network (ϕPT) and those traveling through a homogeneous background medium (ϕ) vary with cell number N at ω=62.1895πc/d.

Download Full Size | PDF

4. Conclusions

We have investigated the singular properties of the finite periodic PT-symmetric network. Based on theoretical derivation, we obtain three expressions for transmission and left and right reflection coefficients. New criteria for the PT-symmetry breaking transition of the periodic PT-symmetric network are analytically derived. These criteria indicate that the exceptional points of the finite periodic PT-symmetric network are only related to the single cell structure, and the cell number has no effect on the frequency location and the number of the exceptional points. In terms of these criteria, the relationship between the transmittance (reflectances) and the cell number is systematically examined. More importantly, we find the conditions for ultrastrong transmission, which are applicable to an arbitrary 1D periodic PT-symmetric network. Besides, we also study the phenomenon of unidirectional and bidirectional transparencies. Using the transmission and reflection coefficient expressions, the conditions forunidirectional and bidirectional transparencies are analytically obtained, and the related optical properties are well introduced. Finally, we find that the finite periodic PT-symmetric network with 1, 27, 53, or 79 unit cells is unidirectionally invisible at the exceptional point. Our work will deepen people’s understanding of the optical properties of PT-symmetric waveguide network. We envision that PT-symmetric networks will open a new route for designing various optical structures or networks with new characteristics and functionalities.

Funding

National Natural Science Foundation of China, Grant Nos. 11674107, 61475049, 11775083, 61774062, and 61771205; Natural Science Foundation of Guangdong Province, Grant No. 2015A030313374.

References

1. J. D. Joannopoulos, P. R. Villeneuve, and S. Fan, “Photonic crystals: putting a new twist on light,” Nature 386, 143 (1997). [CrossRef]  

2. J. C. Knight, J. Broeng, T. A. Birks, and P. S. J. Russell, “Photonic Band Gap Guidance in Optical Fibers,” Science 282, 1476 (1998). [CrossRef]   [PubMed]  

3. J. B. Pendry, “Negative Refraction Makes a Perfect Lens,” Phys. Rev. Lett. 85, 3966–3969 (2000). [CrossRef]   [PubMed]  

4. C. M. Bender and S. Boettcher, “Real Spectra in Non-Hermitian Hamiltonians Having PT Symmetry,” Phys. Rev. Lett. 80, 5243–5246 (1998). [CrossRef]  

5. C. M. Bender, “Making sense of non-Hermitian Hamiltonians,” Rep. Prog. Phys. 70, 947 (2007). [CrossRef]  

6. K. G. Makris, R. El-Ganainy, D. N. Christodoulides, and Z. H. Musslimani, “Beam Dynamics in PT Symmetric Optical Lattices,” Phys. Rev. Lett. 100, 103904 (2008). [CrossRef]  

7. C. E. Rüter, K. G. Makris, R. El-Ganainy, D. N. Christodoulides, M. Segev, and D. Kip, “Observation of parity–time symmetry in optics,” Nat. Phys. 6, 192 (2010). [CrossRef]  

8. M. C. Zheng, D. N. Christodoulides, R. Fleischmann, and T. Kottos, “PT optical lattices and universality in beam dynamics,” Phys. Rev. A 82, 010103 (2010). [CrossRef]  

9. A. Guo, G. J. Salamo, D. Duchesne, R. Morandotti, M. Volatier-Ravat, V. Aimez, G. A. Siviloglou, and D. N. Christodoulides, “Observation of PT-Symmetry Breaking in Complex Optical Potentials,” Phys. Rev. Lett. 103, 093902 (2009). [CrossRef]  

10. Z. Zhang, P. Miao, J. Sun, S. Longhi, N. M. Litchinitser, and L. Feng, “Elimination of Spatial Hole Burning in Microlasers for Stability and Efficiency Enhancement,” ACS Photonics 5, 3016–3022 (2018). [CrossRef]  

11. Z. Lin, H. Ramezani, T. Eichelkraut, T. Kottos, H. Cao, and D. N. Christodoulides, “Unidirectional Invisibility Induced by PT-Symmetric Periodic Structures,” Phys. Rev. Lett. 106, 213901 (2011). [CrossRef]  

12. A. Regensburger, C. Bersch, M. A. Miri, G. Onishchukov, D. N. Christodoulides, and U. Peschel, “Parity-time synthetic photonic lattices,” Nature 488, 167–171 (2012). [CrossRef]   [PubMed]  

13. X. Zhu, L. Feng, P. Zhang, X. Yin, and X. Zhang, “One-way invisible cloak using parity-time symmetric transformation optics,” Opt. Lett. 38, 2821–2824 (2013). [CrossRef]   [PubMed]  

14. L. Feng, X. Zhu, S. Yang, H. Zhu, P. Zhang, X. Yin, Y. Wang, and X. Zhang, “Demonstration of a large-scale optical exceptional point structure,” Opt. Express 22, 1760–1767 (2014). [CrossRef]   [PubMed]  

15. S. Longhi, “PT-symmetric laser absorber,” Phys. Rev. A 82, 031801 (2010). [CrossRef]  

16. Y. D. Chong, L. Ge, and A. D. Stone, “PT-Symmetry Breaking and Laser-Absorber Modes in Optical Scattering Systems,” Phys. Rev. Lett. 106, 093902 (2011). [CrossRef]  

17. S. Longhi, “Bloch Oscillations in Complex Crystals with PT Symmetry,” Phys. Rev. Lett. 103, 123601 (2009). [CrossRef]  

18. S. Longhi, “Spectral singularities and Bragg scattering in complex crystals,” Phys. Rev. A 81, 022102 (2010). [CrossRef]  

19. S. Ding and G. P. Wang, “Nonreciprocal optical Bloch-Zener oscillations in ternary parity-time-symmetric waveguide lattices,” Appl. Phys. Lett. 100, 151913 (2012). [CrossRef]  

20. M. Wimmer, M.-A. Miri, D. Christodoulides, and U. Peschel, “Observation of Bloch oscillations in complex PT-symmetric photonic lattices,” Sci. Rep. 5, 17760 (2015). [CrossRef]   [PubMed]  

21. H. Ramezani, D. N. Christodoulides, V. Kovanis, I. Vitebskiy, and T. Kottos, “PT-Symmetric Talbot Effects,” Phys. Rev. Lett. 109, 033902 (2012). [CrossRef]  

22. H. Ramezani, T. Kottos, R. El-Ganainy, and D. N. Christodoulides, “Unidirectional nonlinear PT-symmetric optical structures,” Phys. Rev. A 82, 043803 (2010). [CrossRef]  

23. N. Bender, S. Factor, J. D. Bodyfelt, H. Ramezani, D. N. Christodoulides, F. M. Ellis, and T. Kottos, “Observation of Asymmetric Transport in Structures with Active Nonlinearities,” Phys. Rev. Lett. 110, 234101 (2013). [CrossRef]   [PubMed]  

24. F. Nazari, N. Bender, H. Ramezani, M. K. Moravvej-Farshi, D. N. Christodoulides, and T. Kottos, “Optical isolation via PT-symmetric nonlinear Fano resonances,” Opt. Express 22, 9574–9584 (2014). [CrossRef]   [PubMed]  

25. C. Yidong, “Nonlinear Optics: Asymmetry from Symmetry,” Nat. Phys. 10, 336 (2014). [CrossRef]  

26. P. Berini and I. De Leon, “Surface plasmon–polariton amplifiers and lasers,” Nat. Photonics 6, 16 (2011). [CrossRef]  

27. M. A. Noginov, G. Zhu, A. M. Belgrave, R. Bakker, V. M. Shalaev, E. E. Narimanov, S. Stout, E. Herz, T. Suteewong, and U. Wiesner, “Demonstration of a spaser-based nanolaser,” Nature 460, 1110 (2009). [CrossRef]   [PubMed]  

28. J. Schindler, A. Li, M. C. Zheng, F. M. Ellis, and T. Kottos, “Experimental study of active LRC circuits with PT symmetries,” Phys. Rev. A 84, 040101 (2011). [CrossRef]  

29. S. Bittner, B. Dietz, U. Günther, H. L. Harney, M. Miski-Oglu, A. Richter, and F. Schäfer, “PT Symmetry and Spontaneous Symmetry Breaking in a Microwave Billiard,” Phys. Rev. Lett. 108, 024101 (2012). [CrossRef]  

30. A. Mostafazadeh, “Spectral Singularities of Complex Scattering Potentials and Infinite Reflection and Transmission Coefficients at Real Energies,” Phys. Rev. Lett. 102, 220402 (2009). [CrossRef]   [PubMed]  

31. S. Alexander, “Superconductivity of networks. A percolation approach to the effects of disorder,” Phys. Rev. B 27, 1541–1557 (1983). [CrossRef]  

32. Z.-Q. Zhang, C. C. Wong, K. K. Fung, Y. L. Ho, W. L. Chan, S. C. Kan, T. L. Chan, and N. Cheung, “Observation of Localized Electromagnetic Waves in Three-Dimensional Networks Of Waveguides,” Phys. Rev. B 81, 5540–5543 (1998).

33. P. Yeh, A. Yariv, and C. S. Hong, “Electromagnetic propagation in periodic stratified media. I. General theory*,” J. Opt. Soc. Am. 67, 423–438 (1977). [CrossRef]  

34. X. Zhu, H. Ramezani, C. Shi, J. Zhu, and X. Zhang, “PT-Symmetric Acoustics,” Phys. Rev. X 4, 031042 (2014).

35. L. Ge, Y. D. Chong, and A. D. Stone, “Conservation relations and anisotropic transmission resonances in one-dimensional PT-symmetric photonic heterostructures,” Phys. Rev. A 85, 023802 (2012). [CrossRef]  

36. H. Ramezani, H.-K. Li, Y. Wang, and X. Zhang, “Unidirectional Spectral Singularities,” Phys. Rev. Lett. 113, 263905 (2014). [CrossRef]  

37. J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling Electromagnetic Fields,” Science 312, 1780 (2006). [CrossRef]   [PubMed]  

38. U. Leonhardt, “Optical Conformal Mapping,” Science 312, 1777 (2006). [CrossRef]   [PubMed]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (9)

Fig. 1
Fig. 1 Schematic of the 1D PT-symmetric periodic optical waveguide network containing N cells. Red lines labeled n1 are the loss segments ( n 1 = 2 + 0.005 ι ), green lines labeled n2 are the gain segments ( n 2 = 2 0.005 ι ), and black lines labeled n0 are the passive segments ( n 0 = 2 ). The lengths of the two loss (gain) segments between adjacent nodes (blue dots) are d1 and d2, and their length ratio is d 1 : d 2 = 1 : 2. The lengths of the two passive segments are d.
Fig. 2
Fig. 2 Logarithm of the absolute values of the eigenvalues of the corresponding scattering matrix for the finite periodic PT-symmetric network with different cell numbers. Red solid lines correspond to N = 1, blue dashed lines correspond to N = 4, and black dashed lines correspond to N = 8.
Fig. 3
Fig. 3 Logarithm of the absolute values of the eigenvalues of the corresponding scattering matrix for the PT-symmetric network with different length ratios. Here the PT-symmetric network contains only one cell. (a) Green solid lines correspond to d 1 : d 2 = 1 : 1, magenta solid lines correspond to d 1 : d 2 = 1 : 3, and blue solid lines correspond to d 1 : d 2 = 1 : 4. (b) Cyan solid lines correspond to d 1 : d 2 = 1 : 5, black solid lines correspond to d 1 : d 2 = 1 : 8, and red dashed lines correspond to d 1 : d 2 = 1 : 17.
Fig. 4
Fig. 4 (a) Transmittance (TN), (b) left reflectance ( R N ( L )), and (c) right reflectance ( R N ( R )) for the finite periodic PT-symmetric network as a function of cell number N. The black dot-solid lines correspond to the PT-symmetric phase ( T 1 < 1 ), the red circle-solid lines correspond to the exceptional point ( T 1 = 1 ), the blue dot-solid lines correspond to the broken phase ( T 1 > 1 ), and the cyan dot-dashed lines correspond to the degeneracy point ( T 1 = 1, not the exceptional point).
Fig. 5
Fig. 5 (a) The imaginary part of 1 / t 1 between frequency 62.10 π c / d and 62.40 π c / d, (b) the values of sin N K with N = 34 in the same frequency range, and (c) the transmittance of the periodic PT-symmetric network with two different cell numbers, where the green line corresponds to N = 1, and the red line corresponds to N = 34
Fig. 6
Fig. 6 Unidirectional transparency characteristics of the finite periodic PT-symmetric network with four cells. (a) Logarithm of the absolute values of the eigenvalues of the corresponding scattering matrix. (b) Transmittance, left reflectance and right reflectance. Phases and delay times of transmitted, left-reflected, and right-reflected waves of the finite periodic PT-symmetric are shown in (c) and (d).
Fig. 7
Fig. 7 Bidirectional transparency characteristics of the finite periodic PT-symmetric network with four cells in the PT-symmetric phase. (a) Logarithm of the absolute values of the eigenvalues of the corresponding scattering matrix. (b) Transmittance, left reflectance and right reflectance. Phases and delay times of transmitted, left-reflected, and right-reflected waves of the PT-symmetric are shown in (c) and (d).
Fig. 8
Fig. 8 Bidirectional transparency characteristics of the finite periodic PT-symmetric network with four cells in the broken phase. (a) Logarithm of the absolute values of the eigenvalues of the corresponding scattering matrix. (b) Transmittance, left reflectance and right reflectance. Phases and delay times of the transmitted, left-reflected, and right-reflected waves of the PT-symmetric are shown in (c) and (d).
Fig. 9
Fig. 9 The absolute value of the phase differences for transmitted waves traveling through the PT-symmetric network ( ϕ PT ) and those traveling through a homogeneous background medium (ϕ) vary with cell number N at ω = 62.1895 π c / d .

Equations (21)

Equations on this page are rendered with MathJax. Learn more.

2 φ i j ( x ) x 2 + ( n i j ω c ) 2 φ i j ( x ) = 0 ,
{ φ 1 , 0 ( x ) = α 1 e ι k ( x d ) + β 1 e ι k ( x d ) , φ 2 N , 2 N + 1 ( x ) = α 2 e ι k x + β 2 e ι k x ,
( α 2 β 2 ) = M N ( α 1 β 1 ) = ( m N ( 11 ) m N ( 12 ) m N ( 21 ) m N ( 22 ) ) ( α 1 β 1 ) .
t N ( L ) = t N ( R ) = t N = 1 m N ( 22 ) , r N ( L ) = m N ( 21 ) m N ( 22 ) , r N ( R ) = m N ( 12 ) m N ( 22 ) .
φ i j ( x ) = α i j e z i j x + β i j e z i j x ,
φ i j ( x ) = φ i sinh  [ z i j ( l i j x ) ] sinh  ( z i j l i j ) + φ j sinh  ( z i j x ) sinh  ( z i j l i j ) ,
j φ i j ( x ) x | x = 0 = 0 .
φ i ε i + φ i 1 κ i 1 , i + φ i + 1 κ i , i + 1 = 0 ,
{ ε i = q coth  ( z i 1 , i d q ) + q coth  ( z i , i + 1 d q ) , κ i 1 , i = q csch  ( z i 1 , i d q ) , κ i , i + 1 = q csch  ( z i , i + 1 d q ) ,
( φ i + 1 φ i ) = ( ε i κ i , i + 1 κ i 1 , i κ i , i + 1 1 0 ) ( φ i φ i 1 ) = Γ i ( φ i φ i 1 ) .
( φ 2 N + 1 φ 2 N ) = ( e ι k d e ι k d 1 1 ) ( α 2 β 2 ) = Q R ( α 2 β 2 ) .
M N = Q R 1 Γ 2 N Γ 2 N 1 Γ 2 N 2 Γ 3 Γ 2 Γ 1 Γ 0 Q L ,
Γ 2 N 1 Γ 2 N 2 = Γ 2 N 3 Γ 2 N 4 = = Γ 3 Γ 2 = Γ = ( Γ 11 Γ 12 Γ 21 Γ 22 ) ,
Γ N 1 = ( Γ 11 U N 2 ( u ) U N 3 ( u ) Γ 12 U N 2 ( u ) Γ 21 U N 2 ( u ) Γ 22 U N 2 ( u ) U N 3 ( u ) ) ,
U N ( u ) = sin  [ ( N + 1 ) arccos  u ] 1 u 2 .
M N = ( m 11 sin  N K sin  K sin  ( N 1 ) K sin  K m 12 sin  N K sin  K     m 21 sin  N K sin  K m 22 sin  N K sin  K sin  ( N 1 ) K sin  K ) ,
{ t N = sin  K ( 1 t 1 sin  N K sin  ( N 1 ) K ) 1 ,   r N ( L ) = r 1 ( L ) ( 1 t 1 sin  ( N 1 ) K sin  N K ) 1 ,   r N ( R ) = r 1 ( R ) ( 1 t 1 sin  ( N 1 ) K sin  N K ) 1 .
( α 2 β 1 ) = S ( ω ) ( α 1 β 2 ) = ( t N r N ( R ) r N ( L ) t N ) ( α 1 β 2 ) ,
s 1 , 2 = sin  K 1 / t 1 × sin  N K sin  ( N 1 ) K ( 1 ± ι sin  N K sin  K 1 T 1 1 ) .
t N = 1 cos  N K + ι Im ( 1 / t 1 ) sin  N K sin  K .
ϕ t N = 1 ι ln  ( | 1 / t 1 sin  N K sin  ( N 1 ) K | 1 / t 1 sin  N K sin  ( N 1 ) K ) .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.