Acta Polytechnica doi:10.14311/AP.2014.54.0106 Acta Polytechnica 54(2):106–112, 2014 © Czech Technical University in Prague, 2014 available online at http://ojs.cvut.cz/ojs/index.php/ap EXCEPTIONAL POINTS IN OPEN AND PT-SYMMETRIC SYSTEMS Hichem Eleucha, Ingrid Rotterb, ∗ a Department of Physics, McGill University, Montreal, Canada H3A 2T8 b Max Planck Institute for the Physics of Complex Systems, D-01187 Dresden, Germany ∗ corresponding author: rotter@pks.mpg.de Abstract. Exceptional points (EPs) determine the dynamics of open quantum systems and cause also PT symmetry breaking in PT symmetric systems. From a mathematical point of view, this is caused by the fact that the phases of the wavefunctions (eigenfunctions of a non-Hermitian Hamiltonian) relative to one another are not rigid when an EP is approached. The system is therefore able to align with the environment to which it is coupled and, consequently, rigorous changes of the system properties may occur. We compare analytically as well as numerically the eigenvalues and eigenfunctions of a 2 × 2 matrix that is characteristic either of open quantum systems at high level density or of PT symmetric optical lattices. In both cases, the results show clearly the influence of the environment on the system in the neighborhood of EPs. Although the systems are very different from one another, the eigenvalues and eigenfunctions indicate the same characteristic features. Keywords: exceptional points, open quantum systems; PT symmetry breaking, dynamical phase transition, non-Hermitian quantum physics, phase rigidity. 1. Introduction Starting with paper [1], it has been shown that a wide class of PT symmetric non-Hermitian Hamilton operators provides entirely real spectra. In the follow- ing years this phenomenon has been studied in many theoretical papers, see the review [2] and the Special Issue [3]. In order to realize complex PT symmetric struc- tures, the formal equivalence of the quantum mechani- cal Schrödinger equation to the optical wave equation in PT symmetric optical lattices [4] can be exploited by involving symmetric index guiding and an anti- symmetric gain/loss profile. Experimental results [5] have confirmed the expectations and have, further- more, demonstrated the onset of passive PT symme- try breaking within the context of optics. This phase transition was found to lead to a loss-induced optical transparency in specially designed pseudo-Hermitian potentials. In another experiment [6], the wave propa- gation in an active PT symmetric coupled waveguide system is studied. Both spontaneous PT symmetry breaking and power oscillations violating left-right symmetry are observed. Moreover, the relation of the relative phases of the eigenstates of the system to their distance from the level crossing point is obtained. The phase transition occurs when this point is approached. The meaning of these results for a new generation of integrated photonic devices is discussed in [7]. Today we have many experimental and theoretical studies related to this topic. On the other hand, non-Hermitian operators are known to describe open quantum systems in a natural manner, see, e.g., [8]. In contrast to the original pa- pers more than 50 years ago, statistical assumptions on the system’s states are not at all necessary today [9] due to the improved accuracy of the experimen- tal as well as theoretical studies. In the present-day papers, the system is assumed to be open due to the fact that it is embedded into the continuum of scatter- ing wavefunctions into which the states of the system can decay. This environment exists always. It can be changed by means of external forces, but cannot be deleted [10]. The states of the system can decay due to their coupling to the environment of scatter- ing wavefunctions but cannot be formed out of the continuum. Hence, the loss is usually nonvanishing, while the gain is zero. The complex eigenvalues of the non-Hermitian Hamiltonian provide both the energy Ei as well as the lifetime τi (inverse proportional to the decay width Γi) of the eigenstate i. Recent studies have shown the important role the singular points in the continuum play for the dynam- ics of open quantum systems, see, e.g., the review [10]. These singular points are usually called exceptional points (EPs) after Kato, who studied their mathemat- ical properties [11] many years ago. The relation of EPs to PT symmetry breaking in optical systems is considered already in the first papers [6, 7]. Neverthe- less, the relation between the dynamical properties of open quantum systems and those of PT symmetric systems has not been considered thoroughly up to now. It is the aim of the present paper to compare di- rectly the influence of EPs onto the dynamics of open quantum systems with that onto PT symmetry break- ing in PT symmetric systems. The comparison is performed on the basis of simple models with only 106 http://dx.doi.org/10.14311/AP.2014.54.0106 http://ojs.cvut.cz/ojs/index.php/ap vol. 54 no. 2/2014 Exceptional Points in Open and PT-Symmetric Systems two levels coupled to one common channel. In both cases, the Hamiltonian is given by a 2 × 2 matrix in the form it is used usually in the literature. We will follow here the representation given for open quantum systems in [10] and for PT symmetric systems used in [12]. In Sect. 2, the non-Hermitian Hamiltonian of an open quantum system is considered. The properties of its eigenvalues and eigenfunctions are sketched, above all in the neighborhood of one or more EPs. In the following section 3, two different non-Hermitian opera- tors that are used in the description of PT symmetric systems, are considered. The similarities and differ- ences to the Hamiltonian of an open quantum system are discussed on the basis of analytical studies (when possible) as well as by means of numerical results. The results are summarized in the last section. 2. Exceptional points in an open quantum system In an open quantum system, the discrete states de- scribed by a Hermitian Hamiltonian HB, are embed- ded into the continuum of scattering wavefunctions, which exists always and cannot be deleted. Due to this fact the discrete states turn into resonance states the lifetime of which is usually finite. The Hamiltonian H of the whole system consisting of the two subsys- tems is non-Hermitian. Its eigenvalues are complex and provide not only the energies of the states but also their lifetimes (being inverse proportional to the widths). The Hamiltonian of an open quantum system reads [10] H = HB + VBCG (+) C VCB, (1) where VBC and VCB stand for the interaction be- tween system and environment and G(+)C is the Green function in the environment. The so-called internal (first-order) interaction between two states i and j is involved in HB while their external (second-order) interaction via the common environment is described by the last term of (1). Generally, the coupling matrix elements of the exter- nal interaction consist of the principal value integral Re〈ΦBi |H|Φ B j 〉−E B i δij = 1 2π P ∫ �′c �c dE′ γ0icγ 0 jc E −E′ , (2) which is real, and the residuum Im〈ΦBi |H|Φ B j 〉 = − 1 2 γ0icγ 0 jc, (3) which is imaginary [10]. Here, the ΦBi and E B i are the eigenfunctions and (discrete) eigenvalues, respectively, of the Hermitian Hamiltonian HB which describes the states in the subspace of discrete states without any interaction of the states via the environment. The γ0ic ≡ √ 2π〈ΦBi |V |ξ E c 〉 are the (energy-dependent) coupling matrix elements between the discrete states i of the system and the environment of scattering wavefunctions ξEc . The γ0kc have to be calculated for every state i and for each channel c (for details see [10]). When i = j, (2) and (3) give the selfenergy of the state i. The coupling matrix elements (2) and (3) (by adding EBi δij in the first case) are often simulated by complex values ωij. In order to study the interaction of two states via one common environment it is convenient to start from two resonance states (instead of two discrete states). Let us consider, as an example, the symmetric 2 × 2 matrix H(2) = ( ε1 ≡ e1 + i2γ1 ω12 ω21 ε2 ≡ e2 + i2γ2 ) , (4) the diagonal elements of which are the two complex eigenvalues εi (i = 1, 2) of a non-Hermitian operator H0. This means that the ei and γi ≤ 0 denote the energies and widths, respectively, of the two states when ωij = 0 (the index c is ignored here for simplicity, c = 1). The ω12 = ω21 ≡ ω stand for the coupling of the two states via the common environment. The selfenergy of the states is assumed to be included into the εi. The two eigenvalues of H(2) are Ei,j ≡ Ei,j + i 2 Γi,j = ε1 + ε2 2 ±Z, Z ≡ 1 2 √ (ε1 −ε2)2 + 4ω2, (5) where Ei and Γi stand for the energy and width, respectively, of the eigenstate i. Resonance states with nonvanishing widths Γi repel each other in energy according to the value of Re(Z) while the widths bifurcate according to the value of Im(Z). The two states cross when Z = 0. This crossing point is an EP according to the definition of Kato [11]. Here, the two eigenvalues coalesce, E1 = E2. According to (5), two interacting discrete states (with γ1 = γ2 = 0) avoid always crossing since ω ≡ ω0 and ε1 − ε2 are real in this case and the condition Z = 0 cannot be fulfilled, (e1 −e2)2 + 4ω20 > 0. (6) In this case, the EP can be found only by analytical continuation into the continuum. This situation is known as avoided crossing of discrete states. It holds also for narrow resonance states if Z = 0 cannot be fulfilled due to the small widths of the two states. The physical meaning of this result has been very well known for many years. The avoided crossing of two discrete states at a certain critical parameter value [13] means that the two states are exchanged at this point, including their populations (population transfer). When ω = iω0 is imaginary, Z = 1 2 ( (e1 −e2)2 + 1 4 (γ1 −γ2)2 + i(e1 −e2)(γ1 −γ2) − 4ω20 )1/2 (7) 107 Hichem Eleuch, Ingrid Rotter Acta Polytechnica is complex. The condition Z = 0 can be fulfilled only when (e1−e2)2 + 14 (γ1−γ2) 2 = 4ω20 and (e1−e2)(γ1− γ2) = 0, i.e., when γ1 = γ2 (or when e1 = e2). In this case, it follows (e1 −e2)2 − 4ω20 = 0 → e1 −e2 = ±2ω0 (8) and two EPs appear. It holds further (e1 −e2)2 > 4ω20 → Z ∈<, (9) (e1 −e2)2 < 4ω20 → Z ∈= (10) independent of the parameter dependence of the ei. In the first case, the eigenvalues Ei = Ei + i2 Γi differ from the original values εi = ei + i/2γi by a contribution to the energies and in the second case by a contribution to the widths. The width bifurcation starts in the very neighborhood of one of the EPs and becomes maximum in the middle between the two EPs. This happens at the crossing point e1 = e2 where ∆Γ/2 ≡ |Γ1/2−Γ2/2| = 4ω0. A similar situation appears when γ1 ≈ γ2 as results of numerical calculations show. The physical meaning of this result is completely different from that discussed above for discrete and narrow resonance states. It means that different time scales appear in the system without any enhancement of the coupling strength to the continuum (for details see [14]). The cross section can be calculated by means of the S matrix σ(E) ∝ |1 − S(E)|2. A unitary repre- sentation of the S matrix in the case of two nearby resonance states coupled to one common continuum of scattering wavefunctions reads [10] S = (E −E1 − i2 Γ1)(E −E2 − i 2 Γ2) (E −E1 + i2 Γ1)(E −E2 + i 2 Γ2) . (11) In this expression, the influence of an EP on the cross section is contained in the eigenvalues Ei = Ei − i2 Γi of H(2). Reliable results can therefore be obtained also when an EP is approached and the S matrix has a double pole. Here, the line shape of the two overlapping resonances is described by S = 1 + 2i Γd E −Ed − i2 Γd − Γ2d (E −Ed − i2 Γd) 2 , (12) where E1 = E2 ≡ Ed and Γ1 = Γ2 ≡ Γd. It devi- ates from the Breit-Wigner line shape of an isolated resonance due to interferences between the two res- onances. The first term of (12) is linear (with the factor 2 in front) while the second one is quadratic. As a result, two peaks with asymmetric line shape appear in the cross section (for a numerical example see Fig. 9 in [15]). The eigenfunctions of the non-Hermitian H(2) are biorthogonal and can be normalized according to 〈Φ∗i |Φj〉 = δij, (13) although 〈Φ∗i |Φj〉 is a complex number (for details see sections 2.2 and 2.3 of [10]). The normalization (13) allows to describe the smooth transition from the regime with orthogonal eigenfunctions to that with biorthogonal eigenfunctions (see below). It follows 〈Φi|Φi〉 = Re〈Φi|Φi〉, Ai ≡〈Φi|Φi〉≥ 1 (14) and 〈Φi|Φj 6=i〉 = i Im〈Φi|Φj 6=i〉 = −〈Φj 6=i|Φi〉, |Bji | ≡ |〈Φi|Φj 6=i| ≥ 0. (15) At an EP Ai → ∞ and |B j i | → ∞. The Ei and Φi contain global features that are caused by many-body forces induced by the coupling ωik of the states i and k 6= i via the environment. They contain moreover the self-energy of the states i due to their coupling to the environment. At the EP, the eigenfunctions Φcri of H (2) of the two crossing states are linearly dependent on one another, Φcr1 →±iΦ cr 2 , Φ cr 2 →∓iΦ cr 1 (16) according to analytical as well as numerical and exper- imental studies, see the appendix of [14] and section 2.5 of [10]. This means that the wavefunction Φ1 of the state 1 jumps, at the EP, via the wavefunction Φ1 ± iΦ2 of a chiral state to ±iΦ2 [16]. The Schrödinger equation with the non-Hermitian operator H(2) is equivalent to a Schrödinger equation with H0 and source term [17] (H0 −εi)|Φi〉 = − ( 0 ωij ωji 0 ) |Φj〉≡ W|Φj〉. (17) Due to the source term, two states are coupled via the common environment of scattering wavefunctions into which the system is embedded, ωij = ωji ≡ ω. The Schrödinger equation (17) with source term can be rewritten in the following manner [17], (H0 −εi)|Φi〉 = ∑ k=1,2 〈Φk|W |Φi〉 ∑ m=1,2 〈Φk|Φm〉|Φm〉. (18) According to the biorthogonality relations (14) and (15) of the eigenfunctions of H(2), (18) is a nonlinear equation. The most important part of the nonlinear contributions is contained in (H0 −εn)|Φn〉 = 〈Φn|W |Φn〉|Φn|2|Φn〉. (19) The nonlinear source term vanishes far from an EP due to 〈Φk|Φk〉 → 1 and 〈Φk|Φl 6=k〉 = −〈Φl 6=k|Φk〉 → 0 according to (13) to (15). Thus, the Schrödinger equation with source term is linear far from an EP, as usually assumed. It is however nonlinear in the neighborhood of an EP. It is meaningful to represent the eigenfunctions Φi of H(2) in the set of basic wavefunctions Φ0i of H 0 Φi = N∑ j=1 bijΦ0j, bij = |bij|e iθij . (20) 108 vol. 54 no. 2/2014 Exceptional Points in Open and PT-Symmetric Systems Also the bij are normalized according to the biorthog- onality relations of the wavefunctions {Φi}. The angle θij can be determined from tan θij = Im(bij)/ Re(bij). From (13) and (16) follows: • When two levels are distant from one another, their eigenfunctions are (almost) orthogonal, 〈Φ∗k|Φk〉≈ 〈Φk|Φk〉 = Ak ≈ 1. • When two levels cross at the EP, their eigenfunc- tions are linearly dependent according to (16) and 〈Φk|Φk〉≡ Ak →∞. These two relations show that the phases of the two eigenfunctions relative to one another change when the crossing point is approached. This can be expressed quantitatively by defining the phase rigidity rk of the eigenfunction Φk, rk ≡ 〈Φ∗k|Φk〉 〈Φk|Φk〉 = A−1k . (21) It holds 1 ≥ rk ≥ 0. The non-rigidity rk of the phases of the eigenfunctions of H(2) follows also from the fact that 〈Φ∗k|Φk〉 is a complex number (unlike the norm 〈Φk|Φk〉, which is a real number) such that the normalization condition (13) can be fulfilled only by the additional postulation Im〈Φ∗k|Φk〉 = 0 (what generally corresponds to a rotation). When rk < 1, an analytical expression for the eigen- functions as a function of a certain control parameter can, generally, not be obtained. The non-rigidity rk < 1 of the phases of the eigenfunctions of H(2) in the neighborhood of EPs is the most important differ- ence between the non-Hermitian quantum physics and the Hermitian one. Mathematically, it causes nonlin- ear effects in quantum systems in a natural manner, as shown above. Physically, it allows the alignment of one of the states of the system to the common environment [10]. Results of numerical calculations are given, e.g., in [18]. The mixing coefficients bij (defined in (20)) of the wavefunctions of the two states due to their avoided crossing are simulated by assuming a Gaus- sian distribution for the coupling coefficients ωi 6=j = ωe−(ei−ej ) 2 (for real ω, the results of the simulation agree with the results [17] of exact calculations). In [18], results of different calculations are shown for illus- tration. Here, the coupling coefficients ω are assumed to be either real or complex or imaginary according to the different possibilities provided by (2) and (3). The main difference of the eigenvalue trajectories with real coupling coefficients ω to those with imagi- nary coupling coefficients ω is related to the relations (6) to (10) obtained analytically. For γ1 6= γ2 and real, complex or even imaginary ω, the results show one EP when the condition Z = 0 is fulfilled. This EP is isolated from other EPs, generally, when the level density is low. In the case of γ1 ≈ γ2 and imaginary ω however, two related EPs appear, see Fig. 1 right panel. Between these two EPs, the widths Γi bifurcate (Fig. 1d) while the energies Ei do not change (Fig. 1b). It is interesting to see that width bifurcation occurs between the two EPs, according to (8) and (10), with- out any enhancement of the coupling strength to the environment. Beyond the two EPs, the eigenvalues approach the original values. In a finite neighborhood of the point at which the two eigenvalue trajectories cross, the eigenfunc- tions are mixed and |bij|→∞ when approaching the EP (Fig. 1f). The phases of all components of the eigenfunctions jump at the EP either by −π/4 or by +π/4 [19]. This means that the phases of both eigen- functions jump in the same direction by the same amount. Thus, there is a phase jump of −π/2 (or +π/2) when one of the eigenfunctions passes into the other one at the EP. This result is in agreement with (16). It holds true for real as well as for imaginary ω. 3. Exceptional points in PT symmetric systems As has been shown in [4], the optical wave equation in PT symmetric optical lattices is formally equiva- lent to a quantum mechanical Schrödinger equation. Complex PT symmetric structures can be realized by involving symmetric index guiding and an antisym- metric gain/loss profile. The main difference between these optical systems and open quantum systems consists in the asymmetry of gain and loss in the first case while the states of an open quantum system can only decay (Im(ε1,2) < 0 and Im(E1,2) < 0 for all states). Thus, the modes involved in the non-Hermitian Hamiltonian in optics appear in complex conjugate pairs while this is not the case in an open quantum system. As a conse- quence, the Hamiltonian for PT symmetric structures in optical lattices may have real eigenvalues in a large parameter range. The 2 × 2 non-Hermitian Hamilto- nian may be written as [4, 12] HPT = ( e− iγ2 w w∗ e + iγ2 ) , (22) where e stands for the energy of the two modes, ±γ describes gain and loss, respectively, and the coupling coefficients w stand for the coupling of the two modes via the lattice. When the PT symmetric optical lat- tices are studied with vanishing gain, the Hamiltonian reads H′PT = ( e− iγ2 w w∗ e ) . (23) In realistic systems, w in (22) and (23) is mostly real (or almost real) [20]. The eigenvalues of the Hamiltonian (22) differ from (5), EPT± = e± 1 2 √ 4|w|2 −γ2 ≡ e±ZPT . (24) A similar expression is derived in [5]. Since e and γ are real, the EPT± are real when 4|w|2 > γ2. Under this 109 Hichem Eleuch, Ingrid Rotter Acta Polytechnica condition, the two levels repel each other in energy, which is characteristic of discrete interacting states. The level repulsion decreases with increasing γ (when the interaction w is fixed). When 4|w|2 = γ2 the two states cross. Here, EPT± = e and γ = ± √ 4|w|2. With further increasing γ and 4|w|2 < γ2 (w fixed for illus- tration), width bifurcation (PT symmetry breaking) occurs and EPT± = e± i 2 √ γ2 − 4|w|2. These relations are in accordance with (8) to (10) for open quantum systems. Two EPs exist according to 4|w|2 = (±γ)2. (25) Further γ2 < 4|w|2 → ZPT ∈<, (26) γ2 > 4|w|2 → ZPT ∈= (27) independent of the parameter dependence γ(a) and of the ratio Re(w)/ Im(w). In the case of the Hamiltonian (23), the eigenvalues read E ′PT ± = e−i γ 4 ± 1 2 √ 4|w|2 − γ2 4 ≡ e−i γ 4 ±Z′PT . (28) We have level repulsion as long as 4|w|2 > γ 2 4 . While level repulsion decreases with increasing γ, loss in- creases with increasing γ. At the crossing point, E ′PT ± = e − i γ 4 . With further increasing γ and 4|w|2 � γ 2 4 E ′PT ± → e− i γ 4 ± i γ 4 = { e e− iγ2 . (29) The two modes (29) behave differently. While the loss in one of them is large, it is almost zero in the other one. Thus, only one of the modes effectively survives. Equation (29) corresponds to high transparency at large γ. Further, two EPs exist according to 4|w|2 = (±γ/2)2 (30) and γ2/4 < 4|w|2 → Z ′ PT ∈<, (31) γ2/4 > 4|w|2 → Z ′ PT ∈=. (32) By analogy with to (25) up to (27), these relations are independent of the parameter dependence of γ and of the ratio Re(w)/ Im(w). Thus, the difference between the eigenvalues Ei of H(2) of an open quantum system and the eigenval- ues of the Hamiltonian of a PT symmetric system consists, above all, in the fact that the Ei depend on the ratio Re(ω)/ Im(ω) while the EPT± and E ′PT ± are independent of Re(w)/ Im(w). There exist however similarities between the two cases. It is interesting to compare the eigenvalues Ei of H(2) obtained for imaginary non-diagonal matrix elements ω, with the eigenvalues of (22) (or (23)) obtained for real w. In both cases, there are two EPs, see Fig. 1. In the first case (right panel), the energies Ei of both states are equal and the widths Γi bifurcate between the two EPs. This situation is characteristic of an open quantum system at high level density with complex (almost imaginary) ω, see Eqs. (8) to (10). In the second case (left panel) however the difference |E1 − E2| of the energies first increases (level repulsion) and then decreases again while the widths Γi of both states vanish in the parameter range between the two EPs in accordance with the analytical results (25) to (27). Between the two EPs, level repulsion causes the two levels to be distant from one another and w is expected to be (almost) real. This result agrees qualitatively with (2) and (3). Similar results are obtained for the eigenvalues of (23). The only difference from those of (22) is that the Γi do not vanish but decrease between the two EPs with increasing a in this case. According to Figs. 1a–d, the role of energy and width is formally exchanged when the eigenvalues of the Hamiltonian (4) are compared with those of (22) (or (23)). In any case, the eigenvalues are influenced strongly by the EPs. Also the eigenfunctions of the Hamiltonian (4) of an open quantum system (with imaginary ω) and those of the Hamiltonians (22) and (23) of a PT sym- metric system (with real w) show similar features. The eigenfunctions ΦPTi of HPT (and Φ ′PT i of H ′ PT ) are biorthogonal with all the consequences discussed in Sect. 2. In contrast to the eigenvalues, they are dependent on the ratio Re(ω)/ Im(ω). The eigenfunctions can be represented in a set of basic wavefunctions in full analogy to the represen- tation of the eigenfunctions Φi of H(2) in (20). They contain valuable information on the mixing of the wavefunctions under the influence of the non-diagonal coupling matrix elements w and w∗ in (22) and (23), respectively, and its relation to EPs. Due to the level repulsion occurring between the two EPs, the coupling coefficients w can be considered to be (almost) real in realistic cases. The phases of the eigenmodes of the non-Hermitian Hamiltonians (22) and (23) are not rigid, generally, in approaching an EP, and spectro- scopic redistribution processes occur in the system under the influence of the environment (lattice). As in the case of open quantum systems, the phase rigidity rk can be defined according to (21). It varies between 1 and 0 and is a quantitative measure for the skewness of the modes when the crossing point is approached. In Figs. 1ef, the eigenfunctions of the Hamiltonian (22) (calculated with real w) are compared to those of the Hamiltonian (4) (calculated with imaginary ω). They show the same characteristic features. As can be seen from Fig. 1e, PT symmetry breaking is accompanied by a mixing of the eigenfunctions in a finite neighborhood of the EPs in PT symmetric systems. This result is in complete analogy to the results shown in Fig. 1f for open quantum systems 110 vol. 54 no. 2/2014 Exceptional Points in Open and PT-Symmetric Systems Figure 1. Energies Ei, widths Γi/2 and wavefunctions |bij| of N = 2 states coupled to K = 1 channel as function of a of a PT symmetric system with Hamiltonian (22) (left panel) and of an open quantum system with Hamiltonian (4) (right panel). Parameters left panel: e = 0.5, γ1 = −γ2 = 0.05a, w = 0.05; right panel: e1 = 1 − 0.5a, e2 = a, γ1/2 = γ2/2 = 0.5, ω = 0.05i. The dashed lines in (a,b) show ei(a). where a hint of width bifurcation can be seen in the mixing of the eigenfunctions around these points. Also the phases of the eigenfunctions jump in both cases by π/4 at the EPs (not shown here). In the parameter region between the two EPs, the eigenfunctions are completely mixed (1:1) in both cases while they are unmixed far beyond the EPs, see Figs. 1ef. 4. Discussion of the results On the basis of 2 × 2 models, we have compared the influence of an EP on the dynamics of an open quantum system with its influence on PT symme- try breaking in a PT symmetric system. In the first case the coupling of the two states via the environ- ment is symmetric (ω12 = ω21 ≡ ω). In the second case however, the formal equivalence of the optical wave equation in PT symmetric optical lattices with a quantum mechanical Schrödinger equation causes the two nondiagonal matrix elements to be complex conjugate (w21 = w∗12). The eigenvalues depend in the first case on the ratio Re(ω)/ Im(ω) while they are independent of Re(w)/ Im(w) in the second case. The eigenfunctions are sensitive to Re(ω)/ Im(ω) and Re(w)/ Im(w), respectively, in both cases. The EPs cause nonlinear effects in their neighbor- hood which determine the evolution of open as well as of PT symmetric systems. Most important for the dynamics of an open quantum system is the regime at high level density where the coupling coefficients are (almost) imaginary. Here, two EPs appear when the decay widths γi of both states are (almost) the same. Approaching the EPs, width bifurcation starts and ends, respectively, while beyond the EPs the widths of both states are equal (or similar) to one another. The energies of the two states show an opposite be- havior: it is Ei = E2 (or Ei ≈ E2) in the parameter range between the two EPs while the states repel each other in energy beyond the EPs. The width bifurca- tion related to the two EPs becomes relevant for the dynamics of an open quantum system at high level density. Here, short-lived and long-lived states are formed which are related to different time scales of the system (for details see [14]). Two EPs appear also in a PT symmetric system, and PT symmetry breaking is directly related to them. From a mathematical point of view however, energy and time are exchanged in comparison with the cor- responding values in an open quantum system. This means that the widths of both states are equal and 111 Hichem Eleuch, Ingrid Rotter Acta Polytechnica vanish in the case of the Hamiltonian (22) with gain and loss in the whole parameter range between the two EPs. In this parameter range, the eigenvalues are real and, furthermore, level repulsion prohibits a small energy distance between the two levels. Therefore the non-diagonal coupling matrix elements w are (almost) real, Re(w) � Im(w). The eigenfunctions of the different 2 × 2 models considered in the present paper show very clearly that the spectroscopic redistribution inside the system is indeed caused by the EPs. However, it shows up in all cases in a finite neighborhood around them. Here the rigidity of the phases of the two eigenfunctions relative to one another is reduced (ri < 1) and an alignment of one of the states to the environment is possible. In the parameter range between the two EPs, the wavefunctions are completely mixed (1:1) as can be seen from the numerical results shown in Fig. 1. Summing up the discussion we state the following. The results obtained by studying PT symmetric opti- cal lattices as well as those received from an investiga- tion of open quantum systems show the characteristic features of non-Hermitian quantum physics. They prove environmentally induced effects that cannot be described convincingly in conventional Hermitian quantum physics. Due to the reduced phase rigidity around an EP, the system is able to align (at least partly) with the environment. This can be seen from PT symmetry breaking occurring in one of the con- sidered systems as well as from the dynamical phase transition taking place at high level density in the other system. References [1] Bender, C.M., and Boettcher, S.: Real Spectra in Non-Hermitian Hamiltonians Having PT Symmetry, Phys. Rev. Lett. 80 (1998) 5243-5246 [2] Bender, C.M.: Making sense of non-Hermitian Hamiltonians, Rep. Progr. Phys. 70 (2007) 947-1018 [3] Special Issue Quantum physics with non-Hermitian operators, J. Phys. A 45, Number 44 (November 2012) [4] Ruschhaupt, A., Delgado, F. and Muga, J.G.: Physical realization of PT symmetric potential scattering in a planar slab waveguide, J. Phys. A 38 (2005) L171-L176; El-Ganainy, R., Makris, K.G., Christodoulides, D.N. and Musslimani, Z.H.: Theory of coupled optical PT- symmetric structures, Optics Lett. 32 (2007) 2632-2634; Makris, K.G., El-Ganainy, R., Christodoulides, D.N. and Musslimani, Z.H.: Beam Dynamics in PT Symmetric Optical Lattices, Phys. Rev. Lett. 100 (2008) 103904 (4pp); Musslimani, Z.H., Makris, K.G., El-Ganainy, R. and Christodoulides, D.N.: Optical Solitons in PT Periodic Potentials, Phys. Rev. Lett. 100 (2008) 030402 (4pp) [5] Guo, A., Salamo, G.J., Duchesne, D., Morandotti, R., Volatier-Ravat, M., Aimez, V., Siviloglou, G.A., and Christodoulides, D.N.: Observation of PT-symmetry breaking in complex optical potentials, Phys. Rev. Lett. 103 (2009) 093902 (4pp) [6] Rüter, C.E., Makris, G., El-Ganainy, R., Christodoulides, D.N., Segev, M., and Kip, D.: Observation of parity-time symmetry in optics, Nature Physics 6 (2010) 192-195 [7] Kottos, T.: Broken symmetry makes light work, Nature Physics 6 (2010) 166-167 [8] Feshbach, H.: Unified theory of nuclear reactions, Ann. Phys. (N.Y.) 5 (1958) 357-390; Feshbach, H.: A unified theory of nuclear reactions. II, Ann. Phys. (N.Y.) 19 (1962) 287-313 [9] I. Rotter: A continuum shell model for the open quantum mechanical nuclear system, Rep. Prog. Phys. 54 (1991) 635-682 [10] Rotter, I.: A non-Hermitian Hamilton operator and the physics of open quantum systems, J. Phys. A 42 (2009) 153001 (51pp) doi: 10.1088/1751-8113/42/15/153001 [11] Kato, T.: Perturbation Theory for Linear Operators Springer Berlin, 1966 [12] Rotter, I.: Environmentally induced effects and dynamical phase transitions in quantum systems, J. Opt. 12 (2010) 065701 (9pp) doi: 10.1088/2040-8978/12/6/065701 [13] Landau, L., Physics Soviet Union 2, 46 (1932); Zener, C.: Non-Adiabatic Crossing of Energy Levels, Proc. Royal Soc. London, Series A 137 (1932) 696-702 [14] Rotter, I.: Dynamical stabilization and time in open quantum systems, Contribution to the Special Issue Quantum Physics with Non-Hermitian Operators: Theory and Experiment, Fortschritte der Physik - Progress of Physics 61 No. 2-3 (2013) 178-193 doi: 10.1002/prop.201200054 [15] Müller, M., Dittes, F.M., Iskra, W. and Rotter, I.: Level repulsion in the complex plane, Phys. Rev. E 52 (1995) 5961-5973 [16] In studies by other researchers, the factor i in (16) does not appear. This difference is discussed in detail and compared with experimental data in the Appendix of [14] and in section 2.5 of [10]. In the present paper, the definitions εi = ei + i2 γi and Ei = Ei + i 2 Γi (with γi ≤ 0 and Γi ≤ 0 for decaying states) are used [17] Rotter, I.: Dynamics of quantum systems, Phys. Rev. E 64 (2001) 036213 (12pp) [18] Eleuch, H. and Rotter, I.: Width bifurcation and dynamical phase transitions in open quantum systems, Phys. Rev. E 87 (2013) 052136 (15pp) doi: 10.1103/PhysRevE.87.052136; Eleuch, H. and Rotter, I.: Avoided level crossings in open quantum systems, Contribution to the Special Issue Quantum Physics with Non-Hermitian Operators: Theory and Experiment, Fortschritte der Physik - Progress of Physics 61 No. 2-3 (2013) 194-204 doi: 10.1002/prop.201200062 [19] Eleuch, H. and Rotter, I., Eur. Phys. J. D 68 (2014) 74 (16pp) doi: 10.1140/epjd/e2014-40780-8 [20] Kottos, T., private communication 112 http://dx.doi.org/10.1088/1751-8113/42/15/153001 http://dx.doi.org/10.1088/2040-8978/12/6/065701 http://dx.doi.org/10.1002/prop.201200054 http://dx.doi.org/10.1103/PhysRevE.87.052136 http://dx.doi.org/10.1002/prop.201200062 http://dx.doi.org/10.1140/epjd/e2014-40780-8 Acta Polytechnica 54(2):106–112, 2014 1 Introduction 2 Exceptional points in an open quantum system 3 Exceptional points in PT symmetric systems 4 Discussion of the results References