• Photonics Research
  • Vol. 7, Issue 8, 862 (2019)
Mario A. Quiroz-Juárez1, Armando Perez-Leija2、3, Konrad Tschernig2、3, Blas M. Rodríguez-Lara4、5, Omar S. Magaña-Loaiza6, Kurt Busch2、3, Yogesh N. Joglekar7、8、*, and Roberto de J. León-Montiel1、9、*
Author Affiliations
  • 1Instituto de Ciencias Nucleares, Universidad Nacional Autónoma de México, Apartado Postal 70-543, 04510 Cd. Mx., Mexico
  • 2Max-Born-Institut, Max-Born-Straße 2A, 12489 Berlin, Germany
  • 3Humboldt-Universität zu Berlin, Institut für Physik, AG Theoretische Optik & Photonik, D-12489 Berlin, Germany
  • 4Tecnológico de Monterrey, Escuela de Ingeniería y Ciencias, Ave. Eugenio Garza Sada 2501, 64849 Monterrey, N.L., Mexico
  • 5Instituto Nacional de Astrofísica, Óptica y Electrónica, Calle Luis Enrique Erro No. 1, Sta. Ma. Tonantzintla, Pue. CP 72840, Mexico
  • 6Department of Physics and Astronomy, Louisiana State University, Baton Rouge, Louisiana 70803, USA
  • 7Department of Physics, Indiana University Purdue University Indianapolis (IUPUI), Indianapolis, Indiana 46202, USA
  • 8e-mail: yojoglek@iupui.edu
  • 9e-mail: roberto.leon@nucleares.unam.mx
  • show less
    DOI: 10.1364/PRJ.7.000862 Cite this Article Set citation alerts
    Mario A. Quiroz-Juárez, Armando Perez-Leija, Konrad Tschernig, Blas M. Rodríguez-Lara, Omar S. Magaña-Loaiza, Kurt Busch, Yogesh N. Joglekar, Roberto de J. León-Montiel, "Exceptional points of any order in a single, lossy waveguide beam splitter by photon-number-resolved detection," Photonics Res. 7, 862 (2019) Copy Citation Text show less

    Abstract

    Exceptional points (EPs) are degeneracies of non-Hermitian operators where, in addition to the eigenvalues, the corresponding eigenmodes become degenerate. Classical and quantum photonic systems with EPs have attracted tremendous attention due to their unusual properties, topological features, and an enhanced sensitivity that depends on the order of the EP, i.e., the number of degenerate eigenmodes. Yet, experimentally engineering higher-order EPs in classical or quantum domains remain an open challenge due to the stringent symmetry constraints that are required for the coalescence of multiple eigenmodes. Here, we analytically show that the number-resolved dynamics of a single, lossy waveguide beam splitter, excited by N indistinguishable photons and post-selected to the N-photon subspace, will exhibit an EP of order N+1. By using the well-established mapping between a beam splitter Hamiltonian and the perfect state transfer model in the photon-number space, we analytically obtain the time evolution of a general N-photon state and numerically simulate the system’s evolution in the post-selected manifold. Our results pave the way toward realizing robust, arbitrary-order EPs on demand in a single device.

    1. INTRODUCTION

    A fundamental postulate of quantum theory is that the Hamiltonian of a (closed) system is Hermitian, which guarantees real energy eigenvalues and a unitary time evolution [1]. This conventional wisdom was upended when Bender and coworkers discovered families of non-Hermitian Hamiltonians with real spectra [2]. The common feature of all such Hamiltonians was that they were invariant under the combined operations of space- and time-reflection, i.e., they were parity and time-reversal (PT) symmetric Hamiltonians. Over the past two decades, it has become clear that non-Hermitian, PT-symmetric Hamiltonians represent classical systems with spatially or temporally separated gain and loss [35]. The spectrum of a PT symmetric Hamiltonian changes from purely real to complex conjugate pairs when the strength of its anti-Hermitian part matches the Hermitian energy scale. This PT-symmetry breaking transition occurs at an exceptional point [68]. The phenomenology of PT-symmetric Hamiltonians with second- and third-order EPs has been extensively explored in optical, mechanical, electrical, and acoustic experimental realizations [5]. Exceptional point degeneracies also occur in mode-selective lossy Hamiltonians, which has enabled investigations of EP-related phenomena in dissipative systems in the classical [9,10] and quantum [11,12] domains, including the realization of a fourth-order EP with single photons [13].

    Many remarkable properties of non-Hermitian systems, such as asymmetric mode switching [14], topological energy transfer [15], robust wireless power transfer [16], and enhanced classical sensitivity [1720] are due to their EP degeneracies. In sharp contrast with Hermitian Hamiltonians, whose eigenmodes continue to span the space irrespective of eigenvalue degeneracies, in the non-Hermitian case, the eigenmodes of the Hamiltonian at an EP do not span the space, and the deficit grows proportional to the order of the EP. This key difference is instrumental to the system sensitivity that scales with the order of the EP [19,20] and has led to the tremendous interest in developing systems with higher-order EPs [2123] and understanding their fundamental quantum limits [24]. However, experimentally realizing classical or quantum systems with higher-order EPs has proven extremely challenging, and EPs beyond the fourth order have not been realized. In particular, integrated platforms where EPs of different orders can be realized are absent.

    In this paper, we propose and theoretically investigate such a platform in a single, lossy waveguide beam splitter in the quantum domain. When excited by a state with N indistinguishable photons and confined to the N-photon subspace, we show that the dynamics of such a beam splitter has an EP of order N+1, which is observable with currently available number-resolving single-photon detectors [2527]. In contrast with the past proposals with multiple waveguides or resonators, where precise parameter tuning is needed to ensure that the higher-order EP does not split into lower-order ones, we show that these EPs are robust due to the bosonic nature of photons and linear nature of the loss at low intensities.

    The paper is structured as follows. In Section 2, we present the formal treatment of the model in the photon-number basis and show analytical results for the time evolution of an arbitrary state. Results from numerical simulations based on the so-called NOON states are presented in Section 3. The paper is concluded in Section 4 with a brief discussion.

    2. LOSSY BEAM SPLITTER IN THE PHOTON-NUMBER BASIS

    The general beam splitter Hamiltonian in second-quantized notation is given by [28] where a^ (a^) and b^ (b^) represent bosonic creation (annihilation) operators for photonic modes in the two waveguides, ω0 is their common propagation constant, the coupling between the two waveguides is given by κ, and Γ is the dissipation coefficient of the lossy waveguide.

    To unveil the link between the waveguide beam splitter and arbitrary-order exceptional points, we represent the Hamiltonian Eq. (1) in the two-mode, N-photon subspace. This subspace is spanned by N+1 orthonormal states |m)|Nm,m=|Nma|mb (0mN) corresponding (Nm) photons in the neutral waveguide and m photons in the lossy waveguide. We emphasize that photon-number-resolving detection is necessary to access different basis states in this subspace [Fig. 1(a)]. It follows from the properties of creation and annihilation operators that the loss-term gives rise to linearly varying potential, i.e., (m|b^b^|m)=mδmm [Fig. 1(b)], while the tunneling term gives rise to index-dependent, nearest-neighbor tunneling amplitude, i.e., (m|a^b^+b^a^|m)=(δm,m+1+δm,m1)2κm(N+1m) [29]. In this basis, Eq. (1) becomes [29,30] where N^=a^a^+b^b^ is the total photon-number operator. The Hermitian generators satisfy the angular momentum algebra [J^z,J^x]=iJ^y, with J^y=i(a^b^a^b^)/2, and its cyclic permutations. Thus, in the N-photon subspace, J^x and J^z are spin S=N/2 representations of the angular momentum operators. In the absence of loss, Γ=0, Eq. (2) reduces to the perfect-state transfer model [3136].

    (a) Schematic of a single, lossy waveguide beam splitter excited with N indistinguishable photons prepared in the state |m)≔|N−m,m⟩=|N−m⟩a|m⟩b, where a represents the neutral (gray) waveguide and b is the lossy (red) waveguide. (b) Mapping onto the N-photon subspace spanned by (N+1) multiphoton states |m), represented as a tight-binding lattice model. The coupling between adjacent “modes” is given by matrix elements of J^x; the linearly increasing loss is also shown. (c) Flow of eigenvalues of H^N for N=4. R(λr) shows level attraction with an EP of order five at Γ=2κ; I(λr) shows the emergence of slow modes past the transition. (d) Intensity I(z) shows the fraction of trials where the system remains in the N-photon subspace, i.e., the post-selection probability. It reflects the order of the exceptional point. The beam splitter parameters are ω0=κ=1 cm−1, and the initial state is |ψ(0)⟩=|0).

    Figure 1.(a) Schematic of a single, lossy waveguide beam splitter excited with N indistinguishable photons prepared in the state |m)|Nm,m=|Nma|mb, where a represents the neutral (gray) waveguide and b is the lossy (red) waveguide. (b) Mapping onto the N-photon subspace spanned by (N+1) multiphoton states |m), represented as a tight-binding lattice model. The coupling between adjacent “modes” is given by matrix elements of J^x; the linearly increasing loss is also shown. (c) Flow of eigenvalues of H^N for N=4. R(λr) shows level attraction with an EP of order five at Γ=2κ; I(λr) shows the emergence of slow modes past the transition. (d) Intensity I(z) shows the fraction of trials where the system remains in the N-photon subspace, i.e., the post-selection probability. It reflects the order of the exceptional point. The beam splitter parameters are ω0=κ=1  cm1, and the initial state is |ψ(0)=|0).

    Within the (N+1)-dimensional subspace, the equidistant eigenvalues of the Hamiltonian are analytically given by the expression where r={S,S+1,,S}. It follows from Eq. (5) that the adjacent difference Δλλrλr1=4κ2Γ2 is purely real when the dissipation coefficient Γ2κ and becomes purely imaginary when the dissipation coefficient is larger, i.e., Γ>2κ. At the transition point Γc=2κ, all eigenvalues become degenerate, and all the eigenmodes coalesce, thus giving rise to an exceptional point of order N+1 [30]. Figure 1(c) shows the (analytical) flow of eigenvalues for a lossy beam splitter excited by N=4 photons, a realistic number that has been achieved in recent experiments [2527]. The beam splitter parameters are set to ω0=κ=1  cm1. The top panel in Fig. 1(c) shows that, as Γ increases, R(λr) undergo level attraction and become degenerate at Γc, remaining constant thereafter. The bottom panel shows that I(λr) increase linearly with Γ and are the same for all eigenmodes for ΓΓc, whereas past the transition point, slowly (and rapidly) decaying eigenmodes emerge. It is worth noting that, in a waveguide beam splitter, the EP of order N+1 appears naturally in the N-photon subspace, and it is always located at Γc=2κ irrespective of N.

    To detect the order of the EP in an experimentally friendly manner [37], we consider the behavior of the intensity I(z) within the N-photon subspace as a function of the propagation distance z, or equivalently, the time. In general, when the lossy beam splitter is excited with an N-photon input, the number-resolving detectors at the output will register any of the (N+1)(N+2)/2 possibilities |pa|qb where 0p,qN with p+qN. Thus, I(z) registers the fraction of trials where the total number of photons detected is exactly N, i.e., we post-select on the manifold where no photons are absorbed in the lossy waveguide [11]. For a normalized initial state |ψ(0), this intensity is given by where G(z)=exp(iH^Nz) is the decaying time evolution operator. At the exceptional point Γ=2κ, the Hamiltonian H^N satisfies the characteristic equation [λr(ω0iκ)N]N+1=0; therefore, the power-series expansion for G(z) terminates at Nth order in z. This implies the post-selection probability I(z)z2Nexp(NΓz) at long distance κz1. The numerically obtained post-selection probability I(z), for input states where N photons are injected into the neutral waveguide for N{5,7,10}, is shown in Fig. 1(d). It clearly shows that the order of the exceptional point is reflected in the results. We note that, under realistic conditions, the ability to resolve EPs of different orders will be adversely affected by noise floor in the number-resolving detectors, nonunit fidelity of the initial state, and the z-variation of the coupling and loss Γ parameters [38,39]. The former will strongly affect subspaces with small N, whereas the noise in the latter does not depend on N.

    To provide visual insight into the EP of order N+1, we consider an (N+1)-dimensional normalized, complex eigenvector |λr of the Hamiltonian [Eq. (2)] and represent it by its three spin-projection components Jα,r=λr|J^α|λr (α=x,y,z). Note that the resulting 3D vectors Jr=(Jx,r,Jy,r,Jz,r) are not orthogonal to each other in the Hermitian limit. Figure 2 shows the evolution of these vectors as a function of loss Γ for N=4 (left column) and N=5 (right column) photon-subspaces. When Γ=Γc/2 [Figs. 2(a) and 2(b)], the spin-projections lie in the xy plane; at Γ=0.99Γc [Figs. 2(c) and 2(d)], the eigenvectors start to coalesce, and the spin-projections point along the positive y axis, both indicating the close proximity to an EP. Past the exceptional point, at Γ=1.5Γc [Figs. 2(e) and 2(f)], the spin-projections now lie in the xz plane.

    Evolution of spin-projections Jr for the N+1 eigenmodes in the post-selected manifold with N=4 (left column) and N=5 (right column) photons, considering different values of the dissipation coefficient: (a), (b) Γ=Γc/2, (c), (d) Γ=0.99Γc, and (e), (f) Γ=1.5Γc. The vector coordinates in the (Jx,Jy,Jz) space are defined by the expectation values of the J^α operators in each eigenstate. When Γ<Γc, the spin projections are in the x−y plane; at the EP, they coalesce along the positive y axis; and when Γ>Γc, they are in the x−z plane.

    Figure 2.Evolution of spin-projections Jr for the N+1 eigenmodes in the post-selected manifold with N=4 (left column) and N=5 (right column) photons, considering different values of the dissipation coefficient: (a), (b) Γ=Γc/2, (c), (d) Γ=0.99Γc, and (e), (f) Γ=1.5Γc. The vector coordinates in the (Jx,Jy,Jz) space are defined by the expectation values of the J^α operators in each eigenstate. When Γ<Γc, the spin projections are in the xy plane; at the EP, they coalesce along the positive y axis; and when Γ>Γc, they are in the xz plane.

    In order to obtain the photon-number-resolved population dynamics of the lossy beam splitter, we note that the time-evolution operator satisfies the Schrödinger equation izG(z)=H^NG(z). Therefore, it can be expressed in terms of the total number operator N^ and the angular momentum operators J^α (with α=x,y,z) by using the Wei–Norman method [40,41], where J^±=J^x±iJ^y are the angular momentum raising and lowering operator, respectively. Note that, because the photon-number operator N^ commutes with J^α, it can be treated as a c-number. f±(z),fz(z) are three complex functions that parameterize the nonunitary time evolution operator and satisfy the following set of coupled, nonlinear differential equations: The solutions of Eq. (8), subject to the initial condition G(0)=1, or equivalently f±(0)=fz(0)=0, are given by We note that the functions f+(z)=f(z) are real, irrespective of whether Δλ is real or purely imaginary, while fz(z) is, in general, complex. It is straightforward to check that, as the system approaches the exceptional point, i.e., Δλ0, the functions f±(z)κz/(1+κz) approach unity at κz1. On the other hand, the function fz(z)2iln(1+κz), in conjunction with the diagonal operator J^z=diag(N/2,,N/2), gives rise to an algebraically growing time evolution operator G(z)zNexp(NΓz/2). Thus, our exact solution, Eq. (7), encodes the order of the exceptional point.

    3. POST-SELECTED DYNAMICS: NUMERICAL RESULTS

    Motivated by the realization of high-order multiphoton entangled states [42,43], we explore the dynamics of the lossy beam splitter excited with NOON-state initial conditions, i.e., |ϕ(0)=(|Na|0b+|0a|Nb)/2=(|0)+|N))/2. Although the post-selection probability I(z)=ϕ(z)|ϕ(z) decreases exponentially with the propagation distance, we will see that, within the post-selected N-photon manifold, signatures of the PT-symmetry breaking transition and the order of the exceptional point are clearly visible. To that end, we consider the normalized, z-dependent occupation function which satisfies m=0NP(|m);z)=1. Figure 3 shows the results for this occupation function for two different input states. The top row in Fig. 3 shows the normalized mode occupations as a function of z for an N=5 state. When Γ=Γc/4 [Fig. 3(a)], we see an asymmetric, oscillatory motion across the six modes with an energy flow from a low-loss region to the high-loss region. At the EP, Γ=Γc, the system reaches a steady state with a weight distributed largely in the low-loss region [Fig. 3(b)]. Past the transition, Γ=1.2Γc, the steady-state is reached slower [Fig. 3(c)], indicating the emergence of slowly decaying eigenmodes for the Hamiltonian Eq. (2). The bottom row in Fig. 3 shows corresponding results for an N=8 NOON state input. Comparing the two rows, it is clear that the period of asymmetric oscillations does not depend on N and the order of the EP is reflected in the post-selected, N-photon manifold results.

    Mode occupation dynamics in the post-selected manifold with NOON state input. (a) For N=5 and small loss, the dynamics show asymmetric oscillations. (b) At the EP, P(|m),z) reaches a steady state with most of the weight localized in the low-loss region. (c) After the transition, the steady-state is reached more slowly. (d)–(f) show qualitatively similar results for an N=8 NOON state input. The waveguide beam splitter parameters are set to ω0=κ=1 cm−1.

    Figure 3.Mode occupation dynamics in the post-selected manifold with NOON state input. (a) For N=5 and small loss, the dynamics show asymmetric oscillations. (b) At the EP, P(|m),z) reaches a steady state with most of the weight localized in the low-loss region. (c) After the transition, the steady-state is reached more slowly. (d)–(f) show qualitatively similar results for an N=8 NOON state input. The waveguide beam splitter parameters are set to ω0=κ=1  cm1.

    4. DISCUSSION

    Despite tremendous interest due to the classical sensitivity enhancement they offer [1720], experimental realizations of exceptional points of higher order have remained elusive. The primary obstacle for such realizations in coupled waveguides, resonators, or other traditional platforms is the fine tuning of system parameters, which is required by a higher-order symmetry necessary for eigenmode degeneracy [44]. When such stringent constraints regarding the ratio of losses or nearest-neighboring coupling amplitudes are not satisfied, a higher-order EP splits into EPs of lower order.

    Here, we have shown that a single, lossy waveguide beam splitter can be used to realize robust EPs of arbitrary order without any fine tuning required. In our proposal, the stringent symmetries required for higher-order EPs are guaranteed by the bosonic nature of input photons, and the linear nature of loss. We have shown that the dynamics observed within the post-selected N-photon subspace has an EP of order N+1. Thus, our analysis passes the burden of fine-tuning the Hermitian and lossy parts of the Hamiltonian onto the dual tasks of creating higher-order NOON states and number-resolving photon detectors, a rapidly maturing technology found in quantum optics laboratories across the globe [2527,45]. Our results, therefore, offer a realistic pathway for realizing EPs of arbitrary order on demand in a single platform.

    Acknowledgment

    Acknowledgment. M. A. Q. J. and R. J. L. M. were supported by CONACyT, and Y. N. J. was supported by the NSF. B. M. R. L. acknowledges financial support from the Marcos Moshinsky Foundation through the 2018 Marcos Moshinsky Young Researcher Chair. O. S. M. L. acknowledges startup funding from Louisiana State University.

    References

    [1] L. Landau, E. Lifshitz. Quantum Mechanics: Non Relativistic Theory, 3(1977).

    [2] C. M. Bender, S. Boettcher. Real spectra in non-Hermitian Hamiltonians having PT symmetry. Phys. Rev. Lett., 80, 5243-5246(1998).

    [3] Y. N. Joglekar, C. Thompson, D. D. Scott, G. Gautam. Optical waveguide arrays: quantum effects and PT symmetry breaking. Eur. Phys. J., 63, 30001(2013).

    [4] L. Feng, R. El-Ganainy, L. Ge. Non-Hermitian photonics based on parity-time symmetry. Nat. Photonics, 11, 752-762(2017).

    [5] R. El-Ganainy, K. G. Makris, M. Khajavikhan, Z. H. Musslimani, S. Rotter, D. N. Christodoulides. Non-Hermitian physics and PT symmetry. Nat. Phys., 14, 11-19(2018).

    [6] M. Mueller, I. Rotter. Exceptional points in open quantum systems. J. Phys. A, 41, 244018(2008).

    [7] W. D. Heiss. The physics of exceptional points. J. Phys. A, 45, 444016(2012).

    [8] T. Kato. Perturbation Theory for Linear Operators, 132(2013).

    [9] A. Guo, G. J. Salamo, D. Duchesne, R. Morandotti, M. Volatier-Ravat, V. Aimez, G. A. Siviloglou, D. N. Christodoulides. Observation of PT-symmetry breaking in complex optical potentials. Phys. Rev. Lett., 103, 093902(2009).

    [10] K. Ding, G. Ma, M. Xiao, Z. Q. Zhang, C. T. Chan. Emergence, coalescence, and topological properties of multiple exceptional points and their experimental realization. Phys. Rev. X, 6, 021007(2016).

    [11] M. Naghiloo, M. Abbasi, Y. N. Joglekar, K. Murch. Quantum state tomography across the exceptional point in a single dissipative qubit(2019).

    [12] J. Li, A. K. Harter, J. Liu, L. de Melo, Y. N. Joglekar, L. Luo. Observation of parity-time symmetry breaking transitions in a dissipative Floquet system of ultracold atoms. Nat. Commun., 10, 855(2019).

    [13] Z. Bian, L. Xiao, K. Wang, X. Zhan, F. A. Onanga, F. Ruzicka, W. Yi, Y. N. Joglekar, P. Xue. Time invariants across a fourth-order exceptional point in a parity-time-symmetric qudit(2019).

    [14] J. Doppler, A. A. Mailybaev, J. Bohm, U. Kuhl, A. Girschik, F. Libisch, T. J. Milburn, P. Rabl, N. Moiseyev, S. Rotter. Dynamically encircling an exceptional point for asymmetric mode switching. Nature, 537, 76-79(2016).

    [15] H. Xu, D. Mason, L. Jiang, J. G. E. Harris. Topological energy transfer in an optomechanical system with exceptional points. Nature, 537, 80-83(2016).

    [16] S. Assawaworrarit, X. Yu, S. Fan. Robust wireless power transfer using a nonlinear parity-time-symmetric circuit. Nature, 546, 387-390(2017).

    [17] J. Wiersig. Enhancing the sensitivity of frequency and energy splitting detection by using exceptional points: application to microcavity sensors for single-particle detection. Phys. Rev. Lett., 112, 203901(2014).

    [18] J. Wiersig. Sensors operating at exceptional points: general theory. Phys. Rev. A, 93, 033809(2016).

    [19] H. Hodaei, A. U. Hassan, S. Wittek, H. Garcia-Gracia, R. El-Ganainy, D. N. Christodoulides, M. Khajavikhan. Enhanced sensitivity at higher-order exceptional points. Nature, 548, 187-191(2017).

    [20] W. Chen, S. Kaya Ozdemir, G. Zhao, J. Wiersig, L. Yang. Exceptional points enhance sensing in an optical microcavity. Nature, 548, 192-196(2017).

    [21] H. Zhao, Z. Chen, R. Zhao, L. Feng. Exceptional point engineered glass slide for microscopic thermal mapping. Nat. Commun., 9, 1764(2018).

    [22] Q. Zhong, D. N. Christodoulides, M. Khajavikhan, K. G. Makris, R. El-Ganainy. Power-law scaling of extreme dynamics near higher-order exceptional points. Phys. Rev. A, 97, 020105(2018).

    [23] S. Wang, B. Hou, W. Lu, Y. Chen, Z. Q. Zhang, C. T. Chan. Arbitrary order exceptional point induced by photonic spin-orbit interaction in coupled resonators. Nat. Commun., 10, 832(2019).

    [24] H.-K. Lau, A. A. Clerk. Fundamental limits and non-reciprocal approaches in non-Hermitian quantum sensing. Nat. Commun., 9, 4320(2018).

    [25] G. Harder, T. J. Bartley, A. E. Lita, S. W. Nam, T. Gerrits, C. Silberhorn. Single-mode parametric-down-conversion states with 50 photons as a source for mesoscopic quantum optics. Phys. Rev. Lett., 116, 143601(2016).

    [26] L. Banchi, W. S. Kolthammer, M. S. Kim. Multiphoton tomography with linear optics and photon counting. Phys. Rev. Lett., 121, 250402(2018).

    [27] O. S. Magaña-Loaiza, R. de J. León-Montiel, A. Perez-Leija, A. B. U’Ren, C. You, K. Busch, A. E. Lita, S. W. Nam, R. P. Mirin, T. Gerrits. Multiphoton quantum-state engineering using conditional measurements(2019).

    [28] W. K. Lai, V. Buek, P. L. Knight. Nonclassical fields in a linear directional coupler. Phys. Rev. A, 43, 6323-6336(1991).

    [29] K. Tschernig, R. de J. León-Montiel, O. S. Magaña-Loaiza, A. Szameit, K. Busch, A. Perez-Leija. Multiphoton discrete fractional Fourier dynamics in waveguide beam splitters. J. Opt. Soc. Am. B, 35, 1985-1989(2018).

    [30] E. M. Graefe, U. Gnther, H. J. Korsch, A. E. Niederle. A non-Hermitian PT symmetric Bose-Hubbard model: eigenvalue rings from unfolding higher-order exceptional points. J. Phys. A, 41, 255206(2008).

    [31] M. Christandl, N. Datta, A. Ekert, A. J. Landahl. Perfect state transfer in quantum spin networks. Phys. Rev. Lett., 92, 187902(2004).

    [32] A. Perez-Leija, R. Keil, A. Kay, H. Moya-Cessa, S. Nolte, L.-C. Kwek, B. M. Rodriguez-Lara, A. Szameit, D. N. Christodoulides. Coherent quantum transport in photonic lattices. Phys. Rev. A, 87, 012309(2013).

    [33] A. Perez-Leija, R. Keil, H. Moya-Cessa, A. Szameit, D. N. Christodoulides. Perfect transfer of path-entangled photons in Jx photonic lattices. Phys. Rev. A, 87, 022303(2013).

    [34] Y. N. Joglekar, A. Saxena. Robust PT-symmetric chain and properties of its Hermitian counterpart. Phys. Rev. A, 83, 050101(2011).

    [35] Y. N. Joglekar, C. Thompson, G. Vemuri. Tunable waveguide lattices with nonuniform parity-symmetric tunneling. Phys. Rev. A, 83, 063817(2011).

    [36] R. J. Chapman, M. Santandrea, Z. Huang, G. Corrielli, A. Crespi, M.-H. Yung, R. Osellame, A. Peruzzo. Experimental perfect state transfer of an entangled photonic qubit. Nat. Commun., 7, 11339(2016).

    [37] R. de J. León-Montiel, M. A. Quiroz-Juárez, J. L. Domínguez-Juárez, R. Quintero-Torres, J. L. Aragón, A. K. Harter, Y. N. Joglekar. Observation of slowly decaying eigenmodes without exceptional points in Floquet dissipative synthetic circuits. Commun. Phys., 1, 88(2018).

    [38] N. A. Mortensen, P. A. D. Gonçalves, M. Khajavikhan, D. N. Christodoulides, C. Tserkezis, C. Wolff. Fluctuations and noise-limited sensing near the exceptional point of parity-time-symmetric resonator systems. Optica, 5, 1342-1346(2018).

    [39] S. Longhi. Loschmidt echo and fidelity decay near an exceptional point(2019).

    [40] J. Wei, E. Norman. Lie algebraic solution of linear differential equations. J. Math. Phys., 4, 575-581(1963).

    [41] W. H. Louisell. Quantum Statistical Properties of Radiation, 7(1973).

    [42] I. Afek, O. Ambar, Y. Silberberg. High-NOON states by mixing quantum and classical light. Science, 328, 879-881(2010).

    [43] J. Zhang, M. Um, D. Lv, J.-N. Zhang, L.-M. Duan, K. Kim. NOON states of nine quantized vibrations in two radial modes of a trapped ion. Phys. Rev. Lett., 121, 160502(2018).

    [44] M. Teimourpour, Q. Zhong, M. Khajavikhan, R. El-Ganainy. Higher order exceptional points in discrete photonics platforms. Parity-Time Symmetry and Its Applications, 261-275(2018).

    [45] A. E. Lita, A. J. Miller, S. W. Nam. Counting near-infrared single-photons with 95% efficiency. Opt. Express, 16, 3032-3040(2008).

    Mario A. Quiroz-Juárez, Armando Perez-Leija, Konrad Tschernig, Blas M. Rodríguez-Lara, Omar S. Magaña-Loaiza, Kurt Busch, Yogesh N. Joglekar, Roberto de J. León-Montiel, "Exceptional points of any order in a single, lossy waveguide beam splitter by photon-number-resolved detection," Photonics Res. 7, 862 (2019)
    Download Citation