Cavity magnomechanics: from classical to quantum

Hybrid quantum systems based on magnons in magnetic materials have made significant progress in the past decade. They are built based on the couplings of magnons with microwave photons, optical photons, vibration phonons, and superconducting qubits. In particular, the interactions among magnons, microwave cavity photons, and vibration phonons form the system of cavity magnomechanics (CMM), which lies in the interdisciplinary field of cavity QED, magnonics, quantum optics, and quantum information. Here, we review the experimental and theoretical progress of this emerging field. We first introduce the underlying theories of the magnomechanical coupling, and then some representative classical phenomena that have been experimentally observed, including magnomechanically induced transparency, magnomechanical dynamical backaction, magnon-phonon cross-Kerr nonlinearity, etc. We also discuss a number of theoretical proposals, which show the potential of the CMM system for preparing different kinds of quantum states of magnons, phonons, and photons, and hybrid systems combining magnomechanics and optomechanics and relevant quantum protocols based on them. Finally, we summarize this review and provide an outlook for the future research directions in this field.


Introduction
Quantum information science and technology, exploiting the principles of quantum mechanics, have great potential to cause revolutionary advances in fields of science and engineering in the twenty-first century.The manipulation and processing of quantum information and the realization of multifunctional quantum technologies require hybrid quantum systems (HQSs) as the implementation platform, which combines different physical components with individual strengths and complementary functionalities [1,2,3].
Up to now, although there have been exhaustive reviews on the interactions between magnons and microwave cavity photons, optical cavity photons, and superconducting qubits [4,8,9], and on the configurations, circuits and materials for realizing the HQS based on magnons [5,7,10,12], there is a lack of a thorough review on the field of cavity magnomechanics (CMM), which studies the interactions among magnons, microwave cavity photons, and vibration phonons [23,24,25,26].The present review exactly fills this gap, which includes underlying theories of the magnomechanical coupling, up-to-date experimental observations, and abundant theoretical studies/protocols, covering both classical and quantum phenomena in the field of CMM.In particular, the nonlinearity provided by the magnetostrictive interaction has been identified as an essential element for the experimental observations of magnomechanically induced transparency, magnomechanical dynamical backaction, magnon-phonon cross-Kerr nonlinearity and magnonic frequency combs, and a prerequisite in many theoretical proposals for preparing macroscopic quantum states of a sub-millimeter sized YIG sphere involving a large number of magnons and phonons, and nonclassical states, e.g., entangled and squeezed states, of microwave fields.All of these make the CMM system a unique platform for exploring rich physics in both classical and quantum regimes, which find broad important applications in macroscopic quantum studies, quantum information processing, quantum technologies, quantum memories, quantum metrology, spectroscopy, nonlinear researches, etc.
The review is organized as follows.Section 2 describes the fundamental theory of the magnomechanical coupling and its fully quantized interaction Hamiltonian, which lays the theoretical foundation for various applications using the magnomechanical systems.Section 3 introduces a series of representative classical phenomena in the field of CMM, of which most have been experimentally demonstrated at room temperature.Section 4 discusses a number of theoretical proposals for preparing different kinds of quantum states in the system, e.g., entangled states, squeezed states and quantum steering of magnons, phonons and photons, fully exhibiting its potential as a HQS.Section 5 introduces hybrid systems combining both magnomechanics and optomechanics and relevant quantum protocols based on them.Finally, Section 6 summarizes this review, outlines important future research directions in the field of CMM by applying the concepts of non-Hermiticity, PT symmetry, and exceptional points to the CMM system, and briefly introduces other magnomechanical platforms beyond the frameworks set by Refs.[4,8,9].

Fundamental theory of the magnomechanical coupling
The magnomechanical coupling describes the interaction between magnetization and elastic strain of a magnetic material.Depending on the distance between magnetic atoms (or ions), there are three kinds of interactions: the spin-orbital interaction, the exchange interaction, and the magnetic dipole-dipole interaction [33].For a cubic crystal, the magnomechanical coupling is described by the magnetoelastic energy density, given by [34,35] where B 1 and B 2 are the magnetoelastic coupling coefficients, M S is the saturation magnetization and M x,y,z are the magnetization components.ϵ i j = 1 2 ∂u i /∂l j + ∂u j /∂l i is the strain tensor of the magnetic crystal, with coordinate indices i, j ∈ {x, y, z}, and u x,y,z are the components of the displacement vector ⃗ u.
A magnon is viewed as a quantized spin wave.According to the Holstein-Primakoff approximation [36], the magnetization can be quantized via and where m (m † ) is the annihilation (creation) operator of the magnon mode, V is the volume of the magnetic crystal, and γ is the gyromagnetic ratio.Note that equation ( 3) is obtained by taking the first-order Taylor series expansion of 1 − (M 2 x + M 2 y )/M 2 S .Substituting equations ( 2) and (3) into the magnetoelastic energy density f me and integrating over the whole volume of the crystal, the semiclassical magnetoelastic Hamiltonian can be obtained as the sum of The magnetoelastic Hamiltonian, H me = H 1 + H 2 , implies different types of magnon-phonon interactions depending on the relation of their frequencies.For specific magnon and phonon frequencies, a certain type of coupling can be dominant in the magnetoelastic coupling, while other coupling mechanisms are negligible.This can be more clearly seen in the fullyquantized interaction Hamiltonian by further quantizing the strain displacement [37].
The magnetoelastic displacement can be decomposed and expressed as the superposition form of ⃗ u(x, y, z) = n,m,k d (n,m,k) ⃗ χ (n,m,k) (x, y, z), (6) where ⃗ χ (n,m,k) (x, y, z) represents the dimensionless displacement eigenmode and d (n,m,k) is the corresponding amplitude with the mode indices (n, m, k).The displacement amplitude d (n,m,k)  can be quantized as where d (n,m,k)   zpm denotes the amplitude of the zero-point motion, and b n,m,k and b † n,m,k are the annihilation and creation operators of the corresponding phonon mode.
The fully-quantized magnetoelastic Hamiltonian can be derived by quantizing the strain displacement in the semiclassical Hamiltonians of equations ( 4) and (5) [37], which accounts for different types of magnon-phonon interactions dependent on their frequencies.Specifically, for the case of the phonon frequencies being much lower than the magnon frequency ω (n,m,k) b ≪ ω m (the case of a large-size magnetic crystal), by neglecting the nonresonant fast-oscillating terms, one obtains the dominant dispersive-type interaction where the dispersive coupling strength Note that the above Hamiltonian is derived under the condition of low-lying magnon excitations, where the second-order term of the magnon excitation m † m is neglected [37].Such a dispersive coupling is the dominant magnon-phonon interaction for sub-millimeter sized YIG spheres [23,24,25,26], which provides essential nonlinearity that is exploited in many proposals for preparing quantum states, as will be discussed in Section 4.
For the case where the phonon mode is nearly resonant with the magnon mode, ω (n,m,k) b ≃ ω m , after neglecting the nonresonant fast-oscillating terms, one gets the following dominant linear interaction where the coupling rate Similarly, the higher-order term has been neglected under the low-lying excitations.This magnon-phonon linear coupling is typically the case for high-frequency phonons, e.g., in the gigahertz range.In this regime, the magnon-phonon strong coupling can be achieved forming magnon-phonon polarons [38].Though strong, such a beam-splitter coupling itself can only realize state-swap interaction between magnons and phonons, but cannot create entanglement and squeezing, as will be discussed in Section 4.
Apart from the aforementioned dispersive and linear interactions, a magnon parametric amplification (PA) interaction becomes dominant when the phonon frequency is about twice of the magnon frequency, ω (n,m,k) b ≃ 2ω m .This leads to the following interaction with the coupling rate Such magnon PA can be used to generate entangled pairs of magnons, in analogy to the generation of entangled photon pairs by optical parametric down-conversion.
To sum up, by fully quantizing the magnetoelastic Hamiltonian, different effective magnomechanical Hamiltonians are derived for different relations between the magnon and phonon frequencies: specifically, i) the dispersive interaction for the phonon frequency being much lower than the magnon frequency, ω (n,m,k) b ≪ ω m ; ii) the linear (beam-splitter) interaction for nearly resonant magnons and phonons, ω (n,m,k) b ≃ ω m ; and iii) the magnon PA interaction when ω (n,m,k) b ≃ 2ω m .Note that all the effective interactions are obtained by neglecting nonresonant fast-oscillating terms, which are good approximations when the magnon frequency and (or) the phonon frequency are (is) much larger than their coupling strength and dissipation rates [37].
Figure 1.The CMM experimental device consists of a three-dimensional microwave cavity made from oxygen-free copper and a YIG sphere supporting both a magnon mode and a mechanical vibration mode.The YIG sphere is placed near the maximum magnetic field of the cavity mode and can be either glued to the end of a silica fiber [23] (a) or placed in a horizontal glass capillary [25,26] (b)-(c).A uniform external magnetic field, provided by a set of permanent magnets (b), is applied to bias the YIG sphere for the magnon-photon coupling.(d) Interactions between the different modes.The magnon mode couples to the cavity magnetic field via the magnetic-dipole interaction with the coupling strength g ma , and to the mechanical vibration mode via the dispersive magnetostrictive interaction with the bare magnomechanical coupling strength g mb .Figures (a)-(d) are adapted from Refs.[23,25,26].

Classical cavity magnomechanics
The exploration of the CMM has started from the experimental studies on classical phenomena, and so far, all the CMM experiments have been operated at room temperature.The typical CMM experimental device is shown in Figure 1, which consists of a threedimensional microwave cavity machined from oxygen-free copper and a highly polished single-crystal YIG sphere supporting both a magnon mode and a mechanical vibration mode.The YIG sphere with the diameter of hundreds of micrometers is placed near the maximum microwave magnetic field of the cavity mode.The YIG sphere can be either glued to the end of a silica fiber [23], or placed (free to move) in a horizontal glass capillary to reduce the mechanical damping [25,26].A uniform external magnetic field is applied to bias the YIG sphere for the magnon-photon coupling.
The Hamiltonian of the CMM system reads [23,24,25,26] where a (a † ) and m (m † ) ([ j, j † ] = 1, j = a, m) are the annihilation (creation) operators of the cavity and magnon modes, respectively, q = (b 2 are the dimensionless position and momentum of the vibration mode (modeled as a mechanical oscillator and [q, p] = i), and ω k (k = a, m, b) are their resonance frequencies.In typical CMM experiments, ω a and ω m are nearly resonant and around ∼10 GHz [23,25,26,39].The magnon frequency can be tuned by varying the external bias magnetic field H via ω m = γH, where the gyromagnetic ratio γ/2π = 28 GHz/T.The magnon mode couples to the cavity magnetic field via the magnetic-dipole interaction with the coupling strength g ma , which can be (much) stronger than the cavity and magnon dissipation rates κ a and κ m , leading to the hybridization of the two modes forming two polaritons [13,14,15,16,17,18].The magnon mode further couples to a mechanical vibration mode via magnetostriction.For a YIG sphere with the diameter in the range of 10 2 -10 3 µm, the mechanical frequency is in the megahertz range [23,25,26,39], which is much lower than the magnon frequency in the gigahertz range, thus promising a dominant dispersive magnomechanical coupling [23,37,40] (Section 2).The bare magnomechanical coupling g mb is typically weak, but the effective coupling can be significantly enhanced by driving either the magnon or the cavity with a strong microwave field.For driving the magnon, e.g., via a loop antenna [39], H dri /ℏ = iΩ m † e −iω 0 t − me iω 0 t , where the Rabi frequency Ω = √ 5 4 γ √ NB 0 [24] denotes the coupling strength between the magnon and the drive magnetic field with frequency ω 0 and amplitude B 0 .N = ρV is the total number of spins with ρ = 4.22 × 10 27 m −3 being the spin density of the YIG and V being the volume of the YIG crystal.For driving the cavity, H dri /ℏ = iE a † e −iω d t − ae iω d t , where E denotes the cavity-drive coupling strength and ω d is the frequency of the drive field.
In what follows, we introduce a series of classical phenomena that have been either experimentally demonstrated or theoretically studied, including magnomechanically induced transparency and absorption, magnomechanical dynamical backaction, magnon-phonon cross-Kerr nonlinearity and mechanical bistability, etc.

Magnomechanically induced transparency and absorption
The radiation-pressure-like dispersive magnomechanical interaction predicts the magnomechanically induced transparency (MMIT) and absorption (MMIA) [23], in analogy to the optomechanically induced transparency (OMIT) and absorption (OMIA) [41,42,43].It describes an interference phenomenon that the transmission/reflection of the weak probe signal is controlled by a strong drive field when the resonance condition is met.In a typical measurement setup of the MMIT (MMIA) [23], the cavity is driven by both a strong microwave field with frequency ω d and amplitude ε d and a weak probe signal with frequency ω s and amplitude ε s .The dynamics of the system is described by the following semiclassical Heisenberg equations: where κ e is the external decay rate of the cavity, and γ b is the mechanical damping rate.Note that the noise terms are removed and the operators are replaced with their expectation values.Equivalently, the position and momentum (q, p) are replaced with (b, b * ) to describe the mechanical mode.
The Heisenberg equations (15) are nonlinear due to the dispersive magnomechanical coupling.They can be solved by linearizing the dynamics of the system around the steadystate values by writing O = Ō + δO (O = a, m, b), where (ā, m, b) are the large mean values and (δa, δm, δb) are the small perturbations.Because the probe signal is much weaker than the drive field, the probe signal can be neglected in deriving the solutions of the mean fields, which are given by By omitting high-order perturbation terms, one obtains the linearized equations for the perturbations (δa, δm, δb), which, in the frame rotating at the drive frequency, are The above equations can be conveniently solved in the frequency domain by taking the Fourier transform of each equation.Having achieved the solution of δa(ω), one then obtains the cavity output field δa out (ω) using the input-output relation δa out (ω) = ε s + i √ 2κ e δa(ω) [23], and thereby the reflection spectrum of the probe signal, i.e., r(ω) ≡ δa out (ω)/ε s .The measured reflection spectra as a function of the probe-drive detuning ∆ sd = ω s − ω d for a red (blue)detuned drive field is shown in Figure 2(a) (Figure 2(b)).The sharp Lorentzian-shaped peak (dip) at ∆ sd = ω b (∆ sd = −ω b ) in the spectra is evidence of the magnomechanical interaction, i.e., the MMIT (MMIA) in Figure 2(a) (Figure 2(b)) [23].Note that in each driving case, the driving frequency ω d is fixed and the mechanical frequency ω b /2π ≈ 11.42 MHz in the experiment [23].
The strong coupling between the cavity and magnon modes leads to two hybridized polariton modes A ± , formed via A + = a cos θ + m sin θ and A − = −a sin θ + m cos θ, with θ = 1 2 arctan 2g ma ω a −ω m .It indicates that the weight of the cavity (magnon) component in each polariton can be varied by altering θ, realized by changing the coupling g ma and/or the detuning ω a − ω m .Since the cavity component has no coupling with the mechanical mode, the interaction between each polariton and the mechanics is fully determined by the interaction between the magnon component and the mechanics, i.e., the polariton-mechanics coupling is also dispersive [44,39].In the spectra of Figure 2 Similarly as the OMIT (OMIA), the MMIT (MMIA) can also cause characteristic changes in dispersion in the range of the transparency (absorption) window, which can be used to manipulate the group velocity of the probe signal.Slow light effects are predicted to occur in the CMM system [45].Long-lived slow light (group delay of the order of millisecond) can be achieved and the group delay can be continuously adjusted by tuning the bias magnetic field [45].Multi-window MMIT, Fano resonances, and slow-to-fast light conversion have been studied in a two-YIG-spheres CMM system [46], where the group delay of the probe signal can be greatly enhanced compared to the single-YIG-sphere case.It is also suggested that magnon squeezing plays an active role in enhancing and controlling the group delay of the probe field [47] and phase control of the slow light effect is possible [48].

Magnomechanical dynamical backaction: mechanical cooling and amplification, and magnonic spring effect
The dispersive magnomechanical interaction also promises dynamical backaction effects on the mechanical mode (similar as those in optomechanics [49]), such as mechanical cooling and amplification, and the "magnonic spring" effect, i.e., the mechanical frequency modified by the magnetostrictive interaction [25].The characterization of these effects requires the measurements of the frequency and linewidth of the mechanical mode.The mechanical frequency and linewidth, as well as the bare magnomechanical coupling g mb , can be extracted by fitting the MMIT signal [23,26], or adopting a microwave homodyne detection scheme, which consists of a microwave generator as the local oscillator, an IQ mixer, and an analogto-digital converter [25].They play a similar role as in the optical homodyne detection of the mechanical state in optomechanics [50].
To enhance the magnomechanical coupling, the system is operated to satisfy the triple- resonance condition [23], where the splitting of the two hybrid modes (i.e., two polaritons) in the spectrum matches the mechanical frequency, ω + − ω − ≡ 4g 2 ma + ∆ 2 am = ω b , with ∆ am ≡ ω a − ω m .In this case, both the drive field and the scattered mechanical sideband can be on resonance with the two polaritons, respectively, leading to drastically enhanced magnomechanical coupling.In the experiment [25], this is implemented by setting the normal mode splitting (NMS) for the resonant case (∆ am = 0) to be slightly smaller than the mechanical frequency, i.e., 2g ma < ω b .Then by slightly detuning the magnon frequency from the cavity frequency (∆ am 0), the NMS can be tuned to exactly match the mechanical frequency.When the triple-resonance condition is met, one can choose to resonantly drive the lower (higher)-frequency polariton, such that the anti-Stokes (Stokes) mechanical sideband resonates with the higher (lower)-frequency polariton, resulting in the cooling (amplification) of the mechanical motion [25].This corresponds to an increased (a reduced) effective mechanical damping rate γ eff b = γ b +δγ b , with δγ b > 0 (δγ b < 0) being the magnomechanically induced additional damping rate, of which the expression is where G mb is the effective magnomechanical coupling strength and , with χ m (ω) (χ a (ω)) being the natural susceptibility of the magnon (cavity) mode.The above mechanism works optimally in the resolved sideband limit κ ± ≪ ω b [51,52,53], with κ ± being the polariton linewidths, which is typically well satisfied in the CMM experiments using YIG [23,25,26,39], benefiting from the low magnonic dissipation rate of the YIG.
Figure 3(a) shows the reduced effective damping rates γ eff b (mechanical linewidths) versus the drive power for three environmental conditions of the YIG sphere placed within a glass capillary, i.e., in air (blue squares), a low pressure (∼15 Torr) of pure helium gas (purple triangles), and vacuum (red circles), corresponding to the natural mechanical damping rate γ b of 247 Hz, 98 Hz, and 59 Hz, respectively.Figure 3(a) is obtained by applying a drive on resonance with the higher-frequency polariton, such that the Stokes scattering is resonantly enhanced.As the drive power increases, a threshold is reached at which the effective damping rate γ eff b approaches zero.By further increasing the power, the mechanical oscillations grow exponentially in time and are ultimately limited by higher-order nonlinear effects.This parametric instability is referred to as phonon lasing.The inset of Figure 3(a) shows the mechanical power spectrum above (orange) and below (green) the threshold power.When the drive is on resonance with the lower-frequency polariton, the anti-Stokes scattering is resonantly enhanced, leading to a broadened mechanical linewidth (or an increased γ eff b ) (green crosses in Figure 3(b)).As a comparison, the red circles in Figure 3(b) correspond to the mechanical amplification process as in Figure 3(a).Due to the small magnomechanical cooperativity (slightly greater than 1) achieved in the experiment, the cooling effect is weak and the effective mean phonon number neff b is only reduced by a few times from the initial large mean thermal phonon number nb .
The mechanical motion can be cooled into its quantum ground state when the cooperativity is significantly improved and the system is at a low bath temperature [54,55].This is a prerequisite for preparing various quantum states in the CMM system [24,56,57,58].Further, it is shown that magnon squeezing can significantly suppress the magnomechanical Stokes scattering and thus becomes particularly useful in realizing ground-state cooling in the unresolved-sideband regime [55], where the sideband cooling mechanism becomes inefficient.This is the case for the CMM system using other ferromagnetic materials, e.g., CoFeB [59], which exhibits much stronger magnetostriction but also much larger magnon dissipation.As for mechanical amplification, phonon laser has also been theoretically studied [60,61,62,63,64], e.g., nonreciprocal [61] and exceptional-point-engineered [64] phonon laser are suggested.
Another magnomechanical backaction effect is called the "magnonic spring" effect, in analogy to the "optical spring" in optomechanics [65].It refers to the phenomenon that the magnomechanical interaction induces a frequency shift of the mechanical mode.The mechanical frequency shift is given by The "magnonic spring" effect is revealed in Figure 3(c) for a mechanical resonance of ω b /2π = 12.6278 MHz in the partial pressure of helium.It shows that the mechanical oscillator is spring-softened (δω b < 0) for a red-detuned drive and spring-hardened (δω b > 0) for a blue-detuned drive.Though weak, this is the first observation of the "magnonic spring" effect [25].Dynamical backaction can, however, be detrimental for some specific applications.A magnomechanical backaction evading experiment has been demonstrated [66].A microwave drive tone, tuned between the two polaritons, simultaneously activates both the Stokes and anti-Stokes scatterings.By carefully balancing the scattering rates of the two processes, dynamical backaction effects are eliminated, confirmed by the measured drive-powerindependent mechanical linewidth [66].

Magnomechanical cross-Kerr nonlinearity and mechanical bistability
As introduced above, the "magnonic spring" effect causes a frequency shift of the mechanical mode.This frequency shift is typically small, at the level of 0.1 kHz as observed in the experiment [25] under a relatively low drive power.However, when the drive power becomes strong, e.g., higher than 10 dBm, a much larger mechanical frequency shift is observed than that is predicted by the dispersive magnomechanical coupling [26].This implies that some higher-order nonlinear coupling between magnons and phonons is activated by the strong drive, resulting in a considerable mechanical frequency shift.
This nonlinear coupling is termed as the magnon-phonon cross-Kerr nonlinearity [26], which originates from the magnetoelastic coupling by including the second-order terms in the nonlinear strain tensor [67,59] where u i are the components of the displacement vector, and l i = i (i = x, y, z).The first-order terms lead to the conventional dispersive (radiation-pressure-like) interaction Hamiltonian, For a not strong drive field, the second-order terms are negligible [23,25], but they can no longer be neglected when the drive becomes sufficiently strong, yielding considerable magnon-phonon cross-Kerr nonlinearity with the interaction Hamiltonian ℏK cross m † mb † b [26], where K cross is the cross-Kerr coefficient.
The dispersive magnomechanical coupling gives rise to an effective susceptibility of the mechanical mode [26] where χ b (ω) is the natural susceptibility of the mechanical mode, but depends on the modified mechanical frequency ωb = ω b + K cross |M| 2 (|M| 2 is the magnon excitation number), which includes the cross-Kerr-induced frequency shift.The effective coupling G mb = g mb M, and as defined in equation ( 18), but here the magnon frequency in χ m (ω) is modified to be ωm = ω m +2K m |M| 2 , which includes the frequency shift dominantly contributed from the magnon self-Kerr nonlinearity [68,69] caused also by the strong drive, and K m is the magnon self-Kerr coefficient.
The effective mechanical susceptibility yields a mechanical frequency shift (i.e., the "magnonic spring" effect) where we write together the frequency shift induced by the cross-Kerr nonlinearity (last term).
In addition, it causes a mechanical linewidth change Note that the above equation is different from equation (18), because χ m (ω) is redefined with the modified magnon frequency ωm .Clearly, the linewidth change is only caused by the dispersive coupling, which can be positive (cooling) or negative (amplification) for a red-or blue-detuned drive, as discussed in Section 3.2.This differs from the frequency shift where both the dispersive and cross-Kerr couplings contribute, and the latter becomes dominant when the drive is strong.This is clearly seen in Figure 4(a) (Figure 4(c)) for a red (blue)-detuned drive.As the power grows, the cross-Kerr-induced frequency shift (green curve) increases rapidly, whereas the dispersive-coupling-induced frequency shift (red curve) increases slowly and is generally very small, in the order of 0.1 kHz.In both the driving situations, the cross-Kerr gives rise to a negative frequency shift because of a negative cross-Kerr coefficient K cross /2π = −5.4pHz in the experiment [26].By contrast, the frequency shifts by the spring effect are opposite under the red-and blue-detuned drives, in agreement with the finding of Ref. [25].
The mechanical linewidth change δγ b versus the drive power for a red (blue)-detuned drive is shown in Figure 4(b) (Figure 4(d)).As expected, a red (blue)-detuned drive results in a dominant anti-Stokes (Stokes) scattering and thus an increased (a reduced) mechanical linewidth.Note that the triple-resonance condition is not fulfilled in the experiment [26].This is because due to the magnon self-Kerr effect, the magnon frequency shifts and thus the splitting of the two polaritons varies by changing the drive power.Therefore, the tripleresonance condition is unlikely to be met as the drive power changes (for fixed cavity frequency ω a and coupling g ma ), but can be achieved for a fixed power.
In Figure 4, both the mechanical frequency and linewidth show a bistable feature with the forward and backward sweeps of the drive power for the two driving situations.This is because both of them originate from the bistable magnon excitation number |M| 2 (by varying the drive power) induced by the magnon self-Kerr effect [26,69].For both the dispersivecoupling-induced frequency shift and linewidth change, the bistability of |M| 2 is mapped to the magnon frequency ωm , then to the polariton susceptibility χ ma (ω), and eventually to the mechanical frequency and linewidth (cf.equations ( 22) and ( 23)).

Strong-coupling CMM: polariton-mechanics normal-mode splitting
Due to the large frequency mismatch between the magnon and mechanical modes, their coupling is typically weak, leaving the CMM experiments [23,25,26] in the weak-coupling regime.This considerably limits diverse applications of the CMM system, such as quantum states preparation [24,56,57,58].Here, the weak coupling can refer to the following two situations, depending on whether the cavity and magnon modes are strongly coupled: i) the effective magnomechanical coupling G mb is smaller than the magnon decay rate κ m ; or ii) the (cavity-magnon) polariton-mechanics coupling G +/−,b is smaller than the polariton decay rate κ +/− .Since the mechanical damping rate is much smaller, the strong coupling is essentially determined by the relation between the coupling G mb (G +/−,b ) and the magnon (polariton) decay rate κ m (κ +/− ).
The natural way of achieving the strong coupling is to improve the effective coupling strength G mb or G +/−,b , e.g., using a strong drive field.However, there are restrictions on such improvement: a strong drive field may cause instability [70] and unwanted nonlinear effects, such as magnon self-Kerr [68,69] and magnon-phonon cross-Kerr [26] nonlinearities.One may consider using other ferromagnetic materials, e.g., CoFeB [59], which exhibit stronger magnetostriction but also significantly increased magnon dissipation.This may end up with a reduced magnomechanical cooperativity and thus may not be an effective means.
An effective approach is offered [39] by significantly reducing the dissipation rate of the polariton mode exploiting coherent perfect absorption (CPA) [71,72,73,74,75,76].Specifically, the cavity is fed by a microwave field to realize an effective gain of the cavity mode.By operating the system at the (cavity) gain-(magnon) loss balance determined by the CPA conditions, i.e., κ ′ a ≡ κ e − κ int = κ m for the resonant case ω m = ω a [39], with κ e (κ int ) being the external (internal) decay rate of the cavity, the decay rate of one of the two polaritons (e.g., the higher-frequency polariton) can be reduced to nearly zero.The experiment [39] realizes the simplest CPA with a single port [71,72,76].Due to the destructive interference between the input and output fields at the cavity port, the amplitude of the cavity output field turns out to be zero, leading to a vanishing dip (−82 dB) in the reflection spectrum at the CPA frequency ω CPA = ω + (Figure 5(a)), which is a typical feature of the CPA.The smallest decay rate of the higher-frequency polariton κ + /2π = 140 Hz is achieved in the experiment [39], which is substantiated by the calculation of the Wigner time delay [77,78] based on the reflection spectrum of the probe field.The polariton decay rate is reduced by four orders of magnitude from the normal value ∼ 1 MHz in the cavity magnonic experiments [13,14,15,16,17,18].
Because of the significantly reduced polariton decay rate and the intrinsically small mechanical damping rate, the polariton-mechanics strong coupling G +,b > κ + , κ b is achieved.Note that the system is actually in the triple cavity-magnon-phonon strong coupling regime, since the polaritons are formed by strongly-coupled cavity photons and magnons.This leads to the NMS at the mechanical sideband in the spectrum being resonant with the higher- frequency polariton [39], as shown in Figure 5(b)-(c).In the experiment, a red-detuned drive field is applied at frequency ω d /2π = 10.04762GHz and the scattered (anti-Stokes) mechanical sideband is thus at ω d + ω b = 2π × 10.05857 GHz, where the mechanical frequency ω b /2π = 10.9485MHz.By varying the drive power and exploiting the magnon self-Kerr effect, the frequency of the polariton can be continuously tuned (Figure 5(c)).When the polariton frequency is approaching and then passing through the mechanical sideband, the normal-mode spectrum shows an anti-crossing feature due to the polariton-mechanics strong coupling.The polariton-mechanics cooperativity, C +,b = G 2 +,b / (κ + κ b ), can be up to C +,b ≈ 4 × 10 3 when the polariton nearly resonates with the mechanical sideband (inset of Figure 5(b)), which is improved by three orders of magnitude than the previous CMM experiments [23,25,26].Similar NMS at the mechanical sideband has been demonstrated in optomechanical systems, as a sign of the optomechanical strong coupling [79,80,81].
The remarkably improved cooperativity can significantly boost the mechanical cooling efficiency and lower the threshold for achieving phonon lasers [25], and makes it possible to prepare magnonic and mechanical quantum states [24,57,58], if further placing the system at a low-temperature environment giving rise to the quantum cooperativity C Q +,b ≡ C +,b /n b much greater than 1.The experiment [39] paves the way towards coherent control and measurement of the quantum states of phonons, photons and magnons.

High-order sidebands and frequency combs
The MMIT and MMIA can be well described by the linearized magnomechanical interaction.Taking into account the nonlinearity of the magnetostriction, it can be used to generate magnonic high-order sidebands (similar to optical high-order sidebands created in optomechanical systems [43]), with frequencies ω d + nω b , where ω d is the driving frequency and n is an integer denoting the order of the magnonic sideband.The efficiency of the magnonic second-order sideband generation [82] can be calculated analytically within the perturbative regime by considering that the generation of the second-order sideband is seeded by the drive field at the first-order sideband.Magnonic high-order sidebands find many applications, e.g., in magnonic information processing and magnonic frequency combs (MFCs).
The concept of the MFC can be regarded as a magnonic counterpart to the optical frequency comb (OFC), which contains a set of discrete modes that are equally spaced in frequency (Figure 6(a)).As is known, OFCs are vital for and have revolutionized the fields of high-precision frequency metrology and spectroscopy.The extension of frequency comb techniques [83,84] to novel physical carriers, such as phonons [85] and magnons [86,87], is likely to open up new vistas for applications, especially in spectroscopy and sensing.Nevertheless, the efficiency of magnonic sideband generation is very weak due to the intrinsically weak magnomechanical interaction.Hence, strategies are highly needed to enhance the generation of magnonic high-order sidebands.
The sideband generation can be enhanced via phonon lasing (or giant mechanical oscillation).Along this line, the generation of MFCs based on resonantly enhanced magnetostrictive interaction has been proposed [88].In Ref. [88], a YIG sphere is placed in the vicinity of a port, which can load a microwave drive field via a loop antenna (Figure 6(b)).A bias magnetic field is applied to establish the coupling between the drive field and the magnon modes.Magnonic and mechanical dynamics can be substantively modified due to the magnetostrictive effect, which results in degenerate and non-degenerate four-wave mixing and frequency conversion of magnons.A threshold for the amplitude of the drive field is obtained, and when the drive field is above the threshold, it leads to unstable dynamics where the magnon number grows over time, yielding sufficiently strong magnetostrictive nonlinearity to generate wave mixing and high-order sidebands (Figure 6(c)).This can be understood from the perspective of phonon lasing.The magnetostrictive interaction leads to the modification of the mechanical damping rate.Below the threshold, the effective mechanical damping rate is positive and the magnonic dynamics is stable.Whereas the drive over the threshold results in a negative effective damping rate and thus the amplification of the mechanical motions manifested as self-induced oscillations (Figure 6(d)).
Experimental generation of the MFCs based on the magnetostrictive interaction has been recently realized [89].In the experiment (Figure 6(e)), a YIG sphere with the diameter of 623.6 µm is placed in a bias magnetic field, and a microwave field with frequency about 4.725 GHz is applied via an antenna to excite and drive the magnon mode.For a weak microwave pump with the power below the threshold of 17.4 dBm, the magnonic spectrum contains only the frequency component of the pump field, with no comb lines.When the pump power reaches the threshold, an MFC emerges with a tooth spacing of ω b = 10.08 MHz (Figure 6(f)) and the number of the comb lines increases steadily with the pump power.
The efficiency of sideband generation can also be enhanced by adopting a two-tone microwave drive [90] or the cavity-magnon dissipative coupling [91].Nonreciprocal sideband generation has been proposed in a spinning microwave magnomechanical system, where a YIG sphere is coupled to a spinning resonator [92].These proposals provide a pathway for the generation of robust MFCs that may be beneficial for on-chip microwave isolation devices and magnon-based precision metrology.

Thermometry
Many quantum information tasks and quantum technologies require low-temperature environments.In cryogenic environments, accurate thermometry can be difficult to implement and often requires calibration to an external reference.To solve this problem, Ref. [93] proposes a thermometry method using the CMM system at a low bath temperature.The thermometry is based on the measurement of two correlation spectra of the cavity output field, i.e., the phase-phase autocorrelation and the amplitude-phase correlation spectra.The temperature of the low-frequency phonon mode (in equilibrium with the bath) is determined by the ratio of the two correlation spectra, i.e., the so-called thermometric relation [93].The thermometric relation as a function of the bath temperature is shown in Figure 7(a).It should be noted that the thermometric relation is independent of experimental parameters, e.g., coupling strengths and decay rates.The thermometric measurement is most accurate at low temperatures when the thermal microwave/magnon occupation is less than unity, implying that higher microwave/magnon resonance frequencies are preferred, cf. Figure 7(b).
Such CMM-based thermometry shows some advantages compared with that employing optomechanical systems.The heating due to optical absorption limits the effective lower temperature range of optical quantum-correlation thermometry to ∼10 K [94,95].By contrast, the use of microwave photons, which cause minimal heating of the mechanical system due to their low energy, allows this protocol to be effectively used at temperatures below ∼500 mK.

Quantum cavity magnomechanics
The nonlinear magnetostrictive interaction provides necessary nonlinearity for preparing various quantum states, e.g., entangled and squeezed states, in the CMM system.By further exploiting, e.g., external quantum fields, the reservoir engineering technique, and measurements on the microwave field, quantum states of the system can be created or enhanced.In general, the preparation of quantum states requires the system to be at a low bath temperature to reduce the detrimental thermal noises.Otherwise, a cooling mechanism for the mechanical motion should be present.

Quantum entanglement
Many proposals indicate that the CMM system holds the potential to prepare entangled states of magnons and phonons of macroscopic YIG crystals.These macroscopic entangled states are relevant to the fundamental studies of macroscopic quantum mechanics [96], the boundary between the quantum and classical worlds [97,98], decoherence theories in macroscopic scale [99,100], etc. Further, the CMM system can also be exploited to prepare entangled states of microwave fields, which find wide applications in quantum information processing [4,8], quantum teleportation [101], and quantum metrology [102].

Magnon-photon-phonon entanglement
The Hamiltonian ( 14) leads to the following quantum Langevin equations (QLEs) by including the dissipation and input noise of each mode, which, in the frame rotating at the drive frequency ω 0 , are given by where ∆ a = ω a − ω 0 , and ∆ m = ω m − ω 0 .a in and m in are input noise operators for the cavity and magnon modes, respectively, which are zero-mean and obey the following correlation functions: ξ is the Langevin force operator, accounting for the Brownian motion of the mechanical oscillator, which is autocorrelated as Here the Markovian approximation has been taken, which is valid for a large mechanical quality factor m, b) are the equilibrium mean thermal photon, magnon, and phonon number, respectively, with k B the Boltzmann constant and T the bath temperature.
The magnon mode is strongly driven to enhance the dispersive magnomechanical coupling strength, leading to a large amplitude of the magnon mode |⟨m⟩| ≫ 1.This allows one to linearize the system dynamics around the large average values.Consequently, the linearized QLEs for the quantum fluctuations can be obtained [24], which can be expressed compactly in the matrix form of u(t) = Au(t) + n(t), where u(t) = [δX a (t), δY a (t), δX m (t), δY m (t), δq(t), δp(t)] T , n Here the quantum fluctuations are written in the quadrature form, i.e., δX O = (δO ), and the noise operators are defined in the same way.The drift matrix A is given by where ∆m = ∆ m + g mb ⟨q⟩ is the effective detuning including the frequency shift due to the magnomechanical interaction, and G mb = i √ 2g mb ⟨m⟩ is the effective magnomechanical coupling strength, with ⟨q⟩ = −(g mb /ω b )|⟨m⟩| 2 and ⟨m⟩ = Ω(i∆ a +κ a ) Because of the linearized dynamics and the Gaussian nature of the quantum noises, the steady state of the quadrature fluctuations is a three-mode Gaussian state, which can be completely characterized by a 6 × 6 covariance matrix (CM) V, with its entries defined as ).The steady-state CM V can be obtained straightforwardly by solving the Lyapunov equation [24,103].With the CM in hand, bipartite and tripartite entanglement of the system can then be quantified by using the logarithmic negativity E N [104,105] and the minimum residual contangle R min τ [106,107], respectively.A precondition for creating entanglement in the system is the low-frequency mechanical mode being significantly cooled (close) to its quantum ground state.To this end, a redtuned microwave field is applied to drive the magnon mode with the detuning ∆m ≃ ω b ≫ κ m (Figure 8(a)).For a relatively weak drive field, this activates the magnomechanical anti-Stokes scattering responsible for cooling the mechanical motion.When the drive field becomes sufficiently strong, the magnomechanical coupling breaks the weak-coupling condition (G mb ≪ ω b ) for taking the rotating wave (RW) approximation to obtain the cooling interaction ∝ δmδb † + δm † δb.Consequently, the counter-RW terms ∝ δm † δb † + δmδb in the linearized magnomechanical interaction start to play the role, i.e., the magnomechanical parametric down-conversion (PDC) interaction is activated, yielding the magnomechanical entanglement.The entanglement is further distributed to the cavity-magnon and cavityphonon subsystems when the cavity resonates with the mechanical Stokes sideband, i.e., ∆ a ≃ −ω b .In this situation, all bipartite subsystems are entangled, as shown in Figure 8(b).
A more fundamental reason lies in the fact that the two mechanical sidebands, i.e., the Stokes and anti-Stokes sidebands, are entangled via the mediation of the mechanical oscillator, which is involved in both the Stokes and anti-Stokes scattering processes.In the Stokes scattering, the mechanics and the Stokes sideband get entangled due to the PDC, whereas in the anti-Stokes scattering, the mechanics and the anti-Stokes sideband realize an effective state-swap interaction.Therefore, the two sidebands become entangled due to the mechanical mediation.Since the cavity and magnon modes are, respectively, resonant with the two sidebands, the cavity and magnon modes thus get entangled.Apart from the simultaneous presence of all bipartite entanglements, the system also shares a genuine tripartite entanglement, as witnessed by a nonzero R min τ in Figure 8(c).The magnon-photon-phonon entanglement can be enhanced by integrating a coherent feedback loop, which contains a controllable beam splitter that determines how much the cavity output field is coherently sent back into the cavity [108].Since the coherent feedback reduces the cavity decay rate [109], yielding an increasement of the magnon amplitude and thus an enhanced magnomechanical coupling, the coherent feedback loop improves the tripartite entanglement.The entanglement can also be enhanced by utilizing the magnon squeezing [110] or the microwave squeezing produced by a parametric amplifier [111,112].In addition, controllable transfer of bipartite entanglements is studied in a parity-timesymmetric-like CMM system by introducing an auxiliary active cavity, which provides a gain to the system [113].Very recently, nonreciprocal entanglement [114,115] and quantum coherence [116] have been studied by adopting a spinning resonator [114,116] or exploiting the chiral cavity-magnon coupling in a torus-shaped microwave cavity [115].

Entanglement between two magnon modes
Two magnon modes of two YIG spheres can get nonlocally entangled by exploiting the nonlinear magnetostriction in one YIG sphere and their linear coupling to a common microwave cavity [58].The system consists of two YIG spheres: one holds a magnon mode and the other supports both a magnon mode and a mechanical vibration mode.This can be realized by adjusting the position of the sphere relative to the direction of the bias magnetic field, since the bare coupling g mb is sensitive to the angle between the direction of the maximum displacement and the bias field [23].The two magnon modes further couple to the same microwave cavity [58], which plays a key role in distributing the entanglement, as discussed later.The Hamiltonian of the system is the sum of equation ( 14) and the following term associated with the second magnon mode (m in equation ( 14) is replaced by m 1 ) where m 2 and m † 2 ([m 2 , m † 2 ] = 1) are the annihilation and creation operators of the second magnon mode, and g m 2 a is the associated cavity-magnon coupling strength.
The mechanism of entangling two magnon modes can be explained by using the partial result of Ref. [24] for the tripartite CMM system, i.e., the cavity and the magnon mode m 1 are entangled under the condition of ∆1 ≃ −∆ a ≃ ω b (here ∆ j ≡ ∆ m j ), cf.Section 4.1.1 and Figure 9(a).Since the second magnon mode m 2 also couples to the cavity via the beam-splitter interaction, which realizes a cavity-magnon state-swap interaction, the two magnon modes thus get entangled because of the mediation of the cavity.The magnon-magnon entanglement is maximized when the cavity and the second magnon mode are resonant, ∆ 2 ≃ ∆ a ≃ −ω b (cf. Figure 9(b)), which is proved to be optimal to realize the cavity-magnon (m 2 ) state transfer [57].To sum up, the optimal condition for getting the magnon-magnon entanglement is thus ∆1 ≃ −∆ a,2 ≃ ω b ≫ κ m [58].
4.1.3.Entanglement between two microwave fields Continuous-variable entanglement of microwave fields can be generated by exploiting the nonlinear magnetostrictive interaction  [58,56,117,118], respectively.and coupling a YIG sphere to two microwave cavities, e.g., of a planar cross-shaped cavity [56].The system consists of two microwave cavities and a YIG sphere supporting a magnon mode and a mechanical mode.The Hamiltonian of the system is the sum of equation ( 14) and the following (a is replaced by a 1 ) where a 2 and a † 2 are the annihilation and creation operators of the second cavity mode, satisfying [a 2 , a † 2 ] = 1, and g ma 2 denotes the coupling rate between the magnon mode with the second cavity mode.
Figure 9(d) shows that two microwave cavity fields get entangled and the maximum entanglement is achieved when the two cavity fields are, respectively, resonant with the two mechanical sidebands, i.e., ∆ 1 = −∆ 2 ≃ ±ω b , where the "±" sign is due to the symmetry of the two cavity fields, and the magnon mode resonant with the blue sideband corresponds to the anti-Stokes process, which significantly cools the phonon mode, thus, removing the main obstacle for observing entanglement.The entanglement originates from the nonlinear magnetostrictive coupling and is distributed to two microwave fields due to the linear magnon-photon coupling.More intuitively, the mechanical motion scatters the microwave driving photons onto two sidebands, which are entangled due to the mechanical mediation as discussed in Section 4.1.1,and if the two cavities are, respectively, resonant with the two sidebands, the two cavity fields get entangled.[117] is proposed to entangle the mechanical vibration modes of two massive YIG spheres in a dual-CMM system, with each subsystem having the same composition as the model described in Section 4.1.1,shown in Figure 9(e).The two cavities are driven by a two-mode squeezed vacuum field, which entangles the two intracavity fields, and, owing to the cavity-magnon beam-splitter interaction, the two magnon modes thus get entangled.Then, each magnon mode is directly driven by a strong red-detuned microwave field, which activates the magnomechanical stateswap interaction, allowing for the transfer of squeezing from the magnon mode to the phonon mode.Consequently, the two phonon modes of two YIG spheres become entangled.However, the above scheme needs an external squeezed vacuum field.By fully exploiting the nonlinear magnetostriction, a more energy-saving protocol [118] is offered without using any external quantum resources.The protocol is based on a hybrid five-mode CMM system, including a microwave cavity mode, two magnon modes, and two mechanical vibration modes, as depicted in Figure 9(f) .The vibration modes and the magnon modes are coupled in a dispersive manner by magnetostriction, and the latter also couple to the microwave cavity via the magnetic-dipole interaction.

Entanglement between two vibrational modes A scheme
The first step is to entangle the vibration mode b 1 and the magnon mode m 2 .From the result of Ref. [24], i.e., the mode b 1 and mode a are entangled, the mode b 1 and mode m 2 are thus entangled due to the cavity-magnon (a-m 2 ) state-swap interaction.Then, switch off the pump on the mode m 1 and, simultaneously, turn on a red-detuned pulsed drive on the mode m 2 .The latter drive activates the magnomechanical beam-splitter interaction, which swaps the magnonic and mechanical states of m 2 and b 2 .Consequently, the previously generated phonon-magnon entanglement is transferred to the two mechanical modes.

Quantum steering
As another form of quantum correlation that lies in between entanglement and Bell nonlocality [119], Einstein-Podolsky-Rosen (EPR) steering characterizes one's capability to control the state of a remote particle via local measurements on his/her particle entangled with the remote one [120], i.e., it essentially embodies quantum nonlocality.In Ref. [24], it is shown that in the CMM system, weak genuine tripartite magnon-photon-phonon entanglement can be achieved in the steady state under a red-detuned (magnon) drive.However, for a red-detuned drive, the steady-state bipartite steerings among magnons, photons and phonons are weak, and therefore genuine tripartite steering cannot be achieved.To enhance the bipartite steerings, one solution is to perform the time-continuous homodyne detection on the output field of the driven microwave cavity [121].The measurement considerably enhances the bipartite steerings, such that robust genuine photon-magnon-phonon tripartite steering can be achieved in the steady state.It is shown that the bipartite and genuine tripartite steering in the bluedetuned regime are stronger than those in the red-detuned regime, and they are optimized when the cavity is nearly resonant with the mechanical sideband (Figure 10(a)-(d)).This is because, in the absence of the measurement, the stability condition of the system severely restricts the strength of the magnomechanical coupling in the blue-detuned regime.Such a restriction, however, is broken by the measurement, leading to stronger bipartite and genuine tripartite entanglement [121].
Distinct from Bell nonlocality and entanglement, EPR steering is essentially asymmetric with respect to two observers.Such an asymmetric and even one-way steering between the cavity mode and the mechanical (magnon) mode in the CMM system can be achieved by the injection of a squeezed (nonclassical) field into the cavity [112].The different damping rates of the cavity, mechanical, and magnon modes lead to asymmetric states with respect to the subsystems, thus giving rise to the asymmetric steerings (Figure 10(e)-(f)).

Squeezed states
Squeezed states represent a kind of nonclassical states with the noise at certain phases below the vacuum fluctuation.Squeezed states of light or microwave fields find a wide range of applications in continuous-variable quantum information [122] and high-precision measurements [123,124,125,126,127].Squeezed states of magnons and phonons of a large-size ferromagnet are macroscopic quantum states, and the magnon squeezing can also be exploited to improve the detection sensitivity of dark-matter axions [128,129].Below we briefly introduce a series of proposals for preparing squeezed states of magnons, phonons, and microwave fields by exploiting different mechanisms.4.3.1.Squeezing transfer from external quantum fields Considering a typical CMM system with the Hamiltonian (14), a broadband squeezed vacuum field is injected to the microwave cavity, which shapes the noise properties of the quantum fluctuations of the cavity field, leading to a squeezed cavity field [57].The squeezing is then transferred to the magnon mode due to the cavity-magnon beam-splitter (state-swap) interaction.A high transfer efficiency requires the resonant condition ∆ a = ∆ m = 0 (in the frame rotating at the squeezed drive frequency ω s , cf. Figure 11(a)-(b)), the strong coupling, and a much smaller magnon decay rate compared to the cavity decay rate, i.e., g ma ≫ κ a ≫ κ m [57,130].
Once the magnon mode is squeezed, by further driving it with a strong reddetuned microwave field with frequency ω 0 (Figure 11(a)), the magnomechanical state-swap interaction is activated, resulting in the transfer of squeezing from magnons to phonons.The state transfer is optimal when the magnon mode is resonant with the anti-Stokes mechanical sideband, ∆ m = ∆ a = ω b (redefined with respect to the strong drive frequency ω 0 , cf. Figure 11(c)).Note that the red-detuned magnon drive also cools the mechanical mode to its quantum ground state, which is a precondition to achieve the squeezing below vacuum noise [57].

Ponderomotive-like squeezing
The dispersive (radiation-pressure-like) magnomechanical interaction can be exploited to generate magnon squeezing [131].Specifically, a deformation displacement of the YIG sample is caused by the magnetostriction, proportional to the magnon excitation number, which in turn modulates the phase of the magnon mode, giving rise to a correlation between the amplitude and the phase of the magnon mode.This correlation results in magnon quadrature squeezing.The mechanism is akin to the ponderomotive squeezing of light induced by radiation pressure in optomechanics [132,133,134,135,136,137], so can be termed as the ponderomotive-like squeezing.The squeezing is then transferred to the cavity mode due to the magnon-cavity beam-splitter coupling, and finally to the output field of the cavity, cf. Figure 11(d), which can be measured via a homodyne detection.
The optimal condition for squeezing requires a strong red-detuned magnon drive field with a small detuning 0 < ∆ m ≡ ω m − ω d < ω b , and a remarkable stationary squeezing appears at the mechanical frequency ω ≃ ω b (in the frame rotating at the drive frequency ω d , Figure 11(e)).Such ponderomotive-like squeezing increases with the magnomechanical cooperativity C m,b = |G mb | 2 /(κ m γ b ), and can be potentially strong [131].

Reservoir-engineered squeezing
The reservoir engineering technique has been successfully applied to the field of optomechanics to prepare squeezed states of mechanical motion [138,139,140,141,142].In view of the similarity between the opto-and magnomechanical interactions, stationary mechanical squeezing can also be obtained by simultaneously driving the magnon mode with both red-and blue-detuned microwave fields [143,144], cf. Figure 12(a).
The two-tone drive can be described by the driving Hamiltonian H dri /ℏ = Ω + e −iω + t + Ω − e −iω − t m † + H.c., where Ω + and Ω − are the Rabi frequencies associated with the two drive fields at frequencies ω ± = ω m ± ω b , respectively.The two-tone drive leads to the steady-state mean value of the magnon mode being approximated as ⟨m⟩ ≈ ⟨m + ⟩e −iω + t + ⟨m − ⟩e −iω − t , where ⟨m ± ⟩ = Ω ± / ω ± − ω m + iκ m − g 2 ma ω ± −ω a +iκ a .The linearized QLEs of the quantum fluctuations, in the interaction picture with respect to ℏω m (a † a+m † m)+ℏω b b † b, are given by where the resonant case ω a = ω m is considered and G ± = g mb ⟨m ± ⟩ are the effective magnomechanical couplings associated with the two drive fields.The squeezing or noise properties can be investigated by solving the dynamical covariance matrix of the system [143].
The system dynamics can be conveniently solved by taking the RW approximation under the condition of κ a(m) , g ma , G ± ≪ ω b , which allows us to neglect the time-dependent coupling terms in equation (28).This then leads to the effective coupling between the magnon and mechanical modes described by δm † G + δb † + G − δb + H.c., which can be rewritten as + , and δB = cosh r δb + sinh r δb † being the Bogoliubov mode, where r = 1 2 ln G − +G + G − −G + is the squeezing parameter.It describes a beam-splitter interaction between the magnon mode and the Bogoliubov mode, and can thus be used to cool the Bogoliubov mode.Since the ground state of the Bogoliubov mode corresponds to the squeezed vacuum state of the mechanical mode, mechanical squeezing can be achieved by cooling the Bogoliubov mode.For the stability reason, G + < G − is required.It is worth noting that for a fixed G − , by increasing G + (via increasing the blue-detuned drive field) the squeezing parameter r grows, but this simultaneously reduces the cooling rate G of the Bogoliubov mode.Consequently, an optimal ratio of G + /G − exists for the mechanical squeezing as a result of the trade-off between the above two effects, cf. Figure 12(b).
Despite the above analysis for the mechanical squeezing, the two-tone drive also induces the magnon squeezing and consequently the intracavity field squeezing, but the degrees of squeezing of both the magnon and intracavity field are very small [143].Instead, Ref. [144] studies the squeezing of the cavity output field and shows that the squeezing can, however, be significantly large (Figure 12(c)).The squeezed output field can be readily accessed and thus can be directly applied to quantum information science and quantum metrology.

Quantum memory
Due to their low damping rates, the mechanical modes of the CMM system can serve as quantum memory for microwave or magnonic quantum states.Further, by exploiting the unique advantages of the CMM system, e.g., tunable magnon frequency and linearized magnomechanical coupling, strong photon-magnon coupling, etc, the CMM system provides a promising platform to efficiently realize the photon-magnon-phonon state transfer.In Ref. [145], the magnon-mediated photon-to-phonon transfer in non-Markovian environments is studied, where the cavity and mechanical modes are assumed to be in structured environments with an Ohmic-and 1/ f -like spectra, respectively.It shows that the fidelity of the photon-to-phonon state transfer can be considerably enhanced due to the presence of the non-Markovian environments.
Differently, Ref. [146] adopts the stimulated Raman adiabatic passage (STIRAP) [147] to achieve high-fidelity state transfer between the cavity and mechanical modes, where the magnon-phonon coupling and the detunings of the cavity and magnon modes with respect to the drive field are designed to be time-dependent.Figure 13(a)-(c) show the instantaneous fidelities for the state transfer and retrieval processes for an initial coherent state |α⟩, Schrödinger's cat state N(|α⟩ + | − α⟩) (N as the normalization factor), and squeezed state e (ra 2 −r * a †2 )/2 |0⟩ (r as the squeezing parameter), respectively, from the cavity to the mechanical mode and back to the cavity mode.For all coherent states, the retrieval fidelity is nearly unit.For cat and squeezed states, even though the retrieval fidelity is rather high, it no longer reaches unit and decreases with the increase of α and r.
In addition, shortcut-to-adiabatic (STA) protocols are provided in Ref. [148] to achieve the high-fidelity state transfer.The STA protocols are based on the counterdiabatic Hamiltonian for transitionless quantum driving (TQD) [149] and the Levis-Riesenfeld (LR) invariant for inverse engineering [150], respectively.It is shown that the TQD approach manifests its power for a larger cat state by achieving a perfect transfer population and the LR-invariant-based approach can further achieve a perfect transfer for any superposed states (Figure 13(d)-(f)).

Optomagnomechanics
It is quite natural to couple the magnetostriction-induced mechanical displacement with an optical cavity via the radiation pressure [49].This then forms a hybrid system consisting of magnons, vibration phonons and optical photons, which is termed as optomagnomechanics (OMM).The OMM system contains both the magno-and optomechanical interactions, which are both a dispersive type (in the case of ω b ≪ ω m ) and thus provide rich nonlinearities for, e.g., preparing various quantum states in the system [37,151,152].A promising physical realization could adopt a micron sized YIG bridge structure [153], which supports a magnon mode with the frequency in gigahertz and a mechanical vibration mode with the frequency ranging from tens to hundreds of megahertz, thus possessing a dominant dispersive magnomechanical coupling (Section 2).To increase the Q factor of the optical cavity, a small high-reflectivity mirror pad could be attached to the surface of the magnomechanical resonator [37] (Figure 14(a)-(d)).The Hamiltonian of the OMM system is given by [37,154].
where c (c † ) is the annihilation (creation) operator of the optical cavity mode with frequency ω c .The bare optomechanical coupling strength g cb is in the range of 10 2 -10 3 Hz for a micronsize mirror [155,156].The bare magnomechanical coupling g mb is about tens of Hz for a micron-size YIG bridge [37] estimated using the theory of Ref. [59].The driving Hamiltonian H dri = iℏϵ(c † e −iω L t − ce iω L t ) + iℏΩ(m † e −iω 0 t − me iω 0 t ), where ϵ denotes the coupling between the cavity and the driving laser at frequency ω L .It is suggested that the OMM configuration can be used to optically measure the magnon population of the ferromagnet [154].Specifically, by resonantly driving the cavity with a weak laser field and under appropriate conditions, a linear relation can be established between the magnon population N m = |⟨m⟩| 2 and the optical phase quadrature ⟨Y c ⟩ via the mediation of the mechanical displacement, i.e., ⟨Y c ⟩ ≃ −4 √ 2 ϵg cb g mb κ 2 c ω b N m [154], with κ c being the cavity decay rate.The excellent linear dependence indicates that the magnon population of the ferromagnet can be read out in the phase quadrature of the optical field, e.g., by homodyning the cavity output field (Figure 14(a)).
The OMM system also enables many quantum protocols, which may find wide applications in quantum information science and technology.For example, by driving the magnon (cavity) mode with a red-detuned microwave (laser) field to activate the magnomechanical (optomechanical) state-swap interaction, the magnon-phonon-photon quantum state transfer can be realized.Consequently, the magnonic state can be read out in the output field of the cavity via the mechanical transduction [37] (Figure 14(b)).Instead, if the red-detuned magnon drive is replaced with a blue-detuned drive to activate the magnomechanical PDC interaction, stationary optomagnonic entanglement can be created [37].Under these conditions, if by further coupling the magnon mode to a microwave cavity via the beam-splitter interaction [14,15] (forming an extended microwave cavity-OMM system, Figure 14(c)), the above generated optomagnonic entanglement [37] can then be distributed to the microwave cavity, yielding stationary microwave-optics entanglement [151].The microwave-optics entanglement finds particularly important applications in building a hybrid quantum network [157].Employing the microwave cavity-OMM system, an entanglement transfer scheme has been suggested to entangle two microwave fields [158].Based on the system used in Ref. [37], if further coupling the optical cavity to an ensemble of two-level atoms (Figure 14(d)), macroscopic entanglement between magnons and atoms and genuine tripartite atom-magnon-phonon entanglement can be generated [152].Recently, an  [37,154,151,152] to realize magnon population detection [154], magnon state readout and optomagnonic entanglement [37], microwave-optics entanglement [151], and magnon-atomic ensemble entanglement [152].(e) Microwave-to-optics transduction protocol using a linear coupling [162].(f) Hybrid system combining magnomechanics with optomechanics to establish coherent couplings among phonons, magnons and photons [163].Figures (a)-(f) are adapted from Refs.[154,37,151,152,162,163], respectively.alternative protocol has been provided for entangling magnons and atoms based on the CMM system by directly coupling the microwave cavity to an atomic ensemble [159].
The above protocols consider the mechanical frequency to be much lower than the magnon frequency.When the size of the magnomechanical system is significantly reduced (considering, e.g., thin YIG films), the mechanical frequency becomes high and close to the magnon frequency [38,160,161].In this regime, the magnomechanical interaction is dominated by the (strong) linear coupling (Section 2).Based on this coupling, a microwaveto-optics conversion protocol is provided [162] exploiting the state-swap interaction in the microwave-magnon, magno-and optomechanical interactions (Figure 14(e)).

Combining magnomechanics with optomechanics via mechanics-mechanics coupling
Another approach to combining magno-and optomechanics is via the direct mechanicsmechanics coupling, implemented by placing the two systems in direct physical contact (Figure 14(f)).This hybrid system can achieve coherent couplings among magnons, vibration phonons and optical photons [163].The Hamiltonian of such a system reads where k = c, m, b 1 , b 2 (k † ) are the annihilation (creation) operators of the optical cavity, magnon, and two mechanical modes, respectively, ω k are the corresponding mode frequencies, and g b denotes the linear coupling between the two mechanical modes of close frequencies.In such a system, when the magnon (optical cavity) mode is driven by a red-detuned microwave (laser) field, the magnomechanical (optomechanical) beam-splitter interaction is activated, which enables the microwave-to-optics conversion via the mechanical mediation.Moreover, the interference effect between the two mechanical motions is also observed by applying a red-(blue-)detuned drive field to the magnon (cavity) mode [163].

Conclusion and outlook
In conclusion, we have reviewed both experimental and theoretical developments in the emerging field of CMM.The experimental observations of the MMIT (MMIA), magnomechanical dynamical backaction, magnon-phonon cross-Kerr nonlinearity, and magnetostrictively induced MFCs are the typical characteristics of the CMM system possessing a nonlinear magnetostrictive interaction.The system has also recently reached the strong-coupling regime and the polariton-mechanics NMS has been observed.As for the theoretical investigations, the progress is much ahead of the experimental studies and many proposals for preparing various quantum states of magnons, phonons, and microwave photons have been offered.These theoretical findings provide beneficial guidance for future experimental development.Despite the above achievements, there is still enormous potential for the development in this field.The quantum-state protocols indicate that the magnomechanical effective coupling and cooperativity should be significantly improved in order to prepare any quantum states of the system, which is likely the most important direction for future experiments.A landmark achievement would be the generation of real macroscopic quantum states in the laboratory.Besides, the CMM system also offers an ideal platform for nonlinear researches by exploiting rich nonlinearities from the magnon-phonon coupling and the YIG.Below, we provide an outlook for future research directions in this field, including non-Hermitian, PT -and anti-PT -symmetry, and exceptional-point-engineered CMM, as well as self-sustaining dynamics and synchronization phenomena by fully exploiting the nonlinearity of the magnetostriction.At last, we introduce some other magnomechanical platforms beyond the frameworks set by Refs.[4,8,9], which, however, are not the focus of this review.

Non-Hermitian and PT -symmetry CMM
Non-Hermitian physics and parity-time (PT ) symmetry have attracted considerable attention over the past decade [164,165,166].Combining non-Hermitian physics with the novel CMM system represents an important research direction and could considerably enrich the research content in this field.
6.1.1.PT -symmetric CMM The PT -symmetric CMM system can be constructed by realizing an effective gain of the microwave cavity mode with the gain rate κ ′ a , e.g., by exploiting the CPA [39,75] (Section 3.4).Due to the intrinsically weak bare magnomechanical coupling g mb , for a not very strong drive, the conditions for realizing the PT -symmetric CMM system are almost identical with those for realizing the PT -symmetric cavity-magnon system, i.e., the cavity and the magnon mode being resonant ω a = ω m and the cavity gain rate and the magnon decay rate being balanced κ ′ a = κ m ≡ κ [75].By diagonalizing the Hamiltonian, one finds that g ma = κ corresponds to the PT phase transition point, i.e., the exceptional point (EP), of the transition from the PT -symmetric phase (g ma > κ) to the broken-PT -symmetric phase (g ma < κ).However, for a sufficiently strong drive field, the above conditions for achieving the PT symmetry are modified due to the appreciable coupling with the mechanical mode, especially when the strong coupling is reached (Section 3.4).Enhanced sideband response [167], controllable output field transmission [168], magnon blockade [169] and magnon chaos [170] have been studied in the PT -symmetric(-like) CMM system.In addition, ultralow-threshold phonon lasing is expected to occur in the PT -symmetric CMM system, enlightened by the results in an analogous optomechanical system [171].
It is worth noting that for the PT -symmetric CMM system, the stability analysis of the system needs to be given sufficient attention, because the introduction of the gain will strongly modify the stability conditions, making the system more easily tend to be unstable.This becomes particularly important when studying some quantum effects, such as mechanical ground-state cooling and entanglement.
As the counterpart of PT symmetry, anti-PT symmetry related phenomena have been demonstrated or predicted in cavity magnonic systems, such as anti-PT -symmetry phase transition [172], enhanced sensing of anharmonicities [173], unconventional singularity [174] and anti-PT -symmetry enhanced microwave-optics interconversion [175], etc.The addition of the mechanical mode by the magnetostrictive effect would make the CMM system a new platform for exploring anti-PT symmetry and topological properties around EPs.
6.1.2.Exceptional-point-engineered CMM Apart from the introduction of the gain into the system as discussed above, the EP can also be achieved in a non-Hermitian CMM system with a dissipative magnon-photon coupling [176].The interaction Hamiltonian of the system includes a dissipative photon-magnon coupling and a coherent magnon-phonon coupling, i.e., H int /ℏ = g ma (a † m + e iΦ am † ) + g mb m † m(b + b † ).(31) The first term describes the novel non-Hermitian coupling caused by the cavity Lenz effect [177,178], and the coupling phase Φ describes the competing coherent and dissipative couplings: Φ = 0 (π) for the photon-magnon coherent (dissipative) coupling.As clearly shown in Figure 15(a)-(b), an EP emerges due to the presence of the dissipative coupling, which can lead to flexible on-off transmission of the probe field (Figure 15(c)).It is also shown that slow-to-fast light switching and enhanced mechanical cooling rate can be achieved by operating the system approaching the EP [176].All of these indicate that the EPassisted CMM system could enable many applications in both classical and quantum regimes.

Nonlinear dynamics: self-sustaining dynamics and synchronization
Most of the studies in the field of CMM are based on an approximate linearized magnomechanical interaction under a strong pump.However, the magnomechanical interaction is intrinsically nonlinear.Including full nonlinearity of the interaction gives rise to rich nonlinear dynamics of the system, such as high-order sidebands (Section 3.5), selfsustaining dynamics [179], and synchronization of mechanical motions [180].In Ref. [179], diverse nonlinear self-sustaining dynamics are predicted, including non-Gaussian phase spreading, amplitude squeezing, and the mixture of multiple limit-cycle states.The nonlinear magnetostriction can also be used to synchronize two mechanical oscillators [180].It shows that by applying a blue-detuned magnon drive, a strong phase correlation can be established between two mechanical oscillators, leading to their synchronization.It also reveals that the strong cavity-magnon linear coupling provides a new degree of freedom to enhance and modulate the synchronization.Apart from the above, many more nonlinear effects remain to be explored, e.g., by exploiting rich nonlinearities of the YIG [26,69].

Other magnomechanical platforms
The preceding sections focus on the CMM system with a dispersive magnetostrictive interaction, typically of a macroscopic ferromagnet (e.g., a sub-mm sized YIG sphere).However, the magnomechanical system can refer to any system that achieves an effective coupling between magnons and phonons, which can be realized in a variety of configurations [8,12,181].
Of particular interest is the coupling between the magnonic mode and the center-of-mass (CoM) motion of a spherical micromagnet, e.g., a YIG sphere [182,183,184,185].The YIG sphere can be trapped in an external harmonic potential (Figure 16(a)), and an oscillating gradient magnetic field B d induces the coupling between the magnon mode and the CoM motion via Gcos(ω d t)(m + m † )(b + b † ) [182], where G is the coupling rate, ω d is the oscillating frequency of the magnetic field, and b denotes the CoM motion.By including also the linear magnetoelastic coupling, strong coupling between the internal acoustic vibration mode and the external low-frequency CoM motion is shown to be possible [182].By further coupling  [182], ground-state cooling of the CoM motion [183], and spin-magnon-phonon strong coupling [184] and giant enhancement of the magnon-spin coupling [185].(d)-(e) Quantum network protocols with magnonic and mechanical nodes connected by optical pulses [186] or microwave photons [187].(f)-(g) SAW-based CMM using an integrated acoustic cavity [192] or a phononic crystal cavity [193].Figures (a)-(g) are adapted from Refs.[182,183,184,186,187,192,193], respectively.the levitated YIG sphere to a microwave cavity (Figure 16(b)), an effective optomechanicallike coupling between the magnon mode and the CoM motion is established, which can be used to cool the CoM motion to its ground state [183].Impressively, this protocol shows the potential to cool the CoM of a sub-centimeter sized sphere to its quantum ground state [183].In addition, by coupling the levitated YIG sphere to a single nitrogen-vacancy spin (Figure 16(c)), strong tripartite spin-magnon-phonon coupling is promising to be reached via applying a parametric (two-phonon) drive to modulate the CoM motion [184].The spinmagnon-phonon coupling can also be realized in the configuration of a micromechanical cantilever holding a YIG sphere [185], where the previous CoM motion of the levitated YIG sphere is replaced with the mechanical motion of the cantilever.By applying a parametric drive to squeeze the mechanical motion, both the spin-and magnon-phonon couplings can be exponentially enhanced, leading to giant enhancement of the effective magnon-spin interaction [185].
The magnon and mechanical modes can also establish a long-range coupling by the aid of optical pulses [186] or a microwave coaxial cable [187,188].As sketched in Figure 16(d), a quantum network protocol consisting of magnonic and mechanical nodes connected by optical pulses is proposed [186].By exploiting both the optomagnonic [9] and optomechanical [49] interactions, it is possible to realize high-fidelity magnon-to-phonon quantum state transfer and entangle two remote magnonic and mechanical nodes [186].In another configuration combining electromechanical and electromagnonic systems (Figure 16(e)), a deterministic scheme is offered to establish entanglement and asymmetric steering between a macroscopic mechanical oscillator and a magnon mode [187].
The phonons can also be the acoustic phonons of the surface acoustic wave (SAW) interacting with ferromagnetic magnons [189,190,191].In such a system, the phonon frequency is in the GHz range and nearly resonant with the ferromagnetic resonance, thus possessing a linear magnon-phonon coupling.In order to enhance the magnon-phonon coupling, an acoustic cavity is designed, formed by either Bragg reflectors [192] (Figure 16(f)) or phononic crystal [193] (Figure 16(g)), which enables the coherent transduction between magnons and acoustic phonons [192] and local manipulation of magnons using phonons strongly confined in the phononic crystal cavity [193].
(a) (Figure 2(b)), the broad Lorentzian resonance dip corresponds to the polariton mode A − (A + ), whose frequency ω − (ω + ) can be tuned by changing the magnon frequency realized by varying the bias magnetic field.When the polariton A − (A + ) is tuned to be resonant with the anti-Stokes (Stokes) mechanical sideband of the drive field, i.e., at ∆ sd = ω b (∆ sd = −ω b ) in Figure 2(a) (Figure 2(b)), the MMIT (MMIA) is most prominent, due to the strong interference between the two resonances.

Figure 2 .
Figure 2. Measured reflection spectra as a function of the probe-drive detuning ∆ sd (with a fixed driving frequency ω d ) for a series of different bias magnetic fields, for (a) a red-detuned drive ω d < ω − ; (b) a blue-detuned drive ω d > ω + .Zoom-in shows detailed spectra of the magnomechanically induced resonances.The figures are adapted from Ref. [23].

Figure 3 .
Figure 3. (a) [(b)] Effective mechanical damping rate γ eff b versus drive power in the mechanical amplification (cooling) process, when the drive is on resonance with the higher (lower)frequency polariton.Inset: Mechanical power spectrum for the green marker at drive power of 0.54 mW (the orange marker at power of 0.68 mW).(c) Mechanical power spectrum versus mechanical frequency ω b and drive-polariton detuning ∆ d .The white dashed line is the fitting curve.Inset shows the drive tone.The figures are adapted from Ref. [25].

Figure 4 .
Figure 4. Bistable mechanical frequency and linewidth under a red (blue)-detuned drive (a)-(b) [(c)-(d)].(a) and (c): The mechanical frequency shift versus the drive power.(b) and (d): The mechanical linewidth variation versus the drive power.The curves are the fitting results.The blue (orange) triangles are the experimental data obtained via forward (backward) sweep of the drive power.The figures are adapted from Ref. [26].

Figure 5 .
Figure 5. (a) Reflection spectra of the higher-frequency polariton.The smallest decay rate of the polariton is achieved at the CPA frequency ω CPA by tuning the magnon frequency.Polariton-mechanics NMS: (b) Reflection spectra of the polariton-mechanics normal modes versus the drive power.Red dots correspond to the frequencies of the dips in the reflection spectra.The inset shows the polariton-mechanics cooperativity versus the frequency of the polariton.(c) Cross sections at different drive powers.The figures are adapted from Ref. [39].

Figure 6 .
Figure 6.(a) Schematic of an OFC with a set of discrete modes that are equally spaced in frequency.(b) Schematic of the generation of MFCs based on the magnetostrictive effect.(c) The magnonic spectrum exhibits a typical frequency-comb structure with the repetition rate f rep = ω b under the drive amplitude over the threshold.(d) Time evolution of the phonon number under a strong drive over the threshold.Inset: phonon number in the logarithmic form versus time.(e) Schematic of the experimental setup for the generation of MFCs via giant mechanical oscillation.(f) The magnonic spectra under the pump power of 20 dBm.Figures (a)-(d) are adapted from Ref. [88] and Figures (e)-(f) are adapted from Ref. [89].

Figure 7 .
Figure 7. (a) Thermometric relation as a function of bath temperature.The solid line represents the analytical result and the solid circles are numerical simulation.(b) Thermometric relation for different values of the microwave cavity frequency (in resonance with the magnon mode): ω m =10, 25, 50 GHz.The figures are adapted from Ref. [93].

Figure 9 .
Figure 9. (a) Frequencies and linewidths of the system adopted to generate magnon-magnon entanglement.(b) Density plot of the magnon-magnon entanglement versus ∆ a and ∆ 2 .(c) Frequencies and linewidths of the system used to entangle two microwave cavity fields.(d) Density plot of the cavity-cavity entanglement vs ∆ 1 and ∆ 2 .(e) Sketch of the dual-CMM system used to entangle two mechanical modes.Two cavities are driven by a two-mode squeezed vacuum field.(f) Sketch of the CMM system used to entangle two mechanical modes of two YIG spheres placed in a common cavity.Figures (a)-(b), (c)-(d), (e) and (f) are adapted from Refs.[58,56,117,118], respectively.

Figure 10 .
Figure 10.(a)-(d) Bipartite steering S i| jk and genuine tripartite steering S cmb versus magnonand cavity-drive detunings ∆ m and δ c via weak continuous measurement.Figures are adapted from Ref. [121].(e) Photon-phonon steering G a⇆b versus G mb ; (f) Photon-magnon steering G a⇆m versus κ a .Figures are adapted from Ref. [112].

Figure 11 .
Figure 11.(a) Frequencies of the system used to generate squeezed states of magnons and phonons.(b) Variance of the magnon amplitude ⟨δx(t) 2 ⟩ vs ∆ m and ∆ a .(c) Variance of the mechanical position ⟨δq(t) 2 ⟩ vs ∆m and ∆ a .(d) Schematic diagram of the ponderomotive-like magnon squeezing, leading to cavity output field squeezing.(e) Noise spectral density of the output field S out W (ω) versus ω and ∆ a .The blank areas in (b), (c) and (e) denote the noise above the vacuum fluctuation.Figures (a)-(c) and (d)-(e) are adapted from Refs.[57, 131], respectively.

Figure 12 .
Figure 12.(a) Sketch of the reservoir-engineered CMM system used to achieve strong squeezing of mechanical motion and cavity output field.The magnon mode is simultaneously driven by two microwave fields at frequencies ω ± .(b) Mechanical squeezing (in units of dB) versus G + /G − at different mean thermal phonon numbers n th .(c) Output field squeezing versus G + /G − at different bath temperatures.Figures are adapted from Refs.[143, 144].

Figure 13 .
Figure 13.Transfer fidelities for a cavity input coherent state (a), cat state (b) and squeezed state (c) via modified STIRAP.Figures are adapted from Ref. [146].The target-state population of the phonon mode for the initial (d) Fock state |n = 4⟩ and (e) cat state of α = 4 in the TQD protocol, and the initial (f) cat states of α = 1, 2, 4 in the LR-invariant approach.Figures are adapted from Ref. [148].

Figure 15 .
Figure 15.Transmission rate of the probe field versus the magnon-cavity detuning δ H = ω m − ω a and the probe-pump detuning ∆ p = ω p − ω 0 − ω b , with ω p (ω 0 ) being the frequency of the probe (pump) field.The figures are adapted from Ref. [176].