Skyrmion stacking in stray field-coupled ultrathin ferromagnetic multilayers

This paper derives a reduced energy functional for ultrathin ferromagnetic multilayers hosting skyrmions and demonstrates that in bilayer systems, the global energy minimizers consist of two antiparallel Néel skyrmions stabilized by stray field coupling.

The Gist

Imagine you are a tiny architect trying to build the most stable, energy-efficient tower possible, but instead of bricks, you are using magnetic spins. These spins are like microscopic compass needles that want to point in a specific direction.

In this paper, the authors are studying a very specific type of magnetic structure called a Skyrmion. Think of a skyrmion not as a solid brick, but as a tiny, swirling tornado of these compass needles. They are stable, particle-like whirls that could one day be used to store data in computers that use almost no electricity.

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

1. The Setup: A Stack of Pancakes

Imagine you have a stack of very thin, magnetic pancakes (ferromagnetic layers). In each pancake, you have placed one of these "magnetic tornadoes" (skyrmions).

• The Problem: Usually, these tornadoes are unstable and want to collapse or fly apart. In single-layer films, we often need special chemical tricks (called DMI) to hold them together.
• The Twist: The authors are looking at what happens when you stack these pancakes on top of each other with tiny gaps in between. They wanted to see if the magnetic fields leaking from one pancake to the next (called "stray fields") could act like glue to hold the tornadoes together without needing those extra chemical tricks.

2. The Discovery: The "Handshake" Effect

The authors did some heavy mathematical modeling (think of it as a super-advanced simulation) to figure out the energy landscape. They found a fascinating phenomenon:

When you have two layers, the skyrmions in the top layer and the bottom layer naturally want to lock hands.

• If the "swirl" in the bottom layer spins clockwise, the one in the top layer wants to spin counter-clockwise.
• More importantly, their "arms" (the parts of the spin pointing sideways) point in opposite directions.

The Analogy: Imagine two people standing on separate trampoline mats. If they both lean toward each other, they create a tension that stabilizes both of them. If they lean away, they might fall off. The paper proves that these magnetic tornadoes naturally want to lean toward each other in a specific, opposite way. This "leaning" creates a magnetic field that acts like a stabilizing force, holding the tornadoes in place even at room temperature.

3. The "Perfect Match" (The Math Part)

The paper uses complex calculus to prove two main things:

1. Existence: These stable, stacked tornadoes can exist. They aren't just a fantasy; the math guarantees that there is a "sweet spot" where the energy is lowest and the structure is happy.
2. The Shape: The most stable shape is when the two tornadoes are concentric (perfectly aligned on top of each other) and have opposite sideways spins.

They also calculated exactly how much energy it takes to keep them apart.

• Close together: They are strongly attracted to each other (like magnets snapping together).
• Far apart: They start to repel each other weakly, like two people who don't want to be too close.

4. Why This Matters for the Future

Why should you care about magnetic tornadoes?

• Super-Efficient Computers: Current computers use electricity to flip bits (0s and 1s), which generates heat and wastes energy. Skyrmions are like tiny, stable marbles that can be moved around with very little energy.
• Room Temperature: Most of these magnetic structures only work at freezing temperatures. This paper suggests that by stacking layers and using this "stray field glue," we might be able to make them stable at room temperature. This is the "holy grail" for making them practical for your phone or laptop.
• No Chemical Tricks Needed: Usually, you need to engineer the materials with specific chemical interfaces to get this stability. This paper shows that the physics of the stack itself might be enough, simplifying manufacturing.

Summary

The authors took a complex 3D magnetic problem and simplified it into a manageable 2D model. They proved that if you stack thin magnetic layers, the magnetic fields between them naturally force the skyrmions to pair up in a specific, stable dance. This "dance" (two concentric, oppositely spinning tornadoes) is the most energy-efficient state, offering a promising new path for building the next generation of ultra-low-power, high-density data storage devices.

In a nutshell: They found a way to use the magnetic "whispers" between stacked layers to glue tiny magnetic tornadoes together, making them stable enough to potentially power our future computers.

Technical Summary

1. Problem Statement and Motivation

The paper addresses the stability and energy landscape of magnetic skyrmions (topologically protected spin textures) in ultrathin ferromagnetic multilayers. While single-layer skyrmions are often stabilized by the Dzyaloshinskii-Moriya Interaction (DMI), this work focuses on stray field-coupled multilayers where the layers are separated by non-magnetic spacers.

• Context: These systems are promising for room-temperature spintronic applications.
• Challenge: In multilayers, the magnetization in adjacent layers interacts via long-range dipolar (stray) fields. This creates a strongly coupled nonlinear system.
• Gap: Existing models often rely on 2D approximations or assume strong interlayer exchange. This paper investigates the regime where interlayer exchange is negligible (due to spacers), and the coupling is purely magnetic (stray field), leading to fully 3D magnetization textures.
• Goal: To derive a rigorous reduced energy functional for $N$ stacked skyrmions, prove the existence of energy minimizers, and completely characterize the ground state for bilayers in the absence of DMI.

2. Methodology

The authors employ a multi-step mathematical approach combining micromagnetic modeling, asymptotic analysis, and variational calculus.

A. Micromagnetic Framework

The study begins with the full 3D micromagnetic energy functional $E(M)$, which includes:

• Intralayer exchange energy.
• Bulk and interfacial magnetocrystalline anisotropy (perpendicular).
• Zeeman energy (external field).
• Interfacial DMI.
• Stray field (demagnetizing) energy: The full non-local dipolar interaction.

The geometry consists of $N$ identical ferromagnetic layers of thickness $d$, separated by non-magnetic spacers, with total stack thickness much smaller than the characteristic variation scale of magnetization in the plane.

B. Asymptotic Reduction

The authors perform an asymptotic expansion in the limit where the layer thickness $d$ is much smaller than the exchange length ($\delta = d/\ell_{ex} \ll 1$).

• Assumption: Magnetization is uniform across the thickness of each layer but varies in the plane.
• Derivation: By expanding the stray field energy kernels (volume-volume, volume-surface, and surface-surface interactions), they derive a reduced 2D energy functional $E_N(\{m_n\})$.
• Key Finding: The expansion reveals a new local dipolar energy term corresponding to interlayer volume-surface interactions.
  • If in-plane magnetizations in adjacent layers are identical, this term vanishes.
  • If they are anti-parallel, this term acts mathematically equivalent to an effective interfacial DMI, favoring Néel-type rotation.

C. Finite-Dimensional Ansatz

Focusing on the regime dominated by intralayer exchange (conformal limit), the authors introduce an ansatz where each layer contains exactly one skyrmion. The magnetization profile in each layer is modeled as a truncated Belavin-Polyakov (BP) skyrmion characterized by:

• Position ($r_n$)
• Radius ($\rho_n$)
• Rotation angle ($\theta_n$)
• Truncation parameter ($L_n$)

This reduces the infinite-dimensional variational problem to a finite-dimensional energy function $F_N(\{\rho_n, \theta_n, r_n\})$.

D. Variational Analysis

The authors analyze the landscape of $F_N$:

1. Existence Proof: They prove that for fixed skyrmion positions, global energy minimizers exist for the radii and angles, preventing the skyrmions from collapsing to zero radius (a common issue in such variational problems).
2. Bilayer Characterization: They explicitly solve the minimization problem for the case $N=2$ (bilayers) with zero DMI ($\kappa=0$).

3. Key Contributions and Results

A. Derivation of the Reduced Energy Functional

The paper provides a rigorous derivation of the effective energy for ultrathin multilayers. A crucial contribution is identifying that stray field coupling between layers with anti-parallel in-plane magnetization generates an effective DMI. This explains why Néel skyrmions can be stabilized in bilayers even without intrinsic DMI.

B. Theorem 1: Existence of Minimizers

The authors prove that for any fixed set of skyrmion centers $\{r_n\}$, there exist optimal radii $\rho_n$ and angles $\theta_n$ that minimize the energy.

• Significance: This resolves the mathematical difficulty that energy could be lowered indefinitely by shrinking a skyrmion to a point. They show that the energy cost of shrinking (due to exchange and anisotropy) outweighs the gain from stray field interactions for small radii.

C. Theorem 2: Complete Characterization of Bilayer Ground State

For a bilayer ($N=2$) with no intrinsic DMI, the global energy minimizer is uniquely characterized as:

1. Concentric: The skyrmions in both layers share the same center ($r_1 = r_2$).
2. Anti-parallel In-Plane Magnetization: The in-plane components are opposite ($m_{\perp,1} = -m_{\perp,2}$).
3. Néel Type: The rotation angle is $\theta_1 = 0$ and $\theta_2 = -\pi$ (or vice versa), corresponding to Néel walls.
4. Chirality: The system exhibits a fixed chirality determined by the direction of magnetization at infinity.
5. Radius: The optimal radius $\rho$ is given explicitly in terms of the dimensionless stray field strength $\bar{\delta}$ using the Lambert W function.

D. Interaction Energy and Stability

• The authors calculate the interaction energy between two skyrmions as a function of their separation distance.
• Result: The interaction is strongly attractive at short distances, stabilizing the concentric Néel configuration.
• Transition: As separation increases, the system undergoes a transition where the skyrmions may switch to Bloch-type configurations (weaker interaction) and eventually exhibit weak repulsion at large distances.

E. Numerical Validation

The theoretical predictions were validated using MuMax3 micromagnetic simulations.

• Parameters: Simulated GdCo-like ferrimagnetic bilayers with negligible spacer thickness.
• Outcome: The simulations confirmed the formation of concentric, anti-parallel Néel skyrmions with a radius (~20.5 nm) matching the theoretical prediction derived from the reduced model.

4. Significance and Implications

• Mechanism for Room-Temperature Skyrmions: The work demonstrates that stray field coupling alone can stabilize compact Néel skyrmions in multilayers without requiring strong DMI. This is critical for designing materials where DMI is weak or difficult to engineer.
• Energy Efficiency: The "flux-closure" like structure (anti-parallel in-plane components) minimizes the magnetostatic energy, suggesting these configurations are energetically favorable for low-power spintronic devices.
• Theoretical Rigor: The paper bridges the gap between full 3D micromagnetics and simplified particle models, providing a mathematically proven framework for analyzing skyrmion stacks.
• Design Guidelines: The explicit formulas for optimal radius and the characterization of the energy landscape provide concrete guidelines for engineering multilayer heterostructures for next-generation memory and logic devices.

In summary, the paper establishes that in ultrathin ferromagnetic multilayers, the interplay of stray fields creates an effective DMI that stabilizes concentric, anti-parallel Néel skyrmions, offering a robust pathway for skyrmion-based technologies independent of intrinsic chiral interactions.