Effects of crystal field and momentum-based frustrated exchange interactions on multiorbital square skyrmion lattice

This paper provides a theoretical investigation into how multiorbital effects, crystal-field anisotropy, and momentum-based frustrated exchange interactions cooperatively stabilize square-shaped skyrmion lattices (S-SkL) in centrosymmetric tetragonal Ce-based magnets.

The Gist

Imagine you are trying to organize a massive, high-speed dance party in a giant square ballroom. Usually, in these "magnetic dance parties," the dancers (which are actually tiny particles called spins) like to move in circles or simple lines.

However, scientists have recently discovered a very rare and special dance move: the Square Skyrmion Lattice (S-SkL). Instead of swirling in messy circles, the dancers form perfectly organized, tiny, swirling "tornadoes" that are arranged in a neat, tight square grid.

This paper explores why these "square tornadoes" happen and how we can "choreograph" them. Here is the breakdown of their discovery using everyday analogies.

1. The Dancers and their "Personalities" (Multiorbital Effects)

In most studies, scientists treat these magnetic dancers as simple, identical robots. But in this paper, the authors look at Ce-based magnets (Cerium), where the dancers are much more complex.

Think of these dancers not as robots, but as people with different "outfits" (these are the orbitals). Some dancers prefer to face the ceiling (easy-axis), while others prefer to face the walls (easy-plane). Because these dancers have different "personalities" and outfits, they don't just react to each other; they interact in complex ways based on which outfit they are wearing. This "multiorbital" complexity is the secret ingredient that allows the square pattern to form.

2. The Music and the Rhythm (Frustrated Exchange)

The "music" in this ballroom is the exchange interaction—the force that tells one dancer how to move based on what their neighbor is doing.

Usually, the music has a simple beat (a single frequency). But in this paper, the music is "frustrated." Imagine a song that has a heavy bass beat, but also a subtle, high-pitched melody playing at a different rhythm. This is what the scientists call higher-harmonic wave vectors.

When the bass beat and the melody clash, the dancers get "frustrated." They can't just follow one simple rhythm. To satisfy both the heavy beat and the high melody at the same time, they settle into the most efficient pattern possible: the Square Skyrmion Lattice.

3. The "Choreography" (The Findings)

The researchers used supercomputers to simulate thousands of different "dance parties" by changing the music and the dancers' outfits. They found three main things:

• The Sweet Spot: The square tornado pattern only appears when the "outfits" (anisotropy) and the "music" (exchange) are perfectly balanced. If the dancers are too stubborn (too much anisotropy) or the music is too simple, the square pattern collapses.
• The "Imperfect" Dance (S-SkL′): Sometimes, if the dancers are too focused on facing the ceiling, the perfect square grid gets slightly squashed or tilted. It’s like a dance troupe that is almost perfectly in formation, but one row is slightly off-center.
• The Bubble Party (MBLs): If the music changes just a little bit, the tornadoes disappear and turn into "bubbles"—clumps of dancers that don't have that special "swirl" (topological charge) but still form a pattern.

Why does this matter? (The Big Picture)

Why do we care about tiny magnetic tornadoes in a square grid? Because these "tornadoes" are incredibly stable and can be moved around with very little energy.

In the future, we want to use these tiny magnetic patterns to store information in computers—much like how a hard drive uses magnetic bits, but much smaller and much faster. By understanding the "choreography" of these complex dancers, scientists are creating a roadmap for the next generation of ultra-powerful, energy-efficient technology.

Technical Summary

Technical Summary: Effects of Crystal Field and Momentum-Based Frustrated Exchange Interactions on Multiorbital Square Skyrmion Lattice

Authors: Yan S. Zha and Satoru Hayami
Date: April 28, 2026 (Preprint/Accepted version)

---

1. Problem Statement

The study addresses the microscopic stabilization mechanism of the square-shaped skyrmion lattice (S-SkL) in centrosymmetric tetragonal magnets, specifically targeting $f$-electron systems like Ce-based compounds.

While skyrmion lattices are well-understood in noncentrosymmetric magnets (stabilized by Dzyaloshinskii–Moriya interaction, DMI), S-SkLs in centrosymmetric systems are more elusive. Previous theoretical models suggested that S-SkLs could emerge from frustrated exchange interactions and higher-harmonic wave vectors. However, most existing research focused on "spin-only" models (e.g., Gd- or Eu-based systems) where orbital angular momentum is quenched. There is a significant knowledge gap regarding how multiorbital effects—specifically the interplay between crystal-field (CF) induced anisotropy, interorbital coupling, and orbital wave functions—contribute to the stability of these topologically nontrivial textures in $f$-electron materials with unquenched orbital moments.

2. Methodology

The authors employ a self-consistent mean-field calculation on a two-dimensional square lattice containing $Ce^{3+}$ ions ($f^1$ configuration). The theoretical framework includes:

• Model Hamiltonian: A localized $4f$ model incorporating:
  • Crystal-Field (CF) Effects: A tetragonal CF Hamiltonian that splits the $J=5/2$ manifold into three Kramers doublets, focusing on two low-lying $\Gamma_{7}$ doublets.
  • Exchange Interaction: A bilinear exchange interaction $J_{ij} \mathbf{J}_i \cdot \mathbf{J}_j$ in momentum space, parameterized by $J_1, J_2, J_3$ to introduce frustration.
  • Zeeman Coupling: Interaction with an external magnetic field $h$.
  • CF Splitting ($\Delta$): A parameter controlling the relative energy of the two $\Gamma_{7}$ doublets.
• Key Parameters for Systematic Investigation:
  • $\alpha$ (Superposition Coefficient): Controls the orbital wave function composition, thereby tuning the competition between intraorbital anisotropy and interorbital coupling.
  • $\gamma$ (Interorbital Weighting): A controlled coefficient used to tune the strength of the off-diagonal (interorbital) coupling.
  • $\xi$ (Higher-Harmonic Ratio): A dimensionless ratio $\xi = J_{Q'\eta} / J_{Q\eta}$ used to quantify the relative weight of higher-harmonic wave vectors (momentum-space frustration).
• Computational Approach: To ensure the global minimum is found, the authors use a massively parallel approach, evolving a diverse pool of initial magnetic-moment configurations (ansatz), including randomized and parameterized states.

3. Key Contributions

• Microscopic Decomposition: The paper provides a detailed energy decomposition ($U_{tot} = U_{\Gamma(1)} + U_{\Gamma(2)} + U_{coupling}$) to disentangle the energetic contributions of individual orbital subspaces and their coupling.
• Anisotropy Mapping: It establishes how the superposition coefficient $\alpha$ dictates whether the system exhibits easy-axis, easy-plane, or isotropic anisotropy within and between orbitals.
• Mechanism Identification: It identifies the "synergy" between crystal-field-induced anisotropy and momentum-based frustration as the primary driver for S-SkL stabilization.
• Phase Classification: The study provides a comprehensive "zoology" of multi-$Q$ states, including new or modified topological and nontopological textures.

4. Results

• S-SkL Stabilization: The S-SkL is robustly stabilized in intermediate magnetic fields. Its stability window is significantly expanded when there is a competition between intraorbital easy-plane anisotropy and interorbital easy-axis anisotropy (occurring at higher $\alpha$ values).
• Role of Multiorbital Effects:
  • The S-SkL is favored when the ground-state doublet is $\Gamma_{7}^{(1)}$.
  • As $|\Delta|$ becomes very large, the system reduces to a single-orbital model, and the S-SkL stability shrinks, proving that interorbital coupling is essential.
• Role of Higher Harmonics ($\xi$): A high $\xi$ (strong higher-harmonic contribution) is mandatory. When $\xi$ is small, the higher-harmonic modes are too "stiff" or energetically costly, and the S-SkL fails to form, resulting in simpler $1Q$ or $2Q$ states.
• Emergent Phases:
  • S-SkL$'$: A topologically nontrivial state with a slight breaking of fourfold rotational symmetry, appearing under strong easy-axis anisotropy.
  • Magnetic Bubble Lattices (MBLs): Both square (S-MBL I) and rectangular (R-MBL I) nontopological soliton lattices were identified.
  • Other $2Q$ States: A variety of states (2Q I–XIV) with varying degrees of scalar chirality and symmetry breaking were mapped.

5. Significance

This work provides the microscopic design principles for realizing skyrmion lattices in a much broader class of $f$-electron materials. By moving beyond spin-only models, the authors demonstrate that the orbital degrees of freedom in Ce-based compounds are not merely secondary but are fundamental to the emergence of complex topological textures. This offers a roadmap for experimentalists to search for S-SkLs in tetragonal Ce-based heavy-fermion systems by targeting specific ranges of crystal-field splitting and orbital compositions.