• Photonics Research
  • Vol. 10, Issue 4, 1134 (2022)
Stefano Biasi†、*, Riccardo Franchi1、†, Filippo Mione, and Lorenzo Pavesi
Author Affiliations
  • Nanoscience Laboratory, Dipartimento di Fisica, Università di Trento, 38123 Trento, Italy
  • show less
    DOI: 10.1364/PRJ.450402 Cite this Article Set citation alerts
    Stefano Biasi, Riccardo Franchi, Filippo Mione, Lorenzo Pavesi, "Interferometric method to estimate the eigenvalues of a non-Hermitian two-level optical system," Photonics Res. 10, 1134 (2022) Copy Citation Text show less

    Abstract

    Non-Hermitian physics has found a fertile ground in optics. Recently, the study of mode coalescence, i.e., exceptional points, has led to the discovery of intriguing and counterintuitive phenomena. Degeneracies are typically modeled through the coupled mode theory to determine the behavior of eigenstates and eigenvalues. However, the complex nature of the eigenvalues makes their characterization from the response spectrum difficult. Here, we demonstrate that a coherent interferometric excitation allows estimation of both the real and imaginary parts of the eigenvalues. We study the clockwise and counter-clockwise modes in optical microresonators both in the case of Hermitian and non-Hermitian intermodal coupling. We show the conditions by which a resonant doublet, due to the dissipative coupling of counter-propagating modes caused by surface roughness backscattering, merges to a single Lorentzian. This permits us to estimate the optimal quality factor of the microresonator in the absence of modal coupling caused by backscattering. Furthermore, we demonstrate that a taiji microresonator working at an exceptional point shows a degeneracy splitting only in one propagation direction and not in the other. This follows from the strongly non-Hermitian intermodal coupling caused by the inner S-shaped waveguide.

    1. INTRODUCTION

    In an open system, energy can be exchanged with the environment. Its description by a Hermitian Hamiltonian is approximated, because the conservation of energy is not a priori satisfied [1]. In this case, a non-Hermitian Hamiltonian has to be used, which leads to many fascinating new phenomena [2]. A radical difference between non-Hermitian (dissipative) and Hermitian (conservative) physics arises in the presence of degeneracies, i.e., when the Hamiltonian eigenvalues coalesce [3]. In the conservative case, these degeneracies are called diabolic points, while in the dissipative one they are called exceptional points. The former are characterized by real eigenvalues that coalesce while maintaining orthogonal eigenvectors [4]. In contrast, the latter are characterized by complex eigenvalues that coalesce simultaneously with the eigenvectors [2]. Close to these exceptional points, counterintuitive and intriguing phenomena occur [3,5]. Since the propagation of light is generally affected by losses, optics constitutes one of the most fertile grounds to study open systems [6].

    In recent years, the realization of optical structures working near exceptional points has allowed the demonstration of non-trivial effects such as unidirectional invisibility [79], loss-induced transparency [10], directional emission [1113], and the PT-symmetric laser [14]. Typically, these optical systems are based on microring resonators, which are described within the temporal coupled mode theory [11,15]. In these systems, knowledge of the eigenvalues represents a fundamental aspect. However, their complex nature makes difficult a simple and direct estimation from their optical response spectra.

    In this work, we describe a method to estimate non-Hermitian eigenvalues based on a coherent interferometric technique to excite the modes of an optical microresonator. The coherent interferometric excitation technique is based on the simultaneous excitation of the clockwise and counter-clockwise modes of a microresonator coupled to a bus waveguide. Its key aspect is that the microresonator response depends on the relative phase difference of the signals input at the two sides of the bus waveguide (Fig. 1). We apply the coherent interferometric excitation to representative Hermitian and non-Hermitian systems which are based on a simple microring or taiji microresonator [8]. In the former case, the coupling between the counter-propagating modes is given by the surface-wall roughness, i.e., the backscattering [16], while in the second, it is related to an S-shaped waveguide embedded in the taiji microresonator [1719].

    Sketch of (a) ring and (b) taiji microresonator coupled to a bus waveguide. The arrows indicate the propagating fields. Specifically, the red (blue) highlights propagation in the clockwise (counter-clockwise) direction. The different symbols are defined in the text.

    Figure 1.Sketch of (a) ring and (b) taiji microresonator coupled to a bus waveguide. The arrows indicate the propagating fields. Specifically, the red (blue) highlights propagation in the clockwise (counter-clockwise) direction. The different symbols are defined in the text.

    In Section 2, we model the bus waveguide/microresonator system response as a function of the relative phase between the signals input in the forward and reverse directions. We show the merging of the resonant doublet, due to the microresonator modes, into a single Lorentzian. We propose an analytical expression for the eigenstates and eigenvectors. We demonstrate that, in both the Hermitian and non-Hermitian cases, the coherent interferometric excitation can be used to determine the microresonator Q-factor in absence of the surface roughness losses. For Hermitian systems, we show that the coupling coefficients can be estimated from the spectral position of the single Lorentzian resonances. For non-Hermitian systems, we demonstrate that the complex nature of the eigenvalues and their characteristic parameters can be studied by changing the relative intensities and phases of the input signals in the coherent interferometric excitation. In this framework, we show that the taiji microresonator works on an exceptional point and we show a direction-dependent response due to the coherent interferometric excitation. A resonance splitting is predicted only for transmission in the propagation direction, which couples to the reflection of the inner S-shaped waveguide.

    In Section 3, we implement the coherent interferometric technique in a simple setup and study two different types of microresonators where non-Hermitian physics can be demonstrated. We extract the different parameters (intrinsic and extrinsic loss coefficients, real and imaginary intermodal coupling coefficients) for a simple microring resonator and a taiji microresonator, and we validate the model by fitting their transmission spectra.

    2. THEORETICAL MODEL

    Let us consider the electric field amplitudes of the clockwise (αCW) and counter-clockwise (αCCW) propagating modes in a simple or taiji microresonator coupled to a bus waveguide [Figs. 1(a) and 1(b)]. We model these systems with the temporal coupled mode theory and the properties of a time-reversal-invariant system [1921]. Therefore, αCW and αCCW are solutions of [22]iddt(αCCWαCW)=(ω0i(γ+Γ)iβ12iβ21ω0i(γ+Γ))(αCCWαCW)2Γ(Ein,LEin,R),where ω0 is the angular frequency of the microresonator resonance, Ein,L (Ein,R) is the electric field amplitude input into the left (right) edge of the bus waveguide, γ and Γ are the intrinsic and extrinsic damping rates, and β12 and β21 are the intermodal coupling coefficients. Precisely, γ describes the losses of the cavity related to intrinsic factors, such as material absorption and bending, while Γ refers to the coupling with the bus waveguide. It is worth noting that γ and Γ are real numbers while β12 and β21 are complex numbers. These latter describe the energy exchange between the counter-propagating modes, which can have either Hermitian (conservative) or non-Hermitian (dissipative) nature. In a simple microring, the intermodal energy exchange is caused by the backscattering, which yields complex coupling coefficients with independent real and imaginary parts [22]. In a taiji microresonator, neglecting the backscattering, β21 reduces to zero while β12 is defined by the inner S-shaped waveguide [8,17,23].

    In the interferometric excitation technique, electric fields Ein,L:=εin,Leiωt and Ein,R:=εin,Reiϕeiωt are coherently input from both sides (left L, right R) of the bus waveguide and the transmitted fields Eout,L and Eout,R are measured to probe the system response. Note the phase difference ϕ between Ein,R and Ein,L. The input and output fields are related by [20,22]:(Eout,REout,L)=(Ein,LEin,R)+i2Γ(αCCWαCW),where in the steady state αCW/CCW:=aCW/CCWeiωt. Equations (1) and (2) lead to the following steady-state output field amplitudes:εout,R=(12Γ(iΔω+γ+Γ)(iΔω+γ+Γ)2β12β21)εin,L+(2Γβ12(iΔω+γ+Γ)2β12β21)eiϕεin,R,εout,L=(12Γ(iΔω+γ+Γ)(iΔω+γ+Γ)2β12β21)eiϕεin,R+(2Γβ21(iΔω+γ+Γ)2β12β21)εin,L,where Δω:=ωω0. Equations (3) and (4) show that the output fields result from the sum of the transmitted and reflected input fields from the microresonator. Note the relevant role played by the phase term ϕ, which allows tuning of the excitation of the two microresonator modes. This role is better exemplified by introducing the super-modes: b1,2:=(αCCW±αCWeiϕ)/2. Here, the sum and difference of the counter-propagating amplitudes are governed by ϕ. In stationary condition, Eq. (1) can be expressed in the super-mode description as ω(b1b2)=(ω0i(γ+Γ+g1)ig2ig2ω0i(γ+Γg1))(b1b2)2Γ(εin,L+εin,R2εin,Lεin,R2),where g1,2:=(β12eiϕ±β21eiϕ)/2. Within the super-mode description, the microresonator excitation results from a linear combination of only the left and right input field amplitudes.

    On the contrary, the coupling terms (g1,2) depend on ϕ, on β12, and on β21. Equation (5) has eigenvaluesλ1,2=ω0±β12β21i(γ+Γ),which do not depent on ϕ since they are an intrinsic property of the system. In contrast, the normalized eigenstates of Eq. (5)v1,2=11+|g1g12g22g2|2(g1g12g22g21)are affected by the interferometric excitation, i.e., ϕ.

    It is customary to define the Hermitian (h) and non-Hermitian (n) intermodal coupling coefficients [12,22]h:=iβ12β21*2,n:=β12+β21*2.

    These definitions allow us to discuss in the following Hermitian and non-Hermitian systems, separately.

    A. Hermitian Coupling and Diabolic Point

    Here, β12=β21*=β, and therefore h0 and n=0. This generates a continuous conservative exchange of energy between the clockwise and the counter-clockwise modes inducing the typical splitting (doublet) in the transmission spectra of a bus waveguide/microring system [22]. Figures 2(a1) and 2(a2) show the computed output spectral intensities for a real coupling coefficient β=33.2  GHz when εin,L=1mW·Hz and εin,R=0 (top spectra) or when εin,L=0 and εin,R=1mW·Hz (bottom spectra), i.e., in a standard, single excitation side experiment. Hence, in the top graph, |εout,R|2 and |εout,L|2 are the transmitted and reflected intensities while the reverse holds for the bottom graph. The output fields are not simple Lorentzians, but show a balanced doublet (same peak intensities) due to the interaction between the counter-propagating modes in the microring [24]. Lorentz reciprocity theorem [25,26] assures equal transmission spectra and the Hermitian coupling yields equal reflection spectra for both input (L or R) directions.

    Theoretical results for a Hermitian coupling. Panels (a1) and (a2) show the transmission and reflection spectra for an excitation from the left and right input, respectively. Maps (b1) and (b2) show, respectively, the right and left output field intensities as a function of the angle (θh) and of the detuning frequency (Δω) for a symmetric interferometric excitation in units of mW·Hz. Panels (c1) and (c2) show the output field spectral lineshapes for fixed values of θh={±π,π/4,π/2,−π/2}. The red (blue) lines denote the left (right) output fields. The graph (d) shows the modulus of the inner product of the eigenstates and their real and imaginary parts as a function of θh. The solid-black, solid-orange, and dashed-magenta lines refer to |⟨v1|v2⟩|, |⟨ℑ[v1]|ℑ[v2]⟩|, and |⟨ℜ[v1]|ℜ[v2]⟩|, respectively. The coupling is Hermitian, and therefore the two eigenstates are orthogonal, i.e., ⟨v1|v2⟩=0. Here, we used the following coefficients: Γ=γ=6.8 GHz, β12=−β21*=33.2 GHz, and ω0=2π·193 THz.

    Figure 2.Theoretical results for a Hermitian coupling. Panels (a1) and (a2) show the transmission and reflection spectra for an excitation from the left and right input, respectively. Maps (b1) and (b2) show, respectively, the right and left output field intensities as a function of the angle (θh) and of the detuning frequency (Δω) for a symmetric interferometric excitation in units of mW·Hz. Panels (c1) and (c2) show the output field spectral lineshapes for fixed values of θh={±π,π/4,π/2,π/2}. The red (blue) lines denote the left (right) output fields. The graph (d) shows the modulus of the inner product of the eigenstates and their real and imaginary parts as a function of θh. The solid-black, solid-orange, and dashed-magenta lines refer to |v1|v2|, |[v1]|[v2]|, and |[v1]|[v2]|, respectively. The coupling is Hermitian, and therefore the two eigenstates are orthogonal, i.e., v1|v2=0. Here, we used the following coefficients: Γ=γ=6.8  GHz, β12=β21*=33.2  GHz, and ω0=2π·193  THz.

    The situation changes when the interferometric excitation is used. In the Hermitian case, as shown in Appendix A, Eq. (7) can be reformulated asv1,2=11+|sin[θh]1cos[θh]|2(isin[θh]1cos[θh]1),where the angle θh=ϕ+arg[β] is introduced to show that the eigenstates do not depend on the strength of the coupling (|β|) but only on its phase. As expected for a Hermitian system, its eigenvectors are orthogonal and, therefore, their inner product |v1|v2|=0 for any ϕ. On the other hand, its eigenvalues depend on |β|:λ1,2=ω0±|β|i(γ+Γ),but not on ϕ. Their imaginary parts (I[λ1,2]=γΓ) are related to the losses, while their real parts (R[λ1,2]=ω1,2=ω0±|β|) describe the splitting of the resonance as a result of the coupling between the counter-propagating (back-reflected) modes. When |β|=0, the system works on a diabolic point [3]. Indeed, the eigenvalues are degenerate (ω1=ω2=ω0), the eigenvectors are orthogonal, and the resonance splitting is linear with |β|.

    Using a symmetric (i.e., same field intensities) interferometric excitation (εin,L=εin,R=1mW·Hz), εout,L and εout,R depend on θh. Figures 2(b1) and 2(b2) show the transmitted field intensities as a function of Δω and θh. A rich dynamics is observed with transmitted intensities which increase or decrease depending on the excitation conditions. Note that since the sum of the input intensities is equal to 2 mW·Hz, it is not surprising to find transmitted intensities larger than 1 mW·Hz. Figures 2(c1)–2(c4) show the spectra of the transmitted intensities for θh={±π,π/4,π/2,π/2}. As a function of θh, different transmission lineshapes at the bus waveguide outputs are observed.

    Figure 2(d) shows the modulus of the inner product of the eigenstates (|v1|v2|) and their real (|R[v1]|R[v2]|) and imaginary (|I[v1]|I[v2]|) parts as a function of θh. Specifically, the inner product of the real and imaginary components follows a periodic function proportional to |cos[θh]|/2, changing their mutual projection as θh varies. When the real and imaginary parts of the eigenstates are parallel, i.e., θh={π,0,π}, the transmission doublet is balanced [see Fig. 2(c1)], and there is a symmetric energy exchange between the counter-propagating modes. On the contrary, when θh=±π/2, the transmission doublet merges into a single resonance and the same spectra are observed from both waveguide ends. In these peculiar cases, both the real and imaginary parts of the eigenstates become orthogonal and their inner products reduce to zero. For these θh values, both the left and right output intensities take the shape of a single Lorentzian.

    Under these conditions, i.e., for θh=ϕ+arg[β]=±π/2+2πm (mZ) and ε0=εin,R=εin,L, Eqs. (3) and (4) reduce to the typical analytical expression of a Lorentzian [27,28]:|εout,R|2=|εout,L|2=(14γΓ(Δω|β|)2+(γ+Γ)2)|ε0|2.

    Let us discuss this result. Equation (11) shows that the transmission of the bus waveguide/microresonator system is the same as the one obtained in the absence of the intermodal coupling (i.e., no backscattering). The position of the Lorentzian is spectrally shifted by the intermodal coupling strength with respect to the microresonator resonance. Thus, the interferometric excitation permits us to extract relevant parameters directly from the measured transmission spectra: The amplitude |β|, by measuring the resonant frequencies of the Lorentzians;The phase arg[β], by measuring the input fields phase difference (ϕ) at the Lorentzians;The intrinsic and extrinsic damping rates γ+Γ, by measuring the full width at half maximum (FWHM) of the Lorentzians.

    These parameters are not easily estimated from the usual (one-side input) spectral response measurements because, the inter-mode interaction strongly modifies the energy stored inside the cavity and, furthermore, the minima of the transmission doublet do not directly relate with the real parts of the eigenvalues.

    In this subsection, we treated the case of a symmetric interferometric excitation. An asymmetric interferometric excitation induces an asymmetric system response, giving rise to a controlled unbalance of the resonant doublet. It is worth noting that at the diabolic point, i.e., |β|=0, the response of the bus waveguide/microresonator system is insensitive to any phase variation in the interferometric excitation. In fact, the counter-propagating modes do not interfere, and the output fields show a Lorentzian lineshape.

    B. Non-Hermitian Coupling and Exceptional Point

    In this case, there is no relation between β12 and β21. Both h and n can be different from zero. Since we are dealing with passive microring, we restrict ourselves to the no-gain case where |n|γ [22]. As a result, n induces a dissipative term which causes an asymmetric energy exchange between the counter-propagating modes. In the case of a single-side input, an unbalanced (i.e., different intensities) transmission doublet is observed. This is shown in Figs. 3(a1) and 3(a2). Again, the Lorentz reciprocity theorem ensures equal transmissions [25,26]. On the other hand, the non-Hermitian intermodal coupling induces different reflections for the different input directions [2].

    In the case of an interferometric excitation, the non-Hermitian intermodal coupling makes the eigenstates of Eq. (7) non-orthogonal (v1|v20). By writing β12/21=|β12/21|eiφ12/21, we can write the eigenvalues as λ1,2=ω0|β12||β21|sin[φ12+φ212]i(Γ+γ|β12||β21|cos[φ12+φ212]).Here, the real part is related to the splitting of the transmission response. In contrast, the imaginary part induces an unbalance in the energy exchange between the counter-propagating modes. In analogy to the Hermitian case, we define the angle θ=ϕφ21φ122π2. This allows for establishing a uniform notation and, especially, treating several cases of non-Hermitian coupling.

    Theoretical results for a non-Hermitian coupling. The transmitted and reflected intensities are shown in panels (a1) and (a2), for left (εin,L=1mW·Hz and εin,R=0) or right (εin,L=0 and εin,R=1mW·Hz) input. The blue (red) curves denote the right (left) output fields. Panels (b1) and (b2) show the intensity of the output fields as a function of the angle θ and frequency detuning Δω for a symmetric interferometric excitation (εin,L=εin,R=1mW·Hz) in units of mW·Hz. (c1)–(c4) Spectral output field intensities for few θ values given in each panel. (d) Modulus of the eigenstate inner product. The non-Hermitian nature of the coupling makes the eigenstates non-orthogonal, i.e., ⟨v1|v2⟩≠0. Here, we used the following coefficients: Γ=γ=6.8 GHz, β12=20.2 GHz, β21=(−20.2+9i) GHz, and ω0=2π·193 THz.

    Figure 3.Theoretical results for a non-Hermitian coupling. The transmitted and reflected intensities are shown in panels (a1) and (a2), for left (εin,L=1mW·Hz and εin,R=0) or right (εin,L=0 and εin,R=1mW·Hz) input. The blue (red) curves denote the right (left) output fields. Panels (b1) and (b2) show the intensity of the output fields as a function of the angle θ and frequency detuning Δω for a symmetric interferometric excitation (εin,L=εin,R=1mW·Hz) in units of mW·Hz. (c1)–(c4) Spectral output field intensities for few θ values given in each panel. (d) Modulus of the eigenstate inner product. The non-Hermitian nature of the coupling makes the eigenstates non-orthogonal, i.e., v1|v20. Here, we used the following coefficients: Γ=γ=6.8  GHz, β12=20.2  GHz, β21=(20.2+9i)  GHz, and ω0=2π·193  THz.

    For a symmetric interferometric excitation, we obtain the maps shown in Figs. 3(b1) and 3(b2). Changing ϕ, i.e., θ, a redistribution of energy between the counter-propagating modes is observed, which results in different doublet shapes of the transmitted output fields. For θ=±π, the dissipative term induces an unbalanced strongly asymmetric doublet [see Fig. 3(c1)]. For θ=±π/2, i.e., Figs. 3(c3) and 3(c4), the doublet merges in a single asymmetric peak, which does not have a simple Lorentzian shape as in the Hermitian case. Moreover, |v1|v2| is different from zero for each θ, while its real and imaginary parts follow a trend similar to the Hermitian case [Fig. 3(d)]. However, |R[v1]|R[v2]| does not decrease to zero at ±π/2, but shows a minimum equal to |v1|v2|. In contrast, |I[v1]|I[v2]| goes to zero as in the Hermitian case.

    The lack of coalescence of the doublet into a simple Lorentzian, when both the real and imaginary parts of the eigenstates reach their inner product minima, is related to the dissipative part of the coupling (n). To get a Lorentzian, we can use an asymmetric interferometric excitation. Indeed, whenεin,R=|β21||β12|εin,Landθ=±π2+2πm,where mZ, the transmission lineshapes reduce to a single Lorentzian lineshape (Fig. 4). The first condition in Eq. (13) compensates for the different strengths of the coupling coefficients, while the second ensures the same interference between the transmitted and reflected fields for both directions. As a result, Eq. (13) establishes the proportionality between the output field intensities. Substituting Eq. (13) in the Eq. (3), we obtain(|εout,R|2|εout,L|2)=(14(γγ˜)Γ(Δωβ˜)2+(γγ˜+Γ)2)(|εin,L|2|εin,R|2),β˜:=|β12||β21|sin[(φ12+φ21)/2],γ˜:=|β12||β21|cos[(φ12+φ21)/2].β˜ and γ˜ are the real and imaginary parts of the shifts of the eigenvalues with respect to ω0 [see Eq. (12)].

    Theoretical results for an asymmetric interferometric excitation which satisfies Eq. (13). Panels (b1) and (b2) show the intensities of the output fields as a function of the angle θ and frequency detuning Δω in units of mW·Hz. The graphs (b1) and (b2) show the spectral output field intensities for θ=π/2 and θ=−π/2, respectively. Here, we used the same coefficients as in Fig. 3, i.e., Γ=γ=6.8 GHz, β12=20.2 GHz, β21=(−20.2+9i) GHz, and ω0=2π·193 THz.

    Figure 4.Theoretical results for an asymmetric interferometric excitation which satisfies Eq. (13). Panels (b1) and (b2) show the intensities of the output fields as a function of the angle θ and frequency detuning Δω in units of mW·Hz. The graphs (b1) and (b2) show the spectral output field intensities for θ=π/2 and θ=π/2, respectively. Here, we used the same coefficients as in Fig. 3, i.e., Γ=γ=6.8  GHz, β12=20.2  GHz, β21=(20.2+9i)  GHz, and ω0=2π·193  THz.

    Figures 4(a1) and 4(a2) show the maps of the transmitted field intensities as a function of Δω and θ for an asymmetric interferometric excitation which satisfies Eq. (13). For θ=±π/2 [Figs. 4(b1) and 4(b2)], the transmitted fields show a simple Lorentzian lineshape, and most importantly, their minima indicate the real part of the eigenvalues. Furthermore, the FWHMs of the Lorentzian lineshapes are related to the imaginary part of the eigenvalues and depend on the intrinsic microresonator losses. Thus, the interferometric excitation allows for the study of both the real and imaginary parts of the eigenvalues. In Fig. 4(b1), the best quality factor of the microresonator in the absence of backscattering corresponds to the value obtained by setting the relative phase ϕ so that θ=π/2.

    When just one of the two coupling coefficients β12/21 tends to zero, the eigenvalues [Eq. (12)] vary following a square root dependence. More importantly, when one of the two coefficients is zero, the eigenvalues become degenerate and the eigenstates [Eq. (7)] become parallel (v1|v2=1). Therefore, the system works on an exceptional point [3]. This situation might be found in an ideal taiji microresonator [Fig. 1(b)]. Neglecting the backscattering, this microresonator has β120 and β21=0. In this case, the energy exchange between the counter-propagating modes is strongly unbalanced in one direction. As can be seen in Figs. 5(a1) and 5(a2), the left- or right-side transmission shows the simple Lorentzian lineshape when a single input excitation is used. On the other hand, the two reflections are completely different due to the non-Hermitian intermodal coupling. The intensity reflected from the right waveguide edge shows a peak, while that from the left waveguide edge is zero. The response of the bus waveguide/taiji microresonator system for a symmetric interferometric excitation is shown in Figs. 5(b1) and 5(b2). It is observed that, unlike the diabolic point case, the exceptional point case is affected by an interferometric excitation. In particular, the left output field is unaffected by a ϕ change and maintains the simple Lorentzian lineshape while the right output field has a lineshape that depends markedly on ϕ [Figs. 5(c1)–5(c4)]. For ϕ=0, the right output field intensity exhibits a balanced doublet [Fig. 5(c4)]. Figure 5(d) shows the inner products of the eigenstates. As expected for a microresonator at the exceptional point, they do not depend on ϕ. In particular, the inner products of the eigenstates and their real parts are equal to 1, while they are 0 for the imaginary ones.

    Theoretical results for a taiji microresonator coupled to a bus waveguide. Panels (a1) and (a2) show the transmitted and reflected intensities as a function of the frequency detuning Δω at the left and right outputs, respectively. The blue (red) curves refer to the output fields from the right (left) waveguide edges. The maps (b1) and (b2) show the intensities of the output fields as a function of the phase ϕ and Δω for a symmetric interferometric excitation in units of mW·Hz. The transmitted spectral intensities are plotted on the graphs (c1)–(c4) for ϕ=−π/2, 3π/4, ±π, and 0, respectively. The graph (d) shows the modulus of the inner product of the eigenstates. As expected, the taiji works on an exceptional point, and therefore the eigenstates are parallel, i.e., ⟨v1|v2⟩=1. Here, we used the following coefficients: Γ=γ=6.8 GHz, β12=12 GHz<2γ, β21=0, and ω0=2π·193 THz.

    Figure 5.Theoretical results for a taiji microresonator coupled to a bus waveguide. Panels (a1) and (a2) show the transmitted and reflected intensities as a function of the frequency detuning Δω at the left and right outputs, respectively. The blue (red) curves refer to the output fields from the right (left) waveguide edges. The maps (b1) and (b2) show the intensities of the output fields as a function of the phase ϕ and Δω for a symmetric interferometric excitation in units of mW·Hz. The transmitted spectral intensities are plotted on the graphs (c1)–(c4) for ϕ=π/2, 3π/4, ±π, and 0, respectively. The graph (d) shows the modulus of the inner product of the eigenstates. As expected, the taiji works on an exceptional point, and therefore the eigenstates are parallel, i.e., v1|v2=1. Here, we used the following coefficients: Γ=γ=6.8  GHz, β12=12  GHz<2γ, β21=0, and ω0=2π·193  THz.

    As a last comment, Eq. (12) establishes that it is possible to reach the exceptional point in several ways. In fact, by vanishing to 0 one of the two intermodal coupling coefficients, e.g., β210 while β120, one can have different situations depending on the linear combination of their relative phases (φ12+φ212). Precisely, when φ12+φ212=π/2+mπ with mZ, the degeneracy is achieved as a result of the change of just the real (resonant) part. Conversely, when φ12+φ212=mπ with mZ, the two eigenvalues coalesce as a consequence of the change of just the imaginary (absorptive) part, whereas, for any other phases, an intermediate situation is obtained, and the degeneracy is achieved in the complex plane of the eigenvalues. Thus, the phase relation between β12 and β21 allows us to study an exceptional point which can be resonant, absorptive or both, in the same way as in the recent work [15]. From Eq. (12), it follows that it is convenient to use the phase relation that makes the perturbation on the eigenvalues purely real in the case of sensing through passive microresonators capable of working on an exceptional point. This permits us to estimate the strength of the perturbation by simply measuring the splitting in the transmission spectra [29,30].

    3. EXPERIMENTAL MEASUREMENTS

    A. Experimental Setup and Samples

    Our setup is sketched in Fig. 6. A fiber-coupled continuous-wave tunable laser (Yenista OPTICS, TUNICS-T100S) is used as the source. It operates at 2 mW, spanning a range of wavelength from 1470 to 1580 nm. Its emission passes through an optical isolator to avoid any spurious back-reflection, and then passes through a 50:50 fiber splitter. This allows for coherent interferometric excitation. After the fiber splitter, the two coherent signals in the different fibers are modulated in intensities by a variable optical attenuator, and are adjusted in polarization by a polarization control stage. We use the transverse electric (TE) polarization in the waveguide. In order to balance the two optical paths, the signal of one fiber is sent to a free space delay line. Two fiber taps of 10% connected to two photodiodes (Thorlabs, PDA10CS2) are used to monitor the input signals. After these, two optical circulators route these signals into two single-mode stripped fibers that couple the light into the sample. As is clear in the figure, the sample is excited from both the left and right sides by two coherent light signals. The output fields are collected by the same stripped fibers and are routed by the circulator into two other photodiode detectors (Thorlabs, PDA10CS2), labeled PDL and PDR in Fig. 6. Finally, an oscilloscope (PicoScope 4000 Series) records simultaneously the signals of all four detectors. The sample is mounted on a thermostat holder whose temperature is controlled by a proportional-integral-derivative controller (SIM960 Analog PID controller) connected to a Peltier cell and a 10 kΩ thermistor. The same setup is also used for standard (single-side excitation) measurements. In these, the signal is input from the left or right only.

    Sketch of the experimental setup. PC, personal computer to control the instruments; CWTL, continuous wave tunable laser to excite the system; OI, optical isolator to protect the laser from back-reflections; FS, fiber splitter to split the signal; VOA, variable optical attenuator to control the amplitudes; FPC, fiber polarization controller to set the transverse electric polarization at the chip grating couplers; DL, delay line to balance the two optical paths; OC, optical circulator to measure and excite simultaneously the counter-propagating fields; PD, photo detector, and PicoScope, PC oscilloscope, to measure the output powers; AS, alignment stage, and IRC, infrared camera, to see and align the stripped fibers with the sample.

    Figure 6.Sketch of the experimental setup. PC, personal computer to control the instruments; CWTL, continuous wave tunable laser to excite the system; OI, optical isolator to protect the laser from back-reflections; FS, fiber splitter to split the signal; VOA, variable optical attenuator to control the amplitudes; FPC, fiber polarization controller to set the transverse electric polarization at the chip grating couplers; DL, delay line to balance the two optical paths; OC, optical circulator to measure and excite simultaneously the counter-propagating fields; PD, photo detector, and PicoScope, PC oscilloscope, to measure the output powers; AS, alignment stage, and IRC, infrared camera, to see and align the stripped fibers with the sample.

    Two different microresonator geometries have been studied, which are based on silicon waveguides with a cross section of 450  nm×220  nm embedded in silica cladding. The first microresonator geometry has a ring shape with a radius of 7 μm, and is point-like coupled to two bus waveguides in the common add and drop configuration with a gap of 300 nm. These devices are described in detail in Ref. [31]. The second is a taiji microresonator with a racetrack geometry. It is obtained by coupling four Euler curves (radius 15 μm and angle 90°) to two pairs of straight sections of length 758 and 3106 nm. The perimeter of the microresonator is about 196.22 μm. The inner S-shaped waveguide is designed using four Euler curves. The central ones have a radius of 10 μm and an angle of 135°, and those at the ends have a radius of 10 μm and an angle of 75°. At both ends of the S-shaped waveguide, an inverse tapering ensures no back-reflections. The gaps between the taiji microresonator and the S-shaped waveguide are 210 nm. The gap between the taiji and the bus waveguide is 241 nm. In both types of microresonators, grating couplers (3.7  dB coupling loss) are used to input and output the light signals. The microresonators have been fabricated in different runs at the IMEC/Europractice facility within the multi-project wafer program.

    B. Method to Control the Phase Shift

    Since the interferometric excitation relies on coupling the light signal in both propagation directions, spurious reflections can induce unwanted interference fringes in the spectral measurements. In order to reduce these, the weak Fabry–Perot oscillations due to the reflections by the grating couplers [32] have been decreased by using glycerol as a matched index fluid. Furthermore, the different optical paths of the fibers after the first 50:50 beam splitter are compensated by a delay line. In this way, the free spectral range (i.e., period) of the interference fringes caused by the large interferometric setup is much wider than a single microresonator resonance doublet.

    In our setup, we do not have an absolute control of the interferometric excitation phase (ϕ). We leave ϕ to vary randomly and extract its value from a fit of the spectra. In fact, the phase shift between the two split signals is affected by fiber relaxation, air fluxes, and room temperature variations whose accurate control is difficult. However, the nearly balanced condition makes the time of these random phase variations much longer than the time of a single spectral acquisition. Specifically, the phase variation rate is measured to be dϕdt0.1  rad/s, while a typical acquisition time for a Δλ=0.3  nm wide spectrum is shorter than Δλ/dλdt<(0.3  nm)/(1nm/s)=0.3  s, where dλdt is the laser wavelength scanning rate. Hence, the phase variation during a single scan is less than 3×102  rad, and therefore, we can safely assume it negligible. A scanning rate of 1  nm/s is permitted by the use of the Picoscope, which synchronizes the laser wavelength scan with the recording of the data from the photodetectors. Furthermore, the symmetrical excitation is made possible by the calibration of all detectors, the measurement of the propagation losses along the two paths, and the estimation of the coupling losses of the gratings (see also Appendix B). The measurement of the resonance lineshapes, in the case of interferometric excitation, has been obtained by acquiring 100 output fields spectra. Each spectrum has a resolution of 0.1 pm. Furthermore, a guard time of 1 s has been set between each acquisition to ensure a sufficiently large relative phase variation.

    C. Results

    1. Microring Resonator

    The experimental results for the microring resonator are reported in Fig. 7. Figures 7(a1) and 7(a2) show the transmitted and reflected intensities for a single left and right excitation, respectively. It can be observed that the backscattering due to the surface-wall roughness gives rise to a slightly unbalanced doublet. The left (blue line in the top panel) and right (red line in the bottom panel) transmission spectra are equal as a consequence of the Lorentz reciprocity theorem. On the contrary, the left (red line) and right (blue line) reflections show similar features but different intensities. This is due to the non-Hermitian intermodal coupling, which induces a different energy exchange between the counter-propagating modes. The fit with the theoretical model [Eqs. (3) and (4) in the single-side excitation] is shown by the dashed-black lines (right output fields) and the dotted-black lines (left output fields). The fit yields the resonance frequency (ω02π·196.593  THz) and intrinsic (γ=4.56±0.01  GHz), extrinsic (Γ=2.662±0.003  GHz), and intermodal coupling [β12=(19.72±0.02)i(0.2±0.4)  GHz and β21=(20.67±0.02)+i(0.8±0.4)  GHz] coefficients. As a result, the non-Hermitian coefficient n=(0.48±0.01)i(0.5±0.3)  GHz is smaller than the Hermitian coefficient h=(0.3±0.3)i(20.20±0.01)  GHz. The effect of n0 is visible in the unbalanced doublet because of the relatively high Q-factor.

    Experimental results for a microring. Panels (a1) and (a2) show the transmitted and reflected intensities as a function of the frequency detuning (Δω) for the right and left excitations. Panels (b1) and (b2) display the doublet splitting (δω=ω</>−ω0) for a symmetric interferometric excitation as a function of the phase θ for the right or left output fields. ω</> are the measured frequency minima of the doublet. The magenta (purple) dots refer to δω=ω<−ω0 for the right (left) output fields. The black (gray) squares refer to δω=ω>−ω0 for the right (left) output fields. The dashed magenta (purple) and dash-dotted black (gray) curves show the theoretical results modeled using the fitting parameters of the one-side excitation. Panels (c1), (c2), (c3), and (c4) show the output field intensities as a function of Δω for θ=0.73π, −0.74π, 0.48π, and −0.52π, respectively. The blue (red) curves refer to the right (left) output field intensities. The dashed-black and dotted-black lines are the fit with the theoretical model for the right and left output fields, respectively. From the fit we obtain the following coefficients: Γ=2.662±0.003 GHz, γ=4.56±0.01 GHz, β12=(−19.72±0.02)−i(0.2±0.4) GHz, and β21=(20.67±0.02)+i(0.8±0.4) GHz.

    Figure 7.Experimental results for a microring. Panels (a1) and (a2) show the transmitted and reflected intensities as a function of the frequency detuning (Δω) for the right and left excitations. Panels (b1) and (b2) display the doublet splitting (δω=ω</>ω0) for a symmetric interferometric excitation as a function of the phase θ for the right or left output fields. ω</> are the measured frequency minima of the doublet. The magenta (purple) dots refer to δω=ω<ω0 for the right (left) output fields. The black (gray) squares refer to δω=ω>ω0 for the right (left) output fields. The dashed magenta (purple) and dash-dotted black (gray) curves show the theoretical results modeled using the fitting parameters of the one-side excitation. Panels (c1), (c2), (c3), and (c4) show the output field intensities as a function of Δω for θ=0.73π, 0.74π, 0.48π, and 0.52π, respectively. The blue (red) curves refer to the right (left) output field intensities. The dashed-black and dotted-black lines are the fit with the theoretical model for the right and left output fields, respectively. From the fit we obtain the following coefficients: Γ=2.662±0.003  GHz, γ=4.56±0.01  GHz, β12=(19.72±0.02)i(0.2±0.4)GHz, and β21=(20.67±0.02)+i(0.8±0.4)  GHz.

    Measurement results by the symmetric interferometric excitation are summarized in Fig. 7. Figures 7(b1) and 7(b2) show the doublet splitting as a function of ϕ. Since our setup does not allow an absolute control of ϕ, θ=ϕφ21φ122π2 was estimated by fitting the 100 experimental spectra of the interferometric excitation with the model of Eqs. (3) and (4). In this, we used the coefficients estimated for the single-side excitation as fixed parameters and the relative phase ϕ as the only free fit parameter. Typically the error on ϕ is of 3×102  rad. In the figure, the doublet splitting is reported as a frequency detuning with respect to the microring resonance δω=ω</>ω0, where ω</> are the measured frequencies of the doublet minima and </> refer to the position with respect to ω0. Figures 7(b1) and 7(b2) report δω for the right and left output fields, respectively. In the plot, we use the following notation: the magenta (purple) dots refer to δω=ω<ω0, the black (gray) squares refer to δω=ω>ω0, and the dashed magenta (purple) and dash-dotted black (gray) curves show the theoretically calculated detuning by Eqs. (3) and (4) with the parameters from the single-side excitation fits. It is worth noticing that the analytical model of Section 2 reproduces the experimental data.

    Figures 7(c1)–7(c4) show the intensities of the output fields as a function of Δω for θ=0.73π, 0.74π, 0.48π, and 0.52π. Notation is the same as for Figs. 7(a1) and 7(a2). It is clear that our theoretical model is in agreement with the experimental data. As expected, the ϕ variation induces a deformation of the doublet which, when θ±π/2, merges into a single resonance. For these θ values, the dots coincide with the square in Figs. 7(b1) and 7(b2). Specifically, Figs. 7(c3) and 7(c4) show that the output fields have a quasi-simple Lorentzian lineshape. Their resonant frequencies correspond to the real part of the eigenvalues with a negative [Fig. 7(c3)] and positive [Fig. 7(c4)] detuning.

    Assuming a pure Hermitian case, we can estimate from the spectra the coefficient parameters following the method described in Section 2. Table 1 shows the comparison between the real and imaginary parts of the eigenvalues λ1,2 extracted by a full fit of the single-side spectra [Figs. 7(a1) and 7(a2)] or estimated from the spectral positions of the collapsed doublet [Figs. 7(c3) and 7(c4)]. In particular, we fit the interferometric excitation spectra with a Lorentzian lineshape to extract the spectral positions and the FWHM of the merged doublets. Table 1 shows that the two methods yield comparable results. Note that the different values and the large error bars for the estimated I[λ1,2] are due to the presence of a non-Hermitian term n0 [see Eq. (12)]. In fact, the slight non-Hermitianity, symmetric interferometric excitation with θ±π/2 causes a weak perturbation to the Lorentzian lineshape, which is reflected in the values given in Table 1 (see also Fig. 3). From the FWHM of the Lorentzian lineshapes, it is still possible to estimate the Q-factor of the microring in absence of surface-wall roughness, which results to be 9×104.

    Real and Imaginary Parts of the Microresonator’s Eigenvaluesa

     Single-Side Excitation [GHz]Interferometric Excitation [GHz]
    |R[λ1,2ω0]|20.2±0.321.0±0.3
    I[λ1]6.73±0.026.8±0.2
    I[λ2]7.71±0.027.51±0.07

    The single-side excitation column shows the values calculated using the parameters obtained from the fit of the spectral transmission and reflection responses [Figs. 7(a1) and 7(a2)]. The column labeled interferometric excitation shows the values estimated by using the experimental data in Figs. 7(c3)–7(c4).

    2. Taiji Microresonator

    The experimental measurements for the taiji microresonator are shown in Fig. 8. As in the case of the microring, we first characterized the taiji microresonator with single-side excitation experiments. Results are reported in Figs. 8(a1) and 8(a2). In agreement with the Lorentz reciprocity theorem, the transmitted intensities are equal. In contrast, the presence of the inner S-shaped waveguide makes the reflection approximately equal to 0 for left excitation and close to 1 for right excitation. Fitting the spectra with the Eqs. (3) and (4) in the single-side excitation yields the values of the resonance frequency (ω02π·194.430  THz) and intrinsic (γ=42.6±0.1  GHz), extrinsic (Γ=75.1±0.1  GHz), and intermodal coupling [β12=(90.8±0.1)i(24.17±0.03)  GHz and β21=(8.7±0.1)i(7.30±0.09)  GHz] coefficients. The coefficient β21 is non-zero due to the presence of surface-wall roughness which causes a weak backscattering. The reflection due to backscattering is about 8.5×103times less intense than the reflection due to the inner S-shaped waveguide. Figures 8(b1) and 8(b2) show the doublet splitting for a symmetric interferometric excitation as a function of the ϕ and δω. As before, ϕ is estimated by fitting 100 spectra with only the relative phase as a free parameter, while all the other coefficients are equal to the ones estimated from the single-side experiments. Noteworthily, the theory shown in Fig. 5 catches the experimental results. The doublet is only observed in the right output field, while the left output field is insensitive to the ϕ variations. Figure 8(b1) shows that the right output field exhibits a clear doublet for a relative phase of about ±π. This doublet is observed until ϕ=±π/2, where it merges into a single peak. In this case, the dots and squares perfectly overlap, and in the region between π/2 and π/2, the peak resonant frequency shifts from the right to the left of ω0. On the contrary, as shown in Fig. 8(b2), the left output field does not exhibit any splitting. The dots and squares overlap for all ϕ values. As a result, the spectrum always shows a single Lorentzian. However, the presence of the small coefficient β210 gives rise to small oscillations of the resonant frequency. As an example, two experimental spectra and their theoretical fits are shown in Figs. 8(c1) and 8(c2) for ϕ=0.93π and 0.35π, respectively. Noteworthily, also in this case the model of Section 2 is in perfect agreement with the experimental results.

    Experimental results for a taiji microresonator. The red (blue) curves refer to the left (right) output field. The dashed and dash-dotted black lines show the fit with the theoretical model. Panels (a1) and (a2) show the output fields as a function of the frequency detuning Δω for a left and right excitation, respectively. Panels (b1) and (b2) display the doublet splitting (δω=ω</>−ω0) for a symmetric interferometric excitation as a function of the phase ϕ for the right or left output fields. ω</> are the measured frequency minima of the doublet. The magenta (purple) dots refer to δω=ω<−ω0 for the right (left) output fields. The black (gray) squares refer to δω=ω>−ω0 for the right (left) output fields. The dashed magenta (purple) and dash-dotted black (gray) curves show the theoretical results modeled by using the fitting parameters of the one-side excitation. Panels (c1) and (c2) show two examples of spectra for an interferometric excitation with a phase ϕ=0.93π and −0.35π, respectively. From the fit, we obtain the following coefficients: Γ=75.1±0.1 GHz, γ=42.6±0.1 GHz, β12=(−90.8±0.1)−i(24.17±0.03) GHz, β21=(8.7±0.1)−i(7.30±0.09) GHz.

    Figure 8.Experimental results for a taiji microresonator. The red (blue) curves refer to the left (right) output field. The dashed and dash-dotted black lines show the fit with the theoretical model. Panels (a1) and (a2) show the output fields as a function of the frequency detuning Δω for a left and right excitation, respectively. Panels (b1) and (b2) display the doublet splitting (δω=ω</>ω0) for a symmetric interferometric excitation as a function of the phase ϕ for the right or left output fields. ω</> are the measured frequency minima of the doublet. The magenta (purple) dots refer to δω=ω<ω0 for the right (left) output fields. The black (gray) squares refer to δω=ω>ω0 for the right (left) output fields. The dashed magenta (purple) and dash-dotted black (gray) curves show the theoretical results modeled by using the fitting parameters of the one-side excitation. Panels (c1) and (c2) show two examples of spectra for an interferometric excitation with a phase ϕ=0.93π and 0.35π, respectively. From the fit, we obtain the following coefficients: Γ=75.1±0.1  GHz, γ=42.6±0.1  GHz, β12=(90.8±0.1)i(24.17±0.03)  GHz, β21=(8.7±0.1)i(7.30±0.09)  GHz.

    4. CONCLUSION

    We theoretically model and experimentally validate the coherent interferometric excitation technique as a method to estimate the real and imaginary parts of the eigenvalues of a Hermitian and non-Hermitian two-level optical system. We apply it to the measurement of the transmission spectra of a bus waveguide/microresonator system where the microresonator is characterized by two degenerate counter-propagating optical modes. The variation of the relative phase and amplitude of the input fields causes the merging of the typical resonant doublet due to the counter-propagating microresonator modes into a single Lorentzian lineshape. This allows a direct estimation of the intermodal coupling coefficients from the transmission spectra in the Hermitian case. In the non-Hermitian case, it allows extrapolating the real part of the eigenvalues from the resonant frequencies and its imaginary part from the resonance shape, i.e., the FWHM. In particular, it is possible to estimate also the microresonator Q-factor in absence of the intermodal coupling between the counter-propagating modes, i.e., without the contribution of surface-wall roughness backscattering. Furthermore, we show that a taiji microresonator works on an exceptional point and that it has an asymmetric transmission under a coherent interferometric excitation. In fact, the transmitted output field from one side remains unaffected by changes in the relative input field phases, while the other transmitted field exhibits rich dynamics with a resonance splitting at characteristic ϕ. These differences are due to the presence of the inner S-shaped waveguide to which the propagating field couples only in one direction.

    These results represent a proof of concept of the coherent interferometric excitation technique. Indeed, the experimental setup is made of discrete optical components and lacks an absolute control of ϕ. However, the coherent interferometric excitation can be easily integrated in a single chip, where phase shifters, splitters, and integrated Mach–Zehnder interferometers can control both the phase and amplitudes of the input fields to excite the microresonator counter-propagating modes; see e.g., Ref. [33]. For this reason, our work paves the way to a family of devices based on the exceptional points for sensing applications where the accuracy of the spectral measurement is improved by the coalescence of the eigenvalues.

    The proposed method and analysis can be also generalized to study other two-level optical systems such as two coupled micoresonators or photonic molecules [34,35] where the complex interplay of different optical modes is difficult to unfold. In particular, the integration of a controlled coherent interferometric excitation into these schemes leads to the study of different kinds of degeneracies, which can be engineered in the spectrum of coupled optical microcavities [15,36]. As a result, the spectral response of an interferometric excitation allows us to directly identify purely absorptive or radiative exceptional points.

    Finally, since the interferometric excitation is based on measuring the linear response of Hermitian and non-Hermitian systems, it can be extended to multi-level structures having more than two input ports and different propagation directions. This allows for the study of higher-order exceptional points [37,38]. As an example, we might apply the coherent interferometric excitation method to three-level systems such as those discussed in Refs. [39,40]. Here, one deals with three input ports that excite three distinct optical cavities in the same propagation direction. It is thus necessary to split the input laser into three beams and to control at least two amplitudes and two phases in order to apply the interferometric excitation. Another example, which involves the excitation of coupled counter-propagating modes, is the two coupled-resonator optical waveguide (CROW) [41,42] in add and drop configuration. In this four-level device, the four ports can be excited simultaneously. Consequently, it is possible to probe the whole phase space in both symmetric and asymmetric interferometric excitation configurations by controlling three relative field amplitudes and three relative phases. These extensions cannot be achieved by the experimental setup discussed here, but are needed to design a dedicated photonic integrated circuit.

    Acknowledgment

    Acknowledgment. We gratefully thank Dr. Iacopo Carusotto and Mr. Alberto Muñoz de las Heras for useful inputs, valuable comments, and interesting discussion, as well as Mr. Enrico Moser for technical support. This work was supported by Q@TN, the joint lab between University of Trento, Fondazione Bruno Kessler (FBK), National Institute for Nuclear Physics (INFN), and the National Research Council (CNR).

    APPENDIX A: HERMITIAN COUPLING

    In this section, we show the Hermitian condition on the coupling coefficients and how it reduces the other characteristic parameters such as g1,2, eigenvalues, and eigenvectors.

    From Eq. (1) of the main text, we define the following backscattering matrix:K=(0iβ12iβ210).This K matrix satisfies the Hermitian condition, i.e., K=K+, if and only ifβ12=β21*=:β.Then, reformulating the backscattering coefficient as β=|β|eiφ is easy to obtain:g1=i|β|sin[θ]R[g1]=0,g2=|β|cos[θ]I[g2]=0,where θ=ϕ+φ=ϕ+arg[β]. Consequently, the eigenvectors changed in v1,2=11+|sin[θ]1cos[θ]|2(isin[θ]1cos[θ]1).

    On the other hand, the eigenvalues reduce to λ1,2=(ω0±|β|)i(γ+Γ). As mentioned in the main text, the eigenvectors depend just on the coefficient θ and are not affected by the strength of the coupling β. Furthermore, the inner product is equal to zero v1|v2=0, and therefore, when β=0, the system works on a diabolic point.

    APPENDIX B: EXPERIMENTAL PROCEDURE

    In this section, we disclose the experimental procedure used to obtain a symmetric interferometric excitation. First, we measured the losses of each component of the experimental setup and we calibrated the responses of the detectors. Then we coupled light into the microresonator and, using the detectors PDC1, PDC2, PDL, and PDR and VOAs, we set the intensity in the two arms almost equal and low enough not to observe nonlinear effects. At this point, the first step is to balance the length of the two optical paths after the first fiber splitter. This is needed to decrease the relative phase variation between the two input fields as a function of wavelength (2πnΔL/λ), and thus to achieve an almost constant relative phase ϕ for a microresonator resonant wavelength. Therefore, we performed repeated wavelength scans, using a scan rate of 100 nm/s, over the entire spectrum and changed the optical length of the right arm to find the condition where the period of the observed oscillations at the PDR and PDL detectors is maximum.

    Once we set the delay line, we used the following procedure to input the same intensities in the bus waveguide sides, i.e., to perform a symmetric interferometric excitation. System.Xml.XmlElementSystem.Xml.XmlElementSystem.Xml.XmlElementSystem.Xml.XmlElementSystem.Xml.XmlElementSystem.Xml.XmlElementSystem.Xml.XmlElementSystem.Xml.XmlElement

    References

    [1] 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).

    [2] Y. Ashida, Z. Gong, M. Ueda. Non-Hermitian physics. Adv. Phys., 69, 249-435(2020).

    [3] M. V. Berry. Physics of nonhermitian degeneracies. Czech. J. Phys., 54, 1039-1047(2004).

    [4] M. Berry, R. Bhandari, S. Klein. Black plastic sandwiches demonstrating biaxial anisotropy. Eur. J. Phys., 20, 1-14(1999).

    [5] Ş. K. Özdemir, S. Rotter, F. Nori, L. Yang. Parity-time symmetry and exceptional points in photonics. Nat. Mater., 18, 783-798(2019).

    [6] M.-A. Miri, A. Alã. Exceptional points in optics and photonics. Science, 363, eaar7709(2019).

    [7] L. Feng, X. Zhu, S. Yang, H. Zhu, P. Zhang, X. Yin, Y. Wang, X. Zhang. Demonstration of a large-scale optical exceptional point structure. Opt. Express, 22, 1760-1767(2014).

    [8] A. Calabrese, F. Ramiro-Manzano, H. M. Price, S. Biasi, M. Bernard, M. Ghulinyan, I. Carusotto, L. Pavesi. Unidirectional reflection from an integrated ‘taiji’ microresonator. Photon. Res., 8, 1333-1341(2020).

    [9] Z. Lin, H. Ramezani, T. Eichelkraut, T. Kottos, H. Cao, D. N. Christodoulides. Unidirectional invisibility induced by PT-symmetric periodic structures. Phys. Rev. Lett., 106, 213901(2011).

    [10] 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).

    [11] B. Peng, Ş. K. Özdemir, M. Liertzer, W. Chen, J. Kramer, H. Ylmaz, J. Wiersig, S. Rotter, L. Yang. Chiral modes and directional lasing at exceptional points. Proc. Natl. Acad. Sci. USA, 113, 6845-6850(2016).

    [12] A. M. de las Heras, I. Carusotto. Unidirectional lasing in nonlinear taiji micro-ring resonators. Photon. Res., 8, 1333-1341(2020).

    [13] L. Jin, Z. Song. Incident direction independent wave propagation and unidirectional lasing. Phys. Rev. Lett., 121, 073901(2018).

    [14] H. Hodaei, M.-A. Miri, M. Heinrich, D. N. Christodoulides, M. Khajavikhan. Parity-time-symmetric microring lasers. Science, 346, 975-978(2014).

    [15] C. Wang, W. R. Sweeney, A. D. Stone, L. Yang. Coherent perfect absorption at an exceptional point. Science, 373, 1261-1265(2021).

    [16] A. Li, T. Van Vaerenbergh, P. De Heyn, P. Bienstman, W. Bogaerts. Backscattering in silicon microring resonators: a quantitative analysis. Laser Photon. Rev., 10, 420-431(2016).

    [17] A. M. de las Heras, R. Franchi, S. Biasi, M. Ghulinyan, L. Pavesi, I. Carusotto. Nonlinearity-induced reciprocity breaking in a single nonmagnetic taiji resonator. Phys. Rev. Appl., 15, 054044(2021).

    [18] Y. G. N. Liu, O. Hemmatyar, A. U. Hassan, P. S. Jung, J.-H. Choi, D. N. Christodoulides, M. Khajavikhan. Engineering interaction dynamics in active resonant photonic structures. APL Photon., 6, 050804(2021).

    [19] S. Fan, W. Suh, J. D. Joannopoulos. Temporal coupled-mode theory for the Fano resonance in optical resonators. J. Opt. Soc. Am. A, 20, 569-572(2003).

    [20] W. Suh, Z. Wang, S. Fan. Temporal coupled-mode theory and the presence of non-orthogonal modes in lossless multimode cavities. IEEE J. Quantum Electron., 40, 1511-1518(2004).

    [21] Q. Li, T. Wang, Y. Su, M. Yan, M. Qiu. Coupled mode theory analysis of mode-splitting in coupled cavity system. Opt. Express, 18, 8367-8382(2010).

    [22] S. Biasi, F. Ramiro-Manzano, F. Turri, P.-E. Larré, M. Ghulinyan, I. Carusotto, L. Pavesi. Hermitian and non-Hermitian mode coupling in a microdisk resonator due to stochastic surface roughness scattering. IEEE Photon. J., 11, 6101114(2019).

    [23] R. Franchi, S. Biasi, A. M. de las Heras, M. Ghulinyan, I. Carusotto, L. Pavesi. Influence of the bus waveguide on the linear and nonlinear response of a taiji microresonator. Opt. Express, 29, 29615-29630(2021).

    [24] Q. Li, Z. Zhang, J. Wang, M. Qiu, Y. Su. Fast light in silicon ring resonator with resonance-splitting. Opt. Express, 17, 933-940(2009).

    [25] R. J. Potton. Reciprocity in optics. Rep. Prog. Phys., 67, 717-754(2004).

    [26] D. Jalas, A. Petrov, M. Eich, W. Freude, S. Fan, Z. Yu, R. Baets, M. Popović, A. Melloni, J. D. Joannopoulos, M. Vanwolleghem, C. R. Doerr, H. Renner. What is-and what is not-an optical isolator. Nat. Photonics, 7, 579-582(2013).

    [27] S. Biasi, P. Guillemé, A. Volpini, G. Fontana, L. Pavesi. Time response of a microring resonator to a rectangular pulse in different coupling regimes. J. Lightwave Technol., 37, 5091-5099(2019).

    [28] J. Heebner, R. Grover, T. Ibrahim. Optical Microresonators: Theory, Fabrication, and Applications, 138(2008).

    [29] J. Wiersig. Review of exceptional point-based sensors. Photon. Res., 8, 1457-1467(2020).

    [30] W. Chen, Å. K. Ãzdemir, G. Zhao, J. Wiersig, L. Yang. Exceptional points enhance sensing in an optical microcavity. Nature, 548, 192-196(2017).

    [31] M. Borghi, D. Bazzanella, M. Mancinelli, L. Pavesi. On the modeling of thermal and free carrier nonlinearities in silicon-on-insulator microring resonators. Opt. Express, 29, 4363-4377(2021).

    [32] R. Marchetti, C. Lacava, A. Khokhar, X. Chen, I. Cristiani, D. J. Richardson, G. T. Reed, P. Petropoulos, P. Minzioni. High-efficiency grating-couplers: demonstration of a new design strategy. Sci. Rep., 7, 16670(2017).

    [33] M. Mancinelli, P. Bettotti, J. M. Fedeli, L. Pavesi. Interferometric switching in coupled resonator optical waveguides-based reconfigurable optical device. Opt. Lett., 38, 217-219(2013).

    [34] S. V. Boriskina. Photonic Molecules and Spectral Engineering, 393-421(2010).

    [35] M. Borghi, A. Trenti, L. Pavesi. Four wave mixing control in a photonic molecule made by silicon microring resonators. Sci. Rep., 9, 408(2019).

    [36] C. Chen, L. Jin, R.-B. Liu. Sensitivity of parameter estimation near the exceptional point of a non-Hermitian system. New J. Phys., 21, 083002(2019).

    [37] M. A. Quiroz-Juárez, A. Perez-Leija, K. Tschernig, B. M. Rodríguez-Lara, O. S. M. Loaiza, K. Busch, Y. N. Joglekar, R. de J. León-Montiel. Exceptional points of any order in a single, lossy waveguide beam splitter by photon-number-resolved detection. Photon. Res., 7, 862-867(2019).

    [38] M. Am-Shallem, R. Kosloff, N. Moiseyev. Parameter estimation in atomic spectroscopy using exceptional points. Phys. Rev. A, 93, 032116(2016).

    [39] 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).

    [40] Y. Wu, P. Zhou, T. Li, W. Wan, Y. Zou. High-order exceptional point based optical sensor. Opt. Express, 29, 6080-6091(2021).

    [41] B. Little, S. Chu, P. Absil, J. Hryniewicz, F. Johnson, F. Seiferth, D. Gill, V. Van, O. King, M. Trakalo. Very high-order microring resonator filters for WDM applications. IEEE Photon. Technol. Lett., 16, 2263-2265(2004).

    [42] T. Lei, A. W. Poon. Modeling of coupled-resonator optical waveguide (crow) based refractive index sensors using pixelized spatial detection at a single wavelength. Opt. Express, 19, 22227-22241(2011).

    Stefano Biasi, Riccardo Franchi, Filippo Mione, Lorenzo Pavesi, "Interferometric method to estimate the eigenvalues of a non-Hermitian two-level optical system," Photonics Res. 10, 1134 (2022)
    Download Citation