Next Article in Journal
An Adaptive Search Algorithm for Multiplicity Dynamic Flexible Job Shop Scheduling with New Order Arrivals
Previous Article in Journal
An Online Review-Driven Picture Fuzzy Multi-Criteria Group Decision-Making Approach for Evaluating the Online Medical Service Quality of Doctors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Transport through Parity–Time Symmetric Lattice Potentials

Department of Physics, Shanghai University, Shanghai 200444, China
*
Author to whom correspondence should be addressed.
Symmetry 2024, 16(6), 640; https://doi.org/10.3390/sym16060640
Submission received: 22 April 2024 / Revised: 13 May 2024 / Accepted: 16 May 2024 / Published: 22 May 2024
(This article belongs to the Special Issue Symmetry-Related Quantum Phases in Exciton-Polariton Condensates)

Abstract

:
We study nonlinear transports of a light field through finite parity–time symmetric lattice potentials. The initial light field is trapped in a source reservoir and is released to expand toward the lattice potentials along the transverse direction due to the nonlinearity. We identify the transports that can be classified into in-band and in-gap transports. In the in-band transport, the light field can tunnel through the lattices into the sink reservoir, and in the in-gap transport, the light field is self-trapped inside the lattices to form a solitary wave.

1. Introduction

Transport properties through potential barriers are fundamentally interesting in quantum physics. They are standard textbook problems. In solid-state systems, a two-dimensional electron gas confined in structured semiconductor heterostructures provides an experimentally controllable means to investigate transport phenomena [1,2,3,4]. Light fields and ultracold atomic gases have been fundamentally physical realizations of quantum simulators [5]. Transports of light fields and ultracold atoms have been widely explored. These two systems also feature the existence of nonlinearity. How nonlinearity affects transport has become an active research field. A spatially localized wave packet of a light beam, arising from the balance between nonlinearity and dispersion, can tunnel through a potential barrier [6]. The nonlinear tunneling of a plane wave through a small potential has been experimentally shown to relate with nonlinearity-induced hysteresis [7]. A wave-packet with an engineered negative effective mass by adjusting nonlinearity has an exotic tunneling behavior through a potential barrier [8]. In both systems of light fields and ultracold atoms, lattice periodic potentials can be easily implemented. For lights, the lattices can be realized with wave-guide arrays [9] or optically induced photonic lattices [10,11]. For ultracold atoms, optical lattices are the standard lattice potentials [12]. Typical transports in lattice potentials are Bloch oscillations [13] and Landau–Zener tunneling [14,15]. After considering nonlinearity, exotic nonlinear transports in lattice potentials include the counter-intuitive self-trapping induced by a repulsive nonlinearity [16,17,18,19]; nonlinear Landau–Zener tunneling, where the adiabaticity is always broken due to nonlinearity no matter how slow the changes of parameters are [20,21,22]; and the breakdown of the superfluidity of nonlinear Bloch waves [23,24,25]. It has been revealed that the nonlinearity can lead to electromagnetically induced transparency amplification [26].
Recently, with the rapid development of non-Hermitian quantum physics, the non-Hermitian periodic potentials with the parity–time reversal (PT) symmetry have attracted great research attention [27]. The PT symmetric lattice potentials generate the Bloch band-gap energy spectrum [28], and all energy may be real-valued [29]. The linear transports of lights in such lattices have been revealed to exhibit exotic features stemming from the non-orthogonality of the Bloch waves [28]. There is a spontaneous phase transition between the PT unbroken and broken phases that can be reflected from whether the Bloch spectrum includes or does not include imaginary modes [30]. After the nonlinearity is introduced to the PT lattices, the nonlinearity can change critical values for the spontaneous PT phase transition [31]. Furthermore, nonlinear Bloch waves without linear analogues can form an asymmetric loop in the Bloch gaps, which plays a fundamental role in the PT phase transition [32]. Bloch oscillations and other transports in the PT symmetric lattice potentials behave differently than in the ordinary lattice potentials [33]. The nonlinear PT symmetric lattice potentials can support the existence of novel solitonic solutions [34]. A great amount of research attention has been paid to the study of the stable transports of these solitons in nonlinear PT symmetric lattices [35,36,37,38].
Here, we studied the nonlinear transports of a light beam through finite PT symmetric lattice potentials. Our study was motivated by possibly experimental realizations of non-Hermitian potentials [39,40,41,42,43] and the existence of nonlinearity in optical systems. However, our results can be easily extended to a similar transport of an ultracold atomic gas in non-Hermitian lattice potentials [16]. Our studied system is similar to that with finite Hermitian lattice potentials in Ref. [16]. It includes two reservoirs behaving as the source and sink for light field and finite PT symmetric lattice potentials placed between the two reservoirs. The initial light field is trapped in one of the reservoirs as the ground state. Once it is released, it shall be allowed to expand toward the lattice potentials due to the kinetic energy that is converted from the nonlinearity. The corresponding infinite lattice potentials feature a nonlinear Bloch band-gap spectrum. Depending on whether the initial spectrum of the light field trapped in the source reservoir is resonant with eigen modes in the nonlinear band or not, the tunneling through the finite lattice potentials can be categorized into two different types. They are in-band transport and in-gap transport. We identified features of these two transports with particular attention to the nonlinear effects. The nonlinear Bloch spectrum of the infinite lattice potentials may have a loop structure that is a pure nonlinear effect. We address that an in-loop transport is an in-band transport.
The rest of this paper is organized as follows. In Section 2, we model our transport system and provide a theoretical frame for the study of nonlinear transport in this system. In Section 3, we show in-band and in-gap transports and identify all features for these two different transports. We describe that the PT symmetrical lattice potentials provide a tunable platform for studying nonlinear transport. Finally, we demonstrate that an in-loop transport belongs to the in-band transport. The conclusion is presented in Section 4.

2. Model

The propagation of light waves in a medium with refractive index n can be described by a nonlinear paraxial equation [28,30,31,33,34]:
i k E Z = 1 2 2 E X 2 + k 0 2 n 0 Δ n E .
Here, E describes the electric field envelope of the light field [28,34], and Z and X are propagation and transverse distances. k 0 = 2 π / λ is the wave vector of the light beam in free space with λ being the wavelength. n 0 is the uniform refractive index of the medium, and k = k 0 n 0 . The total refractive index n is the sum of n 0 and the refractive index change Δ n , i.e., n = n 0 + Δ n . The refractive index change Δ n includes the optically induced photonic potential and the Kerr nonlinearity [31]:
Δ n = V pot + 1 2 α 2 | E | 2 ,
where the potential is categorized into three parts:
V pot = 0 for L A L X < L , V lat ( X ) for L X L , 0 for L < X L + L B .
The scheme of this potential is shown in Figure 1. The lattice potential V lat only exists in a finite region of L X L ; therefore, the lattice is finite, which supports us to study the transport through this finite lattice potential.
V lat ( X ) = n R cos 2 ( π D X ) + i n I sin ( 2 π D X ) ,
Here, n R and n I describe modulation magnitudes of the refraction index by the real and imaginary lattices, respectively, and D relates to the period ofthe lattice potentials [28]. These lattice potentials possess the PT symmetry, V lat * ( X ) = V lat ( X ) . In Equation (2), the Kerr nonlinearity with the first order is represented by the coefficient α 2 , and when α 2 > 0 , the nonlinearity is self-defocusing, and it is self-focusing if α 2 < 0 . In this study, we mainly considered the self-defocusing nonlinearity. For the convenience of numerical calculations and resultant physics analysis, we made the nonlinear equation in Equation (1) dimensionless. It became [31,32]
i ψ z = 1 2 2 ψ x 2 + v lat ( x ) ψ + g | ψ 2 | ψ .
Here, the dimensionless quantities are as follows: x = π D X and z = π 2 k D 2 Z . Therefore, the units of the transverse distance x and the propagation distance z are D / π and k D 2 / π 2 , respectively. ψ = D / ( P 0 π ) E with P 0 is the total power of the incident light, i.e., | E | 2 d X = P 0 . The Kerr nonlinear coefficient is
g = P 0 k 0 2 D n 0 α 2 2 π .
With these units, the field amplitude satisfies the normalization condition | ψ | 2 d x = 1 . The PT symmetric lattice potentials become
v lat ( x ) = V R cos 2 ( x ) + i V I sin ( 2 x ) ,
in the region of π L / D x π L / D . The depths of the real and imaginary lattices are V R = k 0 2 D 2 n 0 n R / π 2 and V I = k 0 2 D 2 n 0 n I / π 2 .
If the PT symmetric lattice potentials are infinite, then the light fields are Bloch waves [28]:
ψ ( x , z ) = e i k x + i b z u ( x ) .
Here, k is the quasimomentum, b is the propagation constant (which corresponds to the chemical potential in ultracold atomic systems), and u ( x ) is a periodic function satisfying u ( x + π ) = u ( x ) . From the Bloch-wave function, we define the spectrum β (which corresponds to the energy in ultracold atomic systems) [32] as follows:
β = d x 1 2 ψ x 2 + v lat ( x ) | ψ | 2 + g 2 | ψ | 4 .
It is noted that the spectrum β is the propagation constant b in the absence of nonlinearity. Generally, they are not equal in the presence of nonlinearity. The relationship β ( k ) forms the Bloch band-gap spectrum. In the absence of nonlinearity, g = 0 , when V I < V R / 2 , all Bloch spectra are real-valued. Correspondingly, the system is said to be in the PT unbroken phase, and when V I > V R / 2 , some imaginary modes appear in the first two bands around the Brillouin zone edges, so the system is in the PT broken phase [28]. V I = V R / 2 is the critical condition for the phase transition between the PT unbroken and broken phases. The transport of light beams in such linear lattice potentials has been studied [28]. In considering the nonlinearity, g 0 , the infinite PT symmetric lattices can still support nonlinear Bloch waves; therefore, the resultant spectrum β ( k ) is a nonlinear Bloch spectrum. Similar to the linear case, whether the nonlinear Bloch spectrum includes imaginary modes can be used to identify the phase transition between the PT unbroken and broken phases [31]. It has been found that the nonlinearity can modify the critical condition for the phase transition, which means that it does not necessarily happen at V I = V R / 2 as in the linear case [31]. Furthermore, there may exist exotic nonlinear Bloch waves that do not have linear analogues. They are located in the energy gaps of the nonlinear Bloch spectrum, forming a symmetric loop structure. These solutions dramatically change the critical condition for the phase transition [32]. The propagation of nonlinear Bloch waves may be destroyed due to dynamical instability [32].
We studied the transport of light beams through the finite PT symmetric lattice potentials. The initial light field is prepared in reservoir A, the ground state [see Figure 1]. This ground state is numerically obtained by calculating the nonlinear equation in Equation (5) with a box potential using the imaginary-time evolution method. Once the initial state is prepared, the right side wall (represented by the dashed vertical line in Figure 1) is removed. The initial field is allowed to expand toward the lattice due to the drivings originating from the release from the reservoir trap and kinetic energy that is converted from the nonlinearity as well. The evolution of the light field is obtained through the numerical calculation of the nonlinear equation in Equation (5).
When the imaginary lattices dominate over the real lattices, the nonlinear spectrum shall include imaginary-energy excitations; therefore, the system is in the PT broken phase. When the light field expands into the lattice region, the imaginary-energy excitations would be excited, leading to the exponential grow of the light field. The field always gains energies from the non-Hermitian lattices and evolves divergently. In the following, we shall study the nonlinear transport of the light field in the PT unbroken phase.

3. Nonlinear Transport

3.1. An In-Band Transport

A typical nonlinear transport is shown in the upper panel in Figure 2 with the parameters g = 198 , V R = 0.28 , and V I = 0.04 . After the light field expands into the lattice potentials, we can define the total power inside the lattices P lat as [16]
P lat ( z ) = L L d x | ψ ( x , z ) | 2 .
The evolution of P lat is recorded in Figure 2(a1). The full distributions | ψ ( x , z ) | 2 around the lattice region of L < x < L are shown in Figure 2(b1). As the light propagates along the z direction, a fraction expands into the lattices [see the red line in Figure 2(c1)]; as a result, P lat increases starting from 0. Soon, it reaches a plateau when z > 400 . This is because before the plateau, the light field starts to tunnel outside the lattices into reservoir B. Around the plateau, the continuous expansion into the lattice region from reservoir A is balanced with leakage on the right side, leading to a constant total power inside the lattice region. As z increases, the other plateau in P lat appears. Between these two plateaus, there is a fast increase. The reason for the discontinuity of the two plateaus is the reflections of the light field around the end of the lattices x = L . As shown in Figure 2(b1), for a very small z, when the light field reaches the end of the lattices, it has the possibility to be reflected back, and the reflected parts travel oppositely toward reservoir A. The reflected parts are interfering with forward-moving parts, causing the dark stripes that generate around the lattice end and travel back. It is noticed that the reflections only happen when the light field is very weak around the lattice end. Figure 2(b1) shows that there is no reflection of a strong light field around the lattice end. In the first plateau of P lat , i.e., 400 < z < 800 [see Figure 2(a1)], there is one dark interference stripe moving backward. In contrast, in the second plateau, z > 800 , there is no reflection [see Figure 2(b1)]. Further, in propagating z, P lat starts to decrease around z = 1300 . This is because the main part of the light field has been transported through the lattice region. From Figure 2(b1), it can be seen that the field on the left side of the lattice region has a dramatic decrease around z = 1500 . The typical distribution (represented by the purple solid line in Figure 2(c1)) also indicates the decrease in the field on the left side of the lattice region. These results indicate that the light field left inside reservoir A becomes weak.
In order to show the nonlinear transport more clearly, we calculate the power flow j as follows [34]:
j ( x , z ) = 1 2 i ψ * x ψ ψ x ψ * .
In the interesting study in [44], an alternative definition of power flow is provided based on a continuity equation from the non-Hermitian dynamical equation. We use the definition in Equation (11), which does not obey the continuity equation; however, it is an experimentally measurable quantity. The power flow at the end of the lattices x = L ,
j end = j ( x = L , z ) ,
is shown in Figure 2(d1). As z increases, j end starts from 0, quickly reaches a steady state, and has a very small and fast oscillation around the steady value (represented by the horizontal dashed line in the figure). The steady value can be explicitly obtained by calculating the change in the power inside reservoir B:
P B = L L + L B d x | ψ | 2 , j ste = d P B d z z = z pla ,
where P B is the power inside reservoir B, and j ste is the change in the power P B at the propagation distance z pla , where the plateaus in P lat in Figure 2(a1) appear. Physically, in the existing regions of the plateaus, the input into the lattices is equal to the output at the end of the lattices. In considering the continuity of the power flow, the steady value calculated from reservoir B naturally reflects the output at the end of the lattices. For a sufficiently large z, the power flow j end decreases from the steady value. This region corresponds to the decreasing of P lat after the plateaus, indicating that the main part of the light field transports through the lattices.
The nonlinear transport shown in the upper panel in Figure 2 indicates an in-band transport. In the plateaus, P lat keeps a constant value; therefore, the average power inside the lattices is P lat / ( 2 L ) . If we assume that the PT symmetric lattices are infinite, the resultant Bloch waves should be normalized over a unit cell as 0 π d x | ψ | 2 / π = 1 . As a result, the light field should be re-normalized | ψ | 2 | ψ | 2 P lat / ( 2 L ) in order to make the Bloch waves re-normalized into 1 over a unit cell of the infinite lattices [45]. Consequently, the nonlinearity coefficient is also re-normalized as g g lat with
g lat = g P lat / ( 2 L ) .
Therefore, the Bloch waves of the infinite lattices are described by the nonlinear equations in Equation (5) by replacing g with g lat ; in this way, the Bloch waves have the average power P lat / ( 2 L ) and are normalized into 1 over a unit cell. With these re-normalized quantities and the re-normalization condition over a unit cell, we numerically solve the nonlinear equation in Equation (5) with the infinite PT symmetric lattice potentials and obtain the nonlinear Bloch-wave solutions. At last, we have the nonlinear Bloch spectrum β ( k ) from Equation (9). We adopt the value of P lat in the second plateau in Figure 2(a1), and the re-normalized g lat becomes g lat = 0.198 . The calculated first two nonlinear Bloch bands are depicted in Figure 2(e1). These two nonlinear bands have a gap opened around the Brillouin zone edge k = 0.5 . The spectrum of the initial field trapped in reservoir A is also shown in this figure as a horizontal dashed line. It can be seen that the initial spectrum just lies inside the first nonlinear band. Physically, we think that the initial spectrum is resonant with the eigen modes of the finite lattice potentials, and the initial light field is allowed to tunnel through the lattice potentials [16]. The decreases of P lat in Figure 2(a1) and j end in Figure 2(d1) for a sufficiently large z provide signatures that the light wave packet transports through the finite lattices. Therefore, the nonlinear transport belongs to the in-band transport.

3.2. An In-Gap Transport

Another typical nonlinear transport is shown in the lower panel in Figure 2 with the parameters g = 398 , V R = 0.28 , and V I = 0.014 . The total power inside the lattice region P lat increases from 0 and reaches a wide plateau when z > 400 . Before this large plateau, P lat experiences two short plateaus as z increases [see Figure 2(a2)]. The existence of these two short plateaus is because there are two dark interference stripes when z > 200 moving backward, as shown from the power distribution | ψ ( x , t ) | 2 in Figure 2(b2). It is interesting to note that these two stripes are not reflected back from the end of the lattices, but from the inside. The formation of the large plateau z > 400 in P lat is also very different from the in-band transport in Figure 2(a1). As can be seen from Figure 2(b2), in the large plateau region, the power distribution keeps the profile without the change as z increases and only occupies the left side of the lattices. A snapshot distribution in the large plateau is indicated in Figure 2(c2) by the red solid line. This clearly shows that the light wave packet localizes on the left side of the lattices and only a very small tail extends into reservoir B. Such a weak extension gives rise to a small power flow into reservoir B. Figure 2(d2) shows that in the large plateau region of P lat , the power flow at the end of the lattice j end oscillates around a steady value that can be obtained from the information in reservoir B using Equation (13). The steady value is smaller in comparison with the in-band transport in Figure 2(d1). Therefore, the formation of the large plateau is due to the fact that there is no further input or output in the lattices. The light wave packet behaves like a solitary wave trapped inside the lattices as z increases [16]. The distribution profile | ψ ( x ) | 2 [see the red line in Figure 2(c2)] is very similar to the solitary-wave solutions existing in the nonlinear infinite lattices that were originally found in Refs. [45,46]. We emphasize that the localization of the light field is the effect of the defocusing nonlinearity, which is very counter-intuitive since the defocusing nonlinearity usually prefers to expand a wave packet. If there is no nonlinearity, the field that penetrates the lattice region will be transported away and cannot be trapped.
In other words, with the help of nonlinearity, the light field cannot transport through the lattices during propagating in the large plateau z region. However, the propagation of such a solitary wave is still accompanied by a very weak leakage at the right end of the lattices [see Figure 2(d1), i.e., j end 0 ]. After a sufficient, long propagation z, the accumulation of the weak leakage becomes serious and makes the solitary wave destabilize. At this moment, the solitary wave begins to collapse. Most of the power shall be transported through the lattices into reservoir B, causing a sharp increase of j end around z = 1000 , as shown in Figure 2(d2). Due to such an increase of j end , the total power inside the lattices P lat decreases dramatically around z = 1000 [see Figure 2(a2)]. The accumulation of the weak leakage that induced transport is similar to the destabilization of the solitary-wave solutions in the infinite lattices that were studied in Ref. [47]. After the collapse, the light field can tunnel through the lattices, and there are continuous inputs. The behavior of the light field is similar as that in the in-band transport shown in the upper panel in Figure 2. In Figure 2(c2), the black dashed line shows the field distribution at z = 1400 . We can see that the distribution shares many similarities with the in-band transport distribution shown by the black dashed line in Figure 2(c1).
The absence of transport due to the formation of the solitary wave and the sharp increase of j end induced by the destabilization of the solitary wave is a signature of the in-gap transport. We use the value of P lat in the large plateau in Figure 2(a1) and calculate the re-normalized nonlinear coefficient as g lat = 0.38 . Finally, we can obtain the nonlinear Bloch spectrum β ( k ) after calculating the nonlinear Bloch-wave solutions for the infinite lattice potentials with the re-normalized nonlinear coefficient. The first two nonlinear bands are shown in Figure 2(d2). These two bands are similar to those of the in-band transport in Figure 2(d1). However, when we place the spectrum of the initial light field trapped in reservoir A in this figure [see the horizontal dashed line], we find that the initial spectrum lies inside the gap of the nonlinear Bloch spectrum. There is no eigen mode of the finite lattice that can be resonant with the initial spectrum; therefore, the initial light field cannot tunnel through the lattices but forms a solitary wave due to the existence of nonlinearity. The breakdown of the solitary wave can trigger the transport through the lattices. This is because the breakdown reduces g lat . Consequently, the nonlinear spectrum moves upward because of the reduced g lat , and then the initial spectrum can enter into the second band of the nonlinear spectrum, generating an in-band transport.

3.3. Tunable Transport

The presence of the imaginary lattices provides a tunable opportunity to study the nonlinear transport of the light field. It is noticed from Figure 2 that both in-band and in-gap transports feature the steady power flow j ste , which can be calculated from Equation (13). For different values of imaginary lattice amplitudes V I , we repeated realizing transports and calculated the steady values, and we show the results in Figure 3. In Figure 3a,b, the nonlinear coefficients g = 198 and g = 398 are the same as the in-band and in-gap transports shown in the upper and lower panels, respectively, in Figure 2. As shown by the red dashed line in Figure 3a, with the increase of V I , the steady value stays the same as the in-band transport in Figure 2(d1), which indicates that the transport is always an in-band transport. When V I > 0.18 , the steady flow becomes divergent. This is because the large imaginary lattice amplitude breaks the PT symmetry and leads to the PT broken phase. Without the nonlinearity, the critical condition for the PT phase transition is V I = V R / 2 = 0.14 . However, the nonlinearity modifies the critical condition, and from the red dashed line, it can be seen that the critical value of V I is larger than the linear prediction. The solid line in Figure 3a is for V R = 0.45 , and in the window of V I the steady flow is always a constant; therefore, the transport is an in-band case. Figure 3b shows the case of a larger nonlinearity. The results are different from those in Figure 3a. For a very small V I , the red dashed line (representing V R = 0.28 ) has the same value as the in-gap transport in Figure 2(d2). Therefore, for the small V I , the transport is still in-gap transport. As V I increases, the steady value soon increases to a value that is the same as the value of the in-band transport. Here, the transport changes into the in-band transport. Further, increasing V I > 0.19 leads it to start to diverge, which means that the PT symmetry is broken. The solid line represents the case of V R = 0.45 [see Figure 3b]. When V I < 0.11 , the steady value indicates the in-gap transport. In increasing the imaginary lattice amplitude, the steady flow increases and enters a plateau with a larger value; therefore, the corresponding transport goes into the in-band transport.

3.4. An In-Loop Transport

In the presence of nonlinearity, the nonlinear Bloch spectrum may have exotic nonlinear solutions that do not have linear analogues. These nonlinear solutions form a loop structure located inside the gaps [32]. We shall study a transport with an initial spectrum of the light field trapped in reservoir A that is resonant with eigen modes in the loop. Such a transport is named an in-loop transport. A typical result is shown in Figure 4. For the chosen parameters, the initial ground state trapped in reservoir A is allowed to expand into the finite lattice potential. The evolution of the total power inside the lattice region is depicted in Figure 4a. As z increases, the total power P lat increases dramatically to reach a plateau. When z > 900 , the plateau begins to collapse, and P lat starts to decrease. In the plateau region, the propagation of the light field along the transverse direction is stationary. This can also be seen from the power flow j end at the end of the lattices shown in Figure 4b. In the plateau region of P lat , j end is mainly around a steady value (which is shown by the dashed line in the figure). However, beyond this region, when z > 900 , j end starts to decrease. The behavior of j end shown in Figure 4b is similar as that of the in-band transport in Figure 2(d1).
Based on the plateau of P lat , we calculated the re-normalized nonlinear coefficient from Equation (14) and obtained g lat = 0.378 . Then, the nonlinear Bloch-wave solutions of the infinite PT symmetric lattice potentials with the re-normalized nonlinear coefficient were systematically examined. Furthermore, from the nonlinear Bloch-wave solutions, we calculated the nonlinear Bloch spectrum. The first two bands in the nonlinear spectrum are shown in Figure 4c. The outstanding feature of the nonlinear spectrum is that there is an interesting loop structure located inside the gap. The nonlinear Bloch-wave solutions in the loop are exotic since they do not have linear analogues. Considering the initial spectrum (shown by the dashed line in the figure), we realized that the initial spectrum is resonant with the modes in the loop. Since the loop structure adheres to the first band, it may be considered a part of the first band. Consequently, the resonance with the loop structure means that the transport is an in-band transport. The transport results shown in Figure 4a,b are similar to those of the in-band transport shown in the upper panel in Figure 2, which is an indication that the in-loop transport essentially becomes the in-band transport.

4. Conclusions

We studied the nonlinear transport of light through finite PT symmetric lattice potentials. The lattice potentials included the real and imaginary lattices. The initial light was trapped in reservoir A as a ground state. It was allowed to expand into the lattice potentials once it was released from reservoir A. The expansion along the transverse direction was derived from the kinetic energy that was converted from the nonlinear energy. After tunneling through the lattice potentials, reservoir B behaved as a sink to hold the transported light field.
We found that the nonlinear transports can be categorized into two types: in-band and in-gap transports. The lattice potentials feature a nonlinear Bloch band-gap spectrum. The in-band transport happens when the initial spectrum of the light field trapped in reservoir A is resonant with eigen modes inside the nonlinear band. The in-band transport is characterized by the existence of a stationary period, during which the input to the lattice region from reservoir A is balanced with the output from the lattice region to reservoir B. After this period, the main part of the light field can tunnel through the lattice potential, and then the light field received by reservoir B will be reduced. The in-gap transport happens when the initial spectrum of the light field trapped in reservoir A lies inside the nonlinear gap of the lattice potentials. The in-gap transport is also characterized by the existence of a stationary period. But the origination of such a stationary period is totally different from that of the in-band transport. This is because the light field cannot tunnel through the lattice potentials but forms a solitary wave localized inside the lattice region. The formed solitary wave can propagate for a long distance leading to the stationary period. However, the propagation of the solitary wave is still accompanied by a very weak leakage into reservoir B. The accumulation of the weak leakage shall finally break the solitary wave and restart the tunneling into reservoir B. The nonlinear Bloch spectrum of the lattice potentials can feature a loop structure that adheres to the first nonlinear band. We identified that the in-loop transport that happens when the initial spectrum of the light field trapped in reservoir A is resonant with modes in the loop belongs to an in-band transport.

Author Contributions

Conceptualization, W.M. and Y.Z.; methodology, Y.Z.; formal analysis, W.M. and Y.Z.; investigation, W.M. and Y.Z.; writing—original draft preparation, W.M.; writing—review and editing, W.M. and Y.Z.; visualization, W.M.; supervision, Y.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China with Grants No. 12374247 and No. 11974235.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. de Aguiar, F.M.; Wharam, D.A. Transport through one-dimensional channels. Phys. Rev. B 1991, 43, 9984–9987. [Google Scholar] [CrossRef] [PubMed]
  2. Shimshoni, E.; Auerbach, A.; Kapitulnik, A. Transport through Quantum Melts. Phys. Rev. Lett. 1998, 80, 3352–3355. [Google Scholar] [CrossRef]
  3. Stangl, J.; Holý, V.; Bauer, G. Structural properties of self-organized semiconductor nanostructures. Rev. Mod. Phys. 2004, 76, 725–783. [Google Scholar] [CrossRef]
  4. Zwanenburg, F.A.; Dzurak, A.S.; Morello, A.; Simmons, M.Y.; Hollenberg, L.C.L.; Klimeck, G.; Rogge, S.; Coppersmith, S.N.; Eriksson, M.A. Silicon quantum electronics. Rev. Mod. Phys. 2013, 85, 961–1019. [Google Scholar] [CrossRef]
  5. Georgescu, I.M.; Ashhab, S.; Nori, F. Quantum simulation. Rev. Mod. Phys. 2014, 86, 153–185. [Google Scholar] [CrossRef]
  6. Barak, A.; Peleg, O.; Stucchio, C.; Soffer, A.; Segev, M. Observation of Soliton Tunneling Phenomena and Soliton Ejection. Phys. Rev. Lett. 2008, 100, 153901. [Google Scholar] [CrossRef] [PubMed]
  7. Wan, W.; Muenzel, S.; Fleischer, J.W. Wave Tunneling and Hysteresis in Nonlinear Junctions. Phys. Rev. Lett. 2010, 104, 073903. [Google Scholar] [CrossRef] [PubMed]
  8. Di Mei, F.; Caramazza, P.; Pierangeli, D.; Di Domenico, G.; Ilan, H.; Agranat, A.J.; Di Porto, P.; DelRe, E. Intrinsic Negative Mass from Nonlinearity. Phys. Rev. Lett. 2016, 116, 153902. [Google Scholar] [CrossRef] [PubMed]
  9. Christodoulides, D.N.; Lederer, F.; Silberberg, Y. Discretizing light behaviour in linear and nonlinear waveguide lattices. Nature 2003, 424, 817–823. [Google Scholar] [CrossRef]
  10. Fleischer, J.W.; Segev, M.; Efremidis, N.K.; Christodoulides, D.N. Observation of two-dimensional discrete solitons in optically induced nonlinear photonic lattices. Nature 2003, 422, 147–150. [Google Scholar] [CrossRef]
  11. Lederer, F.; Stegeman, G.I.; Christodoulides, D.N.; Assanto, G.; Segev, M.; Silberberg, Y. Discrete solitons in optics. Phys. Rep. 2008, 463, 1–126. [Google Scholar] [CrossRef]
  12. Morsch, O.; Oberthaler, M. Dynamics of Bose-Einstein condensates in optical lattices. Rev. Mod. Phys. 2006, 78, 179–215. [Google Scholar] [CrossRef]
  13. Ben Dahan, M.; Peik, E.; Reichel, J.; Castin, Y.; Salomon, C. Bloch Oscillations of Atoms in an Optical Potential. Phys. Rev. Lett. 1996, 76, 4508–4511. [Google Scholar] [CrossRef] [PubMed]
  14. Niu, Q.; Zhao, X.G.; Georgakis, G.A.; Raizen, M.G. Atomic Landau-Zener Tunneling and Wannier-Stark Ladders in Optical Potentials. Phys. Rev. Lett. 1996, 76, 4504–4507. [Google Scholar] [CrossRef] [PubMed]
  15. Kishimoto, T.; Hachisu, H.; Fujiki, J.; Nagato, K.; Yasuda, M.; Katori, H. Electrodynamic Trapping of Spinless Neutral Atoms with an Atom Chip. Phys. Rev. Lett. 2006, 96, 123001. [Google Scholar] [CrossRef] [PubMed]
  16. Rosenkranz, M.; Jaksch, D.; Yin Lim, F.; Bao, W. Self-trapping of Bose-Einstein condensates expanding into shallow optical lattices. Phys. Rev. A 2008, 77, 063607. [Google Scholar] [CrossRef]
  17. Anker, T.; Albiez, M.; Gati, R.; Hunsmann, S.; Eiermann, B.; Trombettoni, A.; Oberthaler, M.K. Nonlinear Self-Trapping of Matter Waves in Periodic Potentials. Phys. Rev. Lett. 2005, 94, 020403. [Google Scholar] [CrossRef]
  18. Trombettoni, A.; Smerzi, A. Discrete Solitons and Breathers with Dilute Bose-Einstein Condensates. Phys. Rev. Lett. 2001, 86, 2353–2356. [Google Scholar] [CrossRef]
  19. Elenewski, J.E.; Chen, H. Real-time transport in open quantum systems from P T -symmetric quantum mechanics. Phys. Rev. B 2014, 90, 085104. [Google Scholar] [CrossRef]
  20. Wu, B.; Niu, Q. Nonlinear Landau-Zener tunneling. Phys. Rev. A 2000, 61, 023402. [Google Scholar] [CrossRef]
  21. Chen, Y.A.; Huber, S.D.; Trotzky, S.; Bloch, I.; Altman, E. Many-body Landau–Zener dynamics in coupled one-dimensional Bose liquids. Nat. Phys. 2011, 7, 61–67. [Google Scholar] [CrossRef]
  22. Guan, Q.; Ome, M.K.H.; Bersano, T.M.; Mossman, S.; Engels, P.; Blume, D. Nonexponential Tunneling due to Mean-Field-Induced Swallowtails. Phys. Rev. Lett. 2020, 125, 213401. [Google Scholar] [CrossRef] [PubMed]
  23. Wu, B.; Niu, Q. Landau and dynamical instabilities of the superflow of Bose-Einstein condensates in optical lattices. Phys. Rev. A 2001, 64, 061603. [Google Scholar] [CrossRef]
  24. Fallani, L.; De Sarlo, L.; Lye, J.E.; Modugno, M.; Saers, R.; Fort, C.; Inguscio, M. Observation of Dynamical Instability for a Bose-Einstein Condensate in a Moving 1D Optical Lattice. Phys. Rev. Lett. 2004, 93, 140406. [Google Scholar] [CrossRef] [PubMed]
  25. Hamner, C.; Zhang, Y.; Khamehchi, M.A.; Davis, M.J.; Engels, P. Spin-Orbit-Coupled Bose-Einstein Condensates in a One-Dimensional Optical Lattice. Phys. Rev. Lett. 2015, 114, 070401. [Google Scholar] [CrossRef] [PubMed]
  26. Iqbal, A.; Khan, N.; Bacha, B.A.; Rahman, A.U.; Ahmad, A. Photon drag enhancement by a slow-light moving medium via electromagnetically-induced transparency amplification. Phys. Lett. A 2017, 381, 3134–3140. [Google Scholar] [CrossRef]
  27. Konotop, V.V.; Yang, J.; Zezyulin, D.A. Nonlinear waves in P T -symmetric systems. Rev. Mod. Phys. 2016, 88, 035002. [Google Scholar] [CrossRef]
  28. Makris, K.G.; El-Ganainy, R.; Christodoulides, D.N.; Musslimani, Z.H. Beam Dynamics in P T Symmetric Optical Lattices. Phys. Rev. Lett. 2008, 100, 103904. [Google Scholar] [CrossRef] [PubMed]
  29. Bender, C.M.; Brody, D.C.; Jones, H.F. Complex Extension of Quantum Mechanics. Phys. Rev. Lett. 2002, 89, 270401. [Google Scholar] [CrossRef]
  30. Makris, K.G.; El-Ganainy, R.; Christodoulides, D.N.; Musslimani, Z.H. P T -symmetric optical lattices. Phys. Rev. A 2010, 81, 063807. [Google Scholar] [CrossRef]
  31. Lumer, Y.; Plotnik, Y.; Rechtsman, M.C.; Segev, M. Nonlinearly Induced PT Transition in Photonic Systems. Phys. Rev. Lett. 2013, 111, 263901. [Google Scholar] [CrossRef]
  32. Zhang, Y.; Chen, Z.; Wu, B.; Busch, T.; Konotop, V.V. Asymmetric Loop Spectra and Unbroken Phase Protection due to Nonlinearities in PT-Symmetric Periodic Potentials. Phys. Rev. Lett. 2021, 127, 034101. [Google Scholar] [CrossRef] [PubMed]
  33. Longhi, S. Bloch Oscillations in Complex Crystals with P T Symmetry. Phys. Rev. Lett. 2009, 103, 123601. [Google Scholar] [CrossRef] [PubMed]
  34. Musslimani, Z.H.; Makris, K.G.; El-Ganainy, R.; Christodoulides, D.N. Optical Solitons in P T Periodic Potentials. Phys. Rev. Lett. 2008, 100, 030402. [Google Scholar] [CrossRef]
  35. Chen, Y.; Yan, Z.; Li, X. One- and two-dimensional gap solitons and dynamics in the PT-symmetric lattice potential and spatially-periodic momentum modulation. Commun. Nonlinear Sci. Numer. Simul. 2018, 55, 287–297. [Google Scholar] [CrossRef]
  36. Tamilselvan, K.; Govindarajan, A.; Inbavalli, I.; Alagesan, T.; Lakshmanan, M. Modulational instability in P T -symmetric Bragg grating structures with saturable nonlinearity. Phys. Rev. A 2023, 107, 053510. [Google Scholar] [CrossRef]
  37. Raja, S.V.; Govindarajan, A.; Lakshmanan, M. Ultra low-power multistability in PT-symmetric periodic structures with saturable coupling, gain and loss. Phys. Lett. A 2024, 509, 129517. [Google Scholar] [CrossRef]
  38. Raja, S.V.; Govindarajan, A.; Mahalingam, A.; Lakshmanan, M. Tailoring inhomogeneous P T -symmetric fiber-Bragg-grating spectra. Phys. Rev. A 2020, 101, 033814. [Google Scholar] [CrossRef]
  39. Guo, A.; Salamo, G.J.; Duchesne, D.; Morandotti, R.; Volatier-Ravat, M.; Aimez, V.; Siviloglou, G.A.; Christodoulides, D.N. Observation of P T -Symmetry Breaking in Complex Optical Potentials. Phys. Rev. Lett. 2009, 103, 093902. [Google Scholar] [CrossRef]
  40. Kottos, T. Broken symmetry makes light work. Nat. Phys. 2010, 6, 166–167. [Google Scholar] [CrossRef]
  41. Rüter, C.E.; Makris, K.G.; El-Ganainy, R.; Christodoulides, D.N.; Segev, M.; Kip, D. Observation of parity–time symmetry in optics. Nat. Phys. 2010, 6, 192–195. [Google Scholar] [CrossRef]
  42. Lin, Z.; Ramezani, H.; Eichelkraut, T.; Kottos, T.; Cao, H.; Christodoulides, D.N. Unidirectional Invisibility Induced by P T -Symmetric Periodic Structures. Phys. Rev. Lett. 2011, 106, 213901. [Google Scholar] [CrossRef] [PubMed]
  43. Sun, Y.; Tan, W.; Li, H.q.; Li, J.; Chen, H. Experimental Demonstration of a Coherent Perfect Absorber with PT Phase Transition. Phys. Rev. Lett. 2014, 112, 143903. [Google Scholar] [CrossRef] [PubMed]
  44. Bagchi, B.; Quesne, C.; Znojil, M. Generalized continuity Equation and modified normalization in pt-symmetric quantum mechanics. Mod. Phys. Lett. A 2001, 16, 2047–2057. [Google Scholar] [CrossRef]
  45. Zhang, Y.; Wu, B. Composition Relation between Gap Solitons and Bloch Waves in Nonlinear Periodic Systems. Phys. Rev. Lett. 2009, 102, 093905. [Google Scholar] [CrossRef] [PubMed]
  46. Alexander, T.J.; Ostrovskaya, E.A.; Kivshar, Y.S. Self-Trapped Nonlinear Matter Waves in Periodic Potentials. Phys. Rev. Lett. 2006, 96, 040401. [Google Scholar] [CrossRef]
  47. Wang, B.; Fu, P.; Liu, J.; Wu, B. Self-trapping of Bose-Einstein condensates in optical lattices. Phys. Rev. A 2006, 74, 063610. [Google Scholar] [CrossRef]
Figure 1. Schematic setup of the transport of a light field through a finite PT symmetric lattice potential. The transverse direction is along x. The light field is initially trapped in reservoir A. Once the wall (represented by a dashed vertical line) is removed, the light expands along the transverse direction x to tunnel through the lattice potential that is located in the region of L < x < L . Reservoir B behaves as a sink to hold the light field that tunnels through the lattice potential. The parameters we use are L A = 145 π ,   L B = 326 π , and L = 5 π (the unit of these distances is D / π , with D being the period of the lattice potential).
Figure 1. Schematic setup of the transport of a light field through a finite PT symmetric lattice potential. The transverse direction is along x. The light field is initially trapped in reservoir A. Once the wall (represented by a dashed vertical line) is removed, the light expands along the transverse direction x to tunnel through the lattice potential that is located in the region of L < x < L . Reservoir B behaves as a sink to hold the light field that tunnels through the lattice potential. The parameters we use are L A = 145 π ,   L B = 326 π , and L = 5 π (the unit of these distances is D / π , with D being the period of the lattice potential).
Symmetry 16 00640 g001
Figure 2. Two different nonlinear transports through the finite PT lattice potentials. The upper panel shows the in-band transport with parameters g = 198 , V R = 0.28 , and V I = 0.04 , and the lower panel shows the in-gap transport with parameters g = 398 , V R = 0.28 , and V I = 0.014 . (a1,a2) The total power inside the lattices P lat as a function of the propagation distance z. (b1,b2) Evolution of the power distribution | ψ ( x , z ) | 2 around the lattice region of L < x < L with L = 5 π . Bright yellow indicates high power, and blue indicates low power. (c1,c2) The power distributions | ψ ( x ) | 2 around the lattice region at particular distances z, which are represented by the horizontal lines in (b1,b2). (d1,d2) The power flow at the lend of the lattices j end as a function of the propagation distance. The dashed horizontal lines indicate the steady values calculated from Equation (13). (e1,e2) Nonlinear Bloch spectra β ( k ) of the infinite PT symmetric lattice potentials with the re-normalized nonlinearity coefficients g lat = 0.198 and g lat = 0.38 , respectively. The horizontal dashed lines represent the spectrum of the initial light field trapped in reservoir A.
Figure 2. Two different nonlinear transports through the finite PT lattice potentials. The upper panel shows the in-band transport with parameters g = 198 , V R = 0.28 , and V I = 0.04 , and the lower panel shows the in-gap transport with parameters g = 398 , V R = 0.28 , and V I = 0.014 . (a1,a2) The total power inside the lattices P lat as a function of the propagation distance z. (b1,b2) Evolution of the power distribution | ψ ( x , z ) | 2 around the lattice region of L < x < L with L = 5 π . Bright yellow indicates high power, and blue indicates low power. (c1,c2) The power distributions | ψ ( x ) | 2 around the lattice region at particular distances z, which are represented by the horizontal lines in (b1,b2). (d1,d2) The power flow at the lend of the lattices j end as a function of the propagation distance. The dashed horizontal lines indicate the steady values calculated from Equation (13). (e1,e2) Nonlinear Bloch spectra β ( k ) of the infinite PT symmetric lattice potentials with the re-normalized nonlinearity coefficients g lat = 0.198 and g lat = 0.38 , respectively. The horizontal dashed lines represent the spectrum of the initial light field trapped in reservoir A.
Symmetry 16 00640 g002
Figure 3. Tunable nonlinear transports represented by the steady power flow in Equation (13) as a function of the imaginary lattice amplitude V I . (a) The nonlinear coefficient g = 198 , which is the same as that of the in-band transport shown in the upper panel in Figure 2. (b) The nonlinear coefficient g = 398 , which is the same as that of the in-gap transport shown in the lower panel in Figure 2. In both (a) and (b), the red dashed lines indicate the real lattice amplitude V R = 0.28 , and the blue solid lines indicate V R = 0.45 .
Figure 3. Tunable nonlinear transports represented by the steady power flow in Equation (13) as a function of the imaginary lattice amplitude V I . (a) The nonlinear coefficient g = 198 , which is the same as that of the in-band transport shown in the upper panel in Figure 2. (b) The nonlinear coefficient g = 398 , which is the same as that of the in-gap transport shown in the lower panel in Figure 2. In both (a) and (b), the red dashed lines indicate the real lattice amplitude V R = 0.28 , and the blue solid lines indicate V R = 0.45 .
Symmetry 16 00640 g003
Figure 4. An in-loop transport with the parameters g = 383 , V R = 0.28 , and V I = 0.066 . (a) The total power inside the lattice region as a function of the propagation distance z. (b) The power flow at the end of the lattices. The horizontal dashed line is the steady value calculated from Equation (13). (c) The nonlinear Bloch spectrum β ( k ) of the infinite PT symmetric lattice potential with the re-normalized nonlinear coefficient g lat = 0.378 . The horizontal dashed line represents the initial spectrum of the light field trapped in reservoir A.
Figure 4. An in-loop transport with the parameters g = 383 , V R = 0.28 , and V I = 0.066 . (a) The total power inside the lattice region as a function of the propagation distance z. (b) The power flow at the end of the lattices. The horizontal dashed line is the steady value calculated from Equation (13). (c) The nonlinear Bloch spectrum β ( k ) of the infinite PT symmetric lattice potential with the re-normalized nonlinear coefficient g lat = 0.378 . The horizontal dashed line represents the initial spectrum of the light field trapped in reservoir A.
Symmetry 16 00640 g004
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mao, W.; Zhang, Y. Nonlinear Transport through Parity–Time Symmetric Lattice Potentials. Symmetry 2024, 16, 640. https://doi.org/10.3390/sym16060640

AMA Style

Mao W, Zhang Y. Nonlinear Transport through Parity–Time Symmetric Lattice Potentials. Symmetry. 2024; 16(6):640. https://doi.org/10.3390/sym16060640

Chicago/Turabian Style

Mao, Wei, and Yongping Zhang. 2024. "Nonlinear Transport through Parity–Time Symmetric Lattice Potentials" Symmetry 16, no. 6: 640. https://doi.org/10.3390/sym16060640

APA Style

Mao, W., & Zhang, Y. (2024). Nonlinear Transport through Parity–Time Symmetric Lattice Potentials. Symmetry, 16(6), 640. https://doi.org/10.3390/sym16060640

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop