pith. sign in

arxiv: 2603.02997 · v1 · submitted 2026-03-03 · ❄️ cond-mat.other

Anisotropic skyrmion liquid phase

Pith reviewed 2026-05-15 17:39 UTC · model grok-4.3

classification ❄️ cond-mat.other
keywords skyrmionmelting transitionanisotropyhexatic phasephase transitiontwo-dimensional systemmagnetic quasiparticles
0
0 comments X

The pith

Anisotropic skyrmion interactions from the atomic lattice lead to a direct solid-liquid transition with persistent orientational order.

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

The paper investigates the melting of two-dimensional skyrmion ensembles under different interaction symmetries. For isotropic interactions, the transition follows the expected path through a hexatic phase to an isotropic liquid. When the atomic lattice imposes anisotropy on the interactions, however, the hexatic phase is suppressed and a direct transition to a liquid occurs where orientational order remains intact up to 30 K. This matters because it shows how microscopic lattice effects can stabilize partial order in what would otherwise be a disordered phase, potentially affecting the behavior of magnetic quasiparticles in real materials.

Core claim

In quasiparticle dynamical simulations of skyrmion systems, isotropic interactions produce a solid phase that melts first into a hexatic phase over a narrow temperature window and then into an isotropic liquid. Forcing the skyrmion-skyrmion interactions to be anisotropic, as induced by the atomic lattice, eliminates the hexatic phase and yields a direct solid-to-liquid transition. Orientational order in this anisotropic liquid phase persists up to temperatures of 30 K.

What carries the argument

Lattice-induced anisotropy in the pairwise interactions between skyrmions, which modifies the sequence of phases during melting according to KTHNY theory.

Load-bearing premise

The quasiparticle dynamical simulations correctly reproduce the anisotropy in skyrmion interactions induced by the atomic lattice without significant higher-order corrections.

What would settle it

Experimental imaging of skyrmion positions at varying temperatures showing either the presence or absence of a hexatic phase and the temperature at which orientational order vanishes in the liquid.

Figures

Figures reproduced from arXiv: 2603.02997 by Daniel Schick, Levente R\'ozsa, Thomas Mutschler, Tim Matthies, Ulrich Nowak.

Figure 1
Figure 1. Figure 1: FIG. 1. Interaction potential between two skyrmions in the [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Examples of simulations of skyrmion ensembles. [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Average number of topologically non-trivial defects as [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Structure factor [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Calculated correlation functions in the system. (a) and (d) Spatial orientational correlation function [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Phase diagrams obtained from the machine-learning [PITH_FULL_IMAGE:figures/full_fig_p007_8.png] view at source ↗
read the original abstract

The nature of the melting transition in two-dimensional systems of particles has attracted considerable research attention since the development of Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) theory. The hexatic phase proposed by this theory has been recently identified experimentally in ensembles of magnetic skyrmions, quasiparticles formed in a magnetically ordered crystal. Here, we use quasiparticle dynamical simulations to study how the anisotropy of the skyrmion-skyrmion interactions induced by the atomic lattice influences the melting transition. For isotropic interactions, we find a transition from a solid phase through a hexatic phase stable in a narrow temperature range to an isotropic liquid phase. However, if the interactions between skyrmions are forced to be anisotropic by the atomic lattice, then a direct solid-liquid transition can be observed with orientational order persisting up to temperatures of 30 K in the liquid phase.

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

1 major / 0 minor

Summary. The manuscript uses quasiparticle dynamical simulations to study the melting transition of 2D magnetic skyrmions. For isotropic skyrmion-skyrmion interactions it recovers the KTHNY sequence (solid-hexatic-isotropic liquid), while lattice-induced anisotropy is reported to produce a direct solid-liquid transition in which orientational order survives in the liquid phase up to 30 K.

Significance. If the modeling assumptions hold, the result would show that lattice anisotropy can eliminate the hexatic window and stabilize long-range orientational order in a skyrmion liquid, offering a concrete mechanism that links microscopic lattice effects to macroscopic phase behavior in 2D particle systems.

major comments (1)
  1. [quasiparticle model and simulation protocol] The central claim of a direct solid-liquid transition and 30 K orientational order in the anisotropic liquid rests on the quasiparticle interaction potential faithfully encoding lattice-induced anisotropy without higher-order corrections. No derivation, parameter values, or comparison to micromagnetic calculations is supplied to establish this, leaving open the possibility that omitted multipole or lattice-relaxation terms restore effective isotropy and reintroduce a hexatic regime.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address the major concern below and will revise the manuscript to incorporate additional details that strengthen the justification of our quasiparticle model.

read point-by-point responses
  1. Referee: The central claim of a direct solid-liquid transition and 30 K orientational order in the anisotropic liquid rests on the quasiparticle interaction potential faithfully encoding lattice-induced anisotropy without higher-order corrections. No derivation, parameter values, or comparison to micromagnetic calculations is supplied to establish this, leaving open the possibility that omitted multipole or lattice-relaxation terms restore effective isotropy and reintroduce a hexatic regime.

    Authors: We agree that a more explicit justification of the interaction potential is required. The anisotropic potential was obtained by integrating out the underlying spin lattice degrees of freedom while retaining the leading lattice-induced anisotropy; in the revised manuscript we will add an appendix containing the full derivation, tabulate all numerical parameters (including the anisotropy strength and cutoff radius), and present a side-by-side comparison of quasiparticle trajectories with micromagnetic simulations on small lattices. These additions will demonstrate that higher-order multipole and relaxation corrections remain small in the relevant temperature range and do not restore an effective isotropy or a hexatic window. revision: yes

Circularity Check

0 steps flagged

Quasiparticle dynamical simulations yield independent results on anisotropic melting without reduction to fitted inputs or self-citation chains

full rationale

The paper derives its central claims (direct solid-liquid transition and persistent orientational order to 30 K under lattice anisotropy) from explicit quasiparticle dynamical simulations that evolve skyrmion positions under an imposed anisotropic interaction potential. No parameter is fitted to the target transition temperatures or hexatic suppression; the anisotropy is introduced directly via the atomic lattice model, and the phase behavior emerges from the time evolution. Existing self-citations to prior skyrmion work are present but not load-bearing for the melting result, which rests on the simulation outputs rather than any imported uniqueness theorem, ansatz, or self-defined mapping. The derivation chain therefore remains self-contained against external benchmarks and does not collapse by construction to its inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The claim rests on the applicability of KTHNY theory to skyrmions and on the accuracy of the chosen interaction model; no new entities are postulated.

axioms (1)
  • domain assumption KTHNY theory describes the melting sequence for skyrmion ensembles
    The paper explicitly invokes the hexatic phase from KTHNY as the baseline for the isotropic case.

pith-pipeline@v0.9.0 · 5459 in / 1120 out tokens · 62534 ms · 2026-05-15T17:39:31.392666+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

44 extracted references · 44 canonical work pages

  1. [1]

    J. M. Kosterlitz and D. J. Thouless, Journal of Physics C: Solid State Physics6, 1181 (1973)

  2. [2]

    B. I. Halperin and D. R. Nelson, Phys. Rev. Lett.41, 121 (1978)

  3. [3]

    D. R. Nelson and B. I. Halperin, Phys. Rev. B19, 2457 (1979)

  4. [4]

    A. P. Young, Phys. Rev. B19, 1855 (1979)

  5. [5]

    Bladon and D

    P. Bladon and D. Frenkel, Phys. Rev. Lett.74, 2519 (1995)

  6. [6]

    E.P.BernardandW.Krauth,Phys.Rev.Lett.107,155704 (2011)

  7. [7]

    S. C. Kapfer and W. Krauth, Phys. Rev. Lett.114, 035702 (2015)

  8. [8]

    Wojciechowski and D

    K. Wojciechowski and D. Frenkel, CMST10, 235 (2004)

  9. [9]

    J. A. Anderson, J. Antonaglia, J. A. Millan, M. Engel, and S. C. Glotzer, Phys. Rev. X7, 021001 (2017)

  10. [10]

    Walsh and N

    L. Walsh and N. Menon, Journal of Statistical Mechanics: Theory and Experiment2016, 083302 (2016)

  11. [11]

    C. J. Guo, D. B. Mast, R. Mehrotra, Y. Z. Ruan, M. A. Stan, and A. J. Dahm, Phys. Rev. Lett.51, 1461 (1983)

  12. [12]

    N. N. Negulyaev, V. S. Stepanyuk, L. Niebergall, P. Bruno, M. Pivetta, M. Ternes, F. Patthey, and W.-D. Schneider, Phys. Rev. Lett.102, 246102 (2009)

  13. [13]

    Cheng, J

    M. Cheng, J. T. Ho, S. W. Hui, and R. Pindak, Phys. Rev. Lett.61, 550 (1988)

  14. [14]

    C. A. Murray and D. H. Van Winkle, Phys. Rev. Lett.58, 1200 (1987)

  15. [15]

    A. J. Armstrong, R. C. Mockler, and W. J. O’Sullivan, Journal of Physics: Condensed Matter1, 1707–1730 (1989)

  16. [16]

    A. H. Marcus and S. A. Rice, Phys. Rev. Lett.77, 2577 (1996)

  17. [17]

    K. Zahn, R. Lenke, and G. Maret, Phys. Rev. Lett.82, 2721 (1999)

  18. [18]

    Guillamón, H

    I. Guillamón, H. Suderow, A. Fernández-Pacheco, J. Sesé, R. Córdoba, J. M. De Teresa, M. R. Ibarra, and S. Vieira, Nature Physics5, 651–655 (2009)

  19. [19]

    I. Roy, S. Dutta, A. N. Roy Choudhury, S. Basistha, I. Mac- cari, S. Mandal, J. Jesudasan, V. Bagwe, C. Castellani, L. Benfatto, and P. Raychaudhuri, Phys. Rev. Lett.122, 047001 (2019)

  20. [20]

    Sharma, D

    R. Sharma, D. Rey, L. Longchambon, A. Perrin, H. Perrin, and R. Dubessy, Phys. Rev. Lett.133, 143401 (2024)

  21. [21]

    A. N. Bogdanov and D. A. Yablonski˘i, Sov. Phys. JETP 68, 101–103 (1989)

  22. [22]

    Bogdanov and A

    A. Bogdanov and A. Hubert, Journal of Magnetism and Magnetic Materials138, 255–269 (1994)

  23. [23]

    Mühlbauer, B

    S. Mühlbauer, B. Binz, F. Jonietz, C. Pfleiderer, A. Rosch, A. Neubauer, R. Georgii, and P. Böni, Science323, 915 (2009)

  24. [24]

    X. Z. Yu, Y. Onose, N. Kanazawa, J. H. Park, J. H. Han, Y. Matsui, N. Nagaosa, and Y. Tokura, Nature465, 901–904 (2010)

  25. [25]

    Seshadri and R

    R. Seshadri and R. M. Westervelt, Phys. Rev. B46, 5150 (1992)

  26. [26]

    Rosch, F

    P.Huang, T.Schönenberger, M.Cantoni, L.Heinen, A.Ma- grez, A. Rosch, F. Carbone, and H. M. Rønnow, Nature Nanotechnology15, 761–767 (2020)

  27. [27]

    Gruber, J

    R. Gruber, J. Rothörl, S. M. Fröhlich, M. A. Brems, F. Kammerbauer, M.-A. Syskaki, E. M. Jefremovas, S. Kr- ishnia, A. Sudbø, P. Virnau, and M. Kläui, Nature Nan- otechnology20, 1405–1411 (2025)

  28. [28]

    Reichhardt, C

    C. Reichhardt, C. J. O. Reichhardt, and M. V. Milošević, Rev. Mod. Phys.94, 035005 (2022)

  29. [29]

    Nishikawa, K

    Y. Nishikawa, K. Hukushima, and W. Krauth, Phys. Rev. B99, 064435 (2019)

  30. [30]

    Schick, M

    D. Schick, M. Weißenhofer, L. Rózsa, J. Rothörl, P. Virnau, and U. Nowak, Phys. Rev. Res.6, 013097 (2024)

  31. [31]

    Dzyaloshinsky, J

    I. Dzyaloshinsky, J. Phys. Chem. Solids4, 241–255 (1958)

  32. [32]

    Moriya, Phys

    T. Moriya, Phys. Rev.120, 91 (1960)

  33. [33]

    S.-Z. Lin, C. Reichhardt, C. D. Batista, and A. Saxena, Phys. Rev. B87, 214419 (2013)

  34. [34]

    A. O. Leonov and M. Mostovoy, Nature Communications 6, 8275 (2015)

  35. [35]

    Rózsa, A

    L. Rózsa, A. Deák, E. Simon, R. Yanes, L. Udvardi, L. Szunyogh, and U. Nowak, Phys. Rev. Lett.117, 157205 (2016)

  36. [36]

    Zázvorka, F

    J. Zázvorka, F. Jakobs, D. Heinze, N. Keil, S. Kromin, S. Jaiswal, K. Litzius, G. Jakob, P. Virnau, D. Pinna, K. Everschor-Sitte, L. Rózsa, A. Donges, U. Nowak, and M. Kläui, Nature Nanotechnology14, 658–661 (2019)

  37. [37]

    A. A. Thiele, Phys. Rev. Lett.30, 230 (1973)

  38. [38]

    Weißenhofer and U

    M. Weißenhofer and U. Nowak, New Journal of Physics22, 103059 (2020)

  39. [39]

    Schick, M

    D. Schick, M. Weißenhofer, L. Rózsa, and U. Nowak, Phys. Rev. B103, 214417 (2021)

  40. [40]

    P. Keim, G. Maret, and H. H. von Grünberg, Phys. Rev. E 75, 031402 (2007)

  41. [41]

    S. Z. Lin, B. Zheng, and S. Trimper, Phys. Rev. E73, 066106 (2006)

  42. [42]

    J. B. McQueen, inProc. of 5th Berkeley Symposium on Math. Stat. and Prob.(1967) pp. 281–297

  43. [43]

    Pedregosa, G

    F. Pedregosa, G. Varoquaux, A. Gramfort, V. Michel, B. Thirion, O. Grisel, M. Blondel, P. Prettenhofer, R. Weiss, V. Dubourg, J. Vanderplas, A. Passos, D. Cour- napeau, M. Brucher, M. Perrot, and E. Duchesnay, Journal of Machine Learning Research12, 2825 (2011)

  44. [44]

    Rózsa, K

    L. Rózsa, K. Palotás, A. Deák, E. Simon, R. Yanes, L. Ud- vardi, L. Szunyogh, and U. Nowak, Phys. Rev. B95, 094423 (2017)