🔬 mesoscale physics

Demagnetization in micromagnetics: magnetostatic self-interactions of bulk chiral magnetic skyrmions

This paper establishes a theoretical and numerical framework demonstrating that long-range magnetostatic dipolar interactions break the energy degeneracy of bulk chiral skyrmions, specifically stabilizing Heusler antiskyrmions into square-lattice crystals while leaving Bloch skyrmions unaffected and slightly shrinking Néel skyrmions.

Original authors: Paul Leask, Martin Speight

Published 2026-02-05

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Paul Leask, Martin Speight

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a giant, invisible ocean made of tiny, spinning tops. In physics, these spinning tops are called spins, and when they line up in a specific, swirling pattern, they create a magnetic shape known as a skyrmion. You can think of a skyrmion as a tiny, stable whirlpool or a magnetic tornado that can move through a material without falling apart.

For a long time, scientists knew that these whirlpools could form in certain materials, but they were mostly looking at them in thin, flat sheets (like a piece of paper). This new paper asks a bigger question: What happens to these whirlpools if we look at a thick, 3D block of material instead?

Here is the story of what the researchers found, broken down into simple concepts:

1. The Three Types of Whirlpools

In the world of these magnetic whirlpools, there are three main "personality types," determined by how the spins twist. The paper calls them:

The Bloch Skyrmion: A smooth, classic swirl.
The Néel Skyrmion: A slightly different twist, like a spiral staircase.
The Antiskyrmion: A more complex, "anti-twist" shape found in special materials called Heusler compounds.

In a perfect, simplified world (without considering the material's own magnetic pull), these three types are essentially twins. They have the exact same energy and would behave identically. It's like having three identical twins who all weigh the same and run at the same speed.

2. The New Ingredient: The "Self-Gravity"

The researchers added a new factor to their simulation: Magnetostatic Self-Interaction.

Think of this as the material's own "self-gravity" or "self-awareness." Every tiny spinning top in the material creates a tiny magnetic field that pushes or pulls on its neighbors. In thin films, this effect is often ignored or treated simply. But in a thick, 3D block, these tiny magnetic forces add up and create a complex "demagnetizing field."

The researchers wanted to see how this "self-gravity" changes the behavior of our three whirlpool twins.

3. The Results: Twins Diverge

When they turned on this "self-gravity," the three identical twins suddenly became very different:

The Bloch Skyrmion (The Classic Swirl): It didn't care at all. The self-gravity had zero effect on it. It stayed exactly the same size and shape. It's like a swimmer in a calm pool who doesn't notice the water moving around them.
The Néel Skyrmion (The Spiral): It got a little bit smaller. The self-gravity squeezed it slightly, making it more compact.
The Antiskyrmion (The Complex Twist): This one had the most dramatic reaction.
- Shape Shift: It lost its perfect circular symmetry. Instead of being a round whirlpool, it squashed into a square shape.
- The Crystal Effect: In the past, scientists thought these whirlpools would just repel each other and drift infinitely far apart. But with the self-gravity turned on, the Antiskyrmions started to attract each other. They didn't just float apart; they grabbed hands and formed a neat, organized square crystal lattice (like a grid of squares).

4. Why This Matters (According to the Paper)

The paper claims that by including this "self-gravity" (the dipole-dipole interaction), they discovered a way to stabilize these magnetic crystals in 3D bulk materials.

Specifically, they found that Heusler antiskyrmions (the complex twist type) naturally want to form a square crystal structure in a 3D block of material, whereas the other types prefer to stay far apart.

The Big Picture Analogy

Imagine you have three types of dancers in a large ballroom:

Dancer A spins in a perfect circle.
Dancer B spins in a spiral.
Dancer C does a complex, twisting move.

If the room is empty, they all dance the same way. But then, imagine the room fills with a thick, sticky gel (the "self-gravity").

Dancer A doesn't feel the gel and keeps spinning perfectly.
Dancer B gets slightly squished by the gel and spins tighter.
Dancer C gets so affected by the gel that they stop spinning alone and start linking up with other Dancers C to form a rigid, square grid, because the gel makes it energetically favorable for them to stick together.

Summary

The paper provides a new mathematical and computer-based framework to study these 3D magnetic whirlpools. Their main discovery is that the material's own magnetic "self-pull" breaks the symmetry between the different types of skyrmions. Most importantly, it reveals that antiskyrmions can naturally form stable, square-shaped crystals in 3D materials, a phenomenon that wouldn't be predicted if you ignored the material's internal magnetic interactions.

Technical Summary: Demagnetization in Micromagnetics of Bulk Chiral Magnetic Skyrmions

Problem Statement
The paper addresses the theoretical and numerical modeling of three-dimensional bulk chiral ferromagnets, specifically focusing on the role of long-range magnetostatic dipole-dipole interactions (DDI) in stabilizing topological spin textures. While previous studies have examined DDI in thin films or specific contexts (e.g., stripe domains or biskyrmions), there is a lack of comprehensive analysis regarding how DDI affects the degeneracy and stability of different skyrmion types in a translationally invariant bulk system.

In the absence of DDI, three distinct Dzyaloshinskii–Moriya interaction (DMI) terms—Dresselhaus, Rashba, and Heusler—are gauge-equivalent. They yield degenerate energy states for their corresponding skyrmion textures: Bloch (Dresselhaus), Néel (Rashba), and antiskyrmion (Heusler). The central problem is to determine how the inclusion of the non-local DDI and its back-reaction on the magnetization breaks this degeneracy and alters the stability and geometry of these textures in a 3D bulk environment.

Methodology
The authors develop a theoretical framework where the magnetization field $\mathbf{m} = M_s \mathbf{n}$ (with $|\mathbf{n}|=1$ ) varies smoothly in a 3D bulk system. To make the problem tractable, they impose translational invariance along the $z$ -direction (the direction of the applied magnetic field), effectively reducing the domain to a 2D plane $\mathbb{R}^2$ while retaining 3D magnetostatic physics.

Energy Functional: The total energy $E$ comprises exchange ( $E_{exch}$ ), DMI ( $E_{DMI}$ ), potential ( $E_{pot}$ including anisotropy and Zeeman terms), and the dipole-dipole interaction ( $E_{DDI}$ ).
Magnetostatic Self-Interaction: The DDI is treated as a non-local term derived from the magnetostatic potential $\psi$ , which satisfies Poisson's equation $\Delta \psi = \mu_0 \rho$ , where the source $\rho = -\nabla \cdot \mathbf{m}$ acts as a "magnetic charge density." The interaction energy is formulated as the electrostatic self-energy of this charge distribution.
Scaling and Stability: The authors demonstrate that the DDI term scales similarly to the potential energy term under spatial rescaling. This allows the system to evade the Hobart–Derrick theorem, which typically prohibits stable static solitons in pure exchange-DMI systems without anisotropy or external fields.
Numerical Approach: The system is solved using a non-local numerical relaxation method.
- Magnetization: An "arrested Newton flow" method (a fourth-order Runge-Kutta scheme with flow-arresting criteria) is used to minimize the energy functional.
- Potential: The Poisson equation is solved at each step using a nonlinear conjugate gradient method (Polak-Ribiere-Polyak).
- Validation: Results are cross-verified against Landau–Lifshitz–Gilbert (LLG) dynamics, confirming consistency despite the LLG method being significantly slower.
Parameters: Simulations utilize micromagnetic parameters representative of soft ferromagnetic multilayers and Heusler alloys (e.g., CoFeB, Mn-based compounds), with specific values for exchange stiffness $J$ , DMI strength $D$ , saturation magnetization $M_s$ , and anisotropy $K_m$ .

Key Contributions

Formalism: The paper establishes a self-consistent framework for calculating the back-reaction of the demagnetizing field on bulk chiral magnets, moving beyond the thin-film approximation.
Gauge Symmetry Breaking: It rigorously demonstrates that while Dresselhaus, Rashba, and Heusler DMI terms are gauge-equivalent in the absence of DDI, the inclusion of DDI breaks this symmetry, leading to distinct energy landscapes for each texture.
Numerical Algorithm: The adaptation of arrested Newton flow combined with conjugate gradient descent for solving the coupled magnetization-Poisson system provides a robust and efficient tool for studying non-local magnetic interactions.

Results
The study analyzes the effect of DDI on isolated skyrmions and skyrmion crystals for the three DMI types:

Isolated Skyrmions:
- Bloch (Dresselhaus): The divergence of the magnetization $\nabla \cdot \mathbf{n}$ is zero (solenoidal). Consequently, the DDI has no effect on the Bloch skyrmion; its energy and size remain unchanged.
- Néel (Rashba): The DDI induces a radially symmetric demagnetizing field. The Néel skyrmion shrinks slightly (from 8.422 nm to 8.415 nm) but retains its axial symmetry.
- Antiskyrmion (Heusler): The DDI breaks the axial symmetry of the antiskyrmion. The magnetic charge density is not rotationally symmetric, leading to a loss of axial symmetry in the relaxed solution.
Skyrmion Crystals (Lattice Stability):
- Bloch and Néel: In the absence of DDI, these skyrmions repel each other, favoring infinite separation. With DDI, they still prefer infinite separation (negative binding energy), though the Néel lattice exhibits slightly higher energy than the Bloch lattice due to the size reduction.
- Heusler Antiskyrmions: The DDI introduces an attractive force that stabilizes a crystalline structure.
  - Symmetry Breaking: The antiskyrmions lose axial symmetry and relax into a square-lattice arrangement (lemniscatic lattice) rather than a triangular one.
  - Binding Energy: Unlike the other types, the Heusler antiskyrmion crystal has a positive binding energy, with an optimal lattice spacing $L^* \approx 39$ nm that minimizes the total energy. The energy of the crystal approaches the energy of two isolated antiskyrmions from below as the lattice size increases.

Significance and Claims
The paper claims that long-range dipolar interactions are a critical factor in the stability of bulk chiral magnetic textures. Specifically, it demonstrates that DDI can:

Lift the energy degeneracy between different DMI-induced skyrmion types.
Break the axial symmetry of antiskyrmions.
Stabilize bulk antiskyrmion crystals in translationally invariant 3D chiral magnets, a phenomenon not observed in the absence of DDI.

The authors note that their findings are specific to the parameter regime chosen (representative of soft ferromagnets and Heusler compounds) and that the inclusion of an external magnetic field yields qualitatively similar results for their specific parameters, though they did not explore the zero-anisotropy, high-field regime where spin-stripe to skyrmion transitions are typically observed. The work provides a theoretical basis for understanding how magnetostatic self-interactions can drive the formation of square-lattice antiskyrmion crystals in bulk materials.