Dispersion and lifetimes of magnons in non-collinear magnets from time dependent density functional theory

Using a novel first-principles approach based on non-collinear KKR Green's functions within time-dependent density functional theory, this study investigates the spin dynamics of Mn$_3$Rh, revealing three linearly dispersing Goldstone modes and characterizing their substantial Landau damping away from the Brillouin zone center.

The Gist

Imagine a dance floor where the dancers aren't just moving in a straight line or spinning in unison. Instead, they are arranged in a triangular pattern, and each dancer is facing a different direction, creating a complex, swirling, non-collinear formation. This is what happens inside a special type of magnetic material called Mn₃Rh, which the scientists in this paper studied.

Here is a breakdown of what the paper is about, using simple analogies:

1. The Setting: A Frustrated Dance Floor

Most magnets are like a crowd of people all facing North (ferromagnets) or alternating North-South (antiferromagnets). But Mn₃Rh is a "frustrated" magnet. Imagine three friends standing in a triangle, each trying to face away from the others. They can't all be perfectly opposite at the same time. This creates a unique, swirling magnetic pattern called a kagome lattice (named after a Japanese woven basket pattern).

The scientists wanted to know: If you poke this magnetic dance floor, how do the waves travel through it, and how long do they last?

2. The Problem: The Old Maps Don't Work

To understand these magnetic waves (called magnons), physicists usually use a simple map called the "Heisenberg model." Think of this like a map of a flat, two-dimensional city. It works great for simple magnets.

However, for these complex, swirling magnets, the old map is like trying to navigate a 3D rollercoaster with a flat paper map. It misses crucial details, specifically:

• The "Landau Damping": In simple magnets, waves travel forever. In these complex metals, the waves crash into the "traffic" of electrons, losing energy and dying out quickly. The old map ignores this traffic.
• The "Polarization": In simple magnets, the dancers just spin up and down. In these complex magnets, the dancers spin in weird, tilted, 3D directions.

3. The New Tool: A High-Definition 4D Camera

The authors developed a brand-new, super-complex computer simulation based on Time-Dependent Density Functional Theory (TD-DFT).

• The Analogy: If the old method was a sketch, this new method is a high-definition, 4D movie that captures every electron, every spin, and every moment in time.
• The Challenge: Doing this math is like trying to solve a puzzle where every piece has 256 different sides. The team had to write special computer code (using "symbolic algebra") to handle the massive amount of data without crashing.

4. The Discovery: Three Waves and a Surprise

When they ran their simulation on Mn₃Rh, they found three main things:

A. Three Distinct Waves (The Goldstone Modes)
Just as a drumhead can vibrate in different ways, this magnetic triangle vibrates in three distinct patterns.

• In the long-wavelength limit (gentle ripples), these waves travel at a constant speed (linear dispersion).
• The Twist: These waves aren't just spinning up and down. They have "non-trivial polarizations," meaning the magnetic moments (the dancers) are doing a complex, coupled precession (wobbling in a specific 3D pattern) that changes depending on where you are on the dance floor.

B. The "Landau Damping" (The Energy Leak)
This is the big discovery. The waves don't just fade away; they get "eaten" by the electrons.

• The Analogy: Imagine a surfer (the magnon) riding a wave. In a calm ocean, they ride forever. In a busy harbor full of motorboats (electrons), the surfer crashes into the wake of the boats and loses energy.
• The Finding: The waves travel fine in the center of the "dance floor" (the Brillouin zone), but as they get to the edges, they crash harder into the electron traffic and die out much faster.

C. The Polarization Surprise (The "Who Gets Hit" Factor)
This is the most fascinating part. The scientists found that how long a wave lasts depends entirely on how it is spinning.

• The Analogy: Imagine two surfers riding waves of the exact same speed and size. One is wearing a bright red wetsuit, the other blue.
  • The Red Surfer (a specific spin pattern) gets hit by the motorboats constantly and crashes immediately.
  • The Blue Surfer (a different spin pattern) somehow slips through the traffic and rides much longer.
• Why? It turns out that the "Red Surfer" is localized on a specific atom (Mn2) where the electron traffic is heavy and chaotic. The "Blue Surfer" is on a different atom (Mn3) where the electron traffic is smoother.

5. Why Does This Matter?

This isn't just about math; it's about the future of computers.

• Spintronics: We are moving from computers that use electric charge to computers that use "spin" (magnetism).
• The Takeaway: If we want to build computers using these materials, we can't just treat all magnetic waves the same. We have to engineer them carefully. By choosing the right "dance move" (polarization), we can make the magnetic signals last longer and travel further, or make them stop quickly where we want them to.

Summary

The paper is like a detective story where the scientists built a super-computer microscope to watch a complex magnetic dance. They discovered that the dancers (spins) move in three unique ways, and that some dance moves are much more durable than others because they avoid crashing into the electron traffic. This gives us a new blueprint for designing faster, more efficient magnetic computers.

Technical Summary

Here is a detailed technical summary of the paper "Dispersion and lifetimes of magnons in non-collinear magnets from time dependent density functional theory" by Eilmsteiner, Ernst, and Buczek.

1. Problem Statement

The paper addresses the challenge of accurately describing spin dynamics in non-collinear (NC) magnetic systems, specifically triangular kagome antiferromagnets like Mn$_3$Rh.

• Limitations of Existing Methods: Traditional approaches, such as the Heisenberg model mapping or adiabatic approximations, fail to capture essential physics in NC systems. Specifically, they cannot account for Landau damping (the decay of spin waves into electron-hole pairs), which is a dominant attenuation mechanism in metallic magnets.
• Complexity of NC Systems: In NC magnets, the single-electron spin is not a good quantum number, and spin and charge degrees of freedom are coupled. This makes the separation of transverse and longitudinal channels (common in collinear magnets) non-trivial.
• Gap in Theory: While Linear Response Time-Dependent Density Functional Theory (LRTDDFT) is the state-of-the-art for capturing Landau damping, its application to complex NC systems has been hindered by algorithmic complexity and the lack of a formalism that treats spin and charge dynamics on an equal footing without adiabatic approximations.

2. Methodology

The authors developed a novel first-principles approach based on LRTDDFT combined with the Korringa-Kohn-Rostoker (KKR) Green's function method.

• Formalism:
  • The true dynamical susceptibility $\chi_{ij}(q, r, r', \omega)$ is calculated by solving the Dyson equation: $\chi = \chi_{KS} + \chi_{KS}(v_C + K_{xc})\chi$.
  • Coupling: Unlike collinear systems, the charge and spin responses are not decoupled. The Coulomb interaction ($v_C$) is retained, and the exchange-correlation kernel ($K_{xc}$) is treated within the Adiabatic Local Spin Density Approximation (ALSDA).
  • Green's Functions: The non-interacting susceptibility $\chi_{KS}$ is computed as a product of two Kohn-Sham Green's functions.
• Computational Innovation:
  • Trace Evaluation: A major hurdle in NC systems is the evaluation of the trace over spin indices in the Green's function product. For KKR, this involves up to 256 matrix elements. The authors overcame this by using automatic FORTRAN code generation based on computational symbolic algebra.
  • Symmetry Handling: The method ensures the correct emergence of Nambu-Goldstone Bosons (NGBs). In triangular antiferromagnets (TAFs), the spontaneous breaking of continuous symmetries (related to a rigid rotor) should yield three linearly dispersing Goldstone modes. The authors implemented strict convergence controls to ensure their numerical scheme respects this symmetry.
• Analysis Tools:
  • Loss Matrix: Collective modes (magnons) are identified as eigenvectors of the anti-Hermitian part of the susceptibility (the loss matrix $L(\omega)$).
  • Landau Maps: To analyze damping origins, the authors constructed "Landau maps" visualizing the spectral density of initial electronic states ($k$) that form Stoner pairs (transitions to $k+q$) contributing to magnon decay.

3. Key Contributions

• First Application to NC Magnets: This is the first application of LRTDDFT to non-collinear magnetic systems, extending the theory beyond collinear limits.
• Unified Treatment: The approach treats spin and charge dynamics simultaneously, avoiding the ambiguities of mapping to effective Heisenberg models.
• Landau Damping in NC Systems: The paper successfully generalizes the concept of Stoner excitations to NC magnets, revealing how specific electronic band structures and sublattice localizations dictate damping rates.
• Polarization-Dependent Lifetimes: The study demonstrates that magnon lifetimes are not solely determined by energy and momentum but are critically dependent on the polarization (spatial form) of the mode.

4. Results

The study focuses on Mn$_3$Rh, a kagome triangular antiferromagnet.

• Electronic Structure: Despite vanishing net magnetization, Mn$_3$Rh exhibits an intrinsically spin-polarized band structure (similar to altermagnets). The Fermi surface features strongly dispersing, weakly spin-polarized bands crucial for Landau damping.
• Magnon Dispersion:
  • Three distinct Goldstone modes are observed, dispersing linearly in the long-wavelength limit ($q \to 0$), consistent with the symmetry breaking of TAFs.
  • At the Brillouin zone (BZ) boundaries (points M and M'), the modes reach energies of ~330 meV.
  • The dispersion agrees qualitatively with previous adiabatic results but includes the crucial non-adiabatic damping effects.
• Mode Polarization and Sublattice Localization:
  • The magnon modes exhibit non-trivial spatial forms. At high-symmetry points (M, M'), the modes split into $\delta$-modes and $\sigma$-modes.
  • Localization: The $\sigma$-mode at point M is localized primarily on the Mn2 sublattice, while $\delta$-modes involve Mn1 and Mn3. This separation is linked to the specific chain of atoms the magnon propagates through.
• Landau Damping and Lifetimes:
  • Energy Dependence: Damping increases with magnon energy, reaching a Full Width at Half Maximum (FWHM) of ~100 meV.
  • Polarization Dependence (Key Finding): Modes with similar energy and momentum but different polarizations exhibit drastically different lifetimes.
    • At intermediate energies (~100 meV), the lifetime difference between modes exceeds a factor of 4.
    • At points M and M', the $\delta$-modes are damped ~50% more strongly than the $\sigma$-mode.
  • Mechanism: The difference arises from Landau hotspots.
    • The $\delta$-modes excite specific electronic transitions that create "hotspots" in the Landau maps (high density of available Stoner pairs).
    • The $\sigma$-mode (localized on Mn2) experiences less damping because the electronic bands associated with the Mn2 site are less polarized along the $\Gamma$-M direction, resulting in fewer available decay channels compared to Mn1 and Mn3.

5. Significance

• Fundamental Physics: The work provides a rigorous, first-principles understanding of how non-collinearity alters spin-wave decay, moving beyond simplified models to capture the interplay between lattice symmetry, electronic structure, and spin dynamics.
• Magnonic Engineering: The discovery that magnon lifetimes can be tuned by polarization and sublattice localization opens new avenues for "magnonic engineering." By designing materials or structures that favor specific polarizations, one could potentially control the propagation distance of spin waves at the single-atom scale.
• Methodological Advancement: The developed framework sets a new standard for calculating spin excitations in complex, frustrated, and non-collinear magnetic materials, which are central to next-generation spintronics and altermagnetism research.