Tunable optical amplification and group delay in cavity magnomechanics
Abstract
In this work, we theoretically investigate the controllable output probe transmission and group delay in a hybrid cavity magnomechanics (CMM) system. The setup comprises a gain (active) cavity and a passive (loss) cavity, which incorporates an optical parametric amplifier (OPA) and two yttrium iron garnet spheres to facilitate magnon-photon coupling. Unlike the single transparency window typically resulting from magnon-photon interactions, we also observe magnomechanically induced transparency due to nonlinear magnon-phonon interactions. Additionally, two absorption dips on either side of the central absorption dip can be asymmetrically modulated into amplification and absorption by varying different system parameters. A -symmetric to broken--symmetric phase transition is observed in both balanced and unbalanced gain-to-loss scenarios. Notably, replacing the second passive cavity with an active one mitigates high absorption and introduces effective gain into the system. Our findings reveal that the group delay of the probe light can be adjusted between positive and negative values by modifying various system parameters. This study provides a robust platform for controlling light propagation in CMM systems, highlighting potential applications in optical communication and signal processing.
I Introduction
Cavity magnomechanics (CMM) [1, 2, 3] has become a rapidly growing field due to its applications in modern quantum technologies [4, 5, 6, 7]. In particular, a new type of hybrid system, primarily relying on collective spin excitations of yttrium iron garnet (YIG) has garnered significant attention and established a platform for quantum information science owing to its robust spin-spin exchange interactions, high spin density, as well as low dissipation rates [8, 9]. Benefiting from the dynamics of YIG results in promising applications, such as spin current manipulation [10, 11], long-time memory [12, 13], microwave-to-optical conversion [14, 15], manipulation of magnons and photons via exceptional points [16, 17, 18], magnon-induced nonreciprocity [19, 20], enhanced tripartite interactions [21], quantum entanglement [22, 23], and precision measurements [24, 25], etc.
One of the fundamental principles in quantum mechanics is the Hermiticity of physical observables, particularly the Hamiltonian operator. This principle ensures that the eigenvalues are real and that probability is conserved. Interestingly, a broad range of non-Hermitian Hamiltonians can exhibit completely real spectra [26, 27]. This is especially true for Hamiltonians in parity-time ()-symmetric systems. -symmetric systems with gain/loss have received considerable interest in the field of quantum optics [28, 29], and shown great potential in optomechanical systems [30, 31]. Remarkably, such systems possess a phase transition [32] from the -symmetric phase to the broken--symmetric phase, which will take place at the exceptional point (EP), where the eigenvalues and the associated eigenvectors coalesce and some intriguing physical phenomena may arise. Optical and quantum optical systems are natural test beds for -symmetric systems [33, 34], which has resulted in showing several exotic phenomena, such as unidirectional invisibility [35, 36], nonreciprocal light propagation [37], optical light stopping, and sensing [38].
New attempts have been made to investigate non-Hermitian physics in magnomechanic quantum systems. The second-order exceptional point is found in a two-mode cavity-magnonic system in which the cavity mode’s gain is produced by coherent perfect absorption [39]. Experiments have been carried out to investigate the engineering of anti--symmetry, focusing on the adiabatic elimination of the cavity field and the dissipative coupling between two magnon modes [40]. In addition to their distinctive spectral responses, these -symmetric systems can control the transmission of the output microwave field [41, 42]. The essential technique underlying such an application ismagnomechanically induced transparency (MMIT) [43, 44], where a narrow spectral hole appears inside the absorption spectrum due to significant magnon-photon interaction. Further research has shown that weak magnon-phonon coupling can result in multiple transmission windows by altering the phase and amplitude of the applied magnetic field [45]. Meanwhile, cavity magnonic system was also utilized to investigate novel phenomena and applications [46, 47], such as, nonreciprocal double-carrier frequency combs [48], bistability or multistability of magnon–polariton [49], magnon-induced chaos [50], and so on.
Optomechanically induced transparency (OMIT) is well-known for its application in investigating the slowing and storing of light within a cavity [51, 52]. Leveraging -symmetry in optomechanical systems improves control over light transmission [53, 54], and enables for light conversion from subluminal to superluminal. However, their work may encounter experimental difficulties because the gain in the adjunct cavity could cause instability in the entire system [55]. The CMM system has significant advantages over optomechanical systems, including strong hybridization between the magnon and photon modes, and a substantially longer lifetime for the phonon mode compared to cavity optomechanics systems. Enhancing the phonon lifetime allows for the transmission of wider spatiotemporal pulses with minimal distortion, enabling them to be stored for longer periods as mechanical excitations. Additionally, CMM systems provide enhanced tunability, as the magnon frequency can be adjusted by varying the external magnetic field. Leveraging these advantages, there has recently been increasing interest in achieving gain in various components of a CMM system [56]. This provides a framework for developing a -symmetry CMM system [39].
Furthermore, it has been demonstrated that the combination of nonlinear optics and CMM systems has led to numerous physical phenomena to strengthen quantum effects [57, 58]. Inserting an optical parametric amplifier (OPA) into the cavity causes an optical amplification [59, 60], enhancement of optical nonlinearity [61], and normal mode splitting [62]. Therefore, the interaction between OPA and the cavity can not only alter the dynamic instabilities [63] but also find applications such as cooling enhancement [64], and magnon blockade [65]. Furthermore, by altering the group delay of the output probe light, one may obtain the slow/fast light effects [66, 67].

In this study, leveraging the MMIT effect, we systematically investigate the tunable optical transmission and its associated group delay in a CMM system. The system consists of two optically coupled microcavities: a gain-assisted auxiliary cavity and a passive cavity, which includes an OPA and two YIG spheres driven coherently by a bichromatic laser field. By applying a perturbation approach to solve the Heisenberg-Langevin equations, we derive an analytical solution for the optical transmission. Our calculations demonstrate that, in addition to the single transparency window typically caused by magnon-photon interactions, we also observe MMIT due to nonlinear magnon-phonon interactions. Furthermore, we find that two additional absorption dips, located on either side of the central absorption dip, can be asymmetrically modulated into amplification and absorption by varying different system parameters. We establish that the transition from a broken -symmetry phase to an unbroken -symmetry phase can be achieved through effective electromechanical coupling, considering both cavity loss and gain. Notably, the transmission rate is significantly enhanced in the active-passive configuration compared to the passive-passive scenario, indicating that the gain from the auxiliary cavity contributes to increasing the transmission of the probe light. Finally, by numerically computing the group delay of the probe field, we show that it is straightforward to achieve the transition between slow and fast light. These attributes suggest that CMM devices can serve as powerful tools for controlling light propagation, with potential applications in optical communication. [68], and long-lifetime quantum memories [52].
This work is organized as follows: In §II we present our proposed theoretical model and derive the steady-state solutions of our system. The results of our research are provided in §III. Finally, in §IV we provide a comprehensive conclusion that encapsulates our research findings. To ensure the completeness of our work, we provide detailed calculations of algebraic equations in §A.1.
II Proposed theoretical model
This section discusses the dynamics of our proposed theoretical model, depicted in Fig. 1. Our model comprises two coupled microwave cavities ( and ). Cavity is treated as a passive (loss) cavity, with the loss rate containing a degenerate OPA and two highly polished single-crystal ferromagnetic YIG spheres (with a 250 m diameter as in Ref [69]). We refer to this cavity as a CMM resonator. Applying an external uniform magnetic field to each YIG sphere activates the magnon mode, which interacts with the microwave cavity mode through magnetic dipole interaction [70, 71, 72]. The external bias magnetic field causes deformation of the lattice structure of the YIG spheres, resulting in magnon-phonon interaction, and vice versa [73]. The strength of the single-magnon magnomechanical interaction is exceedingly small [69] and it is influenced by both the direction of the applied bias field and the diameter of the sphere. It is feasible to disregard the magnomechanical interaction of the YIG1 sphere by either employing a bigger sphere or changing the direction of the applied biased magnetic field [74]. Here, we suppose that the bias field on YIG1 is such that the single-magnon magnomechanical interaction is weak enough to be ignored [69]. The active (gain) cavity () with the gain rate is coupled to the passive (loss) cavity via optical tunneling rate . The CMM resonator is driven by two external optical laser fields: a weak probe and a strong pump field. The symbols and denote the amplitudes, and represent the frequencies and , correspond to the phases of the probe and pump fields, respectively. A weak external microwave field with frequency and the phase is directly drives the YIG2 sphere to enhance the magnon-phonon coupling.
In this study, we assume that the diameters of the YIG spheres are significantly smaller than the wavelengths of the applied fields, allowing us to safely neglect the influence of light radiation pressure in our system [75, 76]. We also ignored the nonlinear Kerr term which may occur due to a strongly driven magnon mode [77, 78]. The total Hamiltonian of the system reads:
(1) | ||||
We assume that both cavities in our CMM system have the same resonance frequency () [79]. Applying the rotating-wave approximation to the system, the interaction terms simplify as , by neglecting the fast oscillating terms . This is valid when which is the case to be considered in the present work ( and are the decay rates of the cavity and magnon modes, respectively). To proceed, it is convenient to shift to a frame rotating at the drive frequency . Following the transformation with , the Hamiltonian in Eq. (1) can be rewritten as [80, 81]:
(2) | ||||
where () shows the frequency detuning between cavity (magnon) and control field. stands for the frequency detuning of probe-control field, and shows the detuning of the magnon drive and control fields. indicates the relative phase between the probe and pump fields. , , and are the frequencies that correspond to the resonances of the photon, phonon, and magnon modes, respectively. The term indicates the photon exchange coupling strength interaction between the two cavities. () signifies the magnon-photon (phonon) coupling rates [82]. represents nonlinear gain, while indicates the phase of OPA, respectively. The term represents the coupling strength between the magnon and the microwave field; in which is the total number of spins (with shows the volume, and is the spin density of YIG sphere) [75, 83], = 28 GHz/T represents the gyromagnetic ratio, with shows the field amplitude. The terms signify the interaction between the cavity field and two input fields. The terms represent the amplitude of the strong pump (weak probe) fields, whereas and are the powers of the corresponding pump and probe fields, respectively. The symbol represents the coupling between the CMM resonator and the output port, defined as , where is the external decay rate of the cavity and is the total decay rate of the cavity fields, consisting of both intrinsic and external decay rates () [84, 85]. In this study, we consider the CMM resonator operating in the critical-coupling regime and set throughout the work [86].
II.1 Quantum dynamics and fluctuations:
In our study, we primarily focus on the average response of the system; thus by considering the Hamiltonian of Eq. (2), and ignoring the input quantum noises [87], the Heisenberg-Langevin equations of our compound CMM system can be expressed as:



(3) | ||||
where and denotes the decay rates of magnons and phonons. Keep in mind that , in Eq. (3) is associated to passive-active compound CMM resonator system, while , is associated to a passive-passive system [88].
Relative to the intensity of the controlled field, we assume that the intensities of the external magnetic driving fields and the weak probe field satisfies the condition and [87, 89]. Thus in this case, we can linearize the set of above dynamical equations of our proposed model by assuming each operator is the sum of its mean value and fluctuation i.e., , where [44, 90], the first term represents the steady-state values and depicts the small fluctuating terms. By using perturbation expansion ansatz in Eq. (3) and considering the time derivative equal to zero we can obtain the steady-state solution as:
(4a) | ||||
(4b) | ||||
(4c) | ||||
(4d) | ||||
(4e) |
Here we define: , , . By performing a series of mathematical calculations and neglecting all higher-order terms the evolution of the perturbation terms of Eq. (3) can be shown as:
(5) | ||||
here, , and signifies the effective strength of the magnon-phonon interaction, and effective magnon detuning respectively.
To determine the amplitudes of the first-order side band, we expand the small fluctuation terms in Eq. (5) by using the following ansatz as [91, 87]:
(6) | ||||
We consider the magnon driving field to become resonant with the probe field frequency which leads us to consider [45, 92], and and denotes the th cavity generated probe field and four-wave mixing field amplitudes, respectively.
Now by employing Eq. (6) into Eq. (5) and comparing the coefficients of similar order, one can get the amplitudes of the first ordered side-bands (detail calculations and the explicit constants are presented in §A.1) as [91, 87]: After solving Eqs. (12) we derive the coefficients for the first lower sidebands. This demonstrates the characteristics of our system. Therefore, the coefficient of the , which indicates the output probe field amplitude of the CMM resonator can be written as:
(7) |
Here, corresponds to the relative phase of the applied fields. Furthermore, by employing the input/output relationship of the cavity [93], we get the normalized output probe field intensity from the CMM resonator as:
(8) |
The real and imaginary components of the () can be quantified as the absorption and dispersion of the probe field at the probe frequency, respectively.


III Results
In this part, we provide our numerical findings to investigate the transmission and group delay of the output probe field in a cavity-magnon system by employing a set of experimentally feasible parameter values [69, 82], unless stated differently. We use the following values, i.e., the degenerate microwave cavity frequency GHz, the first cavity decay rate MHz, and the coupling parameter . The values for the YIG sphere are taken as: the diameter m, the spin density . The frequency of the phonon mode MHz, and decay rate Hz. The magnon decay rate MHz, gyromagnetic ratio = 28 GHz/T, and . In addition, we take the power of the pump field mW, and the amplitude of the probe filed , with .
In §III.1 we study the effect of magnon-photon and magnon-phonon coupling strengths on MMIT window profiles. Next, in §III.2 we examine how the amplitude of the microwave field, its relative phase, and the nonlinear parameter OPA impact the transmission rate of the output probe field. Then, in §III.3 we present our analysis to evaluate the controllable phase transition and MMIT in a passive-active CMM system. Furthermore, in §III.4 we represent our results for optimizing group delays in an active-passive compound CMM system. Finally, in §III.5 we briefly discuss the possible experimental realization of our proposed scheme.
III.1 Effect of magnon-photon and magnon-phonon coupling strengths, and OPA on MMIT window profiles
In this subsection, we will explore the physics behind the formation of MMIT window profiles by systematically analyzing the influence of magnon-photon and magnon-phonon coupling strengths in our model. Figure 2 displays the absorption (dispersion) spectrum of the probe field versus probe detuning () for various coupling strengths. In Fig. 2(a), we consider the scenario where both the magnon-phonon coupling () and the coupling strength of the magnon mode () with the cavity are absent. As a result, only magnon mode is coupled with the cavity. Following this scenario, we observe a typical Lorentzian-type absorption peak of the output spectrum (as shown by a blue solid line) in Fig. 2(a), and the corresponding dispersion curve is indicated by red dotted lines. Interestingly, when the coupling strengths of both magnon modes () are present, and only the magnon-phonon coupling () is absent, we observe that the typical Lorentzian-type absorption peak at resonance transforms into a transparency window. Additionally, two symmetrical absorption dips occur at (left) and (right) as shown in Fig. 2(b). The central peak at represents magnon-induced transparency, while the two sideband absorption peaks indicate magnon-induced absorption (MIA) [66]. The physical reason for these peaks is the generation of magnon-polaritons due to the photon-magnon interaction in our system. As a result, our entire system is effectively reduced to a cavity magnonic (CM) system [94, 95, 96].



It is noteworthy that the width of the transparency window and the amplification in our system can be controlled by adjusting the magnon-photon coupling strength (). However, the results are not shown here. Finally, we consider the scenario where all three couplings are present, i.e., , and Hz. In this case, due to the nonzero magnetostrictive interaction, the single MMIT window at the resonance point () of the Fig. 2(b) splits into double MMIT windows with a central absorption peak, as shown in Fig. 2(c). The physical explanation is that the existence of phonon interactions can create additional pathways for destructive interference, leading to additional transparency windows and our entire system reduced to CMM system [69].
In Fig. 3, we plot the real and imaginary parts of the output probe field as a function of optical detuning for different selected OPA gain parameters . It is evident that changes in the parametric gain influence the transparency window. In the absence of OPA (), we observe clear absorption in our system. However, introducing OPA () converts the positive absorption peak at resonance into a negative peak, indicating gain in our system. The gain provided by the OPA in CMM systems offers substantial benefits, including signal amplification, tunability, sensitivity, and the ability to explore complex physical phenomena. Given these advantages, our work provides greater benefits compared to previously reported studies in [44].


, and . All other parameters remains same as in Fig. 2(c) except for and T.




III.2 Role of microwave field’s amplitude, phase, and OPA on transmission profiles
In the following, we investigate the impact of the microwave field’s amplitude, its relative phase, and the nonlinear parameter OPA on the transmission profiles of the output spectrum. Referring to our discussion about Fig. 2, the positions of the magnetomechanically induced absorption and the MIA peaks are determined, respectively. Here, we are interested in employing an extra magnon driving field to examine the detailed impacts of the magnon driving parameters on the features of the output field in our CMM system. In Fig. 4(a), We graph the transmission spectrum () of the output probe field versus optical detuning () for different magnetic field amplitudes (), keeping the relative phase . It is observed that with an increase in the magnetic field amplitude (), the central absorption peak at the resonance point () becomes slightly deeper. There is a slight enhancement in the transparency spectrum on the left side of the resonance point, coupled with a lower absorption dip on the right side of the resonance point (). Interestingly, when we further increase the magnetic field amplitude (T), the symmetric absorption dips on both sides of the resonance point () become asymmetric. In the blue-sideband regime (), the absorption dip not only transforms into transparency but also exceeds an amplitude of 1, thereby enabling the amplification of the weak probe field. Conversely, in the red-sideband regime (), there is a lower absorption dip as shown by the blue solid line in Fig. 4(a). Controlling the magnitude of the magnon driving field can convert MIA into amplification, which has potential uses in quantum sensing and computation.
Notably, the relative phase of the applied fields has been employed to regulate OMIT in cavity optomechanics devices [97, 98]. We now study the effect of the relative phase of the cavity and magnon driving fields on the transmission rate of the system. In Fig. 4(b), we set the magnetic field amplitude T and display the transmission spectrum () of the output probe field versus the probe detuning with different relative phase (). Figure. 4(b) clearly shows that the absorption resulting from the combined action of magnons and polaritons is sensitive to the relative phase (). Setting leads to an asymmetric transmission profile in the output spectrum. This configuration results in amplification on the right side of the resonance point () and absorption on the left side, as illustrated by the red dotted line in Fig. 4(b). Accordingly, when we set the relative the switching of the asymmetry line shapes becomes noticeable [see blue solid line in Fig. 4(b)]. In contrast, the central absorption dip at remains fairly stable as the relative phase varies from 0 to .
The physical mechanism behind this attribute is that by adjusting the relative phase () of the cavity and magnon driving fields, periodic constructive and destructive interference occurs, significantly altering the cavity-field transmission spectra [97, 99]. As shown in Fig. 4(b), at the exact resonance point (), destructive interference occurs, resulting in strong MIA. In contrast, on either side of the resonance point, varying from 0 to leads to constructive interference and noticeable amplification. This behavior is referred to as magnon-induced amplification (MIAMP). This ability of selectively switching and amplifying the input probe signal could be highly desirable in practical optical communications [97, 100].
To better illustrate the impact of the relative phase on the MMIT profile, we plot the contour map of the transmission spectrum against the optical detuning () and the relative phase as shown in Fig. 4(c). The transmission of the probe light changes periodically. At the exact resonance point, we observe strong absorption. The transmission rate transitions from strong absorption to amplification on either side of the resonance point as the phase is tuned from to .
Therefore, based on the results from Fig. 4, we conclude that the application of an additional magnon driving field in the CMM system allows the MIA and the transparency of the CMM system to be switched by adjusting the magnetic field amplitude and the relative phase while keeping the other system parameters fixed.
To illustrate the significant effect of an OPA on the transmission rate () we depict the graph of transmission rate () versus detuning () as shown in Fig. 5. In Fig. 5(a), in the absence of OPA, specifically when the gain () along with phase difference () are both zero, we observed that our system displays an absorption dip at resonance position () as shown by the red dot curve in Fig. 5 (a). As anticipated, introducing the nonlinear gain () results in the MMIT spectrum displaying an asymmetric Fano line shape and an enhanced transparency window. The primary mechanisms generating such Fano-like peaks are the coupling of magnonic modes within cavities and the interference between direct transmission and resonant tunneling paths. Further increasing the nonlinear gain () of an OPA improves the transmission rate and wider transparency window, indicating amplification in our system, as indicated by the blue solid curve in Fig. 5 (a). This is because increasing the nonlinear gain () of an OPA increases the magnon number, resulting in a broader MMIT window [91] and amplification in our system. In Fig. 5(b), we plot the transmission rate () versus detuning () for different phases of an OPA (, and ) we see that phase shift of an OPA also affects the window width as well as the transmission peak. The maximum transmission peak can be obtained when we set the OPA phase . Hence, we can effectively regulate the transmission rate by modifying the OPA settings.




III.3 Controllable phase transition and MMIT in CMM system
It is commonly established that systems with -symmetry can experience a phase transition, leading to several fascinating and novel phenomena [101, 102].To investigate the occurrence of a phase transition in our -symmetric system, we focus primarily on the optical modes, taking into account both the optical loss and gain. Furthermore, we ignore the influence of OPA, and . This yields an effective non-Hermitian Hamiltonian, as in [103, 88, 104]:
(9) | ||||
Equation (9) can be straightforwardly diagonalized via a Bogoliubov transformation [105], and we obtain the complex frequency as:
(10) |
To investigate the -broken and -symmetric regimes in our system we plot the real and imaginary parts of eigenfrequencies ( against the coupling strength , as depicted in Fig. 6. First, we investigate the equal gain and loss scenario (), as depicted by the blue-solid and red-dashed curves in Figs. 6(ab). In this scenario, two distinct regimes can be observed: (i) when in this scenario the eigenfrequencies of the two supermodes exhibit two distinct frequencies (), and an identical linewidth (, indicating that our system is in a -symmetry regime (as shown by the white region in Figs. 6((a)(b))). (ii) when in this scenario, the two supermodes coalesce and have the identical frequency and different linewidths, indicating that our system is in a broken -symmetry regime (as illustrated by the light pink shaded area in Figs. 6(ab)).
The critical point, i.e., , where the resonance frequencies and the linewidths of the two supermodes become degenerate, is known as the EP or -phase transition point in our hybrid system. As demonstrated in Figs. 6(ab), by enhancing the coupling rate , one can observe the transition from a broken to an unbroken -phase. Unlike the -symmetric system with gain and loss, if both cavities are operating in a passive-passive mode with the same loss rate, i.e., , the supermodes always have different frequencies but identical linewidths. The phase transition point does not occur [see black-solid and green-dotted lines in Figs. 6(ab)].
Now, we examine the scenario where the gain/loss are not perfectly matched i.e., , where shows the gain/loss rate. In Figs. 6(cd) we consider two cases, i.e., () and (). Compared to the equal gain/loss scenario (as in Figs. 6(ab)), one can clearly see that by varying the gain/loss rate ( or ) causes the EP to shift the right or left, which implies that either a weaker or a stronger coupling strength is necessary for -phase transition.
The simplification to gain and loss in the cavities is done to isolate and study the fundamental -symmetric behavior of the system. This approach allows us to focus on the key aspects of -symmetry—gain and loss—without introducing additional complexities that could obscure the desired physical insights. Once we establish a clear understanding of -symmetry in this context, more complex interactions can be systematically added to extend the analysis.
In the above subsections §III.1 and §III.2, we examine the impact of various parameters on the transmission rate. However, we keep , which means we work with a single-loss cavity system. Now, our main task is to incorporate the effect of the tunneling rate and observe its impact on the transmission of our CMM system.


Figure 7 (ab) displays the plot of the optical transmission rate () as a function of optical detuning () in the passive-passive and active-passive CMM systems for different optical tunneling strengths , respectively. In Fig. 7(a) for the passive-passive case, we observed that when the optical tunneling strength is zero, there is a higher absorption dip at the resonance point with a small transmission peak on the left side of the resonance point, as shown by the red dotted curve in Fig. 7(a). We noticed that increasing the optical tunneling strength () for the passive-passive case in our CMM system reduced transmission. This occurs because the increased optical tunneling strength in a passive-passive scenario leads to the formation of two hybrid supermodes, , arising from the coupling between the two cavities [101, 106, 107]. The interference between these hybrid supermodes reduces the efficiency of light transmission through the system, thereby compressing the overall transmission.
Compared to the passive-passive scenario (Fig. 7(a)), the optical transmission rate in the passive-active CMM system (Fig. 7(b)) is significantly boosted. The reason for this enhancement is that as the optical tunneling strength increases, the probe light absorption by the passive cavity weakens at resonance, resulting in a transparency window in the spectrum. The increased photon tunneling strengthens the normal-mode splitting between the passive and active cavities, disrupting resonance absorption. Consequently, Stokes scattering is significantly suppressed while the anti-Stokes field (a two-phonon process) is amplified [108]. When the anti-Stokes field becomes degenerates with the probe field, destructive interference between the light pathways in the coupled cavities suppresses absorption, allowing more light to pass through and thus increasing transmission. Therefore, amplification in our CMM system can be tuned by adjusting the optical tunneling strength [108, 84].
For better comprehension, we provide the analytical results of the optical transmission rate versus optical detuning in two types of systems by setting the optical tunneling strength as in Fig. 7(c). In Fig. 7(c), the black-dotted curve illustrates the conventional passive-passive scenario, whereas the red-solid curve represents the passive-active -symmetric-like system. For the passive-passive scenario, an Autler-Townes-splitting-like spectrum is observed [109, 89]. However, unlike the passive-passive scenario, a strong signal amplification can be observed between the two Autler-Townes absorption dips in a passive-active scenario. The significant amplification peak is centered at , as illustrated by the red solid line in Fig. 7(c). Figure. 7(d) depicts the logarithm of the transmission rate vs optical tunneling strength . For the -symmetric system with gain and loss (passive-active scenario), as increases, the signal transmission gradually changes from absorption to amplification. Significant signal enhancement is achievable at the phase transition point (EP), indicated by the green solid line in Fig. 7(d). Moreover, perfect signal absorption can be achieved in the broken -symmetric regime. Additionally, when the coupling strength is low, no signal amplification occurs in a system without gain, and significant signal absorption is observed. Consequently, the transmission rate is always less than 1 in the passive-passive system, as shown by the dotted-dashed curve in Fig. 7(d). As increases beyond the phase transition point, the logarithm of the transmission rate forms an approximately straight line near zero. Therefore, in the -symmetric CMM system, the absorption and amplification of the probe light can be controlled by adjusting the optical tunneling strength .
Correspondingly, Fig. 8 shows the contour (2D) map of transmission rate vs cavity detuning and relative phase for two different cases. From Fig. 8(a), it is observed that for the passive-passive case, due to the change of relative phase , a periodic change of the transmission rate takes place, and strong absorption is also observed near the point. In Fig. 8(b), the passive-active scenario shows strong amplification at the point. The reason for this is that in the passive-active scenario, our -symmetric CMM system can add nonlinear effects that increase amplification. These nonlinearities can enhance the strength of the signal in the gain region, resulting in greater overall amplification than systems with solely passive components as in Fig.8(a).


III.4 Optimization of group delays in a compound CMM system
In general, the CMM systems can not only achieve the MMIT effect but also cause the light to slow or advance [66, 91]. This characteristic can be explained through the optical group delay. The optical group delay of the output probe light is [110, 44, 111]:
(11) |
A group delay greater than zero () indicates the phenomenon of slow light, whereas a negative group delay () corresponds to the fast light phenomenon. Specifically, the recent advancements in research utilizing slow light have led to a diverse range of applications, including telecommunications and optical data processing.
The optical group delay () of our system’s probing field is currently being examined. To illustrate this, in Fig. 9 we plot the optical group delay () versus optical detuning . In Fig. 9(a), we depict the group delay () of the probe light for differed driving magnetic field amplitudes while keeping the relative phase at . It is evident that the group delay can be adjusted by varying the magnetic field amplitude of the magnon driving field. When the magnetic field amplitude , a negative group delay is observed [as shown by the blue dashed-dotted line in Fig. 9(a)], indicating the fast light effect of the output probe field. Setting the magnetic field amplitude T reduces the negative group delay at the resonance point and increases the positive group delay on the right side of the resonance point [as shown by the black dotted line in Fig. 9(a)], indicating that slow or fast light can be achieved. As the magnetic field amplitude increases to T, the positive group delay rises and reaches a maximum value of near the resonance point (refer to the red solid line in Fig. 9(a)), indicating that the magnon-phonon interaction leads to slow light. Figure. 9(b) illustrates that the group delay can be adjusted by modulating the relative phase of the applied field. By comparing the group delay profiles for the relative phases , and [as shown by the red and black solid lines in Fig. 9(b)], we observe that the positive group delay at the probe detuning corresponding to magnon-induced absorption, significantly decreases and shifts in the case of . This investigation demonstrates that for a fixed value of the magnetic field amplitude (), the dispersion property of the system can be modulated, and the slow and the fast light can be achievable and switchable in the CMM system by using the phase modulation of the applied fields.


To further investigate the group delay of the probe light in our CMM system, we plot the group delay () versus pump power for various values of OPA gains: (blue dashed-dot curve), (black dashed curve), and (red solid curve) as shown in Fig 10. Here, we consider two different scenarios, (i) passive-active system [as shown in Fig 10(a)], (ii) passive-passive system [as shown in Fig 10(b)].
(i) In Fig 10(a), for passive-active system (). In the absence of a degenerate (), the positive group delay of the probe light decreases sharply as the pump power increases, reaching , and then remains constant with further increases in pump power. This behavior corresponds to the slow light effect [see the blue dashed-dot curve in Fig 10(a)]. In sharp contrast, introducing non-linear gain from an OPA (, and ) into the system and increasing the pump power () causes the group delay () of the probe light to become negative [see the light yellow region in Fig 10(a)], resulting in the observation of fast light. One advantage of our CMM system is its ability to gradually produce either fast or slow light by adjusting the OPA gains.
(ii) In Fig 10(b) for passive-passive system (). When compared to the passive-active system, it is evident that the positive group delay has a lower value and decreases sharply as the pump power increases [blue dashed-dot curve in Fig 10(b)], indicating that most of the region exhibits negative group delay and the majority of the probe light demonstrates fast-light behavior. Surprisingly, in the presence of an OPA (, and ), the behavior differs significantly from that of the passive-active system. The positive group delay peak increases with OPA gains, then sharply decreases as the pump power increases up to , after which it remains constant. This suggests that the OPA gains can dramatically affect the dispersion of the system.
In summary, this subsection analyzes the optical group delay of the probe field against probe detuning and pump power, considering different system parameter settings. Our observations suggest that careful adjustment of system parameters can strongly influence the group delay, thereby modulating the speed of light transmission and enabling a transition from slow to fast light within the system. This tunability offers flexibility in controlling the optical properties and can be leveraged for various applications in optical communication and signal processing [52].
III.5 Possible experimental setup
In the following subsection, we present numerical values to demonstrate that the assumptions in our analysis are reasonable and consistent with current state-of-the-art experimental achievements, thereby reinforcing the experimental feasibility of our proposal. We consider a microwave cavity containing YIG spheres and an OPA. This method can be experimentally implemented using a planar cross-shaped cavity or a coplanar waveguide made from high-conductivity copper. The choice of cavity structure can affect the frequency range of the magnomechanical interactions studied, as different designs may support various modes and resonances. Planar structures offer benefits in terms of integration with other components, such as magnetic materials or sensors, potentially enhancing the functionality of the device.
The YIG spheres, with a diameter of m, function as both phonon and magnon resonators, featuring a phonon frequency of MHz and the spin density [69, 82]. Because of its exceptional material as well as geometrical characteristics, the YIG sphere is also a great mechanical resonator. The phonon along magnon modes is connected to the changing magnetization brought about by the magnon excitation, which deforms the YIG sphere’s spherical shape as well as vice versa.
To measure group delay experimentally, techniques like interferometry and pump-probe measurements can be employed. Interferometry compares a reference pulse with the output pulse to determine phase shifts related to group delay, while pump-probe methods use a short pump pulse to excite the system and a delayed probe pulse to analyze the response. High-speed photodetectors can capture the temporal characteristics of the output light, enabling precise determination of group delay and insight into the system’s dynamics. When comparing time-delay values with state-of-the-art setups, recent experiments have demonstrated highly controllable and tunable slow and fast light effects in various cavity magnonic systems [112, 99] achieving group delays on the order of nanoseconds. Notably, the cavity systems in Refs. [113, 114, 115] demonstrate group delays in the microsecond range, which are very close to the values achieved in our system. In our system, by utilizing experimentally realizable parameters, we achieve both pulse advancement and delay, specifically ranging from = -7 and = 9. This demonstrates excellent agreement with the experimental realizations of fast and slow light reported in Refs. [116, 117] and can be further improved with longer lifetimes in magnonic systems. Recent advancements in cavity magnonics make it feasible to realize the characteristics we have reported in state-of-the-art laboratory experiments, and we believe that the system parameters and structure we propose are achievable with current experimental techniques.
IV Conclusions
In our work, we have theoretically investigated the features of MMIT, including the transmission rate, and associated group delay, in a two-optically-coupled CMM system. One passive cavity containing an OPA and two YIG spheres to establish magnon-photon coupling. The photon exchange interaction connected the active and passive microwave cavities. We used a perturbation method to establish the analytical formulas for the optical transmission rate.
Our findings suggest that the transmission spectrum is significantly affected by variations in system parameters. It is worth noting that, unlike the single transparency window observed without magnon-phonon interaction, the presence of magnon-phonon interaction leads to the appearance of two transparency windows and three absorption dips. This phenomenon arises from the combined effects of magnon-phonon and magnon-photon interactions.
Our findings further show that the two absorption dips flanking the central absorption dip may be asymmetrically tuned into amplification and absorption by altering the system parameters. Notably, through the adjustable coupling strength, by modulating both the loss and gain parameters the system can transition from the -symmetry phase to the -symmetry-broken phase, which are characterized by distinct normal mode splitting and linewidths. We also demonstrated that the active cavity generates an effective gain, which significantly boosts the transmission of the probe light in the active-passive case compared to the passive-passive case.
Furthermore, by analyzing the group delay of the probe light and adjusting system settings, a transition from slow to fast light can be achieved. We believe that our findings related to the introduced CMM are significant in the sense that these observations can be applied to real-world experiments where controlling transmission light is an important issue to handle, and thus this work may lay the groundwork for future advances in quantum control associated with CMM systems.
Acknowledgements
Abdul Wahab expresses gratitude for the financial support provided by the China Postdoctoral Research Council and Jiangsu University. This work was supported by the Natural Science Foundation of Jiangsu Province (Grant No. BK20231320), National Natural Science Foundation of China (Grant No. 12174157) and Jiangsu Funding Program for Excellent Postdoctoral Talent No.2024ZB867.
Appendix A Appendix
A.1 Calculations of first-order sideband
In this appendix, we present detailed calculations to obtain the amplitudes of the first-order sideband:
(12) | ||||
with , , , , , , , , , and .
After solving Eq. (12) for () we get the following constants for first-order side-band:
,
,
,
,
,
,
,
,
,
,
,
,
,
,
.
A.2 Stability analysis
In this appendix we determine the stability of the steady states of our system using the Routh-Hurwitz criterion [118]. The fluctuation terms of Eq. (5) in compact matrix form can be rewritten as:
(13) |
where vector :
, in which denotes the transpose of a matrix. The matrix is
given by:
Where , , , , , and . The characteristic equation can be reduced to (with ) where the coefficients can be derived using straightforward but tedious algebra. From the Routh-Hurwitz criterion [118], a solution is stable only if the real part of the corresponding eigenvalue is negative [98]. These analyses confirm that the experimentally accessible parameters in the main text can keep our compound CMM system in a stable zone.
References
- Zare Rameshti et al. [2022] B. Zare Rameshti, S. Viola Kusminskiy, J. A. Haigh, K. Usami, D. Lachance-Quirion, Y. Nakamura, C.-M. Hu, H. X. Tang, G. E. Bauer, and Y. M. Blanter, Phys. Rep. 979, 1 (2022).
- Schlawin et al. [2022] F. Schlawin, D. M. Kennes, and M. A. Sentef, Appl. Phys. Rev. 9, 011312 (2022).
- Qin et al. [2022] Y. Qin, S.-C. Li, K. Li, and J.-J. Song, Phys. Rev. B 106, 054419 (2022).
- Yuan et al. [2022] H. Yuan, Y. Cao, A. Kamra, R. A. Duine, and P. Yan, Phys. Rep. 965, 1 (2022).
- Zheng et al. [2023] S. Zheng, Z. Wang, Y. Wang, F. Sun, Q. He, P. Yan, and H. Y. Yuan, J. Appl. Phys 134, 151101 (2023).
- Wan et al. [2024] Q.-K. Wan, H.-L. Shi, and X.-W. Guan, Phys. Rev. B 109, L041301 (2024).
- Zuo et al. [2024] X. Zuo, Z.-Y. Fan, H. Qian, M.-S. Ding, H. Tan, H. Xiong, and J. Li, New J. Phys. 26, 031201 (2024).
- Lachance-Quirion et al. [2019] D. Lachance-Quirion, Y. Tabuchi, A. Gloppe, K. Usami, and Y. Nakamura, Appl. Phys. Express 12, 070101 (2019).
- Zhang [2023] X. Zhang, Mater. Today Electron. 5, 100044 (2023).
- Chumak et al. [2015] A. V. Chumak, V. I. Vasyuchka, A. A. Serga, and B. Hillebrands, Nat. Phys. 11, 453 (2015).
- Bai et al. [2017] L. Bai, M. Harder, P. Hyde, Z. Zhang, C.-M. Hu, Y. P. Chen, and J. Q. Xiao, Phys. Rev. Lett. 118, 217201 (2017).
- Zhang et al. [2015] X. Zhang, C.-L. Zou, N. Zhu, F. Marquardt, L. Jiang, and H. X. Tang, Nat. Commun. 6, 8914 (2015).
- Shen et al. [2021] R.-C. Shen, Y.-P. Wang, J. Li, S.-Y. Zhu, G. S. Agarwal, and J. Q. You, Phys. Rev. Lett. 127, 183202 (2021).
- Ihn et al. [2020] Y. S. Ihn, S.-Y. Lee, D. Kim, S. H. Yim, and Z. Kim, Phys. Rev. B 102, 064418 (2020).
- Chai et al. [2022] C.-Z. Chai, Z. Shen, Y.-L. Zhang, H.-Q. Zhao, G.-C. Guo, C.-L. Zou, and C.-H. Dong, Photon. Res. 10, 820 (2022).
- Zhang et al. [2019a] X. Zhang, K. Ding, X. Zhou, J. Xu, and D. Jin, Phys. Rev. Lett. 123, 237202 (2019a).
- Yuan et al. [2020] H. Y. Yuan, P. Yan, S. Zheng, Q. Y. He, K. Xia, and M.-H. Yung, Phys. Rev. Lett. 124, 053602 (2020).
- Yang et al. [2020] Y. Yang, Y.-P. Wang, J. W. Rao, Y. S. Gui, B. M. Yao, W. Lu, and C.-M. Hu, Phys. Rev. Lett. 125, 147202 (2020).
- Wang et al. [2019] Y.-P. Wang, J. W. Rao, Y. Yang, P.-C. Xu, Y. S. Gui, B. M. Yao, J. Q. You, and C.-M. Hu, Phys. Rev. Lett. 123, 127202 (2019).
- Zhao et al. [2020a] Y. Zhao, J. Rao, Y. Gui, Y. Wang, and C.-M. Hu, Phys. Rev. Appl. 14, 014035 (2020a).
- Hei et al. [2023] X.-L. Hei, P.-B. Li, X.-F. Pan, and F. Nori, Phys. Rev. Lett. 130, 073602 (2023).
- Zhang et al. [2019b] Z. Zhang, M. O. Scully, and G. S. Agarwal, Phys. Rev. Res. 1, 023021 (2019b).
- Nair and Agarwal [2020] J. M. P. Nair and G. S. Agarwal, Appl. Phys. Lett. 117, 084001 (2020).
- Wolski et al. [2020] S. P. Wolski, D. Lachance-Quirion, Y. Tabuchi, S. Kono, A. Noguchi, K. Usami, and Y. Nakamura, Phys. Rev. Lett. 125, 117701 (2020).
- Lachance-Quirion et al. [2020] D. Lachance-Quirion, S. P. Wolski, Y. Tabuchi, S. Kono, K. Usami, and Y. Nakamura, Science 367, 425 (2020).
- El-Ganainy et al. [2018] R. El-Ganainy, K. G. Makris, M. Khajavikhan, Z. H. Musslimani, S. Rotter, and D. N. Christodoulides, Nature Physics 14, 11 (2018).
- Ji and Yang [2024] X. Ji and X. Yang, J. Phys. Condens. Matter 36, 243001 (2024).
- Miri and Alù [2019] M.-A. Miri and A. Alù, Science 363, eaar7709 (2019).
- Xiao et al. [2021] L. Xiao, T. Deng, K. Wang, Z. Wang, W. Yi, and P. Xue, Phys. Rev. Lett. 126, 230402 (2021).
- Konotop et al. [2016] V. V. Konotop, J. Yang, and D. A. Zezyulin, Rev. Mod. Phys. 88, 035002 (2016).
- Longhi [2018] S. Longhi, Europhys. Lett. 120, 64001 (2018).
- Iorsh et al. [2020] I. Iorsh, A. Poshakinskiy, and A. Poddubny, Phys. Rev. Lett. 125, 183601 (2020).
- Makris et al. [2008] K. G. Makris, R. El-Ganainy, D. N. Christodoulides, and Z. H. Musslimani, Phys. Rev. Lett. 100, 103904 (2008).
- Guo et al. [2009] A. Guo, G. J. Salamo, D. Duchesne, R. Morandotti, M. Volatier-Ravat, V. Aimez, G. A. Siviloglou, and D. N. Christodoulides, Phys. Rev. Lett. 103, 093902 (2009).
- Longhi [2011] S. Longhi, J. Phys. A-Math. Theor. 44, 485302 (2011).
- Feng et al. [2013] L. Feng, Y.-L. Xu, W. S. Fegadolli, M.-H. Lu, J. E. B. Oliveira, V. R. Almeida, Y.-F. Chen, and A. Scherer, Nat. Mater. 12, 108 (2013).
- Peng et al. [2014] B. Peng, Ş. K. Özdemir, F. Lei, F. Monifi, M. Gianfreda, G. L. Long, S. Fan, F. Nori, C. M. Bender, and L. Yang, Nat. Phys. 10, 394 (2014).
- Goldzak et al. [2018] T. Goldzak, A. A. Mailybaev, and N. Moiseyev, Phys. Rev. Lett. 120, 013901 (2018).
- Zhang et al. [2017a] D. Zhang, X.-Q. Luo, Y.-P. Wang, T.-F. Li, and J. Q. You, Nat. Commun. 8, 1368 (2017a).
- Zhao et al. [2020b] J. Zhao, Y. Liu, L. Wu, C.-K. Duan, Y.-x. Liu, and J. Du, Phys. Rev. Appl. 13, 014053 (2020b).
- Wang and Gou [2023] F. Wang and C. Gou, Opt. Lett. 48, 1164 (2023).
- Kong et al. [2019a] C. Kong, H. Xiong, and Y. Wu, Phys. Rev. Appl. 12, 034001 (2019a).
- Bayati et al. [2024] S. Bayati, M. B. Harouni, and A. Mahdifar, Opt. Express 32, 14914 (2024).
- Ullah et al. [2020] K. Ullah, M. T. Naseem, and O. E. Müstecaplıoğlu, Phys. Rev. A 102, 033721 (2020).
- Li et al. [2020] X. Li, W.-X. Yang, T. Shui, L. Li, X. Wang, and Z. Wu, J. Appl. Phys. 128, 233101 (2020).
- Sarma et al. [2021] B. Sarma, T. Busch, and J. Twamley, New J. Phys. 23, 043041 (2021).
- Sun et al. [2021] F.-X. Sun, S.-S. Zheng, Y. Xiao, Q. Gong, Q. He, and K. Xia, Phys. Rev. Lett. 127, 087203 (2021).
- Wang et al. [2023] X. Wang, K.-W. Huang, Q.-Y. Qiu, and H. Xiong, Chaos Solitons Fractals 176, 114137 (2023).
- Shen et al. [2022] R.-C. Shen, J. Li, Z.-Y. Fan, Y.-P. Wang, and J. Q. You, Phys. Rev. Lett. 129, 123601 (2022).
- Xu et al. [2020a] W.-L. Xu, X.-F. Liu, Y. Sun, Y.-P. Gao, T.-J. Wang, and C. Wang, Phys. Rev. E 101, 012205 (2020a).
- Peng et al. [2023] M. Peng, H. Zhang, Q. Zhang, T.-X. Lu, I. M. Mirza, and H. Jing, Phys. Rev. A 107, 033507 (2023).
- Fiore et al. [2011] V. Fiore, Y. Yang, M. C. Kuzyk, R. Barbour, L. Tian, and H. Wang, Phys. Rev. Lett. 107, 133601 (2011).
- Jing et al. [2015] H. Jing, Ş. K. Özdemir, Z. Geng, J. Zhang, X.-Y. Lü, B. Peng, L. Yang, and F. Nori, Sci. Rep. 5, 9663 (2015).
- Li et al. [2016] W. Li, Y. Jiang, C. Li, and H. Song, Sci. Rep. 6, 31095 (2016).
- Vanderhaegen et al. [2021] G. Vanderhaegen, C. Naveau, P. Szriftgiser, A. Kudlinski, M. Conforti, A. Mussot, M. Onorato, S. Trillo, A. Chabchoub, and N. Akhmediev, Proc. Natl. Acad. Sci. U.S.A. 118, e2019348118 (2021).
- Huai et al. [2019] S.-N. Huai, Y.-L. Liu, J. Zhang, L. Yang, and Y.-x. Liu, Phys. Rev. A 99, 043803 (2019).
- Liu et al. [2022] Y. Liu, L. Ling, T. Shui, N. Ji, S. Liu, and W.-X. Yang, J. Opt. Soc. Am. B 39, 1042 (2022).
- Hussain et al. [2022] B. Hussain, S. Qamar, and M. Irfan, Phys. Rev. A 105, 063704 (2022).
- Hu et al. [2019] C.-S. Hu, L.-T. Shen, Z.-B. Yang, H. Wu, Y. Li, and S.-B. Zheng, Phys. Rev. A 100, 043824 (2019).
- Wahab et al. [2024a] A. Wahab, M. Abbas, X. Yang, and Y. Chen, Chaos Solitons Fractals 187, 115436 (2024a).
- Shahidani et al. [2014] S. Shahidani, M. H. Naderi, M. Soltanolkotabi, and S. Barzanjeh, J. Opt. Soc. Am. B 31, 1087 (2014).
- Huang and Agarwal [2009a] S. Huang and G. S. Agarwal, Phys. Rev. A 80, 033807 (2009a).
- Pina-Otey et al. [2016] S. Pina-Otey, F. Jiménez, P. Degenfeld-Schonburg, and C. Navarrete-Benlloch, Phys. Rev. A 93, 033835 (2016).
- Huang and Agarwal [2009b] S. Huang and G. S. Agarwal, Phys. Rev. A 79, 013821 (2009b).
- Zhang et al. [2024] W. Zhang, S. Liu, S. Zhang, and H.-F. Wang, Phys. Rev. A 109, 043712 (2024).
- Kong et al. [2019b] C. Kong, B. Wang, Z.-X. Liu, H. Xiong, and Y. Wu, Opt. Express 27, 5544 (2019b).
- Wahab et al. [2024b] A. Wahab, M. Abbas, N. Akhtar, X. Yang, and Y. Chen, Opt. Quantum Electron. 56, 1971 (2024b).
- Zimmer and Fleischhauer [2004] F. Zimmer and M. Fleischhauer, Phys. Rev. Lett. 92, 253201 (2004).
- Zhang et al. [2016] X. Zhang, C.-L. Zou, L. Jiang, and H. X. Tang, Sci. Adv. 2, e1501286 (2016).
- Osada et al. [2016] A. Osada, R. Hisatomi, A. Noguchi, Y. Tabuchi, R. Yamazaki, K. Usami, M. Sadgrove, R. Yalla, M. Nomura, and Y. Nakamura, Phys. Rev. Lett. 116, 223601 (2016).
- Hou and Liu [2019] J. T. Hou and L. Liu, Phys. Rev. Lett. 123, 107702 (2019).
- Xu et al. [2020b] J. Xu, C. Zhong, X. Han, D. Jin, L. Jiang, and X. Zhang, Phys. Rev. Lett. 125, 237201 (2020b).
- Kittel [1958] C. Kittel, Phys. Rev. 110, 836 (1958).
- Li and Zhu [2019] J. Li and S.-Y. Zhu, New J. Phys. 21, 085001 (2019).
- Li et al. [2018a] J. Li, S.-Y. Zhu, and G. S. Agarwal, Phys. Rev. Lett. 121, 203601 (2018a).
- Li et al. [2019] J. Li, S.-Y. Zhu, and G. S. Agarwal, Phys. Rev. A 99, 021801 (2019).
- Wang et al. [2018a] Y.-P. Wang, G.-Q. Zhang, D. Zhang, T.-F. Li, C.-M. Hu, and J. Q. You, Phys. Rev. Lett. 120, 057202 (2018a).
- Wang et al. [2016] Y.-P. Wang, G.-Q. Zhang, D. Zhang, X.-Q. Luo, W. Xiong, S.-P. Wang, T.-F. Li, C.-M. Hu, and J. Q. You, Phys. Rev. B 94, 224410 (2016).
- Liao et al. [2023] Q. Liao, M. Song, and W. Bao, Chaos Solitons Fractals 166, 112978 (2023).
- Potts et al. [2021] C. A. Potts, E. Varga, V. A. S. V. Bittencourt, S. V. Kusminskiy, and J. P. Davis, Phys. Rev. X 11, 031053 (2021).
- Xu et al. [2021] Y. Xu, J.-Y. Liu, W. Liu, and Y.-F. Xiao, Phys. Rev. A 103, 053501 (2021).
- Nair et al. [2022] J. M. P. Nair, D. Mukhopadhyay, and G. S. Agarwal, Phys. Rev. B 105, 214418 (2022).
- Yu et al. [2020] M. Yu, H. Shen, and J. Li, Phys. Rev. Lett. 124, 213604 (2020).
- Su et al. [2024] Z.-X. Su, A.-D. Zhu, and L. Yu, Phys. Rev. A 110, 023502 (2024).
- Suzuki et al. [2015] H. Suzuki, E. Brown, and R. Sterling, Phys. Rev. A 92, 033823 (2015).
- Weis et al. [2010] S. Weis, R. Rivière, S. Deléglise, E. Gavartin, O. Arcizet, A. Schliesser, and T. J. Kippenberg, Science 330, 1520 (2010).
- Xiong et al. [2012] H. Xiong, L.-G. Si, A.-S. Zheng, X. Yang, and Y. Wu, Phys. Rev. A 86, 013815 (2012).
- Das et al. [2023] S. Das, S. Chakraborty, and T. N. Dey, Phys. Rev. A 108, 033517 (2023).
- Lai et al. [2020] D.-G. Lai, X. Wang, W. Qin, B.-P. Hou, F. Nori, and J.-Q. Liao, Phys. Rev. A 102, 023707 (2020).
- Huang and Agarwal [2011] S. Huang and G. S. Agarwal, Phys. Rev. A 83, 043826 (2011).
- Lu et al. [2023] T.-X. Lu, X. Xiao, L.-S. Chen, Q. Zhang, and H. Jing, Phys. Rev. A 107, 063714 (2023).
- Liu et al. [2024] W. Liu, M. Abbas, S. H. Asadpour, H. R. Hamedi, P. Zhang, and B. C. Sanders, New J. Phys. 26, 093042 (2024).
- Li et al. [2018b] L. Li, W.-X. Yang, Y. Zhang, T. Shui, A.-X. Chen, and Z. Jiang, Phys. Rev. A 98, 063840 (2018b).
- Harder et al. [2018] M. Harder, Y. Yang, B. M. Yao, C. H. Yu, J. W. Rao, Y. S. Gui, R. L. Stamps, and C.-M. Hu, Phys. Rev. Lett. 121, 137203 (2018).
- Yu et al. [2019] C. H. Yu, Y. Yang, J. W. Rao, P. Hyde, Y.-P. Wang, B. Zhang, Y. S. Gui, and C.-M. Hu, AIP Advances 9, 115012 (2019).
- Xu et al. [2019] P.-C. Xu, J. W. Rao, Y. S. Gui, X. Jin, and C.-M. Hu, Phys. Rev. B 100, 094415 (2019).
- Jia et al. [2015] W. Z. Jia, L. F. Wei, Y. Li, and Y.-x. Liu, Phys. Rev. A 91, 043843 (2015).
- Lu et al. [2019] T.-X. Lu, Y.-F. Jiao, H.-L. Zhang, F. Saif, and H. Jing, Phys. Rev. A 100, 013813 (2019).
- Zhao et al. [2021] J. Zhao, L. Wu, T. Li, Y.-x. Liu, F. Nori, Y. Liu, and J. Du, Phys. Rev. Appl. 15, 024056 (2021).
- Li et al. [2017] Y. Li, Y. Y. Huang, X. Z. Zhang, and L. Tian, Opt. Express 25, 18907 (2017).
- Lü et al. [2015] X.-Y. Lü, H. Jing, J.-Y. Ma, and Y. Wu, Phys. Rev. Lett. 114, 253601 (2015).
- Bender et al. [2013] C. M. Bender, M. Gianfreda, i. m. c. K. Özdemir, B. Peng, and L. Yang, Phys. Rev. A 88, 062111 (2013).
- Wang et al. [2022] B. Wang, X. Jia, X.-H. Lu, and H. Xiong, Phys. Rev. A 105, 053705 (2022).
- Wang et al. [2024] Z.-Y. Wang, X.-W. He, X. Han, H.-F. Wang, and S. Zhang, Opt. Express 32, 4987 (2024).
- Walls and Milburn [2012] D. Walls and G. Milburn, Quantum Optics (Springer Berlin Heidelberg, 2012).
- Zhang [2020] Y. Zhang, Phys. Rev. A 101, 023842 (2020).
- Zhang et al. [2017b] X. Y. Zhang, Y. Q. Guo, P. Pei, and X. X. Yi, Phys. Rev. A 95, 063825 (2017b).
- Kundu et al. [2021] A. Kundu, C. Jin, and J.-X. Peng, Ann. Phys. 429, 168465 (2021).
- Wang et al. [2018b] X. Wang, A. Miranowicz, H.-R. Li, F.-L. Li, and F. Nori, Phys. Rev. A 98, 023821 (2018b).
- Zhao et al. [2022] C. Zhao, Z. Yang, R. Peng, J. Yang, C. Li, and L. Zhou, Phys. Rev. Appl. 18, 044074 (2022).
- Lü et al. [2017] H. Lü, Y. Jiang, Y.-Z. Wang, and H. Jing, Photon. Res. 5, 367 (2017).
- Zhang et al. [2014] X. Zhang, C.-L. Zou, L. Jiang, and H. X. Tang, Phys. Rev. Lett. 113, 156401 (2014).
- Qin et al. [2020] T. Qin, J. Yang, F. Zhang, Y. Chen, D. Shen, W. Liu, L. Chen, X. Jiang, X. Chen, and W. Wan, Commun. Phys. 3, 118 (2020).
- Zhu et al. [2024] Y. Zhu, J. Hou, Q. Geng, B. Xue, Y. Chen, X. Chen, L. Ge, and W. Wan, Nat. Commun. 15, 8101 (2024).
- Lake et al. [2020] D. P. Lake, M. Mitchell, B. C. Sanders, and P. E. Barclay, Nat. Commun. 11, 2208 (2020).
- Clark et al. [2014] J. B. Clark, R. T. Glasser, Q. Glorieux, U. Vogl, T. Li, K. M. Jones, and P. D. Lett, Nat. Photonics 8, 515 (2014).
- Stenner et al. [2003] M. D. Stenner, D. J. Gauthier, and M. A. Neifeld, Nature 425, 695 (2003).
- DeJesus and Kaufman [1987] E. X. DeJesus and C. Kaufman, Phys. Rev. A 35, 5288 (1987).