Vortex Gyration Mode in Magnetic Disks

• Vortex gyration mode is a fundamental excitation in soft-magnetic disks where the core precesses due to the interplay of restoring and gyrotropic forces.
• It is described by Thiele’s equation, yielding eigenfrequencies in the 70 MHz to 0.7 GHz range, as validated by high-resolution imaging and simulations.
• Dipolar coupling in disk arrays produces tunable magnonic bands, enabling reconfigurable devices for logic, oscillators, and signal processing.

A vortex gyration mode is a fundamental excitation in patterned soft-magnetic elements, especially thin circular or square disks, where the magnetic ground state forms a vortex structure: the in-plane magnetization curls azimuthally, while a small region at the disk center (the "core") points out of plane. The lowest-frequency excitation of such a vortex is the translational (or gyrotropic) motion of this core about its equilibrium position. This precessional orbit is governed by the interplay of a restoring force (from exchange and magnetostatic energies) and a gyrotropic force (originating from the topological nature of the magnetization). In arrays of such disks, dipolar coupling between vortex cores leads to collective excitations—collective vortex gyration modes—which manifest as distinct magnonic bands in periodic structures. These modes are directly observable via advanced imaging and spectroscopic techniques and underpin a wide range of proposed magnonic, logic, and oscillator devices.

1. Theoretical Formulation: Single-Vortex Gyration

The gyrotropic motion of a vortex core is mathematically described by a collective coordinate (rigid-vortex) approximation, where the magnetization configuration is parameterized by the position $\mathbf{X}(t) = (X, Y)$ of the core center. Projecting the Landau–Lifshitz–Gilbert equation onto this soft mode yields Thiele's equation (neglecting damping for basic dynamics):

$$\mathbf{G} \times \dot{\mathbf{X}} + \frac{\partial U}{\partial \mathbf{X}} = 0$$

where:

• $\mathbf{G} = G\,\hat{\mathbf{z}}$ is the gyrovector, $G = 2\pi\,p\,M_s\,L/\gamma$, with $p = \pm1$ the vortex polarity, $M_s$ the saturation magnetization, $L$ the disk thickness, and $\gamma$ the gyromagnetic ratio.
• $U(\mathbf{X}) = \tfrac12\,\kappa|\mathbf{X}|^2$ is the confinement potential, with $\kappa$ the parabolic stiffness.

Solving the equation gives the gyrotropic eigenfrequency:

$$\mathbf{G} \times \dot{\mathbf{X}} + \frac{\partial U}{\partial \mathbf{X}} = 0$$0

resulting in a circular motion for the core at frequency $$\mathbf{G} \times \dot{\mathbf{X}} + \frac{\partial U}{\partial \mathbf{X}} = 0$$1, with the sense of rotation determined by $$\mathbf{G} \times \dot{\mathbf{X}} + \frac{\partial U}{\partial \mathbf{X}} = 0$$2. Representative values for single-disk eigenfrequencies are in the range 70 MHz to 0.7 GHz depending on material and geometry (Han et al., 2013).

2. Coupled-Vortex Dynamics and Collective Modes

In linear chains or arrays, each vortex core's motion is affected by both its local restoring force and dipolar interactions with neighboring vortices. For a one-dimensional (1D) chain of $$\mathbf{G} \times \dot{\mathbf{X}} + \frac{\partial U}{\partial \mathbf{X}} = 0$$3 identical disks, including nearest-neighbor dipolar coupling and damping, the equations of motion for core $$\mathbf{G} \times \dot{\mathbf{X}} + \frac{\partial U}{\partial \mathbf{X}} = 0$$4 are:

$$\mathbf{G} \times \dot{\mathbf{X}} + \frac{\partial U}{\partial \mathbf{X}} = 0$$5

with:

• $$\mathbf{G} \times \dot{\mathbf{X}} + \frac{\partial U}{\partial \mathbf{X}} = 0$$6 the core polarity and chirality of disk $$\mathbf{G} \times \dot{\mathbf{X}} + \frac{\partial U}{\partial \mathbf{X}} = 0$$7,
• $$\mathbf{G} \times \dot{\mathbf{X}} + \frac{\partial U}{\partial \mathbf{X}} = 0$$8 coupling constants along the chain ($$\mathbf{G} \times \dot{\mathbf{X}} + \frac{\partial U}{\partial \mathbf{X}} = 0$$9) and transverse ($\mathbf{G} = G\,\hat{\mathbf{z}}$0) directions,
• $\mathbf{G} = G\,\hat{\mathbf{z}}$1 the (Gilbert) damping.

Assuming plane-wave solutions, one obtains the collective mode dispersion in an infinite chain:

$\mathbf{G} = G\,\hat{\mathbf{z}}$2

where $\mathbf{G} = G\,\hat{\mathbf{z}}$3 is the 1D wavevector and $\mathbf{G} = G\,\hat{\mathbf{z}}$4 is lattice spacing (Han et al., 2013).

The sign and ordering of $\mathbf{G} = G\,\hat{\mathbf{z}}$5 and $\mathbf{G} = G\,\hat{\mathbf{z}}$6 strongly influence the dispersion curvature and bandwidth: antiparallel $\mathbf{G} = G\,\hat{\mathbf{z}}$7 doubles bandwidth compared to parallel alignment due to stronger dipolar coupling of alternating polarities.

3. Experimental Observations and Mode Analysis

Han et al. (Han et al., 2013) fabricated 5-disk chains of 2 μm diameter Permalloy disks (spacing $\mathbf{G} = G\,\hat{\mathbf{z}}$8 μm, alternating polarity $\mathbf{G} = G\,\hat{\mathbf{z}}$9, uniform chirality $G = 2\pi\,p\,M_s\,L/\gamma$0). Using time-resolved scanning transmission X-ray microscopy (STXM) with 25 nm/35 ps resolution, they directly visualized the propagation and mode structure of collective vortex gyration:

• Excitation of the leftmost disk yields propagating gyration-wave packets.
• Fourier analysis of core trajectories $G = 2\pi\,p\,M_s\,L/\gamma$1 reveals discrete eigenmodes (standing waves) matching theoretical predictions for fixed boundary conditions, with allowed $G = 2\pi\,p\,M_s\,L/\gamma$2, $G = 2\pi\,p\,M_s\,L/\gamma$3.
• Experimental mode frequencies, spatial profiles, and dispersions are in quantitative agreement with both micromagnetic simulations (OOMMF) and coupled Thiele-equation numerics.

A summary table of characteristic features:

Observable	Experimental Value	Interpretation
STXM spatial res.	25 nm	Direct imaging of core trajectories
Mode frequencies	5 peaks, $G = 2\pi\,p\,M_s\,L/\gamma$4	Discrete standing waves, $G = 2\pi\,p\,M_s\,L/\gamma$5 disks
Bandwidth control	via $G = 2\pi\,p\,M_s\,L/\gamma$6, spacing	Tunable by geometry, ordering

4. Influence of Polarization and Chirality

The collective band structure and normal mode profiles depend sensitively on the sequence of vortex polarizations $G = 2\pi\,p\,M_s\,L/\gamma$7 and in-plane chiralities $G = 2\pi\,p\,M_s\,L/\gamma$8. Consequences include:

• For uniform $G = 2\pi\,p\,M_s\,L/\gamma$9 and $p = \pm1$0, the $p = \pm1$1 band is concave-up.
• For alternating $p = \pm1$2, bands are broadened with a possible shift between $p = \pm1$3 and $p = \pm1$4 motion by $p = \pm1$5.
• The bandwidth and dispersion shape can be tailored by engineering the arrangement of $p = \pm1$6 and $p = \pm1$7.

This degree of tunability offers a magnonic analog of the band engineering used in phononic and electronic crystals, but with additional reconfigurability via external field-driven switching of core polarity.

5. Applications: Magnonic Crystals and Signal Processing

The collective vortex-gyration modes in disk arrays serve as a magnonic-crystal platform with the following attributes (Han et al., 2013):

• Reconfigurability: Each vortex core’s polarization $p = \pm1$8 can be electrically or field-set and reset, endowing dynamic functionality (e.g., switching, logic).
• High $p = \pm1$9: Mode frequencies are sharply defined and damping can be minimized using high-quality materials.
• Spectral Engineering: By selecting disk size, lattice spacing, and ordering, one can realize custom band structures, gaps, and group velocities.
• Device Integration: Arrays of such disks are amenable to on-chip integration via standard lithographic techniques.

Potential device applications include delay lines, frequency-selective filters, phase shifters, interference-based logic circuits, and elements for non-volatile information processing where logic states are encoded in mode occupancy or amplitude.

6. Quantitative Confirmation and Engineering Parameters

Key parameters governing vortex-gyration modes include:

• Restoring stiffness $M_s$0 (material, geometry dependent)
• Dipolar couplings $M_s$1 (function of disk spacing)
• Polarization and chirality sequence $M_s$2
• Lattice constant $M_s$3

These quantities all enter analytically in the $M_s$4 dispersion relation and can be experimentally tuned. Both measurements and micromagnetic modeling validate the theoretical framework, with experimental data closely tracing the analytically predicted infinite-chain dispersion.

7. Implications and Outlook

Collective vortex gyration in arrays of dipolar-coupled disks establishes a direct analogy with lattice vibration modes (phonons) in solids, but operating in the GHz regime and with explicit topological (polarity, chirality) degrees of freedom. The established tunability and nonvolatility—combined with high-$M_s$5 operation and potential for on-chip integration—pave the way for magnonic analogs of electronic devices, extending to dynamic, programmable magnonic circuits and logic (Han et al., 2013).

Research continues into optimizing material systems for lower damping, scaling arrays to 2D architectures, and exploiting the topological character of vortex states for robust information encoding and signal transport. A plausible implication is that future spintronic and magnonic computation architectures will leverage reprogrammable vortex-gyration modes as basic signal carriers.