Kerr-Nonlinearity-Triggered Nonclassicality of Magnons in a Photon-Magnon Coupling System

: Hybrid quantum systems have attracted much attention due to the fact that they combine the advantages of different physical subsystems. Cavity QED (cavity quantum electrodynamics) with magnons is a hybrid quantum systems that combines a YIG (Yttrium Iron Garnet) sphere and a 3D (three-dimensional) rectangular microwave cavity. Based on this hybrid photon-magnon system, we obtain an approximate analytic solution by the RWA (rotating wave approximation) with an ingenious transformation. After skillfully diagonalizing the Hamiltonian, we show that the Kerr-nonlinearity interactions could yield a negativity value of the Wigner function, periodic quadrature squeezing effects, antibunching property, and ﬁeld nonclassicality in the magnon. Our work may stimulate the study of nonclassicality of photon-magnon coupling systems and its potential applications in quantum information processing.


Introduction
As we all know, scientists can manipulate various systems at the single quantum level. However, no single type of system can meet all the requirements of quantum information technology. For example, superconducting is recognized as the most competitive solution for scalable solid-state quantum computing, and it is also known to be the nonlinear device with the best performance in the microwave frequency band and ultralow temperature region, but its decoherence time is currently lower than that of natural atoms or ions by three orders of magnitude. Therefore, researchers naturally conceived of utilizing the advantages of solid-state systems, such as microwave resonators, magnon qubits, nanomechanical oscillators and superconducting qubits, to build solid-state hybrid quantum systems with different functions to realize the optimized quantum information processing scheme, which can be explored as a new physical realization method for new quantum computing and meet the development needs of quantum information technology. The purpose of mixing and hybridization is to naturally take advantage of their own strengths and to perform their own duties and functions to achieve more complex and complete information processing, information storage and information transmission functions. In recent years, hybrid quantum systems have attracted much attention due to the fact that they combine the advantages of different physical subsystems [1,2]. Over the past few years, scientists have theoretically proposed [3,4] a hybrid system based on the formation of spin ensemblecoupled microwave cavity modes in single-crystal YIG samples, and obtained experimental verification [5,6]. Unlike spin ensembles in paramagnetic materials with a lower density (such as nitrogen-vacancy color centers in diamonds [7]), spin ensembles in ferromagnetic YIG materials have higher spin densities (~2.1 × 10 22 cm −3 ) and a low dissipation rate, and can be fully polarized below the Curie temperature (~559 K) [8]. Experimentally, a strong coupling between spin ensembles and microwave cavity modes has been achieved in YIG samples [9], which is extremely difficult for spin ensembles in paramagnetic materials. Many exotic phenomena can be observed in cavity-magnon systems, such as cavity spin-wave excitons [10,11], super strong coupling [12,13], mutual conversion between microwave photons and optical photons [13], synchronous spin-photon coupling [14,15], spin-wave dark mode [16], cavity spintronics [17], optical manipulation of cavity spin-wave systems [18], interlayer strong coupling of spin-waves [19], dynamical behavior of collective polaritons [20] and non-Hermitian physics. In addition, spin-wave quanta can also form a strong coupling with other quantum systems, such as superconducting qubits [21], phonons and optical whispering gallery cavities [22]. The above studies all approximate the wave color mode in the YIG sample as a linear harmonic oscillator. However, this approximation only holds under the low-power drive of cavity spin-wave systems. Under high-power driving, the nonlinear effects of the spin-wave mode in YIG will be obviously manifested, such as the bistable effect [23] and the chaotic effect [24].
Enlightened by the above-mentioned works about magnons, we construct a theoretical model of the hybrid quantum structure, which contains a YIG sphere coupling with a 3D microwave cavity, to explore nonclassical behaviors of the magnons' Kerr interaction. In this photon-magnon system, we diagonalize the Hamiltonian of the cavity spin-wave hybrid system with an ingenious transformation. Resorting to the concept of magnon coherent states, we investigate the quantum dynamics of the magnons and find that magnon Kerr nonlinearity can trigger nonclassical properties, such as negativity value of the Wigner function, quadrature squeezing effects, antibunching property and field nonclassicality in the magnon. These hybrid quantum systems may be useful in the preparation of magnon nonclassical states and the realization of a single magnon source.

The Model
Considering a cavity-magnonic system consisting of a cavity mode and a Kittel mode of the magnon, the magnetic dipole-dipole interaction of the Kittle mode with cavity photons could yield a nonlinear model. The Hamiltonian of this cavity-magnonic system can be written as (assume = 1) where H c in Equation (1) is the Hamiltonian of the cavity mode, i.e., The corresponding creation and annihilation operators for cavity mode c (with frequency ω c ) are a † and a, respectively. The interaction between the magnons and photons in the system is introduced by a YIG material sphere placed in the cavity. For the YIG sphere which is uniformly magnetized, its Hamiltonian can be written as two parts [25][26][27][28], i.e., the Zeeman energy and the magneto crystalline anisotropic energy [29]: Here, B 0 represents the applied static magnetic field, and τ and µ 0 are the volume and vacuum permeability of the YIG sphere, respectively. M m is the magnetization of the Kittel mode in the YIG sphere, and B 0 = B 0 e z is the magnetic field along the z direction. H an m represents the energy caused by the anisotropy of the magnetic medium inside the YIG sphere. In case the crystal axis [100] of the YIG sphere is also along the z direction, the anisotropic field H an m is [30][31][32][33][34][35][36][37][38][39] H an m = 2 γ g S m z K m an µ 0 M 2 τ e z , where K m an is the dominant first-order anisotropic coefficient and M is the saturation magnetization. After the Holstein-Primakoff transformation, the simplified Hamiltonian with respect to the magnon becomes where symbol m † (m) represents the magnon creation (annihilation) operator of cavity spin-wave quanta, and the corresponding frequency is . γ g is the gyromagnetic ration. The final term K m † m 2 characterizes the Kerr effect of the magnons. The coefficient K can be calculated by K = − γ 2 g K m an M 2 τ , which represents the intensity of the magnon Kerr nonlinearity. The magnon-photon interaction Hamiltonian is the cavity mode is polarized in the x direction, with V c being the volume of the cavity. Under the RWA in the second quantization, the interaction Hamiltonian can be simplified as Then, the total effective Hamiltonian of the photon-magnon system [30][31][32][33][34][35][36][37][38][39] is In order to diagonalize the Hamiltonian, we use a new pair of bosonic operators A 1 and A 2 [40][41][42], and its expressions are: where A 1 (t) and A 2 (t) are boson operators that satisfy the usual bosonic commutation Next, we need to substitute straightforwardly Equations (7) and (8) into Equation (6) and employ the RWA, i.e., drop the high-frequency fast oscillation term. Thus, the effective Hamiltonian can be written as where ω = (ω c + ω m )/2. For the purpose of simplifying quantum dynamics, let us recur to the Heisenberg equations with respect to A 1 (t) and A 2 (t) as follows where (7)-(11) yield the following expressions For the sake of the subsequent analysis, we adopt the product state of the Glauber coherent state as the initial state in the whole system, i.e., |α a , α m = |α a |α m ; it can also be written as In the following, we will use these initial states to analyze the nonclassical properties of magnons due to Kerr-nonlinearity.

Kerr-Nonlinearity Trigger Nonclassical Properties
I. Negativity value of Wigner function. The Wigner function [43][44][45] has turned out to be remarkably useful in quantum optics, particularly in the characterization and visualization of nonclassical fields. The Wigner function was first introduced by Wigner in 1932 to the calculate quantum correction to a classic distribution function of a quantum-mechanical system [45]. The Wigner function is a real function with a quasiprobability distribution property, invariance of coordinate and momentum exchange, and it can take positive values, which can be analogous to the distribution of the Hamiltonian system in classical statistical mechanics. However, the Wigner function can also take negative values, where the occurrence of negative values is not analogical, which is an important manifestation of nonclassical characteristics. Moreover, the Wigner function does not represent the joint measurement distribution of coordinates and momentum of particles in real phase space; that is, it does not describe the probability distribution of coordinates and momentum of particles at the same time. When the Wigner function presents a negative value, it can clearly reflect the nonclassical property of quantum states. The quantum interference structure of the Wigner function directly reflects the nonclassical characteristics of quantum superposition states. Therefore, the negative part of the Wigner function of a specific quantum state can be used as a symbol of nonclassical behavior to describe the quantum interference effect in a concise and vivid way.
We then investigated the presence of the negative value for the Wigner function using the definition given by [43][44][45] where α is a complex variable, ρ m = |α m α m |, and D(α) = e (αm † −α * m) is the displacement operator for Kittel mode m.
The outcome of the investigation is plotted in Figure 1, which clearly shows the existence of the negative value for the Wigner function when the initial states are ordinary Glauber coherent states. As observed in Figure 1(a-d), the nonlinear effect caused by the Kerr interaction makes the Wigner function of the coherent state negative in some regions of the phase space, thus indicating that the quantum state at this time is a nonclassical state. For the fixed parameter set, the negative volume of the Wigner function may increase, and the nonclassicality of the magnons' quantum state may change sharply when the time t rises. II. Squeezing effect. In quantum theory, the squeezed effect is an important nonclassical phenomenon that is peculiar to the quantum field. The squeezed state reflects the nonclassical characteristics of the bose field by means of a noise level that is even lower than that of the coherent state light field (laser field); that is, the noise fluctuation of an orthogonal phase component in the squeezed light is lower than that of the corresponding orthogonal phase component in the laser field. Therefore, in practical applications, if this component is used to transmit information, the signal-to-noise ratio can be higher than that of the coherent state light field, which has a broad application prospect in optical communication, weak signal detection and quantum nondestructive measurement.
In order to attain the squeezing effect, which delineates the nonclassical property of a boson field in the context of the fluctuations in the quadratures X 1m and X 2m of the field, we have the quadrature components of the magnon defined as [45][46][47] Resorting to Ref. [46], we attain the squeezing coefficients Substituting Equation (15) into Equation (16), one then straightforwardly obtains It is evident that when t = 0, one clearly observes S im (t) = 0, which indicates that there is no orthogonal squeezed phenomenon, since nonlinear photon-magnon interactions are not introduced. For the case of S im (t) < 0, it implies that the quadrature component of magnon X im has been squeezed. In order to show this more clearly, we have drawn a three-dimensional display by choosing g = 0.3, α a = 4 and α m = 2 in Figure 2. When taking into account the ordinary Glauber coherent state as the initial input field state, we find that the output state of the magnon still displays a periodic orthogonal squeezed behavior with the evolution of time t, as shown in Figure 2a,b. In this sense, we see that the squeezing effect can be revealed by calculating either the S 1 (t) or S 2 (t) of the X im component, and one of the two squeeze coefficients S 1 (t) and S 2 (t) always preserves a negative value. As a result, the magnon squeezing behavior exists all the time; it should be mentioned that each has a periodic oscillatory behavior. Simultaneously, as the Kerr coefficient rises, the degree of squeezing is also enhanced. For a fixed time t, the degree of squeezing may increase rapidly, and the fluctuations in the quadrature S im (t) of the magnons may change sharply by boosting the Kerr-nonlinearity strength. Therefore, the degree of the squeezing can be manipulated by selecting a suitable nonlinear strength K.
III. Antibunching phenomenon. The antibunching effect is a remarkable manifestation of nonclassicality. By using the second-order correlation function concept, we can quickly and efficiently identify different quantum statistical behaviors. According to quantum statistical methods, g (2) (0) = 1 implies that magnons are randomly distributed, and we cannot seek out the orientation of a single magnon. g (2) (0) > 1 indicates that magnons are inclined to come in flocks and tend to stay in the same place. g (2) (0) < 1 means that magnons do not like to be together and tend to separate or emerge in different places. In other words, g 2 (0) can reflect the repulsive property or attractive features of two magnons. The second-order correlation function g (2) (0) for a single-mode field is defined as [48][49][50] g (2)  In order to clearly see that the Kerr strength coefficient affects the quantum statistical properties of the single-mode magnon coherent state, we have carried out a direct numerical calculation of g (2) (0) for a single-mode field with respect to the Kerr strength K and dimensionless time t. As shown in Figure 3, we can clearly see the bunching phenomena and the fact that the magnon obeys the super-Poissonian quantum statistics discipline when the Kerr coefficient is smaller than 0.0033. When K > 0.0037, we can obviously observe that the magnon exhibits an antibunching effect that occupies a major area in Figure 3. In the cross regime of 0.0033 < K < 0.0037, the second-order correlation function of the magnon appears more complex when the Kerr strength K varies. It is easy to see that, for a single-mode field, g (2) (0) may decline rapidly when the Kerr-nonlinearity increases. For the sake of confirming this issue, we have plotted g (2) (0) as a function of the Kerr strength K with t = 9 and t = 10, as shown in Figure 4. In order to see this clearly, we draw an auxiliary line (corresponding to g (2) (0) = 1) and also mark the transition point position (K = K t ≈ 0.00332 and 0.00369). In the range 0.003 < K < K t , the second-order correlation function is bigger than 1; however, g (2) (0) of the magnon in the range K t < K < 0.006 is less than 1. Obviously, the magnon number distribution is a Poisson distribution at the transition point K t . This implies that the statistical properties of the magnon can be controlled by varying the Kerr-nonlinearity intensity in the photon-magnon system.  IV. Field nonclassicality. Similar to the light field, the magnon field is also one of the very important bose fields. The nonclassical states of these bose fields play an increasingly important role as a tool for exploring fundamental quantum theory and quantum technologies. Generally speaking, the nonclassical behavior of the bose field includes squeezing, antibunching, sub-Poissonian statistics, particles' number distribution and so on. The nonclassical nature of describing and quantifying quantum states is a central problem in quantum optics. At present, there are various criteria for detecting nonclassicality and several important nonclassical measures to distinguish nonclassical behavior. However, nonclassical ordering may be inconsistent under different metrics. Due to the complex and subtle characteristics of nonclassicality quantification, these nonclassicality quantification methods have different advantages and disadvantages, and their scope of application is also sharply distinct. As a consequence, it is impossible to obtain a universal measure that is applicable to all targets; that is to say, the above nonclassical measures cannot charac-terize all aspects of field nonclassicality either. For the sake of exploring the dynamics of the bose field from as many angles as possible to obtain a more complete picture of the dynamics, it is necessary to seek new nonclassical measures based on different perspectives. Fortunately, by resorting to the concept of generalized Wigner-Yanase skew information, the nonclassicality of the field is a nicely effective way to measure the nonclassical property of magnons. According to references [51,52], its definition can be written as Here, it is necessary to provide a brief description of field nonclassicality. As everyone knows, the Glauber coherent state |α appears to be the closest to the classical behavior of a bose field, and the field nonclassicality of the coherent state is 1/2. In other words, if a state is a pure single-mode quantum state, this pure quantum state's minimum value of field nonclassicality is N(ρ) = 1/2. A natural result is that field nonclassicity N(ρ) ≤ 1/2 for any single-mode classical state, and we can draw the small conclusion that the quantum state is nonclassical when N(ρ) > 1/2. Taking the dimensionless time t and Kerr coefficient K into account, we draw Figure 5 to observe the quantum dynamics of the field nonclassicality. In Figure 5, we can discover that it can exhibit field nonclassicality in most cases. Simultaneously, from the phase diagram, the nonclassical properties of the magnon show a quasiperiodic variation with time and the Kerr coefficients: for a fixed time, the nonclassical properties increase rapidly as the Kerr coefficients become larger. In other words, the increase in time seems to relax the requirement for the Kerr coefficient. It is worth mentioning that we have no choice of t < 20; this is because the Kerr interactions do not accumulate enough, and significant changes will not happen.

Conclusions
In conclusion, we have investigated a hybrid system consisting of a YIG sphere and a three-dimensional rectangular microwave cavity, where the two subsystems are coupled by the fundamental mode magnetic field in the microwave cavity. We diagonalized the Hamiltonian of the photon-magnon system with an ingenious transformation. More interestingly, by using the Wigner function, quadrature squeezing operator, second-order correlation function and field nonclassicality, we showed the Kerr-nonlinearity-triggered nonclassicality of magnons in a photon-magnon coupling system. It has been found that the degree of nonclassicality can be manipulated by adjusting the Kerr nonlinear strength.