---
title: Magnon Kerr Nonlinearity
url: https://www.emergentmind.com/topics/magnon-kerr-nonlinearity
type: topic
---

# Magnon Kerr Nonlinearity

Magnon Kerr nonlinearity is the quartic self-interaction of a magnon mode, usually the Kittel mode of a yttrium-iron-garnet (YIG) sample, appearing in effective Hamiltonians as \(K(m^\dagger m)^2\) or \(-K_m\,m^\dagger m^\dagger m m\). In cavity magnonics and magnomechanics it produces a magnon-population-dependent frequency shift, modifies hybrid magnon–photon–phonon spectra, and, under strong driving, generates effective parametric and squeezing terms. The nonlinearity is ordinarily traced to magnetocrystalline anisotropy in YIG, but later work also realized an effective magnon self-Kerr by dispersively coupling magnons to a superconducting qubit, bringing the Kerr scale into the quantum regime [2401.02251, 1609.07891, 2409.01590, 2602.19671].

## 1. Microscopic origin and model forms

A representative cavity-magnon Hamiltonian in the rotating frame is
\[
H_{\rm sys}
=
\sum_{j=1,2}
\Big[
\Delta_m\,m_j^\dagger m_j
+
g_j\,(m_j^\dagger c+c^\dagger m_j)
+
K_j\,(m_j^\dagger m_j)^2
\Big]
+
\Delta_c\,c^\dagger c
+
\Omega\,(c+c^\dagger),
\]
where \(c\) is the cavity mode, \(m_j\) are magnon modes, \(g_j\) are photon–magnon couplings, and \(K_j(m_j^\dagger m_j)^2\) is the magnon Kerr term. In the single-sphere case the same structure is commonly written as
\[
H=\omega_c a^\dagger a+\omega_m m^\dagger m+K(m^\dagger m)^2+g(a^\dagger m+a m^\dagger)+\Omega_d(m^\dagger e^{-i\omega_d t}+m e^{i\omega_d t}),
\]
so the quartic term directly shifts the magnon resonance by an amount proportional to magnon occupation [2401.02251, 1609.07891].

The microscopic source in YIG is magnetocrystalline anisotropy. In the macrospin description one starts from an anisotropy contribution quadratic in \(S_z\), applies the Holstein–Primakoff transformation \(S_z=S-m^\dagger m\), and obtains a quartic bosonic term. Two equivalent expressions quoted in the literature are
\[
K \simeq \frac{\mu_0\gamma^2 K_1}{2M_s^2V}
\]
and
\[
K=\frac{\mu_0 K_{\rm an}\gamma^2}{M^2V_m},
\]
which make explicit the inverse-volume scaling. In the thin-film formulation the anisotropy energy is written as \(U_{\rm anis}=K_u^*V_m\cos^2\theta\), and expansion to fourth order yields
\[
K=
\begin{cases}
\xi\,(\gamma^2\hbar K_u^*)/(2M_s^2V_m), & \text{in-plane magnetization},\\[4pt]
-(\gamma^2\hbar K_u^*)/(M_s^2V_m), & \text{out-of-plane magnetization},
\end{cases}
\]
so both shape anisotropy and magnetization orientation enter the effective Kerr coefficient [2401.02251, 1609.07891, 2601.22271].

A distinct route to magnon Kerr nonlinearity uses dispersive coupling to a superconducting qubit. After Schrieffer–Wolff elimination, the effective magnon Hamiltonian becomes
\[
H_{\rm eff}=\omega_m' a^\dagger a+\frac{K}{2}a^{\dagger 2}a^2,
\qquad
K=-2\,\frac{g_{qm}^4}{\Delta_{qm}^3},
\]
with \(g_{qm}\) the effective qubit–magnon coupling and \(\Delta_{qm}\) the dispersive detuning. Here the quartic term is not attributed to YIG anisotropy alone, but to virtual qubit processes that generate a sizable self-Kerr in the magnon mode [2602.19671].

## 2. Sign control, scaling, and reported magnitudes

Reported Kerr quantities depend strongly on sample geometry, orientation, and whether one quotes a bare single-magnon coefficient or a drive-induced nonlinear shift.

| Platform | Reported Kerr quantity | Note |
|---|---|---|
| 1 mm YIG sphere in a 3D cavity | \(K/2\pi\sim1\,\mathrm{nHz}\) | Pump-induced shift observable when \(\langle m^\dagger m\rangle\sim10^7-10^8\) [1609.07891] |
| 0.25 mm diameter YIG sphere | \(K/2\pi\sim1-10\,\mathrm{kHz}\), \(|K|/\Gamma\sim10^{-3}-10^{-2}\) | \(g/2\pi\sim10-100\,\mathrm{MHz}\), \(\Gamma/2\pi\sim1\,\mathrm{MHz}\) [2401.02251] |
| 0.28 mm YIG sphere in cavity magnomechanics | \(\chi_m/2\pi\simeq-6.5\,\mathrm{nHz}\), \(\chi_{mp}/2\pi\simeq-5.4\,\mathrm{pHz}\) | Self-Kerr and magnon–phonon cross-Kerr extracted from bistability data [2206.14588] |
| 200 nm YIG film in a 3D loop-gap cavity | \(K_{\rm IP}/2\pi=-57.7\pm0.9\,\mathrm{nHz}\), \(K_{\rm OOP}/2\pi=+170\pm20\,\mathrm{nHz}\) | \(K\approx-3.1\pm0.8\,\mathrm{nHz}\) near \(\theta_H\approx40^\circ\) [2601.22271] |

The sign of the Kerr coefficient is an experimentally tunable quantity. For YIG spheres, rotating the bias field relative to the crystal axes yields \(B_0\parallel[100]\Rightarrow K_j>0\) and \(B_0\parallel[110]\Rightarrow K_j<0\), while the magnitude scales as \(1/V\propto D^{-3}\). In the thin-film system, \(K(\theta_M)\) changes sign continuously, crosses near \(\theta_M\approx40^\circ\), and reaches extrema for out-of-plane and in-plane magnetization. This suggests that geometry and field orientation are coequal control parameters for Kerr engineering rather than secondary perturbations [2401.02251, 2601.22271].

The same point appears in comparative anisotropy estimates. In a bulk YIG sphere \(K_u^*\simeq K_1\), giving \(K_{\rm sphere}\sim10^{-12}-10^{-11}\,\mathrm{Hz}\), whereas in a thin film the demagnetizing field yields \(K_{\rm shape}=\mu_0M_s^2/2\sim10\,\mathrm{kJ/m^3}\), so \(|K_u^*|\gg K_1\). Together with the \(1/V_m\) scaling, this provides a direct route to larger \(|K|\) without abandoning strong magnon–photon coupling [2601.22271].

## 3. Frequency shifts, nonlinear response, and experimental extraction

In the dispersive cavity-magnon experiment of Wang et al., the observable effect of the Kerr term was a pump-induced magnon-frequency shift and a smaller cavity-frequency shift. In the dispersive limit \(|\Delta|=|\omega_c-\omega_m|\gg g\), the steady-state shift obeys
\[
[\Delta\omega_m^2+(\gamma_m/2)^2]\Delta\omega_m=cP,
\]
equivalently
\[
\Delta\omega_m^3+(\gamma_m/2)^2\Delta\omega_m-cP=0.
\]
At very small \(P\), \(\Delta\omega_m\) is linear in \(P\); at high power, \(\Delta\omega_m\propto(cP)^{1/3}\). The cavity shift is
\[
\Delta\omega_c=\frac{2g^2}{\Delta^2}\,K\,\langle m^\dagger m\rangle
\simeq
\frac{2g^2}{\Delta^2}\,\Delta\omega_m.
\]
For \(P=11\,\mathrm{dBm}\) the cavity central frequency shifted by \(\simeq0.7\,\mathrm{MHz}\), while the magnon shift was \(\simeq7\,\mathrm{MHz}\) [1609.07891].

In cavity magnomechanics, the Kerr coefficient was extracted by fitting the reflection spectrum and tracking a power-dependent magnon detuning
\[
\tilde\Delta_m=\omega_m-2\chi_m|\langle m\rangle|^2-\omega_d.
\]
The mechanical-frequency shift was then decomposed into the conventional radiation-pressure spring shift and a magnon–phonon cross-Kerr term \(\chi_{mp}|\langle m\rangle|^2\). Pump sweeps upward and downward produced hysteresis in both \(\omega_m\) and \(\omega_b\), identifying a Kerr-modified bistable regime rather than a purely linear magnomechanical response [2206.14588].

In the thin-film experiment, the coupled equations
\[
\dot c=-i\Delta_c c-(\kappa_{\rm tot}/2)c-ig m+\sqrt{\kappa_{\rm ext}}\,c_{\rm in},
\qquad
\dot m=-i[\Delta_m+2K\langle m^\dagger m\rangle]m-(\kappa_m/2)m-igc+\sqrt{\kappa_m}\,m_{\rm in}
\]
were fitted to power-dependent \(S_{21}\). At \(P=25\,\mathrm{dBm}\), \(N_m\sim10^{14}\) and the induced magnon shift reached several MHz. In the qubit-assisted quantum-regime experiment, the dispersive coefficients were determined by direct spectroscopy and by fitting Kerr-induced phase-space shearing: \(\chi/2\pi=4.9\pm0.1\,\mathrm{MHz}\), \(\delta/2\pi=0.36\pm0.02\,\mathrm{MHz}\), and therefore \(K/2\pi\approx-0.72\,\mathrm{MHz}\) [2601.22271, 2602.19671].

## 4. Interference physics, blockade, and nonreciprocity

One of the most direct uses of magnon Kerr nonlinearity is unconventional photon blockade. In the weak-drive truncation,
\[
g^{2}(0)\simeq \frac{2|C_{20}|^2}{|C_{10}|^4},
\]
and photon antibunching follows from destructive interference between a direct two-photon path, \(|0,0\rangle\to|1,0\rangle\to|2,0\rangle\), and a magnon-assisted path, \(|0,0\rangle\to|1,0\rangle\to|0,1\rangle\to|1,1\rangle\to|2,0\rangle\). Because the intermediate states shift by \(+2K\) for \(K>0\) and by \(-2K\) for \(K<0\), the cancellation condition depends on the Kerr sign:
\[
\frac{g^2K}{\Delta_m'+2K}+(\Delta_c'+\Delta_m'+K)(\Delta_m'+K)=\Omega^2,
\qquad
\Delta'=\Delta-i\Gamma/2.
\]
Solving this gives
\[
\Delta_{\rm opt}\simeq \pm(\sqrt{3}/6)\Gamma,
\]
with the sign set by \(\operatorname{sign}(K)\). In the one-sphere proposal this produces nonreciprocal unconventional photon blockade; in the two-sphere version, symmetric \(g_j\) and \(K_j\) yield reciprocal blockade, whereas asymmetric couplings or asymmetric Kerr strengths produce nonreciprocity [2401.02251].

The same quartic nonlinearity also drives high-order frequency conversion. In bichromatically driven cavity–magnon systems the output field acquires sidebands at \(\omega_\ell\pm n\Delta\), with exact spacing \(\Delta=\omega_p-\omega_\ell\). The effective nonlinear shift is \(\delta\omega_b=2\mathcal K|\beta|^2\), and the sideband amplitudes scale approximately as \((\delta\omega_b)^n\). Raising the microwave drive to an optimal \(\sim12\,\mathrm{mW}\) populates many sidebands, while changing \(\Delta\) directly tunes the comb spacing. In the nonreciprocal two-cavity configuration, the effective Kerr scale is \(K_{\rm eff}\simeq2K n_m\), and the reported isolation ratios near optimal detuning \(\Delta_p\simeq-0.3\Delta_c\) are \(I_t\approx5.6\,\mathrm{dB}\) for the first-order sideband, \(I_\tau\approx9.9\,\mathrm{dB}\) for the second-order sideband, and \(I_\ell\approx30\,\mathrm{dB}\) for the third-order sideband [1806.08289, 2106.09542].

In hybrid magnomechanics, Kerr-induced frequency shifts reshape transparency and dispersion. The effective susceptibility
\[
\chi_{n_1}(\omega)=\big[\kappa_{n_1}+i(\Delta_{n_1}+2K|n_{1s}|^2-\omega)\big]^{-1}
\]
moves the transparency windows differently for opposite propagation directions, producing asymmetric Fano line shapes, multiple transparency dips, and large directional group-delay contrasts. Under the stated strong-coupling and mechanical-\(Q\) conditions, forward and backward group delays can differ by \(>100\,\mu\mathrm{s}\) [2606.13412].

Kerr nonlinearity can also be leveraged through non-Hermitian engineering. In the \(\mathcal PT\)-symmetric cavity-magnonic proposal, balancing gain and loss near the exceptional point strongly amplifies the effective nonlinearities, so that perfect magnon blockade is obtained even for \(\chi/\kappa_m\approx10^{-2}\), and photon blockade appears simultaneously although there is no bare photonic Kerr term. In a polaritonic formulation, the Kerr-shifted spectrum
\[
E_n=n\Delta_+ + (K/4)n^2
\]
supports \(2\)–\(5\) polariton blockade with \(F_p(n)>0.95\) for \(K/g=0.05\), \(\kappa/g=0.001\), \(\Omega/g=0.005\)–\(0.01\), and \(G/g=0.005\)–\(0.01\) [2209.04228, 2603.16089].

## 5. Kerr dressing, squeezing, entanglement, and engineered couplings

Under strong driving, the magnon Kerr term does more than shift a resonance: it generates effective parametric interactions. In the photon–phonon squeezing protocol, linearization around a large magnon amplitude and the squeezing transformation
\[
U(r,\theta)=e^{\frac r2(e^{-i\theta}m^2-e^{i\theta}m^{\dagger2})},
\qquad
r=\tfrac12\operatorname{arctanh}\!\Bigl(\tfrac{2K}{\Delta_m}\Bigr),
\qquad
K=K_m|\langle m\rangle|^2,
\]
convert the Kerr term into a dressed detuning \(\Delta_m'=(\Delta_m-2K)/\cosh2r\) and renormalized couplings. After eliminating magnon fluctuations, the effective photon–phonon Hamiltonian becomes
\[
H_{\rm eff}=\delta\,a^\dagger a+\delta\,b^\dagger b+g_{\rm eff}(a^\dagger b^\dagger+ab),
\]
with
\[
g_{\rm eff}
=
-\,gGe^r
\Bigl[
\frac{\sinh r}{\Delta'_m-\omega_b}
+
\frac{\cosh r}{\Delta'_m+\omega_b}
\Bigr].
\]
Numerically, moderate \(r\sim0.2\)–\(0.3\) maximizes \(g_{\rm eff}\), and realistic parameters allow \(\sim10\,\mathrm{dB}\) of two-mode squeezing for \(r\approx0.25\) and \(g\approx G\approx0.1\,\omega_b\) [2409.01590].

A closely related construction appears in magnomechanical cooling. Two-tone driving generates a squeezing rate \(\lambda e^{i\theta}=K B_+B_-\), and the Bogoliubov mode
\[
\beta=b\cosh r+b^\dagger e^{i\phi}\sinh r
\]
diagonalizes the magnon sector. The optical damping rates
\[
\Gamma_\pm=2\mathcal G^2\,\Re\chi_\pm
\]
then give
\[
\Gamma_{\rm opt}=\Gamma_- - \Gamma_+,
\qquad
n_{\rm final}=\frac{\gamma_b n_{\rm th}+\Gamma_+}{\gamma_b+\Gamma_{\rm opt}}.
\]
In the sideband-unresolved regime, the analytic optimum is
\[
\Delta_m'=\tfrac12\sqrt{\gamma_m^2+4\omega_b^2},
\qquad
\lambda_{\rm opt}=\frac{\gamma_m}{2},
\]
and the representative YIG-sphere parameters yield \(\Gamma_{\rm opt}/2\pi\sim10\,\mathrm{kHz}\) with \(n_{\rm final}\lesssim1\) [2504.06484].

Kerr nonlinearity also mediates steady-state entanglement. In the two-magnon cavity proposal, strong blue-detuned driving creates effective parametric couplings
\[
G_j=2\Delta_j\,\Re\langle m_j\rangle^2,
\qquad
F_j=2\Delta_j\,\Im\langle m_j\rangle^2,
\]
and the resulting quadratic fluctuation Hamiltonian supports magnon–photon squeezing and cavity-mediated magnon–magnon entanglement. For the quoted realistic parameters, numerical solutions of the Lyapunov equation yield logarithmic negativities \(E_{m_1m_2}\sim0.1\)–\(0.2\) at \(T=10\,\mathrm{mK}\) [1904.04167]. In a single ferrimagnetic crystal containing a Kittel mode and a higher-order magnetostatic mode, self-Kerr and cross-Kerr coexist through
\[
H_{\rm int}=K_b\,b^\dagger b^\dagger bb+K_c\,c^\dagger c^\dagger cc+G\,b^\dagger b\,c^\dagger c,
\]
and strong drives raise the effective nonlinearities to \(\tilde K_b/2\pi\approx15\,\mathrm{MHz}\), \(\tilde K_c/2\pi\approx24\,\mathrm{MHz}\), and \(\tilde G/2\pi\approx19.4\,\mathrm{MHz}\), with entanglement surviving up to \(T_e\lesssim0.15\)–\(0.2\,\mathrm{K}\) [2205.11865].

More elaborate hybridizations use Kerr magnons as squeezed intermediaries. In the NV-center proposal, the squeezed-basis tripartite coupling
\[
G_s
=
G_k\bigl(\cosh\xi_l\cosh\xi_r+\sinh\xi_l\sinh\xi_r\bigr)
=
\frac{G_k}{4}e^{\xi_l+\xi_r}
\]
gives, for \(\xi_l=\xi_r=\xi\),
\[
G_s/G_k=\tfrac14 e^{2\xi},
\qquad
C_s/C_k=\tfrac1{64}e^{6\xi},
\]
and perfect magnon blockade emerges when \(G_s\) is sufficiently large, numerically for \(\xi\gtrsim2\) [2509.01884]. In the driven spin–magnon problem, the Kerr-induced squeezing parameter \(r_k\) enhances the coupling as \(\mathcal G_k\approx g_k e^{r_k}\), invalidates the single-Kittel-mode approximation under strong pumping, and leads to bound-state-induced population trapping or persistent Rabi-like oscillation [2308.05927].

The quantum-regime experiment provides an explicit demonstration that a sufficiently large magnon Kerr can directly generate nonclassical magnon states. There, a 1-mm YIG sphere coupled dispersively to a superconducting qubit exhibited \(K/2\pi\approx-0.72\,\mathrm{MHz}\), and Kerr evolution produced quadrature variances of \(\sim0.8\), corresponding to \(\sim1.0\,\mathrm{dB}\) of squeezing, with mean magnon number less than one [2602.19671].

## 6. Bistability, soft modes, sensitivity enhancement, and conceptual boundaries

Beyond antibunching and squeezing, magnon Kerr nonlinearity reorganizes nonlinear dynamics at the semiclassical level. In cavity magnomechanics, mechanical bistability arises when magnetostriction, magnon self-Kerr, and magnon–phonon cross-Kerr act simultaneously; the experiment identified hysteresis in the magnon resonance, the mechanical resonance, and the mechanical linewidth under \(20\)–\(30\,\mathrm{dBm}\) pumping [2206.14588]. In easy-axis ferromagnets coupled to a microwave cavity, the Kerr term stabilizes the soft mode by opening a finite gap,
\[
\Delta_{\rm gap}=\sqrt{4g^2+(2K|b_0|^2)^2},
\]
and, in the simplest zero-drive picture,
\[
\Delta_{\rm gap}(K,g)\simeq 2^{1/3}g^{2/3}K^{1/3}.
\]
The same analysis predicts chaos and comb-like behavior near the mode crossing, with comb spacing
\[
\delta\omega \simeq \frac{2K|b|^2}{2\pi},
\]
and reports agreement with \(0.34\,\mathrm{GHz}\) sidebands [2602.12550].

Sensitivity to very small Kerr shifts can be increased through exceptional-point physics. In the coherent-perfect-absorption proposal, imposing the CPA condition yields an effective non-Hermitian Hamiltonian that can host an EP3. Near that point, a tiny Kerr-induced shift \(\Delta_K\ll\Gamma\) produces cube-root eigenvalue splitting,
\[
\Omega_\ell=\Omega_{\rm EP3}+\lambda_\ell(\Delta_K/\gamma_2)^{1/3}\gamma_2+\cdots,
\]
so the observable dip separation obeys
\[
\delta\omega_p\propto(\Delta_K)^{1/3}.
\]
The paper’s stated implication is that a Kerr shift far below the linewidth can still be converted into a much larger spectral splitting in the cavity transmission [2211.08922].

A recurrent boundary condition in this literature is the distinction between bare Kerr strength and effective nonlinear response. Bulk YIG can have an intrinsic single-magnon Kerr that is extremely small, yet the same platform can display large pump-induced shifts, strong blockade, or measurable squeezing once one exploits large occupations, shape anisotropy, \(\mathcal PT\)-symmetric enhancement, exceptional-point amplification, or qubit-mediated dispersive engineering [1609.07891, 2601.22271, 2209.04228, 2211.08922, 2602.19671]. This suggests that “magnon Kerr nonlinearity” in current research denotes not only a microscopic quartic coefficient, but also a design principle for converting weak intrinsic anharmonicity into experimentally resolvable nonlinear magnon dynamics.

Source: https://www.emergentmind.com/topics/magnon-kerr-nonlinearity