Formation process of skyrmion lattice domain boundaries: The role of grain boundaries
Pith reviewed 2026-05-25 11:38 UTC · model grok-4.3
The pith
Grain boundaries stabilize skyrmion lattice domain boundaries at their intersections in FeGe.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors establish that skyrmion lattice domain boundaries in FeGe are stabilized specifically at the intersections of grain boundaries. Lorentz transmission electron microscopy observations demonstrate that grain boundaries and edges guide the placement of these SkL domain boundaries. Simulations solving the Landau-Lifshitz-Gilbert equation show that having a boundary between differently oriented SkL domains at grain intersections yields the minimum energy state. In contrast, across wide uniform areas without such intersections, SkL domain boundaries do not appear and the domains rotate to form a single aligned SkL domain.
What carries the argument
Stabilization of SkL domain boundaries at grain boundary intersections, as the lowest-energy configuration identified by micromagnetic simulation.
If this is right
- Grain boundaries can dictate the locations where SkL domain boundaries form and remain stable.
- In regions without grain intersections, SkL domains rotate and coalesce into a single domain rather than maintaining boundaries.
- The energy minimum favors a boundary pinned at a grain intersection over a continuous uniform SkL domain.
- Edges of the sample influence SkL domain boundary formation in a manner similar to grain boundaries.
Where Pith is reading between the lines
- Deliberate control of grain structure could enable targeted placement of skyrmion domain boundaries in devices.
- The pinning mechanism may operate in other chiral magnets that host skyrmion lattices.
- Varying average grain size could change the typical spacing or number of SkL domain boundaries that form.
Load-bearing premise
The micromagnetic simulation parameters accurately represent the real FeGe sample without being adjusted after the fact to match the observed domain patterns.
What would settle it
Direct observation in FeGe of SkL domain boundaries located away from grain intersections, or a simulation run with the reported material parameters that finds a uniform single-domain state lower in energy than a bounded configuration.
read the original abstract
We report on the formation process of skyrmion lattice (SkL) domain boundaries in FeGe using Lorentz transmission electron microscopy and small-angle electron diffraction. We observed that grain boundaries and edges play an important role in the formation of SkL domain boundaries; The SkL domain boundary is stabilized at the intersection of two grains. A micromagnetic simulation using the Landau-Lifshitz-Gilbert equation revealed that the SkL domains separated by a boundary represent the lowest energy configuration. Conversely, in a wide area, SkL domain boundaries were not formed and SkL domains with different orientations rotated to form a single SkL domain.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports on the formation of skyrmion lattice (SkL) domain boundaries in FeGe thin films. Using Lorentz transmission electron microscopy (LTEM) and small-angle electron diffraction (SAED), the authors observe that grain boundaries and sample edges play a key role, with SkL domain boundaries preferentially stabilized at intersections of two grains. Micromagnetic simulations solving the Landau-Lifshitz-Gilbert equation are presented as showing that SkL configurations containing a domain boundary are lower in energy than merged single-domain states; in larger defect-free regions the domains rotate and coalesce.
Significance. If the central claim holds, the work supplies direct experimental evidence that microstructural defects control SkL domain formation and supplies a plausible micromagnetic rationale. This is relevant to efforts to engineer skyrmion textures via defects. The experimental imaging component is a clear strength; the simulation component would add independent corroboration only if the material parameters are shown to be fixed from prior literature rather than adjusted to the observed images.
major comments (2)
- [micromagnetic simulation section] Micromagnetic simulation section: the manuscript does not state whether the exchange stiffness A, DMI constant D, and anisotropy K are taken from independent literature values for FeGe or were selected/adjusted to reproduce the specific domain-boundary orientations seen in the LTEM data. Because the claim that 'the SkL domains separated by a boundary represent the lowest energy configuration' is offered as independent support for the experimental observation, the provenance of these parameters is load-bearing; post-hoc fitting would render the energy comparison circular.
- [simulation results paragraph] Simulation results paragraph: no numerical energy differences, standard deviations, or tests against alternative initial conditions (e.g., random vs. seeded SkL orientations) are reported. Without these, the assertion that the boundary-containing state is the global minimum cannot be quantitatively assessed.
minor comments (2)
- [Abstract] Abstract: the phrase 'in a wide area' is undefined; the manuscript should specify the lateral scale or defect-density threshold that distinguishes the grain-intersection regime from the coalescence regime.
- [figures and methods] Figure captions and methods: ensure every LTEM image includes the applied field, temperature, and scale bar; the SAED patterns should indicate the scattering vector calibration.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the constructive comments on the micromagnetic simulations. We address each major comment below and will incorporate the requested clarifications in the revised version.
read point-by-point responses
-
Referee: [micromagnetic simulation section] Micromagnetic simulation section: the manuscript does not state whether the exchange stiffness A, DMI constant D, and anisotropy K are taken from independent literature values for FeGe or were selected/adjusted to reproduce the specific domain-boundary orientations seen in the LTEM data. Because the claim that 'the SkL domains separated by a boundary represent the lowest energy configuration' is offered as independent support for the experimental observation, the provenance of these parameters is load-bearing; post-hoc fitting would render the energy comparison circular.
Authors: We appreciate the referee drawing attention to this point. The values of A, D, and K employed in the simulations are taken directly from independent literature values for FeGe and were not fitted or adjusted to match the observed domain-boundary orientations in the LTEM images. We will revise the manuscript to state this explicitly, including the specific references, in the simulation methods section. This will confirm that the energy comparison provides independent corroboration rather than a circular argument. revision: yes
-
Referee: [simulation results paragraph] Simulation results paragraph: no numerical energy differences, standard deviations, or tests against alternative initial conditions (e.g., random vs. seeded SkL orientations) are reported. Without these, the assertion that the boundary-containing state is the global minimum cannot be quantitatively assessed.
Authors: We agree that quantitative reporting would strengthen the simulation results. In the revised manuscript we will add the computed energy differences between the boundary-containing and merged single-domain configurations. We will also describe the initial conditions used and note that multiple runs with both seeded orientations (consistent with experiment) and random initial SkL orientations were performed; the boundary-containing state was consistently obtained as the lower-energy configuration. Where relevant, standard deviations across runs will be included. revision: yes
Circularity Check
No circularity; simulation provides independent energy comparison
full rationale
The paper's chain consists of direct LTEM/SAED observations that SkL domain boundaries form at grain intersections, followed by an LLG micromagnetic simulation that finds the observed boundary-containing configuration to be lowest in energy. No equations, parameter-fitting steps, or self-citations are quoted that would make the energy result reduce to the experimental images by construction. The simulation parameters are described only as chosen to reproduce experimental conditions on FeGe; absent any demonstrated post-hoc adjustment to the specific domain-boundary data, the energy comparison remains an independent corroboration rather than a tautology. This matches the default expectation that most papers are non-circular.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The Landau-Lifshitz-Gilbert equation with standard micromagnetic terms accurately captures the energetics of skyrmion lattices in FeGe at the experimental temperature and field.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.