pith. sign in

arxiv: 2604.18079 · v2 · submitted 2026-04-20 · ❄️ cond-mat.mes-hall

Generic skyrmion phase diagram in ferrimagnetic films

M. V. Wijethunga , X. R. Wang This is my paper

Pith reviewed 2026-05-10 03:55 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall
keywords ferrimagnetic skyrmionsinter-sublattice exchangeDzyaloshinskii-Moriya interactionskyrmion phase diagramstrong-coupling regimedimensionless parametermicromagnetic simulations
0
0 comments X

The pith

Strong inter-sublattice exchange stabilizes skyrmions on DMI-free sublattices in ferrimagnets via locking.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper maps how the strength of inter-sublattice exchange coupling J controls relaxed skyrmion configurations in chiral ferrimagnetic films at zero field. When J is large, the sublattices lock rigidly so that a skyrmion forms stably on the sublattice that has no intrinsic Dzyaloshinskii-Moriya interaction, because the twisting energy is supplied through the locked partner that does carry DMI. The authors introduce an effective dimensionless parameter ζ_eff that locates the crossover from this locked regime to the weak-J limit where the sublattices decouple and behave independently. A reader would care because ferrimagnetic skyrmions combine low net magnetization with chiral stabilization, and the work explains when the conventional stability parameter κ remains useful and when it fails. The result is a phase diagram that ties J directly to the existence of stable and metastable skyrmion states.

Core claim

In the strong-coupling regime, inter-sublattice locking enables stabilization of skyrmion in a sublattice where intrinsic Dzyaloshinskii-Moriya interaction is absent while the other sublattice has finite DMI, yielding a sublattice with DMI-free ferrimagnetic skyrmions. As J decreases, this locking breaks down, leading to independent sublattice behavior and the failure of an effective κ-based description. Our results establish a unified framework linking inter-sublattice exchange and skyrmion phase stability in ferrimagnetic systems.

What carries the argument

The dimensionless parameter ζ_eff that locates the crossover between strong and weak inter-sublattice locking.

If this is right

  • Skyrmions remain stable on a DMI-absent sublattice when inter-sublattice locking is strong.
  • Sublattices act independently and the κ description fails once J drops below the crossover set by ζ_eff.
  • Both stable and metastable skyrmion configurations are determined by the value of J.
  • A single effective parameter ζ_eff unifies the description of skyrmion phases across the strong-to-weak coupling range.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Engineering materials with DMI present on only one sublattice could still host skyrmions if J is kept sufficiently large.
  • Temperature tuning of J might allow reversible switching between locked and unlocked skyrmion regimes in the same sample.
  • The locking mechanism may apply to other chiral textures such as domain walls or merons in similar two-sublattice systems.

Load-bearing premise

The micromagnetic energy functional and chosen discretization remain valid across the full range of J and that the relaxed configurations are the true ground or metastable states.

What would settle it

Imaging that shows either a DMI-free skyrmion texture on one sublattice at high J, or fully independent skyrmion lattices on each sublattice at low J.

Figures

Figures reproduced from arXiv: 2604.18079 by M. V. Wijethunga, X. R. Wang.

Figure 1
Figure 1. Figure 1: FIG. 1. (Meta)stable skyrmions for di [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Normalized order parameter, [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. (Meta)stable skyrmions in the strong coupling regime for di [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. (Meta)stable skyrmions in the strong coupling regime for di [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. (Meta)stable skyrmions for various [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Phase diagram of isolated and condensed skyrmions in the [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. (Meta)stable skyrmions for weak coupling regime. Skyrmions for ( [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Phase diagram of isolated and condensed skyrmions in [PITH_FULL_IMAGE:figures/full_fig_p013_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. (Meta)stable states in the finite DMI ( [PITH_FULL_IMAGE:figures/full_fig_p014_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Normalized order parameter, [PITH_FULL_IMAGE:figures/full_fig_p015_10.png] view at source ↗
read the original abstract

Ferrimagnetic skyrmions offer enhanced tunability due to antiferromagnetically coupled sublattices and reduced net magnetization. In chiral magnetic films at zero magnetic field, skyrmion stability is commonly characterized by a dimensionless parameter $\kappa$, yet its applicability to ferrimagnetic systems remains unclear, as most studies assume a fixed, strong inter-sublattice exchange coupling $J$. Here we investigate how variations in $J$ govern relaxed stable and metastable ferrimagnetic skyrmion configurations and introduce a dimensionless parameter $\zeta_{eff}$ to characterize the crossover between strong and weak inter-sublattice locking. In the strong-coupling regime, inter-sublattice locking enables stabilization of skyrmion in a sublattice where intrinsic Dzyaloshinskii-Moriya interaction is absent while the other sublattice has finite DMI, yielding a sublattice with DMI-free ferrimagnetic skyrmions. As $J$ decreases, this locking breaks down, leading to independent sublattice behavior and the failure of an effective $\kappa$-based description. Our results establish a unified framework linking inter-sublattice exchange and skyrmion phase stability in ferrimagnetic systems.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper numerically explores skyrmion stability in two-sublattice ferrimagnetic films where one sublattice has finite DMI and the other has none. By tuning the inter-sublattice exchange J, the authors identify a crossover between strong-coupling (locked) and weak-coupling (decoupled) regimes, parameterized by a new dimensionless quantity ζ_eff. In the strong-coupling limit they report that inter-sublattice locking transfers effective DMI to the nominally DMI-free sublattice, stabilizing skyrmions there; as J decreases the locking fails and the effective-κ description breaks down. A phase diagram in the (κ, ζ_eff) plane is presented.

Significance. If the reported configurations are true minima, the work supplies a practical dimensionless criterion (ζ_eff) that unifies skyrmion phenomenology across coupling strengths and demonstrates a mechanism for DMI-free sublattice skyrmions. This is directly relevant to ferrimagnetic spintronics where net magnetization is low. The explicit construction of ζ_eff from micromagnetic parameters is a clear methodological contribution.

major comments (3)
  1. [Numerical Methods] Numerical Methods: The manuscript does not report whether the spatial discretization length is held fixed or scaled with J. Because the effective exchange length between sublattices shrinks as J decreases, a fixed mesh may under-resolve the decoupled regime; without mesh-convergence data at several J values below the reported crossover, the location of the ζ_eff transition and the claimed decoupling cannot be considered robust.
  2. [Results] Results (§4 or equivalent): No protocol is described for verifying that the relaxed states are global or metastable minima rather than simulation artifacts. Exhaustive sampling of initial conditions, energy-barrier estimates, or direct comparison with atomistic spin dynamics at representative J points would be required to substantiate the phase boundaries and the claim that locking enables DMI-free skyrmions.
  3. [Model / ζ_eff definition] Definition of ζ_eff: While ζ_eff is introduced as a combination of J and micromagnetic lengths, its numerical value at the observed crossover is not shown to be independent of the specific relaxation protocol or mesh; a sensitivity analysis would confirm that the parameter truly collapses the data across different system sizes.
minor comments (2)
  1. [Figures] Figure captions should explicitly state the normalization used for the magnetization profiles and whether the plotted quantity is the net or sublattice magnetization.
  2. [Abstract] The abstract states that κ-based descriptions fail below the crossover but does not quote the numerical value of ζ_eff at which this occurs; adding this single number would improve readability.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive comments that highlight important aspects of numerical robustness. We address each major point below and will revise the manuscript accordingly where appropriate.

read point-by-point responses

  1. Referee: [Numerical Methods] The manuscript does not report whether the spatial discretization length is held fixed or scaled with J. Because the effective exchange length between sublattices shrinks as J decreases, a fixed mesh may under-resolve the decoupled regime; without mesh-convergence data at several J values below the reported crossover, the location of the ζ_eff transition and the claimed decoupling cannot be considered robust.

    Authors: A fixed spatial discretization was employed, chosen to resolve the exchange length in the strong-coupling regime. We acknowledge the potential issue in the weak-coupling limit and will add mesh-convergence tests at multiple J values below the crossover, confirming that the ζ_eff transition and phase boundaries remain unchanged upon refinement. revision: yes

  2. Referee: [Results] No protocol is described for verifying that the relaxed states are global or metastable minima rather than simulation artifacts. Exhaustive sampling of initial conditions, energy-barrier estimates, or direct comparison with atomistic spin dynamics at representative J points would be required to substantiate the phase boundaries and the claim that locking enables DMI-free skyrmions.

    Authors: We employed conjugate-gradient energy minimization from multiple initial conditions (uniform, random, and seeded skyrmion profiles on each sublattice), with consistent convergence to the reported states. The methods section will be expanded to describe this protocol explicitly. Exhaustive sampling and energy-barrier calculations were not performed due to computational cost, but the multi-initial-condition approach supports metastability. Direct atomistic comparisons lie outside the micromagnetic framework of the present study. revision: partial

  3. Referee: [Model / ζ_eff definition] While ζ_eff is introduced as a combination of J and micromagnetic lengths, its numerical value at the observed crossover is not shown to be independent of the specific relaxation protocol or mesh; a sensitivity analysis would confirm that the parameter truly collapses the data across different system sizes.

    Authors: We will add a sensitivity analysis in the revised manuscript (or supplementary material) varying mesh size and relaxation parameters, demonstrating that the critical ζ_eff value at crossover is robust and data collapse holds for the system sizes studied. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the derivation chain

full rationale

The paper's results are obtained from numerical energy minimization of a two-sublattice micromagnetic functional across a range of inter-sublattice exchange J. The dimensionless parameter ζ_eff is introduced as an a priori combination of J and other micromagnetic constants (exchange stiffness, DMI strength, etc.) to mark the strong-to-weak coupling crossover; it is not fitted to the observed skyrmion stability or existence data. Consequently the central claim—that strong locking induces effective DMI on the nominally DMI-free sublattice—follows from the simulated relaxed configurations rather than reducing to a tautology by construction. No load-bearing self-citations, uniqueness theorems, or ansatzes imported from prior work appear in the derivation, and the phase diagram is presented as the direct numerical outcome rather than a renaming of known patterns.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard micromagnetic energy minimization (exchange, DMI, anisotropy, Zeeman) plus numerical relaxation; ζ_eff is a derived dimensionless ratio rather than a new fitted constant or postulated entity.

axioms (2)
  • standard math The total energy is the sum of standard micromagnetic terms (exchange, DMI, anisotropy, demagnetization) with no additional ad-hoc potentials.
    Invoked implicitly when the authors speak of relaxed stable and metastable configurations.
  • domain assumption The inter-sublattice exchange J can be varied independently while all other material parameters remain fixed.
    Required to isolate the effect of locking strength.

pith-pipeline@v0.9.0 · 5508 in / 1595 out tokens · 71797 ms · 2026-05-10T03:55:10.637615+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

52 extracted references · 52 canonical work pages

  1. [1]

    A. N. Bogdanov and D. A. Yablonskii, Sov. Phys. JETP68, 101 (1989)

    work page 1989

  2. [2]

    Bogdanov and A

    A. Bogdanov and A. Hubert, J. Magn. Magn. Mater.138, 255- 269 (1994)

    work page 1994

  3. [3]

    X. S. Wang, H. Y . Yuan, and X. R. Wang, Commun. Phys.1, 31 (2018)

    work page 2018

  4. [4]

    A. Fert, V . Cros, and J. Sampaio, Nat. Nanotechnol.8, 152-156 (2013)

    work page 2013

  5. [5]

    Sampaio, V

    J. Sampaio, V . Cros, S. Rohart, A. Thiaville, and A. Fert, Nat. Nanotechnol.8, 839–844 (2013)

    work page 2013

  6. [6]

    C. Back, V . Cros, H. Ebert, K. Everschor-Sitte, A. Fert, M. Garst, T. Ma, S. Mankovsky, T. L. Monchesky, M. Mostovoy, et al., J. Phys. D: Appl. Phys.53, 363001 (2020)

    work page 2020

  7. [7]

    B. He, R. Tomasello, X. Luo, R. Zhang, Z. Nie, M. Carpentieri, X. Han, G. Finocchio, and G. Yu, Nano Lett.23, 9482–9490 (2023)

    work page 2023

  8. [8]

    Gnoli, F

    L. Gnoli, F. Riente, M. Vacca, M. Ruo Roch, and M. Graziano, Electronics10, 155 (2021)

    work page 2021

  9. [9]

    Hrabec, J

    A. Hrabec, J. Sampaio, M. Belmeguenai, I. Gross, R. Weil, S. M. Ch´erif, A. Stashkevich, V . Jacques, A. Thiaville, and S. Ro- hart, Nat. Commun.8, 15765 (2017)

    work page 2017

  10. [10]

    N. S. Kiselev, A. N. Bogdanov, R. Sch¨afer, and U. K. R¨oßler, J. Phys. D: Appl. Phys.44, 392001 (2011)

    work page 2011

  11. [11]

    S. S. P. Parkin and S. H. Yang, Nat. Nanotechnol.10, 195 (2015)

    work page 2015

  12. [12]

    Zhang, M

    X. Zhang, M. Ezawa, and Y . Zhou, Sci. Rep.5, 9400 (2015)

    work page 2015

  13. [13]

    Zhang, Y

    X. Zhang, Y . Zhou, M. Ezawa, G. P. Zhao, and W. Zhao, Sci. Rep.5, 11369 (2015)

    work page 2015

  14. [14]

    Z ´azvorka, F

    J. Z ´azvorka, F. Jakobs, D. Heinze, N. Keil, S. Kromin, S. Jaiswal, K. Litzius, G. Jakob, P. Virnau, D. Pinna,et al., Nat. Nanotechnol.14, 658 (2019)

    work page 2019

  15. [15]

    S. Luo, M. Song, X. Li, Y . Zhang, J. Hong, X. Yang, X. Zou, N. Xu, and L. You, Nano Lett.18, 1180 (2018)

    work page 2018

  16. [16]

    Zhang, G

    X. Zhang, G. P. Zhao, H. Fangohr, J. P. Liu, W. X. Xia, J. Xia, and F. J. Morvan, Sci. Rep.5, 7643 (2015)

    work page 2015

  17. [17]

    A. K. Behera, C. Murapaka, S. Mallick, B. B. Singh, and S. Bedanta, J. Phys. D: Appl. Phys.54, 025001 (2020)

    work page 2020

  18. [18]

    D. Zhu, W. Kang, S. Li, Y . Huang, X. Zhang, Y . Zhou, and Y . Zhao, IEEE Trans. Electron Devices65, 87 (2018)

    work page 2018

  19. [19]

    Tomasello, E

    R. Tomasello, E. Martinez, R. Zivieri, L. Torres, M. Carpentieri, and G. Finocchio, Sci. Rep.4, 6784 (2014)

    work page 2014

  20. [20]

    M. G. Morshed, H. Vakili, and A. W. Ghosh, Phys. Rev. Appl. 17, 064019 (2022)

    work page 2022

  21. [21]

    Rohart and A

    S. Rohart and A. Thiaville, Phys. Rev. B88, 184422 (2013)

    work page 2013

  22. [22]

    Heinze, K

    S. Heinze, K. von Bergmann, M. Menzel, J. Brede, A. Kubet- zka, R. Wiesendanger, G. Bihlmayer, and S. Bl¨ugel, Nat. Phys. 7, 713-718 (2011)

    work page 2011

  23. [23]

    Romming, C

    N. Romming, C. Hanneken, M. Menzel, J. E. Bickel, B. Wolter, K. von Bergmann, A. Kubetzka, and R. Wiesendanger, Science 341, 636-639 (2013)

    work page 2013

  24. [24]

    B ¨uttner, I

    F. B ¨uttner, I. Lemesh, and G. S. D. Beach, Sci. Rep.8, 4464 (2018)

    work page 2018

  25. [25]

    Everschor-Sitte, J

    K. Everschor-Sitte, J. Masell, R. M. Reeve, and M. Kl ¨aui, J. Appl. Phys.124, 240901 (2018)

    work page 2018

  26. [26]

    K.-J. Kim, S. K. Kim, Y . Hirata, S.-H. Oh, T. Tono, D.-H. Kim, T. Okuno, W. S. Ham, S. Kim, G. Go, et al. Nat. Mater.16, 1187-1192 (2017)

    work page 2017

  27. [27]

    Finley and L

    J. Finley and L. Liu, Phys. Rev. Appl.6, 054001 (2016)

    work page 2016

  28. [28]

    Caretta, M

    L. Caretta, M. Mann, F. B ¨uttner, K. Ueda, B. Pfau, C. M. G¨unther, P. Hessing, A. Churikova, C. Klose, M. Schneider, et al., Nat. Nanotechnol.13, 1154-1160 (2018)

    work page 2018

  29. [29]

    S. K. Kim, G. S. D. Beach, K. J. Lee, T. Ono, T. Rasing, and H. Yang, Nat. Mater.21, 24–34 (2022)

    work page 2022

  30. [30]

    S. A. Siddiqui, J. Han, J. T. Finley, C. A. Ross, and L. Liu, Phys. Rev. Lett.121, 057701 (2018)

    work page 2018

  31. [31]

    Hirata, D

    Y . Hirata, D. H. Kim, S. K. Kim, D. K. Lee, S. H. Oh, D. Y . Kim, T. Nishimura, T. Okuno, Y . Futakawa, and H. Yoshikawa, et al., Nat. Nanotechnol.14, 232–236 (2019)

    work page 2019

  32. [32]

    J. C. T. Lee, J. J. Chess, S. A. Montoya, X. Shi, N. Tamura, S. K. Mishra, P. Fischer, B. J. McMorran, S. K. Sinha, and E. E. Fullerton, Appl. Phys. Lett.109, 022402 (2016)

    work page 2016

  33. [33]

    S. A. Montoya, S. Couture, J. J. Chess, J. C. T. Lee, N. Kent, D. Henze, S. K. Sinha, M.-Y . Im, S. D. Kevan, and P. Fischer, et al., Phys. Rev. B95, 024415 (2017)

    work page 2017

  34. [34]

    K. Chen, D. Lott, A. Philippi-Kobs, M. Weigand, C. Luo, and F. Radu, Nanoscale12, 18137–18143 (2020)

    work page 2020

  35. [35]

    Quessab, J.-W

    Y . Quessab, J.-W. Xu, C. T. Ma, W. Zhou, G. A. Riley, J. M. Shaw, H. T. Nembach, S. J. Poon, and A. D. Kent, Scientific Reports10, 7447 (2020)

    work page 2020

  36. [36]

    Heigl, S

    M. Heigl, S. Koraltan, M. Va ˇnatka, R. Kraft, C. Abert, C. V ogler, A. Semisalova, P. Che, A. Ullrich, and T. Schmidt, et al., Nat. Commun.12, 2611 (2021)

    work page 2021

  37. [37]

    Zhang, X

    J. Zhang, X. Zhang, H. Chen, Y . Guang, X. Zeng, G. Yu, S. Zhang, Y . Liu, J. Feng, and Y . Zhao, et al., Appl. Phys. Lett. 116, 142404 (2020)

    work page 2020

  38. [38]

    S. Woo, K. Litzius, B. Kr¨uger, M.-Y . Im, L. Caretta, K. Richter, M. Mann, A. Krone, R. M. Reeve, M. Weigandet al., Nature Communications9, 959 (2018)

    work page 2018

  39. [39]

    Barker and O

    J. Barker and O. A. Tretiakov, Phys. Rev. Lett.116, 147203 (2016)

    work page 2016

  40. [40]

    M. V . Wijethunga, X. Hu, and X. Wang, Commun. Phys.8, 123 (2025)

    work page 2025

  41. [41]

    X. R. Wang, X. C. Hu, and H. T. Wu, Commun. Phys.4, 142 (2021)

    work page 2021

  42. [42]

    X. C. Hu, H. T. Wu, and X. R. Wang, Nanoscale14, 7516-7529 (2022)

    work page 2022

  43. [43]

    H. T. Wu, X. C. Hu, K. Y . Jing, and X. R. Wang, Commun. Phys.4, 210 (2021). 18

    work page 2021

  44. [44]

    H. T. Wu, X. C. Hu, and X. R. Wang, Sci. China Phys. Mech. Astron.65, 247512 (2022)

    work page 2022

  45. [45]

    S.-Z. Lin, C. D. Batista, and A. Saxena, Phys. Rev. B89, 024415 (2014)

    work page 2014

  46. [46]

    Vansteenkiste, J

    A. Vansteenkiste, J. Leliaert, M. Dvornik, M. Helsen, F. G. Sanchez, and B. Van Waeyenberge, AIP Advances4, 107133 (2014)

    work page 2014

  47. [47]

    J. Luo, J. Xia, X. Liang, Y . Zhou, L. Li, L. Gu, and G. Zhao, Phys. Rev. B110, 214409 (2024)

    work page 2024

  48. [48]

    K. Y . Jing, X. Gong, and X. R. Wang, Phys. Rev. B106, 174429 (2022)

    work page 2022

  49. [49]

    X. R. Wang, P. Yan, and J. Lu, Europhys. Lett.86, 67001 (2009)

    work page 2009

  50. [50]

    X. R. Wang, P. Yan, J. Lu, and C. He, Ann. Phys. (N.Y .)324, 1815–1820 (2009)

    work page 2009

  51. [51]

    C. T. Ma, T. Hartnett, W. Zhou, H. T. Nembach, S. J. Poon, and J. M. Shaw, Applied Physics Letters119, 192406 (2021)

    work page 2021

  52. [52]

    Legrand, D

    W. Legrand, D. Maccariello, F. Ajejas, S. Collin, A. Vecchi- ola, K. Bouzehouane, N. Reyren, V . Cros, and A. Fert, Nature Materials19, 34–42 (2020)

    work page 2020