Universal theory of domain-wall width in multi-sublattice Heisenberg magnets

This paper proposes a universal expression for the domain-wall width in multi-sublattice Heisenberg magnets by establishing an exact connection between the domain-wall profile and long-wavelength spin-wave dispersion, a framework that accurately predicts widths across various magnetic orders and lattice structures while providing a microscopic foundation for their temperature dependence.

The Gist

Imagine a magnet not as a solid, uniform block, but as a vast crowd of tiny, spinning tops (atoms) all trying to point in the same direction. Sometimes, this crowd splits into two groups: one group pointing "up" and another pointing "down." The invisible line where these two groups meet is called a domain wall.

Think of a domain wall like a transition zone or a "ramp" on a highway. On one side, all cars (spins) are driving North; on the other, they are driving South. The domain wall is the curved section of the road where the cars gently turn around. The width of this wall is simply how many cars it takes to make that turn.

The Problem: A One-Size-Fits-All Rule That Broke

For simple magnets (like a standard fridge magnet), scientists had a perfect, simple recipe to calculate how wide this turn would be. It was like a rule saying: "The width depends on how tightly the cars hold hands (exchange) versus how strongly they want to stay in their lanes (anisotropy)."

However, the real world is messy. Many advanced magnets are made of multiple subgroups (sublattices) of atoms that interact in complex ways. Some might be heavy, some light; some might pull, others push. In these complex "multi-sublattice" magnets, the old simple rule stopped working. Scientists didn't have a universal way to predict the width of the turn in these complicated crowds.

The Solution: A Universal "Traffic Map"

The authors of this paper propose a universal formula that works for any type of magnetic order—whether it's a simple crowd, a split crowd (ferromagnet), a fighting crowd (antiferromagnet), or a mixed crowd (ferrimagnet).

Here is the core idea using an analogy:

The "Spin-Wave" Analogy:
Imagine the magnetic atoms are dancers.

• Spin Waves: If you nudge the dancers slightly, they ripple through the crowd like a wave. These ripples are called "spin waves."
• The Domain Wall: A domain wall is like a giant, static ripple frozen in place.

The paper's big discovery is that you can predict the size of the frozen ripple (the wall) by studying the tiny ripples (the waves).

The authors found that if you look at the "energy map" of how these tiny waves move (specifically, how fast they move and how much energy it takes to start them), you can mathematically calculate the width of the domain wall.

How They Proved It

They didn't just guess; they built a massive digital simulation of these atomic crowds. They tested their new formula on:

1. Rock-salt magnets: Complex 3D structures with two types of atoms.
2. Honeycomb magnets: Flat, 2D structures (like graphene) that look like a beehive.
3. Kagome magnets: Flat structures with a pattern of triangles and stars.

In every single case, from simple to highly complex, their new "universal formula" matched the computer simulations perfectly. It worked whether the temperature was near absolute zero or getting close to the point where the magnetism disappears.

The "Temperature" Twist

The paper also explains what happens when you heat things up.

• Cold: The atoms are stiff and hold their positions tightly. The formula works easily.
• Hot: The atoms start shaking and dancing wildly. This changes the "rules" of how they hold hands.
• The Fix: The authors showed that their formula can be "renormalized" (adjusted) to account for this shaking. By measuring how the tiny waves change as the temperature rises, the formula can still accurately predict how the domain wall width changes, all the way up to the point where the magnet stops working.

The Takeaway

In simple terms, this paper provides a master key for understanding magnetic walls. Before, scientists needed a different key for every different type of complex magnet. Now, they have one universal key that works for all of them, based on the simple idea that the shape of a frozen wave (the wall) is determined by the behavior of the tiny ripples (the spin waves).

This allows scientists to predict the behavior of complex magnetic materials without needing to simulate every single atom every time, bridging the gap between the tiny atomic world and the larger devices we might use in the future.

Technical Summary

Technical Summary: Universal Theory of Domain-Wall Width in Multi-Sublattice Heisenberg Magnets

Problem Statement
Magnetic domain walls (DWs) are fundamental topological textures in condensed matter physics, governing phenomena ranging from chirality-dependent transport to ultrafast angular-momentum flow. In spintronic applications, such as racetrack memory, the DW width ($\delta_{DW}$) is a critical parameter that determines energy, stability, mobility, and dynamic response. While the width in simple ferro- and antiferromagnets is well-described by the compact expression $\delta_{DW} = \sqrt{A/K}$ (where $A$ is exchange stiffness and $K$ is uniaxial anisotropy), this description fails in complex multi-sublattice systems. In such systems, including ferrimagnets, antiferromagnets, and low-dimensional magnets, the presence of multiple spin-wave (SW) branches, sublattice-dependent couplings, and mode hybridization obscures the identification of effective parameters linking the SW spectrum to the DW width. Consequently, a general, universal relation connecting measurable SW properties to DW width in multi-sublattice Heisenberg magnets has been lacking.

Methodology
The authors propose a theoretical framework establishing an exact connection between the long-wavelength spin-wave dispersion and the domain-wall profile in generic multi-sublattice Heisenberg magnets with uniaxial anisotropy. The core of the approach involves:

1. Theoretical Derivation: The authors derive a universal expression for $\delta_{DW}$ based on the approximation of the spin-wave spectrum as $\varepsilon^\beta_q = \varepsilon^\beta_0 + A_\beta q^2$, where $\beta$ indexes the SW branches. The proposed universal formula is:
   $$\delta_{DW} = \sqrt{\sum_\beta \frac{A_\beta}{\varepsilon^\beta_0}}$$
   where $A_\beta$ is the SW stiffness and $\varepsilon^\beta_0$ is the energy gap of the $\beta$-th mode.

2. Atomistic Spin Dynamics (ASD) Simulations: To validate the theory, the authors perform large-scale ASD simulations on various lattice structures. These simulations include:

   • Rock-salt type magnets: Covering ferrimagnetic (FI), ferromagnetic (FM), and antiferromagnetic (AFM) orderings.
   • Two-dimensional (2D) magnets: Specifically honeycomb and kagome lattices.
   • Temperature Effects: The simulations incorporate thermal fluctuations up to the critical temperature ($T_c$) to account for the thermal renormalization of spin-model parameters.

3. Validation Protocol: The DW width extracted from ASD simulations is obtained by fitting the sublattice-resolved magnetization profiles to a hyperbolic tangent function. These empirical widths are compared against the values predicted by the universal relation using parameters derived from the simulated SW spectra.

Key Contributions and Results

• Universal Expression: The paper establishes that the DW width in multi-sublattice systems is the square root of the sum of the ratios of SW stiffness to energy gap for all active branches ($\delta_{DW} = \sqrt{\sum A_\beta/\varepsilon^\beta_0}$). This unifies the description of DWs across ferro-, antiferro-, and ferrimagnetic orders.
• Ferrimagnetic Generalization: For a two-sublattice ferrimagnet, the authors derive a simplified expression in the exchange-dominated regime. They demonstrate that the standard ferromagnetic relation is a special case of their universal formula, with effective exchange stiffness and anisotropy density emerging from the contributions of both high-energy and low-energy SW branches.
• Temperature Dependence: The framework successfully describes the temperature dependence of $\delta_{DW}$ up to the magnetic ordering transition. By utilizing thermally renormalized parameters from ASD simulations, the theory captures the evolution of the DW width as SW gaps and stiffnesses change with temperature. Notably, near the compensation temperature in ferrimagnets, the contribution from the high-energy branch becomes comparable to the low-energy branch, altering the width dynamics.
• Validation Across Lattices:
  • Rock-Salt Structures: Excellent quantitative agreement is found between the universal relation and ASD simulations for FM, AFM, and FI systems. In the AFM case, the necessity of summing contributions from two degenerate branches is highlighted.
  • 2D Honeycomb Lattice: The relation holds for a 2D honeycomb ferromagnet, accurately predicting the width across a range of exchange interactions ($J_1, J_2, J_3$), including regimes where the system approaches instability.
  • 2D Kagome Lattice: For a three-sublattice kagome ferromagnet, the theory correctly accounts for three SW branches, noting that flat bands (where stiffness $A \to 0$) do not contribute to the width.

Significance
The paper claims to provide a "universal expression" and a "unified framework" for describing magnetic textures across distinct ordering types. By establishing a direct link between microscopic spin-wave properties and the continuum description of domain walls, the work offers a powerful tool for understanding DW profiles in complex spin systems. The approach bridges the gap between atomistic spin interactions and continuum micromagnetics, enabling the prediction of DW widths in multi-sublattice magnets without relying on ad-hoc parameter fitting. Furthermore, the establishment of a general microscopic foundation for temperature dependence allows for a self-consistent description of DWs from zero temperature to the ordering transition, facilitating the engineering of magnetic textures in advanced spintronic materials.