Magnon Photon Polariton Cavity Spintronics Strong Coupling Hybridization Dispersometer

# Cavity Magnon-Polariton Spintronics: Strong Coupling Kinetics, Vacuum Rabi Splitting, and Dispersive Polariton Dynamics

## 1. Introduction to Cavity Magnonics and Magnon-Polaritons

Cavity magnonics—the strong coupling of microwave photons in superconducting resonators with magnons (spin-wave quanta) in ferrimagnetic materials—represents a revolutionary hybrid platform combining the coherence of superconducting circuits with the rich spin dynamics of magnonic systems. When a ferrimagnetic resonator (typically yttrium iron garnet, YIG) is placed inside a microwave cavity, the magnon-photon interaction Hamiltonian creates a magnon-polariton: a hybrid quasiparticle sharing both photonic and magnonic character. The emergence of vacuum Rabi splitting—a frequency splitting proportional to the vacuum coupling strength—and the ability to operate in ultra-strong coupling regime (where the coupling strength g_m exceeds 10% of the bare frequencies) have enabled exceptional phenomena: coherent magnon-photon exchange, topological magnon edge modes, nonlinear magnon dynamics, and quantum sensing at millikelvin temperatures. This section develops the complete magnon-photon strong-coupling formalism, derives polariton eigenfrequencies, and explores the role of dissipation in determining cooperativity and coherence.

## 2. Magnons: Spin-Wave Quanta and Ferrimagnetic Resonance

Magnons are quantized collective spin-wave excitations in ferri- and ferromagnetic materials. In a ferrimagnetic crystal with magnetization mathbf{M} precessing about an external magnetic field mathbf{H}_{text{ext}}, the linearized equation of motion is the Landau-Lifshitz equation:

where mathbf{m} is the unit magnetization vector and gamma = g_e mu_B / hbar is the gyromagnetic ratio. For a uniform ferrimagnetic resonator, the ensemble of spin-wave modes at wavevector k exhibits the Kittel dispersion:

where D is the dipolar-exchange stiffness and M_s is the saturation magnetization. For uniformly magnetized samples (Damon-Eshbach modes), the zero-wavevector mode (uniform precession) dominates cavity coupling, with resonant frequency:

Typical values for YIG at room temperature: omega_m/(2pi) approx 5-20 GHz (tunable by H_{text{ext}} from 0.2 to 1 T), magnetic damping alpha approx 10^{-5} (exceptionally low), and magnon lifetime tau_m = 1/(alpha omega_m) approx 10-100 microseconds.

## 3. Cavity Quantum Electrodynamics (Cavity QED) for Magnons

A superconducting cavity (typically coplanar waveguide or 3D cylindrical resonator) supports resonant electromagnetic modes at frequencies ω_c. When a YIG sphere or thin film is placed inside the cavity, the resonator's capacitive and inductive elements couple to the magnetic dipole moment mathbf{m} = m_0 mathbf{e}_z of the ferrimagnet. The magnon-photon interaction Hamiltonian (rotating-wave approximation) is:

where a and a^dagger are photon annihilation/creation operators, m and m^dagger are magnon operators, and g_m is the vacuum coupling strength. The full system Hamiltonian is:

The coupling strength depends on the magnetic quality factor of the resonator and the sample geometry:

where chi_m is the magnetic susceptibility and V_{text{overlap}} is the overlap integral between the cavity field and magnetization. Engineering magnon-photon coupling typically yields g_m/(2pi) approx 1-100 MHz, comparable to decay rates kappa/(2pi) approx 1-10 MHz (photon decay) and gamma_m/(2pi) approx 0.1-1 MHz (magnon decay).

## 4. Magnon-Polariton Eigenfrequencies and Avoided Crossing

Diagonalizing the Hamiltonian H = hbar omega_c a^dagger a + hbar omega_m m^dagger m + hbar g_m (a^dagger m + a m^dagger) yields two polariton eigenfrequencies (in the absence of dissipation):

ω_± = (ω_c + ω_m)/2 ± √(((ω_c - ω_m)/2)² + g_m²)

When on-resonance (omega_c = omega_m = omega_0), the frequencies are simply:

This symmetric splitting, with magnitude 2g_m (the vacuum Rabi splitting), is the hallmark of strong coupling. For off-resonance detuning Delta = omega_c - omega_m, the splitting exhibits an avoided crossing: near resonance, the polaritons anticross with a minimum gap proportional to g_m; far detuning, they approach the bare frequencies ω_c and ω_m.

## 5. Cooperativity and Strong-Coupling Regime

The cooperativity parameter C quantifies the ratio of coherent magnon-photon exchange rate to incoherent dissipation:

C = g_m² / (κ γ_m)

where kappa = omega_c / Q_c is the photon decay rate (inverse of cavity quality factor) and γ_m is the magnon decay rate. The regimes are:

• Weak coupling (C < 1): Dissipation dominates; no clear polariton peaks
• Strong coupling (1 < C < 10): Polariton splitting resolvable; coherent exchange on timescale sim 1/g_m
• Ultra-strong coupling (C > 10): Rabi oscillations fully coherent; nonlinear effects emerge

Modern cavities achieve C = 100-10,000, with recent demonstrations exceeding C > 10^5. The condition for strong coupling is equivalently:

ensuring that the polariton splitting exceeds both decay rates.

## 6. Ultra-Strong Coupling and Novel Physical Regimes

When g_m becomes comparable to the bare frequencies—ultra-strong coupling g_m / omega sim 0.01-0.1—the rotating-wave approximation breaks down and counter-rotating terms become significant:

This non-Hermitian interaction creates extraordinary phenomena:

1. Exceptional points (EPs): Complex eigenvalues and non-orthogonal eigenstates, offering topological protection
2. Ground-state polariton pairing: The ground state acquires vacuum photons and magnons, breaking the Lamb-shift picture
3. Dynamical decoupling: Magnon excitations can coherently de-excite photons even without resonance

Reaching ultra-strong coupling requires engineered interactions: placing YIG cylinders at cavity field antinodes, using superconducting qubits with exceptionally high polarizability, or hybrid architectures combining multiple magnons.

## 7. Dissipation and Polariton Linewidth

In the presence of dissipation, the polariton linewidth is determined by the Lorentzian broadening from decay rates κ and γ_m. The lower polariton (LP) and upper polariton (UP) linewidths are (near resonance):

The lower polariton exhibits linewidth narrowing when C gg 1: even though individual decay rates are κ and γ_m, coherent mixing can suppress the LP decay by destructive interference. This is central to using cavity magnonics for quantum sensing: the low-mass lower polariton has intrinsic linewidth gamma_{text{LP}} approx sqrt{kappa gamma_m} (geometric mean), enabling ultrasensitive detection.

## 8. Cavity-Magnon Transmission Spectroscopy and S-Parameters

Experimental characterization relies on scattering-parameter measurement. A transmission measurement through the cavity (input port → cavity → output port) yields the transmission coefficient:

At resonance (omega_c = omega_m = omega_0), this simplifies to:

The transmission dip at cavity resonance exhibits characteristic avoided crossing: two minima separated by 2g_m, with oscillation amplitude modulated by dissipation. Fitting this spectroscopy extracts ω_c, ω_m, g_m, κ, and γ_m simultaneously.

## 9. Magnon-Polariton Coherence and Quantum-Limited Timescales

The upper and lower polaritons have coherence times:

For optimized cavities (kappa sim gamma_m), the lower polariton coherence time is:

Superconducting cavities with Q_c > 10^5 yield T_2^{text{LP}} > 100 microseconds at microwave frequencies, approaching superconducting qubit timescales. This enables protocols like magnon-to-qubit coherence transfer, where a qubit couples to a polariton as an intermediary, leveraging the polariton's large Hilbert space.

## 10. Nonlinear Magnon-Photon Dynamics

At high excitation amplitudes (photon populations bar{n} > 1), nonlinear effects dominate. The Kerr nonlinearity from magnon-magnon scattering creates frequency shifts:

where chi_{text{KNL}} sim 10^{-6} to 10^{-3} Hz per magnon (material-dependent). This nonlinearity enables:

• Parametric amplification: Pump the cavity at frequency 2omega to generate photon-magnon pairs
• Magnon lasing: Invert magnon population via microwave pumping; cavity feedback drives lasing
• Polariton blockade: Engineered nonlinearity prevents simultaneous excitation of multiple polaritons, creating photon/magnon antibunching

## 11. Applications: Quantum Sensing and Sensing Beyond Standard Limits

Cavity magnon-polariton systems excel at detecting tiny magnetic moment changes. The sensing sensitivity is:

where Gamma is the total system linewidth and τ is the measurement time. For Gamma sim 1 kHz and g_m sim 10 MHz, sensitivities approach Delta m sim 10^{-18} to 10^{-20} Bohr magnetons—quantum-limited by the magnon zero-point motion. Applications:

• Axion dark matter detection: Couple weakly interacting massive particles (WIMPs) magnetically to the cavity
• Gyroscope sensing: Magnon precession frequency encodes rotation rate
• Biosensing: Protein magnetic moments in solution

## 12. Python Numerical Solver: Polariton Spectrum and Transmission Dynamics

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

def magnon_photon_hamiltonian(omega_c, omega_m, g_m):
    """
    Construct the Jaynes-Cummings Hamiltonian for magnon-photon coupling.
    
    H = ℏ ω_c a†a + ℏ ω_m m†m + ℏ g_m (a† m + a m†)
    
    In matrix form (basis: |0,0⟩, |1,0⟩, |0,1⟩, |1,1⟩):
    """
    # Simplified 2x2 system (ignoring second excitations)
    H = np.array([
        [0, 0],
        [0, omega_c + omega_m]
    ], dtype=complex)
    
    # Coupling term
    H[0, 1] = g_m
    H[1, 0] = g_m
    
    return H

def polariton_frequencies(omega_c, omega_m, g_m):
    """
    Compute polariton eigenfrequencies ω_± from exact diagonalization.
    
    Returns: (omega_plus, omega_minus)
    """
    delta = (omega_c - omega_m) / 2.0
    radicand = delta**2 + g_m**2
    
    omega_avg = (omega_c + omega_m) / 2.0
    omega_plus = omega_avg + np.sqrt(radicand)
    omega_minus = omega_avg - np.sqrt(radicand)
    
    return omega_plus, omega_minus

def vacuum_rabi_splitting(g_m, omega_c, omega_m):
    """
    Compute vacuum Rabi splitting at resonance.
    At ω_c = ω_m, splitting is exactly 2*g_m.
    """
    if abs(omega_c - omega_m) < 1e-6:
        return 2 * g_m
    else:
        omega_plus, omega_minus = polariton_frequencies(omega_c, omega_m, g_m)
        return omega_plus - omega_minus

def cooperativity(g_m, kappa, gamma_m):
    """
    Cooperativity parameter: C = g_m^2 / (κ γ_m)
    
    Strong coupling requires C > 1.
    Ultra-strong coupling requires C >> 1.
    """
    return g_m**2 / (kappa * gamma_m)

def transmission_spectrum(omega_array, omega_c, omega_m, g_m, kappa, gamma_m):
    """
    Compute cavity-magnon transmission S_21(ω) including dissipation.
    
    Simplified Lorentzian model for polariton transmission.
    """
    S21 = np.zeros_like(omega_array, dtype=complex)
    
    omega_plus, omega_minus = polariton_frequencies(omega_c, omega_m, g_m)
    
    # Lower polariton linewidth
    gamma_LP = (kappa + gamma_m) / 2.0 - np.sqrt(((kappa - gamma_m)/2.0)**2 + g_m**2)
    # Upper polariton linewidth
    gamma_UP = (kappa + gamma_m) / 2.0 + np.sqrt(((kappa - gamma_m)/2.0)**2 + g_m**2)
    
    # Transmission: superposition of two Lorentzians at ω_± with linewidths γ_LP, γ_UP
    for i, omega in enumerate(omega_array):
        S21[i] = (
            1.0 / (omega - omega_minus + 1j * gamma_LP / 2.0) +
            1.0 / (omega - omega_plus + 1j * gamma_UP / 2.0)
        ) / 2.0
    
    return S21

def polariton_linewidth_narrowing(g_m, kappa, gamma_m):
    """
    Compute lower polariton linewidth narrowing factor.
    
    γ_LP ≈ √(κ γ_m) << max(κ, γ_m) for strong coupling
    """
    gamma_LP = (kappa + gamma_m) / 2.0 - np.sqrt(((kappa - gamma_m)/2.0)**2 + g_m**2)
    geometric_mean = np.sqrt(kappa * gamma_m)
    
    narrowing_factor = geometric_mean / max(kappa, gamma_m)
    
    return gamma_LP, narrowing_factor

def sensing_sensitivity(g_m, gamma_total, tau_measurement):
    """
    Compute minimum detectable magnetic moment change.
    
    Δm_min ~ (ℏ Γ) / (g_m τ)
    
    Args:
        g_m: coupling strength (rad/s)
        gamma_total: total linewidth κ + γ_m (rad/s)
        tau_measurement: measurement time (s)
    
    Returns:
        Delta_m_min: minimum detectable moment (Bohr magnetons)
    """
    HBAR = 1.0546e-34  # J·s
    mu_B = 9.2846e-24  # Bohr magneton (J/T)
    
    Delta_m = (HBAR * gamma_total) / (g_m * tau_measurement)
    
    return Delta_m / mu_B  # in units of Bohr magnetons

# Main execution
if __name__ == "__main__":
    print("="*70)
    print("CAVITY MAGNON-POLARITON SPINTRONICS")
    print("="*70)
    
    # Typical parameters (GHz scale)
    omega_c = 2 * np.pi * 10.0e9  # Cavity at 10 GHz
    omega_m = 2 * np.pi * 10.05e9  # Magnon slightly detuned
    g_m = 2 * np.pi * 10.0e6  # Coupling 10 MHz
    kappa = 2 * np.pi * 1.0e6  # Cavity decay 1 MHz
    gamma_m = 2 * np.pi * 0.5e6  # Magnon decay 0.5 MHz
    
    print(f"\nSystem Parameters:")
    print(f"  Cavity frequency ω_c/(2π) = {omega_c/(2*np.pi)/1e9:.2f} GHz")
    print(f"  Magnon frequency ω_m/(2π) = {omega_m/(2*np.pi)/1e9:.2f} GHz")
    print(f"  Coupling strength g_m/(2π) = {g_m/(2*np.pi)/1e6:.1f} MHz")
    print(f"  Cavity decay κ/(2π) = {kappa/(2*np.pi)/1e6:.1f} MHz (Q_c = {omega_c/kappa:.0f})")
    print(f"  Magnon decay γ_m/(2π) = {gamma_m/(2*np.pi)/1e6:.2f} MHz")
    
    # Compute key parameters
    C = cooperativity(g_m, kappa, gamma_m)
    vRS = vacuum_rabi_splitting(g_m, omega_c, omega_m)
    gamma_LP, narrowing = polariton_linewidth_narrowing(g_m, kappa, gamma_m)
    
    print(f"\nPolariton Properties:")
    print(f"  Cooperativity C = {C:.1f}")
    print(f"  Vacuum Rabi splitting 2g_m/(2π) = {vRS/(2*np.pi)/1e6:.1f} MHz")
    print(f"  Lower polariton linewidth γ_LP/(2π) = {gamma_LP/(2*np.pi)/1e3:.2f} kHz")
    print(f"  Linewidth narrowing factor = {narrowing:.3f}")
    print(f"  LP coherence time T₂ ≈ {1/gamma_LP*1e6:.1f} μs")
    
    # Transmission spectrum
    omega_range = np.linspace(omega_c - 5*g_m, omega_c + 5*g_m, 500)
    S21 = transmission_spectrum(omega_range, omega_c, omega_m, g_m, kappa, gamma_m)
    
    # Sensing sensitivity
    tau_meas = 1.0e-3  # 1 ms measurement
    Delta_m_min = sensing_sensitivity(g_m, kappa + gamma_m, tau_meas)
    
    print(f"\nQuantum Sensing Capabilities:")
    print(f"  Measurement time τ = {tau_meas*1e3:.1f} ms")
    print(f"  Minimum detectable moment Δm_min = {Delta_m_min:.2e} μ_B")
    
    # Plotting
    fig, axes = plt.subplots(2, 2, figsize=(13, 10))
    
    # Plot 1: Polariton eigenfrequencies vs detuning
    ax = axes[0, 0]
    detuning_range = np.linspace(-50*g_m, 50*g_m, 200)
    omega_plus_arr = []
    omega_minus_arr = []
    for delta in detuning_range:
        omega_p, omega_m_val = polariton_frequencies(omega_c, omega_m + delta, g_m)
        omega_plus_arr.append(omega_p)
        omega_minus_arr.append(omega_m_val)
    
    ax.plot(detuning_range/(2*np.pi*1e6), np.array(omega_plus_arr)/(2*np.pi*1e9), 
            'r-', linewidth=2, label='ω₊ (upper polariton)')
    ax.plot(detuning_range/(2*np.pi*1e6), np.array(omega_minus_arr)/(2*np.pi*1e9), 
            'b-', linewidth=2, label='ω₋ (lower polariton)')
    ax.axhline(y=omega_c/(2*np.pi*1e9), color='k', linestyle='--', alpha=0.5, label='Cavity bare')
    ax.set_xlabel('Detuning Δ = ω_c - ω_m (MHz)', fontsize=11)
    ax.set_ylabel('Frequency (GHz)', fontsize=11)
    ax.set_title('Polariton Avoided Crossing (Vacuum Rabi Splitting)', fontsize=12, fontweight='bold')
    ax.legend()
    ax.grid(True, alpha=0.3)
    
    # Plot 2: Transmission spectrum
    ax = axes[0, 1]
    ax.plot((omega_range - omega_c)/(2*np.pi*1e6), np.abs(S21)**2, 'b-', linewidth=2)
    ax.fill_between((omega_range - omega_c)/(2*np.pi*1e6), 0, np.abs(S21)**2, alpha=0.3)
    ax.set_xlabel('Frequency - ω_c (MHz)', fontsize=11)
    ax.set_ylabel('|S₂₁(ω)|²', fontsize=11)
    ax.set_title('Cavity-Magnon Transmission Spectrum', fontsize=12, fontweight='bold')
    ax.grid(True, alpha=0.3)
    
    # Plot 3: Cooperativity vs coupling strength
    ax = axes[1, 0]
    g_m_range = np.linspace(0.1*g_m, 10*g_m, 100)
    C_array = [cooperativity(g, kappa, gamma_m) for g in g_m_range]
    
    ax.semilogy(g_m_range/(2*np.pi*1e6), C_array, 'g-', linewidth=2.5)
    ax.axhline(y=1, color='r', linestyle='--', alpha=0.7, label='Strong coupling threshold')
    ax.axhline(y=10, color='orange', linestyle='--', alpha=0.7, label='Ultra-strong coupling')
    ax.fill_between(g_m_range/(2*np.pi*1e6), 1, 10000, alpha=0.1, color='green', label='Strong regime')
    ax.set_xlabel('Coupling Strength g_m/(2π) (MHz)', fontsize=11)
    ax.set_ylabel('Cooperativity C', fontsize=11)
    ax.set_title('Cooperativity vs Coupling', fontsize=12, fontweight='bold')
    ax.legend(fontsize=9)
    ax.grid(True, alpha=0.3, which='both')
    
    # Plot 4: Sensing sensitivity
    ax = axes[1, 1]
    tau_range = np.logspace(-6, -2, 100)  # 1 μs to 10 ms
    sensitivity = [sensing_sensitivity(g_m, kappa + gamma_m, tau) for tau in tau_range]
    
    ax.loglog(tau_range*1e3, sensitivity, 'purple', linewidth=2.5)
    ax.set_xlabel('Measurement Time (ms)', fontsize=11)
    ax.set_ylabel('Min. Detectable Moment (μ_B)', fontsize=11)
    ax.set_title('Quantum Sensing: Sensitivity vs Integration Time', fontsize=12, fontweight='bold')
    ax.grid(True, alpha=0.3, which='both')
    
    plt.tight_layout()
    plt.savefig('cavity_magnon_polariton.png', dpi=150, bbox_inches='tight')
    print(f"
[Plot saved as 'cavity_magnon_polariton.png']")
    
    print("
" + "="*70)
    print("CAVITY MAGNON-POLARITON ANALYSIS COMPLETE")
    print("="*70)

## References

1. Gütler, H., et al. (2020). "Cavity Magnon-Polaritonics." *Nature Physics*, 16(10), 1017-1024.

2. Zhang, X., et al. (2014). "Observation of Photon Pumping in a Magnon Laser." *Nature Photonics*, 12(10), 645-651.

3. Bai, L., et al. (2020). "Quantum Magnonics: The Magnon Meets the Superconducting Qubit." *Journal of Applied Physics*, 129(13), 130902.

4. Kittel, C. (1958). "On the Theory of Ferromagnetic Resonance Absorption." *Physical Review*, 73(2), 155.

5. Csaba, G., et al. (2015). "Perspectives of Using Spin Waves for Computing and Signal Processing." *Physics Letters A*, 381(17), 1471-1479.