cavity magnonics magnon polariton strong coupling hybrid quantum systems

# Cavity Magnonics: Magnon-Polariton Strong Coupling, Hybrid Quantum Systems, and Exceptional Points

## 1. Introduction to Cavity Magnonics

Cavity magnonics is the emerging field studying coherent coupling between magnons (collective spin-wave excitations in ferromagnetic materials) and photons confined in microwave cavities, creating hybrid magnon-photon quasiparticles (magnon-polaritons). This hybrid platform uniquely combines the tunability and long coherence times of ferromagnetic spin waves with the ease of control and measurement of microwave photons, enabling exceptional sensitivity for sensing, efficient frequency conversion, and exploration of non-Hermitian physics including exceptional points and topological phenomena.

A typical cavity magnonics system consists of a yttrium iron garnet (YIG) sphere or thin film placed inside a 3D or 2D microwave cavity (typically operating at 5-10 GHz). When the magnon frequency (tunable via applied magnetic field) is brought into resonance with the cavity resonance, strong magnon-photon coupling occurs, producing a characteristic avoided-crossing (anticrossing) in the microwave transmission spectrum—the signature of coherent hybrid oscillations.

This article develops cavity magnonics from first principles, covering the magnon and photon Hamiltonians, magnon-photon coupling mechanisms, Tavis-Cummings and Jaynes-Cummings frameworks, magnon-polariton eigenmodes, cooperativity and strong coupling criteria, dissipative effects and non-Hermitian physics including exceptional points, and applications in quantum transduction and sensing. A Python solver computes transmission spectra and polariton anticrossings.

## 2. Ferromagnetic Resonance and Magnon Modes

In a ferromagnetic material with saturation magnetization M_s, spin-wave excitations (magnons) emerge as collective precessions of the magnetization around an equilibrium direction. The precession frequency in an applied magnetic field B₀ (Zeeman field) is given by the Larmor precession:

$$\omega_m = \gamma g_s \mu_B B_0 / \hbar = \gamma_m B_0,$$

where γ_m is the gyromagnetic ratio for magnons (~1.76 × 10¹¹ rad/(T·s) for electrons). For YIG, a ferrite with exceptional magnetic properties (low damping, large M_s), typical ferromagnetic resonance frequencies are 5-10 GHz at fields B₀ ~ 0.2-0.4 T.

The magnon Hamiltonian for a single mode in the cavity is:

$$H_m = \hbar \omega_m (m^\dagger m + 1/2),$$

where m and m† are the magnon annihilation and creation operators obeying bosonic commutation relations [m, m†] = 1, analogous to photon operators but for magnons.

## 3. Microwave Cavity Photons

A 3D microwave cavity (e.g., a copper box with dimensions ~cm) confines electromagnetic radiation with resonant frequencies determined by boundary conditions:

$$\omega_c = \frac{c}{n} \sqrt{\left(\frac{m\pi}{L_x} ight)^2 + \left(\frac{n\pi}{L_y} ight)^2 + \left(\frac{p\pi}{L_z} ight)^2},$$

where c is the speed of light, n is the refractive index (~1 for vacuum/air), m, n, p are mode integers, and L_x, L_y, L_z are cavity dimensions. In practice, cavities are designed to have a single resonant mode (the fundamental mode) at the desired frequency, with higher-order modes suppressed.

The photon Hamiltonian is:

$$H_c = \hbar \omega_c (a^\dagger a + 1/2),$$

where a and a† are photon operators with [a, a†] = 1.

## 4. Magnon-Photon Coupling Mechanism

When a YIG sample is placed inside a microwave cavity, the time-varying magnetic field of the cavity photons exerts a torque on the YIG magnetization, creating a coupling between magnons and photons. The coupling arises from the interaction energy:

$$H_{ ext{int}} = -\mu_0 M_s V_{ ext{YIG}} \mathbf{m} \cdot \mathbf{B}_{ ext{cav}},$$

where V_YIG is the YIG sample volume and B_cav is the cavity magnetic field, proportional to a + a†. The rotating-wave approximation leads to:

$$H_{ ext{int}} = \hbar g_m (a^\dagger m + a m^\dagger),$$

where g_m is the single-magnon-single-photon coupling strength:

$$g_m = \frac{\mu_0 M_s V_{ ext{YIG}}}{2V_{ ext{cavity}}} \sqrt{\omega_c \omega_m}.$$

Coupling strength scales as the square root of cavity and magnon frequencies, and inversely with cavity volume. Typical values: g_m ~ 1-100 MHz for reasonably sized samples and cavities.

## 5. Tavis-Cummings Hamiltonian and Magnon-Polariton Eigenmodes

The total Hamiltonian for the magnon-photon system is:

$$H = \hbar \omega_c a^\dagger a + \hbar \omega_m m^\dagger m + \hbar g_m (a^\dagger m + a m^\dagger),$$

known as the Tavis-Cummings Hamiltonian (or Jaynes-Cummings for single magnons). Diagonalizing this Hamiltonian yields the eigenmodes—magnon-polaritons:

$$\omega_{\pm} = \frac{\omega_c + \omega_m}{2} \pm \sqrt{\left(\frac{\omega_c - \omega_m}{2} ight)^2 + g_m^2}.$$

When ω_c = ω_m (resonance), the polariton frequencies split symmetrically:

$$\omega_{\pm}( ext{resonance}) = \frac{\omega_c + \omega_m}{2} \pm g_m = \omega_0 \pm g_m,$$

where ω₀ = ω_c = ω_m. The splitting 2g_m is the vacuum Rabi splitting, analogous to cavity QED but with magnons replacing atoms.

As the detuning |ω_c - ω_m| increases from zero, the avoided crossing becomes progressively smaller. Far from resonance (|ω_c - ω_m| >> g_m), the eigenmodes revert to approximately pure cavity and magnon modes.

## 6. Cooperativity and Strong Coupling

The regime of strong magnon-photon coupling is characterized by the cooperativity:

$$C = \frac{g_m^2}{\kappa \gamma_m},$$

where κ is the cavity photon decay rate and γ_m is the magnon damping rate. In the strong coupling limit (C >> 1), coherent oscillations between magnons and photons dominate over dissipation, and the magnon-polariton anticrossing is clearly resolved in transmission spectra.

Strong coupling requires:
- High coupling: Large g_m (achieved with larger samples, higher frequencies, or smaller cavities)
- Low cavity loss: High-quality-factor cavities (Q_c = ω_c / κ >> 1000)
- Low magnon damping: Minimal dissipation in YIG (damping parameter α ~ 10⁻⁴ for high-quality YIG, giving γ_m = α ω_m)

Typical values for strong coupling: C ~ 1-100 in well-designed systems.

## 7. Experimental Signatures: Magnon-Polariton Anticrossing

In transmission spectroscopy, the cavity response S₂₁(ω, B₀) is measured as a function of frequency ω and applied magnetic field B₀. The magnon frequency ω_m ∝ B₀ shifts linearly with field; the cavity frequency ω_c is independent of B₀. Sweeping B₀ brings ω_m through resonance with ω_c, producing a characteristic avoided crossing pattern:

• Far detuning (B₀ away from resonance): Two separate transmission peaks at ω ≈ ω_c and ω ≈ ω_m
• Near resonance: Transmission peaks split and anticross, with minimum transmission at the crossing point (where the upper and lower polaritons are equally excited)
• Strong coupling (C >> 1): Clear anticrossing with minimal transmission dip; Rabi splitting 2g_m easily resolved
• Weak coupling (C << 1): Weak anticrossing; peaks nearly touch

The width and shape of the anticrossing encode the coupling strength and dissipation rates, allowing extraction of g_m, κ, and γ_m from experimental data.

## 8. Dissipative Coupling and Non-Hermitian Physics

Real systems experience energy loss from both cavity decay (radiation into free space) and magnon damping (scattering into phonons, impurities). The effective Hamiltonian becomes non-Hermitian:

$$H_{ ext{eff}} = \hbar \omega_c a^\dagger a + \hbar \omega_m m^\dagger m + \hbar g_m (a^\dagger m + a m^\dagger) - i\frac{\hbar \kappa}{2} a^\dagger a - i\frac{\hbar \gamma_m}{2} m^\dagger m,$$

where the imaginary terms represent decay. Non-Hermitian Hamiltonians exhibit exotic phenomena absent in Hermitian systems:

1. Exceptional points (EPs): Coalescence points where two or more eigenvectors become identical (complete non-orthogonality). At an EP, the Hamiltonian can be non-diagonalizable, leading to superradiance, chiral behavior, and unidirectional transmission.

2. Chiral eigenvalue trajectories: Eigenvalues form closed loops (non-analytic behavior) as parameters are varied—impossible in Hermitian systems.

3. Enhanced sensitivity: Near exceptional points, small parameter changes produce large eigenvalue shifts, enabling ultrasensitive sensing.

For the magnon-photon system, an exceptional point occurs when:

$$\frac{\kappa - \gamma_m}{2} = \sqrt{g_m^2 - \left(\frac{\kappa + \gamma_m}{4} ight)^2},$$

which requires tuning the detuning ω_c - ω_m and damping rates appropriately.

## 9. Magnon-Photon Transduction

Cavity magnonics provides a platform for quantum transduction—converting quantum information between different physical domains. A microwave photon absorbed by the cavity-magnon system can, under proper conditions, be converted into a magnon, which then couples to other degrees of freedom (e.g., phonons, spin currents, optical photons).

Applications include:
- Magnon-to-microwave transduction: Converting low-frequency magnons to microwave photons for readout
- Microwave-to-optical transduction: Via magnons coupled to optical cavities
- Hybrid quantum processors: Combining superconducting qubits (coupled to photons) with magnonic systems for enhanced processing

## 10. Magnon-Polariton Energy Landscape and Dynamics

The eigenenergies of the magnon-polariton system are:

$$E_{\pm}(B_0) = \frac{\hbar(\omega_c + \gamma_m B_0)}{2} \pm \frac{\hbar}{2}\sqrt{(\omega_c - \gamma_m B_0)^2 + 4g_m^2}.$$

As B₀ varies, the eigenenergy landscape traces out two branches (upper and lower polaritons) that anticross at resonance. The separation at resonance (Rabi splitting) is 2ℏg_m.

Time-domain Rabi oscillations can be observed when a magnon or photon is excited, causing oscillations between the two polariton states at rate g_m.

## 11. Numerical Solver: Magnon-Polariton Transmission Spectra

import numpy as np
import matplotlib.pyplot as plt
from scipy.linalg import eigh

def magnon_photon_hamiltonian(omega_c, omega_m, g_m, kappa, gamma_m):
    """
    Non-Hermitian effective Hamiltonian for magnon-photon system.
    Returns eigenvalues (complex) and eigenvectors.
    """
    H = np.array([
        [omega_c - 1j*kappa/2, g_m],
        [g_m, omega_m - 1j*gamma_m/2]
    ], dtype=complex)
    
    evals, evecs = np.linalg.eig(H)
    return evals, evecs, H

def transmission_spectrum(freq_range, B_0_range, omega_c, gamma_m, g_m, kappa, gamma_m_base):
    """
    Compute S21 transmission spectrum vs frequency and magnetic field.
    """
    S21 = np.zeros((len(B_0_range), len(freq_range)))
    
    for ib, B_0 in enumerate(B_0_range):
        omega_m = gamma_m * B_0  # Magnon frequency proportional to field
        
        for iw, omega in enumerate(freq_range):
            # Admittance matrix (cavity-magnon coupling)
            Y = np.array([
                [omega - omega_c + 1j*kappa/2, -g_m],
                [-g_m, omega - omega_m + 1j*gamma_m_base/2]
            ], dtype=complex)
            
            # Transmission coefficient (simplified model)
            det = np.linalg.det(Y)
            S21[ib, iw] = np.abs(-g_m / det)
    
    return S21

def plot_magnon_polaritons():
    """
    Plot magnon-polariton anticrossing and transmission spectra.
    """
    fig, axes = plt.subplots(2, 2, figsize=(14, 10))
    
    # Parameters
    omega_c = 2 * np.pi * 10e9  # Cavity frequency 10 GHz (normalized)
    g_m = 2 * np.pi * 50e6      # Coupling 50 MHz
    kappa = 2 * np.pi * 10e6    # Cavity decay 10 MHz
    gamma_m_base = 2 * np.pi * 5e6   # Magnon damping 5 MHz
    
    # Magnetic field sweep
    B_0_range = np.linspace(0.1, 0.5, 100)
    gamma_m = 2.8e6  # Gyromagnetic ratio (arbitrary units)
    omega_m_range = gamma_m * B_0_range
    
    # Compute polariton frequencies
    omega_plus = []
    omega_minus = []
    
    for omega_m in omega_m_range:
        omega_p = (omega_c + omega_m) / 2 + np.sqrt(((omega_c - omega_m)/2)**2 + g_m**2)
        omega_m_eig = (omega_c + omega_m) / 2 - np.sqrt(((omega_c - omega_m)/2)**2 + g_m**2)
        omega_plus.append(omega_p)
        omega_minus.append(omega_m_eig)
    
    omega_plus = np.array(omega_plus)
    omega_minus = np.array(omega_minus)
    
    # Panel 1: Magnon-polariton anticrossing (energy landscape)
    ax = axes[0, 0]
    ax.plot(B_0_range, omega_plus / (2*np.pi*1e9), 'r-', linewidth=2.5, label='Upper polariton')
    ax.plot(B_0_range, omega_minus / (2*np.pi*1e9), 'b-', linewidth=2.5, label='Lower polariton')
    ax.axhline(omega_c / (2*np.pi*1e9), color='k', linestyle='--', linewidth=1, alpha=0.5, label='Cavity')
    ax.plot(B_0_range, omega_m_range / (2*np.pi*1e9), 'g--', linewidth=1, alpha=0.5, label='Magnon')
    ax.set_xlabel('Magnetic field B₀ (T)', fontsize=11)
    ax.set_ylabel('Frequency (GHz)', fontsize=11)
    ax.set_title('Magnon-Polariton Anticrossing: Avoided Level Crossing', fontsize=12)
    ax.legend(fontsize=10)
    ax.grid(True, alpha=0.3)
    
    # Panel 2: Rabi splitting vs coupling
    ax = axes[0, 1]
    g_m_vals = np.linspace(10, 500, 100) * 2*np.pi*1e6
    rabi_splitting = 2 * g_m_vals / (2*np.pi*1e9)
    ax.plot(g_m_vals / (2*np.pi*1e6), rabi_splitting, 'purple', linewidth=2.5)
    ax.fill_between(g_m_vals / (2*np.pi*1e6), 0, rabi_splitting, alpha=0.3, color='purple')
    ax.set_xlabel('Coupling strength g_m (MHz)', fontsize=11)
    ax.set_ylabel('Rabi splitting 2g_m (GHz)', fontsize=11)
    ax.set_title('Vacuum Rabi Splitting vs Coupling', fontsize=12)
    ax.grid(True, alpha=0.3)
    
    # Panel 3: Cooperativity regimes
    ax = axes[1, 0]
    C_values = np.logspace(-1, 3, 100)
    ax.loglog(C_values, C_values, 'r-', linewidth=2.5, label='C = 1 (boundary)')
    ax.fill_between(C_values, 1e-2, 1e2, where=(C_values < 1), alpha=0.2, color='red', label='Weak coupling (C<1)')
    ax.fill_between(C_values, 1e2, 1e4, where=(C_values > 1), alpha=0.2, color='green', label='Strong coupling (C>1)')
    ax.set_xlabel('Cooperativity C = g_m² / (κ γ_m)', fontsize=11)
    ax.set_ylabel('Regime indicator', fontsize=11)
    ax.set_title('Strong vs Weak Coupling Regimes', fontsize=12)
    ax.legend(fontsize=10, loc='upper left')
    ax.set_xlim([0.1, 1000])
    ax.set_ylim([0.01, 100])
    
    # Panel 4: Summary and exceptional points
    ax = axes[1, 1]
    ax.axis('off')
    summary_text = """
    Cavity Magnonics Summary

    Host material: YIG (Yttrium Iron Garnet)
    • Saturation magnetization: ~200 kA/m
    • Damping parameter α: ~10⁻⁴ (very low)
    • Gyromagnetic ratio: 1.76×10¹¹ rad/(T·s)

    Typical system parameters:
    • Cavity frequency: 5-10 GHz
    • Magnon frequency: 5-10 GHz (tunable via B₀)
    • Coupling strength: 10-200 MHz
    • Cavity Q-factor: ~1000-10000
    • Magnon damping: 1-10 MHz

    Coupling mechanism:
    • Magnetic interaction between YIG spin waves
      and cavity photons
    • Strength ∝ √(ω_c ω_m) V_YIG / V_cavity

    Magnon-polariton regimes:
    • Weak coupling (C << 1): Modes decouple
    • Strong coupling (C >> 1): Clear anticrossing
    • Near resonance: Maximum Rabi splitting

    Non-Hermitian physics:
    ✓ Exceptional points (EP)
    ✓ Chiral behavior
    ✓ Unidirectional transmission
    ✓ Enhanced sensitivity

    Applications:
    ✓ Quantum transduction
    ✓ Hybrid quantum systems
    ✓ Magnonic circuitry
    ✓ Sensitive magnetometry
    ✓ Frequency converters
    """
    ax.text(0.05, 0.95, summary_text, transform=ax.transAxes, fontsize=9,
            verticalalignment='top', family='monospace',
            bbox=dict(boxstyle='round', facecolor='lightcyan', alpha=0.8))
    
    plt.tight_layout()
    plt.show()

print("Cavity Magnonics Analysis")
plot_magnon_polaritons()

## 12. Multi-Magnon and Multi-Photon Interactions

Beyond the single-magnon single-photon coupling, systems with multiple magnons or photons exhibit richer dynamics:

• Magnon number >1: Collective magnon modes in different regions of the YIG sample couple to a single cavity photon mode (Tavis-Cummings regime with N magnons). The coupling strength scales as √N.
• Photon number >1: Higher Fock states of the cavity field interact with magnons, leading to photon blockade (suppression of multi-photon states due to nonlinear coupling).

## 13. Applications: Quantum Transduction and Sensing

### 13.1 Magnon-to-Microwave Quantum Transduction
Converting quantum information from low-frequency magnonic systems (GHz range) to microwave photons (3-10 GHz) enables efficient quantum state transfer between otherwise incompatible platforms.

### 13.2 Magnonic Sensing
The high sensitivity of magnon-polariton systems near exceptional points enables detection of magnetic field changes at extreme sensitivity (pT/√Hz ranges), surpassing conventional magnetometers.

### 13.3 Hybrid Quantum Processors
Coupling magnon systems to superconducting qubits via microwave cavities creates hybrid processors combining advantages of both platforms.

## 14. Challenges and Future Directions

Current challenges include:
- Magnon damping: Minimizing loss to maintain coherence
- Coupling efficiency: Maximizing g_m without increasing losses
- Scalability: Building arrays of coupled magnon-photon systems
- Integration: Incorporating YIG with superconducting circuits on-chip

Future directions:
- Non-Hermitian topological effects in magnon systems
- Magnon-based quantum computing
- Room-temperature magnonic devices

## 15. Conclusion

Cavity magnonics represents a fertile platform for exploring strong light-matter coupling in a new regime—magnon-photon interactions. The hybrid magnon-polariton eigenmodes exhibit rich physics including avoided crossings, Rabi splittings, and exotic non-Hermitian phenomena like exceptional points. The coupling strength, dissipation rates, and cooperativity fundamentally determine the system behavior. Applications span quantum transduction, sensing, and hybrid quantum information processing. The presented numerical solver demonstrates magnon-polariton anticrossing spectra and transmission characteristics. As materials and fabrication techniques continue to advance, cavity magnonics will enable new quantum technologies at the interface of magnetism and microwave photonics.