Abstract
In non-Hermitian settings, the particular position at which two eigenstates coalesce in the complex plane under a variation of a physical parameter is called an exceptional point. An open disordered system is a special class of non-Hermitian system, where the degree of scattering directly controls the confinement of the modes. Herein a non-perturbative theory is proposed which describes the evolution of modes when the permittivity distribution of a 2D open dielectric system is modified, thereby facilitating to steer individual eigenstates to such a non-Hermitian degeneracy. The method is used to predict the position of such an exceptional point between two Anderson-localized states in a disordered scattering medium. We observe that the accuracy of the prediction depends on the number of localized states accounted for. Such an exceptional point is experimentally accessible in practically relevant disordered photonic systems.
© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement
1. Introduction
Most physical systems are subject to inherent losses, either because of dissipation or openness, and are described mathematically by a non-Hermitian Hamiltonian. This generally implies that the eigenstates of the Hamiltonian are complex as well as non-orthogonal. In such systems, the interrelation between pairs of eigenstates under the variation of a set of external parameters is essentially dictated by the existence of exceptional points (EPs). At an EP not only eigenvalues but also eigenstates degenerate, i.e. eigenvectors have the same amplitudes and a constant phase shift. The eigenvalue surfaces display a singular topology in the presence of an EP [1] while the eigenstates accumulate a residual geometric phase upon encircling the EP in parameter space [2,3]. Since their introduction by Kato in 1966 [4], EPs have been associated with numerous often counterintuitive physical effects [5–8] that include mode hybridization [9], quantum phase transitions [10], reversal of the pump dependence of a laser [11], Parity-Time (PT) symmetry breaking [12,13], increased sensitivity to perturbations [14,15] or spectrally broad coherent perfect absorption [16]. They have been observed experimentally in a variety of different systems such as microwave billiards [17], chaotic optical microcavities [18], optomechanical arrangements [19], Fabry-Pérot resonators [20], two level atoms in high-Q cavities [8], coupled photonic-crystals [21], coupled acoustic cavities [22], optical fibres [23] or electronic circuits [24].
Open random media constitute a special class of non-Hermitian systems. Here, modal confinement depends almost entirely on the degree of scattering. In the regime of strong scattering, modes can experience Anderson-localization—a multi-scattering interference phenomenon where the spatial extension of the modes becomes smaller than the system size resulting in transport inhibition [25]. These disorder-induced localized states, for example, provide a natural optical cavity in random lasers [26–28] and showed to be good candidates for cavity quantum electrodynamics (QED) [29,30], with the main advantage of being inherently disorder-robust. By locally tuning the disorder, e.g. by gradually adjusting the dielectric permittivity, these states can be engineered to spatially and spectrally overlap to form so-called necklace states [31–33], which represent open channels in a nominally localized system [34,35]. The formation of necklace states is considered a key mechanism in the transition from localization to the diffusive regime [36]. While EPs and, in general, PT-symmetry have been studied in the context of disordered media [21,28,37–41], coalescing localized modes and their associated degeneracies have not been studied so far.
In this work, the coalescence of two Anderson-localized optical modes in a two-dimensional (2D) dielectric random system is invoked by manipulating the dielectric permittivity at two different locations in the system. We first propose a general analytical analysis to follow the spectral and spatial evolution of modes in 2D dielectric open media. Then, this approach is applied to the specific case of Anderson-localized modes to identify the position of an EP in the parameter space. The prediction is tested numerically with finite element method (FEM) simulations. Remarkably, in this open disordered system a conventional two-mode model, prominently used in a variety of EP-related studies, would fail badly. The modification of a few single scatterers affects a multitude of modes and eventually influences the position of the degeneracies. This emphasizes the complexity of eigenstate interactions in disordered media. We believe that our approach opens the way to a controlled local manipulation of the permittivity and the possibility to engineer hybridized and coupled localized modes. Furthermore, we think this approach can be easily extended to other kinds of resonant systems, e.g. coupled arrays of cavities [42,43].
2. Theory and applications
2.1 Generalized eigenvalue problem
First, we consider the general case of a finite-size dielectric medium in 2D space with an inhomogeneous dielectric constant distribution $\epsilon (\mathbf {r})$. The speed of light is set to $c$ = 1 such that frequencies are hence measured in multiples of $c$ having a unit of $\text {m}^{-1}$. The electric field is assumed transverse and thus its polarization is neglected such that it satisfies the scalar Helmholtz equation
where $E(\mathbf {r},\omega )$ stands for the scalar electric field in the frequency domain. Eigensolutions of Eq. (1) define the modes or eigenstates of the problemNow, consider two locations $\mathbf {R}_1$ and $\mathbf {R}_2$ where the permittivity is varied
The permittivity distribution $\tilde {\epsilon }(\mathbf {r})$ describes a new disordered system associated with a new set of eigenstates $(\tilde {\Omega }_i, |{\tilde {\Psi }_i}\rangle )_{i\in \mathbb {N}}$. The electric field of the modified system written in terms of the basis of the original random system reads
where $b_i(\omega )$ are the new expansion coefficients. Inserting Eq. (7) into Eq. (6) givesProjecting $\langle {\Psi _j^{*}}|$ onto Eq. (8) using Eq. (2) and the biorthogonal product in Eq. (4) leads to
If we consider a finite set of $N$ modes, Eq. (9) can be written conveniently in the form of a generalized eigenvalue problem
Equation (11) extends this result to any number of interacting modes $N>2$ essentially describing a network of linearly coupled oscillators.
2.2 Identifying an EP between two Anderson-localized states
The formalism presented above is now applied to a 2D random collection of 896 circular dielectric scatterers (radius $60~\text {nm}$) with dielectric permittivity $\epsilon = 4$ embedded in a host material of index $\epsilon _{\text {mat}} = 1$, with a filling fraction of $40\%$ [Fig. 1(a)]. The system dimensions are $L\times L = 5.3\,\times 5.3 \,\mu \text {m}^{2}$. The interference effects, originating from multiple scattering leads to spatial confinement by the disorder known as Anderson localization and characterized by a localization length $\xi$. In the spectral range considered in the following, $\xi$ is computed through an averaged correlation between modes [49] and estimated to be around $1 \mu \text {m} \ll L$. Thus, the energy of the modes is mostly confined in the system: The modes are localized. The two circular regions with $340~\text {nm}$ diameter centered at $\mathbf {R}_1$ and $\mathbf {R}_2$ are shown in Fig. 1(a) as red and green circles, respectively. Within these regions the dielectric permittivity of the scatterers is varied from $\epsilon$ to $\epsilon + \Delta \epsilon _1$ and $\epsilon + \Delta \epsilon _2$, respectively.
The initial modes $(\Omega _i,\, |{\Psi _i}\rangle )$, $i \in 1,\dots,N$, which are the only input requested by Eq. (11), are computed using FEM simulation [50,51] with absorbing boundary conditions that are placed $0.4\, {\mu }\text {m}$ away from each side of the system. A large number of modes ($N = 90$) are computed for the initial system [Fig. 1(b)] in a narrow spectral range and we check that none of them are degenerate. Two localized states $|{\Psi _1}\rangle$ and $|{\Psi _2}\rangle$, corresponding to $\Omega _1$ and $\Omega _2$ respectively, are selected for being spectrally close [Fig. 1(b)] but spatially distinct [Fig. 1(c),(d)]. The spectral distance between mode $i$ and mode $1$ is defined as $d(1,i) = |\Omega _1 - \Omega _i|$ and is visualized in the color-coding of Fig. 1(b). Here, mode $2$ has the largest spectral overlap with mode $1$ but we will see that the influence of other nearby modes cannot be neglected in the modal interaction. The initial modes obtained via a FEM are spatially defined in a finite spatial domain $V$ while the biorthogonal product of Eq. (4) is an integral over all $\mathbb {R}^{2}$. However, it can be split into an integral over $V$ and a second term at the boundary of $V$ which replaces outside propagation [52]. Because the modes are localized their amplitudes along the boundary of $V$ are very small. Neglecting the edge term in the biorthogonal product leads to an inaccuracy of around $0.8\%$ in the ensuing calculation of the EP position. Boundary terms can therefore be safely neglected in the biorthogonal product.
Now, the permittivity distribution $\epsilon (\mathbf {r})$ is altered by $\left (\Delta \epsilon _1,\Delta \epsilon _2\right )$ and the corresponding eigensystem is calculated from Eq. (11) for $N=60$ interacting modes for each pair $\left (\Delta \epsilon _1,\Delta \epsilon _2\right )$ in this two-dimensional parameter space, from which a new set of modes $(\tilde {\Omega }_i,\, |{\tilde {\Psi }_i}\rangle )$ is obtained. The spectral distance $\tilde {d}(1,2) = |\tilde {\Omega }_1 - \tilde {\Omega }_2|$ is shown in Fig. 2(a) for all values of the pair $\left (\Delta \epsilon _1,\Delta \epsilon _2\right )$. The sharp drop of $\tilde {d}(1,2)$ to zero observed at $(\Delta \epsilon _1,\Delta \epsilon _2)_{\text {EP}} = (0.939,0.90)$ reveals the existence of an EP, where both eigenvalues become degenerate. The evolution in the parameter space of the real and imaginary part of the eigenvalues, $\tilde {\Omega }_1$ and $\tilde {\Omega }_2$, is represented separately in Fig. 2(c),(d). The intricate topology of self-intersecting Riemann sheets at precisely the position $(\Delta \epsilon _1,\Delta \epsilon _2)_{\text {EP}}$, corroborates the existence of an EP, where both the real and imaginary part of the eigenvalues cross. At this singular position, the corresponding eigenfunctions become collinear. This is demonstrated in Fig. 2(b), where the function $(1-|\langle {\tilde {\Psi }_1}|{\tilde {\Psi }_2}\rangle |)^{-1}$ is found to diverge.
2.3 Testing the theoretical prediction
To substantiate the analytical prediction, the first two eigenvectors ($|{\tilde {\Psi }_1}\rangle$, $|{\tilde {\Psi }_2}\rangle$) of the modified permittivity distribution $\tilde {\epsilon }(\mathbf {r})$ are calculated via numerical FEM simulations for each point $(\Delta \epsilon _1,\Delta \epsilon _2)$ on a grid in a range that encloses $(\Delta \epsilon _1,\Delta \epsilon _2)_{\text {EP}}$. The sampling is set to $0.04$ along the two parameters $\Delta \epsilon _1$ and $\Delta \epsilon _2$. The numerically computed eigenvalues merge at $(\Delta \epsilon _1,\Delta \epsilon _2) = (0.92,0.88) \pm (0.04,0.04)$, such that the predicted value from Eq. (11) is within the error bars. The collinearity of the eigenvectors is also numerically confirmed (not shown). Finally, the coalescing mechanism of modes 1 and 2 is described in Fig. 2(e) where modes with identical real (imaginary) parts are plotted in the parameter space, $(\Delta \epsilon _1,\Delta \epsilon _2)$, as blue (red) dots. At the EP where real and imaginary parts are equal simultaneously, the red and blue trajectories meet. Near this position, we compute the real (imaginary) part of the two eigenvectors, $\tilde {\Psi }_1(x,y)$ and $\tilde {\Psi }_2(x,y)$. Their distribution is shown in the insets of Fig. 2(e)]. Starting from the unperturbed modes $\Psi _1(x,y)$ and $\Psi _2(x,y)$ [shown in Fig. 1(c),(d)], both eigenvectors are progressively modified when approaching the EP. In the vicinity of the EP, $\text {Re}(\tilde {\Psi }_1)$ converges to $\text {Im}(\tilde {\Psi }_2)$ and $\text {Im}(\tilde {\Psi }_1)$ to $-\text {Re}(\tilde {\Psi }_2)$. Therefore, the eigenvectors display the same amplitude and a phase shift of $\pi /2$, $\tilde {\Psi }_1 = -i \tilde {\Psi }_2$, i.e. they become collinear at the EP. The amplitude of the degenerate eigenstate at the EP, $|{\tilde {\Psi }_{\text {EP}(1,2)}}\rangle$, is shown in the inset of Fig. 3(a). It forms a beaded chain which connects both ends of the system and is similar to necklace states studied in [31].
2.4 Influence of other localized modes
Here, in contrast to higher-order EPs [53], the EP under study characterizes the coalescence of only two modes. Yet, when the permittivity distribution gets modified to converge to this EP, the others modes play a role in the deformation of the two coalescing modes. To investigate how much different localized modes of the system contribute to the degenerate state $|{\tilde {\Psi }_{\text {EP}(1,2)}}\rangle$, we measure their spatial overlap with the degenerate eigenstate. The norm of their biorthogonal projection, $| \langle {\Psi _{i}^{*}\, }|{\tilde {\Psi }_{\text {EP}(1,2)}}\rangle |$, is shown in Fig. 3(a). Remarkably enough, it demonstrates the fading, though not negligible, influence of nearby modes with $i>2$. Both modes 1 and 2 contribute to $60\%$ to the degenerate EP state, while other modes have a vanishing contribution with increasing index, though modes 4 and 5, still contribute significantly for $25\%$. Their influence is also highlighted in Fig. 3(b), where the position of the EP, $(\Delta \epsilon _1,\Delta \epsilon _2)_{\text {EP}}$, is calculated using the analytical prediction, Eq. (11), for an increasing number of contributing modes $N$. Fluctuations of the EP position are large for small $N$, while the convergence is reached for value of $N$ larger than $55$. This observation emphasizes a fundamental point, namely that the spatial amplitude distribution of the degenerate state at the EP results from the coupling of a large number of modes.
3. Conclusion
In conclusion, we proposed a general analytical theory to study the evolution of modes in open media when the permittivity is varied. This non-perturbative approach relies on the linearity of $\epsilon (\mathbf {r})$ in the Helmholtz equation: For any variation $(\Delta \epsilon _1, \, \Delta \epsilon _2)$ a new set of modes can be computed from Eq. (11) simply from the knowledge of the initial system. This approach describes the system by an infinite set of modes acting as oscillators coupled via the modification of the permittivity. We considered the specific case of random media in the localization regime and show that our theory can be used to investigate the mode coupling and hybridization resulting from a local perturbation. In particular, by changing the local index at two different small areas, two modes are brought to coalescence and give rise to an EP. Remarkably, the accuracy of the theoretical prediction is shown to strongly depend on the number of modes considered. A large number of initial modes is required, which shows that the mode evolution is dictated by multimode interactions. In terms of experimental realization, such a manipulation of the disorder can be easily implemented on existing setups [27,54,55]. The permittivity landscape can be shaped reversibly, e.g. using laser illumination to switch on nonlinear index changes [56–58] or by moving one or several scatterers [59]. EPs can be calculated for any pair of modes and even generalized to hybridization of three or more eigenstates, which opens the way to the control of light-matter interaction in random media. The effective creation of necklace states, for instance, allows the formation of open transmission channels in opaque systems. Alternatively, using mode repulsion in the vicinity of an EP, the disorder can be engineered to increase the spatial confinement of the modes and consequently their Q-factor. Finally, we believe this approach can be extended to others types of complex optical systems, e.g. photonic crystal cavity arrays used for quantum simulation [42,43,60] where fabrication inaccuracy could be compensated a posteriori by such an external control.
Funding
Agence Nationale de la Recherche (10-IDEX-0001-02 PSL*, 12-BS09-003-01); Israel Science Foundation (1871/15, 2074/15, 2630/20); United States - Israel Binational Science Foundation (2015694); H2020 Marie Skłodowska-Curie Actions (840745).
Acknowledgments
We thank S. Rotter for enlightening and fruitful discussions.
Disclosures
The authors declare no conflicts of interest.
Data availability
Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.
References
1. W. D. Heiss and A. L. Sannino, “Avoided level crossing and exceptional points,” J. Phys. A: Math. Gen. 23(7), 1167–1178 (1990). [CrossRef]
2. R. S. Whitney and Y. Gefen, “Berry Phase in a Nonisolated System,” Phys. Rev. Lett. 90(19), 190402 (2003). [CrossRef]
3. A. Carollo, I. Fuentes-Guridi, M. F. Santos, and V. Vedral, “Geometric Phase in Open Systems,” Phys. Rev. Lett. 90(16), 160402 (2003). [CrossRef]
4. T. Kato, Perturbation theory for linear operators, vol. 132 of Classics in Mathematics (Springer-Verlag Berlin Heidelberg, Berlin, Heidelberg, 1995), 2nd ed.
5. W. D. Heiss, “Phases of wave functions and level repulsion,” Eur. Phys. J. D 7(1), 1–4 (1999). [CrossRef]
6. M. Liertzer, L. Ge, A. Cerjan, A. D. Stone, H. E. Türeci, and S. Rotter, “Pump-induced exceptional points in lasers,” Phys. Rev. Lett. 108(17), 173901 (2012). [CrossRef]
7. J. Feilhauer, A. Schumer, J. Doppler, A. A. Mailybaev, J. Böhm, U. Kuhl, N. Moiseyev, and S. Rotter, “Encircling exceptional points as a non-Hermitian extension of rapid adiabatic passage,” Phys. Rev. A 102(4), 040201 (2020). [CrossRef]
8. Y. Choi, S. Kang, S. Lim, W. Kim, J.-R. Kim, J.-H. Lee, and K. An, “Quasieigenstate Coalescence in an Atom-Cavity Quantum Composite,” Phys. Rev. Lett. 104(15), 153601 (2010). [CrossRef]
9. M. Kammerer, F. Merz, and F. Jenko, “Exceptional points in linear gyrokinetics,” Phys. Plasmas 15(5), 052102 (2008). [CrossRef]
10. W. D. Heiss and M. Müller, “Universal relationship between a quantum phase transition and instability points of classical systems,” Phys. Rev. E 66(1), 016217 (2002). [CrossRef]
11. M. Brandstetter, M. Liertzer, C. Deutsch, P. Klang, J. Schöberl, H. E. Türeci, G. Strasser, K. Unterrainer, and S. Rotter, “Reversing the pump dependence of a laser at an exceptional point,” Nat. Commun. 5(1), 4034 (2014). [CrossRef]
12. S. Klaiman, U. Günther, and N. Moiseyev, “Visualization of Branch Points in $\mathcal {P}\mathcal {T}$-Symmetric Waveguides,” Phys. Rev. Lett. 101(8), 080402 (2008). [CrossRef]
13. A. Guo, G. J. Salamo, D. Duchesne, R. Morandotti, M. Volatier-Ravat, V. Aimez, G. A. Siviloglou, and D. N. Christodoulides, “Observation of $\mathcal {P}\mathcal {T}$-Symmetry Breaking in Complex Optical Potentials,” Phys. Rev. Lett. 103(9), 093902 (2009). [CrossRef]
14. M. P. Hokmabadi, A. Schumer, D. N. Christodoulides, and M. Khajavikhan, “Non-Hermitian ring laser gyroscopes with enhanced Sagnac sensitivity,” Nature 576(7785), 70–74 (2019). [CrossRef]
15. Y.-H. Lai, Y.-K. Lu, M.-G. Suh, Z. Yuan, and K. Vahala, “Observation of the exceptional-point-enhanced Sagnac effect,” Nature 576(7785), 65–69 (2019). [CrossRef]
16. C. Wang, W. R. Sweeney, A. D. Stone, and L. Yang, “Coherent perfect absorption at an exceptional point,” Science 373(6560), 1261–1265 (2021). [CrossRef]
17. B. Dietz, H. L. Harney, O. N. Kirillov, M. Miski-Oglu, A. Richter, and F. Schäfer, “Exceptional Points in a Microwave Billiard with Time-Reversal Invariance Violation,” Phys. Rev. Lett. 106(15), 150403 (2011). [CrossRef]
18. S.-B. Lee, J. Yang, S. Moon, S.-Y. Lee, J.-B. Shim, S. W. Kim, J.-H. Lee, and K. An, “Observation of an Exceptional Point in a Chaotic Optical Microcavity,” Phys. Rev. Lett. 103(13), 134101 (2009). [CrossRef]
19. H. Xu, D. Mason, L. Jiang, and J. G. E. Harris, “Topological energy transfer in an optomechanical system with exceptional points,” Nature 537(7618), 80–83 (2016). [CrossRef]
20. Y. Zhu and L. Zhu, “Accessing the Exceptional Points in Coupled Fabry–Perot Resonators through Hybrid Integration,” ACS Photonics 5(12), 4920–4927 (2018). [CrossRef]
21. K. H. Kim, M. S. Hwang, H. R. Kim, J. H. Choi, Y. S. No, and H. G. Park, “Direct observation of exceptional points in coupled photonic-crystal lasers with asymmetric optical gains,” Nat. Commun. 7(1), 13893 (2016). [CrossRef]
22. K. Ding, G. Ma, M. Xiao, Z. Q. Zhang, and C. T. Chan, “Emergence, Coalescence, and Topological Properties of Multiple Exceptional Points and Their Experimental Realization,” Phys. Rev. X 6(2), 021007 (2016). [CrossRef]
23. A. Bergman, R. Duggan, K. Sharma, M. Tur, A. Zadok, and A. Alù, “Observation of anti-parity-time-symmetry, phase transitions and exceptional points in an optical fibre,” Nat. Commun. 12(1), 486 (2021). [CrossRef]
24. Y. Choi, C. Hahn, J. W. Yoon, and S. H. Song, “Observation of an anti-PT-symmetric exceptional point and energy-difference conserving dynamics in electrical circuit resonators,” Nat. Commun. 9(1), 2182 (2018). [CrossRef]
25. P. W. Anderson, “Absence of Diffusion in Certain Random Lattices,” Phys. Rev. 109(5), 1492–1505 (1958). [CrossRef]
26. V. Milner and A. Z. Genack, “Photon Localization Laser: Low-Threshold Lasing in a Random Amplifying Layered Medium via Wave Localization,” Phys. Rev. Lett. 94(7), 073901 (2005). [CrossRef]
27. J.-K. Yang, H. Noh, M. J. Rooks, G. S. Solomon, F. Vollmer, and H. Cao, “Lasing in localized modes of a slow light photonic crystal waveguide,” Appl. Phys. Lett. 98(24), 241107 (2011). [CrossRef]
28. B. Kumar, R. Homri, Priyanka, S. K. Maurya, M. Lebental, and P. Sebbah, “Localized modes revealed in random lasers,” Optica 8(8), 1033–1039 (2021). [CrossRef]
29. L. Sapienza, H. Thyrrestrup, S. Stobbe, P. D. Garcia, S. Smolka, and P. Lodahl, “Cavity quantum electrodynamics with anderson-localized modes,” Science 327(5971), 1352–1355 (2010). [CrossRef]
30. J. Gao, S. Combrie, B. Liang, P. Schmitteckert, G. Lehoucq, S. Xavier, X. Xu, K. Busch, D. L. Huffaker, A. De Rossi, and C. W. Wong, “Strongly coupled slow-light polaritons in one-dimensional disordered localized states,” Sci. Rep. 3(1), 1994 (2013). [CrossRef]
31. L. Labonté, C. Vanneste, and P. Sebbah, “Localized mode hybridization by fine tuning of two-dimensional random media,” Opt. Lett. 37(11), 1946–1948 (2012). [CrossRef]
32. J. Bertolotti, S. Gottardo, D. S. Wiersma, M. Ghulinyan, and L. Pavesi, “Optical Necklace States in Anderson Localized 1D Systems,” Phys. Rev. Lett. 94(11), 113903 (2005). [CrossRef]
33. P. Sebbah, B. Hu, J. M. Klosner, and A. Z. Genack, “Extended Quasimodes within Nominally Localized Random Waveguides,” Phys. Rev. Lett. 96(18), 183902 (2006). [CrossRef]
34. J. B. Pendry, “Quasi-extended electron states in strongly disordered systems,” J. Phys. C: Solid State Phys. 20(5), 733–742 (1987). [CrossRef]
35. J. B. Pendry, “Symmetry and transport of waves in one-dimensional disordered systems,” Adv. Phys. 43(4), 461–542 (1994). [CrossRef]
36. C. Vanneste and P. Sebbah, “Complexity of two-dimensional quasimodes at the transition from weak scattering to Anderson localization,” Phys. Rev. A 79(4), 041802 (2009). [CrossRef]
37. A. F. Tzortzakakis, K. G. Makris, and E. N. Economou, “Non-Hermitian disorder in two-dimensional optical lattices,” Phys. Rev. B 101(1), 014202 (2020). [CrossRef]
38. S. Kalish, Z. Lin, and T. Kottos, “Light transport in random media with PT symmetry,” Phys. Rev. A: At., Mol., Opt. Phys. 85(5), 055802 (2012). [CrossRef]
39. O. Bendix, R. Fleischmann, T. Kottos, and B. Shapiro, “Exponentially Fragile $\mathcal {P}\mathcal {T}$ Symmetry in Lattices with Localized Eigenmodes,” Phys. Rev. Lett. 103(3), 030402 (2009). [CrossRef]
40. H. Vemuri, V. Vavilala, T. Bhamidipati, and Y. N. Joglekar, “Dynamics, disorder effects, and $\mathcal {PT}$-symmetry breaking in waveguide lattices with localized eigenstates,” Phys. Rev. A 84(4), 043826 (2011). [CrossRef]
41. O. Vázquez-Candanedo, J. C. Hernández-Herrejón, F. M. Izrailev, and D. N. Christodoulides, “Gain- or loss-induced localization in one-dimensional $\mathcal {PT}$-symmetric tight-binding models,” Phys. Rev. A 89(1), 013832 (2014). [CrossRef]
42. M. J. Hartmann, F. G. S. L. Brandão, and M. B. Plenio, “Strongly interacting polaritons in coupled arrays of cavities,” Nat. Phys. 2(12), 849–855 (2006). [CrossRef]
43. A. D. Greentree, C. Tahan, J. H. Cole, and L. C. L. Hollenberg, “Quantum phase transitions of light,” Nat. Phys. 2(12), 856–861 (2006). [CrossRef]
44. P. T. Leung, S. Y. Liu, S. S. Tong, and K. Young, “Time-independent perturbation theory for quasinormal modes in leaky optical cavities,” Phys. Rev. A 49(4), 3068–3073 (1994). [CrossRef]
45. P. T. Leung, S. Y. Liu, and K. Young, “Completeness and orthogonality of quasinormal modes in leaky optical cavities,” Phys. Rev. A 49(4), 3057–3067 (1994). [CrossRef]
46. P. T. Leung, W. M. Suen, C. P. Sun, and K. Young, “Waves in open systems via a biorthogonal basis,” Phys. Rev. E 57(5), 6101–6104 (1998). [CrossRef]
47. P. M. Morse and H. Feshbach, Methods of theoretical physics (McGraw-Hill, New York, 1953).
48. N. Moiseyev, P. R. Certain, and F. Weinhold, “Resonance properties of complex-rotated hamiltonians,” Mol. Phys. 36(6), 1613–1630 (1978). [CrossRef]
49. D. Laurent, O. Legrand, P. Sebbah, C. Vanneste, and F. Mortessagne, “Localized modes in a finite-size open disordered microwave cavity,” Phys. Rev. Lett. 99(25), 253902 (2007). [CrossRef]
50. R. Touzani, “OFELI, An Object Oriented Finite Element Library,” http://ofeli.org/ (2013).
51. V. Hernandez, J. E. Roman, and V. Vidal, “SLEPc: A Scalable and Flexible Toolkit for the Solution of Eigenvalue Problems,” ACM Transactions on Mathematical Software 31(3), 351–362 (2005). [CrossRef]
52. M. B. Doost, W. Langbein, and E. A. Muljarov, “Resonant state expansion applied to two-dimensional open optical systems,” Phys. Rev. A 87(4), 043827 (2013). [CrossRef]
53. H. Hodaei, A. U. Hassan, S. Wittek, H. Garcia-Gracia, R. El-Ganainy, D. N. Christodoulides, and M. Khajavikhan, “Enhanced sensitivity at higher-order exceptional points,” Nature 548(7666), 187–191 (2017). [CrossRef]
54. F. Riboli, P. Barthelemy, S. Vignolini, F. Intonti, A. De Rossi, S. Combrie, and D. S. Wiersma, “Anderson localization of near-visible light in two dimensions,” Opt. Lett. 36(2), 127–129 (2011). [CrossRef]
55. F. Riboli, N. Caselli, S. Vignolini, F. Intonti, K. Vynck, P. Barthelemy, A. Gerardino, L. Balet, L. H. Li, A. Fiore, M. Gurioli, and D. S. Wiersma, “Engineering of light confinement in strongly scattering disordered media,” Nat. Mater. 13(7), 720–725 (2014). [CrossRef]
56. R. Espinola, M. Tsai, J. T. Yardley, and R. Osgood, “Fast and low-power thermooptic switch on thin silicon-on-insulator,” IEEE Photonics Technol. Lett. 15(10), 1366–1368 (2003). [CrossRef]
57. C. Ropp, N. Bachelard, D. Barth, Y. Wang, and X. Zhang, “Dissipative self-organization in optical space,” Nat. Photonics 12(12), 739–743 (2018). [CrossRef]
58. S. Schönhuber, N. Bachelard, B. Limbacher, M. A. Kainz, A. M. Andrews, H. Detz, G. Strasser, J. Darmo, S. Rotter, and K. Unterrainer, “All-optical adaptive control of quantum cascade random lasers,” Nat. Commun. 11(1), 5530 (2020). [CrossRef]
59. K. Y. Bliokh, Y. P. Bliokh, V. Freilikher, A. Genack, and P. Sebbah, “Coupling and level repulsion in the localized regime: From isolated to quasiextended modes,” Phys. Rev. Lett. 101(13), 133901 (2008). [CrossRef]
60. A. Majumdar, A. Rundquist, M. Bajcsy, V. D. Dasika, S. R. Bank, and J. Vučković, “Design and analysis of photonic crystal coupled cavity arrays for quantum simulation,” Phys. Rev. B 86(19), 195312 (2012). [CrossRef]