Spin Waves &amp; Magnons in a Ferromagnet

3D Simulations

Article

Electromagnetism · Condensed Matter · ⏱ ~12 min read · Last updated: 9 July 2026

Spin Waves & Magnons — Collective Excitations in a Ferromagnet

Below its Curie temperature (the threshold above which heat disrupts magnetic order entirely), a ferromagnet's spins align, but not perfectly and not statically — the lowest-energy excitations above the fully aligned ground state are not single flipped spins but collective, wave-like precessions that ripple coherently through the whole lattice. Felix Bloch quantised these spin waves in 1930, calling the quantum of excitation a magnon: a boson carrying one unit of spin-flip and a well-defined momentum and energy, in direct analogy to phonons for lattice vibrations and photons for the electromagnetic field. Magnons now underpin an entire subfield, magnonics, aiming to carry and process information with spin currents instead of charge currents.

TL;DR: A ferromagnet's spins support collective wave-like ripples called spin waves, quantised as magnons — bosons with definite energy and momentum, analogous to phonons. Their dispersion is quadratic, giving Bloch's T^3/2 law for magnetisation loss with temperature, while antiferromagnetic magnons disperse linearly instead, following a T² law.

Written by MySimulator Team · Reviewed by MySimulator Editorial Review

1. The Heisenberg Exchange Model

Localised spins on a lattice, nearest-neighbour exchange interaction: H = -J Σ_<ij> S_i \cdot S_j - gμ_B B Σ_i S_i^z J > 0 \to ferromagnetic coupling (parallel spins favoured, ground state: all spins aligned along +z) J < 0 \to antiferromagnetic coupling (antiparallel neighbours favoured) Exchange energy J arises from the Pauli exclusion principle plus Coulomb repulsion (NOT classical dipole-dipole interaction, which is far too weak to produce room-temperature ferromagnetism) — a purely quantum- mechanical effect. Typical exchange energies: J/k_B ~ 100s of Kelvin for 3d transition-metal ferromagnets (Fe, Co, Ni), setting the Curie temperature scale (T_C ~ zJS(S+1)/3k_B in mean-field theory, z = coordination number).

2. Classical Spin Precession Picture

Semiclassical picture: think of each spin as a classical magnetic moment vector precessing about its local effective field. Landau- Lifshitz equation of motion: dS_i/dt = -(γ) S_i \times B_eff,i B_eff,i = -(1/gμ_B) \partialH/\partialS_i = (J/gμ_B) Σ_j S_j + B (external + exchange field from neighbours) In the ground state all spins point along z and are static. A spin wave is a small transverse deviation δS_i that precesses with phase advancing smoothly from site to site: S_i^x, S_i^y \propto cos(k\cdotr_i - ωt), forming a coherent wave of precession travelling through the lattice — no single spin fully flips, each just tips slightly and the tip direction rotates around the z-axis as the wave passes.

3. Magnon Dispersion Relation

1D ferromagnetic chain, lattice spacing a, spin S, nearest- neighbour exchange J: Exact dispersion: ℏω(k) = 4JS[1 - cos(ka)] Long-wavelength limit (ka ≪ 1): ℏω(k) \approx 2JSa²k² — QUADRATIC dispersion (unlike phonons, which are LINEAR ω \propto k at small k). This quadratic dispersion is the direct signature of exchange- dominated spin waves and comes from the underlying Heisenberg Hamiltonian's rotational (not translational) symmetry breaking. 3D simple-cubic lattice: ℏω(k) = 2JS[3 - cos(k_x a) - cos(k_y a) - cos(k_z a)] Magnon group velocity: v_g = dω/dk = 4JSa²k/ℏ (long- wavelength limit) — vanishes at k=0 (uniform precession/ ferromagnetic resonance mode) and grows linearly with k for small k, unlike the constant group velocity of acoustic phonons.

4. Bloch's T^3/2 Law

At finite temperature T, thermally excited magnons reduce the net magnetisation from its T=0 saturation value M(0). Number of thermally excited magnons per site (Bose-Einstein statistics, quadratic dispersion, 3D): n_magnon(T) \propto T^(3/2) Bloch's T^3/2 law: [M(0) - M(T)] / M(0) = A \cdot T^(3/2) A depends on exchange stiffness D = 2JSa² and lattice geometry; A \propto [k_B/(4πD)]^(3/2) \cdot ζ(3/2) (ζ = Riemann zeta function, from the Bose-Einstein integral) Physical origin: quadratic dispersion \to magnon density of states g(ε) \propto \sqrtε at low energy (same as free non-relativistic particles in 3D) \to thermal population integral yields the characteristic 3/2 power, distinct from the T³ Debye law for phonons (linear dispersion, different density of states). Experimentally verified in Fe, Ni, and many ferromagnetic insulators (e.g. YIG) at low temperature.

5. Magnon Quantisation (Holstein-Primakoff)

The Holstein-Primakoff transformation maps spin operators to boson creation/annihilation operators a†, a (exact for a single site): S_i^z = S - a_i†a_i S_i^+ = \sqrt(2S) \sqrt(1 - a_i†a_i/2S) a_i \approx \sqrt(2S) a_i (large-S / low-excitation approximation) S_i^- = (S_i^+)† Substituting into the Heisenberg Hamiltonian and Fourier-transforming to momentum space (a_k = N^(-1/2) Σ_i e^(-ik\cdotr_i) a_i) diagonalises H into a sum of independent harmonic oscillators: H \approx E_0 + Σ_k ℏω(k) \cdot a_k†a_k Each mode k is now literally a quantum harmonic oscillator — a magnon is one quantum of that oscillator, a boson with energy ℏω(k), momentum ℏk, and spin angular momentum ℏ (it reduces total S_z by exactly one unit per magnon created). Magnons obey Bose-Einstein statistics, can Bose-condense (magnon BEC observed at room temperature in YIG under microwave pumping, Demokritov et al. 2006), and interact weakly at low density (magnon-magnon scattering appears only at higher order in the 1/2S expansion).

6. Antiferromagnetic Magnons

Antiferromagnet: two interpenetrating sublattices A, B with opposite spin orientation (Néel state), J < 0 exchange. Because neighbouring spins are already antiparallel, spin-wave theory needs a Holstein-Primakoff transformation on EACH sublattice separately (with an extra sign flip on B), followed by a Bogoliubov transformation to diagonalise the resulting Hamiltonian (which mixes a_k and b_{−k} creation/annihilation operators — the ground state is a squeezed vacuum, not the naive Néel state). Resulting dispersion is LINEAR at small k (unlike the ferromagnetic quadratic law): ℏω(k) ≈ 2|J|Sza · |k| (z = coordination number) — antiferromagnetic magnons behave like relativistic massless particles or acoustic phonons at long wavelength. Consequence: the antiferromagnetic analogue of Bloch's law is T² (not T^3/2), matching the linear-dispersion density of states, just as acoustic phonons give the Debye T³ specific heat.

7. JavaScript Spin-Chain Simulator

// Classical Landau-Lifshitz integration of a 1D Heisenberg spin chain
// (small-angle spin-wave regime, periodic boundary conditions)

function createChain(N, S = 1) {
  // spins start aligned along z, with a small transverse perturbation
  const spins = [];
  for (let i = 0; i < N; i++) {
    spins.push({ x: 0.05 * Math.sin((2 * Math.PI * i) / N), y: 0, z: S });
  }
  return spins;
}

function effectiveField(spins, i, J, Bz) {
  const N = spins.length;
  const left = spins[(i - 1 + N) % N];
  const right = spins[(i + 1) % N];
  return {
    x: J * (left.x + right.x),
    y: J * (left.y + right.y),
    z: J * (left.z + right.z) + Bz,
  };
}

function cross(a, b) {
  return {
    x: a.y * b.z - a.z * b.y,
    y: a.z * b.x - a.x * b.z,
    z: a.x * b.y - a.y * b.x,
  };
}

// Advance the chain by one Landau-Lifshitz Euler step (gamma absorbed into dt)
function step(spins, J, Bz, dt) {
  const N = spins.length;
  const torques = spins.map((s, i) => {
    const Beff = effectiveField(spins, i, J, Bz);
    const t = cross(s, Beff);
    return { x: -t.x, y: -t.y, z: -t.z };
  });
  return spins.map((s, i) => ({
    x: s.x + dt * torques[i].x,
    y: s.y + dt * torques[i].y,
    z: s.z + dt * torques[i].z,
  }));
}

// Predicted magnon dispersion (long-wavelength limit) for comparison
function magnonFrequency(k, J, S, a = 1) {
  return 4 * J * S * (1 - Math.cos(k * a)); // hbar = 1 units
}

// Example: N=64 chain, J=1, Bz=0, integrate 500 steps
let chain = createChain(64, 1);
for (let t = 0; t < 500; t++) chain = step(chain, 1, 0, 0.01);
console.log('Predicted omega(k=2pi/64):', magnonFrequency((2 * Math.PI) / 64, 1, 1));

8. Applications

Magnonics & Spin-Wave Logic

Spin waves carry information without moving electrons, so magnonic logic gates and interferometers can process data with far lower Joule heating than CMOS charge-based circuits.

YIG-Based Microwave Devices

Yttrium iron garnet (YIG) has extremely low magnon damping, making it the standard material for tunable microwave filters, oscillators, and magnon Bose-Einstein condensate experiments.

Spintronics & Spin Pumping

Magnons injected at a ferromagnet/normal-metal interface convert to a pure spin current via the inverse spin Hall effect — a key mechanism for spin-based memory and logic devices.

Magnon-Photon & Magnon-Phonon Coupling

Strong coupling between magnons and microwave photons in cavities, or magnons and acoustic phonons in magnetostrictive films, enables hybrid quantum systems for coherent information transfer.

💡 Open Maxwell's Equations Simulation →