Next Article in Journal
Application of Terahertz Time-Domain Spectroscopy to Study the Microheterogeneities of Solutions: A Case Study of Aqueous Sugar Solutions
Next Article in Special Issue
Measurement-Based Control of Quantum Entanglement and Steering in a Distant Magnomechanical System
Previous Article in Journal
Minimisation of Parasitic Capacitance in Lumped-Element Electro-Absorption Modulators for High-Speed Optical Components
Previous Article in Special Issue
Proposal for Optomagnonic Teleportation and Entanglement Swapping
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Generation of Second-Order Sideband through Nonlinear Magnetostrictive Interaction

1
School of Information Engineering, Nanyang Institute of Technology, Nanyang 473004, China
2
School of Physics, Huazhong University of Science and Technology, Wuhan 430074, China
*
Author to whom correspondence should be addressed.
Photonics 2023, 10(8), 886; https://doi.org/10.3390/photonics10080886
Submission received: 13 July 2023 / Revised: 27 July 2023 / Accepted: 30 July 2023 / Published: 1 August 2023
(This article belongs to the Special Issue Hybrid Quantum Magnonics)

Abstract

:
Nonlinear interaction between the magnon mode and the mechanical mode in a magnomechanical system is treated analytically where the magnon mode is coherently driven by a bichromatic microwave drive field consisting of a strong pumping field and a weak probe field and that works within a perturbative regime. Using experimentally achievable parameters, we show that the magnonic second-order sideband is generated and can be considerably enhanced by increasing the power of the pumping field. The suppression of the magnonic second-order sideband generation at the resonance point is discussed. Furthermore, the efficiency of magnonic second-order sideband generation can be well controlled by adjusting the applied bias magnetic field strength, which is a particular feature compared to the optical second-order sideband. In addition to offering insights into the magnomechanical nonlinearity, the present results have the potential to pave the way for exploring practical applications for achieving high-precision measurement in magnonics.

1. Introduction

A spin wave is the collective excitation of magnetization in ordered magnets, and a magnon is the spin wave quanta [1,2,3]. The ferrimagnetic insulator yttrium iron garnet (YIG) with high spin density, very low magnetic damping, and a high Curie temperature provides an alternative platform for quantum information processing [4,5,6,7]. The utilizaton of a magnon as a carrier to transmit and process information is an emerging research field in spintronics [1]. Coherent coupling between magnons and microwave photons through magnetic dipole interaction has been demonstrated experimentally in a three-dimensional microwave cavity coupled with a YIG sphere [6], where the interplay can achieve a strong and even an ultrastrong coupling regime, leading to cavity-magnon polaritons [7]. Many important physical phenomena, such as Rabi oscillation, magnetically induced transparency, the Purcell effect [4], and magnon gradient memory [8], have been reported in cavity magnonics based on coherent coupling. Moreover, due to the powerful compatibility of magnons, optical photons [9,10,11,12,13,14,15] and superconducting qubits [16,17,18] can be coherently coupled to magnons, which offers the prospect of fabricating new hybrid quantum systems.
The magnetostrictive effect describes the phenomenon that, when the magnon mode in the YIG sphere is excited, the varying magnetization induced by the magnon excitation leads to geometrical deformation of the YIG sphere, which can be treated as a mechanical resonator [19]. Currently, the magnomechanical system [20,21,22,23,24], which describes the coherent coupling between magnon modes and mechanical modes (phonons) in the YIG sphere via the magnetostrictive effect, has attracted extensive attention in the magnonic field and many exciting advances have been achieved, These include achieving magnon–photon–phonon entanglement [25,26,27], mechanical bistability [28], nonreciprocal transmission [29,30,31,32], squeezing of the microwave field [33], and magnonic frequency combs [34,35,36]. On the other hand, non-Hermitian cavity magnonics has also aroused considerable interest and has emerged as a new frontier in the magnonic field [37,38,39,40]. Non-Hermitian cavity magnomechanical systems have been proposed to generate magnon chaos [41], enhance sideband responses [42], and manipulate optical transmission [43]. It is well known that nonlinear phenomena can offer insights into understanding of the magnon–phonon interaction and can help promote the development of cavity magnomechanics.
In the present work, we theoretically investigate an important nonlinear effect, the second-order sideband in a magnomechanical system. We obtain the analytical expression for the amplitude of the magnonic second-order sideband by employing a perturbation technique. We show that the efficiency of the second-order sideband generation exhibits a strong dependence on the microwave driving field. Moreover, we show that the second-order sideband control can be achieved by adjusting the external magnetic field. All the relevant parameters are well within current experimental reach. Our work establishes a theoretical framework for generating a magnonic second-order sideband and provides a way to study other nonlinear phenomena induced by the magnetostrictive effect and to explore corresponding applications.

2. Physical Setup and Dynamical Equation

We consider a magnomechanical system, as shown in Figure 1a, where a YIG sphere is placed near a port, which is designed for inputting a microwave field to drive the YIG sphere via a superconducting microwave line with a loop antenna. The YIG sphere is magnetized by a uniform bias static magnetic field (along the z direction, i.e., H = H z ^ ) to induce the magnon mode. Due to the finite geometry, the magnetized YIG sphere confines the magnon modes and acts as a magnonic resonator. As a fundamental magnon mode, the Kittel mode is a uniform magnetization precession. The frequency of the Kittel mode is determined by the external bias magnetic field H , viz. ω m = γ H , where γ / 2 π = 28 GHz / T is the gyromagnetic ratio and H denotes the magnetic field strength [6]. For the bias magnetic field ranging from 0 from 1 T , the resonance frequency of the Kittel magnon mode corresponds to 0 to 28 GHz . The quantum mechanical Hamiltonian of the Kittel mode can be expressed as H ^ m = ω m m ^ m ^ via the Holstein–Primakoff transformation [44], where m ^ = V m 2 γ M s ( M x i M y ) is the annihilation operator of the Kittel mode, with V m the YIG sphere volume, M s is the saturation magnetization, and M x , y are the magnetization components. Here, we assume that the Kittel magnon mode of the YIG sphere is excited by a bichromatic microwave drive field, and setting of the magnetic component of the microwave drive field in the x y plane. The Hamiltonian of the driving of the Kittel mode by the microwave drive fields consisting of the pumping and probe fields is of the form H ^ d = i Ω 0 ( m ^ e i ω 0 t m ^ e i ω 0 t ) + i Ω p ( m ^ e i ω p t m ^ e i ω p t ) within the condition of the low-lying excitations [25], where Ω i = 5 4 γ N B i ( i = 0 , p ) are the Rabi frequencies denoting the coupling strength between the drive magnetic fields and the magnon mode, where N = ρ V m is the total number of spins in the YIG sphere, with ρ = 4.22 × 10 27 m 3 the spin density of the YIG sphere, and B 0 and B p are the amplitudes of the drive magnetic fields, with frequencies ω 0 and ω p corresponding to the pumping field and the probe field, respectively. The power of the microwave drive field P i = B i 2 2 μ 0 S c ( i = 0 , p ) , where μ 0 is the vacuum magnetic permeability, c is the speed of an electromagnetic wave propagating through the vacuum, and S = π R 2 is the cross-sectional area, with R the radius of the YIG sphere; therefore, B i = 1 R 2 P i μ 0 π c .
According to the magnetostrictive effect, the varying magnetization induced by the microwave magnetic field causes deformation of the YIG sphere and the deformation can also change the magnetization, which leads to a nonlinear interaction between the Kittel magnon mode and the vibrational mode. The Hamiltonian of the magnetostrictive interaction can be written as H ^ i n t = g m b m ^ m ^ ( b ^ + b ^ ) [20], where b ^ is the annihilation operator of the phonon (the quantum of the mechanical vibration mode) and g m b is the single magnon–phonon coupling strength, which is simmlar to the optomechanical interaction [45,46]. In sum, the Hamiltonian of the system includes three parts, i.e., H ^ = H ^ m + H ^ d + H ^ i n t , which is the starting point of the following discussion.
Based on the above Hamiltonian and by introducing the dissipation with the Markov approximation [47] in a frame rotating with frequency ω 0 , i.e., a unitary transformation U ^ ( t ) = exp i ω 0 m ^ m ^ t is applied, the Heisenberg–Langevin equations ( o ^ ˙ = i [ H ^ , o ^ ] ) [48] describing the system dynamics are given by:
d m ^ d t = ( i Δ m κ m ) m ^ i g m b m ^ ( b ^ + b ^ ) + Ω 0 + Ω p e i δ t + 2 κ m m ^ i n ,
d b ^ d t = ( i ω b κ b ) b ^ i g m b m ^ m ^ + 2 κ b b ^ i n ,
where Δ m = ω m ω 0 is the detuning between the pumping field and the magnon mode, δ = ω p ω 0 is the detuning between the pumping field and the probe field, ω b is the frequency of the mechanical mode, and κ m and κ b are, respectively, the decay rate of the magnon and the mechanical mode; m i n and b i n are the noise operators for the magnon and mechanical modes, respectively, which are zero mean, i.e., m ^ i n ( t ) = 0 and b ^ i n ( t ) = 0 , and are characterized by the correlation functions m ^ i n ( t ) m ^ i n ( t ) = N m ( ω m ) + 1 δ ( t t ) , m ^ i n ( t ) m ^ i n ( t ) = N m ( ω m ) δ ( t t ) , and b ^ i n ( t ) b ^ i n ( t ) = N b ( ω b ) + 1 δ ( t t ) , b ^ i n ( t ) b ^ i n ( t ) = N b ( ω b ) δ ( t t ) [47], where N m ( b ) ( ω m ( b ) ) = exp ( ω m ( b ) / k B T ) 1 1 are the equilibrium mean thermal magnon (phonon) number, with k B the Boltzmann constant and T the environmental temperature. Within the semiclassical approximation, the magnon and phonon operators are reduced to their expectation values [34,49], viz., m = m ^ , m * = m ^ , b = b ^ , and b * = b ^ , and the quantum noise terms can be safely dropped. In addition, the mean-field approximation by factorizing averages, i.e., m ^ m ^ = m * m , is used to deal with the nonlinear term i g m b m ^ m ^ . Then, the Heisenberg–Langevin equations read
d m d t = ( i Δ m κ m ) m i g m b m ( b + b * ) + Ω 0 + Ω p e i δ t ,
d b d t = ( i ω b κ b ) b i g m b m * m .
Assuming that the pumping field is much stronger than the probe field, we can employ the perturbation method for solving the evolution Equations (2a) and (2b); that is, the solution of the Equations (2a) and (2b) can be written as the sum of the steady-state values and the small perturbations. Physically, the strong pumping field provides the steady-state solutions ( m 0 , b 0 ) of the system, and the weak probe field provides the perturbations around the stable states. Accordingly, the total solutions of Equations (2a) and (2b) can be described by m = m 0 + δ m and b = b 0 + δ b . Setting d m 0 / d t = 0 and d b 0 / d t = 0 , the steady-state solutions of Equations (2a) and (2b) can be obtained as
m 0 = Ω 0 i Δ m + κ m i g m b ( b 0 + b 0 * ) ,
b 0 = i g m b | m 0 | 2 i ω b + κ b .
| m 0 | 2 and | b 0 | 2 denote, respectively, the steady-state average magnon number and the phonon number provided by the pumping field Ω 0 . According to Equations (3a) and (3b), we can obtain the following function relationship of | m 0 | 2 :
4 g m b 4 ω b 2 ( ω b 2 + κ b 2 ) 2 ( | m 0 | 2 ) 3 + 4 Δ m g m b 2 ω b ω b 2 + κ b 2 ( | m 0 | 2 ) 2 + ( κ m 2 + Δ m 2 ) | m 0 | 2 = 5 γ 2 N μ 0 P 0 8 S c .

3. Results and Discussion

We display the steady-state magnon number | m 0 | 2 as a function of the power of the microwave drive field P 0 for different values of the magnomechanical coupling strength g m b in Figure 2. As shown by the red line in Figure 2, without the magnomechanical coupling, i.e., g m b = 0 , the steady-state magnon number | m 0 | 2 undergoes linear growth with increase in the power of the microwave drive field P 0 . Interestingly, the steady-state magnon number | m 0 | 2 experiences a bistable behavior of P 0 in a certain parameter region when we consider the magnomechanical coupling, implying that the magnonic response exhibits two different stable states. As shown by the blue line (corresponding to the magnomechanical coupling strength g m b / 2 π = 9.88 mHz ) in Figure 2, for the case where P 0 < 143.1 mW , only one solution exists and the system has no bistability. For the case where 143.1 mW < P 0 < 8.867 W , the cubic Equation (4) for the steady-state magnon number | m 0 | 2 yields three real roots. The largest and smallest roots are stable, and the middle one is unstable, which is indicated by the dashed line. The system gives rise to bistability in this region [6,50]. The magnonic bistability can be seen from the hysteresis loop. Specifically, the steady-state magnon number | m 0 | 2 initially lies in the lower stable branch (corresponding to the smallest root). When the power of the microwave drive field is gradually increased, the steady-state magnon number | m 0 | 2 scans the lower stable branch. When the steady-state magnon number | m 0 | 2 arrives at the end of the lower stable branch, i.e., the first critical point, it jumps to the upper stable branch. After jumping to the upper stable branch, the steady-state magnon number | m 0 | 2 starts decreasing and still follows the upper stable branch when the power of the microwave drive field is decreased. When the steady-state magnon number | m 0 | 2 reaches to the second critical point, it jumps down to the lower stable branch. Furthermore, when we increase the magnomechanical coupling strength g m b to 1.1 × 9.88 mHz , the power of the microwave drive field needed to observe the bistability is relatively lower, as shown by the black line in Figure 2. Here, the power of the pumping field is kept below 10 mW throughout this work to avoid bistability in the system. Specifically, the low-lying excitations indicate m ^ m ^ 2 N s , where s = 5 / 2 is the spin number of the ground state Fe 3 + ion in the YIG sphere. For a 250- μ m -diameter YIG sphere, the total number of spins N 3.5 × 10 16 [20]. The magnon number for P 0 = 10 mW is about 1.1 × 10 14 , which satisfies the assumption of the low-lying excitations | m | 2 5 N = 17.5 × 10 16 .
By substituting the ansatzs m = m 0 + δ m and b = b 0 + δ b into the Equations (2a) and (2b), we can obtain the evolution equation of the perturbations δ m and δ b made by the probe field:
ϕ = M ϕ + N ϕ * + Ω p e i δ t σ ,
where ϕ = ( δ m , δ b ) T , σ = ( 1 , 0 ) T , and
= i Δ m + κ m + d / d t 0 0 i ω b + κ b + d / d t ,
M = i g m b ( b 0 + b 0 * + δ b + δ b * ) i g m b m 0 i g m b ( m 0 * + δ m * ) 0 ,
N = 0 i g m b m 0 i g m b m 0 0 .
In the perturbative regime, we assume that the perturbations δ m and δ b can be expressed in the following series form [50]:
δ m = m 1 + e i δ t + m 1 e i δ t + m 2 + e 2 i δ t + m 2 e 2 i δ t ,
δ b = b 1 + e i δ t + b 1 e i δ t + b 2 + e 2 i δ t + b 2 e 2 i δ t ,
where the coefficients m 1 ± and m 2 ± are the amplitude of the magnonic first- and second-order sideband components with the frequency ω 0 δ and ω 0 2 δ , respectively. m 1 + ( m 2 + ) and m 1 ( m 2 ) correspond to the lower and upper sidebands, respectively. The physical picture of such ansatzes Equations (6a) and (6b) can be understood as follows: When the Kittel mode in the YIG sphere is excited by the pumping field and the probe field, there are sidebands generated with a series of frequencies ω 0 ± n δ , where the integer n is the order of the sideband. These frequency components are generated by the parametric frequency conversion stemmed from the nonlinear terms i g m b δ m ( δ b + δ b * ) and i g m b δ m * δ m , which is quite similar to second- and high-order harmonics generation in a nonlinear medium [51,52]. The first-order upper sideband ω 0 + δ is referred to as the anti-Stokes field, and the first-order lower sideband ω 0 δ is known as the Stokes field. The output field with a frequency ω 0 + 2 δ is the second-order upper sideband, and the frequency ω 0 2 δ is the second-order lower sideband. The third-order sideband and the other higher-order sideband components are very weak in contrast to the second-order sideband and are validly ignored in the perturbative regime. By substituting the ansatzes Equations (6a) and (6b) into the Equation (5), we can obtain
m 1 = Ω p Ξ ( δ ) i g m b λ m 0 ,
m 2 = i g m b m 0 ( Θ + Υ ) ϖ ( ψ , θ ) ρ ( δ )
where Ξ ( δ ) = i Δ m + κ m i δ i g m b ( b 0 * + b 0 ) , λ = i g m b m 0 * [ 1 / ( i ω b i δ + κ b ) 1 / ( i ω b i δ + κ b ) ] , m 1 + = i g m b m 0 λ * ( m 1 ) * / Ξ ( δ ) , ϖ ( x , y ) = x i g m b ( b 0 + b 0 * ) i g m b m 0 y , θ = π m 0 * + ( π + m 0 ) * , π ± = i g / ( i ω b + κ b ± 2 i δ ) , ψ ± = i Δ m + κ b ± 2 i δ , τ 1 = ( π * + π + ) m 0 * , τ 2 = ( π + π + ) m 0 , Υ = ( π + π + * ) ( m 1 + ) * m 1 + i g m b λ ( m 1 ) 2 , and
ρ ( δ ) = g m b 2 | m 0 | 2 τ 2 ( π * m 0 + π + m 0 ) ϖ ( ψ + , τ 1 ) ,
Θ = i τ 2 g m b m 0 * Υ i g m b λ ( m 1 + ) * m 1 ϖ * ( ψ + , τ 1 ) .
Using the input–output relation [48] s o u t = ( Ω 0 + Ω p e i δ t ) / κ m κ m m , we can obtain the output fields of the magnons:
s o u t = Ω 0 / κ m κ m m 0 + ( Ω p / κ m κ m m 1 ) e i δ t κ m m 1 + e i δ t κ m m 2 e 2 i δ t κ m m 2 + e 2 i δ t ,
where the coefficients Ω 0 / κ m κ m m 0 , Ω p / κ m κ m m 1 , and κ m m 1 + denote the amplitude of the output fields at the frequencies ω 0 , ω p , and 2 ω 0 ω p , respectively, and the amplitude of the output fields at the second-order lower sideband with the frequency 3 ω 0 2 ω p and the second-order upper sideband with the frequency 2 ω p ω 0 can be obtained as κ m m 2 + and κ m m 2 , respectively. Here, we focus on the generation of the second-order upper sideband. To reveal substantially the second-order upper sideband generation process, we define the transmission of the probe field as t p = κ m ( Ω p / κ m κ m m 1 ) / Ω p = 1 κ m m 1 / Ω p , and the dimensionless quantity η = | κ m m 2 / Ω p | as the efficiency of the second-order upper sideband generation. Specifically, we scale the output intensity of the second-order upper sideband generation with respect to the intensity of the probe field.
In order to better understand the physical process of the second-order sideband induced by the magnomechanical nonlinearity, in Figure 3a, we plot the transmission ratio of the probe field | t p | 2 as a function of the detuning δ and the power of the pumping field P 0 . As shown in Figure 3b, a transmission peak appears at δ / ω b = 1 , and we find that the transmission characteristics of the probe field exhibit a strong dependence on the pumping field. In detail, when we increase the power of the pumping field to 0.5 mW and 5 mW , the peak values reach about 0.8 and 1, as shown in Figure 3c,d. This implies that the generation efficiency of the first-order upper sideband becomes higher with increase in the power of the pumping field P 0 . The first-order upper sideband is derived from the one-phonon upconversion of the pumping field, which is called an anti-Stokes process [46,53,54]. In addition, Figure 3a shows that the transmission peak of the probe field becomes wider and wider with increase in the power of the pumping field. The width of the transmission peak of the probe field is proportional to the steady-state magnon number, i.e., Γ m = κ m + 4 g m b 2 | m 0 | 2 / κ m [45], and the steady-state magnon number | m 0 | 2 can be enhanced by increasing the power of the pumping field.
Figure 4 shows that the efficiency η of the magnonic second-order sideband generation varies with the detuning δ for different magnomechanical coupling strengths g m b . We can see that the value of η is zero in the case of g m b = 0 , which means that no magnonic second-order sideband is generated due to the nonlinear terms [i.e., i g m b δ m ( δ b + δ b * ) and i g m b δ m * δ m ] are absent. When the magnomechanical interaction is considered, i.e., g m b / 2 π = 9.88 mHz , one can find that an obvious magnonic second-order sideband generation can be observed, i.e., η = 1.3 × 10 3 , which means that the amplitude of the output magnonic second-order sideband is 1.3 × 10 3 of the amplitude of the probe field. The linewidth of the second-order sideband is about 4.2 × 10 3 Hz , arising from the intrinsic mechanical damping rate and the microwave drive field induced damping rate, which can be obtained as κ m + 4 g m b 2 | m 0 | 2 / κ m [45] and is consistent with the first-order sideband spectrum. The increase in the magnomechanical coupling strength g m b enhances the generation efficiency of the magnonic second-order sideband. By experimentation, the improvement in the magnomechanical coupling strength can be realized by scaling down the YIG sphere size [20].
Figure 5a plots the efficiency η of the magnonic second-order sideband generation varying with the detuning δ and the power of the pumping field. We can see that the increased pumping field intensity can result in considerable enhancement of the efficiency η . In detail, as shown in Figure 5b–d, the maximum values of the efficiency η are about 4 × 10 3 , 6 × 10 3 , and 6.4 × 10 3 , which correspond to the power of the pumping field P 0 = 0.1 mW , 1 mW , 5 mW , respectively. On the other hand, one can find that the efficiency η becomes very small at the resonance point δ / ω b = 1 . The physical reason for this phenomenon can be understood as follows: The fundamental mechanism for second-order sideband generation is that the anti-Stokes field and the probe field undergo a one-phonon upconverted process. However, destructive interference occurs between the probe field and the anti-Stokes field, which leads to a transparency window at δ / ω b = 1 , as shown in Figure 3a. Then the generation process of the second-order sideband is subdued at δ / ω b = 1 [50]. The linewidth of the second-order sideband becomes wider with increase in the power of the pumping field P 0 , while it is much less than the linewith of the magnon resonance mode (∼ κ m ), such that the magnonic second-order sideband may provide a tool to achieve high-precision measurement.
Finally, we discuss the influence of the applied bias magnetic field on the magnonic second-order sideband generation. Figure 6 shows that the efficiency η can be tuned by adjusting the bias magnetic field strength H. For a certain pumping power, we find that the efficiency of second-order sideband generation can be enhanced by tuning the bias magnetic field. The maximum value of the efficiency η is about 15.4 × 10 3 when the bias magnetic field strength H 280.67 mT , while the efficiency η is about 6 × 10 3 when P 0 = 1 mW , as shown in Figure 5c. This implies that the nonlinear magnetostrictive interaction between the magnon mode and the mechanical mode in the YIG sphere is extremely enhanced. It can be seen that the sideband efficiency is asymmetric, and the critical value H = 280.71 mT corresponds to the frequency of the magnon mode ω m / 2 π = 7.86 GHz . The generation of the magnonic second-order sideband can be well controlled by only tuning the applied bias magnetic field without changing the other system parameters, which provides an effective and feasible way to manipulate the magnetostrictive effect and to realize convenient magnomechanical control.

4. Conclusions

In conclusion, we have theoretically studied the nonlinear interaction between the magnon mode and the mechanical mode and discussed the generation and control of the magnonic second-order sideband in a magnomechanical system by using the perturbation method. We have shown that the magnomechanical nonlinearity makes the system bistable, and that the magnonic second-order sideband can be tuned by the microwave drive field and the external magnetic field. Based on the low damping and the flexible controllability of the system, the magnonic second-order sideband presented here possesses the features of narrow linewidth and convenient control. Beyond fundamental scientific significance, our results may have potential applications in high-precision measurement, such as force sensors [55], magnetostrictive transducers [56] and actuators [57].

Author Contributions

L.Y. carried out the calculations, wrote the main manuscript text and prepared all the figures; B.W. reviewed the manuscript and contributed to the interpretation of the work and the writing of the manuscript; H.X. participated in the discussions. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number: 12004202.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chumak, A.V.; Vasyuchka, V.I.; Serga, A.A.; Hillebrands, B. Magnon spintronics. Nat. Phys. 2015, 11, 453–461. [Google Scholar] [CrossRef]
  2. Soykal, O.O.; Flatte, M.E. Strong field interactions between a nanomagnet and a photonic cavity. Phys. Rev. Lett. 2010, 104, 077202. [Google Scholar] [CrossRef] [Green Version]
  3. Tabuchi, Y.; Ishino, S.; Ishikawa, T.; Yamazaki, R.; Usami, K.; Nakamura, Y. Hybridizing ferromagnetic magnons and microwave photons in the quantum limit. Phys. Rev. Lett. 2014, 113, 083603. [Google Scholar] [CrossRef] [Green Version]
  4. Zhang, X.; Zou, C.-L.; Jiang, L.; Tang, H.X. Strongly coupled magnons and cavity microwave photons. Phys. Rev. Lett. 2014, 113, 156401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Bai, L.; Harder, M.; Chen, Y.P.; Fan, X.; Xiao, J.Q.; Hu, C.-M. Spin pumping in electrodynamically coupled magnon-photon systems. Phys. Rev. Lett. 2015, 114, 227201. [Google Scholar] [CrossRef] [Green Version]
  6. Wang, Y.-P.; Zhang, G.-Q.; Zhang, D.; Li, T.-F.; Hu, C.-M.; You, J.Q. Bistability of cavity magnon-polaritons. Phys. Rev. Lett. 2018, 120, 057202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Zhang, D.; Wang, X.-M.; Li, T.-F.; Luo, X.-Q.; Wu, W.; Nori, F.; You, J.Q. Cavity quantum electrodynamics with ferromagnetic magnons in a small yttrium-iron-garnet sphere. NPJ Quantum Inf. 2015, 1, 15014. [Google Scholar] [CrossRef] [Green Version]
  8. Zhang, X.; Zou, C.-L.; Zhu, N.; Marquardt, F.; Jiang, L.; Tang, H.X. Magnon dark modes and gradient memory. Nat. Commun. 2015, 6, 8914. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Osada, A.; Hisatomi, R.; Noguchi, A.; Tabuchi, Y.; Yamazaki, R.; Usami, K.; Sadgrove, M.; Yalla, R.; Nomura, M.; Nakamura, Y. Cavity Optomagnonics with Spin-Orbit Coupled Photons. Phys. Rev. Lett. 2016, 116, 223601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Kusminskiy, S.V.; Tang, H.X.; Marquardt, F. Coupled spin-light dynamics in cavity optomagnonics. Phys. Rev. A 2016, 94, 033821. [Google Scholar] [CrossRef] [Green Version]
  11. Zhang, X.; Zhu, N.; Zou, C.-L.; Tang, H.X. Optomagnonic Whispering Gallery Microresonators. Phys. Rev. Lett. 2016, 117, 123605. [Google Scholar] [CrossRef] [Green Version]
  12. Sharma, S.; Rameshti, B.Z.; Blanter, Y.M.; Bauer, G.E.W. Optimal mode matching in cavity optomagnonics. Phys. Rev. B 2019, 99, 214423. [Google Scholar] [CrossRef] [Green Version]
  13. Xu, Y.-J.; Song, J. Nonreciprocal magnon laser. Opt. Lett. 2021, 46, 5276–5279. [Google Scholar] [CrossRef]
  14. Liu, Z.-X.; Li, Y.-Q. Optomagnonic frequency combs. Photon. Res. 2022, 10, 2786–2793. [Google Scholar] [CrossRef]
  15. Fan, Z.; Zuo, X.; Qian, H.; Li, J. Proposal for optomagnonic teleportation and entanglement swapping. Photonics 2023, 10, 739. [Google Scholar] [CrossRef]
  16. Tabuchi, Y.; Ishino, S.; Noguchi, A.; Ishikawa, T.; Yamazaki, R.; Usami, K.; Nakamura, Y. Coherent coupling between a ferromagnetic magnon and a superconducting qubit. Science 2015, 349, 405. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Lachance-Quirion, D.; Tabuchi, Y.; Ishino, S.; Noguchi, A.; Ishikawa, T.; Yamazaki, R.; Nakamura, Y. Resolving quanta of collective spin excitations in a millimeter-sized ferromagnet. Sci. Adv. 2017, 3, e1603150. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Liu, Z.-X.; Xiong, H.; Wu, Y. Magnon blockade in a hybrid ferromagnet-superconductor quantum system. Phys. Rev. B 2019, 100, 134421. [Google Scholar] [CrossRef] [Green Version]
  19. Spencer, E.G.; LeCraw, R.C. Magnetoacoustic resonance in yttrium iron garnet. Phys. Rev. Lett. 1958, 1, 241–243. [Google Scholar] [CrossRef]
  20. Zhang, X.; Zou, C.-L.; Jiang, L.; Tang, H.X. Cavity magnomechanics. Sci. Adv. 2016, 2, 1501286. [Google Scholar] [CrossRef] [Green Version]
  21. Li, J.; Zhu, S.-Y. Entangling two magnon modes via magnetostrictive interaction. New J. Phys. 2019, 21, 085001. [Google Scholar] [CrossRef]
  22. Xu, G.-T.; Zhang, M.; Wang, Z.-Y.; Wang, Y.; Liu, Y.-X.; Shen, Z.; Guo, G.-C.; Dong, C.-H. Ringing spectroscopy in the magnomechanical system. Fundam. Res. 2023, 3, 45–49. [Google Scholar] [CrossRef]
  23. Li, J.; Gröblacher, S. Entangling the vibrational modes of two massive ferromagnetic spheres using cavity magnomechanics. Quantum Sci. Technol. 2021, 6, 024005. [Google Scholar] [CrossRef]
  24. Potts, C.A.; Varga, E.; Bittencourt, V.A.S.V.; Kusminskiy, S.V.; Davis, J.P. Dynamical backaction magnomechanics. Phys. Rev. X 2021, 11, 031053. [Google Scholar] [CrossRef]
  25. Li, J.; Zhu, S.-Y.; Agarwal, G.S. Magnon-photon-phonon entanglement in cavity magnomechanics. Phys. Rev. Lett. 2018, 121, 203601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Yu, M.; Shen, H.; Li, J. Magnetostrictively induced stationary entanglement between two microwave fields. Phys. Rev. Lett. 2020, 124, 213604. [Google Scholar] [CrossRef]
  27. Fan, Z.-Y.; Qiu, L.; Gröblacher, S.; Li, J. Microwave-optics entanglement via cavity optomagnomechanics. arXiv 2022, arXiv:2208.10703. [Google Scholar]
  28. Shen, R.-C.; Li, J.; Fan, Z.-Y.; Wang, Y.-P.; You, J.Q. Mechanical bistability in Kerr-modified cavity magnomechanics. Phys. Rev. Lett. 2022, 129, 123601. [Google Scholar] [CrossRef]
  29. Xu, Y.; Liu, J.-Y.; Liu, W.; Xiao, Y.-F. Nonreciprocal phonon laser in a spinning microwave magnomechanical system. Phys. Rev. A 2021, 103, 053501. [Google Scholar] [CrossRef]
  30. Yang, Z.-B.; Liu, J.-S.; Zhu, A.-D.; Liu, H.-Y.; Yang, R.-C. Nonreciprocal transmission and nonreciprocal entanglement in a spinning microwave magnomechanical system. Ann. Phys. 2020, 532, 2000196. [Google Scholar] [CrossRef]
  31. Kong, C.; Liu, J.; Xiong, H. Nonreciprocal microwave transmission under the joint mechanism of phase modulation and magnon Kerr nonlinearity effect. Front. Phys. 2023, 18, 12501. [Google Scholar] [CrossRef]
  32. Wang, X.; Huang, K.-W.; Xiong, H. Nonreciprocal sideband responses in a spinning microwave magnomechanical system. Opt. Express 2023, 31, 5492–5506. [Google Scholar] [CrossRef]
  33. Li, J.; Wang, Y.-P.; You, J.Q.; Zhu, S.-Y. Squeezing microwaves by magnetostriction. Natl. Sci. Rev. 2023, 10, nwac247. [Google Scholar] [CrossRef]
  34. Xiong, H. Magnonic frequency combs based on the resonantly enhanced magnetostrictive effect. Fundam. Res. 2023, 3, 8–14. [Google Scholar] [CrossRef]
  35. Liu, Z.-X.; Peng, J.; Xiong, H. Generation of magnonic frequency combs via a two-tone microwave drive. Phys. Rev. A 2023, 107, 053708. [Google Scholar] [CrossRef]
  36. Xu, G.-T.; Zhang, M.; Wang, Y.; Shen, Z.; Guo, G.-C.; Dong, C.-H. Magnonic frequency comb in the magnomechanical resonator. arXiv 2023, arXiv:2306.07985. [Google Scholar]
  37. Zhao, J.; Liu, Y.; Wu, L.; Duan, C.-K.; Liu, Y.-X.; Du, J. Observation of anti-PT-symmetry phase transition in the magnon-cavity-magnon coupled system. Phys. Rev. Appl. 2020, 13, 014053. [Google Scholar] [CrossRef] [Green Version]
  38. Wang, B.; Liu, Z.-X.; Kong, C.; Xiong, H.; Wu, Y. Magnon-induced transparency and amplification in PT-symmetric cavity-magnon system. Opt. Express 2018, 26, 20248–20257. [Google Scholar] [CrossRef] [PubMed]
  39. Wang, B.; Jia, X.; Lu, X.-H.; Xiong, H. PT-symmetric magnon laser in cavity optomagnonics. Phys. Rev. A 2022, 105, 053705. [Google Scholar] [CrossRef]
  40. Zhang, G.-Q.; Chen, Z.; Xu, D.; Shammah, N.; Liao, M.; Li, T.-F.; Tong, L.; Zhu, S.-Y.; Nori, F.; You, J.Q. Exceptional point and cross-relaxation effect in a hybrid quantum system. PRX Quantum 2021, 2, 020307. [Google Scholar] [CrossRef]
  41. Wang, M.; Zhang, D.; Li, X.H.; Wu, Y.Y.; Sun, Z.Y. Magnon chaos in PT-symmetric cavity magnomechanics. IEEE Photon. J. 2019, 11, 5300108. [Google Scholar] [CrossRef]
  42. Huai, S.-N.; Liu, Y.-L.; Zhang, J.; Yang, L.; Liu, Y.-X. Enhanced sideband responses in a PT-symmetric-like cavity magnomechanical system. Phys. Rev. A 2019, 99, 043803. [Google Scholar] [CrossRef] [Green Version]
  43. Lu, T.-X.; Zhang, H.; Zhang, Q.; Jing, H. Exceptional-point-engineered cavity magnomechanics. Phys. Rev. A 2021, 103, 063708. [Google Scholar] [CrossRef]
  44. Holstein, T.; Primakoff, H. Field dependence of the intrinsic domain magnetization of a ferromagnet. Phys. Rev. 1940, 58, 1098. [Google Scholar] [CrossRef]
  45. Aspelmeyer, M.; Kippenberg, T.J.; Marquardt, F. Cavity optomechanics. Rev. Mod. Phys. 2014, 86, 1391. [Google Scholar] [CrossRef] [Green Version]
  46. Xiong, H.; Wu, Y. Fundamentals and applications of optomechanically induced transparency. Appl. Phys. Rev. 2018, 5, 031305. [Google Scholar] [CrossRef]
  47. Gardiner, C.W.; Zoller, P. Quantum Noise; Springer: Berlin, Germany, 2000. [Google Scholar]
  48. Walls, D.F.; Milburn, G.J. Quantum Optics; Springer: Berlin, Germany, 1994. [Google Scholar]
  49. Wang, B.; Lu, X.-H.; Jia, X.; Xiong, H. Coherent stimulated amplification of the skyrmion breathing. Chaos Solitons Fractals 2023, 171, 113484. [Google Scholar] [CrossRef]
  50. Xiong, H.; Si, L.-G.; Zheng, A.-S.; Yang, X.; Wu, Y. Higher-order sidebands in optomechanically induced transparency. Phys. Rev. A 2012, 86, 013815. [Google Scholar] [CrossRef] [Green Version]
  51. Ruzicka, B.A.; Werake, L.K.; Xu, G.; Khurgin, J.B.; Sherman, E.Y.; Wu, J.Z.; Zhao, H. Second-harmonic generation induced by electric currents in GaAs. Phys. Rev. Lett. 2012, 108, 077403. [Google Scholar] [CrossRef] [Green Version]
  52. Velotta, R.; Hay, N.; Mason, M.B.; Castillejo, M.; Marangos, J.P. High-order harmonic generation in laser-aligned molecules. Phys. Rev. A 2002, 65, 053805. [Google Scholar]
  53. Chen, H.-J. Multiple-Fano-resonance-induced fast and slow light in the hybrid nanomechanical-resonator system. Phys. Rev. A 2021, 104, 013708. [Google Scholar] [CrossRef]
  54. Chen, H.-J. The fast–slow light transitions induced by Fano resonance in multiple nanomechanical resonators. Opt. Laser Technol. 2023, 161, 109242. [Google Scholar] [CrossRef]
  55. Yang, Q.; Yan, R.; Fan, C.; Chen, H.; Liu, F.; Liu, S. A magneto-mechanical strongly coupled model for giant magnetostrictive force sensor. IEEE Trans. Magn. 2007, 43, 1437–1440. [Google Scholar] [CrossRef]
  56. Kim, Y.Y.; Kwon, Y.E. Review of magnetostrictive patch transducers and applications in ultrasonic nondestructive testing of waveguides. Ultrasonics 2015, 62, 3–19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Chen, L.; Zhu, Y.; Ling, J.; Feng, Z. Theoretical modeling and experimental evaluation of a magnetostrictive actuator with radial-nested stacked configuration. Nonlinear Dyn. 2022, 109, 1277–1293. [Google Scholar] [CrossRef]
Figure 1. (Color online) (a) Schematic diagram of a magnomechanical system. A uniform bias static magnetic field (H) is applied along the z-direction to magnetize the YIG sphere. The magnon mode is directly driven via a microwave antenna. The varying magnetization causes the deformation of the YIG sphere. (b) Frequency spectrum of the output fields. The yellow line indicates the frequency of the pumping field ω 0 , the red line is the frequency of the Stokes field ω 0 δ , and the blue line is the frequency of the anti-Stokes field (the probe field) ω 0 + δ . The frequency of the probe field is ω p , and the detuning between the pumping field and the probe field is δ = ω p ω 0 . The violet line is the second-order sideband ω 0 + 2 δ . The frequency of the magnon mode is ω m .
Figure 1. (Color online) (a) Schematic diagram of a magnomechanical system. A uniform bias static magnetic field (H) is applied along the z-direction to magnetize the YIG sphere. The magnon mode is directly driven via a microwave antenna. The varying magnetization causes the deformation of the YIG sphere. (b) Frequency spectrum of the output fields. The yellow line indicates the frequency of the pumping field ω 0 , the red line is the frequency of the Stokes field ω 0 δ , and the blue line is the frequency of the anti-Stokes field (the probe field) ω 0 + δ . The frequency of the probe field is ω p , and the detuning between the pumping field and the probe field is δ = ω p ω 0 . The violet line is the second-order sideband ω 0 + 2 δ . The frequency of the magnon mode is ω m .
Photonics 10 00886 g001
Figure 2. (Color online) Magnon bistability as seen through the steady-state magnon number | m 0 | 2 varied with the power of the microwave drive field P 0 for the three different magnomechanical coupling strengths g m b / 2 π = 0 , 9.88 mHz , and 1.1 × 9.88 mHz , corresponding to the red, blue, and black lines, respectively. The other system parameters we use are ω b / 2 π = 11.42 MHz , κ m / 2 π = 0.56 MHz , κ b / 2 π = 150 Hz , and Δ m = ω b , respectively, which are chosen from the recent experiment [20].
Figure 2. (Color online) Magnon bistability as seen through the steady-state magnon number | m 0 | 2 varied with the power of the microwave drive field P 0 for the three different magnomechanical coupling strengths g m b / 2 π = 0 , 9.88 mHz , and 1.1 × 9.88 mHz , corresponding to the red, blue, and black lines, respectively. The other system parameters we use are ω b / 2 π = 11.42 MHz , κ m / 2 π = 0.56 MHz , κ b / 2 π = 150 Hz , and Δ m = ω b , respectively, which are chosen from the recent experiment [20].
Photonics 10 00886 g002
Figure 3. (Color online) (a) The transmission spectrum | t p | 2 of the probe field as a function of δ and the power of the pumping field P 0 . (bd) plot the transmission spectrum | t p | 2 in the case of P 0 = 10 μ W , 0.5 mW , and 5 mW , respectively. The magnomechanical coupling strength g m b / 2 π = 9.88 mHz . We use Ω p / Ω 0 = 0.05 , and the other parameters are the same as in Figure 2.
Figure 3. (Color online) (a) The transmission spectrum | t p | 2 of the probe field as a function of δ and the power of the pumping field P 0 . (bd) plot the transmission spectrum | t p | 2 in the case of P 0 = 10 μ W , 0.5 mW , and 5 mW , respectively. The magnomechanical coupling strength g m b / 2 π = 9.88 mHz . We use Ω p / Ω 0 = 0.05 , and the other parameters are the same as in Figure 2.
Photonics 10 00886 g003
Figure 4. (Color online) The efficiency η of the magnonic second-order sideband generation as a function of the detuning δ for different values of the magnomechanical coupling strength g m b . The power of the pumping field is P 0 = 10 μ W . The other parameters are the same as in Figure 2.
Figure 4. (Color online) The efficiency η of the magnonic second-order sideband generation as a function of the detuning δ for different values of the magnomechanical coupling strength g m b . The power of the pumping field is P 0 = 10 μ W . The other parameters are the same as in Figure 2.
Photonics 10 00886 g004
Figure 5. (Color online) (a) The contour map of the efficiency η as a function of the detuning δ and the power of the pumping field P 0 . (bd) plot the spectrum of the magnonic second-order sideband in the case of P 0 = 0.1 mW , 1 mW , and 5 mW , respectively. The magnomechanical coupling strength g m b / 2 π = 9.88 mHz . The other parameters are the same as Figure 2.
Figure 5. (Color online) (a) The contour map of the efficiency η as a function of the detuning δ and the power of the pumping field P 0 . (bd) plot the spectrum of the magnonic second-order sideband in the case of P 0 = 0.1 mW , 1 mW , and 5 mW , respectively. The magnomechanical coupling strength g m b / 2 π = 9.88 mHz . The other parameters are the same as Figure 2.
Photonics 10 00886 g005
Figure 6. (Color online) The efficiency η as a function of the detuning δ and the bias magnetic field strength H. The power of the pumping field is P 0 = 1 mW . The other parameters are the same as Figure 5.
Figure 6. (Color online) The efficiency η as a function of the detuning δ and the bias magnetic field strength H. The power of the pumping field is P 0 = 1 mW . The other parameters are the same as Figure 5.
Photonics 10 00886 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yang, L.; Wang, B.; Xiong, H. Generation of Second-Order Sideband through Nonlinear Magnetostrictive Interaction. Photonics 2023, 10, 886. https://doi.org/10.3390/photonics10080886

AMA Style

Yang L, Wang B, Xiong H. Generation of Second-Order Sideband through Nonlinear Magnetostrictive Interaction. Photonics. 2023; 10(8):886. https://doi.org/10.3390/photonics10080886

Chicago/Turabian Style

Yang, Lei, Bao Wang, and Hao Xiong. 2023. "Generation of Second-Order Sideband through Nonlinear Magnetostrictive Interaction" Photonics 10, no. 8: 886. https://doi.org/10.3390/photonics10080886

APA Style

Yang, L., Wang, B., & Xiong, H. (2023). Generation of Second-Order Sideband through Nonlinear Magnetostrictive Interaction. Photonics, 10(8), 886. https://doi.org/10.3390/photonics10080886

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop