PTsymmetry breaking and maximal chirality in a nonuniform PTsymmetric ring
Abstract
We study the properties of an site tightbinding ring with parity and timereversal () symmetric, Hermitian, sitedependent tunneling and a pair of nonHermitian, symmetric, loss and gain impurities . The properties of such lattices with open boundary conditions have been intensely explored over the past two years. We numerically investigate the symmetric phase in a ring with a positiondependent tunneling function that, in an open lattice, leads to a strengthened symmetric phase, and study the evolution of the symmetric phase from the open chain to a ring. We show that, generally, periodic boundary conditions weaken the symmetric phase, although for experimentally relevant lattice sizes , it remains easily accessible. We show that the chirality, quantified by the (magnitude of the) average transverse momentum of a wave packet, shows a maximum at the symmetric threshold. Our results show that although the wavepacket intensity increases monotonically across the breaking threshold, the average momentum decays monotonically on both sides of the threshold.
1 Introduction
One of the founding tenets of quantum theory is that the Hamiltonian of a ”closed” system is Hermitian with respect to the standard inner product; this property of the Hamiltonian implies that the energy spectrum of the system is purely real, and that its eigenvectors form a complete orthonormal set [1, 2, 3]. The quantum theory, based on this and other axioms, has been phenomenally successful in explaining experimental results, and predicting various nontrivial and counterintuitive phenomena. Traditionally, nonHermitian Hamiltonians, with eigenvalues that have an imaginary part, have been used to model open, dissipative quantum systems such as a resistor [4]. Since the seminal work of Bender and coworkers in the late nineties, over the past decade, a new class of nonHermitian Hamiltonians has been extensively explored [5, 6, 7]. Although these Hamiltonians are not Hermitian under the standard inner product, they are invariant under the combined parity () and timereversal () operation, and thus are called symmetric Hamiltonians. The eigenvalues of such a Hamiltonian are real over a range of parameters, although the corresponding eigenvectors are not orthogonal under the standard inner product. This parameter space where all eigenvalues are real is called the symmetric phase, and the emergence of the first pair of complex eigenvalues that occurs when the Hamiltonian leaves this parameter space is called symmetry breaking. Since timereversal is an antilinear operator, it follows that when the symmetry is broken, the eigenvectors are no longer simultaneous eigenfunctions of the combined operation.
Over the past three years, symmetric Hamiltonians for tightbinding lattice models have been extensively investigated [8, 9, 10, 11]. They have included one dimensional chains with nonHermitian, symmetric tunneling [12], nonHermitian symmetric pair of onsite impurity potentials [13], twodimensional lattices with symmetric potential [14], onedimensional chain of dimers with symmetric impurities on each dimer [15], and symmetric spin models [16] where the effects of boundary conditions are claimed to be negligible. These theoretical investigations have been accompanied by the experimental exploration of consequences of symmetry breaking in coupled optical waveguides [17, 18, 19] and electrical circuits [20]. These ongoing studies have hinted at the deep and rich set of phenomena that follow from nonHermitian Hamiltonians with balanced loss and gain. Although most of these studies have focused on open chains, some recent works have shown that the effect of boundary conditions on the symmetric phase is nontrivial [21, 22].
The current and emerging promising candidates for symmetric systems  coupled optical waveguides and coupled electrical circuits  are presently focused on the realization of symmetric open chains. For electrical circuits, since the spatial circuit geometry can be easily manipulated, introducing inductive coupling between the “end sites of an open chain” is relatively straightforward [20]. On the other front, two dimensional (2D) lattices of coupled optical waveguides have been extensively experimentally explored [23, 24] in the context of Anderson localization. Thus, creating a onedimensional system with periodic boundary conditions  a 2D lattice with only boundary waveguides and no interior waveguides  is feasible with current fabrication technology. Therefore, realization and experimental investigation of one dimensional symmetric rings is likely to occur in near future.
In this paper, we numerically investigate the symmetric phase diagram in a one dimensional system with periodic boundary condition (a symmetric ring). We consider an site lattice with a positiondependent, paritysymmetric tunneling function , and a pair of balanced gain and loss impurities . In the Hermitian limit, , such nonuniform open chain shows tunable energy spectra, wave packet evolution, and HanburyBrownTwiss correlations [25]. Our primary results are a follows: i) For , we find that the symmetric phase is generally weakened in comparison with its open chain counterpart. However, for experimentally relevant system sizes , the critical impurity strength is of the same order of magnitude as its openchain counterpart. ii) For , for a ring is substantially enhanced from its open chain value for all impurity locations. iii) The average transverse momentum of a wave packet shows a maximum at the symmetry breaking threshold, and decays monotonically on the two sides of it.
The plan of the paper is as follows. In Sec. 2, we introduce the tightbinding model and discuss the differences between properties of an open chain and a ring that become relevant for small . In Sec. 3, we present the results for the symmetric phase diagram as a function of the scale factor that determines the tunneling between end points of an open chain. In Sec. 4, we discuss the time evolution of a wave packet across the symmetric threshold. We show that the chirality, which encodes the average clockwise or counterclockwise motion of the wave packet around the ring, is quantified by the average momentum and shows a maximum at the threshold. The paper is concluded with a brief discussion and open questions in Sec. 5. Throughout this paper, we consider coupled optical waveguides with complex refractive index and gain as the prototype realization of a symmetric ring [17, 18, 19].
2 TightBinding Model
We start with a onedimensional, tightbinding, site lattice with two nonHermitian impurities located at mirrorsymmetric sites , where and . The Hermitian, positiondependent tunneling Hamiltonian for the open chain is
(1) 
where is the tunneling between sites and , and () represents the creation (annihilation) operator for a singleparticle state localized on site . The difference equation for the coefficients of an eigenfunction of the Hamiltonian Eq.(1) with eigenenergy is given by
(2) 
where open boundary conditions imply . We remind the reader that since the tunneling amplitudes is site dependent, traditional methods for solving Eq.(2) such as the Bethe ansatz are not applicable. Apart from a few exceptions [10], the eigenvalues and eigenvectors of Eq.(1) must be obtained numerically. The impurity potential is given by
(3) 
The action of the parity operator is given by , where and the antilinear timereversal operator implies . From these definitions, it follows that the nonHermitian potential is odd under parity and timereversal individually, and but obeys . The symmetry breaking in such a nonuniform open chain has been extensively studied, albeit only numerically due to the absence of analytical methods that are applicable to a system with sitedependent tunneling [26]. The Hamiltonian for a symmetric ring is given by
(4) 
where we choose the tunneling between the end points of the open chain, sites and , as and the scalefactor allows us to continuously extrapolate from an open chain to a symmetric ring. Before discussing the results for a ring, , we briefly recall the results for an open chain and establish the terminology [26]. In the symmetric phase, the eigenvalue spectrum of an site open chain is nondegenerate and symmetric about zero, and the eigenfunctions for energies are related to each other. The bandwidth of the energy spectrum scales as for and for because the bandwidth is determined by the largest tunneling amplitude. We use quarterbandwidth, as the energy scale and note that when , the threshold impurity strength is given by .
Note that although the distinction between an open chain and a ring is expected to vanish [16] in the limit , for small the differences between the two can be substantial. As an extreme case, let us consider a uniform 3site lattice with nearestneighbor tunneling . The nondegenerate, particlehole symmetric spectrum of such an open chain is given by . On the other hand, the spectrum of a threesite ring is given by and is, in general, asymmetric about zero and degenerate.
The lefthand panel in Fig. 1 shows typical spectra for Hamiltonian 1 with and . The energy is normalized by its maximum value and since , the spectrum is virtually identical for an open chain and a ring. When (black line) we obtain the expected cosine spectrum, for (green line) and (red line), the spectrum is quasilinear, and when (blue line), the spectrum bandedges represent eigenstates localized at the end of the chain that are generically present when [25]. The righthand panel in Fig. 1 shows the difference between energies of a ring and an open chain for a lattice with sites and . We find that for the difference is greatest near the center of the band, whereas for the spectral difference is greatest at the band edges. This is expected since the localized edgestates that exist for are most influenced by the introduction of periodic boundary condition. The spectral differences vanish with increasing but remain pertinent for experimentally relevant lattice sizes [27, 28, 29].
3 symmetric Phase Diagram
We now discuss the numerically obtained symmetric phase diagram for the Hamiltonian as a function of increasing loss and gain impurity strength , fractional location of the gain impurity, , and the dimensionless scale factor that extrapolates between an open chain () and a ring (). The symmetric phase is called robust if the critical impurity strength measured in units of the energy scale, is nonzero as ; it is called fragile if as . In an open chain corresponds to farthest impurities whereas the in a ring, farthest impurities correspond to . We remind the reader that in an open chain, when the critical value of impurity strength decreases as the distance between the impurities increases whereas for , the symmetric phase is vanishingly small for almost all values of impurity location [26].
Figure 2 shows the typical evolution of the symmetric phase diagram in the plane as a function of the scale parameter . These results are for a lattice with sites, and (lefthand panel) and (righthand panel). For (black circles), the symmetric phase is robust when the loss and gain impurities are closest to each other, . Then the critical impurity strength is given by when is even and when is odd [26]. For , we see that the critical impurity strength is, in general, suppressed relative to its value for an open chain for all values of impurity positions . Thus, the symmetric phase in a ring with nonuniform tunneling is weaker than its counterpart in an open chain. However, results in Fig. 2 also show that the critical impurity strength is still an appreciable fraction of its value in an open chain. The lefthand panel shows that for a lattice with , the critical impurity strength is strongly suppressed except for a few specific values of ; this feature is robust irrespective of and is related to the exactly linear spectrum of an open chain [10, 12]. The righthand panel shows that for lattice, with a nonlinear spectrum, the critical impurity strength reaches a plateau for most impurity locations . As is increased from zero (red squares, green diamonds) to one (blue triangles) the symmetric Hamiltonian, Eq.(4), is perturbed only slightly. The tunneling introduced between sites and , , is the smallest among all links for . However, it has a dramatic effect on even when the loss and gain impurities are away from this tunneling link, .
When , the tunneling introduced between sites and is the largest among all links and is proportional to the bandwidth . Therefore, we can expect that the fragile symmetric phase in an open chain will be strengthened in a ring. Fig. 3 shows the symmetric phase diagram for an , lattice. When (black circles), we obtain the extremely fragile phase diagram of an open chain. As the tunneling is increased from (red squares), (green diamonds) to (blue triangles), we see that the symmetric phase in a ring is substantially strengthened for all impurity positions. Since symmetric rings () with and can be mapped onto each other with an appropriate redefinition of the impurity location , the critical impurity strength decreases with when (Fig. 3) whereas it increases with for (Fig. 2).
We remind the reader that in a ring with constant tunneling, accessible to analytical treatment, the critical impurity threshold is zero [22]. These numerical results show that in a ring with nonuniform tunneling there is a large region, below the bluetrianglescurve in Figs. 2 and 3, where the symmetry is exact.
4 Chirality across the symmetric Phase Boundary
The timeevolution of a wave packet that is initially localized to a single site has been traditionally used to probe the degrees and signatures of symmetry breaking in coupled optical waveguides [18, 19]. In these systems, the wave function denotes the singletransversemode electric field in each waveguide [27], the symmetric impurities which absorb or emit the corresponding electromagnetic radiation are engineered, and the time evolution of an initially normalized wave packet is given by where is the scaled Planck’s constant. The timeevolution operator is not unitary irrespective of whether the symmetry is exact or broken; therefore, in general, the net intensity is not one, where denotes the site and timedependent intensity. This violation of unitarity is determined by the parity of number of lattice sites [26], and the localized or extended nature of eigenstates [30]. In all cases, however, the net intensity increases monotonically across and exponentially past the symmetric threshold [31].
The typical timeevolution of an initially localized state shows that apart from spreading across different sites, the wave packet undergoes a preferential clockwise or anticlockwise motion around the ring. We quantify this tendency, chirality, by a dimensionless, Hermitian, discretemomentum operator on a ring [22],
(5) 
where is a timeevolved wave function, and the normalization factor in the denominator is required due to its nonunitary time evolution. Note that due to the CauchySchwartz inequality, the dimensionless momentum satisfies .
When the Hamiltonian is Hermitian and the initial state is localized to a single site, the momentum symmetrically oscillates about zero. As the impurity strength approaches the threshold value, , remains constant over long time intervals where the time interval is measured in units of [26]. We use the timeaveraged momentum to denote this steady state value, , and choose the time interval such that the steadystate value is independent of it.
Figure 4 shows the evolution of the dimensionless momentum across the symmetric phase boundary for an ring with tunneling functions (lefthand panel) and (righthand panel). The initial positions of the wave packet are at (blue circles dotted line) and (red squares solid line). Note that when , the threshold impurity strengths for and differ by an order of magnitude (Fig 2). In all cases, the average momentum is zero when . Its linear increase with at small and the different slopes for different initial wave packet locations are both expected from a firstorder perturbation theory. Remarkably, in all cases, reaches the same maximum possible value (of unity) at the symmetric threshold and decreases monotonically beyond it even though the net intensity increases exponentially past the threshold.
At this point, it is worthwhile to recall the corresponding results for a symmetric ring with two (constant) tunnelings between the loss and gain impurities. In that case, the threshold strength is independent of the impurity location , the sign of the momentum is determined by the path with the higher tunneling amplitude, and it reaches a universal maximum value of one at the threshold [26].
Figure 5 shows that for a ring with nonuniform, positiondependent tunneling profile , the behavior of the steadystate momentum is not as straightforward. Both panels present results for an ring with and two initial wave packet locations, (blue circles dotted line) and (red squares solid line). The lefthand panel has the gain impurity at location () and it shows that the sign of the steadystate momentum is now negative. The righthand panel has gain impurity at location () and it shows that the maximum value attained by the steadystate momentum is not unity. These results imply that, contrary to the twotunneling ring, sign of the momentum and the maximum value it attains at the threshold are both dependent upon the impurity location. We emphasize, however, that these results are independent of the initial wave packet location, and the qualitative behavior of across the symmetric threshold is identical in all cases.
5 Discussion
In this paper, we have numerically explored the symmetric phase diagram for a finite, nonuniform symmetric lattice as a function of its boundary conditions. We have shown that for experimentally relevant lattice sizes, the differences between properties of an open chain and a ring are nontrivial, particularly for tunneling profiles with . Generically, we found that for , the symmetric phase in a ring is weaker than its counterpart in an open chain. In contrast, when the symmetric phase in a ring is substantially strengthened. Since for the tunneling perturbation that is required to change an open chain into a ring is comparable to the bandwidth of the open chain, the strengthening of the symmetric phase is reasonable.
We have shown that the symmetry breaking is accompanied by a qualitatively universal behavior of dimensionless, average momentum : the momentum is zero when , increases linearly with , and its magnitude reaches maximum at the breaking threshold, accompanied by monotonic decay on both sides of the threshold. We have also found that, in contrast with the twotunneling model [22], here, the sign of the momentum and its maximum value at the threshold are dependent upon the impurity locations, but not on the initial wave packet location.
These numerical results raise several questions. For an open chain, the location of the pair of eigenvalues that become degenerate and then complex is uniquely determined by the impurity location irrespective of the value of [26]; no such claim seems possible for the symmetry breaking in a ring. The universal presence of the peak in at the threshold suggests that it may be driven by the exceptional point at the threshold where two eigenvalues become degenerate and the corresponding eigenvectors become parallel [32]. However, lacking analytical methods for the nonuniform tunneling ring, a systematic numerical investigation of the dependence of the sign and the maximum value of the chirality remains an open problem.
6 Acknowledgment
This work was supported by NSF grant DMR1054020.
References
References
 Dirac P A M 1996 The Principles of Quantum Mechanics (New York, NY: Oxford University Press)
 Sakurai J J 1995 Modern Quantum Mechanics (Reading, MA: AddisonWesley)
 Dennery P and Krzywicki A 1995 Mathematics for Physicists (Mineola, NY: Dover Publications)
 Bruus H and Flensberg K 2010 Manybody Quantum Theory in Condensed Matter Physics (New York, NY: Oxford University Press)

Bender C M and Boettcher S 1998 Phys. Rev. Lett. 80 5243
Bender C M, Brody D C, and Jones H F 2002 Phys. Rev. Lett. 89 270401  For a review, see Bender C M 2007 Rep. Prog. Phys. 70 947
 For a review, see Mostafazadeh A 2010 Int. J. Geom. Meth. Mod. Phys. 7 1191 (arXiv:08105643)

Znojil M 2010 Phys. Rev. A 82 052113
Znojil M 2007 Phys. Lett. A 40 13131  Bendix O, Fleischmann R, Kottos T, and Shapiro B 2009 Phys. Rev. Lett. 103 030402
 Longhi S 2010 Phys. Rev. B 82 041106(R)
 Jin L and Song Z 2009 Phys. Rev. A 80 052107
 Joglekar Y N and Saxena A 2011 Phys. Rev. A 83 050101(R)
 Joglekar Y N, Scott D D, Babbey M, and Saxena A 2010 Phys. Rev. A 82 030103(R)
 Szameit A, Rechtsman M C, BahatTreidel O, and Segev M 2011 Phys. Rev. A 84 021806(R)
 Bendix O, Fleischmann R, Kottos T, and Shapiro B 2010 J. Phys. A 43 265305
 CastroAlvared O A and Fring A 2009 J. Phys. A 42 465211
 Guo A et al. 2009 Phys. Rev. Lett. 103 093902
 Rüter C E, Makris K G, ElGanainy R, Christodoulides D N, Segev M, and Kip D 2010 Nat. Phys. 6 192
 Feng L et al. 2011 Science 333 729
 Schindler J, Li A, Zheng M C, Ellis F M, and Kottos T 2011 Phys. Rev. A 84 040101(R)
 Znojil M 2011 Phys. Lett. A 375 3435
 Joglekar Y N and Scott D D 2012 Phys. Rev. A 85 062105
 Schwartz T, Bartal G, Fishman S, and Segev M 2007 Nature 446 52
 Szameit A and Nolte S 2010 J. Phys. B 43 163001
 Joglekar Y N, Thompson C, and Vemuri G 2011 Phys. Rev. A 83 063817
 Scott D D and Joglekar Y N 2011 Phys. Rev. A 83 050102(R)
 Christodoulides D N, Lederer F, Silberberg Y 2003 Nature 424 817
 Lahini Y, Avidan A, Pozzi F, Sorel M, Morandotti R, and Silberberg Y 2008 Phys. Rev. Lett. 100 013906
 PerezLeija A, MoyaCessa H, Szameit A, and Christodoulides D N 2010 Opt. Lett. 35 2409
 Vemuri H, Bhimadipati T, Vavilala V, and Joglekar Y N 2011 Phys. Rev. A 84 043826
 Zheng M C, Christodoulides D N, Fleischmann R, and Kottos T 2010 Phys. Rev. A 82 010103(R)
 Heiss W D and Harney H L 2001 Eur. J. Phys. D 17 149