Pith. sign in

REVIEW 4 major objections 6 minor 50 references

A quasicrystalline lattice makes the skyrmion density — and with it the topological Hall conductivity — tunable continuously to zero by magnetic field.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-02 07:58 UTC pith:LB6WN77J

load-bearing objection A well-executed multiscale study whose headline claim—quasi-continuous skyrmion density tuning—is a property of the fitted effective model, not yet demonstrated in the full spin dynamics. the 4 major comments →

arxiv 2607.07948 v2 pith:LB6WN77J submitted 2026-07-08 cond-mat.str-el cond-mat.mes-hall

Tunable Emergent Gauge Fields from Skyrmions in a Quasicrystalline Lattice

Leandro M. Chinellato , Flavia A. G\'omez Albarrac\'in , Cristian D. Batista , Pablo S. Cornaglia This is my paper
classification cond-mat.str-el cond-mat.mes-hall
keywords skyrmionsquasicrystalDzyaloshinskii-Moriya interactiontopological Hall effectpinning hierarchyskyrmion densityemergent gauge fieldclassical spin model
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper studies skyrmions — small, topologically protected swirls of magnetic moments — placed on a two-dimensional quasicrystal made of square and triangular tiles. It claims that the quasicrystal's local environments create a hierarchy of pinning energies with many nearly degenerate minima, so skyrmions can sit at many different sites with only tiny energy differences. Because of this, as an external magnetic field approaches the saturation value, skyrmions are removed one pinning site at a time rather than all at once; the skyrmion density falls continuously to zero with a logarithmic dependence on field, unlike the abrupt first-order collapse seen on periodic lattices. Since the topological Hall conductivity of an electron bath coupled to the texture is controlled by the skyrmion density, this gives a magnetic-field knob that smoothly turns the Hall response off near saturation. The paper's claim, if correct, identifies quasicrystalline magnets as a platform for analog, field-controlled topological textures and a large magnetoelectric response.

Core claim

On the quasicrystalline lattice the fully polarized state is never an exact ground state; the imbalance of nearest-neighbor vectors leaves a small Dzyaloshinskii–Moriya torque, producing a quasi-fully polarized background. Single skyrmions stabilize at specific local patches, and their energies form a singular spectrum of quasi-degenerate minima because each patch recurs in the quasicrystal with exponentially shrinking energy splittings. Near the critical field H_c, the field acts as a chemical potential, and the exponentially short-ranged skyrmion–skyrmion repulsion, V(r)=V_0 e^{-r/λ_V}, balances the linear energy gain; minimizing the resulting point-particle energy functional yields a dens

What carries the argument

The central machinery is the effective point-particle model for skyrmions near saturation. Skyrmions are treated as occupying discrete pinning sites {R_i} of the quasicrystal, each with a field-dependent energy ε_s(R_i;H)=ε_0(R_i)+α(H-H_c), and interacting through an isotropic exponential repulsion V(|r|)=V_0 e^{-|r|/λ_V} extracted from two-skyrmion relaxations. The field H plays the role of a chemical potential; balancing it against the exponential repulsion in the dilute limit produces the singular equation of state n_sk ~ [ln(Λ/(H_c-H))]^{-2}. The underlying microscopic input is the exponentially localized skyrmion profile, whose decay length is set by the gapped magnon spectrum of the qu

Load-bearing premise

The load-bearing premise is the fitted effective model in which skyrmions are point particles with an isotropic exponential repulsion and a single, site-independent field slope α; if the repulsion is anisotropic or α varies strongly from site to site, the predicted smooth density law and Hall tuning break down.

What would settle it

Compute the full two-skyrmion interaction V(R_i-R_j) for all inequivalent pinning pairs near H_c, and the single-skyrmion energy ε_s(R_i;H) for each site as a function of field. If the interaction is strongly anisotropic or the slope α varies by more than a few percent across the A, B, and C patches, the inverse-square-logarithm density law and the smooth field tuning of σ_xy fail. Alternatively, a large-scale Monte Carlo simulation at system sizes where finite-size effects are controlled should show whether the skyrmion number follows 1/[ln(Λ/(H_c-H))]^2 or exhibits a jump.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

Share X Bluesky LinkedIn Reddit HN

If this is right

  • Approaching saturation on a quasicrystal, the equilibrium skyrmion density vanishes continuously instead of jumping, so the total topological charge can be set to any small value by choosing the field.
  • The topological Hall conductivity σ_xy tracks the skyrmion density and can be smoothly reduced to zero with small magnetic-field changes near H_c, giving a field-controlled Hall response.
  • The hierarchical pinning spectrum means many nearly degenerate skyrmion configurations exist, so small perturbations can select different configurations — a natural multistability for memory elements.
  • The effective point-particle model reproduces the low-energy skyrmion arrangements from the full spin simulations, supporting the view that pinning plus exponential repulsion governs the low-density phase.
  • The exponentially weak repulsion and linear chemical potential imply the density law n_sk ~ 1/[ln(Λ/(H_c-H))]^2, a direct corollary that can be tested in larger simulations.

Where Pith is reading between the lines

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

  • If the density law holds, near H_c a tiny field change produces a large relative change in skyrmion spacing; this amplification could be used in sensors or switches that respond to skyrmion-density-dependent transport or optical properties.
  • The same singular pinning hierarchy should also govern skyrmion dynamics: driving skyrmions across such a substrate is likely to show directional locking and field/density-dependent ordering transitions, analogous to vortices on quasiperiodic arrays — a testable prediction beyond the static results.
  • At finite temperature, thermal fluctuations will round the logarithmic divergence; identifying the crossover temperature at which the quasi-continuous behavior is smeared would give a concrete experimental target.
  • If the single-slope α and isotropic V(r) assumptions are relaxed, one expects the ideal inverse-square-logarithm law to be replaced by a broader distribution of site occupancies; measuring the variance in local skyrmion energies would show how close real quasicrystals are to the idealized regime.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 6 minor

Summary. The paper studies a classical Heisenberg model with Dzyaloshinskii–Moriya interactions on a two-dimensional quasicrystalline (square-triangle) lattice. It reports that the quasiperiodic lattice generates a hierarchy of single-skyrmion pinning energies, that this hierarchy leads to a quasi-continuous suppression of the skyrmion density as the saturation field H_c is approached, and that this density tunability translates into continuous control of the topological Hall conductivity of an attached itinerant electron bath. The quantitative centerpiece is an effective point-particle model in which skyrmions occupy discrete pinning sites with a linear-in-field site energy and an exponential repulsion; from these forms the paper derives the asymptotic law n_sk ~ 1/[ln(Λ/(H_c-H))]^2 (Eq. 15) and uses it to explain smooth Hall response. Full-spin Monte Carlo, single-skyrmion relaxations, and KPM magnon calculations are used as supporting evidence.

Significance. If the central claim were established, the paper would identify a new and potentially useful mechanism — intrinsic quasiperiodic pinning — for field-controlled continuous tuning of topological charge and Hall response. The paper has clear strengths: it uses multiple complementary methods (BFGS relaxation, full-spin Monte Carlo with error bars, effective-particle simulated annealing, Kubo-formula transport), it makes specific falsifiable predictions (Eq. 15, smooth σ_xy(H)), and the observation of a three-peak hierarchy in single-skyrmion energies (Fig. 5) is interesting and well documented. However, the load-bearing quantitative claims rest on assumptions in the effective model that are not validated by the full spin model, and the derivation of Eq. (15) does not in fact use the quasicrystalline energy hierarchy. These issues make the central conclusion, as stated in the abstract and Sec. VIII, currently unsupported.

major comments (4)
  1. [§VI, Eq. (15)] The density law n_sk ~ 1/[ln(Λ/(H_c-H))]^2 is derived solely from balancing a linear chemical potential α(H-H_c) against an exponential repulsion V_0 e^{-r/λ_V}. It does not use the singular hierarchy of pinning energies that the paper emphasizes; any dilute system with exponential repulsion and a linear chemical potential would give the same law. To support the claim that quasiperiodicity 'enables' this quasi-continuous suppression, the authors should either derive Eq. (15) from the actual distribution of ε_0(R_i) and show that the discrete site set realizes the law, or present a periodic-lattice analog showing a first-order transition under the same interaction and chemical potential. As it stands, Eq. (15) is a property of the assumed exponential interaction, not of the quasicrystalline hierarchy, and the contrast with periodic lattices is asserted rather than derived.
  2. [§V–§VI, Fig. 10] The central physical issue is equilibration across topological sectors. The full-spin Monte Carlo (green points, Fig. 10) shows skyrmions persisting above H_c up to H≈1.1, and Sec. V states that 'once formed, skyrmions are protected by a finite topological energy barrier.' The effective point-particle model (orange points) instead permits skyrmion creation and annihilation without any barrier, by construction, and the paper reports that it yields lower energies than the spin MC. Thus the 'equilibrium' branch with the sharp drop to zero at H_c is a property of the barrier-free effective model, not of the spin Hamiltonian. The abstract and Sec. VIII claim that the system 'enables' quasi-continuous suppression and smooth Hall control; this is not demonstrated for the actual spin model, whose dynamics and simulated-annealing thermodynamics retain a metastable skyrmion population. The authors
  3. [§VI, Eq. (13)] The universal linear field dependence ε_s(R_i;H)=ε_0(R_i)+α(H-H_c) with a single slope α is supported by only three sites (A, B, C) in Fig. 8(a). Given that the paper's main qualitative point is the singular, hierarchical distribution of ε_0, it is not justified to assume that α is site-independent across the entire hierarchy. If α varies from site to site, or if ε_s is nonlinear in H-H_c for some sites, the effective chemical potential picture changes and Eq. (15) may not hold. The authors should extract ε_s(R_i;H) over a statistically meaningful sample of pinning sites and verify linearity and universality of α, or assess how a distribution of slopes affects the density law.
  4. [§VII, Fig. 11] The Hall conductivity shown in Fig. 11 is computed from 'MC samples of the periodic approximant' (Sec. VII), which are the metastable spin-MC configurations (green branch of Fig. 10), not the effective-model equilibrium configurations that display the predicted quasi-continuous density drop. The smooth variation of σ_xy(H) near H_c is therefore inherited from the metastable branch, not from the equilibrium branch on which the central claim rests. The statement that 'the continuous tunability of the skyrmion number under an applied magnetic field translates into continuous control of the Hall conductivity' is not supported for the equilibrium branch unless the calculation is repeated on the effective-model configurations (orange branch) or on spin configurations with a controlled topological charge. Please clarify which configurations enter Fig. 11 and, if the equilibrium claim is intende
minor comments (6)
  1. [Fig. 4 caption] Typo: 'The lower panel shows shows the decay...' — remove the duplicated 'shows.'
  2. [Sec. II] 'spin S^z=1 boundary conditions' is unclear. Do the boundary spins have S^z fixed to +1? Please state explicitly.
  3. [Sec. VI, Eq. (14) and Supplemental Eq. (S3)] Notation is inconsistent: the main text uses λ_V for the interaction range (≈1.1) and λ for the single-skyrmion profile decay (≈0.51), but the Supplemental Material's Eq. (S3) uses λ for the interaction range. Please unify notation.
  4. [Sec. III] Section heading reads 'QUASI-FULL Y POLARIZED STATE' — likely a spacing typo for 'QUASI-FULLY.' Also, 'at finite H' in Sec. III is ambiguous (magnetic field vs. Hamiltonian); clarify.
  5. [Sec. VI] The claim that the effective model 'reproduces the main features' of the spin MC would be stronger if accompanied by a quantitative comparison (e.g., radial distribution functions or site-occupancy correlations) between the effective-model configurations and the spin-MC configurations at the same field. Currently the comparison is qualitative.
  6. [Eq. (19)] In the Kubo formula, the denominator should presumably be (E_n-E_m)^2+η^2 for the Lorentzian broadening; the expression as written has the η^2 term added outside the square. Please check the standard form.

Circularity Check

0 steps flagged

No significant circularity: Eq. (15) is a derived consequence of independently computed energetic inputs, not a refit of the density.

full rationale

The derivation chain is self-contained. The spin Hamiltonian Eq. (1) is the input; single-skyrmion energies and the pair interaction V(r) are computed by minimizing this Hamiltonian (BFGS), and V(r) is then fitted to an exponential form (Eq. 14) while ε_s(R_i;H) is taken linear in H-H_c (Eq. 13). The dilute-limit law n_sk ~ 1/[ln(Λ/(H_c-H))]^2 (Eq. 15) follows mathematically from those inputs; it is not obtained by fitting n_sk(H) data, so it is a model prediction rather than a tautology or a fitted parameter renamed as a result. The exponential form is also microphysically motivated by the gapped magnon spectrum. The Hall conductivity is computed from actual spin configurations through Eq. (19), not from Eq. (15) alone. Self-citations ([30], [34], [47]) are used as analogies, implementation details, or prior particle-model support and are not load-bearing; no uniqueness theorem or ansatz is imported from the authors' prior work. The paper explicitly flags its main physical limitation: Sec. V states that 'once formed, skyrmions are protected by a finite topological energy barrier,' and Fig. 10 shows the full-spin MC retaining skyrmions above H_c. This is an acknowledged equilibration/metastability concern that weakens the practical claim of continuous tunability, but it is not a circularity in the derivation.

Axiom & Free-Parameter Ledger

5 free parameters · 7 axioms · 0 invented entities

The free parameters are the microscopically fitted inputs of the effective model (interaction strength/range, chemical-potential slope, and the model parameters D and H_c). The axioms are standardness of the quasicrystal construction, classical-spin and no-back-reaction domain assumptions, and the ad hoc point-particle/exponential-repulsion assumption that carries the central density law.

free parameters (5)
  • DM interaction strength D/J = 1.08
    Chosen to make skyrmions a few lattice constants in size; all derived quantities (λ, λ_V, V_0) depend on it. The paper argues small-D/J physics is representative but shows no D scan.
  • Single-skyrmion profile fit A, λ = A=0.0337, λ=0.513
    Fit to 1-S_z(r) in Fig. 4; used to characterize exponential localization and to connect λ_V≈2λ.
  • Two-skyrmion interaction fit V_0, λ_V = V_0=671, λ_V≈1.1
    Fit to BFGS two-skyrmion energies (Fig. 8d), V(r)≈671 e^{-0.91r}. This exponential form drives the singular equation of state (15).
  • Chemical-potential slope α = not quoted (α>0)
    Slope of ε_s(R_i;H) vs H fitted to Fig. 8(a); assumed universal across all pinning sites in Eq. (13), a load-bearing assumption.
  • Critical field H_c = 0.9957
    Numerically determined as the field where a single skyrmion on the best (A-type) patch becomes unstable; the expansion Eq. (13) is anchored at H_c.
axioms (7)
  • standard math The square-triangle substitution tiling is quasiperiodic with Pisot eigenvalue, self-similar, and admits finite periodic approximants.
    Used in Sec. II to construct the lattice and the N=2911 periodic approximant with PBC.
  • domain assumption Spins are classical unit vectors; quantum fluctuations are neglected except for linear spin-wave analysis.
    Hamiltonian (1) is classical; LSWT is used only for the magnon spectrum, not for the skyrmion energetics.
  • domain assumption Above H_c the quasi-fully polarized state is well described by dropping the residual DM torque linear terms in the HP expansion.
    Sec. III argues the torque is 'parametrically small' but gives no quantitative bound; this underlies the magnon spectrum and the exponential interaction argument.
  • ad hoc to paper Skyrmions behave as point particles on discrete pinning sites with pairwise additive, isotropic, exponentially decaying repulsion.
    Assumed in Sec. VI effective model; the exponential form is fitted, not derived, and the pairwise additivity is not tested against three-skyrmion configurations.
  • ad hoc to paper Single-skyrmion energy is linear in field with a universal slope α for all pinning sites.
    Eq. (13): site-independence of α is not checked; a skyrmion on a different local environment could have a different Zeeman cost.
  • domain assumption Equilibrium (annealed) skyrmion densities are physically reachable despite topological protection barriers.
    Sec. V shows metastable skyrmions persist above H_c; the tunability claims assume these barriers are overcome in practice.
  • domain assumption Itinerant electrons do not back-react on the spin texture.
    Sec. VII: feedback neglected; strong-coupling Hund coupling is assumed, and only the pure spin texture generates the Hall response.

pith-pipeline@v1.3.0-alltime-deepseek · 18591 in / 18556 out tokens · 174195 ms · 2026-08-02T07:58:12.504620+00:00 · methodology

0 comments
Cite this review

Pith. "Pith review of Tunable Emergent Gauge Fields from Skyrmions in a Quasicrystalline Lattice." pith.science (2026). https://pith.science/paper/LB6WN77J

@misc{pith2026260707948,
  author       = {Pith},
  title        = {Pith review of: Tunable Emergent Gauge Fields from Skyrmions in a Quasicrystalline Lattice},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/LB6WN77J}},
  note         = {Machine review of arXiv:2607.07948}
}
Share X Bluesky LinkedIn Reddit HN
read the original abstract

We study magnetic skyrmions in a two-dimensional quasicrystalline lattice using a classical Heisenberg model with Dzyaloshinskii-Moriya interactions and an external magnetic field. The competition between the skyrmion-skyrmion repulsion and an emergent quasiperiodic pinning landscape gives rise to a sequence of distinct skyrmion lattice configurations as a function of field. The resulting hierarchy of quasiperiodic pinning potentials, characterized by closely spaced quasi-degenerate minima, enables a quasi-continuous suppression of the skyrmion density as the saturation field is approached, in sharp contrast to the strongly first-order collapse of skyrmion crystals on periodic lattices. This provides a direct mechanism for controlling the topological charge and, consequently, the emergent gauge field for itinerant electrons. As a consequence, the Hall conductivity can be strongly modified with small changes in the magnetic field and driven smoothly to zero near saturation. This field-controlled tunability, rooted in the underlying multistability, identifies quasicrystalline magnets as a platform for tunable topological textures, with potential applications in magnetic memory and magnetoelectronic response.

Figures

Figures reproduced from arXiv: 2607.07948 by Cristian D. Batista, Flavia A. G\'omez Albarrac\'in, Leandro M. Chinellato, Pablo S. Cornaglia.

Figure 1
Figure 1. Figure 1: FIG. 1. Inflation rules used to generate the quasicrystalline [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Dynamical spin structure factor in the fully polarized [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 2
Figure 2. Figure 2: IV. SINGLE-SKYRMION SOLUTIONS In addition to the quasi-fully polarized ground state, the system supports metastable single-skyrmion solu￾tions on the QFP state at fields H ≈ Hc [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Squared probability amplitude [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Single skyrmion on the quasicrystalline lattice at [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. (a) Single skyrmion energy as a function of the posi [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. a) Probability distribution of the angles between [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: (c) shows the interaction potential as a func￾tion of the relative position r. The potential is deter￾mined at the discrete relative positions at which pairs of skyrmions can be stabilized. The absence of data near (0, 0) reflects the fact that the interaction is too strong at these short distances for the skyrmions to re￾main pinned at their original positions. As expected, the interaction decreases expon… view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Spin configurations near the critical field, obtained [PITH_FULL_IMAGE:figures/full_fig_p009_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Average total skyrmion number [PITH_FULL_IMAGE:figures/full_fig_p009_10.png] view at source ↗
Figure 10
Figure 10. Figure 10: While metastable skyrmion configurations per [PITH_FULL_IMAGE:figures/full_fig_p010_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Hall conductivity [PITH_FULL_IMAGE:figures/full_fig_p011_11.png] view at source ↗

Review history (2 revisions) →

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

50 extracted references · 1 linked inside Pith

  1. [1]

    Y. Wu, B. Francisco, Z. Chen, W. Wang, Y. Zhang, C. Wan, X. Han, H. Chi, Y. Hou, A. Lodesani,et al., A van der waals interface hosting two groups of magnetic skyrmions, Advanced Materials34, 2110583 (2022)

    2022

  2. [2]

    Huang, D.-F

    K. Huang, D.-F. Shao, and E. Y. Tsymbal, Ferroelectric control of magnetic skyrmions in two-dimensional van der Waals heterostructures, Nano Letters22, 3349 (2022)

    2022

  3. [3]

    Kurumaji, T

    T. Kurumaji, T. Nakajima, M. Hirschberger, A. Kikkawa, Y. Yamasaki, H. Sagayama, H. Nakao, Y. Taguchi, T.-h. Arima, and Y. Tokura, Skyrmion lattice with a giant topological hall effect in a frustrated triangular-lattice magnet, Science365, 914 (2019)

    2019

  4. [4]

    J. A. M. Paddisonet al., Magnetic interactions of the centrosymmetric skyrmion material gd2pdsi3, Physical Review Letters129, 137202 (2022)

    2022

  5. [5]

    Hirschberger, T

    M. Hirschberger, T. Nakajima, S. Gao, L. Peng, A. Kikkawa, T. Kurumaji, M. Kriener, Y. Yamasaki, H. Sagayama, H. Nakao, K. Ohishi, K. Kakurai, Y. Taguchi, X. Yu, T.-h. Arima, and Y. Tokura, Skyrmion phase and competing magnetic orders on a breathing kagome lattice, Nature Communications10, 5831 (2019)

    2019

  6. [6]

    J. Mo, L. M. Chinellato, F. Williams, A. Kikkawa, J. A. Paddison, M. D. Frontzek, G. Sala, C. Pasco, K. Barros, T. Nakajima,et al., Skyrmion and meron crystals in in- termetallic gd3ru4al12: Microscopic model insights into chiral phases, arXiv preprint arXiv:2512.24071 (2025)

    arXiv 2025

  7. [7]

    Khatua and K.-Y

    J. Khatua and K.-Y. Choi, Frustration and chirality in three-dimensional trillium lattices: insights and perspec- tives, Journal of Physics: Condensed Matter37, 483001 (2025)

    2025

  8. [8]

    Y. Arai, K. Nakayama, A. Honma, S. Souma, D. Shiga, H. Kumigashira, T. Takahashi, K. Segawa, and T. Sato, Origin of multiple skyrmion phases in eual4, Nature Communications17, 3162 (2026)

    2026

  9. [9]

    J.-U. Lee, S. Lee, J. H. Ryoo, S. Kang, T. Y. Kim, P. Kim, C.-H. Park, J.-G. Park, and H. Cheong, Ising- type magnetic ordering in atomically thin FePS 3, Nano letters16, 7433 (2016)

    2016

  10. [10]

    K. S. Burch, D. Mandrus, and J.-G. Park, Magnetism in two-dimensional van der Waals materials, Nature563, 47 (2018)

    2018

  11. [11]

    Huang, M

    B. Huang, M. A. McGuire, A. F. May, D. Xiao, P. Jarillo- Herrero, and X. Xu, Emergent phenomena and proximity effects in two-dimensional magnets and heterostructures, Nature Materials19, 1276 (2020)

    2020

  12. [12]

    W. Chen, Z. Sun, Z. Wang, L. Gu, X. Xu, S. Wu, and C. Gao, Direct observation of van der waals stacking–dependent interlayer magnetism, Science366, 983 (2019)

    2019

  13. [13]

    Gibertini, M

    M. Gibertini, M. Koperski, A. F. Morpurgo, and K. S. Novoselov, Magnetic 2D materials and heterostructures, Nature nanotechnology14, 408 (2019)

    2019

  14. [14]

    Y. Xu, A. Ray, Y.-T. Shao, S. Jiang, K. Lee, D. We- ber, J. E. Goldberger, K. Watanabe, T. Taniguchi, D. A. Muller,et al., Coexisting ferromagnetic– antiferromagnetic state in twisted bilayer cri3, Nature nanotechnology17, 143 (2022)

    2022

  15. [15]

    Q. Tong, F. Liu, J. Xiao, and W. Yao, Skyrmions in the moir´ e of van der waals 2d magnets, Nano Letters18, 7194 (2018), pMID: 30285450, https://doi.org/10.1021/acs.nanolett.8b03315

    doi:10.1021/acs.nanolett.8b03315 2018

  16. [16]

    S. J. Ahn, P. Moon, T.-H. Kim, H.-W. Kim, H.-C. Shin, E. H. Kim, H. W. Cha, S.-J. Kahng, P. Kim, M. Koshino, et al., Dirac electrons in a dodecagonal graphene qua- sicrystal, Science361, 782 (2018)

    2018

  17. [17]

    W. Yao, E. Wang, C. Bao, Y. Zhang, K. Zhang, K. Bao, C. K. Chan, C. Chen, J. Avila, M. C. Asensio,et al., Quasicrystalline 30 twisted bilayer graphene as an incom- mensurate superlattice with strong interlayer coupling, Proceedings of the National Academy of Sciences115, 6928 (2018)

    2018

  18. [18]

    C.-Y. Hao, Z. Zhan, P. A. Pantale´ on, J.-Q. He, Y.-X. Zhao, K. Watanabe, T. Taniguchi, F. Guinea, and L. He, Robust flat bands in twisted trilayer graphene moir´ e qua- sicrystals, Nature Communications15, 8437 (2024)

    2024

  19. [19]

    Tamura, Y

    R. Tamura, Y. Muro, T. Hiroto, K. Nishimoto, and T. Takabatake, Long-range magnetic order in the qua- sicrystalline approximant cd 6 tb, Physical Review B—Condensed Matter and Materials Physics82, 220201 (2010)

    2010

  20. [20]

    M. G. Kim, G. Beutier, A. Kreyssig, T. Hiroto, T. Ya- mada, J. W. Kim, M. de Boissieu, R. Tamura, and A. I. Goldman, Antiferromagnetic order in the quasicrystal approximant cd 6tb studied by x-ray resonant magnetic scattering, Phys. Rev. B85, 134442 (2012)

    2012

  21. [21]

    Woodland, F

    L. Woodland, F. Labib, D. Khalyavin, P. Manuel, F. Orlandi, H. Luetkens, and R. Tamura, Magnetic ground state of a prototype quasicrystal approximant: a candidate for octahedral spin ice physics (2026), arXiv:2605.26044 [cond-mat.str-el]

    Pith/arXiv arXiv 2026

  22. [22]

    Mukherjee, A

    A. Mukherjee, A. B. Sanyal, and E. Dagotto, Unconven- tional skyrmions in an interfacial frustrated triangular lattice, Phys. Rev. B108, 014408 (2023)

    2023

  23. [23]

    Paredes, J

    R. Paredes, J. L. Arag´ on, and R. A. Barrio, Nonperiodic hexagonal square-triangle tilings, Phys. Rev. B58, 11990 (1998)

    1998

  24. [24]

    H. Lane, H. Zhang, D. Dahlbom, S. Quinn, R. D. Somma, M. Mourigal, C. D. Batista, and K. Barros, Kernel poly- nomial method for linear spin wave theory, SciPost Phys. 17, 145 (2024)

    2024

  25. [25]

    Dahlbom, H

    D. Dahlbom, H. Zhang, C. Miles, S. Quinn, A. Niraula, B. Thipe, M. Wilson, S. Matin, H. Mankad, S. Hahn, D. Pajerowski, S. Johnston, Z. Wang, H. Lane, Y. W. Li, X. Bai, M. Mourigal, C. D. Batista, and K. Barros, Sunny.jl: A julia package for spin dynamics, Journal of Open Source Software10, 8138 (2025)

    2025

  26. [27]

    Each magnetic moment is represented as a unit-length classical vector parametrized using a skew-symmetric matrix [48]

  27. [28]

    R. R. Curtin, M. Edel, R. G. Prabhu, S. Basak, Z. Lou, and C. Sanderson, The ensmallen library for flexible nu- merical optimization, Journal of Machine Learning Re- search22, 1 (2021)

    2021

  28. [29]

    Sanderson and R

    C. Sanderson and R. Curtin, Armadillo: a template- based C++ library for linear algebra, Journal of Open Source Software1, 26 (2016)

    2016

  29. [30]

    P. S. Cornaglia, M. Nu˜ nez, and D. Garcia, Unveiling ex- 13 otic magnetic phases in fibonacci quasicrystals through machine learning, Physical Review B108, 144429 (2023)

    2023

  30. [31]

    The skyrmion position is computed as the centroid weighted by 1−S z

  31. [32]

    Berg and M

    B. Berg and M. L¨ uscher, Definition and statistical dis- tributions of a topological number in the lattice o(3)σ- model, Nuclear Physics B190, 412 (1981)

    1981

  32. [33]

    Baake and U

    M. Baake and U. Grimm,Aperiodic order, Vol. 1 (Cam- bridge University Press, 2013)

    2013

  33. [34]

    P. S. Cornaglia, L. M. Chinellato, and C. D. Batista, Quasicrystalline chiral soliton lattices in a fibonacci heli- magnet, Physical Review B110, 054430 (2024)

    2024

  34. [35]

    Creutz, Overrelaxation and monte carlo simulation, Physical Review D36, 515 (1987)

    M. Creutz, Overrelaxation and monte carlo simulation, Physical Review D36, 515 (1987)

    1987

  35. [36]

    Bogdanov and A

    A. Bogdanov and A. Hubert, Thermodynamically stable magnetic vortex states in magnetic crystals, Journal of Magnetism and Magnetic Materials138, 255 (1994)

    1994

  36. [37]

    X. Z. Yu, Y. Onose, N. Kanazawa, J. H. Park, J. H. Han, Y. Matsui, N. Nagaosa, and Y. Tokura, Real-space ob- servation of a two-dimensional skyrmion crystal, Nature 465, 901 (2010)

    2010

  37. [38]

    Nagaosa and Y

    N. Nagaosa and Y. Tokura, Topological properties and dynamics of magnetic skyrmions, Nature Nanotechnol- ogy8, 899 (2013)

    2013

  38. [39]

    A. O. Leonov and M. Mostovoy, Multiply periodic states and isolated skyrmions in an anisotropic frustrated mag- net, Nature Communications6, 8275 (2015)

    2015

  39. [40]

    C. D. Batista, S.-Z. Lin, S. Hayami, and Y. Kamiya, Frus- tration and chiral orderings in correlated electron sys- tems, Reports on Progress in Physics79, 084504 (2016)

    2016

  40. [41]

    M. Raju, A. Yagil, A. P. Petrovi´ c, X. Zhang, H. Ishizuka, Y. Li, Y. Okamura, F. Kagawa, S. Seki, Y. Otani, N. Na- gaosa, and C. Panagopoulos, The evolution of skyrmions in ir/fe/co/pt multilayers and their topological hall sig- nature, Nature Communications10, 696 (2019)

    2019

  41. [42]

    Mandal, S

    A. Mandal, S. Satpathy, and B. R. K. Nanda, Tuning the band topology and topological hall effect in skyrmion crystals via the spin-orbit coupling, Physical Review B 112, 064428 (2025)

    2025

  42. [43]

    Blundell,Magnetism in condensed matter(OUP Ox- ford, 2001)

    S. Blundell,Magnetism in condensed matter(OUP Ox- ford, 2001)

    2001

  43. [44]

    De Nittis and M

    G. De Nittis and M. Lein,Linear response theory: an analytic-algebraic approach(Springer, 2017)

    2017

  44. [45]

    Reichhardt and C

    C. Reichhardt and C. Olson Reichhardt, Dynamical or- dering and directional locking for particles moving over quasicrystalline substrates, Physical review letters106, 060603 (2011)

    2011

  45. [46]

    C. O. Reichhardt and C. Reichhardt, Vortex dynamics and symmetry locking on quasiperiodic and periodic sub- strates, Physica C: Superconductivity479, 45 (2012)

    2012

  46. [47]

    S.-Z. Lin, C. Reichhardt, C. D. Batista, and A. Sax- ena, Particle model for skyrmions in metallic chiral mag- nets: Dynamics, pinning, and creep, Physical Review B—Condensed Matter and Materials Physics87, 214419 (2013)

    2013

  47. [48]

    A. V. Ivanov, V. M. Uzdin, and H. J´ onsson, Fast and ro- bust algorithm for energy minimization of spin systems applied in an analysis of high temperature spin configu- rations in terms of skyrmion density, Computer Physics Communications260, 107749 (2021). Supplemental Material: Tunable Emergent Gauge Fields from Skyrmions in a Quasicrystalline Lattice ...

    2021

  48. [49]

    Kramer and A

    B. Kramer and A. MacKinnon, Localization: theory and experiment, Reports on Progress in Physics56, 1469 (1993)

    1993

  49. [50]

    J.-B. Suck, M. Schreiber, and P. H¨ aussler,Quasicrystals: An introduction to structure, physical properties and ap- plications, Vol. 55 (Springer Science & Business Media, 2013)

    2013

  50. [51]

    A. Koga, S. Sakai, Y. Matsushita, and T. Ishimasa, Hype- runiform properties of square-triangle tilings, Phys. Rev. B110, 094208 (2024)

    2024