This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Tunable optical amplification and group delay in cavity magnomechanics

Abdul Wahab [email protected] Department of Physics, Jiangsu University, Zhenjiang, 212013, China    Muqaddar Abbas Ministry of Education Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, Shaanxi Province Key Laboratory of Quantum Information and Quantum Optoelectronic Devices, School of Physics, Xi’an Jiaotong University, Xi’an 710049, China    Xiaosen Yang [email protected]    Yuanping Chen [email protected] Department of Physics, Jiangsu University, Zhenjiang, 212013, China
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 𝒫𝒯\mathcal{PT}-symmetric to broken-𝒫𝒯\mathcal{PT}-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 (𝒫𝒯\mathcal{PT})-symmetric systems. 𝒫𝒯\mathcal{PT}-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 𝒫𝒯\mathcal{PT}-symmetric phase to the broken-𝒫𝒯\mathcal{PT}-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 𝒫𝒯\mathcal{PT}-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-𝒫𝒯\mathcal{PT}-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 𝒫𝒯\mathcal{PT}-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 𝒫𝒯\mathcal{PT}-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  𝒫𝒯\mathcal{PT}-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].

Refer to caption
Figure 1: Schematic configuration of hybrid nonlinear CMM system. The system is made up of two connected microwave cavities. A passive (loss) cavity (a1a_{1}) that contains an OPA and two highly polished single-crystal ferromagnetic YIG spheres (m1,m2m_{1},m_{2}), and an active (gain) cavity (a2a_{2}). The κ1\kappa_{1} and κ2\kappa_{2} represent the loss and gain rates, respectively. Both the cavities can be interconnected through the optical tunneling rate JJ. The loss cavity is driven by a strong external pump and a weak probe laser field. The magnon m2m_{2} is directly controlled by a microwave control field with frequency Ωm\Omega_{m}. Detailed information about the schematic configuration can be shown in  §II.

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 𝒫𝒯\mathcal{PT}-symmetry phase to an unbroken 𝒫𝒯\mathcal{PT}-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 (a1a_{1} and a2a_{2}). Cavity a1a_{1} is treated as a passive (loss) cavity, with the loss rate κ1\kappa_{1} containing a degenerate OPA and two highly polished single-crystal ferromagnetic YIG spheres (with a 250 μ\mum 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 (a2a_{2}) with the gain rate κ2\kappa_{2} is coupled to the passive (loss) cavity via optical tunneling rate JJ. The CMM resonator is driven by two external optical laser fields: a weak probe and a strong pump field. The symbols εp\varepsilon_{p} and εl\varepsilon_{l} denote the amplitudes, ωp\omega_{p} and ωl\omega_{l} represent the frequencies and ϕp\phi_{p}, ϕl\phi_{l} correspond to the phases of the probe and pump fields, respectively. A weak external microwave field with frequency Ωm\Omega_{m} and the phase ϕm\phi_{m} 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 Km^jm^jm^jm^jK\hat{m}_{j}^{\dagger}\hat{m}_{j}^{\dagger}\hat{m}_{j}\hat{m}_{j} which may occur due to a strongly driven magnon mode [77, 78]. The total Hamiltonian of the system reads:

H^t=\displaystyle\hat{H}_{\text{t}}= ω1a^1a^1+ω2a^2a^2+ωbb^b^+J(a^1a^2+a^2a^1)\displaystyle\hbar\omega_{1}\hat{a}_{1}^{\dagger}\hat{a}_{1}+\hbar\omega_{2}\hat{a}_{2}^{\dagger}\hat{a}_{2}+\hbar\omega_{b}\hat{b}^{\dagger}\hat{b}+\hbar J\left(\hat{a}_{1}^{\dagger}\hat{a}_{2}+\hat{a}_{2}^{\dagger}\hat{a}_{1}\right) (1)
+j=12[ωmjm^jm^j+gj(a^1+a^1)(m^j+m^j)]\displaystyle+\hbar\sum_{j=1}^{2}\left[\omega_{mj}\hat{m}_{j}^{\dagger}\hat{m}_{j}+g_{j}\left(\hat{a}_{1}+\hat{a}_{1}^{\dagger}\right)\left(\hat{m}_{j}+\hat{m}_{j}^{\dagger}\right)\right]
+gmbm^2m^2(b^+b^)+iG(eiθa^12e2iωlt H.c.)\displaystyle+\hbar g_{mb}\hat{m}_{2}^{\dagger}\hat{m}_{2}\left(\hat{b}^{\dagger}+\hat{b}\right)+i\hbar G\left(e^{i\theta}\hat{a}_{1}^{\dagger 2}e^{-2i\omega_{l}t}-\text{ H.c.}\right)
+i2ηaκ1εl(a^1eiωltiϕl H.c.)\displaystyle+i\hbar\sqrt{2\eta_{a}\kappa_{1}}\varepsilon_{l}\left(\hat{a}_{1}^{\dagger}e^{-i\omega_{l}t-i\phi_{l}}-\text{ H.c.}\right)
+i(εmm2eiΩmtiϕmH.c.)\displaystyle+i\hbar\left(\varepsilon_{m}{m}_{2}^{\dagger}e^{-i\Omega_{m}t-i\phi_{m}}-\text{H.c.}\right)
+i2ηaκ1εp(a^1eiωptiϕp H.c.).\displaystyle+i\hbar\sqrt{2\eta_{a}\kappa_{1}}\varepsilon_{p}\left(\hat{a}_{1}^{\dagger}e^{-i\omega_{p}t-i\phi_{p}}-\text{ H.c.}\right).

We assume that both cavities in our CMM system have the same resonance frequency (ω1=ω2=ωa\omega_{1}=\omega_{2}=\omega_{a}[79]. Applying the rotating-wave approximation to the system, the interaction terms simplify as gj(a^1+a^1)(m^j+m^j)gj(m^ja^1+m^ja^1)g_{j}(\hat{a}_{1}+\hat{a}_{1}^{\dagger})(\hat{m}_{j}+\hat{m}_{j}^{\dagger})\rightarrow g_{j}(\hat{m}_{j}^{\dagger}\hat{a}_{1}+\hat{m}_{j}\hat{a}_{1}^{\dagger}), by neglecting the fast oscillating terms gj(a^1m^j+a^1m^j)g_{j}(\hat{a}_{1}\hat{m}j+\hat{a}1^{\dagger}\hat{m}j^{\dagger}). This is valid when ωagj,κ1,κmj\omega_{a}\gg g_{j},\kappa_{1},\kappa_{mj} which is the case to be considered in the present work (κ1\kappa_{1} and κmj\kappa_{mj} 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 ωl\omega_{l}. Following the transformation H^rot=RH^tR+i(Rt)R\hat{H}_{\text{rot}}=R\hat{H}_{\text{t}}R^{\dagger}+i\hbar(\frac{\partial R}{\partial t})R^{\dagger} with R=eiωl(a^1a^1+a^2a^2+m^jm^j)tR=e^{i\omega_{l}(\hat{a}_{1}^{\dagger}\hat{a}_{1}+\hat{a}_{2}^{\dagger}\hat{a}_{2}+\hat{m}_{j}^{\dagger}\hat{m}_{j})t}, the Hamiltonian in Eq. (1) can be rewritten as [80, 81]:

H^rot=\displaystyle\hat{H}_{\text{rot}}= Δa(a^1a^1+a^2a^2)+ωbb^b^+J(a^1a^2+a^2a^1)\displaystyle\hbar\Delta_{a}\left(\hat{a}_{1}^{\dagger}\hat{a}_{1}+\hat{a}_{2}^{\dagger}\hat{a}_{2}\right)+\hbar\omega_{b}\hat{b}^{\dagger}\hat{b}+\hbar J\left(\hat{a}_{1}^{\dagger}\hat{a}_{2}+\hat{a}_{2}^{\dagger}\hat{a}_{1}\right) (2)
+j=12[Δmjm^jm^j+gj(m^ja^1+mja^1)]\displaystyle+\hbar\sum_{j=1}^{2}\left[\Delta_{mj}\hat{m}_{j}^{\dagger}\hat{m}_{j}+g_{j}\left(\hat{m}_{j}^{\dagger}\hat{a}_{1}+m_{j}\hat{a}_{1}^{\dagger}\right)\right]
+gmbm^2m^2(b^+b^)+iG(eiθa^12eiθa^12)\displaystyle+\hbar g_{mb}\hat{m}_{2}^{\dagger}\hat{m}_{2}\left(\hat{b}^{\dagger}+\hat{b}\right)+i\hbar G\left(e^{i\theta}\hat{a}_{1}^{\dagger 2}-e^{-i\theta}\hat{a}_{1}^{2}\right)
+i2ηaκ1εl(a^1a^1)+i(εmeiδmltiϕmm^H.c.)\displaystyle+i\hbar\sqrt{2\eta_{a}\kappa_{1}}\varepsilon_{l}\left(\hat{a}_{1}^{\dagger}-\hat{a}_{1}\right)+i\hbar\left(\varepsilon_{m}e^{-i\delta_{ml}t-i\phi_{m}}\hat{m}^{\dagger}-\text{H.c.}\right)
+i2ηaκ1εp(a^1eiδpltiϕpl H.c.),\displaystyle+i\hbar\sqrt{2\eta_{a}\kappa_{1}}\varepsilon_{p}\left(\hat{a}_{1}^{\dagger}e^{-i\delta_{pl}t-i\phi_{pl}}-\text{ H.c.}\right),

where Δa:=ωaωl\Delta_{a}:=\omega_{a}-\omega_{l} (j=12Δmj:=ωmjωl\sum_{j=1}^{2}\Delta_{mj}:=\omega_{mj}-\omega_{l}) shows the frequency detuning between cavity (magnon) and control field. δpl:=ωpωl\delta_{pl}:=\omega_{p}-\omega_{l} stands for the frequency detuning of probe-control field, and δml:=Ωmωl\delta_{ml}:=\Omega_{m}-\omega_{l} shows the detuning of the magnon drive and control fields. ϕpl=ϕpϕl\phi_{pl}=\phi_{p}-\phi_{l} indicates the relative phase between the probe and pump fields. ωa\omega_{a}, ωb\omega_{b}, and ωm\omega_{m} are the frequencies that correspond to the resonances of the photon, phonon, and magnon modes, respectively. The term J(a^1a^2+a^2a^1)\hbar J\left(\hat{a}_{1}^{\dagger}\hat{a}_{2}+\hat{a}_{2}^{\dagger}\hat{a}_{1}\right) indicates the photon exchange coupling strength interaction between the two cavities. j=12gj\sum_{j=1}^{2}g_{j} (gmbg_{mb}) signifies the magnon-photon (phonon) coupling rates [82]. GG represents nonlinear gain, while θ\theta indicates the phase of OPA, respectively. The term εm=5N/4B0γ\varepsilon_{m}=\sqrt{5N/4}B_{0}\gamma represents the coupling strength between the magnon and the microwave field; in which N=ρVN=\rho V is the total number of spins (with VV shows the volume, and  ρ=4.22×1027cm3\rho=4.22\times 10^{27}\textup{cm}^{-3} is the spin density of YIG sphere) [75, 83], γ/2π\gamma/2\pi= 28 GHz/T represents the gyromagnetic ratio, with B0B_{0} shows the field amplitude. The terms i2ηaκ1εl(a^1a^1)+i2ηaκ1(a^1eiδpltiϕplεp H.c.)i\hbar\sqrt{2\eta_{a}\kappa_{1}}\varepsilon_{l}\left(\hat{a}_{1}^{\dagger}-\hat{a}_{1}\right)+i\hbar\sqrt{2\eta_{a}\kappa_{1}}\left(\hat{a}_{1}^{\dagger}e^{-i\delta_{pl}t-i\phi_{pl}}\varepsilon_{p}-\text{ H.c.}\right) signify the interaction between the cavity field and two input fields. The terms εl=Pl/ωl(εp=Pp/ωp)\varepsilon_{l}=\sqrt{P_{l}/\hbar\omega_{l}}(\varepsilon_{p}=\sqrt{P_{p}/\hbar\omega_{p}}) represent the amplitude of the strong pump (weak probe) fields, whereas PlP_{l} and PpP_{p} are the powers of the corresponding pump and probe fields, respectively. The symbol ηa\eta_{a} represents the coupling between the CMM resonator and the output port, defined as ηa=κc1/κ1\eta_{a}=\kappa_{c_{1}}/\kappa_{1}, where κc1\kappa_{c_{1}} is the external decay rate of the cavity and κ1\kappa_{1} is the total decay rate of the cavity fields, consisting of both intrinsic and external decay rates (κ1=κi+κc1\kappa_{1}=\kappa_{i}+\kappa_{c_{1}}[84, 85]. In this study, we consider the CMM resonator operating in the critical-coupling regime and set ηa=0.5\eta_{a}=0.5 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:

Refer to caption
Refer to caption
Refer to caption
Figure 2:  Absorption Re[α][\alpha] (blue solid line) and dispersion Im[α]Im[\alpha] (red dotted line) profiles of the output probe field are shown versus optical detuning δ/ωb\delta/\omega_{b}. We select the coupling strengths as (a) g1=2κ1g_{1}=2\kappa_{1}g2=gmb=0g_{2}=g_{mb}=0 (b) g1=g2=2κ1g_{1}=g_{2}=2\kappa_{1}gmb=0g_{mb}=0 (c) g1=g2=2κ1g_{1}=g_{2}=2\kappa_{1}gmb=2πg_{mb}=2\pi Hz. The parameter values are considered as: ωc=2π×7.86\omega_{c}=2\pi\times 7.86 GHz, κ1=π×3.35\kappa_{1}=\pi\times 3.35 MHz, κ2=κ1\kappa_{2}=\kappa_{1}D=250μD=250\mum, ρ=4.22×1027cm3\rho=4.22\times 10^{27}\textup{cm}^{-3}=1.054×1034Js\hbar=1.054\times 10^{-34}\ \mathrm{J\cdot s}ωb=2π×11.42\omega_{b}=2\pi\times 11.42 MHz, κb=300π\kappa_{b}=300\pi Hz, κm1=κm2=π×1.12\kappa_{m_{1}}=\kappa_{m_{2}}=\pi\times 1.12 MHz, B0=0B_{0}=0N=3×1016N=3\times 10^{16}γ\gamma= 2π×2\pi\times 28 GHz/T, Pl=10P_{l}=10 mW, εp=0.0001εl\varepsilon_{p}=0.0001\varepsilon_{l}ηa=0.5\eta_{a}=0.5Δa=Δm1=Δm2=ωb\Delta_{a}=\Delta_{m_{1}}=\Delta_{m_{2}}^{\prime}=\omega_{b}G=0G=0θ=0\theta=0Φm=0\Phi_{m}=0, and J=0J=0.
a˙1=\displaystyle\left\langle\dot{a}_{1}\right\rangle= (iΔaκ1)a1ij=12gjmjiJa2\displaystyle\left(-i\Delta_{a}-\kappa_{1}\right)\left\langle a_{1}\right\rangle-i\sum_{j=1}^{2}g_{j}\langle m_{j}\rangle-iJ\left\langle a_{2}\right\rangle (3)
+2ηaκ1εl+2ηaκ1εpeiδpltiϕpl+2Geiθa1,\displaystyle+\sqrt{2\eta_{a}\kappa_{1}}\varepsilon_{l}+\sqrt{2\eta_{a}\kappa_{1}}\varepsilon_{p}e^{-i\delta_{pl}t-i\phi_{pl}}+2Ge^{i\theta}a_{1}^{\dagger},
m1˙=\displaystyle\langle\dot{m_{1}}\rangle= (iΔm1κm1)m1ig1a1,\displaystyle\left(-i\Delta_{m1}-\kappa_{m1}\right)\langle m_{1}\rangle-ig_{1}\left\langle a_{1}\right\rangle,
m2˙=\displaystyle\langle\dot{m_{2}}\rangle= (iΔm2κm2)m2ig2a1\displaystyle\left(-i\Delta_{m2}-\kappa_{m2}\right)\langle m_{2}\rangle-ig_{2}\left\langle a_{1}\right\rangle
+εmeiδmltiϕmigmbm2(b+b),\displaystyle+\varepsilon_{m}e^{-i\delta_{ml}t-i\phi_{m}}-ig_{mb}\langle m_{2}\rangle\left(\left\langle b^{\dagger}\right\rangle+\langle b\rangle\right),
b˙=\displaystyle\langle\dot{b}\rangle= (iωbκb)bigmbm2m2,\displaystyle\left(-i\omega_{b}-\kappa_{b}\right)\langle b\rangle-ig_{mb}\left\langle m_{2}^{\dagger}\right\rangle\langle m_{2}\rangle,
a˙2=\displaystyle\left\langle\dot{a}_{2}\right\rangle= (iΔa+κ2)a2iJa1,\displaystyle\left(-i\Delta_{a}+\kappa_{2}\right)\left\langle a_{2}\right\rangle-iJ\left\langle a_{1}\right\rangle,

where j=12κmj\sum_{j=1}^{2}\kappa_{mj} and κb\kappa_{b} denotes the decay rates of magnons and phonons. Keep in mind that κ1>0\kappa_{1}>0, κ2>0\kappa_{2}>0 in Eq. (3) is associated to passive-active compound CMM resonator system, while κ1>0\kappa_{1}>0, κ2<0\kappa_{2}<0 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 |εp|εl|\varepsilon_{p}|\ll\varepsilon_{l} and |ε|mεl|\varepsilon|_{m}\ll\varepsilon_{l} [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., O=Os+O+eiδt+OeiδtO=O_{s}+O_{+}e^{-i\delta t}+O_{-}e^{i\delta t}, where O=a1,a2,b,m1,m2O=a_{1},a_{2},b,m_{1},m_{2} [44, 90], the first term OsO_{s} represents the steady-state values and δO=O+eiδt+Oeiδt\delta{O}=O_{+}e^{-i\delta t}+O_{-}e^{i\delta t} 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:

a1s\displaystyle a_{1s} =(λ2+2Geiθ)2ηaκ1εlλ2λ3,\displaystyle=\frac{(\lambda_{2}+2Ge^{i\theta}){\sqrt{2\eta_{a}\kappa_{1}}\varepsilon_{l}}}{\lambda_{2}\lambda_{3}}, (4a)
m1s=ig1a1siΔm1+κm1,m_{1s}=\frac{-ig_{1}a_{1s}}{i\Delta_{m1}+\kappa_{m1}}, (4b)
m2s=ig2a1siΔm2+κm2,m_{2s}=\frac{-ig_{2}a_{1s}}{i\Delta_{m2}^{\prime}+\kappa_{m2}}, (4c)
bs=igmb|m2s|2iωb+κb,\displaystyle b_{s}=\frac{-ig_{mb}\left|m_{2s}\right|^{2}}{i\omega_{b}+\kappa_{b}}, (4d)
a2s=iJa1s(iΔa+κ2).\displaystyle a_{2s}=\frac{iJa_{1s}}{\left(-i\Delta_{a}+\kappa_{2}\right)}. (4e)

Here we define: λ1=(κ1+iΔa)((κ2iΔa))J2(κ2iΔa)\lambda_{1}=\frac{{(\kappa_{1}+i\Delta_{a})((\kappa_{2}-i\Delta_{a}))-J^{2}}}{{(\kappa_{2}-i\Delta_{a})}}λ2=(κ1iΔa)+(g12iΔm1+κm1)+(g22iΔm2+κm2)(J2κ2+iΔa)\lambda_{2}=(\kappa_{1}-i\Delta_{a})+(\frac{g_{1}^{2}}{-i\Delta_{m1}+\kappa_{m1}})+(\frac{g_{2}^{2}}{-i\Delta_{m2}+\kappa_{m2}})-(\frac{J^{2}}{\kappa_{2}+i\Delta a})λ3=(λ1λ24G2λ2)+(g12iΔm1+κm1)+(g22iΔm2+κm2)\lambda_{3}=(\frac{\lambda_{1}\lambda_{2}-4G^{2}}{\lambda_{2}})+(\frac{g_{1}^{2}}{i\Delta_{m1}+\kappa_{m1}})+(\frac{g_{2}^{2}}{i\Delta_{m2}+\kappa_{m2}}). By performing a series of mathematical calculations and neglecting all higher-order terms (δoδo)(\delta o\delta o) the evolution of the perturbation terms of Eq. (3) can be shown as:

δa˙1=\displaystyle\delta\dot{a}_{1}= (iΔa+κ1)δa1iJδa2ij=12gjδmj\displaystyle-\left(i\Delta_{a}+\kappa_{1}\right)\delta a_{1}-iJ\delta a_{2}-i\sum_{j=1}^{2}g_{j}\delta m_{j} (5)
+2ηaκ1εpeiδpltiϕpl+2Geiθδa1,\displaystyle+\sqrt{2\eta_{a}\kappa_{1}}\varepsilon_{p}e^{-i\delta_{pl}t-i\phi_{pl}}+2Ge^{i\theta}\delta a_{1}^{\dagger},
δm1˙=\displaystyle\delta\dot{m_{1}}= (iΔm1+κm1)δm1ig1δa1,\displaystyle-\left(i\Delta_{m_{1}}+\kappa_{m_{1}}\right)\delta{m_{1}}-ig_{1}\delta a_{1},
δm2˙=\displaystyle\delta\dot{m_{2}}= (iΔm2+κm2)δm2ig2δa1iFδb\displaystyle-\left(i\Delta_{m_{2}}^{\prime}+\kappa_{m_{2}}\right)\delta{m_{2}}-ig_{2}\delta a_{1}-iF\delta b
iFδb+εmeiδmltiϕm,\displaystyle-iF\delta b^{\dagger}+\varepsilon_{m}e^{-i\delta_{ml}t-i\phi_{m}},
δb˙=\displaystyle\delta\dot{b}= (iωb+κb)δbiFδm2iFδm2,\displaystyle-\left(i\omega_{b}+\kappa_{b}\right)\delta b-iF\delta m_{2}^{\dagger}-iF^{*}\delta m_{2},
δa˙2=\displaystyle\delta\dot{a}_{2}= (iΔaκ2)δa2iJδa1,\displaystyle-\left(i\Delta_{a}-\kappa_{2}\right)\delta a_{2}-iJ\delta a_{1},

here, F=gmbm2sF=g_{mb}m_{2s}, and Δm2=Δm2+gmb(bs+bs)\Delta_{m_{2}}^{\prime}=\Delta_{m_{2}}+g_{mb}(b_{s}+b_{s}^{*}) 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]:

δa1=\displaystyle\delta{a}_{1}= A1+eiδt+A1eiδt,\displaystyle A_{1+}e^{-i\delta t}+A_{1-}e^{i\delta t}, (6)
δm1=\displaystyle\delta{m}_{1}= M1+eiδt+M1eiδt,\displaystyle M_{1+}e^{-i\delta t}+M_{1-}e^{i\delta t},
δm2=\displaystyle\delta{m}_{2}= M2+eiδt+M2eiδt,\displaystyle M_{2+}e^{-i\delta t}+M_{2-}e^{i\delta t},
δb=\displaystyle\delta b= B1+eiδt+B1eiδt,\displaystyle B_{1+}e^{-i\delta t}+B_{1-}e^{i\delta t},
δa2=\displaystyle\delta a_{2}= X1+eiδt+X1eiδt.\displaystyle X_{1+}e^{-i\delta t}+X_{1-}e^{i\delta t}.

We consider the magnon driving field to become resonant with the probe field frequency which leads us to consider δ=δpl=δml\delta=\delta_{pl}=\delta_{ml} [45, 92], and Ai+A_{i+} and AiA_{i-} denotes the iith 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 A1+A_{1+}, which indicates the output probe field amplitude of the CMM resonator can be written as:

A1+=(β4εmeiΦ2ηaκ1εpα13).\displaystyle A_{1+}=\left(\frac{\beta_{4}\varepsilon_{m}e^{i\Phi}-\sqrt{2\eta_{a}\kappa_{1}}\varepsilon_{p}}{\alpha_{13}}\right). (7)

Here, Φ=ϕmϕpl\Phi=\phi_{m}-\phi_{pl} corresponds to the relative phase of the applied fields. Furthermore, by employing the input/output relationship of the cavity sout=sin2ηκ1a1s_{\text{out}}=s_{\text{in}}-\sqrt{2\eta\kappa_{1}}a_{1} [93], we get the normalized output probe field intensity from the CMM resonator as:

T=|tp|2=|2ηaκ1A1+εp1|2.\displaystyle T=|t_{p}|^{2}=\left|\frac{\sqrt{2\eta_{a}\kappa_{1}}A_{1+}}{\varepsilon_{p}}-1\right|^{2}. (8)

The real and imaginary components of the α\alpha (α=tp+1\alpha=t_{p}+1) can be quantified as the absorption and dispersion of the probe field at the probe frequency, respectively.

Refer to caption
Refer to caption
Figure 3: (a) Absorption Re[α][\alpha] and (b) dispersion Im[α]Im[\alpha] profiles of the output probe field are shown versus optical detuning δ/ωb\delta/\omega_{b} for the different selected parameters of OPA gains as G=0G=0G=0.5κ1G=0.5\kappa_{1}G=1κ1G=1\kappa_{1} and θ=0\theta=0. All other parameters remains same as in Fig. 2(c) except for B0=5×1010B_{0}=5\times 10^{-10}T.

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 ωc=2π×7.86\omega_{c}=2\pi\times 7.86 GHz, the first cavity decay rate κ1=π×3.35\kappa_{1}=\pi\times 3.35 MHz, and the coupling parameter ηa=0.5\eta_{a}=0.5. The values for the YIG sphere are taken as: the diameter D=250μD=250\mum, the spin density ρ=4.22×1027cm3\rho=4.22\times 10^{27}\textup{cm}^{-3}. The frequency of the phonon mode ωb=2π×11.42\omega_{b}=2\pi\times 11.42 MHz, and decay rate κb=300π\kappa_{b}=300\pi Hz. The magnon decay rate κm1=κm2=π×1.12\kappa_{m_{1}}=\kappa_{m_{2}}=\pi\times 1.12 MHz, gyromagnetic ratio γ\gamma= 2π×2\pi\times 28 GHz/T, and N=3×1016N=3\times 10^{16}. In addition, we take the power of the pump field Pl=10P_{l}=10 mW, and the amplitude of the probe filed εp=0.0001εl\varepsilon_{p}=0.0001\varepsilon_{l}, with Δa=Δm1=Δm2=ωb\Delta_{a}=\Delta_{m_{1}}=\Delta_{m_{2}}^{\prime}=\omega_{b}.

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 (δ/ωb\delta/\omega_{b}) for various coupling strengths. In Fig. 2(a), we consider the scenario where both the magnon-phonon coupling (gmbg_{mb}) and the coupling strength of the magnon mode m2m_{2} (g2g_{2}) with the cavity are absent. As a result, only magnon mode m1m_{1} 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 m1,m2m_{1},m_{2} (g1,g2g_{1},g_{2}) are present, and only the magnon-phonon coupling (gmbg_{mb}) is absent, we observe that the typical Lorentzian-type absorption peak at resonance δ=ωb\delta=\omega_{b} transforms into a transparency window. Additionally, two symmetrical absorption dips occur at δ=0.7ωb\delta=0.7\omega_{b}(left) and δ=1.3ωb\delta=1.3\omega_{b}(right) as shown in Fig. 2(b). The central peak at δ=ωb\delta=\omega_{b} 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].

Refer to caption
Refer to caption
Refer to caption
Figure 4:  Transmission |tp|2|t_{p}|^{2} spectrum of the output probe field versus optical detuning δ/ωb\delta/\omega_{b} is shown for different values of the magnetic field amplitude (B0B_{0}) and the relative phase (Φ\Phi) of magnon driving field εm\varepsilon_{m}. (a) The magnetic field amplitude B0=0B_{0}=0B0=5×1010B_{0}=5\times 10^{-10}T, and  B0=10×1010B_{0}=10\times 10^{-10}T. (b) The relative phase Φ=0\Phi=0, Φ=π\Phi=\pi and B0=12×1010B_{0}=12\times 10^{-10}T (c) Contour map of transmission rate |tp|2|t_{p}|^{2} vs optical detuning δ/ωb\delta/\omega_{b} and relative phase Φ\Phi. All other parameters remains same as in Fig. 2(c) except for B0=5×1010B_{0}=5\times 10^{-10}T.

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 (g2g_{2}). However, the results are not shown here. Finally, we consider the scenario where all three couplings are present, i.e., g1=g2=2κ1g_{1}=g_{2}=2\kappa_{1}, and gmb=2πg_{mb}=2\pi Hz. In this case, due to the nonzero magnetostrictive interaction, the single MMIT window at the resonance point (δ=ωb\delta=\omega_{b}) 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 δ/ωb\delta/\omega_{b} for different selected OPA gain parameters GG. It is evident that changes in the parametric gain influence the transparency window. In the absence of OPA (G=0G=0), we observe clear absorption in our system. However, introducing OPA (G0G\neq 0) 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].

Refer to caption
Refer to caption
Figure 5:  Transmission rate |tp|2|t_{p}|^{2} spectrum of the output probe field versus optical detuning δ/ωb\delta/\omega_{b} for different values of OPA gains and phases. (a) The OPA gains parameters are chosen as: G=0G=0G=0.5κ1G=0.5\kappa_{1}G=1κ1G=1\kappa_{1}, and θ=0\theta=0. (b) The phases of OPA are set as: θ=0\theta=0,  θ=0.5π\theta=0.5\piθ=1.5π\theta=1.5\pi

, and G=0.5κ1G=0.5\kappa_{1}. All other parameters remains same as in Fig. 2(c) except for Φ=0\Phi=0 and B0=5×1010B_{0}=5\times 10^{-10}T.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 6:  (a) Real and (b) Imaginary components of the eigenfrequencies (λΔa)/κ1\lambda-\Delta_{a})/\kappa_{1} are depicted versus coupling strength J/k1J/k_{1} for identical gain and loss (k1=k2k_{1}=k_{2}). (c) Real and (d) imaginary components of eigenfrequencies (λΔa)/κ2\lambda-\Delta_{a})/\kappa_{2} are depicted versus coupling strength J/k2J/k_{2} for different gain/loss ratios RR, where red dashed and blue solid lines corresponds to R=0.5R=0.5, and green (black) dashed lines correspond to R=2R=2, respectively. Distinct background colors are utilized to differentiate between the unbroken and broken 𝒫𝒯\mathcal{PT}-symmetric regimes.

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 (|tp|2|t_{p}|^{2}) of the output probe field versus optical detuning (δ/ωb\delta/\omega_{b}) for different magnetic field amplitudes (B0B_{0}), keeping the relative phase Φ=0\Phi=0. It is observed that with an increase in the magnetic field amplitude (B0B_{0}), the central absorption peak at the resonance point (δ=ωb\delta=\omega_{b}) 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 (δ=ωb\delta=\omega_{b}). Interestingly, when we further increase the magnetic field amplitude (B0=10×1010B_{0}=10\times 10^{-10}T), the symmetric absorption dips on both sides of the resonance point (δ=ωb\delta=\omega_{b}) become asymmetric. In the blue-sideband regime (δ=1.3ωb\delta=1.3\omega_{b}), 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 (δ=0.7ωb\delta=0.7\omega_{b}), 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 B0=12×1010B_{0}=12\times 10^{-10}T and display the transmission spectrum (|tp|2|t_{p}|^{2}) of the output probe field versus the probe detuning with different relative phase (Φ\Phi). Figure. 4(b) clearly shows that the absorption resulting from the combined action of magnons and polaritons is sensitive to the relative phase (Φ\Phi). Setting Φ=0\Phi=0 leads to an asymmetric transmission profile in the output spectrum. This configuration results in amplification on the right side of the resonance point (δ=ωb\delta=\omega_{b}) and absorption on the left side, as illustrated by the red dotted line in Fig. 4(b). Accordingly, when we set the relative Φ=π\Phi=\pi the switching of the asymmetry line shapes becomes noticeable [see blue solid line in Fig. 4(b)]. In contrast, the central absorption dip at δ=ωb\delta=\omega_{b} remains fairly stable as the relative phase varies from 0 to π\pi.

The physical mechanism behind this attribute is that by adjusting the relative phase (Φ\Phi) 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 (δ=ωb\delta=\omega_{b}), destructive interference occurs, resulting in strong MIA. In contrast, on either side of the resonance point, varying Φ\Phi from 0 to π\pi 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 Φ\Phi on the MMIT profile, we plot the contour map of the transmission spectrum against the optical detuning (δ/ωb\delta/\omega_{b}) and the relative phase Φ\Phi 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 2π-2\pi to 2π2\pi.

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 (|tp|2|t_{p}|^{2}) we depict the graph of transmission rate (|tp|2|t_{p}|^{2}) versus detuning (δ/ωb\delta/\omega_{b}) as shown in Fig. 5. In Fig. 5(a), in the absence of OPA, specifically when the gain (GG) along with phase difference (θ\theta) are both zero, we observed that our system displays an absorption dip at resonance position (δ=ωb\delta=\omega_{b}) as shown by the red dot curve in Fig. 5 (a). As anticipated, introducing the nonlinear gain (G=0.5κ1G=0.5\kappa_{1}) 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 (G=1κ1G=1\kappa_{1}) 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 (GG) 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 (|tp|2|t_{p}|^{2}) versus detuning (δ/ωb\delta/\omega_{b}) for different phases of an OPA (θ=0.5π\theta=0.5\pi, θ=π\theta=\pi and θ=1.5π\theta=1.5\pi) 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 θ=0.5π\theta=0.5\pi. Hence, we can effectively regulate the transmission rate by modifying the OPA settings.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 7:  (a-b) The analytical results of the optical transmission rate (|tp|2|t_{p}|^{2}) versus optical detuning (δ/ωb\delta/\omega_{b}) for various optical tunneling strength JJ. (a) passive-passive (κ1=κ2\kappa_{1}=-\kappa_{2}) (b) passive-active (κ1=κ2\kappa_{1}=\kappa_{2}). (c) Plot of transmission rate (|tp|2|t_{p}|^{2}) versus optical detuning (δ/κ1\delta/\kappa_{1}) with J/κ1=0.8J/\kappa_{1}=0.8. (d) The logarithm of the transmission rate (|tp|2|t_{p}|^{2}) versus optical tunneling strength J/κ1J/\kappa_{1} for passive-passive and passive-active CMM system. With J/κ1=0.8J/\kappa_{1}=0.8 and δ=ωb\delta=\omega_{b}. All other parameters remains same as in Fig. 2(c) except for  G=0.5κ1G=0.5\kappa_{1}Φ=0\Phi=0θ=π/2\theta=\pi/2 and B0=5×1010B_{0}=5\times 10^{-10}T.

III.3 Controllable phase transition and MMIT in CMM system

It is commonly established that systems with 𝒫𝒯\mathcal{PT}-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 𝒫𝒯\mathcal{PT}-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 g1=g2=gmb=0g_{1}=g_{2}=g_{mb}=0. This yields an effective non-Hermitian Hamiltonian, as in [103, 88, 104]:

Heff=\displaystyle H_{eff}= (Δaik1)a^1a^1+(Δa+ik2)a^2a^2\displaystyle\hbar\left(\Delta_{a}-{ik_{1}}\right)\hat{a}_{1}^{\dagger}\hat{a}_{1}+\left(\Delta_{a}+{ik_{2}}\right)\hat{a}_{2}^{\dagger}\hat{a}_{2} (9)
+J(a^1a^2+a^1a^2).\displaystyle+\hbar J\left(\hat{a}_{1}^{\dagger}\hat{a}_{2}+\hat{a}_{1}\hat{a}_{2}^{\dagger}\right).

Equation (9) can be straightforwardly diagonalized via a Bogoliubov transformation [105], and we obtain the complex frequency as:

λ±=Δa+i(κ1κ2)2±4J2(κ1+κ2)22.\lambda_{\pm}=\Delta_{a}+\frac{i(\kappa_{1}-\kappa_{2})}{2}\pm\frac{\sqrt{4J^{2}-\left({\kappa_{1}+\kappa_{2}}\right)^{2}}}{2}. (10)

To investigate the 𝒫𝒯\mathcal{PT}-broken and 𝒫𝒯\mathcal{PT}-symmetric regimes in our system we plot the real and imaginary parts of eigenfrequencies (λΔa)/κ1\lambda-\Delta_{a})/\kappa_{1} against the coupling strength J/k1J/k_{1}, as depicted in Fig. 6. First, we investigate the equal gain and loss scenario (k1=k2k_{1}=k_{2}), as depicted by the blue-solid and red-dashed curves in Figs. 6(a-b). In this scenario, two distinct regimes can be observed: (i) when J>|κ1+κ2|/2J>|\kappa_{1}+\kappa_{2}|/2 in this scenario the eigenfrequencies of the two supermodes exhibit two distinct frequencies (Δa±4J2(k1+k2)2/2\Delta_{a}\pm\sqrt{4J^{2}-(k_{1}+k_{2})^{2}/2}), and an identical linewidth (k1k2)/2k_{1}-k_{2})/2, indicating that our system is in a 𝒫𝒯\mathcal{PT}-symmetry regime (as shown by the white region in Figs. 6((a)-(b))). (ii) when J<|κ1+κ2|/2J<|\kappa_{1}+\kappa_{2}|/2 in this scenario, the two supermodes coalesce and have the identical frequency and different linewidths, indicating that our system is in a broken 𝒫𝒯\mathcal{PT}-symmetry regime (as illustrated by the light pink shaded area in Figs. 6(a-b)).

The critical point, i.e., J=|κ1+κ2|/2J=|\kappa_{1}+\kappa_{2}|/2, where the resonance frequencies and the linewidths of the two supermodes become degenerate, is known as the EP or 𝒫𝒯\mathcal{PT}-phase transition point in our hybrid system. As demonstrated in Figs. 6(a-b), by enhancing the coupling rate JJ, one can observe the transition from a broken to an unbroken 𝒫𝒯\mathcal{PT}-phase. Unlike the 𝒫𝒯\mathcal{PT}-symmetric system with gain and loss, if both cavities are operating in a passive-passive mode with the same loss rate, i.e., κ1=κ2\kappa_{1}=-\kappa_{2}, 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(a-b)].

Now, we examine the scenario where the gain/loss are not perfectly matched i.e., k1=Rk2(R1)k_{1}=Rk_{2}~{}(R\neq 1), where RR shows the gain/loss rate. In Figs. 6(c-d) we consider two cases, i.e., R=0.5R=0.5 (k1=0.5k2k_{1}=0.5k_{2}) and R=2R=2 (k1=2k2k_{1}=2k_{2}). Compared to the equal gain/loss scenario (as in Figs. 6(a-b)), one can clearly see that by varying the gain/loss rate (R=2R=2 or R=0.5R=0.5) causes the EP to shift the right or left, which implies that either a weaker or a stronger coupling strength JJ is necessary for 𝒫𝒯\mathcal{PT}-phase transition.

The simplification to gain and loss in the cavities is done to isolate and study the fundamental 𝒫𝒯\mathcal{PT}-symmetric behavior of the system. This approach allows us to focus on the key aspects of 𝒫𝒯\mathcal{PT}-symmetry—gain and loss—without introducing additional complexities that could obscure the desired physical insights. Once we establish a clear understanding of 𝒫𝒯\mathcal{PT} -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 J=0J=0, which means we work with a single-loss cavity system. Now, our main task is to incorporate the effect of the tunneling rate JJ and observe its impact on the transmission of our CMM system.

Refer to caption
Refer to caption
Figure 8:  Contour (2D) map of transmission rate |tp|2|t_{p}|^{2} vs cavity detuning Δa/ωb\Delta_{a}/\omega_{b} and relative phase Φ\Phi. (a) passive-passive (κ1=κ2\kappa_{1}=-\kappa_{2}) (b) passive-active (κ1=κ2\kappa_{1}=\kappa_{2}). All other parameters remains same as in Fig. 2(c) except for J=1.1κ1J=1.1\kappa_{1}δ=ωb\delta=\omega_{b}G=0.5κ1G=0.5\kappa_{1}Φ=0\Phi=0θ=0\theta=0 and B0=5×1010B_{0}=5\times 10^{-10}T.

Figure 7 (a-b) displays the plot of the optical transmission rate (|tp|2|t_{p}|^{2}) as a function of optical detuning (δ/ωb\delta/\omega_{b}) in the passive-passive and active-passive CMM systems for different optical tunneling strengths JJ, respectively. In Fig. 7(a) for the passive-passive case, we observed that when the optical tunneling strength JJ 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 (J=0.5κ1,0.8κ1J=0.5\kappa_{1},0.8\kappa_{1}) 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, a^±=(a^1±a^2)/2\hat{a}_{\pm}=(\hat{a}_{1}\pm\hat{a}_{2})/\sqrt{2}, 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 JJ 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 JJ [108, 84].

For better comprehension, we provide the analytical results of the optical transmission rate |tp|2|t_{p}|^{2} versus optical detuning δ/κ1\delta/\kappa_{1} in two types of systems by setting the optical tunneling strength J=0.8κ1J=0.8\kappa_{1} 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 𝒫𝒯\mathcal{PT}-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 δ=7κ1\delta=7\kappa_{1}, as illustrated by the red solid line in Fig. 7(c). Figure. 7(d) depicts the logarithm of the transmission rate |tp|2|t_{p}|^{2} vs optical tunneling strength J/κ1J/\kappa_{1}. For the 𝒫𝒯\mathcal{PT}-symmetric system with gain and loss (passive-active scenario), as JJ 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 𝒫𝒯\mathcal{PT}-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 JJ increases beyond the phase transition point, the logarithm of the transmission rate forms an approximately straight line near zero. Therefore, in the 𝒫𝒯\mathcal{PT}-symmetric CMM system, the absorption and amplification of the probe light can be controlled by adjusting the optical tunneling strength JJ.

Correspondingly, Fig. 8 shows the contour (2D) map of transmission rate |tp|2|t_{p}|^{2} vs cavity detuning Δa/ωb\Delta_{a}/\omega_{b} and relative phase Φ\Phi for two different cases. From Fig. 8(a), it is observed that for the passive-passive case, due to the change of relative phase Φ\Phi, a periodic change of the transmission rate takes place, and strong absorption is also observed near the Δa=ωb\Delta_{a}=\omega_{b} point. In Fig. 8(b), the passive-active scenario shows strong amplification at the Δa=ωb\Delta_{a}=\omega_{b} point. The reason for this is that in the passive-active scenario, our 𝒫𝒯\mathcal{PT}-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).

Refer to caption
Refer to caption
Figure 9: Optical group delay τg\tau_{g}(μs\mu s) of the probe light plotted as a function of optical detuning δ/ωb\delta/\omega_{b}: We set (a) The magnetic field amplitude B0=0B_{0}=0B0=6×1010B_{0}=6\times 10^{-10}T,  B0=12×1010B_{0}=12\times 10^{-10}T and relative phase Φ=0\Phi=0. (b) The relative phase Φ=0\Phi=0Φ=π\Phi=\pi and B0=12×1010B_{0}=12\times 10^{-10}T . All other parameters remain same as in Fig. 2(c). except for J=0J=0θ=0\theta=0G=0G=0 and g1=0g_{1}=0.

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]:

τg=darg(tp)dδ.\tau_{g}=\frac{d\arg\left(t_{p}\right)}{d\delta}. (11)

A group delay greater than zero (τg>0\tau_{g}>0) indicates the phenomenon of slow light, whereas a negative group delay (τg<0\tau_{g}<0) 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 (τg\tau_{g}) of our system’s probing field is currently being examined. To illustrate this, in Fig. 9 we plot the optical group delay τg\tau_{g}(μs\mu s) versus optical detuning δ/ωb\delta/\omega_{b}. In Fig. 9(a), we depict the group delay τg\tau_{g}(μs\mu s) of the probe light for differed driving magnetic field amplitudes B0B_{0} while keeping the relative phase at Φ=0\Phi=0. 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 B0=0B_{0}=0, 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  B0=6×1010B_{0}=6\times 10^{-10}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 B0=12×1010B_{0}=12\times 10^{-10}T, the positive group delay rises and reaches a maximum value of 3.8μs3.8\mu s 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 Φ\Phi of the applied field. By comparing the group delay profiles for the relative phases Φ=0\Phi=0, and Φ=π\Phi=\pi [as shown by the red and black solid lines in Fig. 9(b)], we observe that the positive group delay at the probe detuning δ1.002ωb\delta\approx 1.002\omega_{b} corresponding to magnon-induced absorption, significantly decreases and shifts in the case of Φ=π\Phi=\pi. This investigation demonstrates that for a fixed value of the magnetic field amplitude (B0B_{0}), 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.

Refer to caption
Refer to caption
Figure 10: Optical group delay τg\tau_{g}(μs\mu s) of the probe light plotted as a function of pump power Pl(mW)P_{l}(mW), for various values of OPA gains: G=0G=0, G=0.2κ1G=0.2\kappa_{1}, and G=0.3κ1G=0.3\kappa_{1}. We set (a) passive-active system (κ1=κ2\kappa_{1}=\kappa_{2}) (b) passive-passive system (κ2=κ1\kappa_{2}=-\kappa_{1}). All other parameters remains same as in Fig. 2(c). except for J=1.06k1J=1.06k_{1}θ=0\theta=0Φ=0\Phi=0 and B0=5×1010B_{0}=5\times 10^{-10}T.

To further investigate the group delay of the probe light in our CMM system, we plot the group delay τg\tau_{g}(μs\mu s) versus pump power Pl(mW)P_{l}(mW) for various values of OPA gains: G=0G=0 (blue dashed-dot curve), G=0.2κ1G=0.2\kappa_{1} (black dashed curve), and G=0.3κ1G=0.3\kappa_{1} (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 (κ1=κ2\kappa_{1}=\kappa_{2}). In the absence of a degenerate (G=0G=0), the positive group delay of the probe light decreases sharply as the pump power Pl(mW)P_{l}(mW) increases, reaching τg=0\tau_{g}=0, 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 (G=0.2κ1G=0.2\kappa_{1}, and G=0.3κ1G=0.3\kappa_{1}) into the system and increasing the pump power (Pl(mW)P_{l}(mW)) causes the group delay (τg\tau_{g}) 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 (κ1=κ2\kappa_{1}=-\kappa_{2}). 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 (G=0.2κ1G=0.2\kappa_{1}, and G=0.3κ1G=0.3\kappa_{1}), 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 Pl=4(mW)P_{l}=4(mW), 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 D=250μD=250\mum, function as both phonon and magnon resonators, featuring a phonon frequency of ωb=2π×11.42\omega_{b}=2\pi\times 11.42 MHz and the spin density ρ=4.22×1027cm3\rho=4.22\times 10^{27}\textup{cm}^{-3} [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 τg\tau_{g} = -7μs\mu s and τg\tau_{g} = 9μs\mu s. 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 𝒫𝒯\mathcal{PT}-symmetry phase to the 𝒫𝒯\mathcal{PT}-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:

A1+h1=\displaystyle A_{1+}{h_{1}}= iJX1++ig1M1++ig2M2+\displaystyle iJX_{1+}+ig_{1}M_{1+}+ig_{2}M_{2+} (12)
ηaκ1εpeiϕpl2GeiθA1,\displaystyle-\sqrt{\eta_{a}\kappa_{1}}\varepsilon_{p}e^{-i\phi_{pl}}-2Ge^{i\theta}A_{1-}^{*},
A1h2=\displaystyle A_{1-}{h_{2}}= iJX1+ig1M1+ig2M22GeiθA1+,\displaystyle iJX_{1-}+ig_{1}M_{1-}+ig_{2}M_{2-}-2Ge^{i\theta}A_{1+}^{*},
X1+h3=\displaystyle X_{1+}{h_{3}}= iJA1+,\displaystyle iJA_{1+},
X1h4=\displaystyle X_{1-}{h_{4}}= iJA1,\displaystyle iJA_{1-},
B1+h5=\displaystyle B_{1+}{h_{5}}= iF(M2+M2+),\displaystyle iF(M_{2-}^{*}+M_{2+}),
B1h6=\displaystyle B_{1-}{h_{6}}= iF(M2++M2),\displaystyle iF(M_{2+}^{*}+M_{2-}),
M2+h7=\displaystyle M_{2+}{h_{7}}= ig2A1++iF(B1++B1)εmeiϕm,\displaystyle ig_{2}A_{1+}+iF(B_{1+}+B_{1-}^{*})-\varepsilon_{m}e^{-i\phi_{m}},
M2h8=\displaystyle M_{2-}{h_{8}}= ig2A1+iF(B1+B1+),\displaystyle ig_{2}A_{1-}+iF(B_{1-}+B_{1+}^{*}),
M1+h9=\displaystyle M_{1+}{h_{9}}= ig1A1+,\displaystyle ig_{1}A_{1+},
M1h10=\displaystyle M_{1-}{h_{10}}= ig1A1,\displaystyle ig_{1}A_{1-},

with h1=iΔa+iδκ1h_{1}=-i\Delta_{a}+i\delta-\kappa_{1}, h2=iΔaiδκ1h_{2}=-i\Delta_{a}-i\delta-\kappa_{1}, h3=iΔa+iδ+κ2h_{3}=-i\Delta_{a}+i\delta+\kappa_{2}, h4=iΔaiδ+κ2h_{4}=-i\Delta_{a}-i\delta+\kappa_{2}, h5=iωb+iδκbh_{5}=-i\omega_{b}+i\delta-\kappa_{b}, h6=iωbiδκbh_{6}=-i\omega_{b}-i\delta-\kappa_{b}, h7=iΔm2+iδκm2h_{7}=-i\Delta_{m2}^{\prime}+i\delta-\kappa_{m2}, h8=iΔm2iδκm2h_{8}=-i\Delta_{m2}^{\prime}-i\delta-\kappa_{m2}h9=iΔm1+iδκm1h_{9}=-i\Delta_{m1}+i\delta-\kappa_{m1}, and h10=iΔm1iδκm1h_{10}=-i\Delta_{m1}-i\delta-\kappa_{m1}.

After solving Eq. (12) for (A1+A_{1+}) we get the following constants for first-order side-band:
α1=h5h6h7+F2(h6h5)\alpha_{1}=h_{5}h_{6}^{*}h_{7}+F^{2}(h_{6}^{*}-h_{5}),
α2=h2+J2h4+g12h10\alpha_{2}={h_{2}}^{*}+\frac{J^{2}}{{h_{4}}^{*}}+\frac{g_{1}^{2}}{{h_{10}^{*}}}
α3=h8α2+g22\alpha_{3}=h_{8}^{*}\alpha_{2}+g_{2}^{2},
α4=2iGeiθg2+ig2α2\alpha_{4}=2iGe^{-i\theta}g_{2}+ig_{2}\alpha_{2},
α5=ig2α42Geiθα3\alpha_{5}=-ig_{2}\alpha_{4}-2Ge^{-i\theta}\alpha_{3},
α6=ig2h7α2\alpha_{6}=ig_{2}h_{7}\alpha_{2},
α7=ig2α5+ig2α2α3\alpha_{7}=-ig_{2}\alpha_{5}+ig_{2}\alpha_{2}\alpha_{3},
α8=ig2α6+α2α3h7\alpha_{8}=ig_{2}\alpha_{6}+\alpha_{2}\alpha_{3}h_{7},
α9=h8α2α3\alpha_{9}=h_{8}^{*}\alpha_{2}\alpha_{3},
α10=α1α9F2(h6h5)α8\alpha_{10}=\alpha_{1}\alpha_{9}-F^{2}(h_{6}^{*}-h_{5})\alpha_{8},
α11=ig2h5h6α9F2(h6h5)α7\alpha_{11}=ig_{2}h_{5}h_{6}^{*}\alpha_{9}-F^{2}(h_{6}^{*}-h_{5})\alpha_{7},
α12=(α5α10+α6α11α2α3α10)\alpha_{12}=\left(\frac{\alpha_{5}\alpha_{10}+\alpha_{6}\alpha_{11}}{\alpha_{2}\alpha_{3}\alpha_{10}}\right),
α13=h1+J2h3ig1h9ig2α11α10+2Geiθα12,\alpha_{13}=h_{1}+\frac{J^{2}}{h_{3}}-\frac{ig_{1}}{h_{9}}-\frac{ig_{2}\,\alpha_{11}}{\alpha_{10}}+2Ge^{-i\theta}\alpha_{12},
β1=g12α2α2α3\beta_{1}=g_{1}^{2}\alpha_{2}-\alpha_{2}\alpha_{3},
β2=F2(h6h5)β1h5h6α9\beta_{2}=F^{2}(h_{6}^{*}-h_{5})\beta_{1}-h_{5}h_{6}^{*}\alpha_{9},
β3=(α6β2α2α3α10ig1α2α2α3)\beta_{3}=\left(\frac{\alpha_{6}\beta_{2}}{\alpha_{2}\alpha_{3}\alpha_{10}}-\frac{ig_{1}\alpha_{2}}{\alpha_{2}\alpha_{3}}\right),
β4=(ig1β2α102Geiθ)\beta_{4}=\left(\frac{ig_{1}\beta_{2}}{\alpha_{10}}-2Ge^{i\theta}\right).

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:

δ𝐱˙=𝐀𝐱,\delta\dot{\mathbf{x}}=\mathbf{A}\mathbf{x}, (13)

where vector :
𝐱=(δa1,δa1,δm1,δm1,δm2,δm2,δb,δb,δa2,δa2)T\mathbf{x}=(\delta a_{1},\delta a_{1}^{\dagger},\delta m_{1},\delta m_{1}^{\dagger},\delta m_{2},\delta m_{2}^{\dagger},\delta b,\delta b^{\dagger},\delta a_{2},\delta a_{2}^{\dagger})^{T}, in which TT denotes the transpose of a matrix. The matrix 𝐀\mathbf{A} is given by:

𝐀=(ξ12Geiθig10ig2000iJ02Geiθξ20ig10ig2000iJig10ξ300000000ig10ξ4000000ig2000ξ50iFiF000ig2000ξ6iFiF000000iFiF(iωb+κb)0000000iFiF0(iωb+κb)00iJ0000000(iΔaκ2)00iJ0000000(iΔaκ2)).\mathbf{A}=\begin{pmatrix}\xi_{1}&2Ge^{i\theta}&-ig_{1}&0&-ig_{2}&0&0&0&-iJ&0\\ 2Ge^{-i\theta}&\xi_{2}&0&ig_{1}&0&ig_{2}&0&0&0&iJ\\ -ig_{1}&0&\xi_{3}&0&0&0&0&0&0&0\\ 0&ig_{1}&0&\xi_{4}&0&0&0&0&0&0\\ -ig_{2}&0&0&0&\xi_{5}&0&-iF&-iF&0&0\\ 0&ig_{2}&0&0&0&\xi_{6}&iF^{*}&iF^{*}&0&0\\ 0&0&0&0&-iF^{*}&-iF&-\left(i\omega_{b}+\kappa_{b}\right)&0&0&0\\ 0&0&0&0&iF&iF^{*}&0&-\left(-i\omega_{b}+\kappa_{b}\right)&0&0\\ -iJ&0&0&0&0&0&0&0&-\left(i\Delta_{a}-\kappa_{2}\right)&0\\ 0&iJ&0&0&0&0&0&0&0&-\left(-i\Delta_{a}-\kappa_{2}\right)\end{pmatrix}.

Where ξ1=(iΔa+κ1)\xi_{1}=-\left(i\Delta_{a}+\kappa_{1}\right), ξ2=(iΔa+κ1)\xi_{2}=-\left(-i\Delta_{a}+\kappa_{1}\right), ξ3=(iΔm1+κm1)\xi_{3}=-\left(i\Delta_{m_{1}}+\kappa_{m_{1}}\right)ξ4=(iΔm1+κm1)\xi_{4}=-\left(-i\Delta_{m_{1}}+\kappa_{m_{1}}\right)ξ5=(iΔm2+κm2)\xi_{5}=-\left(i\Delta_{m_{2}}^{\prime}+\kappa_{m_{2}}\right), and ξ6=(iΔm2+κm2)\xi_{6}=-\left(-i\Delta_{m_{2}}^{\prime}+\kappa_{m_{2}}\right). The characteristic equation |AΛI|=0|A-\Lambda I|=0 can be reduced to Λn+A1Λn1+.+An1Λ+An=0\Lambda^{n}+A_{1}\Lambda^{n-1}+....+A_{n-1}\Lambda+A_{n}=0 (with n=10n=10) 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 Λ\Lambda 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