Dynamic self-organisation and pattern formation by magnon-polarons

Fast, efficient, and non-volatile reversal of magnetisation is a central challenge in condensed matter physics with technological relevance. Large-angle, high-amplitude precessional switching differs fundamentally from conventional ferromagnetic resonance (FMR), as it is accompanied by spin-wave instabilities (Suhl instability) that funnel energy into a broad spectrum of incoherent spin waves, effectively increasing dissipation far beyond the small damping typical of FMR. As the switching is driven faster and harder, more energy is injected and redistributed into these modes, promoting chaotic magnetisation dynamics. Beyond intrinsic magnetostatic nonlinearities, coupling to the lattice provides an additional nonlinear channel: excitation of phonons can interact with magnetisation dynamics, yet its influence on large-cone-angle switching has been unclear. To address this, the authors employ multi-scale single-shot pump–probe experiments spanning nanoseconds to milliseconds. Mid-infrared pulses resonantly drive optical phonons, transiently modifying magnetic anisotropy and initiating large-amplitude precession and switching. The ensuing instabilities heavily populate low-energy magnon states; in the presence of magnetoelastic coupling this repopulation collapses into the hybridised magnon–polaron region. The phase-synchronisation and amplification of these magnon-polarons ultimately produce periodic ripple patterns and even full magnetisation reversal at the wave maxima, revealing an unexpected route from chaos to order in ultrafast magnetic switching.

Prior work established that high-power FMR leads to spin-wave instabilities (Suhl instability) that redistribute energy into a broad spectrum of magnons, causing enhanced effective damping and chaotic dynamics during large-angle precession. The Landau–Lifshitz–Gilbert framework applies, but effective dissipation can be orders of magnitude larger than in small-angle FMR due to multi-magnon processes and spin-wave generation. Additional nonlinear behavior can arise from magnetoelastic coupling: interactions with acoustic phonons can be amplified by magnetisation oscillations at comparable frequencies. YIG-based systems are noted for low magnetic and elastic damping, enabling strong magnon–phonon coupling and the formation of magnon-polarons when dispersions cross. Previous Brillouin light scattering and theory have shown that an overpopulated magnon gas can self-accumulate at magnetoelastic anti-crossings (forming hybrid magnetoelastic bosons), and parametric excitation can resonantly enhance magnons at these crossing points. However, how such coupling shapes large-angle, ultrafast magnetisation switching and pattern formation had remained unclear.

Materials: A 7.5 µm-thick doped iron-garnet film Lu1.69Y0.65Bi0.66Fe3.85Ga1.15O12 (Lu:YIG) grown on (001) GGG with trace Pb impurities was used. Material parameters include saturation magnetisation Ms ≈ 300 kA/m and weak in-plane fourfold anisotropy of 4 kA/m, yielding four in-plane easy axes and large domains at zero field.

Pump–probe imaging: Experiments were conducted at the FELIX facility. Single transform-limited mid-infrared pump micropulses (λ = 13 µm, energy ≈ 80 µJ, duration ≈ 2 ps) were focused to an elliptical spot (FWHM ≈ 300 µm × 130 µm) on the sample. Magneto-optical imaging used defocused 532 nm probe pulses of 5 ns duration (defining temporal resolution) with ≈2 µm spatial resolution. To enhance in-plane magneto-optical contrast, the sample was tilted by ≈30° relative to the probe beam. Single-shot images were recorded at variable pump–probe delays from nanoseconds to milliseconds to track spatially resolved dynamics and switching. Background subtraction isolated switched domains and early-time ripple patterns.

Time-resolved measurements: Sequential single-shot images at delays from −75 ns to 1005 ns captured pattern nucleation and propagation. Fourier transforms of selected image regions within the first 30 ns quantified ripple periodicities and were compared to the pump’s spatial spectrum.

Micromagnetic simulations: Simulations used MuMax3 with sample size 10 × 10 × 0.5 µm3 and cell size 9.76 × 9.76 × 500 nm3, including exchange interactions. Parameters: exchange A ≈ 3.7 pJ/m, Ms = 300 kA/m, biaxial anisotropy implemented with Ku = 750 J/m3 and easy axes along (1,1,0) and (1,−1,0). Gilbert damping α = 0.12, increasing quadratically up to α = 10 in outer 10% to form absorbing boundaries. Strain was introduced following a transient modification of the strain tensor to mimic nonlinear phononics-driven strain from the IR pulse. The excitation had an elliptical Gaussian spatial profile with standard deviations corresponding to FWHM 0.784 µm and 1.80 µm along (1,1,0) and (1,−1,0), respectively (aspect ratio 2.3), and temporal FDHM τ = 10 ps, peaking 50 ps after simulation start. The simulations included magnetoelastic energy to model coupling to strain and reproduced switching patterns.

Analysis and controls: The spectral dependence of switching (Supplementary Information) supports a phononic mechanism. Consideration of elastic parameters and excitation frequency ruled out Faraday-wave origins for the ripples. Estimated magnon–polaron wavelengths were computed from FMR frequency (~0.5 GHz near zero field) and acoustic velocities (TA ≈ 3.8 km/s, LA ≈ 7.2 km/s).

• A single mid-IR pump pulse transforms an initially uniform in-plane magnetisation into a four-domain, triangular switching pattern emanating diagonally from the excitation center; switching is absent at the very center. The switched domains exhibit twofold rotational symmetry.
• Early-time (first ≈30 ns) imaging reveals robust, outward-propagating quasi-periodic ripple patterns of reversed domains with visible periods on the order of ~10 µm, despite the much larger (≈300 µm) pump spot. Ripples decay within ≈300 ns, whereas large triangular domains persist up to hundreds of microseconds to milliseconds.
• Micromagnetic simulations including magnetoelastic coupling reproduce the observed four-domain switching pattern, supporting a phonon-driven anisotropy mechanism (nonlinear phononics).
• The ripple periodicity matches magnon–polaron wavelengths estimated from sample parameters: with FMR ≈ 0.5 GHz and sound velocities of ~3.8 km/s (TA) and ~7.2 km/s (LA), the inferred wavelengths are ~8–15 µm, in good agreement with the smallest observed ripple spacings.
• The results indicate that spin-wave instabilities generated by large-angle precession feed a broad magnon spectrum that, via magnetoelastic coupling, collapses and accumulates near magnon–phonon anti-crossings (magnon-polarons). These hybrid waves phase-synchronise, reaching amplitudes sufficient to cause complete magnetisation reversal at wave maxima.
• Alternative mechanisms such as Faraday waves in elastic media are excluded based on the elastic properties of YIG and the excitation conditions.
• Fourier analysis of early-time images shows spectral peaks corresponding to the ripple periodicities that are not present in the pump’s spatial spectrum, highlighting self-organisation by the magnetic–elastic system.

The study addresses how nonlinear spin-wave instabilities and magnetoelastic coupling affect large-angle, ultrafast magnetisation switching. Rather than producing only chaotic dynamics, the instability-driven broad magnon population self-organises via strong coupling to acoustic phonons, forming magnon-polarons that accumulate near dispersion anti-crossings. This leads to emergent spatial order: quasi-periodic ripple domains and a robust four-domain triangular switching pattern. Phase synchronisation of magnon-polaron wave packets explains how sub-spot-size periodic reversal can arise, with wavelengths set by the hybridised magnon–phonon modes, not by the optical excitation profile. These findings reveal an alternative, coherent mechanism for magnetisation reversal mediated by short-wavelength magnetoelastic waves and suggest that coupling to the lattice can funnel and organise magnetic energy into specific modes with enhanced group velocity, enabling efficient angular momentum transfer. The work thus reframes ultrafast switching in low-damping garnets as a route from chaos to order through hybrid quasiparticles, with broader relevance to synchronisation phenomena observed across physical and biological systems.

The work demonstrates that large-amplitude, ultrafast excitation of a low-damping magnetic garnet leads to self-organisation of magnetisation into periodic ripple patterns and macroscopic triangular domains via the formation and phase synchronisation of magnon-polarons. Simulations and spectral analyses support a phonon-driven mechanism wherein a broad magnon spectrum collapses to magnon–phonon anti-crossings, yielding coherent, high-amplitude magnetoelastic waves capable of inducing complete magnetisation reversal at their maxima. These insights provide an alternative pathway for fast, energy-efficient magnetisation control using coherent, short-wavelength magnetoelastic excitations. Future research should resolve the sub-nanosecond spatiotemporal development of magnon-polarons and synchronisation dynamics to deepen understanding of phonon-induced switching and general wave entrainment phenomena.

The time resolution of the single-shot magneto-optical imaging is limited by the 5 ns probe pulse duration, preventing direct observation of sub-nanosecond formation and phase-locking of magnon-polarons; the authors explicitly note the need for sub-ns studies. Spatial resolution (~2 µm) and the defocused imaging approach may under-resolve the finest features. The analysis infers magnon–polaron condensation from imaging, simulations, and spectral estimates rather than direct frequency-resolved measurements of hybrid modes within the first nanoseconds.