Issue 
EPJ Appl. Metamat.
Volume 9, 2022
Metamaterials for Novel Wave Phenomena in Microwaves, Optics, and Mechanics



Article Number  14  
Number of page(s)  15  
DOI  https://doi.org/10.1051/epjam/2022006  
Published online  01 July 2022 
https://doi.org/10.1051/epjam/2022006
Research Article
How to achieve exceptional points in coupled resonators using a gyrator or PTsymmetry, and in a timemodulated single resonator: high sensitivity to perturbations^{★}
^{1} Department of Electrical Engineering and Computer Science, University of California, Irvine, CA 92697, USA
^{2} Department of Mathematics, University of California, Irvine, CA 92697, USA
^{*} email: f.capolino@uci.edu
Received:
6
December
2021
Accepted:
25
February
2022
Published online: 1 July 2022
We study the rise of exceptional points of degeneracy (EPD) in various distinct circuit configurations such as gyratorbased coupled resonators, coupled resonators with PTsymmetry, and in a single resonator with a timevarying component. In particular, we analyze their high sensitivity to changes in resistance, capacitance, and inductance and show the high sensitivity of the resonance frequency to perturbations. We also investigate stability and instability conditions for these configurations; for example, the effect of losses in the gyratorbased circuit leads to instability, and it may break the symmetry in the PTsymmetrybased circuit, also resulting in instabilities. Instability in the PTsymmetry circuit is also generated by breaking PTsymmetry when one element (e.g., a capacitor) is perturbed due to sensing. We have turned this instability “inconvenience” to an advantage, and we investigate the effect of nonlinear gain in the PTsymmetry coupledresonator circuit and how this leads to an oscillator with oscillation frequency very sensitive to perturbation. The circuits studied in this paper have the potential to lead the way for a more efficient generation of highsensitivity sensors that can detect very small changes in chemical, biological, or physical quantities.
Key words: Coupled resonators / exceptional points of degeneracy (EPDs) / gyrator / perturbation theory / sensor / timemodulation
© A. Nikzamir et al., Published by EDP Sciences, 2022
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
In order to improve the detection limit for small perturbation sensors, an operation based on the exceptional point of degeneracy concept (EPD) can be a valuable option. EPD is a point at which two or more eigenvalues and the corresponding eigenvectors of the system coalesce (i.e., they are degenerate) [1–14]. The main feature of an exceptional point is the strong full degeneracy of the relevant eigenmodes (including their eigenvectors), justifying the presence of “D” in EPD that stands for “degeneracy” [15]. Exceptional points have been observed in various circuits as in coupled resonators [16–22], temporallyperiodic systems [23–25], and spatially periodic structures [26–28].
Sensors are used widely and extensively in many industrial, automotive, and medical applications. In recent years, numerous approaches have been used to sense variations of physical, biological, or chemical changes, e.g., to sense pressure [16,29,29], temperature [30], humidity [31], electron beam velocity [32], and chemical or biological quantities [33–36]. Since sensitivity is an essential feature of measurement systems, various types of sensors are needed to sense. Thus, lowcost, simple, and highlysensitive sensors are desirable to measure different quantities. In conventional sensors, the perturbation by a small value Δ results in shifting the system’s eigenfrequency by an amount in the same order of Δ (linear behavior), like perturbing a simple LC resonator (or tank). It means that in conventional sensing applications, the perturbation and the measurable changes, like frequency, follow the same trend. To increase the sensitivity, EPDs have offered a new method. Rather than typical detecting systems, where the eigenfrequency changes are of the same order of the perturbation Δ, the change of an eigenfrequency Δω in a system working at an EPD of order two follows the behavior Δω∝ Δ^{1/2} as shown in [37–40]. EPDs are found in many physical settings under PTsymmetry properties [7,16–19,41]. However, EPDs are also found in more general forms, which do not require a system to satisfy PTsymmetry [23–25,37,42]. Breaking the symmetry due to the external perturbance eigenvalues split from the degenerated value can be used as a detecting system.
This paper treats EPDs using four methods: First, by using the circuit’s characteristic equation and showing the degeneracy of the eigenmodes. Second, the coalescence of eigenvectors, which is observed analytically. Third, timedomain results obtained from simulations show the linear growth revealing the secondorder EPD. The related frequency response is also studied, and it is shown how it is associated with the EPD. Fourth, the approximate fractional power expansion series using the Puiseux series shows the bifurcation and square rootlike behavior of the eigenvalues with respect to perturbations.
To address the effectiveness of the discussed EPD circuits, we consider the systems’ sensitivity to perturbations of capacitance, inductance, and resistance, which are used as sensing elements, depending on the application. In Section 2, we provide an overview of the three circuit configurations leading to secondorder EPDs. In Section 3, we describe the secondorder EPD in the gyratorbased circuit with parallel configuration in lossless and lossy cases. Also, we study the stability of the circuit, manifested in its eigenfrequencies. Details and analysis of the series configuration, not shown here, are discussed in [21,22,43,44]. In Section 4, we study the secondorder EPD in two coupled resonators with balanced gain and loss satisfying PTsymmetry. In this circuit, we consider nonlinear gain, and by perturbing only one side of the coupled resonator, we break PTsymmetry making the system unstable and causing oscillation. In [16], when a perturbation occurs on one side of the coupled resonator circuit, the other side is tuned accordingly to keep PTsymmetry. This procedure made the operational regime difficult to implement since the exact value of the changes should be known a priori. Moreover, the sign of the perturbation was consistent with the bifurcation direction, which means that only either positive or negative changes in the circuit’s quantities could be sensed, not both. Instead, working as in the oscillator regime enables sensing of perturbations with both positive and negative signs. Also, when using nonlinear gain and the oscillatory regime, there is no need to tune a circuit’s side to keep symmetry, as we show in Section 4. Nonlinearity plays a critical role in this circuit operational regime and helps us find the EPD more easily. While the circuit with EPD and nonlinearity is sensitive to any perturbation sign, the nonlinearity keeps the circuit at an EPD even with a 1% missmatch between the gain and loss. Finally, in Section 5, we study a single LC resonator with a timevarying element and show how this simple configuration leads to an EPD by just tuning the modulation frequency. In order to find the EPD, we assume a timevarying capacitor connected to a fixed inductance. Then, the loss effect is taken into account in the analysis. We study the eigenfrequency dispersion by varying modulation frequency, where the LC tank is connected to loss or gain. Also, EPD is found in the timemodulated circuit when the LC tank is terminated with timevarying gain or loss.
2 Secondorder EPD in three distinct possible configurations
We discuss here three different methods to obtain a secondorder EPD. First, the degeneracy of the eigenvalues (i.e., eigenfrequencies) occurs in a gyratorbased circuit where two LC tanks are coupled with a gyrator, as shown in Figure 1a. Second, we study two mutually coupled LC circuits with balanced gain and loss, satisfying PTsymmetry, shown in Figure 1b.
In Figures 1a and 1b, both circuits support square root sensitive behavior to perturbations of the resistance, capacitance, and inductance. However, some differences must be observed: for example, the conventional PTsymmetry circuit in Figure 1b needs the tuning on both sides to keep the PTsymmetry and avoid selfoscillations (arising from nonzero Im (ω)). In other words, it is not possible to change one side only and observe two purely real eigenfrequencies. Instead, the gyratorbased circuit shows that the perturbation on only one side leads to two purely real eigenfrequencies. Thus, there is no need to tune the other side after a perturbation in the gyratorbased circuit. However, any gain or loss in the gyratorbased circuit will cause oscillation and instability, as we will discuss later in more detail.
We show how to turn the instability of a circuit (e.g., the one in Figure 1b) to our advantage. Using nonlinear gain, we make the circuit oscillate and saturate, and by perturbing the capacitance on one side, a square rootlike change in oscillation frequency is detected. Working in an oscillator regime has certain advantages compared to the conventional PTsymmetry regime of operations in the previous literature [16]. For instance, there is no need to tune the circuit to reach PTsymmetry again after the perturbed values (of a capacitor, for example) are measured in this circuit. Also, nonlinearity helps to finetune the circuit to the EPD in an easier way. We discuss these concepts in more detail later on in Section 4.
The third circuit we discuss here, shown in Figure 1c, is a linear timevarying (LTV) system in which an inductor is connected to the timevarying capacitor. This configuration does not need any negative components to realize the EPD, like a negative capacitance and inductance in the gyratorbased circuit or an active gain element in the PTsymmetric coupledresonator circuit. In addition, we need only one timevarying resonator in this third scheme rather than two fixed resonators. Similarly, an EPD is found in the LTV circuit when an inductor and capacitor are connected in parallel to the timevarying loss or gain.
Fig. 1 Three different configurations where support EPD (a) Two resonators coupled through gyrator (b) Two mutually coupledresonator terminated with balanced gain and loss (PTsymmetry circuit) (c) Timevarying circuit which just one component in the circuit needed to change over time periodically. The timevarying component here could be capacitor or inductor or loss. 
3 EPD in gyratorbased circuit
In this section, we study the first scheme to obtain EPD by using two coupled LC tanks connected through a gyrator. An ideal gyrator is a nonreciprocal linear twoport device whose current on one port is related to the voltage on the other port. More details about the gyrator and various realization methods are discussed in references [45–48]. The instantaneous relations between voltages and currents on the gyrator are described by (1)
where the gyration resistance R_{g} has a unit of Ohm with the direction indicated by an arrow in the circuit.
We find the eigenvalues (i.e., the eigenfrequencies) and demonstrate the condition for obtaining an EPD at the desired frequency. Finally, we show the sensing potentials by applying a perturbation, and we study the effects of losses on the stability of eigenfrequencies. In addition, by using a timedomain circuit simulator, we verify the circuit behavior predicted by the theoretical calculations and also show that the eigenfrequencies can be predicted by using the Puiseux fractional power series expansion.
3.1 Lossless configuration
Two parallel LC tanks are coupled by a gyrator, as shown in Figure 2a. By writing the circuit equations and defining the state vector as leads to (2)
and is the circuit matrix. The eigenfrequencies of this circuit are calculated by solving the characteristic equation [22,49], (4)
The characteristic equation is quadratic in ω^{2} and all the coefficients are real, so both ω and −ω and ω and ω^{*} are solutions. The angular eigenfrequencies are determined as (5) (6) (7)
According to equation (5), a necessary condition for an EPD to occur is b = 0, which results in an EPD angular frequency of . So, we rewrite equation (7) as (8)
In order to obtain an EPD with real angular frequency, we consider the case with purely real value for ω_{01} and ω_{02}, so the value of either C_{1} or C_{2} should be negative. As a result, to have a real value for ω_{01} and ω_{02}, one resonator needs to be composed of both negative C and L, and more details are discussed in [43,44]. Another scenario with an unstable uncoupled resonator is conceivable, which was studied for a series configuration in [21]. As an example, here we use the following values for the components shown in Figure 2a: L_{1} = 100 μH, L_{2} = − 100 μH, C_{2} = − 100 nF, and R_{g} = 50 Ω. Then, the positive capacitance C_{1} is found by solving the quadratic equation obtained from the EPD condition. Since the equation of the EPD condition is quadratic, it would yield two answers for C_{1}. In this paper, we select the value that leads to real EPD frequency, which is C_{1} = 13.51 nF. Then the corresponding value for EPD angular frequency is calculated as ω_{e} = 5.22 × 10^{5} rad/s. In this circuit, an opampbased invertor could realize the negative components. The designed circuit to obtain negative impedance is shown in Figure 2b, which converts the impedance Z_{l} (ω) to Z_{invert} (ω) = − Z_{l} (ω). Thus, we can provide the required negative capacitance and inductance by employing that configuration.
The real and imaginary parts of perturbed eigenfrequencies normalized to the EPD angular frequency by varying R_{g} are shown in Figures 2c and 2d, and analogous results by perturbing the positive capacitance C_{1} are shown in Figures 2e and 2f. In addition, the eigenfrequencies are well approximated by using the Puiseux fractional power series expansion (green dashed lines). Appendix A provides the concept and formulas for this method. The approximated results obtained by the Puiseux series show an excellent agreement with the “exact” values calculated directly from the eigenvalue problem. The coefficients of the Puiseux series up to secondorder for the mentioned example are calculated as, α_{1} = 3.85 × 10^{5} rad/s, and α_{2} = 1.42 × 10^{5} rad/s when perturbing R_{g}, and α_{1} = j2.07 × 10^{5} rad/s, and α_{2} = − 1.72 × 10^{5} rad/s when perturbing C_{1}. The bifurcation of the real part of the eigenfrequencies, which indicate the stable sensing region, is observed when R_{g} > R_{g,e}, and C_{1} < C_{1,e}.
Timedomain simulation result for the voltage v_{1} (t) is obtained using the Keysight Advanced Design System (ADS) timedomain circuit simulator, and the result is plotted in Figure 2g. The frequency spectrum corresponding to the simulated timedomain voltage is found by taking the Fast Fourier Transform (FFT), as shown in Figure 2h. The results are obtained using the initial voltage of 1 mV on the left capacitor C_{1}. We observe that the voltage increases linearly with time. As well known in circuit theory, this linear growth indicates that two system eigenvalues collided, and a double pole describes the system response. This is a specific property of a secondorder EPD. The oscillation frequency is ω_{o} = 5.22 × 10^{5} rad/s, which is exactly equal to EPD angular frequency ω_{e}.
We now observe the EPD in a frequencydomain analysis by calculating the circuit’s total input admittance Y_{total} (ω), as shown in Figure 2a. We define the two admittances of two LC tanks as Y_{1} = jωC_{1} + 1/(jωL_{1}) and Y_{2} = jωC_{2} + 1/(jωL_{2}). Then the transferred admittance of the Y_{2} to the left side is . The total admittance observed from the circuit input port is calculated by (9)
The resonant angular frequencies are achieved by imposing Y_{total} (ω) = 0; the normalized resonance frequencies to the EPD angular frequency by varying the gyration resistance are shown in Figure 2i. Two zeros of total admittance coincide exactly at EPD.
Fig. 2 (a) The gyratorbased circuit with the ideal gyrator connected two parallel LC tanks. (b) Opampbased circuit configuration to obtain negative inductance and capacitance. The sensitivity of the (c), (e), real and (d), (f), imaginary parts of the eigenfrequencies to (c), (d) gyration resistance, (e), (f) positive capacitance C_{1}. Solid lines: solution of eigenvalue problem of Eq. (2); greendashed lines: Puiseux series approximation truncated to its second term. Voltage of positive capacitance v_{1} (t) at EPD (g) timedomain, and (h) frequencydomain. The frequencydomain result is calculated by applying an FFT with 10^{6} samples in the time window of 0 to 3 ms. (i) Root locus of zeros of Y_{total} (ω) = 0 showing the real versus imaginary parts of resonance frequencies by perturbing gyration resistance. At the EPD, two zero collide at ω_{e} and the system’s total admittance has the form of Y_{total} (ω) ∝ (ω − ω_{e}) ^{2}. 
3.2 Lossy configuration and stability
In this section, we study the loss effects on the eigenfrequencies of the gyratorbased circuit. Two parallel resistors R_{1} and R_{2} are connected to both resonators, as shown in Figure 3a. By writing down the circuit equations and assuming the same state vector as introduced in the lossless case, the associated Liouvillian formalism reads as (10) (11)
where γ_{1} = 1/(R_{1} C_{1}) and γ_{2} = 1/(R_{2} C_{2}) represent the losses made by resistors. The eigenfrequencies are found by solving the characteristic equation, (12)
Here, eigenfrequencies ω and −ω^{*} are both roots of the characteristic equation. In order to have a stable circuit, eigenfrequencies should be purely real, but the characteristic equation in equation (12) has some imaginary coefficients. Here, to have purely real eigenfrequencies in the lossy circuit, the oddpower terms of the angular eigenfrequency in the characteristic equation should vanish. Otherwise, a complex eigenfrequency is needed to satisfy the characteristic equation. There is no condition to make both ω and ω^{3} coefficients equal to zero [21]. Hence, eigenfrequencies are always complex, leading to instabilities that cause oscillations in the circuit.
By considering the same value for components as already used in the lossless case, the evolution of eigenfrequencies is shown in Figures 3b and 3c. In these plots, loss on the first resonator is changed, and loss on the second resonator is eliminated. Moreover, in Figures 3d and 3e, we perturb the loss on the second resonator while the loss on the first resonator is removed. When both losses are zero, the system has the same EPD frequency of a lossless configuration, but perturbed eigenfrequencies are complexvalued for any amount of losses. So, the lossy circuit oscillates at the frequency associated with the real part of the unstable eigenfrequency. Also, the eigenfrequency is extremely sensitive to either positive or negative variations in the parallel resistances (square root behavior due to the perturbation). A working option is based on preventing the circuit from reaching saturation by switching off the circuit and operating on the circuit’s transient response, as was done in [23] for an EPD based on a time modulated circuit.
Fig. 3 (a) Schematic view of the lossy gyratorbased circuit, with a resistor in each resonator. Variation of (b) real and (c) imaginary parts of the angular eigenfrequencies to a resistor perturbation on the left resonator, i.e., when γ_{1} changes and γ_{2} = 0. Variation of (d) real and (e) imaginary parts of the angular eigenfrequencies to a resistor perturbation on the right resonator, i.e., when −γ_{2} changes and γ_{1} = 0. 
4 EPD in PTsymmetric coupled resonators and nonlinearity effects
This section discusses the EPD in two mutually coupled resonators based on PTsymmetry. This is the circuits studied so far by most of the researchers in the last decade [7,16–18,20]. We show the occurrence of an EPD by using the concept of the eigenvector coalescence parameter. Moreover, we study the resonance condition when the total admittance of the circuit is equal to zero (i.e., the double zero condition). The negative conductance in the analyzed circuit could be achieved via crosscoupled or opampbased circuits. The negative conductance obtained from these transistorbased circuits has nonlinearity effects due to the saturation. Thus, the nonlinearity in negative conductance would alter the circuit operation, as discussed later on. We model the nonlinearity with a cubic iv characteristic and show the timedomain analysis and frequency responses by using timedomain simulations that are the right tool when nonlinearity is present.
Moreover, we stress the EPD sensitivity characteristic and provide an example where the sensing scheme involves the perturbation of a capacitance. It means that the PTsymmetry is broken. We demonstrate the high sensitivity behavior of the eigenfrequencies for either positive or negative changes in capacitance and show that the system becomes unstable. In [16], the authors discussed sensitivity using two PTsymmetric coupled resonators. They demonstrated the sensitivity of the eigenfrequencies of the circuit due to the capacitance perturbation while they kept the PTsymmetry configuration. Thus, they needed to tune the other (nonsensing) side of the circuit to have balanced capacitance on both sides to keep the PTsymmetry even after each sensing operation. It means that the exact value of the changes in the sensing capacitance should be somehow known to tune the other side, which is not possible in practical sensing scenarios. Also, they could only measure the perturbation in the bifurcation direction (i.e., only the negative (or positive) capacitance changes, based on the design). Both positive and negative capacitance perturbation sensing ranges should be desirable, and the tuning process should be made easier since there is no priori knowledge of the sensing capacitance variation. Hence, it is not possible to keep the system PTsymmetric while sensing unless possible iterative schemes are researched that guess the unknown capacitance value. Finally, we confirm the eigenfrequency’s sensitivity and square root behavior to the perturbation by using the Puiseux fractional power series expansion.
4.1 EPD in mutual coupled resonators with PTsymmetry
Two coupled LC tanks terminated on the left side with a gain given by the negative conductance −G_{1} and terminated on the right side with loss G_{2} are illustrated in Figure 4a. By writing Kirchhoff’s current law, we obtain the two equations for the circuit (13)
where Q_{1} is the capacitor charge on the gain side (left resonator), Q_{2} is the capacitor charge on the loss side (right resonator), and k = M/L is the transformer coupling coefficient. In addition, , , , and are the first and the second time derivatives of the capacitors’ charge. We define the state vector as where superscript T denotes the transpose operation. Therefore, the circuit evolution is described by (14) (15)
Assuming signals in the form of Q_{n} ∝ e^{jωt}, C_{1} = C_{2} = C_{0} and G_{1} = G_{2} = G; we get the eigenfrequencies of the circuit by solving the characteristic equation, , leading to (16)
In the above equations, and According to equation (17), the required condition to obtain an EPD is b = 0, which leads to an EPD angular frequency of (18)
In the presented example, we use L = 1 μH, C_{1} = C_{2} = 1 nF, where inductors are mutually coupled via k = 0.2, and terminated with balanced gain and loss G_{1} = G_{2} = G. After solving the eigenvalue problem, the evolution of real and imaginary parts of the eigenfrequencies are shown in Figures 4b and 4c. The eigenfrequencies of the circuit coalesce at a specific balanced gain/loss value of γ = γ_{e} = 0.205, where γ_{e} that leads to an EPD is derived as (19)
To validate the results, we assume γ = 0, which means there is no gain or loss in the system and it results in two simple mutual LC tanks. The mentioned circuit has two pairs of eigenfrequencies and . If we remove the coupling, i.e., k = 0, the eigenfrequencies are equal to the independent LC tank circuits ±ω_{0}. For the values, 0 < γ < γ_{e}, the system’s eigenfrequencies are purely real, and the system has two fundamentals real eigenfrequencies. For the values, γ_{e} < γ, the two eigenfrequencies are complex conjugate, and system solutions grow or damp depending on the sign of the imaginary part of the angular eigenfrequencies; the system exhibits self (unstable) oscillations at the frequency associated with the real part of the eigenfrequency. The eigenvector coalescence factor is defined to evaluate how the circuit’s operation point is close to an EPD and measure the coalescence of two eigenvectors. It is defined as C. F. = sin(θ), where cos(θ) is (20)
In the determined equation ⟨, ⟩ is the inner product and ∥ ∥ denotes the eigenvector norm. The coalescence factor for the presented example is shown in Figure 4d. As we observe in this plot, two eigenvectors have coalesced at the corresponding value for EPD.
Fig. 4 (a). Two mutually coupled LC tanks terminated with balanced gain on the left and loss on the right side G_{1} = G_{2} = G (b) real and (c) imaginary parts of evaluated eigenfrequencies by perturbing G. In the illustrated plots, only eigenvalues with the positive real parts are drawn. (d) The coalescence factor corresponding to calculated eigenvectors which showing the degeneracy of two eigenvectors exactly at G_{e}(e) Root locus of zeros of Y_{in} (ω) − G = 0 showing the real and imaginary parts of resonance frequencies of the circuit when perturbing both load resistance and gain G. At the EPD, the system’s total admittance is Y_{in} (ω) − G_{e} ∝ (ω − ω_{e}) ^{2}; hence it shows a double zero at ω_{e}. 
4.2 Root locus of zeros of admittance
In this section, the resonance condition based on the vanishing of the total admittance is studied. We find the admittance Y_{in}, shown in Figure 4a, and demonstrate its double zero at the EPD. The resonance condition for this circuit is expressed as (21)
Here, the circuit is PTsymmetric, assuming linear gain and loss with G_{1} = G_{2} = G. We calculate the eigenfrequency by finding the zeros of the Y_{in} (ω) − G, which results in the same eigenfrequencies obtained from . From the zeros trajectory, both ω and −ω and ω and ω^{*} are solutions of equation (21), and we only show the eigenfrequencies with positive real value in Figure 4e.
4.3 Nonlinear gain and oscillator characteristics
In this section, we discuss the oscillator characteristics and nonlinear gain effects in two wireless coupled resonators, as in Figure 5. The transient and frequency response of the system are discussed when using a cubic (nonlinear) negative conductance model of the gain element. The system’s parameters are the same as in the previous section where G_{1} = G_{2} = G_{e} = 6.49 ms, and the EPD angular frequency is ω_{e} = 3.19 × 10^{6} rad/s. The relation between voltage and current of the nonlinear negative conductance is [50] (22)
where −G_{1} is the smallsignal negative conductance and α = G_{1}/3 is a thirdorder nonlinearity that is related to the active device’s saturation. We now assume the smallsignal nonlinear gain G_{1} to be slightly bigger than the balanced loss, as G_{1} = 1.001G_{2} = 1.001G_{e} to make the circuit slightly unstable (slightly breaking PTsymmetry). The timedomain response and frequency response obtained from Keysight ADS timedomain circuit simulator are shown in Figures 5b and 5c, where the circuit operates in the proximity of the EPD. To show the sensitivity of this oscillator, we perturb the capacitor C_{2} on the lossy side by 0.5%. The oscillation frequency at each perturbation, shown with the black dots, is found by taking the FFT of the timedomain voltage signal at the capacitor C_{1} (on the gain side) after reaching saturation. The FFT is calculated by using 10^{6} samples in the time window of 1000 periods after saturation, for each perturbed case. The system shows a distinct saturated oscillation frequency at each perturbed capacitor value. The frequency shift from the case without perturbation could be easily measured. For the sake of comparison with the linear case, solid blue and red lines in Figures 5d and 5e show the eigenfrequency evolution versus load capacitance perturbation as ΔC_{2} = (C_{2} − C_{2,e})/C_{2,e} (C_{2,e} is the capacitor’s value at EPD) by assuming linear gain for −G_{1}. The green dashed lines show the Puiseux fractional power series expansion truncated to its second order, which exhibits the squarerootlike sensitivity of the eigenvalues to a perturbation. Using both linear and nonlinear gain, the resonance frequency behavior shows the square root variation with capacitance perturbation. The difference in the frequency values between the nonlinear timedomain simulation and theoretical eigenvalue solutions arises from the nonlinearity and the subsequent saturation regime. The Puiseux series coefficients are calculated as, α_{1} = 5.35 × 10^{6} − j4.84 × 10^{6} rad/s, and α_{2} = − 7.90 × 10^{6} − j1.62 × 10^{6} rad/s.
The use of nonlinear gain in the circuit and the saturation effects make the EPD sensing regime robust. An errorcorrection method is discussed to enhance the robustness of sensing using nonlinearity in [51]. Also, the nonlinearity works as a selfcorrecting process in two coupled optical ring resonators in [52]. Nonlinearity in our proposed circuit helps maintain the oscillation frequency at the EPD frequency, within a range of small mismatches between gain and loss. The results obtained from the simulation shown in Figure 5f demonstrate that even with a 1% mismatch between gain and loss, the circuit oscillates at the same frequency as the case with balanced loss and gain. The red curve shows the oscillation frequency for the system with an exact match between the (nonlinear) smallsignal gain and loss (i.e., G_{1} = G_{2}), the dashed green is related to the case with G_{1} = 1.001G_{2}, and the black points denote the case with even larger smallsignal gain, G_{1} = 1.01G_{2}. They all show the same results in terms of shifted oscillation frequency versus C_{2} perturbation.
Fig. 5 (a). Two mutually coupled LC tanks terminated with nonlinear gain on the left and linear loss resistance on the right side where always G_{1} = 1.001 × G_{2}. (b) Time domain response. (c) Frequency response of the gainside capacitor voltage at the EPD point. (d) Real and (e) imaginary parts of the eigenfrequencies versus load capacitance perturbation ΔC_{2}; solid line: result from ; green dashed line: Puiseux fractional power series expansion terminated to its second order; black dots are obtained from the nonlinear time domain simulation. (f) Oscillation frequency versus load capacitance perturbation for three different mismatches between gain and loss δ = 0, 0.001, 0.01: the three different gains provide the same saturated oscillation frequencies. 
5 EPD in a timevarying single resonator
We now discuss a completely different way to obtain an EPD in a single resonator containing a timevarying element. It can be used as a highly sensitive circuit. As in the PTsymmetry system with balanced gain and loss discussed in the previous section, the EPD’s highly sensitive characteristics are also found in the timevarying single resonator, without the need for a gain component. For instance, we show an EPD in a periodic timevarying simple LC circuit in Figure 6a. We summarize the general formulation that can be applied to electronic and optical circuits. By assuming the state vector Ψ (t) = [Ψ_{1} (t), Ψ_{2} (t)] ^{T} where Tis the transpose operator. The differential equation describing the state vector time evolution is (23)
where (t) is the 2 × 2 timevariant system matrix. Here, the system matrix changes periodically in time, unlike the cases discussed in Sections 3 and 4. Thus, we do not have only one system matrix, and we employ timeperiodic analysis to achieve eigenvalues using the transition matrix. For LTV systems with period T_{m}, the state vector evolution from the time instantt to t + T_{m} is given by (24)
where is the state transition matrix that transfers the state vector Ψ from t_{1} to t_{2} [53]. The eigenvalue problem is (25)
where is a twobytwo identity matrix and λ represents an eigenvalue. The eigenvalues are found by solving the characteristic equation , leading to (26)
where Φ_{12} and Φ_{11} are elements of the twobytwo matrix . For the illustrated circuit in Figure 6a, the eigenvalues are λ_{p} = e^{j2πfp Tm}, with p = 1, 2, where f_{p} are the two resonance frequencies, with all f_{p} ± nf_{m} harmonics (n is the integer number with modulation frequency f_{m} = 1/T_{m}).
We now demonstrate the degeneracy in an LTV−LC tank shown in Figure 6a. The capacitance C (t) varies between two values C_{1} = 1.5C_{0} and C_{2} = 0.5C_{0} with period T_{m}, where C_{0} = 20 nF.
Defining the state vector Ψ (t) = [q (t), i (t)] ^{T} with capacitor’s charge q (t) and the inductance current i (t), we find the system matrix as (27)
The resonant frequencies versus modulation frequency f_{m} are shown in Figures 6b and 6c. We restrict the plot to frequencies with positive real value, in the range of 0 < f/f_{m} < 1, which could be identified as the fundamental Brillouin Zone (BZ) in a timevarying system. EPD happens at f_{m,e} = 71.72 kHz and f_{m,e} = 63.95 kHz, where the subscript e denotes the corresponding value at the EPD. At an EPD, two eigenvectors and eigenvalues collide, corresponding to a nondiagonalizable transition matrix with a degenerate eigenvalue λ_{e} which is related to the resonance frequency f_{e}. In this configuration, two scenarios may happen to have the EPD (i.e., when the state transition matrix is equivalent to a secondorder JordanBlock matrix). First, when the degenerate eigenvalue is λ_{e} = − 1, which is related to a resonance frequency f_{e} = f_{m}/2, and due to the time periodicity, it also happens at harmonics (f_{e} = f_{m}/2 ± nf_{m}). Second, when λ_{e} = 1, which is related to f_{e} = 0 and to the harmonics f_{e} = ± nf_{m}. Note that here we assume a lossless LC tank (besides the energy injection due to time variation), and we consider the loss effects later on. For modulation values such that 63.95 kHz < f_{m} < 71.72 kHz, the system has two real resonance frequencies, whereas for modulation frequency such that 71.72 kHz < f_{m} < 74.13 kHz, the system experiences complex resonance frequencies, which cause instability and oscillation (rising signal associated with the resonance frequency’s negative imaginary part). Here, we tend to work on the stable part, which has two different real resonance frequencies. At the EPD frequency associated with a modulation frequency of f_{m} = 63.95 kHz, the capacitor’s voltage grows linearly, considering the initial condition of v_{c} (0^{−}) = 50 mV, which indicates that two eigenfrequencies have coalesced.
Fig. 6 (a) Time periodic capacitor connected to the inductor in parallel. The capacitor varies between two values C_{1} and C_{2} with T_{m} period, as shown as an inset. (b) Real and (c) imaginary parts of resonanse frequency evolution varying modulation frequency f_{m}. (b) The time domain signal revealing the second order EPD due to the capacitor’s voltage linear growth, with initial condition of v_{c} (0^{−}) = 50mV, and f_{m} = 63.95kHz. 
5.1 Loss effects on LTV circuit
We shall consider and study loss effects. To validate the occurrence of the EPDs in temporally LC resonator with losses, we assume an LC tank where the constant conductance G is associated to the losses or gain is connected in parallel, as shown in Figure 7a. In this circuit, the capacitance C changes between twolevel capacitance (C_{1} = 1.5C_{0} and C_{2} = 0.5C_{0}) with period T_{m}, as in the previous lossless case. We consider two scenarios where we connect the system to a loss (G > 0) or gain (G < 0) element. Defining the state vector Ψ (t) = [q (t), i (t)] ^{T} with capacitor’s charge q (t) and the inductance current i (t), we find the system matrix as (28)
Fig. 7 (a) Time periodic capacitor C (t) connected to the inductor L_{0} and associated loss/gain Gin parallel. The capacitor varies between two values C_{1} and C_{2} with T_{m} period, as shown in the subset. (b) Real and (c) imaginary parts of resonance frequency evolution varying modulation frequency f_{m} where the LC tank is connected to the lossy conductance G = 1mS. (d) The time domain signal of the second order EPD which is indicated the decaying signal associated to the positive imaginary part of eigenfrequencies. (e) Real and (f) imaginary parts of resonance frequency evolution varying modulation frequency f_{m} where the LC tank is connected to the negative conductance G = − 1mS. (g) Capacitor’s voltage obtained from Keysight ADS circuit simulator, which shows the rising signal. In all time domain simulation, the capacitor has an initial condition of v_{c} (0^{−}) = 50mV. 
5.1.1 Timevarying capacitor: lossy case (G > 0)
The eigenfrequencies’ dispersion diagram in Figures 7b and 7c show real and imaginary parts of the eigenfrequencies versus modulation frequency shown. The system parameters are the same as those as in the previous section: L_{0} = 33 μH, C_{0} = 20 nF, G = 1 mS. Figure 7d shows the capacitor’s voltage at the EPD associated with f_{m} = 63.95 kHz due to the initial condition of v_{c} (0^{−}) = 50 mV obtained from Keysight ADS timedomain circuit simulator. In this configuration, the LTV circuit with loss shows that the imaginary part of the eigenfrequency at every EPD is positive. The state vector, which contains the capacitor’s voltage and inductor current, is proportional to the e^{jωt}. An eigenfrequency with a positive imaginary part leads to a decaying signal (exponential decay in system state vector), as shown for the case in Figure 7d.
5.1.2 Timevarying capacitor: gain case (G < 0)
Figures 7e and 7f show the real and imaginary parts of the system’s eigenfrequencies where the system’s parameters are selected as L_{0} = 33 μH, C_{0} = 20 nF, G = − 1 mS. The negative conductance G here shows the gain which results in EPD frequencies with a negative imaginary part. Thus, the state vector is rising exponentially and making the system unstable. Figure 7g shows the capacitor’s voltage at an EPD associated to f_{m} = 63.95 kHz due to the initial condition of v_{c} (0^{−}) = 50 mV obtained from Keysight ADS circuit simulator. The signal shows a rising behavior, which makes the system unstable and oscillating.
5.2 Timevarying conductance
In this section, we consider the timevarying loss/gain element with the time periodicity of T_{m} in the LTV circuit in Figure 8a. We show the occurrence of EPDs. The parallel conductance is set to G_{1} for half period, and to G_{2} for the other half. Generally, G_{1} and G_{2} could be positive and negative values acting as loss or gain in this scheme. By defining the state vector Ψ (t) = [q (t), i (t)] ^{T} with capacitor’s charge q (t) and the inductance current i (t), we find the system matrix as (29)
We have three scenarios where the time average of the conductance, (30)
is positive (lossy), negative (gain), and zero. The average Avg (G) > 0 means that loss is dominant, whereas Avg (G) < 0 means that gain is dominant, and when Avg (G) = 0 the system has balanced gain and loss. The solution for the system’s state vector Ψ (contains the solution of the charge on the capacitor and current on the inductor) is proportional to the eigenfrequency as Ψ ∝ e^{jωt}. The signal at an EPD, i.e., when the system experiences the coalescence of the eigenvalues at a real ω_{e}, voltages and currents grow linearly as te^{jωt}. This is due to the double pole in the Laplace transform of a signal of a system at the EPD. Moreover, when the eigenfrequencies are complex, signals in the system (currents and voltages) experience exponential growth or decay. In this section, the value of components in the LC tank is set as L_{0} = 33 μH and C_{0} = 20 nF.
Fig. 8 (a) Circuit scheme including LC tank connected to the timevarying loss. (b) real and (c) imaginary parts of eigenfrequency versus modulation frequency f_{m}, where the loss average is zero. (d) Timedomain signal capacitance’s voltage v_{c} (t), which shows the linear growth at EPD. (e) real and (f) imaginary parts of eigenfrequency versus modulation frequency f_{m}, where the loss average is positive (G works as a lossy component in average) (g) Timedomain signal capacitance voltage v_{c} (t) is decaying related to the positive imaginary part of eigenfrequency at EPD. (h) real and (i) imaginary parts of eigenfrequency versus modulation frequency f_{m}, where the loss average is negative (G works as a gain component in average) (j) Timedomain signal capacitance voltage v_{1} (t) is rising related to the negative imaginary part of eigenfrequency at EPD. In all timedomain simulation the capacitor has an initial condition of v_{c} (0^{−}) = 50mV. 
5.2.1 Timevarying conductance: zero average Avg (G) = 0
Figures 8b and 8c show the complex dispersion diagram, eigenfrequencies versus frequency modulation, with a zero timeaverage conductance. The conductance for half a period T_{m}/2 is G_{1} = 4 mS while for the other T_{m}/2 is G_{2} = − 4 mS. The EPDs eigenfrequencies are realvalued here, and the dispersion diagram looks the same as the dispersion diagram of the lossless system. To validate it, we calculate the determinant of the transition matrix as (31)
Thus, for zero timeaverage conductance (Avg (G) = 0) we have which leads to λ_{e} = ± 1. Thus, under the mentioned conditions, the system has a realvalued EPD frequency f_{e}. Figure 8d shows the timedomain signal obtained from Keysight ADS that indicates the secondorder EPD behavior associated to f_{m} = 56.2 kHz, which exhibits a linear growth of the capacitor’s voltage as v_{c} ∝ t cos(ωt + θ).
5.2.2 Timevarying conductance: positive average Avg (G) > 0
The real and imaginary parts of the eigenfrequencies for positive time average conductance (Avg (G) > 0) are shown in Figures 8e and 8f. The system’s parameters are, G_{1} = 4mS and G_{2} = − 2mS, hence the loss is dominant in the system. The EPDs are complex frequencies with a positive imaginary part, which is associated with a damping signal in the circuit, as exhibited in Figure 8g, where modulation frequency is f_{m} = 56.2 kHz.
5.2.3 Timevarying conductance: negative average Avg (G) < 0
Figures 8h and 8i show the real and imaginary parts of the eigenfrequencies for negative time average conductance (Avg (G) < 0). The conductances are G_{1} = 2 mS and G_{2} = − 4 mS, hence the gain is dominant in the system. The EPDs frequencies have a negative imaginary part corresponding to an exponential rise of the signal making the system unstable, as shown in Figure 8j, where modulation frequency is f_{m} = 56.2 kHz.
6 Conclusion
We considered three different circuit configurations supporting an EPD of order two: gyratorbased, PTsymmetry based, and linear timevarying systems. All the configurations exhibit ultrasensitive responses to perturbations, though their operational regimes differ. Each design has some advantages compared to the others. For example, in a gyratorbased circuit, the system has purely real perturbed eigenfrequencies when perturbing one component (e.g., a capacitor), while negative capacitance and inductance are needed to realize such a circuit, which require active components. Small losses or gains in this circuit cause instability. Though it seems to be a complication at first sight and may require working in the transient regime before reaching saturation, instability offers the possibility to work in the unstable oscillatory regime. An EPD is also present in two coupled resonators with balanced gain and loss, i.e., satisfying PT symmetry. The presence of gain in one side of the circuit requires active electronic components to provide a precise gain value. To achieve the sensitive feature of the EPD, a component is varied in a standard sensing scheme (e.g., a capacitor), and the perturbed resonant frequency is detected. However, the sensing scheme proposed in the PTsymmetry regime required tuning the capacitance on the nonsensing part of the circuit to keep the circuit under PTsymmetry while the sensing component is varied [16]. This was done to keep the two shifted frequencies realvalued and avoid instabilities. (However, in a sensing scheme, the value of the varied component is usually the one to be measured; hence it is not known a priori). This complication is not needed in the circuit using the gyrator, and it is also not needed in the circuit based on the single LTV resonator.
Compared to the conventional PTsymmetric circuit where people observed shifted resonance frequencies [16,18], here instead, we have shown that we can work in an oscillatory regime generated by the instability and the nonlinear behavior of the circuit. In other words, we have turned the instability due to broken PTsymmetry (due to a perturbation) to our advantage. The oscillation frequency after reaching saturation is very sensitive to perturbations. Still, it shows the squarerootlike dependency with respect to perturbations and the possibility to measure both signs of an element perturbation (this is not possible with the other schemes in the linear regime discussed in this paper). This nonlinear oscillator scheme is also robust in terms of bringing the system near the EPD, independent of the amount of (nonlinear) small gain used.
Finally, EPDs are found in linear timevarying single resonator circuits where a timevarying capacitor is connected to an inductor. There is no need for gain and lossy elements, though the time modulation requires active components. The effect of additional loss and gain has been discussed for this circuit. We have also shown that timevarying gain or loss connected to a stationary LC tank is another method that leads to EPDs. The simple tuning procedure is one important advantage of the LTV circuit compared to gyratorbased and PTsymmetric circuits. In PTsymmetric and gyratorbased circuits, a tuning process is needed to obtain an EPD, e.g., variable capacitors, gain, or resistors are deemed necessary. On the other hand, in the LTV circuit, the EPD is found by simply changing the modulation frequency, which is done easily in electronics. Note that to obtain the system’s resonance frequency, we need a long enough time interval to measure the signal to derive its frequency response with good accuracy. Thus, the rise/fall time of the signal is important. In this paper, we used the practical quantities as modulation frequency discussed in [23]. There are two work regimes for the system to study the resonance frequencies evolution. First, the circuit needs to be reset periodically, like the circuit in [21], and resonance frequencies are found by looking at the transient response. Second, let the system saturate and then study the system in an oscillatory regime. Note that the system must be designed properly to have a rational rise/fall time in signals. At the same time, fast rising or damping signals cause losing the frequency response resolution.
Acknowledgments
This material is based upon work supported by the National Science Foundation (NSF) awards ECCS1711975 and by the Air Force Office of Scientific Research Award No. FA95501910103.
Appendix A Puiseux fractional power series expansion
The sensitivity of a system due to the perturbation of a system’s component is detected by measuring changes in the system’s observables, like the system’s resonance frequency. The Puiseux fractional power series expansion helps us find the eigenvalues ω_{p} related to the perturbations when working at an EPD of order p. We consider a small perturbation Δ_{X} of a system parameter X as (A1)
where X_{e} is the parameters’ value at the EPD, and X is the parameter’s value after applying perturbation.
Away from an EPD of order p, the system matrix is diagonalizable, and there are p independent eigenvectors. Whereas at an EPD, is nondiagonalizable, and the system has only one eigenvector and p − 1 generalized eigenvectors, and p repeated eigenvalues (i.e., eigenfrequencies). Therefore, the system matrix is similar to a matrix that contains a p × p dimension Jordan block.
For a system with the characteristic equation of under the condition at the EPD, the matrix is made of a p × p Jordan block. Thus, the perturbed eigenfrequencies could be expressed by the Puiseux fractional power series expansion, including powers of such as (A2)
where and the series is a kind of convergent Taylor series of power of . Note that the eigenvalues follow the Puiseux fractional power series expansion at and very near the EPD frequency is a way to validate the existence of an EPD (bifurcation of the eigenvalues). Equation (A2) shows that the eigenfrequency shift from an EPD ω (Δ_{X}) − ω_{e} is proportional to for small Δ_{X}. For second order EPDs, a perturbation Δ_{X} results in the perturbed eigenvalues ω_{p} (Δ_{X}) with p = 1, 2, and the Puiseux fractional power series expansion of ω_{p} (Δ_{X}) is given by [3,54] (A3)
The first two coefficients are expressed (A4) (A5)
where . The coefficients are calculated at the EPD, i.e., at Δ_{X} = 0 and ω = ω_{e}.
References
 M.I. Vishik, L.A. Lyusternik, The solution of some perturbation problems for matrices and selfadjoint or nonselfadjoint differential equations I, Russian Math. Surv. 15, 1 (1960) [CrossRef] [Google Scholar]
 P. Lancaster, On eigenvalues of matrices dependent on a parameter, Numerisc. Math. 6, 377 (1964) [CrossRef] [Google Scholar]
 T. Kato, Perturbation Theory for Linear Operators (SpringerVerlag New York Inc., New York, 1966) [Google Scholar]
 W.D. Heiss, A.L. Sannino, Avoided level crossing and exceptional points, J. Phys. A: Math. General 23, 1167 (1990) [CrossRef] [Google Scholar]
 W.H. Steeb, W.D. Heiss, Energy eigenvalue levels and singular point analysis, Phys. Lett. A 152, 339 (1991) [CrossRef] [Google Scholar]
 A.P. Seyranian, Sensitivity analysis of multiple eigenvalues, J. Struct. Mech. 21, 261 (1993) [Google Scholar]
 C.M. Bender, S. Boettcher, Real spectra in nonHermitian Hamiltonians having PT symmetry, Phys. Rev. Lett. 80, 5243 (1998) [Google Scholar]
 W.D. Heiss, Repulsion of resonance states and exceptional points, Phys. Rev. E 61, 929 (2000) [CrossRef] [Google Scholar]
 C.M. Bender, M.V. Berry, A. Mandilara, Generalized PT symmetry and real spectra, J. Phys. A 35, L467 (2002) [CrossRef] [Google Scholar]
 A. Figotin, I. Vitebskiy, Oblique frozen modes in periodic layered media, Phys. Rev. E 68, 036609 (2003) [Google Scholar]
 A. Figotin, I. Vitebskiy, Gigantic transmission bandedge resonance in periodic stacks of anisotropic layers, Phys. Rev. E 72, 036619 (2005) [Google Scholar]
 A. Guo et al., Observation of PTsymmetry breaking in complex optical potentials, Phys. Rev. Lett. 103, 093902 (2009) [CrossRef] [Google Scholar]
 W.D. Heiss, The physics of exceptional points, J. Phys. A 45, 444016 (2012) [Google Scholar]
 B. Peng et al., Parity–timesymmetric whisperinggallery microcavities, Nat. Phys. 10, 394 (2014) [Google Scholar]
 F. Alex, V. Ilya, Slow light in photonic crystals, Waves Random Complex Media 16, 293 (2006) [CrossRef] [Google Scholar]
 F. Alex, V. Ilya, Frozen light in photonic crystals with degenerate band edge, Phys. Rev. E. 74, 066613 (2006) [CrossRef] [Google Scholar]
 M.V. Berry, Physics of nonhermitian degeneracies, Czechoslovak J. Phys. 54, 1039 (2004) [Google Scholar]
 P.Y. Chen et al., Generalized parity–time symmetry condition for enhanced sensor telemetry, Nat. Electr. 1, 297 (2018) [CrossRef] [Google Scholar]
 M. Sakhdari, M. Farhat, P.Y. Chen, PTsymmetric metasurfaces: wave manipulation and sensing using singular points, New J. Phys. 19, 065002 (2017) [CrossRef] [Google Scholar]
 J. Schindler, A. Li, M.C. Zheng, F.M. Ellis, T. Kottos, Experimental study of active LRC circuits with PT symmetries, Phys. Rev. A 84, 040101 (2011) [Google Scholar]
 H. Hodaei, M.A. Miri, M. Heinrich, D.N. Christodoulides, M. Khajavikhan, Paritytimesymmetric microring lasers, Science 346, 975 (2014) [CrossRef] [Google Scholar]
 T. Stehmann, W.D. Heiss, F.G. Scholtz, Observation of exceptional points in electronic circuits, J. Phys. A: Math. General 37, 7813 (2004) [CrossRef] [Google Scholar]
 K. Rouhi, A. Nikzamir, A. Figotin, F. Capolino, Exceptional point in a degenerate system made of a gyrator and two unstable resonators, Phys. Rev. A 105, 032214 (2022) [CrossRef] [Google Scholar]
 A. Nikzamir, K. Rouhi, A. Figotin, F. Capolino, Demonstration of exceptional points of degeneracy in gyratorbased circuit for highsensitivity applications., arXiv preprint arXiv:2107.00639 (2021) [Google Scholar]
 H. Kazemi, M.Y. Nada, A. Nikzamir, F. Maddaleno, F. Capolino, Experimental demonstration of exceptional points of degeneracy in linear time periodic systems and exceptional sensitivity, J. Appl. Phys. 131, 144502 (2022) [Google Scholar]
 H. Kazemi, M.Y. Nada, T. Mealy, A.F. Abdelshafy, F. Capolino, Exceptional points of degeneracy induced by linear timeperiodic variation, Phys. Rev. Appl. 11, 014007 (2019) [Google Scholar]
 H. Kazemi et al., Ultrasensitive radio frequency biosensor at an exceptional point of degeneracy induced by time modulation, IEEE Sensors J. 21 (2021) [Google Scholar]
 A. Figotin, I. Vitebskiy, Oblique frozen modes in periodic layered media, Phys. Rev. E 68, 036609 (2003) [Google Scholar]
 M.Y. Nada, M.A.K. Othman, F. Capolino, Theory of coupled resonator optical waveguides exhibiting highorder exceptional points of degeneracy, Phys. Rev. B 96, 184304 (2017) [Google Scholar]
 A.F. Abdelshafy, T. Mealy, E. Hafezi, A. Nikzamir, F. Capolino, Exceptional degeneracy in a waveguide periodically loaded with discrete gain and radiation loss elements, Appl. Phys. Lett. 118, 224102 (2021) [CrossRef] [Google Scholar]
 P.J. Chen, S. Saati, R. Varma, M.S. Humayun, Y.C. Tai, Wireless intraocular pressure sensing using microfabricated minimally invasive flexiblecoiled LC sensor implant, J. Microelectromech. Syst. 19, 721 (2010) [CrossRef] [Google Scholar]
 T.Q. Trung, S. Ramasundaram, B.U. Hwang, N.E. Lee, An allelastomeric transparent and stretchable temperature sensor for bodyattachable wearable electronics, Adv. Mater. 28, 502 (2016) [CrossRef] [Google Scholar]
 Y. Feng, L. Xie, Q. Chen, L.R. Zheng, Lowcost printed chipless RFID humidity sensor tag for intelligent packaging, IEEE Sensors J. 15, 3201 (2015) [CrossRef] [Google Scholar]
 Y. Qian, N. Sahar, D. Prativa, N. Sangjun, H. Emily, K. Taeil, E. Rahim, All3DPrinted, Flexible, and Hybrid Wearable Bioelectronic Tactile Sensors Using Biocompatible Nanocomposites for Health Monitoring, Adv. Mater. Technol. 7, 2101034 (2021) [Google Scholar]
 K. Rouhi, R. Marosi, T. Mealy, A.F. Abdelshafy, A. Figotin, F. Capolino, Exceptional degeneracies in traveling wave tubes with dispersive slowwave structure including spacecharge effect, Appl. Phys. Lett. 118, 263506 (2021) [Google Scholar]
 P.J. Chen, D.C. Rodger, S. Saati, M.S. Humayun, Y.C. Tai, Microfabricated implantable parylenebased wireless passive intraocular pressure sensors, J. Microelectromech. Syst. 17, 1342 (2008) [CrossRef] [Google Scholar]
 S.R. Corrie, J.W. Coffey, J. Islam, K.A. Markey, M.A.F. Kendall, Blood, sweat, and tears: developing clinically relevant protein biosensors for integrated body fluid analysis, The Analyst 140, 4350 (2015) [CrossRef] [Google Scholar]
 P. Tseng, B. Napier, L. Garbarini, D.L. Kaplan, F.G. Omenetto, Functional, RF‐trilayer sensors for tooth‐mounted, wireless monitoring of the oral cavity and food consumption, Adv. Mater. 30, 1703257 (2018) [CrossRef] [Google Scholar]
 Y.J. Zhang et al., Noninvasive glucose sensor based on paritytime symmetry, Phys. Rev. Appl. 11, 044049 (2019) [CrossRef] [Google Scholar]
 Z.P. Liu et al., Metrology with PTsymmetric cavities: enhanced sensitivity near the PTphase transition, Phys. Rev. Lett. 117, 110802 (2016) [CrossRef] [Google Scholar]
 J. Wiersig, Enhancing the sensitivity of frequency and energy splitting detection by using exceptional points: application to microcavity sensors for singleparticle detection, Phys. Rev. Lett. 112, 203901 (2014) [CrossRef] [Google Scholar]
 J. Wiersig, Sensors operating at exceptional points: general theory, Phys. Rev. A 93, 033809 (2016) [Google Scholar]
 M.A. Miri, A. Alù, Exceptional points in optics and photonics, Science 363, 6422 (2019) [Google Scholar]
 C.E. Rüter, K.G. Makris, R. ElGanainy, D.N. Christodoulides, M. Segev, D. Kip, Observation of paritytime symmetry in optics, Nat. Phys. 6, 192 (2010) [Google Scholar]
 L. Feng, Z.J. Wong, R.M. Ma, Y. Wang, X. Zhang, Singlemode laser by paritytime symmetry breaking, Science 346, 972 (2014) [CrossRef] [Google Scholar]
 A. Figotin, Perturbations of circuit evolution matrices with Jordan blocks, J. Math. Phys. 62, 042703 (2021) [Google Scholar]
 A. Figotin, Synthesis of lossless electric circuits based on prescribed Jordan forms, J. Math. Phys. 61, 122703 (2020) [Google Scholar]
 B.D.H. Tellegen, The gyrator, a new electric network element, Philips Res. Rep. 3, 81 (1948) [Google Scholar]
 A. Antoniou, Gyrators using operational amplifiers, Electr. Lett. 3, 350 (1967) [Google Scholar]
 M. Ehsani, I. Husain, M.O. Bilgic, Power converters as natural gyrators, IEEE Trans. Circ. Syst I 40, 946 (1993) [CrossRef] [Google Scholar]
 W.H. Holmes, S. Gruetzmann, W.E. Heinlein, Highperformance directcoupled gyrators, Electr. Lett. 3, 45 (1967) [Google Scholar]
 A. Nikzamir, K. Rouhi, A. Figotin, F. Capolino, Exceptional points of degeneracy in gyratorbased coupled resonator circuit, in 2021 Fifteenth International Congress on Artificial Materials for Novel Wave Phenomena (Metamaterials) (2021), pp. 302–304 [Google Scholar]
 D. Oshmarin, A.F. Abdelshafy, A. Nikzamir, M.M. Green, F. Capolino, Experimental demonstration of a new oscillator concept based on degenerate band edge in microstrip circuit, arXiv preprint arXiv:2109.07002 (2021) [Google Scholar]
 S. Kananian, G. Alexopoulos, A.S.Y. Poon, Couplingindependent realtime wireless resistive sensing through nonlinear PT symmetry, Phys. Rev. Appl. 14, 064072 (2020) [CrossRef] [Google Scholar]
 S. Ramezanpour, A. Bogdanov, Tuning exceptional points with Kerr nonlinearity, Phys. Rev. A 103, 043510 (2021) [CrossRef] [Google Scholar]
 J.A. Richards, Analysis of Periodically TimeVarying Systems (Springer, Berlin Heidelberg, 1983) [CrossRef] [Google Scholar]
 A. Welters, On explicit recursive formulas in the spectral perturbation analysis of a Jordan block, SIAM J. Matrix Anal. Appl. 32, 1 (2011) [Google Scholar]
Cite this article as: Alireza Nikzamir, Kasra Rouhi, Alexander Figotin, Filippo Capolino, How to achieve exceptional points in coupled resonators using a gyrator or PTsymmetry, and in a timemodulated single resonator: high sensitivity to perturbations, EPJ Appl. Metamat. 9, 14 (2022)
All Figures
Fig. 1 Three different configurations where support EPD (a) Two resonators coupled through gyrator (b) Two mutually coupledresonator terminated with balanced gain and loss (PTsymmetry circuit) (c) Timevarying circuit which just one component in the circuit needed to change over time periodically. The timevarying component here could be capacitor or inductor or loss. 

In the text 
Fig. 2 (a) The gyratorbased circuit with the ideal gyrator connected two parallel LC tanks. (b) Opampbased circuit configuration to obtain negative inductance and capacitance. The sensitivity of the (c), (e), real and (d), (f), imaginary parts of the eigenfrequencies to (c), (d) gyration resistance, (e), (f) positive capacitance C_{1}. Solid lines: solution of eigenvalue problem of Eq. (2); greendashed lines: Puiseux series approximation truncated to its second term. Voltage of positive capacitance v_{1} (t) at EPD (g) timedomain, and (h) frequencydomain. The frequencydomain result is calculated by applying an FFT with 10^{6} samples in the time window of 0 to 3 ms. (i) Root locus of zeros of Y_{total} (ω) = 0 showing the real versus imaginary parts of resonance frequencies by perturbing gyration resistance. At the EPD, two zero collide at ω_{e} and the system’s total admittance has the form of Y_{total} (ω) ∝ (ω − ω_{e}) ^{2}. 

In the text 
Fig. 3 (a) Schematic view of the lossy gyratorbased circuit, with a resistor in each resonator. Variation of (b) real and (c) imaginary parts of the angular eigenfrequencies to a resistor perturbation on the left resonator, i.e., when γ_{1} changes and γ_{2} = 0. Variation of (d) real and (e) imaginary parts of the angular eigenfrequencies to a resistor perturbation on the right resonator, i.e., when −γ_{2} changes and γ_{1} = 0. 

In the text 
Fig. 4 (a). Two mutually coupled LC tanks terminated with balanced gain on the left and loss on the right side G_{1} = G_{2} = G (b) real and (c) imaginary parts of evaluated eigenfrequencies by perturbing G. In the illustrated plots, only eigenvalues with the positive real parts are drawn. (d) The coalescence factor corresponding to calculated eigenvectors which showing the degeneracy of two eigenvectors exactly at G_{e}(e) Root locus of zeros of Y_{in} (ω) − G = 0 showing the real and imaginary parts of resonance frequencies of the circuit when perturbing both load resistance and gain G. At the EPD, the system’s total admittance is Y_{in} (ω) − G_{e} ∝ (ω − ω_{e}) ^{2}; hence it shows a double zero at ω_{e}. 

In the text 
Fig. 5 (a). Two mutually coupled LC tanks terminated with nonlinear gain on the left and linear loss resistance on the right side where always G_{1} = 1.001 × G_{2}. (b) Time domain response. (c) Frequency response of the gainside capacitor voltage at the EPD point. (d) Real and (e) imaginary parts of the eigenfrequencies versus load capacitance perturbation ΔC_{2}; solid line: result from ; green dashed line: Puiseux fractional power series expansion terminated to its second order; black dots are obtained from the nonlinear time domain simulation. (f) Oscillation frequency versus load capacitance perturbation for three different mismatches between gain and loss δ = 0, 0.001, 0.01: the three different gains provide the same saturated oscillation frequencies. 

In the text 
Fig. 6 (a) Time periodic capacitor connected to the inductor in parallel. The capacitor varies between two values C_{1} and C_{2} with T_{m} period, as shown as an inset. (b) Real and (c) imaginary parts of resonanse frequency evolution varying modulation frequency f_{m}. (b) The time domain signal revealing the second order EPD due to the capacitor’s voltage linear growth, with initial condition of v_{c} (0^{−}) = 50mV, and f_{m} = 63.95kHz. 

In the text 
Fig. 7 (a) Time periodic capacitor C (t) connected to the inductor L_{0} and associated loss/gain Gin parallel. The capacitor varies between two values C_{1} and C_{2} with T_{m} period, as shown in the subset. (b) Real and (c) imaginary parts of resonance frequency evolution varying modulation frequency f_{m} where the LC tank is connected to the lossy conductance G = 1mS. (d) The time domain signal of the second order EPD which is indicated the decaying signal associated to the positive imaginary part of eigenfrequencies. (e) Real and (f) imaginary parts of resonance frequency evolution varying modulation frequency f_{m} where the LC tank is connected to the negative conductance G = − 1mS. (g) Capacitor’s voltage obtained from Keysight ADS circuit simulator, which shows the rising signal. In all time domain simulation, the capacitor has an initial condition of v_{c} (0^{−}) = 50mV. 

In the text 
Fig. 8 (a) Circuit scheme including LC tank connected to the timevarying loss. (b) real and (c) imaginary parts of eigenfrequency versus modulation frequency f_{m}, where the loss average is zero. (d) Timedomain signal capacitance’s voltage v_{c} (t), which shows the linear growth at EPD. (e) real and (f) imaginary parts of eigenfrequency versus modulation frequency f_{m}, where the loss average is positive (G works as a lossy component in average) (g) Timedomain signal capacitance voltage v_{c} (t) is decaying related to the positive imaginary part of eigenfrequency at EPD. (h) real and (i) imaginary parts of eigenfrequency versus modulation frequency f_{m}, where the loss average is negative (G works as a gain component in average) (j) Timedomain signal capacitance voltage v_{1} (t) is rising related to the negative imaginary part of eigenfrequency at EPD. In all timedomain simulation the capacitor has an initial condition of v_{c} (0^{−}) = 50mV. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.